The violin mode harmonics overlap with the ~15070Hz test mass bulk modes. This may explain the mystery modes that have been rung up in the vicinity at Hanford alog27659 15063Hz and Livingston LLO alog20100 15085Hz.
This spectrum was taken from these measurements alog27743. Two violin mode harmonics are visible between 15kHz to 15.1kHz and 15.5kHz and 15.6kHz. Several test mass modes are also visible. The violin modes were elevated and we were actuating on 15077Hz ITMX mode and 15072 ITMY mode (which are very large and far off the top of the plot).
Jenne, Sheila, Peter, Evan, Stefan
We had been using AS45I for SRM signal in the last two days, and they seemed to not work well for yaw today. Peter and Jenne noticed that the violin modes showed up more strongly in some quadratns of the AS45 WFS than others, so we drove a line in DARM to check the balancing of the AS45 WFS. We drove DARM at 21 Hz, and looked at the transfer functions to the individual quadrants before the phasing. We set the gains so that the signal amplitudes would all match quadrant 1. For AS A, we ended up with gains as different from 0 as 0.69
We then attempted to phase these signals in the same way that the LSC AS45 detector is phased, by turning on and off the DARM offset and making sure that all the DARM signal was in the Q phase. To do this we changed some phases by up to 20 degrees. We had to use a step size of 0.5degrees on AS A which is in loop for DHARD to avoid loosing lock.
After doing this we tried moving the SRM alingment to see if the AS45 signals were better for SRM. We saw that pit and yaw were cross coupled, and we didn't have much pitch signal at all. We reverted these changes since we don't really think that the gains in the electronics can be this badly matched.
C. Cahillane I have revamped the uncertainty budget to include covariances between all stages of actuation and all time-dependent parameters. I computed each parameter's covariances in real and imaginary coordinates to provide a consistent basis. I then compiled an6 x 6
Actuation Covariance MatrixC_A
, a2 x 2
Sensing Covariance MatrixC_S
, and an8 x 8
Kappa Covariance MatrixC_K
. Then I compile them into a giant covariance matrixC
:_ _ | C_A 0 0 | C = | 0 C_S 0 | |_ 0 0 C_K _|
Then, I multiply by some conspicuous Jacobian vectorsJ(f)
to get the final2 x 2
uncertainty matrix σ_R^2(f):σ_R^2 = J * C * J'
whereJ
looks like:_ _ | d Re(R) d Im(R) | | --------- --------- .... | | d Re(p_i) d Re(p_i) | J(f) = | | | d Re(R) d Im(R) | | --------- --------- .... | |_d Im(p_i) d Im(p_i) _|
(I was able to use complex differentiation and Cauchy-Riemann here to make the derivatives easier. Recall that R = 1/C + D*A. Now I can compute dR/dA = D and dR/dC = -1/C^2 to formJ(f)
, thanks to 200 year old mathematics) Finally, to make the uncertainties readable by humans, I divideσ_R^2(f)
by|R(f)|^2
, rotateσ_R^2(f)
byangle(R(f))
via a rotation matrix, and read off the square roots of the diagonal of the rotatedσ_R^2(f)
to get the magnitude and phase uncertainties plotted below. I have plotted the uncertainty at GPSTime = 1135136350, the time of the Boxing Day Event. The plot shows an overall increase in magnitude uncertainty of about 1% at low frequency. Phase uncertainty increased by about 0.5 degrees at low frequency. The effects are more dramatic at Livingston. Check out LLO aLOG 26542.
C. Cahillane I have reproduced the uncertainties including covariance for GW150914 for the calibration companion paper. We will have to update the associated uncertainty calculation sections of the paper. I have also attached two .txt files for the R_C01/R_C03 response comparison and the associated uncertainty. Something I failed to emphasize above: Our uncertainties in the response function are now fully covariant... the plots I show of the magnitude and phase are only approximations to the true uncertainty. I have looked at the 3D plots of the covariant ellipses, and it's a fairly good approximation in this case.
C. Cahillane
I have attached and printed my relative covariance matrix. Please see DCC T1600227 for an explanation of the relative covariance matrix.
Basically, the below is percentage covariances.
Re(A_U) Im(A_U) Re(A_P) Im(A_P) Re(A_T) Im(A_U) Re(C_R) Im(C_R) Re(K_T) Im(K_T) Re(K_P) Im(K_P) Re(K_C) Im(K_C) Re(f_C) Im(f_C)
Re(A_U) 0.0166 0.0083 0.0139 0.0079 0.0146 0.0067 0 0 0 0 0 0 0 0 0 0
Im(A_U) 0.0083 0.0209 0.0091 0.0169 0.0071 0.0178 0 0 0 0 0 0 0 0 0 0
Re(A_P) 0.0139 0.0091 0.0163 0.0052 0.0157 0.0066 0 0 0 0 0 0 0 0 0 0
Im(A_P) 0.0079 0.0169 0.0052 0.0181 0.0057 0.0156 0 0 0 0 0 0 0 0 0 0
Re(A_T) 0.0146 0.0071 0.0157 0.0057 0.0251 0.0047 0 0 0 0 0 0 0 0 0 0
Im(A_T) 0.0067 0.0178 0.0066 0.0156 0.0047 0.0187 0 0 0 0 0 0 0 0 0 0
Re(C_R) 0 0 0 0 0 0 0.0207 0.0079 0 0 0 0 0 0 0 0
Im(C_R) 0 0 0 0 0 0 0.0079 0.0208 0 0 0 0 0 0 0 0
Re(K_T) 0 0 0 0 0 0 0 0 0.0025 -0.0002 0.0019 -0.0018 -0.0004 0 0.0004 0
Im(K_T) 0 0 0 0 0 0 0 0 -0.0002 0.0025 0.0017 0.0019 0.0001 0 0.0001 0
Re(K_P) 0 0 0 0 0 0 0 0 0.0019 0.0017 0.0035 -0.0003 0.0002 0 -0.0003 0
Im(K_P) 0 0 0 0 0 0 0 0 -0.0018 0.0019 -0.0003 0.0035 0.0006 0 -0.0005 0
Re(K_C) 0 0 0 0 0 0 0 0 -0.0004 0.0001 0.0002 0.0006 0.0037 0 -0.0036 0
Im(K_C) 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
Re(f_C) 0 0 0 0 0 0 0 0 0.0004 0.0001 -0.0003 -0.0005 -0.0036 0 0.0054 0
Im(f_C) 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
Jenne, Peter, Jeff, Terra
We took charge measurements on the ITMs by driving 20.1 Hz into H1:SUS-ITMX/Y_L3_DRIVEALIGN_L2L_EXC with 100k cts and looked at the coupling to DARM. We stepped up and down the ESD bias voltage from zero and found the bias that gave zero coupling to get charge measurements, where Vcharge= (20/218) * bias * 40. ITMX zeroed with bias = 3k, ITMY zeroed with bias = 2.5k. See bias stepping ITMX and ITMY spectra attached.
ITMX charge: 9.15 V, ITMY charge: 7.6 V
ITMX:
BIAS offset | RMS AMPL (m) | PEAK AMPL (m) | PHASE (deg) |
---|---|---|---|
0 | 4.9x10-15 | 6.9x10-15 | 37 |
50K | 7.6x10-14 | 1.1x10-13 | -142 |
100K | 1.6x10-13 | 2.6x10-13 | -142 |
-50K | 8.5x10-14 | 1.2x10-13 | 38 |
-100K | 1.7x10-13 | 2.4x10-13 | 38 |
ITMY:
BIAS offset | RMS AMPL (m) | PEAK AMPL (m) | PHASE (deg) |
---|---|---|---|
0 | 5.4x 10-15 | 7.6x10-15 | -142 |
50K | 9.3x 10-14 | 1.3x10-13 | 38 |
100K | 1.9 x 10-13 | 2.7x10-13 | 38 |
-50K | 1.0 x 10-13 | 1.4x10-13 | -142 |
-100K | 2.0x 10-13 | 2.8x10-13 | -142 |
Approximating alpha: The first term in the expression for the force produced by the ESDs is the attractive force between the ESD fringe fields and the test mass: F = alpha(Vbias - Vsignal)2, where alpha is a constant of proportionality. With the known bias voltage Vb, signal drive voltage Vs, drive frequency f, and now the peak amplitude xpp , we used the largest bias offset (100k) to approximate alpha for both ITMs. Work is attached.
ITMX: alpha = 1.78 x 10-11 N/V2, ITMY: alpha = 1.85 x 10-11 N/V2
These LHO ITM force coefficients agree with LLO's. Using Valera and Den's 2015 measurements (and assuming an 80 Hz ESD drive), I calculated alpha for LLO ITMX: alpha = 1.46 x 10-11 N/V2
Stefan, Nutsinee
So we found out that the ITMY vstop filter has been turned off. That's likely the cause of mysterious ring up of ITMY violin modes last night and tonight. I was able to damp the three highest modes (ITMY MODE2, MODE3, MODE5) with the settings on the table I confirmed on June 13. I haven't had a chance to re-comfirm all the settings. So until then I would either comment out IY violin mode damping guardian lines or just skip it for now.
Nevermind. I commented out every IY line except for MODE2, 3, and 5.
Evan, Sheila
We tried the settings in the guardian (and in the new wiki) for ITMY mode 1 and mode 6, they worked so they are back in the guardian now. ETMX mode7 (507.195Hz) is the largest violin mode now, and I don't see settings for this mode on the wiki, so I tried FM1, FM4, and a +30 gain (BP, 100dB, 0 degrees of phase). This seemed to ring it up very slowly, so if we hadn't lost lock for unrelated reasons I would have tried a 60 degrees phase shift.
Wanted to get H1 to atleast DRMI before today's Press Conference, so jumped into locking first thing this morning, but no DRMI locks, and PRMI would not stay lock for more than a second, so opted for an Initial Alignment, but it wasn't trivial. A few notes from locking & alignment:
Another IFO Note: After the alignment, Krishna mentioned conditions were quiet so, we took the ITMs to the 90 blends and NO Sensor Correction.
BSC ISI Guardian Note:
Krishna noticed that Sensor Correction was turned off for ETMy & ETMx ISI's at around 1am last night by taking a gain from 1 to 0 by hand. We want to get in the habit of not doing this (because Guardian does not monitor these gain channels -at the moment-), and use our fancy new Guardian medm (sitemap/O-1/ISI Blend Filters.adl)
Daily activity log attached as pdf.
I checked the PSL chillers, per FAMIS work order 4155, assigned to me today, and found the chiller was filled yesterday, so did not need any water.
First post is in the ground and concrete is curing.
1505 - 1520 hrs. local -> To and from Y-mid Opened exhaust check-valve bypass valve. Opened LLCV bypass valve 1/2 turn -> LN2 at exhaust in < 1 minute -> Restored valves to as found Next over-fill to be Friday, June 17th.
Bad PEM corner station MAG and MIC channels
Richard, Fil, Robert, Dave:
We found that all the MAG and MIC corner stations channels have not been correct since the power recovery Mon 6/6. These channels span several ADC cards, so it did not look like a software issue. However to rule this out we power cycled h1oaf0, its IO Chassis and the PEM AA chassis with no change to the signals. Problem was identified as two DC power supplies which needed to be turned on after the outages (18V for MAG, 36V for MIC). Robert confirms these channels are now good. We are still unsure why all the MAG and MIC channels read about 3000 counts on the IOP and not zero, we may investigate further next Tuesday.
h1calcs new code
Joe, Jeff, Dave:
new h1calcs model went in. Some channels were promoted from standard single precision floating point to the new double precision floating point data type. The RCG appends the suffix "_DBL_DQ" to these channels, compared with "_DQ" for standard chans. The channels changed are:
[H1:CAL-DARM_CTRL_WHITEN_OUT_DBL_DQ]
[H1:CAL-DARM_ERR_WHITEN_OUT_DBL_DQ]
[H1:CAL-DELTAL_CTRL_DBL_DQ]
[H1:CAL-DELTAL_CTRL_PUM_DBL_DQ]
[H1:CAL-DELTAL_CTRL_TST_DBL_DQ]
[H1:CAL-DELTAL_CTRL_UIM_DBL_DQ]
[H1:CAL-DELTAL_RESIDUAL_DBL_DQ]
Initially H1:CAL-DELTAL_EXTERNAL was made double, but this meant the DTT display on the control room TV stopped working, this channel was then reverted to single precision.
ETMY Hardware Watchdog
Jim, Dave:
To verify the HWWD software is running, I pulled out the top monitor cable on the front panel of the EY HWWD unit. The LED fault light came on, showing that the software did start after the power outage.
DMT broadcaster channel configuration
Joe B, Maddie, Dave
Joe and Maddie requested some channel additions/removals from the DMT channel list, these were installed.
QUAD_MASTER change
Jeff
h1susetmx,y were restarted against a change to the QUAD_MASTER.mdl file. Next Tuesday QUAD_ITM_MASTER.mdl will be modified and the itm models restarted.
Resync of atomic clock to NTP GPS receiver
Jim:
Following the power outage, Jim re-synchronized the atomic clock to the 1PPS signal coming from the NTP GPS receiver. This provides an independent timing of the Timing Master.
The Timing Solutions time code generator in the MSR has been resync'ed to the Symmetricom NTP server following the power outage of June 4,5. The following steps were performed: Sync the Symmetricom 4310B Cesium Frequency Standard to the Symmetricom S350 NTP Server. Attach the 1PPS output of the Symmetricom S350 NTP server to the Sync input of the Symmetricom 4310B Frequency Standard (aka atomic clock). Attach matched coax cables from the Sync input of the 4310B and the 1PPS Out of the 4310B using coax T's to an oscilloscope to measure the difference between the two signals. Difference between sync pulse and 1PPS out was 178mS. Send sync arm command to the 4310B via the RS-232 input to the 4310B. Observe that the difference between the sync pulse and the 1PPS out of the 4310B is now in the range of 30nS. Note that the jitter of the S350 is much greater than the output of the 4310B, so precise measurement is not really possible to the 1nS range. Also, it's only possible to sync the 4310B to a 1PPS input to +/- 100nS. Sync the Timing Solutions timing distribution system to the 4310B. Next, attach the 1PPS out of the 4310B to the Timing Solutions TCG 1PPS Sync Input. Using the oscilloscope, measure the difference between the TCG 1PPS Sync Input and the TCG 1PPS Output. This difference was in excess of 430mS. Adjust the slew using the Timing Solutions MCA user interface to bring the difference as small as possible. This involved several iterations as it is only possible to adjust the slew by a maximum of 100mS for each adjustment. Also note that the minimum adjustment is 10nS. The ending best difference obtained was 4nS between the 4310B 1PPS output and the TCG 1PPS output. The 4310B 1PPS output is connected to the TCG Advance 1PPS input to monitor coincidence to +/- 20nS between the 4310B and TCG 1PPS. No adjustment of the Timing Solutions TCT units at the end station is needed, as the only adjustment possible is to compensate for the length of the optical fiber between the corner and end stations which is a fixed value. The 1PPS output of the Timing Solutions TDS is connected to the corner station timing comparator input 1 while the Symmetricom S350 NTP server is connected to input 2. The MEDM screen for the timing comparator shows the difference between the timing master 1PPS derived from a separate GPS receiver and the various inputs.
Subsytem Updates
REMINDERS:
Carl, Terra, Rich A.
1. We used the never-before-tested LVLN ITMY ESD driver to ring up and damp two mechanical modes of ITMY, 15072 Hz and 14979 Hz. Below is the amplitude evolution of the 15072 Hz peak as we rang it up, allowed it to ring down naturally, rang it up again, and then damped it down fully with a gain sign flip.
We drove and damped similarly for the 14979 Hz peak. 15072 Hz is the vertical differential drumhead mechanical mode; 14979 Hz is the horizontal version. For both cases, we tightly bandpassed the H1:OMC-PI_DCPD_64KZ_A/B signal, added a +60deg damping filter, and added gain to the damping filter until saturation. Positive gain drove up, negative damped down.
I've fit the natural ring down with f(x) = a * exp(bx), where tau = - (1/b). Then Q = pi * 15072 * tau = 31.5 million.
2. Interestingly, we realized after the above tests that we had not turned the ESD bias on for either ITM. After turning on both to 100K cts (to DC offset), we just had time to ring back up the ITMY 15072 Hz mode before a lockloss. Below is a comparison of the ring ups (note we lost lock and did not damp for the ring down portion below). Green trace is the 15072 Hz ring up without bias, blue trace is with bias (time shifted for ease of comparison).
A slight slope difference is visible but we'll look into this more. Thoughts are that we could use the difference in responses to measure test mass charge coefficients as discussed here and here.
3. We turned ring heaters off and on for future mode-mass identification analysis. For the record (since RH messing with violin modes was a concern tonight):
We were also driving the ITMX mode at 15077Hz before and after the bias was turned on. The data is a little more confusing. As the phase was being varied before hand to try find the optimum damping phase. In the plot the amplitude is made to coincide when the damping phase with no bias had the largest response. There was no attempt to optimise the phase in the case where the bias was on other that to try the positive and negative of the previous 'best' phase. Interesting points are:
For the 15077Hz mode the phase that excited the mode was flipped whent he bias was engaged.
The response with the bias engaged is relatively larger when compared to the ITMY 15072Hz mode.
For the 15077Hz mode the amplitude response with bias is about twice the amplitude response without bias.
Mode identification.
In the attached plot the relative change in frequency of the four likely drum head modes around 15200Hz is plotted as a function of time. The ring heaters were adjusted as follows:
ETMX 0.5W per segment to 0W 02:29
ETMY 0.5W per segment to 0W 03:00
ITMX 0W to 0.5W per segment 03:33
ITMX 0.5W per segment to 0W 04:02
ITMY 0W to 0.5W per segment 04:06
ITMY 0.5W per segment to 0W 04:36
The response in frequency shows that the 15218Hz mode is a cooling ETMX, 15219Hz is a cooling ETMY, 15197Hz a heating ITMX and 15192 a heating ITMY mode.
As we were only operating at 2W input the signal to noise ratio of modes is a lot lower and many fewer modes are visible compared to Livingston measurements T1600141.
The following is the list of modes identified.
Measured Frequencies | Simulated Frequency | ||||||||
ITMX | ITMY | ETMX | ETMY | ITMX | ITMY | ETMX | ETMY | Shape description | |
6044 | 6042 | 6055 | 6054 | 6057 | 6054 | 6053 | 6053 | Butterfly | |
8162 | 8160 | 8161 | 8158 | 8194 | 8190 | 8189 | 8188 | Drumhead | |
9812 | 9809 | 9829 | 9831 | 9827 | 9827 | 9832 | 9832 | Drumhead | Vertical |
9860 | 9863 | 9881 | 9882 | 9879 | 9881 | 9885 | 9885 | Drumhead | Horizontal |
10415 | 10437 | 10434 | 10433 | 10432 | 3 point ripple | ||||
10423 | 10419 | 10428 | 10426 | 10462 | 10460 | 10463 | 10462 | Drumhead | |
12992 | 12988 | 12999 | 13000 | 13022 | 13026 | 13035 | 13036 | Drumhead | Differential |
15077 | 15072 | 15093 | 15096 | 15102 | 15103 | Drumhead | Horizontal | ||
15197 | 15192 | 15218 | 15219 | 15227 | 15227 | 15232 | 15231 | Drumhead | |
15540 | 15545 | 15541 | 15544 | 15544 | Drumhead | Vertical | |||
15627 | 15635 | 15632 | 15634 | 15634 | Drumhead | Horizontal | |||
19555 | 19553 | 19589 | 19588 | 19595 | 19594 | Drumhead |
Stefan, Nutsinee
A very quick summary: ITMY modes rung up tonight. All the damp settings has been changed. I suspect ETMX damp phase might have been changed as well. ITMX didn't seen to have a problem. We suspect the ring heater changed the violin mode frequency.
Below I attached a screenshot of what works on ITMY modes tonight.
We looked at the frequency of mode 8 last night, two nights before around GW151226 and after a O1 power ouitage (20151025 9UTC). (attached)
Two days ago the mode was heavily damped, so the frequency estimate might be a bit off. At any rate, all frequencies are within 3mHz (or 22deg of feed-back phase in the current band-pass filter). So we ough to be able to damp this mode consistently.
Chandra, Gerardo -Adjusted CP5 Dewar pressure regulator 1/4 turn CW to raise Dewar pressure (nominally at 15 psi). Will wait hours/days to see results. -Measured Dewar vacuum jacket pressure = 20 microns, satisfactory. -Magnehelic readout is noisy - jumps from 35-40 in. of H20, even when CP5 is in manual mode. -Manually tested full range of stroke of electronic actuator - looks OK. -Tightened actuator support/brace nuts along unistrut. -Noticed needle of valve has wiggle (more than CP6). Verified packing nut is adjusted properly and coupling nut and needle are fully threaded. So now we wait to see if adjusting pressure of Dewar has an impact.
Not able to trend - was sitting around 100% but maybe just big swing?
Sheila, Haocun
We had lots of locklosses today, and some were probably caused by the side OSEM driving of the beamsplitter M1, as shown in the figures attached.
(H1:SUS-BS_M1_MASTER_OUT_SD_DQ & H1:SUS-BS_M1_DAMP_T_IN1_DQ)
It seems to have happened again at 22:20:12 UTC June 15th
Tagging SUS and CDS.
Tonight I had a look at several other of these locklosses, and other than the one at 23:57:59 that Hoacun plotted, it seems mostly like the glitch in the side osem comes after the lockloss. Many of these locklosses seem to happen durring the CM board switching.
It does seem like it would be worth investigating this side osem issue, since we know that it causes at least some locklosses.
New today is the the SEI_CONF Guardian node that will manage all of the other configuration nodes that we have made the past few weeks. This will hopefully make it much easier to choose a configuration for all of SEI. Sensor Correction for the entire SEI can be turned OFF from here, and then be brought back, of course.
The states and their names are very likely to change in the near future. So keep an eye out for updates, and we will also try to keep documentation up to date.
I also move the GUARD_OVERVIEW.adl around a bit since it was getting crowded in spots.
I changed the Blend_PAGE_MAIN.adl to include all of the config nodes and removed the SC ON/OFF switches. I have also included a guide to help choose states, but take it lightly because I'm sure it will change again soon. The right corner shows if the MATCH bank gains are 1 (green) to allow for sensor correction (Enabled) in that chamber.