Hello, I am Varun Kelkar. I will be working at the 40m lab as a SURF student this summer with Eric Quintero on Audio processing for real time control system signals. This week I will mostly be working on implementing basic DSP C-code offline. Currently I am trying to write a code for noise whitening.
Varun has received 40m specific basic safety training today.
Local EQ 3.5 mag at 2:28 UTC May 24, 2016 Rancho Cucamonga, Ca.....no damage
WFS locking point seemed degraded; I hand aligned and reset the WFS offsets as usual.
ITMX oplev recentered. While doing so, I noticed an ETMX excursion rear its head for the first time in a long while :
There was no active length control on ETMX, only OSEM damping + oplevs. Afterwards, its still moving around with only local damping on. I'm leaving the oplevs off for now.
I think you should use the current actual PRC & SRC cavity lengths as measured, as it would be simplest to simply replace the folding mirror optics without changing the macroscopic lengths / optic positions. (EDIT: Gautam rightly points out that we have to move things around regardless, since our current lengths include propagation through the folding mirror subtrates)
Moreover, the recycling cavity lengths you posted are not the right "ideal" lengths to use, as they do not account for the complex reflectivities of the sidebands off of the arm cavities (I have made this mistake myself). See this 40m wiki page for details.
In short, given our current modulation frequency, the ideal lengths to use would be:
These are the lengths that the recycling cavity optics were positioned for (though we did not achieve them perfectly). If you do a finer PRC/SRC length scan around the DRFPMI resonance of your model, you would presumably see some undesired sideband splitting.
Having played around with a toy finesse model, I went about setting up a model in which the RC folding mirrors are not flipped. I then repeated the low-level tests detailed in the earlier elog, after which I ran a few spatial mode overlap analyses, the results of which are presented here. It remains to do a stability analysis.
Overview of model parameters (more details to follow):
Results (general note: positive RoC in these plots mean a concave surface as seen by the beam):
Next step is to carry out a stability analysis...
Contacted Charles regarding use of Altium. Got to know that Altium is installed on cit40m iMac in Win7 on VirtualBox. Had to update Virtualbox to get it working. Altium now works for sometime, but then fails, saying that it is unlicensed.
I did a quick measurement of the ITMX oplev loops, both pitch and yaw have about the same upper UGF as previous measurements with the previous laser; about 4 Hz.
Tested the AGC today with LSC cavity transmission signal and error signal. Not in real time still.
Key to attachments:
cav_tr-eps-converted-to.pdf: LSC cavity transmission signal input
cav_tr_out-eps-converted-to.pdf: LSC cavity transmission signal, output of the AGC.
2005 ALL oplev servos use Coherent DIODE LASERS # 31-0425-000, 670 nm, 1 mW
Sep. 28, 2006 optical lever noise budget with DC readout in 40m, LIGO- T060234-00-R, Reinecke & Rana
May 22, 2007 BS, SRM & PRM He Ne 1103P takes over from diode
May 29, 2007 low RIN He Ne JDSU 1103P selected, 5 purchased sn: T8078254, T8078256, T8078257, T8078258 & T8077178 in Sep. 2007
Nov 30, 2007 Uniphase 1103P divergence measured
Nov. 30, 2007 ETMX old Uniphase 1103P from 2002 dies: .............., running time not known......~3-5 years?
May 19, 2008 ETMY old Uniphase 1103P from 1999 dies;.....................running time not known.....~ ?
Oct. 2, 2008 ITMX & ITMY are still diodes, meaning others are converted to 1103P earlier
JDSU 1103P were replaced as follows:
May 11, 2011 ETMX replaced, life time 1,258 days or 3.4 years
May 13, 2014 ETMX , LT 1,098 days or 3 y
May 22, 2012 ETMY, LT 1,464 days or 4 y
Oct. 5, 2011 BS & PRM, LT 4 years, laser in place at 1,037 days or 2.8 y
Sep. 13, 2011 ITMY old 1103P & SRM diode laser replaced by 1125P ..........old He life time is not known, 1125P in place 1,059 days or 2.9 y
June 26, 2013 ITMX 622 days or 1.7 y note: we changed because of beam quality.........................laser in place 420 days or 1.2 y
Sep. 27, 2013 purchased 3 JDSU 1103P lasers, sn: P893516, P893518, P893519 ......2 spares ( also 2 spares of 1125P of 5 mW & larger body )
May 13, 2014 ETMX, .............laser in place 90 d
May 22, 2012 ETMY,
Oct. 7, 2013 ETMY, LT 503 d or 1.4 y............bad beam quality ?
Aug. 8, 2014 ETMY, .............laser in place 425 days or 1.2 y
Sept. 5, 2014 new 1103P, sn P893516 installed at SP table for aLIGO oplev use qualification
May 23, 2016 ITMX dead laser sn P845648 replaced after 1062 days [2.9 yrs] by 1103P, sn P859884, with output output 2.6 mW, nicely round beam quality at 15 meters.
May 23, 2016 ITMX dead He/Ne laser sn P845648 replaced after 1062 days [2.9 yrs] by 1103P, sn P859884, with output 2.6 mW, nicely round beam quality at 15 meters.
Power just before viewport 1 mW, returning light on qpd 154 microW = 7,500 counts
As we realized during the EX table switch, the transmitted beam height from the arm is not exactly 4" relative to the endtable, it is more like 4.75" at the X-end (yet to be investigated at the Y-end). As a result, the present configuration involves the steering optics immediately before the Oplev and TransMon QPDs sending the beam downwards at about 5 degrees. Although this isn't an extremely large angle, we would like to have things more level. For this purpose, Steve has ordered some Aluminium I-beams (1/2 " thick) which we can cut to size as we require. The idea is to have the QPD enclosures mounted on these beams and then clamped to the table. One concern was electrical isolation - but Steve thinks Delrin washers between the QPD enclosure and the mount will suffice. We will move ahead with getting these machined once I investigate the situation at the Y end as well.. The I beams should be here sometime next week...
I have written a basic version of AGC, and have done some tests with a data file. will do tests on whitening and agc today. Also, today I have to go to the SSN office. Hence will be late.
Finished writing the code on whitening. I have to still test it. uploaded on github noise cancellation repo. @eric could you give me some data of noise power spectral density for testing the code?
I've been working on putting together a Finesse model for the current 40m configuration. The idea was to see if I could reproduce a model that is in agreement with what we have been seeing during the recent DRFPMI locks. With Antonio and EricQs help, I've been making slow progress in my forays into Finesse and pyKat. Here is a summary of what I have so far.
Having put together the .kat file (code attached, but this is probably useless, the new model with RC folding mirrors the right way will be what is relevant), I was able to recover a power recycling gain of ~7.5. The arm transmission at full lock also matches the expected value (125*80uW ~ 10mW) based on a recent measurement I did while putting the X endtable together. I also tuned the arm losses to see (qualitatively) that the power recycling gain tracked this curve by Yutaro. EricQ suggested I do a few more checks:
Conclusion: It doesn't look like I've done anything crazy. So unless anyone thinks there are any further checks I should do on this "toy" model, I will start putting together the "correct" model - using RC folding mirrors that are oriented the right way, and using the "ideal" RC cavity lengths as detailed on this wiki page. The plan of action then is
Sidenote to self: It would be nice to consolidate the most recent cavity length measurements in one place sometime...
We have good RGA scan now. There was no scan for 3 months.
3-4 hrs ago we run out of nitrogen. We are back to Vacuum Normal
It looks like the hardware reset did the trick. Previously, I had just tried ssh-ing into c0rga and rebooting it. At the time, however, Steve and I noticed that the various LEDs on the RGA unit weren't on, as they are supposed to be in the nominal operating state. Today, Steve reported that all LEDs except the RS232 one were on today, so I just tried following the steps in this elog again, looks like things are back up and running. I'm attaching a plot of the scan generated using plotrgascan MATLAB script, it looks comparable to the plot in elog 11697, which if I remember right, was acceptable.
Unless there is some reason we want to keep this c0rga machine, I will recommission one of the spare Raspberry Pis lying around to interface with the RGA scanner when I get the time...
Our last RGA scan is from February 14, 2016 We had a power outage on the 15th
Gautom has not succeded reseting it. The old c0rga computer looks dead. Q may resurrect it, if he can?
The c0rga computer was off, I turned it on via front panel button. After running RGAset.py, RGAlogger.py seems to run. However, there are error messages in the output of the plotrgascan MATLAB script; evidiently there are some negative/bogus values in the output.
I'll look into it more tomorrow.
This is a cold scan.
ALSX noise is solidly within past acceptable performance levels. The DRFPMI was locked on four out of six attempts.
Some housekeeping was done:
The recombination of the QPD signals to common / differential is imperfect, and limited how well we could keep the interferometer aligned, since the QPD at X has changed. This needs some daytime work.
Some sensing matrix measurements were made, to be meditated upon for how to 1F the DRMI.
As an aside, Gautam and I noticed numerous green beams coming from inside the vacuum system onto the PSL table. They exist only when green is locked to the arms. Some of them come out at very non-level angles and shine in many places. This doesn't make me feel very happy; I suppose we've been living with it for some time.
Indeed. This is why the LSC PDs have a 2f notch in addition to the 1f resonance. In recent versions, we also put a 2f notch in the feedback of the preamp which comes after the diode but before the mixer. The overall 1f to 2f ratio that we get is in the 50-60 dB region. I don't think we have to go that far with this thing; having a double LC already seems like it should be pretty good, or we could have a single LC bandpass with a 2f notch all in one Pomona box.
As I was looking at filter designs, it seemed difficult to get 40dB of supression at 2F with a bandpass without going to a pretty high order, which would mean a fair number of lossy inductors.
I'll keep working on it. Maybe we don't need 40dB...
Seems weird to design a PD lowpass with a corner at the modulation frequency. Recall what our strategy is with the other photodetectors we use for PDH servos: bandpass, not low-pass, and the band has to be wide enough to not effect the phase of the servo.
I've build the filter, and it seems to have the desired TF shape.
I also re-purposed the 70k lowass to a ~120k lowpass by changing the 68nF caps to 22nF caps, since we still want some post-mixer rolloff.
However, putting the ELPF in the chain caused some weird shapes in the OLG. I still need to get to the bottom of it. However, just with the post-mixer LPF modification, here's what the OLG looks like:
As Rana surmises, we definitely still add a boost and maintain a 10k UGF. I still need to look into the state of the remote boost....
I suggested in an earlier elog that after the repair of the NPRO, the PZT capacitance may have changed dramatically. This seems unlikely - I measured the PZT capacitance with the BK Precision LCR meter and found it to be 2.62 nF, which is in excellent agreement with the numbers from elogs 3640 and 4354 - but this makes me wonder how the old setup ever worked. If the PZT capacitance were indeed that value, then for the Pomona box design in elog 4354, and assuming the PM at ~216kHz which was the old modulation frequency was ~30rad/V as suggested by the data in this elog, we would have had a modulation depth of 0.75 if the Function Generator were set to output a Signal at 2Vpp (2Vpp * 0.5 * 0.05 * 30rad/V = 1.5rad pp)! Am I missing something here?
Instead of using an attenuator, we could instead change the capacitor in the pomona box from 47pF mica to 5pF mica to realize a modulation depth of ~0.2 at the new modulation frequency of 231.25 kHz. In any case, as elog 4354 suggests, the phase introduced by this high-pass filter is non-zero at the modulation frequency, so we may also want to install an all-pass filter which will allow us to control the demodulation phase. This should be easy enough to implement with an Op27 and passive components we have in hand...
All seems very fishy. Its not good to put attenuators and filters in nilly-willy.
This morning I poked around with the green layout a bit. I found that the iris immediately preceding the viewport was clipping the ingoing green beam too much, opening it up allowed for better coupling to the arm. I also tweaked the positions of the mode matching lenses and did some alignment, and have since been able to achieve GTRX values of around 0.5.
I also removed the 20db attenuator after the mixer, and turned the servo gain way down and was able to lock easily. I then adjusted the gain while measuring the CLG, and set it where the maximum gain peaking was 6dB, which worked out to be a UGF of around 8kHz. On the input monitor, the PDH horn-to-horn voltage going into the VGA is 2.44V, which shouldn't saturate the G=4 preamp stage of the AD8336, which seems ok.
The ALS sensitivity is now approaching the good nominal state:
There remains some things to be done, including comprehensive dumping of all beams at the end table (especially the reflections off of the viewport) and the new filters to replace the current post-mixer LPF, but things look pretty good.
We took an OLG measurement of the green PDH loop. It seems consistent with past measurements. I've added a trace for the the post-mixer lowpass, to show its contribution to the phase loss. (EDIT: updated with measured LPF TF)
I used this measured OLG and the datasheet laser PZT conversion factor to calibrate the control signal monitor into the AUX laser frequency noise, it looks consistent with the frequency noise measured via the PSL PLL (300 Hz/rtHz @ 100Hz). Above a few tens of kHz, the control signal measurement is all analyzer noise floor, due to the fourth order 70kHz lowpass after the mixer (the peaks change height significantly depending on the analyzer input range, so I don't think they're on the laser). Gautam will follow up with more detailed measurements of both the error and control signals as he noisebudgets, this was just intended as a quick consistency check.
Today we spent some time looking into the PDH situation at the X end. A summary of our findings.
With Steve's help, I installed the Oplev earlier today. I adjusted the positions of the two lenses until I deemed the spot size on the QPD satisfactory by eye. As a quick check, I verified using the DTT template that the UGF is ~5Hz for both pitch and yaw. There is ~300uW of power incident on the QPD (out of ~2mW from the HeNe). In terms of ADC counts, this is ~13,000 counts which is about what we had prior to taking the endtable apart. There are a couple of spots from reflections off the black glass plate in the vacuum chamber, but in general, I think the overall setup is acceptable.
This completes the bulk of the optical layout. The only bits remaining are to couple the IR into the fiber and to install a power monitoring PD. Pictures to follow shortly.
Now that the layout is complete, it remains to optimize various things. My immediate plan is to do the following:
I will also need to upload the layout drawing to reflect the layout finally implemented.
Not directly related:
The ETMx oplev servo is now on. I then wanted to see if I could lock both arms to IR. I've managed to do this successfully - BUT I think there is something wrong with the X arm dither alignment servo. By manually tweaking the alignment sliders on the IFOalign MEDM screen, I can get the IR transmission up to ~0.95. But when I run the dither, it drives the transmission back down to ~0.6, where it plateaus. I will need to investigate further.
GV Edit: There was some confusion while aligning the Oplev input beam as to how the wedge of the ETM is oriented. We believe the wedge is horizontal, but its orientation (i.e. thicker side on the right or left?) was still ambiguous. I've made a roughly-to-scale sketch (attachment #1) of what I think is the correct orientation - which turns out to be in the opposite sense of the schematic pinned up in the office area.. Does this make sense? Is there some schematic/drawing where the wedge orientation is explicitly indicated? My search of the elog/wiki did not yield any..
Antonio/Gautam are now developing a more up to date Finesse model of our recycling cavities to see what we can have there before our power recycling gain or cavity geometric stability is compromised. Expect that we will here a progress report on the model on Wednesday.
Steve sent 4 of our 1" diameter G&H HR mirrors to Josh Smith at Fullerton for scatter testing. Attached photo is our total stock before sending.
We can get as much, if not more, attenuation of the 1F line in the mixer output that we get from the post-mixer LPF from using the following passive filter between the PD and mixer RF input:
There should still be some kind of LPF after the mixer, but I haven't yet determined what it should be; this will determine how much phase the PDH loop wins. At most, this should win around 25 degrees at 10kHz.
The filter was designed by referencing the "Handbook of Filter Synthesis" by Zverev, looking for an elliptic filter for matched source and load impedences, 40dB min attenuation in the stopband, a stopband frequency that starts at twice the corner frequency, and minimizing the VSWR between the PD and filter in the passband.
In terms of the tables in the book, this means: n=5, rho=2%, theta=30deg, K**2 = 1.0. The dimensionless component values were scaled by the corner frequency of 200kHz, and reference impedence of 50 Ohm. (The corner is a little lower than the real modulation frequency, since the nonzero resistance of the inductors pushes the frequency up a bit)
The ideal capactior values do not correspond to things we have in hand, so I checked our stock and chose the closest value to each one.Unsurprisingly, due to these component substitutions, and the fact that the coilcraft inductors have a resistance of about 7 Ohms, the predicted TF of the realizable filter does not match the design filter exactly. However, the predicition still looks like it will meet the requirement of 40dB of supression of the 2F line in the PD signal. (Since we have tunable inductors, I've used the ideal inductor values in generating the TF. In practice I'll inspect the TF while I tune them)
[In this TF plot, I've multiplied the real response by 2 to account for the voltage division that occurs with ideal 50 Ohm impedance matching, to make 0dB the reference for proper matching]
The filter's phase delay at the modulation frequency is just about 180, which as a time delay of 5usec works out to 9 degrees of phase loss at 10kHz in the PDH loop. According to some old measurements, the current LPF costs something like 35 degrees at 10k, so this wins at most around 25 degrees, depedent on what LPF we put after the mixer.
LISO source both traces is attached!
After a second round of F.C. application, I think the window is clean enough and there are no residual F.C. pieces anywhere near the central parts of the window (indeed I think we got most of it off). So I am going to go ahead and install the Oplev.
It looks very promising.
Using the modulation frequency suggested here, I hooked up the PDH setup at the X-end and succeeded in locking the green to the X arm. I then rotated the HWP after the green Faraday to maximize TRX output, which after a cursory alignment optimization is ~0.2 (I believe we were used to seeing ~0.3 before the end laser went wonky). Obviously much optimization/characterization remains to be done. But for tonight, I am closing the PSL and EX laser shutters and applying first contact to the window once more courtesy more PEEK from Koji's lab in W Bridge. Once this is taken care of, I can install the Oplev tomorrow, and then set about optimizing various things in a systematic way.. MC autolocker has also been disabled...
Side note: for the IR Transmon QPD, we'd like a post that is ~0.75" taller given the difference in beam height from the arm cavity and on the endtable. I will put together a drawing for Steve tomorrow..
OK - but give us a circuit diagram and the expected before/after loop plots. Got to make sure we keep the right impedance from PD to mixer. Some of the Thorlabs PDs have a 50 Ohm instead of 0 Ohm source impedance. Maybe you can try it out now since the green arm is ready.
The 2F product out of the mixer is a natural concern when demodulating. However, I think this isn't so big of a deal in our green PDH servos; 420kHz isn't so high of a frequency that the servo amplifiers are bandwidth or slew-rate limited. Furthermore, the amplitude of this line is supressed by the loop somewhat, since it arises from the same field product that the loop is acting on. Measuring the Y end mixer output with a high impedance probe and the AG4395 shows it to be something like -50dBm.
In fact, the main thing that the pomona LPFs are accomplishing right now is filtering the 1F content of the mixer output that arises from the second order sideband creating a signal at 2F, and beating with the LO at (2F-1F)=1F. This line is something like -30dBm (5mVrms) at the mixer output; I can reproduce this amplitude with a back-of-the envelope calculation using a modulation depth of 0.3, 8V out of the PD at DC when unlocked, the mixer datasheet, and the nominal cavity parameters.
The nice thing about this is that we don't need to filter this after the mixer, we can use a [bandpass/lowpass/notch] filter before the mixer (as is done in the LSC demod boards) to filter out the 2F (420kHz) content of the PD signal, which will only introduce some small amount of linear time delay to the PDH loop, instead of the wicked phase loss from the current post-mixer LPF. We can then replace that 70kHz filter with something of lower order or higher corner frequency to win a good deal of phase in the PDH loop.
Building: Campus (see attached Map)
Date: Manhole 1 - May 3 through May 5
Manhole 2 – May 6 through May 10
Manhole 2 - May 16 through May 19
Manhole 3 – May 11 through May 19
Time: Noise: 7:00 a.m. To 5:00 p.m.
Access: 24 Hours a day
Interruption: Noise/Vehicular & Pedestrian Access
Storm Drain Manholes
*In order to repair 3 manholes associated with a large storm drain that runs north-south through the campus, work will take place at the
3 manholes shown on the map. This work will interrupt vehicular and pedestrian access on the paths adjacent to the manholes. Though the work at Manholes 1 and 2 will allow vehicular and pedestrian access around the manholes, the work at Manhole 3 will completely block the driveway running south from the southeast corner of Parking Lot 11. Noise will also be created by the repair
The IR Transmon system is almost completely laid out, only the QPD remains to be installed. Some notes:
I feel like once the above are resolved, the next step would be to PDH lock the green to the arm and see what sort of transmission we get on the PSL table. It may be the polarization or just alignment, but for some reason, the transmitted green light from the X arm is showing up at GTRY now (up to 0.5, which is the level we are used to when the Y arm has green locked!). So a rough plan of action:
Bah, we need ruby slippers for all future suspensions. Prism with curved backside and smooth grooves.
No aluminum, no cry.
Earth quake stops need viton tips.
Wirestandoffs are still aluminum.
The anti-symmetric port
spider webs fly in the wind
Layout as of today. Most of the green path is done. The Green REFL PD + PZT mirrors have not been hooked up to their respective power sources yet (I wonder if it's okay to start laying cables through the feedthroughs on either end of the table already, or if we want to put whatever it is that makes it airtight eventually in first?). A rough power budget has been included (with no harmonic separator just before the window), though some optimization can be done once the table is completely repopulated.
A zoomed-in version of the REFL path.
Some general notes:
I am closing the PSL shutter and the EX laser shutters for the night as I have applied a layer of first contact to the window for cleaning purposes, and we don't want any laser light incident on it. It may be that the window is so dirty that we may need multiple F.C. cleaning rounds, we will see how the window looks tomorrow...
Gautom is progressing with the layout nicely. The X-arm transmission window have not seen cleaning for decades. This should be the time to do it. Here is picture of dirtiness.
It is not that simple... How much effort should we put in it? The hole table with 1W inno laser plus... set up now about ~500 lbs We can pull it off carefully, but it is not risk free.
We should look at our other signal port windows! Gautom's long reach able him to do the first contact cleaning without moving anything. It is great!
Dan sealed the leak today.
Steve pointed out that in the aftermath of the Nitrogen running out a couple of times last week, the RGA had shut itself off thinking that there was a leak and so it was not performing the scheduled scans once a day. So the data files from the scheduled scans were empty in the /opt/rtcds/caltech/c1/scripts/RGA/logs directory. The wiki page for getting it up and running again is up-to-date, but the script RGAset.py did not exist on the c0rga machine, which the RGA is communicating with via serial port. I copied over the script RGAset.py from rossa to c0rga and ran the script on that machine - but the error flags it returned were not all 0 (indicating some error according to the manual) - so I edited the script to send just the initialize command ('IN0') and commented out the other commands, after which I got error flags which were all 0. After this, I ran a manual scan using 'RGAlogger.py', and it appears that the RGA is now able to take scans again - I'm attaching a plot of the scan results. We've saved this scan as a reference to compare against after a few days.
I've made progress on the new layout up to the doubling oven. After doing the coarse alignment with the diode current to the NPRO at ~1A, I turned it back up to the nominal 2A. I then rotated the HWP before the IR Faraday such that only ~470mW of IR power is going into the doubler (the rest is being dumped on razor beam dumps). After tuning the alignment of the IR into the doubling oven using the steering mirror + 4 axis translation stage on which the doubling oven is mounted, I get ~3.2mW of green after the harmonic separator and a HR mirror for green. The mode looks pretty good to the eye (see attachment #1), and the conversion efficiency is ~1.45%/W - which is somewhat less than the expected 2%/W but in the ballpark. It may be that some fine tweaking of the alignment + polarization while monitoring the green power can improve the situation a little bit (I think it may go up to ~4mW, which would be pretty close to 2%/W conversion efficiency). The harmonic separator also seems to be reflecting quite a bit of green light along with IR (see attachment #2) - so I'm not sure how much of a correction that introduces to the conversion efficiency.
While doing the alignment, I noticed that some amount of IR light is actually transmitted through the HR mirrors. With ~500mW of incident light at ~45 degrees, this transmitted light amounts to ~2mW. Turns out that this is also polarization dependant (see attachment #3) - for S polarized light, as at the first two steering mirrors after the NPRO, there is no transmitted light, while for P-polarized light, which is what we want for the doubling crystal, the amount transmitted is ~0.5%. The point is, I think the measured levels are consistent with the CVI datasheet. We just have to take care find all these stray beams and dump them.
I will try and optimize the amount of green power we can get out of the doubler a little more (but anyway 3mW should still be plenty for ALS). Once I'm happy with that, I will proceed with laying out the optics for mode-matching the green to the arm.