At 00:42:03 UTC Verbals announced an incoming 5.5M Earthquake from Tonga.
At 00:56:16:UTC H1 dropped from Observing due to a SQZr Unlock that was very likely caused by the increased ground motion.
We went back to Observing at 01:01:39 UTC after the SQZ Man relocked itself.
I did put the OPS mode into earthquake for that moment.
At 01:02:10 UTC Earthquake mode was activated.
TITLE: 08/22 Eve Shift: 2330-0500 UTC (1630-2200 PST), all times posted in UTC
STATE of H1: Observing at 154Mpc
OUTGOING OPERATOR: Oli
CURRENT ENVIRONMENT:
SEI_ENV state: CALM
Wind: 16mph Gusts, 9mph 3min avg
Primary useism: 0.03 μm/s
Secondary useism: 0.19 μm/s
QUICK SUMMARY:
I inherited a 15+ hour lock from the day shift operator.
Everything looks like its running well.
Wind forcast looks great as well.
TITLE: 08/22 Day Shift: 1430-2330 UTC (0730-1630 PST), all times posted in UTC
STATE of H1: Observing at 154Mpc
INCOMING OPERATOR: Tony
SHIFT SUMMARY: Observing at 155Mpc and have been Locked for over 15 hours. Nothing happened during my shift at all.
LOG:
14:30 UTC Observing and have been Locked for 6.5 hours
Start Time | System | Name | Location | Lazer_Haz | Task | Time End |
---|---|---|---|---|---|---|
15:37 | FAC | Kim | OpticsLab | n | Tech clean | 15:56 |
17:17 | SPI | Jeff | OpticsLab | n | SPI Inventory | 21:06 |
Leo, Camilla, Jennie W.
Over the summer, a robust mode-matching simulation for the SQZ beam to the OMC was built using the a la mode library. This report serves to summarize the current state of the simulation and the relevant information one can extract from it. Link to git code: https://git.ligo.org/leendert.schrader/alm-beam-simulation-for-sqz/-/tree/main SQZ_simulation_v4.m is currently the most recent version. #1 A detailed description of the squeezed beam path. Distances, curvature radii, incidence angles, refractive indices, etc. are all within the code file with comments on where they were obtained. This includes all information necessary to propagate the beam starting at the SQZ output (BM4), through ZM4-6, the OFI, SRM, and into the OMC. These can also be found in the google doc in T2500228, though the curvature radii for ZM4/5 are not what are used in the simulation. Related aLogs: 85282, 85339 #2 The beam characteristics between ZM5 and ZM6. Using the beam profiling data from SQZT7 one can plot a "q manifold" which assigns a q parameter with respect to strain gauge values for ZM4 and ZM5. That is, if one knows the strain gauge for ZM4/5, they can know the beam after ZM5. Attached is a sample of the q manifold obtained by fitting the horizontal beam width data. Related aLogs: 85917, 86365, 86519, 85775 #3 Simulated and experimental mode-matching at the OMC. Simulated: By Propagating the q manifold forward to the OMC using the beam path constructed in a la mode, we can simulate mode-matching with respect to ZM4/5 strain gauge values. The simulation was found to be accurate for high mode-matching, but underestimates matching on the fringes of the strain gauge ranges (compared to collected data). Experimental: Multiple OMC scans for mode-match were obtained, which allows us to characterize mode-matching across any strain gauge values for ZM4/5. Attached is a picture of both the simulated and experimental matches against each other. The surface with a grid fits the experimental data, while the surface without grid is the simulated data. Related aLogs: 86445, 86467, 86520 #4 ZM4 and ZM5 curvature estimates. We estimate ZM4 curvature by interpolating data in E2100289. Note that preload was increased from 46 in lb to 75 in lb, hence why we had to "interpolate." We predict the curvature of ZM4 acts as R_ZM4 = 2/(-0.026296*V_ZM4 - 0.310465) where V_ZM4 is the strain gauge value for ZM4. This equation is achieved by assuming a linear voltage - to - optical power relationship (hence the '2 divided by a linear term'). The q manifold is then used to predict the ZM5 curvature values. Therefore, by ranging the strain gauge in the simulation GUI, one can find simulated ranges for ZM4 and ZM5 curvature. Note: ZM4 was selected for the above estimate, as we found it behaved "more linearly" than ZM5. Miscellaneous: the simulation contains the unique SRM ray transfer matrix in the M_SRM variable (see the first interim progress report in T2500228 for details on computation). For all reports and presentations related to the summer project, see the main DCC document T2500228.
We have now been Locked for almost 13 hours and have been in Observing since I got in this morning. Nothing to report.
Fri Aug 22 10:08:42 2025 INFO: Fill completed in 8min 38secs
Closes FAMIS#26550, last checked 86343
Laser Status:
NPRO output power is 1.865W
AMP1 output power is 69.91W
AMP2 output power is 140.3W
NPRO watchdog is GREEN
AMP1 watchdog is GREEN
AMP2 watchdog is GREEN
PDWD watchdog is GREEN
PMC:
It has been locked 9 days, 22 hr 41 minutes
Reflected power = 23.7W
Transmitted power = 105.3W
PowerSum = 129.0W
FSS:
It has been locked for 0 days 7 hr and 57 min
TPD[V] = 0.8324V
ISS:
The diffracted power is around 4.0%
Last saturation event was 0 days 12 hours and 31 minutes ago
Possible Issues:
PMC reflected power is high
TITLE: 08/22 Day Shift: 1430-2330 UTC (0730-1630 PST), all times posted in UTC
STATE of H1: Observing at 153Mpc
OUTGOING OPERATOR: TJ
CURRENT ENVIRONMENT:
SEI_ENV state: CALM
Wind: 15mph Gusts, 10mph 3min avg
Primary useism: 0.02 μm/s
Secondary useism: 0.19 μm/s
QUICK SUMMARY:
Observing at 153Mpc and have been Locked for 6.5 hours.
Our normal LVK teamspeak is back up, so I've moved us back over!
TITLE: 08/22 Eve Shift: 2330-0500 UTC (1630-2200 PST), all times posted in UTC
STATE of H1: Observing at 149Mpc
INCOMING OPERATOR: TJ
SHIFT SUMMARY:
Uneventful shift up until a M7.5 EQ took over! Shift's almost over and I'm guessing we have another 45+min until we get out of LARGE_EQ mode. Will be switching over to automatic operations shortly and gave TJ a status text earlier tonight. Picket Fence continues to flash yellow.
LOG:
H1's been down about 45min so far due to a M7.5 EQ....waiting for the earth to calm down, but the R-Waves are barely rolling through now.
Just giving an update---It's been about almost 2hrs since the earthquake rolled thorugh (took down H1/L1/V1). Picket Fences has been flashing between yellow & red the last hour or so. Gave TJ (owl shift) a heads up. In the last 30min, EQ band has finally started finally slowly turning downward.
Additional Notes:
Jeff, Oli
Today we tested out the SR3 estimators for both Pitch and Yaw (86491). The ASC loops where we see a difference are DHARD, MICH, SRC1, and SRC2, since these all go through the SRC to their output ports. As expected, we don't see any effect on CHARD, PRC, or INP.
I've plotted comparisons of the ASC control and error signals during those times along with a reference time of earlier that morning.
The ultimate conclusion is that the SR3 P and Y estimators do work to make some of the ASC noise better, and the good thing is that in the places where we see that the estimator makes the noise worse, they all seem like easy improvements/fixes! The amount of improvement that we see doesn't look like much, but that's also expected since it's just SR3 that has the estimator out of the several suspensions in each ASC path. We expect that as we add more estimators along these paths, we will keep picking away at the suspension noise. We already see really impressive broadband improvement between 1-5 Hz on SRC1 and SRC2!
Main areas of concern:
- For all ASC loops and both P and Y, we see heightened SR3 resonances on around 0.7 Hz when the SR3 P estimator is on. This is because the OSEM P blend filter was made without including the 0.7-0.8 Hz mode. That is being added in and that should solve that issue
- Between 0.1 and 0.2 Hz for some loops, turning on the estimators puts in some extra noise from the GS13s, which aren't as sensitive at those lower frequencies. To solve this, we will be adding a high pass filter above those frequencies
- A few peaks aren't being damped as well with the estimator as they were with normal damping; however, now that we've shown that the estimator decreases our noise(86455, 86503), we are hoping to be able to increase gain somewhere in the estimator system (possibly in SR3_M1_EST_Y_FUSION_MEAS_BP) to better damp our resonances while still seeing improvements in noise as compared to normal damping
Here are the times used to grab the traces:
Estimator ON for SR3 P and Y (RED) 2025-08-21 16:35:00 UTC
Estimator ON for SR3 Y (BLUE) 2025-08-21 15:34:06 UTC
Estimator ON for SR3 P (GREEN) 2025-08-21 16:05:00 UTC
Estimator OFF (PINK) 2025-08-21 10:27:00 UTC
ASC | Pitch | Pitch Notes/Comments (besides 0.7 Hz SR3 P estimator resonances) |
Yaw | Yaw Notes/Comments (besides 0.7 Hz SR3 P estimator resonances) |
DHARD | Error | 0.6-0.8 Hz - SR3 P resonance elevated whenever P estimator is ON | Error | 0.2-0.4 Hz - Transient noise with only Y Est ON? 0.6-0.8 Hz - SR3 P resonance elevated whenever P estimator is ON |
Control | 0.6-0.8 Hz - SR3 P resonance elevated whenever P estimator is ON | Control | 0.2-0.4 Hz - Transient noise with only Y Est ON? 0.6-0.8 Hz - SR3 P resonance elevated whenever P estimator is ON |
|
MICH | Error | 0.6-0.8 Hz - SR3 P resonance elevated whenever P estimator is ON 1.1 Hz - Every estimator ON configuration lowers noise |
Error | 0.1-0.2 Hz - Estimators ON elevates noise; need to high pass GS13 0.3-0.4 Hz - Transient noise with only Y Est ON? |
Control | 0.6-0.8 Hz - SR3 P resonance elevated whenever P estimator is ON 1.0 Hz - Y estimator ON alone adds some noise, but P estimator fixes it 1.1-1.5 Hz - Every estimator ON configuration lowers noise |
Control | 0.1-0.2 Hz - Estimators ON elevates noise; need to high pass GS13 0.3-0.4 Hz - Transient noise with only Y Est ON? 0.6-0.8 Hz - SR3 P resonance elevated whenever P estimator is ON |
|
SRC1 | Error | 0.6-0.8 Hz - SR3 P resonance elevated whenever P estimator is ON 1.1 Hz - Every estimator ON configuration lowers noise |
Error | 0.15-0.4 Hz - Transient noise with only Y Est on? 0.6 Hz - Every estimator configuration lowers noise 0.6-0.8 Hz - SR3 P resonance elevated whenever P estimator is ON 1.0 Hz - Estimators ON elevate noise 1.2-3.1 Hz - Every estimator ON configuration lowers noise |
Control | 0.6-0.8 Hz - SR3 P resonance elevated whenever P estimator is ON 1.1 Hz - Every estimator ON configuration lowers noise 3.0-10.0 Hz - P and Y estimator damping lowers noise; this is good but we don't know why it's doing that there |
Control | 0.25-0.4 Hz - Transient noise with only Y Est on? 0.6-0.8 Hz - SR3 P resonance elevated whenever P estimator is ON 1.0 Hz - Estimators ON elevate noise 1.2-3.0 Hz - Every estimator ON configuration lowers noise 3.0-10.0 Hz - Just P estimator ON gives us the best noise performance, and P and Y estimators ON gives us the worst noise performance. This is NOT good, and also doesn't match with what we see with the SRC1 P Control signals, so we don't know why it looks like this at all |
|
SRC2 | Error | 0.1-0.2 Hz - Estimators ON elevates noise; need to high pass GS13 0.6-0.8 Hz - SR3 P resonance elevated whenever P estimator is ON 1.1-3.5 Hz - Every estimator ON configuration lowers noise |
Error | 0.6-0.8 Hz - SR3 P resonance elevated whenever P estimator is ON 1.0 Hz - Estimator ON elevates noise 1.1-2.1 Hz - Every estimator ON configuration lowers noise 2.1 Hz - Estimator ON elevates noise 2.5-5.0 Hz - Every estimator ON configuration lowers noise |
Control | 0.1-0.2 Hz - Estimators ON elevates noise; need to high pass GS13 0.6-0.8 Hz - SR3 P resonance elevated whenever P estimator is ON 1.1-3.5 Hz - Every estimator ON configuration lowers noise |
Control | 0.6-0.8 Hz - SR3 P resonance elevated whenever P estimator is ON 1.0 Hz - Estimator ON elevates noise 1.1-2.1 Hz - Every estimator ON configuration lowers noise 2.1 Hz - Estimator ON elevates noise 2.5-5.0 Hz - Every estimator ON configuration lowers noise |
I created an updated blend for the SR3 pitch estimator. These are in the SUS SVN next to the first version. (see LHO log 84452)
The design script is blend_SR3_pitchv2.m
the Foton update script is make_SR3_Pitch_blend_v2.m
these are both in {SUS_SVN}/HLTS/Common/FilterDesign/Estimator/ - revision 12608
The update script will install the new blends into FM2 of SR3_M1_EST_P_FUSION_{MEAS/MODL}_BP with the name pit_v2. pit_v1 should still be in FM1. Turn on FM2, turn FM1 off.
Jeff and Oli tried the first one, and see that the first 2 modes (about 0.65 and 0.75 Hz) are seeing more motion with the estimator damping than the normal damping. To correct this, _v2 add OSEM signal to the estimator for those modes. See plots below - First 2 are the _v2 blend and a zoom of the _v2 blend. figure 3 shows the measured pitch plant vs. OSEM path - the modes line up pretty well. There is a bit of shift because the peaks are close together. I expect hope this will not matter. Figure 4 shows the plant vs. the model path. Now all 4 modes are driven by the measured OSEM signal instead of the model.
It is interesting to see that the model was not doing a good job of predicting the motion at the first 2 peaks. This is (I guess) because either (a) the model and the plant are different - or - (b) there are unmodeled drives pushing the plant (the suspension) that the model doesn't know about.
I'm guessing the answer is (b - unmodeled drives) and is likely from DAC noise. I think this because
1 - The plant fit is smooth and really good.
2 - In the yaw analysis that Edgard and Ivey are doing (not yet posted) the first mode of the yaw plant can be seen with the OSEM, but the ISI motion is much too small to excite that level of motion. But the OSEM can see motion, so something is exciting that motion.
3 - The DAC noise is the only thing I can think of.
Quick chat with Jeff indicates that the DAC noise models at those frequencies are not well trusted. We'll try something anyway and see if it is close. I don't see how to use the estimator to deal with that noise - we'd need to have an accurate realtime measurement.
Updated make_SR3_Pitch_blend_v2.m to r12610 after fixing the filter name and the subblock it writes to. Will be loading these in the next time we lose lock.
Thanks Oli!
Also - as a note to myself - I've attached 1 sec of drive signal from the SR3 outputs at the time when both estimators were on (2025-08-21 16:40 UTC, see LHO log 86491 ). The pitch drive is about 30 to 50 counts pk-pk, and not particularly high frequency, compared to the 16384 model rate. This suggests that the low frequency DAC noise is worth following up, and also that it could theoretically be improved with whitening filters. However, since the DC levels are ~ 10k counts to hold the alignments, a simple gain probably wont work. plz note - I am NOT suggesting any changes here, just logging some observations for followup.
Overdue comment - The real problem seems to be that the L to P path was not installed, but it is now, see aLOG 86567.
We do need to look into the DAC noise, however. The extra motion has been fixed by adding the L2P path AND changing the blends.
At the beginning of this shift, Tony noticed that we had errors for the LHO Control Room Teamspeak. He mentioned this to Jonathan, and they found out that all the LIGO/MIT Teamspeak servers are down and that we should use the Virgo server to maintain communications. So---You can find the "LHO Control Room" teamspeak account at:
Of course you can always phone the Control Room as well.
TITLE: 08/21 Eve Shift: 2330-0500 UTC (1630-2200 PST), all times posted in UTC
STATE of H1: Observing at 150Mpc
OUTGOING OPERATOR: TJ
CURRENT ENVIRONMENT:
SEI_ENV state: SEISMON_ALERT
Wind: 12mph Gusts, 6mph 3min avg
Primary useism: 0.05 μm/s
Secondary useism: 0.17 μm/s
QUICK SUMMARY:
H1's been locked 3.25+hrs. TJ passed on what he needed to do for relocking earlier (I've not had to lock H1 for any of my last 6-evening shifts in a row this week...we'll see how tonight goes for my last (7th EVE) shift of the week!). Environmentally, we had a small EQ from Japan roll through, winds are low, and µseism has slowly popped above the 50th percentile in the last 24hrs.
Leo, Jennie, Camilla
Followed setup instructions from 80010, with 75mW injected in the SEED beam, we had 1mW on OPO_IR_PD_DC slightly lower than last time. Also took SQZ_FC to FC_MISALIGNED. And opened SQZ beam divertor and fast shutter. Have counts of ~60 on ASC-AS_A and B and 7e-4 on ASC-OMC-A and B NSUMs. Jennie took the OMC PZT down to zero to start and then ran a template userapps/../omc/h1/templates/OMC_scan_sqz_beam.xml.
Repeated mode scans with different ZM PSAMs offsets, the DC centering loops could account for the pitch change for PSAMs, for the extreme PSAMS values, I increased the servo limiter from 85 to 95.
Jennie and Leo have the data to analyze.
See the ndscope with the channels we used to monitor attached.
This template is saved in /ligo/home/jennifer.wright/Documents/OMC_scan/20250819_SQZ_beam_scan_with_ZM_changes.yaml
Following up on these OMC scans: attached is the table with all computed mode-mismatch values.
Leo, Jennie W., Camilla Below is a plot of the OMC scans fitted with a surface polynomial. The plot is from the presentation in T2500228, so the labels on axes will be different. This can be plotted on Matlab simply using the following code block (requires Curve Fitting Toolbox to use fit()). OMCX = [5.5, 4, 2.3, 3, 8.1, 2.3, 2.3, 6, 9.5, 6, 6, 9, 3, 8, 4, 9.5, 9.5]; OMCY = [-0.8, 0.34, 0.2, 0.85, -0.4, -0.1, 0.2, -0.4, -0.6, -4.5, -5.3, 2, 2, -2, -2, -5, -0.6]; OMCdata = [1-4.048/100, 1-3.412/100, 1-3.074/100, 1-3.234/100, 1-5.073/100, ... 1-3.176/100, 1-2.884/100, 1-2.174/100, 1-3.343/100, 1-9.490/100, 1-11.865/100, 1-6.282/100, 1-2.313/100, ... 1-3.013/100, 1-3.021/100, 1-9.179/100, 1-3.243/100]; omc = fit([transpose(OMCX), transpose(OMCY)], transpose(OMCdata), 'poly22'); p = plot(omc);
Leo, Jennie, Camilla, WP 12694
Jennie followed instructions for set up in 85775. We removed the SQZ beam iris at the bottom of the LPM (added for alignment capture during OFI vent work). Then took beam profiles in this SQZ path of the SEED beam with various PSAMS settings, adjusted PSAMS settings as in 85775 and used the servo for nominal settings. All data is attached. Jennie then reverted settings back to nominal.
Leo, Jennie W., Camilla The attached pdf contains all the beam q parameters fitted to the collected beam width data. Only the 13.5% data was fitted, as the D4S data was too inconsistent to obtain confident q values. Fitting was performed with the a la mode beamPath.fitBeamWidth function. The attached q parameters were individually plotted using a la mode and verified for their data-fitting accuracy. As mentioned in the document, all q parameters are located immediately after the interaction with ZM5 (through the view of BM4 -> ZM4 -> ZM5 beam travel).
Leo, Jennie W., Camilla Attached is a plot of the q manifold from the q parameter data, which allows for characterizing the beam smoothly with respect to ZM4/5 strain gauge voltage values. The image is taken from the presentation uploaded to T2500228. The real plot will likely have slightly different labels to axes. Link to git code for plotting: https://git.ligo.org/leendert.schrader/alm-beam-simulation-for-sqz/-/tree/main