J. Kissel The message: I tried to improve the Bounce mode damping, I've installed some new filters; no change in the amount of damping for the bounce modes, and they haven't been cooled any more than when I got started. I've made no changes the guardian code, so in the next lock stretch my work will be erased / cleaned up (except for the new filters), and what used to be turned on turned on before today will be turned on again. After spending an hour or so trying to get better results out of the installed bounce mode damping filters to no avail, I tried installing much more narrow band-pass filters. I did this because of Shiela's entry (LHO aLOG 18440) seemed to indicate we know the frequencies really well, and I could see what I thought was ETMX at 9.77 [Hz] interacting with what I thought was ITMX at 9.83 [Hz]. However, I found after exploring the parameter space (with gains and +/- 30 or 60 [deg] filters) with the more narrow filters, the situation got more confusing because the top mass of ETMY, ETMX, and ITMX would come in and out of showing 9.83 [Hz]. See for example the attached .pdf of the top mass vertical OSEMs for each. Yuck! After a while I convinced myself that having both IX and EX using the narrow band bp9.7-9.8 (in FM10 of IX and FM5 in EX) with a small gain of +0.1 and +0.05 and their normal -60 [deg] and +60 [deg] filters wasn't making it worse, but I couldn't get any improvement. This mode interaction are stable / non-existent at low power, but as soon as we elevate to 23 [W], the MICH oscillation creeps in mentioned earlier today (LHO aLOG 18478) creeps in and breaks the lock.
I'm leaving the ISC_LOCK requested to be at LSC_FF, so hopefully the IFO will return to it's automatic locking of ~20-30 minute, ~55 [Mpc] lock stretches. The guardian will turn off the undisturbed bit when it looses lock, so I won't bother setting it as I leave this time. Note that Robert's still here, and other than a trip to the Y-arm beam tubes between 12:30p and 1:30p PDT, he hasn't performed any injections (to my knowledge).
J. Kissel I've processed the DARM OLG TFs from Thursday (see 18443), and found that the DARM Coupled Cavity Pole Frequency seems relatively stable at 350 +/- 10 [Hz] with respect to IFO power. Admittedly, the data at 23 [W] is only over a small frequency range and is varying by ~5% (even though, as usual, I applied at 0.99% coherence threshold cut on the data points), so that's why I put a +/- 10 [Hz] on the number I quote. I think we'll either have to (a) get more DARM sweeps, (b) use Kiwamu's LSC DCCP tracker, and/or (b) get the PCAL demodulation up and running to check it this is consistent from lock-stretch to lock stretch. The ITMs are under control for all of these data sets, and the SRC1 yaw loop is tuned as described in LHO aLOG 18442 for the latter three. For all four measurements SRCL to DARM subtraction, ("FF") was off.
Just a quick report; following Jeff's analysis of the DARM open loops, I looked into the DARM cavity pole tracker signals from the successive 7.5 hr lock stretch (alog 18455).
Enjoy some interesting plots:
IFO back up; bounce mode looks better (maybe a slightly different alignment?); 2.5 [kHz] violin harmonic is still quite rung up. Turns out our DMT session had just frozen on the wall -- the IFO has been locking itself for the past 8 hours!! The lock stretches have only been ~1/2 hours long though, I'll see if its these poorly damped modes. Robert has joined the party, and he plans to do injections later in the morning, but for now I've checked the undisturbed bit at 2015 May 16 16:10:00 UTC.
Wow, that is awesome. Re-locked itself all the way from DOWN to LSC_FF 8 times overnight. Very very cool.
Lost lock at 16:23 UTC. Looking at some lock loss plots and ASAIR 90, the symptom seems to be a slow, 0.3-0.4 [Hz] oscillation in MICH that slowly builds up over time. Maybe this the is 0.4 [Hz], first longitudinal mode of the ITMs getting slowly rung up from radiation pressure? #wildspeculation
Stable lock attempts continue, I've been able to measure a high-resolution ASD, and identify that it's a competition between the ETMX and ITMX's bounce mode at 9.77 and 9.84 [Hz] being rung up the most. I'm trying to hold the ISC_LOCK guardian in DC readout (before the coil drivers get switched to low range) and noodle the with top mass bounce mode damping settings on these SUS to see if I can get better action. But it's not working out well at all. I'll keep exploring the parameter space... The remarkable thing is how consistent the lock-acquisition sequence is -- nice work commissioning vanguard!
J. Kissel I came in ~10 minutes ago, and found that the IFO had just finished getting to low noise, in the LSC_FF state. ITMX's highest vertical (bounce) mode and the 4th harmonic (2.5 [kHz]) of somebody's violin mode are very rung up, but the IFO appears to be stable(-ish). I briefly turned on the undisturbed bit, but when I noticed that DMT wasn't computing the inspiral range, I found there was a state above LSC_FF, "NOISE_TUNINGS" so I've turned off the undisturbed bit, changed the nominal state to NOISE_TUNINGS, and went for it (since all it appeared to do was turn on a SRCL offset). This broke the lock after a minute or so. Hurumph. I've now changed the nominal state back to LSC_FF, and I'm gunna leave it there is time. Once it gets back up (hopefully), I'll reengage the undisturbed bit.
J. Kissel, for B. Lantz and J. Harms There's been some interesting discussions on the DetChar mailing list about seasonal variation of ground motion around the world's network of gravitational wave observatories. Two things (of many) that have come out of it that might be of particular interest to LHO: - Jan Harm's data set from 2009/2010 T1500224, showing what I'm calling "meta" ASDs showing the percentile ASDs for each month over a year. - Brian Lantz has used that data, shown it in a slightly different way which can be found in SEI aLOG 766. More data is to come. Why do we care? It's merely to call attention to the fact that for the active isolation feed-back, (a) We shouldn't expect a "one (blend filter / sensor correction) filter fixes all" solution for both LHO and LLO. (b) Seasonal variations of the secondary microseismic may explain why we've needed to work harder at the microseism this past winter (i.e. use the 45 [mHz] DeRosa, LLO blend) and now we're doing just fine with a 90 [mHz] blend everywhere.
Dan, Keita, Stefan - Updated the OMC model with correct math and enough monitor points to be useful. - We tried commissioning the model sitting off fringe at 23W, but confused ourself a bit. Everything started to make sense when Keita pointed out that we have a fringe offset when sitting at zero digital offset, which manifests itself as different OMC trans DC power levels on the left and right side of the fringe. - With the nominal 1e-5 DARM ERR counts offset (about 14pm), we observed 18mAmp, while -1e-5 gives 15mAmp. Thus we have about 0.6pm fringe offset when the RF length signal is zeroed. - With that we laid out the strategy for setting up the parameters: 0) Set H1:OMC-READOUT_SIMPLE_GAIN to 0 (Turns off simple path.) 1) Estimate the minimum power at zero DARM offset to get a number for the contrast defect. It turns out this is close to 0, so 0 will do. 2) Pick the transition point. In our case 14pm. Thus set both XF (fringe length) and X0 (fringe offset) to 0. 3) Adjust H1:OMC-READOUT_PREF_OFFSET until H1:OMC-READOUT_STEP2 is equal to one on average. (H1:OMC-READOUT_STEP2 is the power normalized by the power at the transition point.) 4) Set H1:OMC-READOUT_X0OFFSET_GAIN to 1 now the .H1:OMC-READOUT_ERR_INMON signal is zero on average. 5) Use H1:OMC-READOUT_ERR_GAIN to match the gains in the RF sensor (H1:LSC-DARM_IN1) and the OMC path (H1:OMC-READOUT_ERR_OUT). 6) Set H1:LSC-OMC_DC_OFFSET to minus the H1:LSC-DARM_OFFSET. Now the offsets are matched. 7) Transition the input matrix. 8) Ramp H1:LSC-OMC_DC_OFFSET and H1:LSC-DARM_OFFSET to 0 (at the same time) - We had high winds for about two hours. But after that successfully tested the script that follows this logic.
Here's the MEDM screen Stefan built for the new handoff code.
The OMC guardian and the ISC_lOCK guardian have been updated to follow the above script. After the transition to DC readout, the DARM offset is offloaded to the "fringe offset" in the OMC-READOUT path.
We were able to transition to DC readout at 2.2W, 3e-5 DARM offset, and increase the power to 22W while decreasing the DARM offset to 1e-5. This was done with only a casual regard for ramp times; the process was very robust. We have coded a step in INCREASE_POWER that changes the fringe offset to maintain 20mA in DCPD sum. Eventually, we will juggle the states in the ISC_LOCK guardian so we transition to DC before the power up, and then smoothly adjust the fringe offset during the power up sequence. We tested this once and it works well, but we didn't have the energy to perform a shuffling of guardian states. For now, we have left the guardian as usual, but the handoff from the OMC to DARM is performed using the new normalization path. It worked once, before we left.
We never reached low noise tonight, we lost lock a few times in the coil driver switching step. At least one of these locklosses was due to what appeared to be a guardian error, in which the ISC_LOCK guardian skipped the COIL_DRIVERS state, went to the ETMY ESD handoff, and then ran the ETMY ESD handoff *again* before losing the lock. See the second screenshot attached.
Here's my interpretation of what happend with the ISC_LOCK guardian (see comments in the screen shot below):
This is all kosher (i.e. no bug) assuming there was indeed a second request for LOWNOISE_ESD_ETMY at 7:20. If there wasn't, and the second request truly came out of nowhere, then that's weird. Where could it have come from?
If the second request was legitimate, i.e. it came from the operator, it might have been accidental, or the operator was maybe un-aware of the intended behavior (not unlikely given that behavior is probably not widely known). However, it might also indicate that the behavior is non-intuitive and/or problematic and should maybe be re-examined.
The same-state redirect behavior was added early in commissioning the SUS guardians. (It was actually added in response to some of the problems we've had handling the DC alignment biases, which is actually slated to change (ECR 746).) I'm not sure how widely this behavior is know or taken advantage of. It's not actually being used in the SUS guardians anymore. Given the current crappy MEDM UI, it's also quite easy to accidentally re-request the same state.
My intention was that all REQUESTable states be idempotent, e.g. it would never hurt anything to simply re-run them again. If we were to make sure that all requestable states really were idempotent, then the behavior is probably fine as is. This is actually quite easy to accomplish in this case by simply making the relevant non-idempotent action state request=False
, and adding a new requestable idle state right after it.
If this is too cumbersome, the we need to do some or all of the following:
Over the lock period last night, the CPS ASD appears to experience changes as time passes.
Just to focus on the BS, the ST2 Vertical between 1 and 5 hz see the largest variation. See the two attached images, the lower right panel is the ST2 Vertical. The lighter thicker colors are references before the CPS timing change on Tuesday. The finer darker traces are current running. I've look at the ITMY & ITMX too over the few hours of the lock loast night, and at some time or another, all the traces are at or below the reference, suggesting the CPS are okay.
Brute Force coherence report is available here for last night lock:
https://ldas-jobs.ligo.caltech.edu/~gabriele.vajente/bruco_1115716006/
Some highlights, which I find very interesting:
Not much is happening at higher frequency, at least for broad band coherence. However there are a lot of single frequency bins with significant coherence, and I can't list all of them here (I'm too lazy). So look into the table at your favorite line frequency, you might be rewarded with some coherence.
After the time that Gabriele looked at, the coherence of CARM with DARM went up to nearly 1 above 500 Hz (and 0.5 at 100 Hz). This seemed to continue for the rest of the lock. All of the loudest glitches, which had peak frequencies at about 200 Hz, showed up in CARM as well as DARM. What is going on? The first plot is a reproduction of Gabriele's coherence plot. The second is about a half hour later, and the third near the end of the lock. The last is a coherence spectrogram of CARM with DARM during this very loud glitch near the end of the lock, which showed up as strong in CARM as DARM.
LSC-CARM is not used anymore for any feedback to the IFO; we use it sometimes to make an audio channel out of DARM so that we can listen to the IFO. Stefan was tuning the audio filters after we hit the intent bit last night.
The actual feedback to MC2 from the CARM slow loop goes through the LSC-MCL filter bank. If you want an out of loop CARM sensor you can look at LSC-REFL_A_9I, for the error and control signals you can look at the common mode board channels, eg LSC-REFL_SERVO_ERR_DQ or CTRL_DQ
By the way, instead of stealing the CARM filter bank to process audio, you can apply the filter directly in the audio player:
$ cdsutils audio H1:CHANNEL-NAME -f 'filterspec'
where 'filterspec' is a filter design string, as you would type into Foton. (This was requested in CDS bug 797 and added in cdsutils r444.)
It seems unneccesarily confusing to use a named channel with a very clear symantic meaning for something other than it's intended purpose. Can we not just make some separate filter banks just for filtering audio output? Either that or use the built-in cdsutils audio filtering that Chris mentioned.
I spent a bit of time this morning clearing out remaining SDF diffs in SUS, LSC, and ASC. This involved the usual alog hunt followed by some clarifications from commissioners. Things still showing diffs at the moment and why:
- We need to capture a new snap file for the OMC due to yesterday's model change which will clear out it's large SDF diff list.
- The ASCIMC still has one diff due to a 10e-17 size change in one value. I think this feature is due for an upgrade at some point, so until then, this 1 diff will be lingering.
- Same with ISCEX - there are 2 diffs with very small diffs which won't go away with the accept button. Will linger until update.
- We need some assistance with clearing the PSL diffs from team PSL.
- I've not paid much attention to the bottom monitors, namely ODC and CAL.
Comment posted to the wrong aLOG. Deleted!
Sheila, Evan, Stefan - Reduced the whitening gain on both AS 36 ASC diodes. The new value is 24dB. We increase the input matrix values from 1 to 2 to compensate (for MICH and SRC1, i.e. the SRM loop). - We tried to phase the AS_A_36 quadrants by wiggling SRM in yaw and pitch, maximizing the I signal. Not sure how meaningful that was though, because the actual signal we ended up using was Q. - We switched the the sensor for SRC1_YAW (SRM) to 0.3*AS_A_36_Q + 0.3*AS_B_36I. This signal clearly reacted to our mysterious SRC cavity drifts, and had zero offset at the good locking point. - We increased the power to 23W. This is the maximum power currently available. With this new SRM yaw sensor the 23W was very stable. - The DARM offset was 14pm (or 1e-5 cts). - This gave us a recycling gain of 38. - Measured the open loop gain and verified the calibration. - This configuration was stable for 1h. - We took SRCL noise injections at: 8:46:15UTC and 8:59:40UTC (SRCL cut-off filter off and on) - We took MICH noise injections at: 8:52:15UTC and 8:55:40UTC (SRCL cut-off filter on and off) - Pushed the intent bit at 9:06:30 UTC. - The inspiral range seems quite remarkably flat at 57 Mpc. A stron indication that we are now purley electronics noise limited at low frequencies. To do tomorrow: - Explore DARM offset: the new OMC code is ready to install, and will allow on-the-fly OMC offset tuning. - Add the new ASC INMATRIX to Guardian. - Add the new DARM offset to Guardian. - Generate noise budget.
Attached is a new spectrum, before and after screenshots of the AS A 36 WFS phases, and 3 OLG measurements made at 3 different input powers.
Great progress, congratulations!!
I will take a closer look at the noise budget over the weekend, but it seems (roughly) that above 20 Hz we are currently limited by DAC→ESD noise and quantum noise. As usual, there is some uncertainty in coupling strength of the ESD to DARM, which may explain the discrepancy in the plot from 20 Hz to 100 Hz.
This is fantastic -- perfect timing for the low noise Electro-Static Driver installation next week!
Excellent news! Well done.
Looks like discharging and lower noise driver will do some good. Nice going.
Truly great progress. Nice work, LHO commissioning team!
Beautiful! Nice and stable as well! I was not able to change my slides for the invited LIGO talk at APS Northwest Section Meeting, but I did announce your success. Great timing too.
Well done all at LHO. That's great progress.
Bounce | Roll | ||
ETMY | 9.730 Hz | 13.816 Hz | 14854 |
ETMX | 9.77 | ||
ITMY | 9.8135 | 13.934Hz | 18395 |
ITMX | 9.8469 | 15400 |
the remaining roll modes are at 13.89 and 13.98 Hz, but we don't know which is from ETMX and which is from ITMX
The ITMY bounce mode is 9.83 Hz.
This is a follow up entry of LHO ALOG 17601.
A couple of days ago, the discrepancy of the response for DCPDA and DCPDB were found. This was basically caused by misadjusted filter modules for the anti-whitening filters. Some of them were using design values (like Z10:P1) and some others were just left as they had been imported from the LLO setup.
In order to correctly take the whitening transfer functions into account, the wiring of the in-vacuum and in-air connections were necessary to be tracked down. The 1st attachment shows the sufficiently detailed wiring chain for this task. Using the test data (links indicated in the diagram), we can reconstruct what the correct anti-whitening filters should be. The summary can be found below.
[Trivia for Rich: DCPD1 (Transmission side of the OMC BS) is connected to HEAD2, and DCPD2 (Relfection side of the OMC BS) is connected to HEAD1. This is because of twisted D1300369. This cable has J2 for HEAD2 and J3 for HEAD1. This twist exists in LLO and LHO consistently, as far as I know]
=======
Characteristics of the DCPD electronics chain
Complex poles/zeros are expressed by f0 and Q
DCPD A
(DCPD at the transmission side of the OMC DCPD BS)
- Preamp D060572 SN005
Transimpedance: Z_LO = 100.2, Z_HI = 400.0
Voltage amplification ZPK: zeros: 7.094, 7.094, (204.44 k, 0.426), poles: 73.131, 83.167, 13.71k, 17.80k, gain: 1.984
- Whitening filter D1002559 S1101603
(This document defines the gain not at the DC but at a high frequency. The gain below is defined as a DC gain.)
CH5 Whitening
Filter 1: zero 0.87, pole 10.07, DC gain 10.36/(10.07/0.87)
Filter 2: zero 0.88, pole 10.15, DC gain 10.36/(10.15/0.88)
Filter 3: zero 0.88, pole 10.20, DC gain 10.36/(10.20/0.88)
Gain: “0dB”: -0.051dB (nominal), “3dB”: 2.944dB, “6dB”: 5.963dB, “12dB”: 11.84dB, “24dB”: 24.04dB
DCPD B
(DCPD at the reflection side of the OMC DCPD BS)
- Preamp D060572 SN004
Transimpedance: Z_LO = 100.8, Z_HI = 400.9
Voltage amplification ZPK: zeros: 7.689, 7.689, (203.90 k, 0.429), poles: 78.912, 90.642, 13.69k, 17.80k, gain: 1.983
- Whitening filter D1002559 S1101603
CH6 Whitening
Filter 1: zero 0.88, pole 10.13, DC gain 10.41/(10.13/0.88)
Filter 2: zero 0.87, pole 9.96, DC gain 10.40/( 9.96/0.87)
Filter 3: zero 0.88, pole 10.15, DC gain 10.41/(10.15/0.88)
Gain: “0dB”: -0.012dB (nominal), “3dB”: 2.982dB, “6dB”: 6.007dB, “12dB”: 11.87dB, “24dB”: 24.04dB
=======
Now we put these transfer functions into the model and check if we can reproduce the observed relative difference (Attachment 2). In deed, the measurement is well explained by the model below 30Hz where the measurement S/N was good. As we saw in the previous entry, the difference of the DCPDA and DCPDB after the whtening compensation is 20% max. Note that further inspection revealed that this 20% difference is, in fact, mostly coming from the difference of the preamp transfer functions rather than the miscompensation.
So this was the relative calibration between DCPDA and DCPDB. How is the compensation performance of each one? The 3rd attachment shows how much of current we get at the output as H1:OMC-DCPD_A_OUT, H1:OMC-DCPD_B_OUT, and H1:OMC-DCPD_SUM_OUT, if we give 1mA of photocurrent to DCPD_A, DCPD_B, or both (half and half). Ideally, this should be the unity. The plot shows how they have not been adjusted. For the our main GW channel we take sum of two DCPDs. The individual deviations were averaged and thus the sum channel has max 10% deviation from the ideal compensation. This shows up in the GW channel.
=======
So let’s implement correct compensation. Basically we can place the inverse filter of the each filters. The preamplifier, however, includes some poles and zeros whose frequency are higher than the nyquist frequency. Here we just ignore them and assess how the impact is.
The result is shown as the 4th attachment. Upto 1kHz, the gain error is less than 1%. This increases to 5% above 3kHz. The phase error is 7deg at 1kHz. This increases to 20deg above 3kHz. These are the effect of the ignored pole/zeros. Note that these are static error. In fact, the phase error is quite linear to the frequency. Thus this behaves as a time delay of ~18.5us. Since the phase delay at 100Hz is small, the impact to the DARM feedback servo is minimal. For the feedforward subtraction, however, this might cause some limitation of the subtraction performance. In practice, we measure the coupling transfer function in order to adjust the subtraction, in any case. Therefore this delay would not be a serious problem.
The filter bank to implement the new compensation was already configured. The filter file is attached as foton_DCPDfilters.txt
.
Once we lock the full IFO, we measure the DARM OLTF and give it to Kiwamu for recalibration.
With the new filters, the balance is extremely good now.
This indirectly suggests that the individual compensations are done pretty well.
J. Kissel Since the front-end calibration did not account for this whitening compensation mis-match, i.e. it assumed perfect compensation, the calibration of the sensing function was simply *wrong* (inaccurate) at these frequencies were there was a mis-match. (Recall the DARM UGF is ~40 [Hz], so the mismatch began influencing the calibration only above ~40 [Hz]) As such, now that the whitening and preamps have been more accurately compensated the calibration as it stands has now simply become *more correct*. Therefore we will not need to change or correct anything in the front end calibration filters. Stay tuned for further study.
Jeff -- don't be so hasty. The absolute DC gain of the sensing function (or the inverse sensing function in the CAL CS model) is set by scaling an open loop gain TF measurement to a model. Thus far, open loop gain TFs have only been taken between ~10 and ~100 [Hz], exactly where this discrepancy occurs. Thus, the IFO's DC sensing function is likely off in overall scale factor by the ~10-20% caused by this discrepancy. So, once we get the IFO back up, we'll take another open loop gain transfer function, compare it against the prior, determine a new DC gain for optical gain / sensing function, and update the calibration accordingly.
At the section "Characteristics of the DCPD electronics chain", I wrote something inconsistent with the other part of the entry.
DCPD A is the DCPD at the reflection side of the OMC DCPD BS
DCPD B is the DCPD at the transmission side of the OMC DCPD BS
My hand written cartoon is correct.
I wish I could correct the aLOG entry that is older than 24 hours.