40m QIL Cryo_Lab CTN SUS_Lab TCS_Lab OMC_Lab CRIME_Lab FEA ENG_Labs OptContFac Mariner WBEEShop
  Cryo Lab eLog  Not logged in ELOG logo
Entry  Fri Apr 30 12:07:25 2021, aaron, DailyProgress, Laser, marconi locking 252D9563-0122-418B-9B97-29D44FEF83E9.jpeg
    Reply  Fri Apr 30 15:15:33 2021, rana, DailyProgress, Laser, marconi locking 
       Reply  Mon May 3 23:28:56 2021, rana, DailyProgress, Laser, delay line frequency discriminator DFD.jpg
          Reply  Wed May 5 11:54:19 2021, aaron, DailyProgress, Laser, delay line frequency discriminator 9x
             Reply  Wed May 5 22:21:17 2021, rana, DailyProgress, General, delay line frequency discriminator 
             Reply  Thu May 6 10:48:40 2021, aaron, DailyProgress, Laser, delay line frequency discriminator 3EFFFA27-AAF6-434F-A1AB-FE7C657A464C.jpegDL_3CH.pdf
                Reply  Tue May 11 14:37:48 2021, aaron, DailyProgress, Laser, delay line frequency discriminator, current noise measurement 
                   Reply  Wed May 12 15:35:09 2021, aaron, DailyProgress, Laser, delay line in a box 69A8E6A2-7F85-464E-BDAA-8C6BBD2910B5.jpeg1D05A31E-7DBC-448B-A2F4-A8805EFAC393.jpeg25C2DD1E-7776-491E-BBCE-5F2D1EA935E3.jpeg
                   Reply  Wed May 12 16:52:08 2021, aaron, DailyProgress, Laser, delay line frequency discriminator 3878D102-C5BB-46AE-B90B-540CFDE3100B.jpegSPSR785_NX_12-05-2021_173959.pdfSPSR785_SX_12-05-2021_183434.pdfSPSR785_NS_12-05-2021_191513.pdf
                      Reply  Thu May 13 16:31:40 2021, aaron, DailyProgress, Laser, delay line frequency discriminator 6x
                         Reply  Thu May 13 19:09:16 2021, rana, DailyProgress, Laser, delay line frequency discriminator 
                         Reply  Fri May 14 18:43:25 2021, aaron, DailyProgress, Laser, delay line frequency discriminator 994ECC7F-C5FE-4E7B-8F19-B1C118DCE6CE.png
                            Reply  Thu May 20 12:50:56 2021, aaron, DailyProgress, Laser, delay line frequency discriminator 
                               Reply  Fri May 21 14:15:21 2021, aaron, DailyProgress, Laser, current noise of custom current drivers 5012D147-8B13-478C-BB83-57A71A2B304A.jpeg467DEC59-01B1-40CD-9F4E-4D9D87183A78.jpegABC2A51B-6726-4F9C-AD88-551ED7C90EEF.jpeg
                               Reply  Tue Aug 17 14:04:27 2021, aaron, DailyProgress, Laser, delay line frequency discriminator 3A1E242C-F4F6-4A87-B2E8-E74839BE542C.jpeg
                                  Reply  Wed Aug 18 09:55:13 2021, aaron, DailyProgress, Laser, delay line frequency discriminator 0A96102F-5465-456B-93BF-E2B9D4A7E96C.jpegtrim.2553D656-98A0-442D-9185-4C5252B2B727.MOVtrim.787F04FC-8DED-4501-B67D-FF98EC5D5D56.MOV
                                     Reply  Thu Aug 19 14:21:35 2021, aaron, DailyProgress, Laser, delay line frequency discriminator 007FC2D7-F29B-4636-9CE7-04CBC3B789F7.jpeg1F1E68F0-199F-4E1C-B909-883A521E06CF.jpegScreenshot_from_2021-08-19_17-53-58.png
                                        Reply  Fri Aug 20 11:03:44 2021, aaron, DailyProgress, Laser, delay line frequency discriminator 30A43669-5DFF-400A-BEB4-7EF71BBAE1A6.jpegA3BF5DA5-32E5-42A5-A1D8-9270C3D88518.pngScreenshot_from_2021-08-20_16-11-26.pngScreenshot_from_2021-08-20_16-20-48.png
                                           Reply  Mon Aug 23 13:38:57 2021, aaron, DailyProgress, Laser, delay line frequency discriminator Screenshot_from_2021-08-23_16-42-48.png
                                              Reply  Tue Aug 24 16:44:39 2021, aaron, DailyProgress, Laser, delay line frequency discriminator beats_regions.pdfTeraxion_noise.pdfEW_timeseries.pdf
                                                 Reply  Thu Aug 26 11:58:50 2021, aaron, DailyProgress, Laser, delay line frequency discriminator 9x
                                                    Reply  Fri Aug 27 15:27:38 2021, rana, DailyProgress, Laser, delay line frequency discriminator 
                                                    Reply  Tue Aug 31 15:14:28 2021, aaron, DailyProgress, Laser, quieting the Rio W and Teraxion beat notes, but not much help beats_regions_base2.pdfbeats_regions_base2.pdfTeraxion_noise_base2.jpgDLFD_foam.pdf
                                                       Reply  Wed Sep 1 12:20:18 2021, rana, DailyProgress, Laser, quieting the Rio W and Teraxion beat notes, but not much help 
                                     Reply  Thu Aug 19 16:32:47 2021, rana, DailyProgress, Laser, delay line frequency discriminator 
       Reply  Tue May 4 10:50:22 2021, aaron, DailyProgress, Laser, marconi locking 6DC5CB12-BCDF-4550-AFB8-42D3DC354AED.jpegSPSR785_04-05-2021_152326.pdf
          Reply  Wed May 5 00:27:10 2021, rana, DailyProgress, Laser, marconi locking 
Message ID: 2802     Entry time: Thu Aug 26 11:58:50 2021     In reply to: 2801     Reply to this: 2804   2807
Author: aaron 
Type: DailyProgress 
Category: Laser 
Subject: delay line frequency discriminator 

Yesterday (August 25), I measured another 10 min at 488 Hz of Moku phasemeter data for the three pairs of beat notes. All three lasers had been on overnight, so there was no longer low frequency drift of the Teraxion laser. Rather than amplifying then picking off the beat note, I sent the RF output of the 1611 directly to the Moku's phasemeter input.

what kind of spectral density estimate to use?

Today, I've been figuring out how to get more averages out of our data. One approach (the one used above) is a modification of Welch's method:

  1. Start Welch's method with as large a window as possible, given the desired number of averages. For example, if there are 2**16 data points and we want at least 4 averages, then use Welch's method with 50% overlap and 2**15 points per segment. 
  2. From the first application of Welch's method, save the first N frequency bins (including the DC value). Our spectral estimate runs from 0 to (N-1)*f_0 in steps of size f_0.
  3. Next, decrease the number of samples per segment by a factor of N, and repeat Welch's method. In our example above, we are now able to take 4*N averages. Assuming the rounding is handled correctly, the frequency resolution of the second Welch's periodogram is f_1=N * f_0. 
  4. From the second estimate, save the N-1 samples from f_1 to (N-1)*f_1.
  5. Continue applying Welch's method to successively smaller window sizes and saving the 'low frequency' data from each. The iteration terminates when the window size is unnacceptably small (for example, nperseg <= 2), at which point you can save the remaining spectrum up to the Nyquist frequency.

The above procedure sacrifices some frequency resolution at the higher frequencies in exchange for additional averaging. The tradeoff with resolution is necessary, because the window size determines not only the smallest resolvable frequency, but also the spacing of frequencies in the spectrum. For the spectra from the previous elog in this thread where N=10, the total measurement time is 10 minutes, and the sampling rate 488 Hz, there is evidently more noise higher in the frequency decade (7-9e^n) than lower (1-3e^n). More consistent averaging can be achieved by setting N=2, but at the expense of most of the high frequency resolution (only 33 frequency bins survive the procedure). 

One workaround is to modify the procedure so the frequency binning is mostly set at the beginning, by the 'highest resolution' available. Then, perform Welch's method with as small a window as possible while still resolving the frequency bins. Care must be taken at high frequency: eventually, the 'FSR' of our Welch's method cannot resolve an f_0 difference in frequency into an integer change in the number of samples per segment. The spectrum can either be cut off at that frequency, or the procedure can continue while accepting nonstandard bin widths at high frequency.

Other workarounds are perhaps less desirable. One could accept nonuniform frequency binning, and simply compute Welch's method for every available choice of nperseg. This would maximize the number of averages in each bin, but especially at low frequency, there will be substantial correlation between adjacent frequency bins. Another workaround is to save the entire spectrum wherever evaluated, then combine the data later. One must again worry about correlations between the measurements: at high frequency, we would be combining coarse data with many averages with fine data with fewer averages. 

Another approach entirely is to do something smarter than Welch's method. In our meeting today, Chris suggested I look into multitapering. Spectral estimates can reduce bias due to leakage by introducing a tapered window, at the cost of increased measurement variance. Welch's method heals the variance relative to standard tapering by overlapping the windowed segments, at the cost of some frequency resolution. Multitapering instead minimizes loss of information by increasing the number of degrees of freedom of the estimates. The Here are a few resources on the topic:

While checking out Percival and Walden, I stumbled across parametric methods for spectral estimation -- those where an early spectral estimate is used to refine the procedure and spectral estimate iteratively. Perhaps up the alley of some recent discussions at our group meeting.

error propagation

Simply applying Gaussian error propagation is not quite right, because the PSD is exponential distributed (the ASD is Rayleigh distributed, see Evan's note T1500300). Each ASD is Rayleigh distributed with

  • mode \sigma
  • mean \left(\frac{\pi}{2}\right)^{1/2}\sigma \equiv \mu
  • rms \sqrt{2}\sigma = \frac{2}{\sqrt{\pi}}\mu 
  • median \sqrt{\ln 4}\sigma
  • variance of \frac{4-\pi}{2}\sigma^2 \equiv \delta(X_{ij})^2

For a large number N of averages, the central limit theorem lets us estimate the mean of each PSD, x_{ij}^2, with normal distributed uncertainty and variance \delta(X_{ij}^2)^2/N. Our three corner hat estimates of the ASD Y are based on the scaled, root mean-squared sum of three such PSD estimates, so for each frequency bin we can estimate the variance of the laser's ASD by

Y = \frac{X_{ab}^2 + X_{ac}^2 - X_{bc}^2}{2} \\ \delta(X_{ij}^2)=2x_{ij}\delta(X_{ij}) \\ \delta(X_{ab}^2 + X_{ac}^2 - X_{bc}^2) = \sqrt{\delta(X_{ab}^2)^2+\delta(X_{ac}^2)^2 + \delta(X_{bc}^2)^2} \\ \implies \delta(Y)^2 = x_{ab}^2\delta(X_{ab})^2+x_{bc}^2\delta(X_{bc})^2+x_{ac}^2\delta(X_{ac})^2\\

I'll use the final equation above, along with the number of averages, to estimate the uncertainty in each frequency bin of the final frequency noise ASD of the individual lasers. In particular, the filled region is \sqrt{Y \pm \delta(Y)}

Attachments

[OK, I'm having a lot of trouble uploading pdfs to the elog this week, even with rasterizing. I've dropped these figures along with one set for the case of 'Welch's with no averaging' onto gaston under /home/controls/cryo_lab/Figures/3CH ]

  1. Time series of Rio E x Rio W beat frequency
  2. "" Rio E x Teraxion ""
  3. "" Rio W x Teraxion ""
  4. "" Marconi ""
  5. Welch's estimate of the frequency noise on the beat notes, with increased averages every decade
  6. Three corner hat estimate of the frequency noise on the Teraxion laser, from the ASD in (5)
  7. Welch's estimate of the frequency noise ont he beat notes, with increased averages every factor of 4 in frequency
  8. Three corner hat estimate of the frequency noise on the Teraxion laser, from the ASD in (7)
  9. Welch's estimate of the frequency noise on the beat notes, with averages increasing every factor of 2 in frequency
  10. Three corner hat estimate of the frequency noise on the Teraxion laser, from the ASD in (9)

I think these results warrant a more careful measurement, especially in the decade around 1 Hz. Also, the error bars are obviously way underestimated. 

Attachment 1: EW_timeseries.pdf  118 kB  Uploaded Fri Aug 27 18:53:25 2021  | Hide | Hide all
EW_timeseries.pdf
Attachment 2: EX_timeseries.pdf  114 kB  Uploaded Fri Aug 27 18:53:42 2021  | Hide | Hide all
EX_timeseries.pdf
Attachment 3: WX_timeseries.pdf  102 kB  Uploaded Fri Aug 27 18:53:55 2021  | Hide | Hide all
WX_timeseries.pdf
Attachment 4: Ma_timeseries.pdf  95 kB  Uploaded Fri Aug 27 18:54:07 2021  | Hide | Hide all
Ma_timeseries.pdf
Attachment 5: beats_regions_base10.pdf  29 kB  Uploaded Sat Aug 28 21:19:32 2021  | Hide | Hide all
beats_regions_base10.pdf
Attachment 6: Teraxion_noise_base10.jpg  595 kB  Uploaded Sat Aug 28 21:19:43 2021  | Hide | Hide all
Teraxion_noise_base10.jpg
Attachment 7: beats_regions_base4.pdf  24 kB  Uploaded Sat Aug 28 21:19:52 2021  | Hide | Hide all
beats_regions_base4.pdf
Attachment 8: beats_regions_base2.pdf  22 kB  Uploaded Sat Aug 28 21:20:08 2021  | Hide | Hide all
beats_regions_base2.pdf
Attachment 9: Teraxion_noise_base2.jpg  535 kB  Uploaded Sat Aug 28 21:20:25 2021  | Hide | Hide all
Teraxion_noise_base2.jpg
ELOG V3.1.3-