Daniel, Evan
We looked at the DARM spectrum in the vicinity of the OMC dither line (at 4100 Hz) and its second harmonic (at 8200 Hz) using the OMC DCPD IOP channels.
The spectrum around the line itself appears to be bilinear, perhaps as one would expect (if the OMC is exactly on resonance, there should be no first-order modulation of the transmitted power). The second harmonic of the line is quite clean, with no obvious upconversion. Is this what we should expect?
That's what you expect (when OMC transmission is quadratic to the length).
Power propto L^2 = (L0+dL+dither)^2 = L0^2 + 2L0(dL+dither) + dL^2+ 2dL*dither + dither^2
where L0 and dL are the DC- and AC-ish component of the OMC length and dither is the dither.
dL*dither^2 will only appear in L^3 response, and (dL*dither)^2 in L^4.
Summary:
I was baffled to find that each and every omega scan I have seen so far shows some wandering lines in H1:IMC-PZT_PIT_OUT_DQ and YAW that extends to hundreds of Hz, e.g.
https://ldas-jobs.ligo-wa.caltech.edu/~nairwita/wdq/H1_1126294545/#H1:IMC
These wandaring lines seem to be irrelevant of any of the glitches, and they're there even when there is no fast feedback to the PZT (indeed, since Sept. 09 the fast feedback has been disabled https://alog.ligo-wa.caltech.edu/aLOG/index.php?callRep=21300).
This is an artefact of single precision rounding error that happens when the data is saved on frame files (in the frontend the data is double).
That tells us that we need to be careful when looking at signals with big DC and small AC.
Analysis 1: Measured long data spectrum VS predicted single precision rounding error spectrum agree very well.
When you look at the first attachment, you'll see that PZT_OUT and PZT_IN1 are different now (blue and green) though the only difference should be that the former has a fixed offset added in. Also, the ratio of the noise level above 20Hz is approximately equal to the ratio of the DC value.
Though the rounding error spectrum depends on the wave form before the casting, if you take a long (relative to the frequency analyzed) data and make an average, like I did here, you can expect that it becomes white-ish. Since the RMS of the signals are totally dominated by the DC, and since I know that the exponents did not change during the measurement, you can calculate the expected white-ish noise spectrum as
2^(floor(log2(DC level))) / 2^23 / sqrt(12) / sqrt(1024),
i.e. 1.1E-6 /sqrtHz and 2.2E-6 for PZT PIT and YAW OUT respectively, which perfectly agree with the measured noise floor. 1/2^23 is a "1-bit" fraction, and 1/sqrt(12) because the RMS of +-0.5 bit fraction is 1/sqrt(12) bits when you assume a uniform distribution, and sqrt(1024) because this is first down sampled into 2kHz and then cast into single.
Analysis 2: Short data spectrum shows wandering structures
Since what we really want to know is short term variability at around 100Hz, I looked at the 6Hz BW spectrum.
In the video attachment, on the top is IMC-DOF_3_Y_IN1_DQ signal, which is the error point of one of ASC sensors and that's the only signal that goes to PZT YAW. On the bottom is the IMC-PZT_YAW_OUT_DQ.
Red are short data spectrum (2 averages exponential), blue are the long spectrum (0.047HZ, 10 averages fixed). This is from around the same GPS time as the example omega scan time.
Look how variable the bottom is, you can easily see the same thing as Omega scan!
This is not the case with the top plot. This is because the DOF3 signal is NOT limited by the single precision casting due to small DC value relative to the noise floor. DOF3 signal is zero on average, and RMS is almost always smaller than 0.05, so the expected rounding error noise floor is in the vicinity of 10^-10 to 10^-11/sqrtHz.
Analysis 3: How large is the "true" noise going to PZTs?
Third attachment shows DOF_3_PIT_IN1_DQ and YAW (these are the only source of PZT feedback after Sept. 09) at the example omega scan time, projected on the PZT output using the DOF3 filter. This can be regarded as the "true" noise as the high frequency noise of DOF3 is not limited by the rounding error.
Anyway, as you can see, PZT_OUT spectrum is at least 5 orders of magnitude larger than the "true" noise level at 100Hz, and there's nothing to worry about.
I just realized that I have had two LSC test point channels open since the beginning of this lock acquisition and in Observation Mode. I just turned them off. 20:52 UTC Stopped running H1:LSC-TR_X_NORM_OUT and H1:LSC-TX_Y_NORM_OUT test points in dataviewer. Checked that they are no longer present on the H1LSC_GDS_TP medm.
Whoa. How do we know if we have test points up? Heck, I might have been looking at these channels, too (or did I have non-test point channel versions of these channels?). Is it bad to look at Test Points? (if it is, is there a monitor to show us test points are being looked at?)
This is ok as long as we don't open way too many.
ER8 Day 32. No restarts reported.
ER8 Day 31. No restarts reported.
15:41 UTC Joe D. called from VPW. Turning on a couple of water pumps and a water tap. 16:32 UTC lockloss, earthquake, either 5.4 in Northern Mid-Atlantic Ridge (15:59:42 UTC) or 4.5 63km W. of Ovalle, Chile (15:58:37 UTC) 16:36 UTC Reset OMC SUS WD ALS is not staying locked, still riding out earthquake 16:46 UTC Kyle and Gerardo leaving control room to start replacement of ion pump at mid Y (WP 5500) 16:52 UTC Attempting to lock DRMI 16:56 UTC Richard reports he turned on AC for LDAS room in staging building, it had not been on 16:58 UTC Locked DRMI, lost DRMI lock, but guardian did not go to down, continued to attempt to lock DRMI 17:00 UTC Kyle and Gerardo at mid Y 17:16 UTC Searching for ALS DIFF IR resonance by hand 17:18 UTC Done searching for IR resonance by hand, the peak went by on the strip tool, but guardian says it is moving on. I guess it caught it. 17:32 UTC Back to Low Noise. Seismic 0.03-0.1 Hz band is around 10^-1 um/s. 17:38 UTC Back to Observing. Range is low ~ 54 MPc. 17:43 UTC Kyle reports seeing a helicopter flying low over mid Y about 5 minutes ago. Range seems to be rising as seismic comes down. 18:52 UTC Kyle and Gerardo back. Ion pump has been replaced. Aux cart attached to new ion pump left running. They left mid Y at 18:40 UTC.
A reminder to those on site. Please do not touch thermostats without talking to John, Bubba, Richard or the occupant of a room. This morning the cleaning crew called to let me know the new DCS room in the warehouse was warm. I investigated and the room was 88 degF and not in cooling mode. This was not a good state. I have turned the unit to cooling and it was immediately bringing the temperature down. Luckily we do not have a large heat load yet. Also on Wednesday when doing a safety tour we found the warehouse VPW fan set to auto instead of on. This is a problem as we loose pressure in the room. Just because we do not have "DO NOT TOUCH" signs on everything, does not mean it should be fiddled with.
DCS switched Disk2Disk copying of raw and commissioning frames, and rd and aggregated hoft generation from ER8 directories to O1 directories this morning before 8 am.
We also have 2 extra Disk2Disk processes running to archive commissioning frames from Sept 1 to the end of ER8 and for the beginning of O1
TITLE: 9/18 OWL Shift: 7:00-15:00UTC (00:00-8:00PDT), all times posted in UTC
STATE OF H1: In Observation Mode with a range which was trending down.
SUPPORT: Most of the shift solo.
SHIFT SUMMARY: Other than the 2hrs of downtime due to low frequency seismic, the shift went fairly smoothly, and was nice to take H1 into O1 in Observation Mode!
Incoming DAY Operator: Patrick T.
DAY'S ACTIVITIES:
I was curious what caused the lock loss at 12:04 (1126613091.1). There is an oscillation that starts in ASC-AS_B_RF36_I 0.5s before lockloss, and SRC2 and the SRM start to wiggle at ~4Hz. An 8Hz oscillation starts in DHArd and CHard 0.25s later, and ITMY followed by ITMX L2 stages saturate just before the lockloss.
TITLE: 09/18 [DAY Shift]: 15:00-23:00 UTC (08:00-16:00 PDT), all times posted in UTC STATE Of H1: Locked at Low Noise, Observing Mode, range around 68 MPc. OUTGOING OPERATOR: Corey G. QUICK SUMMARY: Lights are off in the LVEA, PSL enclosure, end X, end Y and mid X. I can not tell from the camera if they are off at mid Y. Seismic in 0.03 - 0.1 Hz band has recovered from a peak to ~1 um/s about 5 hours ago, now around 10^-2 um/s. Seismic in 0.1 - 0.3 Hz band has been trending slowly down over the last 18 hours, is now below 10^-1 um/s. Winds are less than 20 mph, may be starting to come up. Range has been trending down since the beginning of the lock about an hour and a half ago. Started around 75 MPc, now around 68 MPc. O1 HAS BEGUN!
~10:50 Left to cook salmon and bagel/cream cheese, and set ISC_LOCK to DRMI
11:22 Came back to H1 locked on DRMI (since 11:18)
NOTE: So DRMI finally locked up, and on the 0.03-0.1 seismic band, the velocity was down to 0.1um/s. Recap:
Observation Mode was set to Earthquake from 10:03-11:34 (should be more like 9:30-11:18, where 11:18-11:33 would be Lock Acquisition).
11:35 Back to Observation Mode (with a range trending toward 80Mpc).
ITM ISI GS13 SDFs ACCEPTED with HI gain.
TITLE: 9/18 OWL Shift: 7:00-15:00UTC (00:00-8:00PDT), all times posted in UTC
STATE OF H1: Smooth sailing at a steady ~74Mpc. O1 begins in T-Minus 6hrs!
OUTGOING OPERATOR: Ed M.
SUPPORT: Kissel hung out for an hour, but flying solo now.
QUICK SUMMARY: H1 has been locked at NOMINAL_LOW_NOISE since 9/17 at 13:49UTC (that's ~19hrs) & the current Observation Segment is going on 15.5hrs. Range is averaging ~74Mpc.
Seismically, we look decent. Terramon is warning of an impending 5.1 Chilean quake (0.27um/s) at 9:05utc. Low winds.
Ed has left me with a task list for the next lockloss:
TITLE: Sep 17 EVE Shift 23:00-07:00UTC (16:00-00:00 PDT), all times posted in UTC
STATE Of H1: Observing
LOCK DURATION: Entire Shift
SUPPORT: Jeff K
INCOMING OPERATOR: Corey
Activity log:
23:00 IFO locked and In Observing mode. 70Mpc
00:03 Reset a timing glitch at ETMX
03:23 ETMY Saturation
03:45 ETMY Saturation
04:45 ETMY Saturation
04:46 ETMY Saturation
05:30 ETMY Saturation
End-of-Shift Summary: IFO Locked entire shift 72-76Mpc. No obtrusive Seismic or wind activity. 5 ETMY glitches. Approx. 12 hours coincidence with LLO.No Lockloss, therefore the ITM GS13s are still in LO GAIN state. Handing of to Corey.
J. Kissel Though I have not yet flipped the sign on te actuation filter as I suspect I need to (see LHO aLOG 21627, and more to come), I've heard on several calls today that in order to resume, hardware injections need to go though an offline test showing that the wave-form must be passed through each sites actuation filter offline to confirm it will not cause saturations in the IFO. In order to facilitate that study, I attach the H1CALCS.txt foton file currently running. Recall that you want the produce FMs 3 and 4 in the H1:CAL-INJ_HARDWARE bank (see LHO aLOG 21487).
C. Cahillane Today I have fit phase delays on the various actuation stage residuals (L1, L2, L3) and a zero pair on the magnitude residual of the L1 actuation stage to appease the calibration gods. The L2 and L3 actuation magnitude weighted mean residuals were already flat, so I did not have to fit them. This gives us a good idea of what our systematics look like for actuation. Slight phase delays on all at 100 Hz, and large deviation for the L1 stage magnitude. I will use these results to "correct" systematic error in our model by making a puesdo-model with the systematics included and comparing the two. To do: 1) Get realistic values for kappa uncertainty 2) Apply systematic "corrections" to a pseudo-model and compare with actual model 3) Make correct A_pu calculation corrections from Mathematica notebook 4) LLO
Summary:
Not much to report. IFO locked and Observing - 74MPC. Wind < 10mph only a mild rise in earthquake graph. ~ 10.5 hours of coincidental locking with LLO. 5 ETMY glitches.
This is a budget for the DARM noise on 2015-09-12. Compared to the previous budget, the only major difference (that we know of) is the elimination of the excess 45 MHz AM.
For the DARM estimation I use the DCPD sum, which is calibrated into milliamps. Calibration into freerunning displacement requires the DARM OLTF (I used the CalSVN measurement from 2015-09-10) and an estimate of the optical plant. I use Kiwamu's measurement of the optical plant along with Sudarshan's corrections to arrive at an optical gain of 3.39 mA/pm, and a DARM pole of 349 Hz.
Quantum noise is taken from GWINC, using the same parameters as previously: recycling gain 37 W/W, 107 kW circulating power, 755 kW BS power, 87 % quantum efficiency, 14 % other readout losses. The dark noise is taken from a measurement made several months ago.
The null stream is higher than this shot noise prediction by about 7 %. This was the case previously, but this is almost certainly a coincidence, since the previous budget used data with two stages of whitening on (which we already know has a small calibration discrepancy). Previously I had used pcal directly to estimate the optical plant (without the necessary pcal correction factors), so this could explain part of the issue.
The attachment shows sum, null, and quantum noise curves, along with 20.0 mA shot noise. It could be that we are missing some readout loss, or the GWINC curve could still require some more tuning of the power.
The seismic and Newtonian curves are the vanilla ones from GWINC, as is the coating Brownian noise.
DAC noise (using Peter's model) is propagated forward to displacement noise at the quads. The elevated noise from Chris's pringle measurements is not yet included.
Intensity noise and frequency noise are as described previously. Intensity noise was measured with the ISS outer loop off. I injected into the ISS inner loop error point and used the outer loop's out-of-loop PD array (sum 5–8) to estimate the RIN. Frequency noise is a lower bound, including only the noise of the first stages of the CARM loop (PD dark noise, PD shot noise, and electronics noise of the CARM board).
LSC control noises include PRCL, MICH, and SRCL. ASC control noises include dHard, BS, and SR2 loops (these are our high-bandwidth loops). I injected broadband noise for each of these and estimated the coupling via ratio of excess power (not coherent TFs).
Compared to the last budget, the SRCL coupling has increased slightly (see attached injection spectra), which is responsible for the slightly increased LSC noise around 50 Hz. Unclear whether this was a short-term excursion or a long-term drift of the SRCL coupling.
Squeeze film damping and residual gas are calculated as before. End station pressures are assumed to be 1×10−8 torr of molecular hydrogen. For residual gas, I have used 5×10−9 torr of molecular hydrogen.
IMC input jtter was measured by driving the IMC PZT in pitch and in yaw. On Keita's suggestion, as a sensor for the PZT motion I use WFS B I for pitch and WFS A I for yaw. Again the coupling is estimated by excess power ratio. Jitter at the OMC is not yet incorporated.
Based on Bruco, there is still some coherence with the endstation QPDs around 78 Hz, as Gabriele has already noted. But other than that, no channel seems to jump out as the culprit for the noise from 50 Hz to 150 Hz. So either we are looking at some nonlinear process, or we are seeing some kind of displacement noise that is not captured by the digital system.
In the above entry I have added a zoomed plot in the region where we have unexplained noise.
In this entry I have attached the residual of the measured noise and the budgeted noise, in two versions. The first uses the total noise trace as given in the above entry. The second uses a total noise trace in which the quantum noise has been increased by 7 %, since (as mentioned above) the quantum noise curve seems to underestimate the shot noise by this amount.
[Of course, simply scaling the overall quantum noise is not correct, but since the radiation pressure part of the curve is buried below the total noise, it doesn't matter too much here.]
Injection times, 2015-09-12 Z: