40m QIL Cryo_Lab CTN SUS_Lab TCS_Lab OMC_Lab CRIME_Lab FEA ENG_Labs OptContFac Mariner WBEEShop
  40m Log, Page 81 of 335  Not logged in ELOG logo
ID Date Author Type Category Subjectup
  16499   Fri Dec 10 15:59:23 2021 PacoUpdateBHDFinished Coil driver (even serial number) units tests

[Paco, Anchal]

We have completed modifications and testing of the HAM Coil driver D1100687 units with serial numbers listed below. The DCC tree reflects these changes and tests (Run/Acq modes transfer functions).

SERIAL # TEST result
S2100608 PASS
S2100610 PASS
S2100612 PASS
S2100614 PASS
S2100616 PASS
S2100618 PASS
S2100620 PASS
S2100622 PASS
S2100624 PASS
S2100626 PASS
S2100628 PASS
S2100630 PASS
S2100632 PASS
S2101648** FAIL (Ch1, Ch3 run mode)
S2101650** FAIL (Ch3 run mode)
S2101652** PASS
S2101654** PASS

** A fix had to be done on the DC power supply for these. The units' regulated power boards were not connected to the raw DC power, so the cabling had to be modified accordingly (see Attachment #1)

Attachment 1: dc_fail.jpg
dc_fail.jpg
  16514   Thu Dec 16 15:32:59 2021 AnchalUpdateBHDFinished Coil driver (odd serial number) units tests

We have completed modifications and testing of the HAM Coil driver D1100687 units with serial numbers listed below. The DCC tree reflects these changes and tests (Run/Acq modes transfer functions).

SERIAL # TEST result
S2100609 PASS
S2100611 PASS
S2100613 PASS
S2100615 PASS
S2100617 PASS
S2100619 FAIL (CH2 phase)
S2100621 PASS
S2100623 PASS
S2100625 PASS
S2100627 PASS
S2100629 PASS
S2100631 PASS
S2100633 Waiting for more components
S2101649** PASS
S2101651** PASS
S2101653** PASS
S2101655** PASS

** A fix had to be done on the DC power supply for these. The units' regulated power boards were not connected to the raw DC power, so the cabling had to be modified accordingly.

Further, Paco fixed the two even serial number units (S2101648, S211650) that failed the test. The issues were minor soldering mistakes that were easily resolved.

  16517   Thu Dec 16 17:57:17 2021 AnchalUpdateBHDFinished Coil driver (odd serial number) units tests

S2100619 was fixed by Koji and it passed the test after that.

Quote:
SERIAL #  
S2100619 FAIL (CH2 phase)

 

  12329   Mon Jul 25 10:54:55 2016 PrafulUpdateComputer Scripts / ProgramsFinished MEDM Tab on Summary Pages

The MEDM screen capture tab is now working and up on the summary pages: https://nodus.ligo.caltech.edu:30889/detcharsummary/day/20160725/medm/

Please let me know if you have any suggestions or notice any issues.

  12330   Mon Jul 25 12:22:18 2016 ranaUpdateComputer Scripts / ProgramsFinished MEDM Tab on Summary Pages

Looks pretty great. However, there's two problems:

1) Some of the MEDM screens don't show the time. You can fix this by editing the screens and copy/paste from screens which have working screens.

2) The snapshot script seems to not grab the full MEDM screen sometimes.

These are not a very big deal, so you can get the microphones working first and we can take care of this afterwards.

  12439   Wed Aug 24 23:47:30 2016 PrafulUpdateElectronicsFinished Prototype Box

Gautam helped me drill holes in a metal box and I set up my circuit inside. Everything seems to be working so far. Tomorrow I'll be suspending the box near the PSL and setting up a data channel. Attached are some pictures of the box- sorry some of the angles turned out weird.

Attachment 1: out1.pdf
out1.pdf
Attachment 2: out2.pdf
out2.pdf
Attachment 3: out3.pdf
out3.pdf
Attachment 4: in1.pdf
in1.pdf
Attachment 5: in2.pdf
in2.pdf
  3116   Thu Jun 24 16:59:24 2010 josephbUpdateVACFinished restoring the access connector and door

[Jenne,  Kiwamu, Steve, Sharmila, Katherine, Joe]

We finished bolting the door on the new ITMX (old ITMY) and putting the access connector section back into place.  We finished with torquing all the bolts to 40 foot-pounds.

  2597   Fri Feb 12 13:56:16 2010 josephbUpdateComputersFinishing touches on IP switch over

The GPIB interfaces have been updated to the new 192.168.113.xxx addresses, with Alberto's help.

Spare ethernet cables have been moved into a cabinet halfway down the x-arm.

The illuminators have a white V error on the alarm handler, but I'm not sure why.  I can turn them on and off using the video screen controls (except for the x arm, which has no computer control, just walk out and turn it on).

There's a laptop or two I haven't tracked down yet, but that should be it on IPs. 

At some point, a find and replace on 131.215.xxx.yyy addresses to 192.168.xxx.yyy should be done on the wiki.  I also need to generate an up to date ethernet IP spreadsheet and post it to the wiki.

 

  4056   Wed Dec 15 12:46:18 2010 KojiSummaryIOOFinishing up the vac work

What else?

v: Edit on Dec 15 10PM
v: Edit on Dec 16 10PM

JD:  We should check OSEMs for all optics *after* table leveling.  Some of them (esp. BS and ITMX) are currently close to their limits right now.

KA: Check green alignment.

Take photos of the tables.

Fix the leveling weights



Location    Optics            Action
--------------------------------------------------------------
@ITMX -     v POX             alignment
            v POP1/POP2       alignment
            v Table Leveling

@ITMY -     POY               mirror replacement (45deg->0deg) / alignment
            v SR2-TT          alignment
            v SRM Tower       alignment / EQ-stop release
            v SRM             alignment
            v SRM OSEM
            vvSRM OPLEV (X2)  install (VIS)/ alignment
            v ITMY OPLEV (X2)   install (VIS)/ alignment
            v OM1/OM2         install (DLC 45deg)/ alignment       
            v Table Leveling

@BSC -      v OM3             install (DLC 45deg/ alignment)
            v OM4(PZT)          neutralize, adjustment
            IPPOS steering    alignment
            v BS OPLEV        alignment
           
v PRM OPLEV(x2)     alignment
            Beam dumps
            Table Leveling

@IMC -      v REFL              mirror replacement (45deg->0deg)

@ETMX -     Al foil removal
            Table Leveling

@ETMY -     ETMY damping
            OSEM
            OPLEV
            Al foil removal
            Table Leveling

@OMC -      v OM5(PZT)        neutralize, adjustment

@ITM/ETM -  Mirror Wiping

  15490   Thu Jul 16 14:41:22 2020 gautamUpdateGeneralFire extinguisher inspection

The (masked) tech accessed all areas in the lab (office area, control room, VEA) between ~230pm-3pm. The laser safety goggles he used have been kept aside for appropriate sanitaiton.

  15530   Mon Aug 17 21:24:43 2020 gautamUpdateGeneralFire extinguisher inspection

A technician came to the lab today at ~4pm. He entered the VEA (with booties and googles), and also the clean and bake lab. The whole procedure lasted ~10 minutes. I did not follow him around, but was available in the control room throughout the process. I think the whole episode went without incident.

BTW, this guy didn't ring the doorbell, I just happened to be here when he came by. I don't know if this is usual practise - are we happy with the technicians entering the VEA and/or clean and bake labs without supervision? AFAIK, this wasn't scheduled.

  1940   Tue Aug 25 02:37:53 2009 ranaConfigurationComputer Scripts / ProgramsFirefox 3.5 installed for 64 bit linux in apps/
Attachment 1: DSC_0620.JPG
DSC_0620.JPG
  5863   Thu Nov 10 16:26:46 2011 MirkoUpdateComputersFirefox kills elog

Had to restart the elog many times. For some reason firefox 8 on Win 7 kills the elog pretty consistently when trying to make a new entry. IE9 works fine ....

  14010   Sat Jun 23 13:08:41 2018 JonUpdateAUXFirst Coherent AUX Scan of PRC Using AM Sidebands

[Jon, Keerthana, Sandrine]

Thu.-Fri. we continued with PRC scans using the AUX laser, but now the "scanned" parameter is the frequency of AM sidebands, rather than the frequency of the AUX carrier itself. The switch to AM (or PM) allows us to coherently measure the cavity transfer as a function of modulation frequency.

In order to make a sentinel measurement, I installed a broadband PDA255 at an unused pickoff behind the first AUX steering mirror on the AS table. The sentinel PD measures the AM actually imprinted on the light going into the IFO, making our measurement independent of the AOM response. This technique removes not only the (non-flat) AOM transfer function, but also any non-linearities from, e.g., overdriving the AOM. The below photo shows the new PD (center) on the AS table.

With the sentinel PD installed, we proceeded as follows.

  • Locked IFO in PRMI on carrier.
  • Locked AUX PLL to PSL.
  • Tuned the frequency of the AUX laser (via the RF offset) to bring the carrier onto resonance with the PRC.
  • Swept the AOM modulation frequency 0-60 MHz while measuring the AUX reflection and injection signals.

The below photo shows the measured transfer function [AUX Reflection / AUX Injection]. The measurement coherence is high to ~55 MHz (the AOM bandwidth is 60 MHz). We clearly resolve two FSRs, visible as Lorentzian dips at which more AUX power couples into the cavity. The SURFs have these data and will be separately posting figures for the measurements.

With the basic system working, we attempted to produce HOMs, first by partially occluding the injected AUX beam with a razor blade, then by placing a thin two-prong fork in the beam path. We also experimented with using a razor blade on the output to partially occlude the reflection beam just before the sensor. We were able to observe an apparent secondary dip indicative of an HOM a few times, as shown below, but could not repeat this deterministically. Besides not having fine control over the occlusion of the beams, there is also large few-Hz angular noise shaking the AS beam position. I suspect from moment to moment the HOM content is varying considerably due to the movement of the AS beam relative to the occluding object. I'm now thinking about more systematic ways to approach this.

 

  14011   Sat Jun 23 20:54:35 2018 KojiUpdateAUXFirst Coherent AUX Scan of PRC Using AM Sidebands

How much was the osc freq of the marconi? And then how much was the resulting freq offset between PSL and AUX?

Are we supposed to see two dips with the separation of an FSR? Or four dips (you have two sidebands)?

And the distance between the dips (28MHz-ish?) seems too large to be the FSR (22MHz-ish).
cf https://wiki-40m.ligo.caltech.edu/IFO_Modeling/RC_lengths

  14017   Tue Jun 26 10:06:39 2018 keerthanaUpdateAUXFirst Coherent AUX Scan of PRC Using AM Sidebands

(Jon, Keerthana, Sandrine)

I am attaching the plots of the Reflected and transmitted AUX beam. In the transmission graph, we are getting peak corresponding to the resonance frequencies, as at that frequency maximum power goes to the cavity. But in the Reflection graph, we are obtaining dips corresponding to the resonance frequency because maximum power goes to the cavity and the reflected beam intensity becomes very less at those points.

 

Attachment 1: TRANS.pdf
TRANS.pdf
Attachment 2: REFL.pdf
REFL.pdf
  7860   Wed Dec 19 21:35:33 2012 KojiSummaryGeneralFirst Contact Training with Margot

[Koji, Steve]

First Contact Training with Margot

  12223   Tue Jun 28 20:43:23 2016 KojiSummaryCOCFirst Contact cleaning practice

Made a dry run of the in-situ cleaning for a 3inch optic.

Attachment 1: The Al dummy mass is clamped in the suspension cage.
Attachment 2: The front surface was painted. The nominal brush with the FC bottle was used.
Attachment 3: Zoom in of the front surface.
Attachment 4: The back surface was painted.
Attachment 5: The back surface was peeled.
Attachment 6: The front surface was peeled too.
Attachment 7: The peeled layers.

Findings:

1. To paint a thick layer (particlarly on the rim) is the key to peel it nicely.

2. It was helpful for easier peeling to have mutiple peek tabs. Two tabs were sufficient for ~1" circle.

3. The nominal brush with the bottle was OK although one has to apply the liquid many times to cover such a large area. A larger brush may cause dripping.

4. The nominal brush was sufficiently long once the OSEMs are removed. In any case it is better to remove the OSEMs.

Attachment 1: IMG_20160628_170335196.jpg
IMG_20160628_170335196.jpg
Attachment 2: IMG_20160628_171547769.jpg
IMG_20160628_171547769.jpg
Attachment 3: IMG_20160628_171607802.jpg
IMG_20160628_171607802.jpg
Attachment 4: IMG_20160628_172328190.jpg
IMG_20160628_172328190.jpg
Attachment 5: IMG_20160628_174541960.jpg
IMG_20160628_174541960.jpg
Attachment 6: IMG_20160628_174556004.jpg
IMG_20160628_174556004.jpg
Attachment 7: IMG_20160628_174617198.jpg
IMG_20160628_174617198.jpg
  7388   Fri Sep 14 16:39:14 2012 ericq, jenneUpdateGeneralFirst In Vac Picture

After much fussing, we got a picture of MMT1 with the beam.

Using the iris doesn't seem feasible. Since it has to be significantly separated from the optic, it is hard to judge whether it is centered, especially in yaw.

It took ~30 min to get this picture. Comments on whether this kind of picture is good enough are welcomed, since there are many more to be taken.

Attachment 1: mmt1.jpg
mmt1.jpg
  7389   Fri Sep 14 18:15:43 2012 ranaUpdateASCFirst In Vac Picture

 

 Looks good. Any way that you can tell in an unambiguous way, where the beam is, is very good. Ideally we want to have1-2 mm accuracy.

  2964   Fri May 21 00:51:06 2010 JenneUpdateIOOFirst MC mode measuring (hopefully) done

[Jenne, Kiwamu, Steve]

Round 1 of measuring the MC mode is pretty much done.  Yay.

Earlier today, Steve and I launched the MC beam off the flat mirror just after the Faraday, and sent it down toward ETMY(new convention). We ended up not being able to see it all the way at the ETM because we were hitting the beam tube, but at the ITM chamber we could see that the beam looked nice and circle-y, so wasn't being clipped in the Faraday or anywhere else.  To do this we removed 2 1inch oplev optics.  One was removed from the BS table, and wrapped in foil and put in a plastic box.  The other was just layed on its' side on the BS table. 

I then took the beam out of the BS chamber, in order to begin measuring the mode.  I left the flat fixed mirror in the place of what will be PZT SM1, and instead used the PZT mirror to turn the beam and get it out the BS chamber door.  (Thoughts of getting the beam to the BS oplev table were abandoned since this was way easier, since Kiwamu and Steve had made the nifty table leg things.)  Kiwamu and I borrowed an 2inch 45P Y1 optic from the collection on Koji's desk (since we have ZERO 2inch optics on the random-optic-shelf....no good), to shoot the beam down the hallway of the Yarm (new convention).  We used the beam scan on a rolling cart to measure the beam at various distances.  I made some sweet impromptu plum bobs to help make our distance measurements a bit more accurate.

We stopped at ~25 feet from the BS chamber, since the spot was getting too big for the beam scanner.  If it turns out that I can't get a good fit with the points I have, I'll keep everything in-chamber the same, and do the farther distances using the good ol' razor blade technique.

I have measurements for the distances between the beam scan head and the opening of the BS chamber.  Tomorrow, or very soon after, I need to measure the distances in-chamber between the MC and the BS chamber opening.  Plots etc will come after I have those distances.

Next on the to-do list:

1.  Measure distances in-chamber for first mode scan.

2.  Plot spot size vs. distance, see if we need more points. Take more points if needed.

3.  Put in MMT1, repeat measurements.

4.  Put in MMT2, rinse and repeat.

5.  Move the PZT mirror to its new place as SM1, and figure out how to connect it.  Right now the little wires are hooked up on the BS table, but we're going to need to make / find a connector to the outside world from the IOO table. This is potentially a pretty big pain, if we don't by happenstance have open connectors on the IOO table.

  5727   Fri Oct 21 18:20:54 2011 MirkoUpdateCDSFirst OAF version running

[Jenne, Jamie, Mirko]

We got the first version of the oaf code based on Matt"s code running!! :-)
Produces already data for e.g. MICH DOF. But don"t trust that. It's only 10 taps long and delay is not adjusted.

  16007   Thu Apr 8 17:04:43 2021 Anchal, PacoUpdateSUSFirst Successful Coil Balancing

Today, we finally crossed the last hurdle and got a successful converging coil balancing run. laugh


What was the issue with POS?

  • Position of the MC2 mirror is being sensed using C1:IOO-MC_F_DQ channel which is proportional to the resonant frequency of the locked IMC.
  • However, this sensor is always 180 degrees out of phase of our actuator, the coils.
  • When the coils push the mirror forward, the length of the cavity actually decreases.
  • We added an extra option of providing a sign to the sensors such that -1 will be multiplied to sensed values for sensors which measure in opposite direction to the actuation.
  • This is important, because the feedback is applied to the coil output matrix assuming a particular direction of acctuation.
  • When we gave negative sign for the position sensor, it all started making sense and the algorithm started converging.

First run parameters:

  • We used binwidth of 0.25 Hz and duration of excitation as 41s. This would give welch and csd averaging of 19. We used median averaging to ignore outliers.
  • This iteration was run after PIT and YAW were separetly uncoupled before. We'll post a clean start to end run results in near future.
  • The iteration works in following manner:
    • Define a constant coil matrix C = [[1, 1, 1], [1, 1, -1], [1, -1, 1], [1, -1, -1]] which is ideal coil output matrix.
    • In each iteration, the output matrix Ok is defined as (note @ is the matmul operator):
      Ok = C @ Ak
      where Ak is a 3x3 matrix. A-1 is identity matrix.
    • At the end of each iteration, a sensing matrix is calculated in dimensions sensedDOF x excitedDOF, Sk
    • For next iteration, Ak+1 is calcualted by:
      Ak+1 = Ak - b * (Sk - I)
      where I is the identity matrix.
    • At convergence, the sensing matrix would become same as identity and matrix A will stop updating.
  • For this run, we kept the parameter b to be 0.05. This is similar to the KP parameter in PID loops. It should be between 0 and 1.
  • Since b value was small enough to allow for convergence from the inital point, but later it slowed down the process a lot.
  • Ideally, we should figure out a way to increase this paramter when the coil has been balanced somewhat, to increase the speed of the algorithm.
  • Secondly, we have a code which excites all DOFs at different frequencies directly using excitation channels in coil output matrix using awg.py. But for some reason, the excitation channel for 4th row in the output matrix column only connects intermittantly. Because of this, we can't use this method reliably yet. We can investigate more into it if suggested.

Balancing characteristics:

  • Attachment 1 shows how the distance of sensing matrix falls as iterations increase. We only ran for 50 iterations.
  • Attachment 2 shows how different off-diagonal terms of sensing matrix decreased.
    • Note that POS -> PIT, POS -> YAW and PIT-YAW have settled down to the noise floor.
    • The noise floor can be improved by increasing the excitation amplitude and/or increasing the duration of measurement.
  • Attachment 3 shows the evolution of sensing matrix as iterations move.

Final balanced output matrix:

Final balanced output coil matrix for MC2
POS PIT YAW COILS
1.02956 1.13053 1.19116 UL
1.01210 1.09188 -0.74832 UR
0.98737 -0.85502 0.70485 LR
0.96991 -0.89366 -1.23463 LR
Final Sensing Matrix
  Exc POS Exc PIT Exc YAW
Sens POS 1 -2.96e-2 8.00e-3
Sens PIT 8.58e-4 1 -4.84e-3
Sens YAW 5.97e-4 -1.15e-3 1

Code features and next:

  • Majority of the code is in two files: scripts/SUS/OutMatCalc/MC2crossCoupleTest.py and scripts/SUS/OutMatCalc/crossCoupleTest.py .
  • The code runs from start to end without human involevement and restores the state of channels in any case (error, kyboard interrupt, end of code) using finally statement.
  • Currently, each excitation is done one at a time through LockIn1. As mentioned above, this can be sped up 3 times if we get the awg.py to work reliably.
  • The complete code is in python3 and currently is run through native python3 on allegra (a new debian10 workstation with latest cds-workstation installed).
  • The code can be easily generalized for balancing any optic. Please let us know if we should work on making the generalized optic.
  • We're also working on thinking about increasing b as iterations move forward and the error signal becomes smaller.
  • We can also include the uncertainty in the Sensing matrix measurement to provide a weighted feedback. That way, we can probably increase b more.
Attachment 1: SDistanceFromIdentity.pdf
SDistanceFromIdentity.pdf
Attachment 2: SmatIterations.pdf
SmatIterations.pdf
Attachment 3: SmatrixPlots.pdf
SmatrixPlots.pdf SmatrixPlots.pdf SmatrixPlots.pdf SmatrixPlots.pdf SmatrixPlots.pdf SmatrixPlots.pdf SmatrixPlots.pdf SmatrixPlots.pdf
  565   Wed Jun 25 11:36:14 2008 Max JonesUpdateComputer Scripts / ProgramsFirst Week Update
For the first week I have been modifying the noise budget script in caltech/NB to run with 40 m parameters and data. As per Rana's instructions, I have tried to run the script with only seismic and Darm sources. This involves identifying and changing channel names and altering paramter files (such as NB/ReferenceData/C1IFOparams.m). To supply the parameter files, I have copied the H1 files with (as yet only) slight modification. The channel name changes have been made to mirror the sites for the most part. Two figures are attached which show the current noise budget. The Day plot was taken 6/23/08 at ~10:30 am. The night plot was taken 6/22/08 at ~11:00 pm . Note that the SRD curve is for the sites and not for the 40 m (I hope to change that soon). Also in one of the plots the DARM noise signal is visible. Obviously this needs work. A list of current concerns is

1) I am using a seismic transfer function made by previous SURF student Ryan Kinney to operate with channels of the form C1:PEM_ACC-ETMY_Y (should I be using C1:DMF-IX_ACCY?) and the channels I am currently using are the acceleraometers for the mode cleaner with names of the form C1:PEM_ACC-MC1_X. Rana said that he thinks these may be the same but I need to be sure.

2) We don't have a DARM_CTRL channel but the code requires it, currently I am using DARM_ERR as a substitute which is probably partly responsible for the obvious error in DARM noise.

Any suggestions are appreciated. Max.
Attachment 1: C1_NoiseBudgetPlot_Day.eps
C1_NoiseBudgetPlot_Day.eps
Attachment 2: C1_NoiseBudgetPlot_Night.eps
C1_NoiseBudgetPlot_Night.eps
  6879   Wed Jun 27 11:33:28 2012 MashaUpdateGeneralFirst Week Update

This week I wrote Matlab code, most of which can be found in /users/masha

First, I wrote a simulation seismicFilter.m which filters noisy seismic noise with a desired signal of non-seismic noise. The signals are purely simulated, so I played around with zero-pole-gain generation of transfer functions to obtain them. The function takes the number of taps, the filter type (Wiener or adaptive nlms) as well as an iteration step size and number of iterations, and generates PSD plots of the witness signal, the desired signal, the estimated (filtered) signal, and the error. I'm not sure that I am properly implementing the Wiener part of the code, and I assume the line "[W, R, P] = miso_firlev(TAPS, noisySeismicSignal1, seismicSignal2); " generates W, a filter with TAPS number of weights, but  then "[y, error] = filter(W, 1, noisySeismicSignal1);" generates an error signal of size TAPS rather than N, the size of the original signal. Perhaps I should calculate error using e(t) = d(t + a) - w(t)*x(t), where "a" is the delay.

I have various screenshots in my directory of what seismicFilter.m generates, and I will take a larger screenshot, as well as generate a learning curve (for error vs. number of taps) when I can use Sasha's computer for a bit, since it both has more computing power and a larger screen.

The funciton filterConvergence.m, meanwhile, is similar, except it takes two file names as real data, and uses realDataFilter.m to run the filtering. Currently, I am working with data from C1:IOO-MC_F_DQ-Online  and C1:PEM-SEIS_GUR1_X_IN1_DQ-Online, and I will include screenshots of these once I get on Sasha's computer.

In order to generate the data, meanwhile, I had to modify the python script, and thus wrote mashaImportingData.py for myself. Likewise, plotSignalFromFile.m visualizes this data, both in the time domain and in the frequency domain.

On the side, I wrote an NLMS filter in adaptiveFilterSimulationNLMS.m, and compared is to Matlab's NLMS filter in NLMStest.m (using generated data) and adaptiveFilterSimulation.m using twn input signals. Right now, it's faster on smaller inputs and smaller tap sizes, but then begins to choke and run slower than the Matlab one when these get to big. In order to improve it, I have to develop a better method of generating the initial weights.

As far as machine learning goes, once I find the number of taps for the convergence of both the Wiener filter and the NLMS filter, I will email Denis for further instructions. At some point, however, I should generate learning sets from the seismometers and the MCL (or the DARM), and thus have to find adequate times at which I can take data (probably not from the DARM, however, because it was rarely on).

Thanks for reading!

  6878   Wed Jun 27 11:27:49 2012 LizUpdate First Week Update!

This week, the other SURF students and I got acquainted with the caltech campus, LIGO 40m lab and the expectations of the SURF program.  We went to a lot of safety meetings and lectures that established a framework for the jobs we will be doing over the course of the summer.  I went on several tours of the 40m interferometer (one each with Jenne, Jamie and Steve) to get an overview of the layout and specifics of the setup.  I read parts of R. Ward and A. Parameswaran's theses and Saulson's book in order to prepare myself and gain a broader understanding of the purpose of LIGO.

I also began working in Python this week, primarily graphing PSDs of data from the C1:SUS-ETMY_SENSOR_LR, C1:SUS-ETMY_SENSOR_LL, C1:SUS-ETMY_SENSOR_UR, and C1:SUS-ETMY_SENSOR_UL channels.  I will eventually be using Python to generate the plots for the summary pages, so this is good practice.  The code that I have been working on can be found in /users/elizabeth.davison/script5.py.  Additionally, I have been going through the G1 summary pages and attempting to understand the plots available on them and the code that is available.

My plans for the upcoming week begin with modifying my code and potentially calibrating the channel data so that it is in units of length instead of counts.  I will also access the code from the G1 pages and go over it in depth, hopefully gaining insight into the structure of the website.

  3099   Tue Jun 22 20:07:08 2010 JenneUpdate40m UpgradingFirst attempt at Tip Tilt hanging

[Jenne, Steve, Nancy, Gopal]

We made an attempt at hanging some of the Tip Tilt eddy current dampers today. 

Photo 1 shows the 2 ECDs suspended.

Procedure:

(1) Loosen the #4-40 screws on the side of the ECDs, so the wire can be threaded through the clamps.

(2) Place the ECDs in the locator jigs (not shown), and the locator jigs in the backplane (removed from main TT structure), all laying flat on the table.

(3) Get a length of Tungsten wire (0.007 inch OD = 180um OD), wipe it with acetone, and cut it into 4 ~8cm long segments (long enough to go from the top of the backplane to the bottom).

(4) Thread a length of wire through the clamps on the ECDs, one length going through both ECDs' clamps.

(5) One person hold the wire taught, and straight, and as horizontal as possible, the other person tightened the clamping screws on the ECDs.

(6) Again holding the wire in place, one person put the clamps onto the backplane (the horizontal 'sticks' with 3 screws in them).

(7) The end. In the future, we'll also clip off extra pieces of wire.

When we held up the backplane to check out our handy work, it was clear that the bottom ECD was a much softer pendulum than the top one, since the top one has the wire held above and below, while the bottom one only has the wire held on the top.  I assume we'll trim the wire so that the upper ECD is only held on the top as well?

Lessons learned:

* This may be a 3 person job, or a 2 people who are good at multitasking job.  The wire needs to be held, the ECDs need to be held in place so they don't move during the screwing/clamping process, and the screws need to be tightened.

* Make sure to actually hold the wire taught. This didn't end up happening successfully for the leftmost wire in the photo, and the wire is a bit loose between the 2 ECDs.  This will need to be redone.

* We aren't sure that we have the correct screws for the clamps holding the wire to the backplane.  We only have 3/16" screws, and we aren't getting very many threads into the aluminum of the backplane.  Rana is ordering some 316 Stainless Steel (low magnetism) 1/4" #4-40 screws.  We're going for Stainless because Brass (the screws in the photo), while they passed their RGA scan, aren't really good for the vacuum.  And titanium is very expensive.  

The 2nd photo is of the magnet sticking out of the optic holder.  The hole that the magnet is sitting in has an aluminum piece ~2/3 of the way through.  A steel disk has been placed on one side, and the magnet on the other.  By doing this, we don't need to do any press-fitting (which was a concern whether or not the magnets could withstand that procedure), and we don't need to do any epoxying.  We'll have to wait until the ECDs are hung, and the optic holder suspended, to see whether or not the magnet is sticking out far enough to get to the ECDs. 

Attachment 1: 2_ECDs_small.jpg
2_ECDs_small.jpg
Attachment 2: MagnetStickingOutFar_small.jpg
MagnetStickingOutFar_small.jpg
  5639   Sun Oct 9 17:13:46 2011 kiwamuUpdateLSCFirst attempt to estimate mode matching efficiency using interferometer

The efficiency of the mode matching (MM) to PRC (Power-Recycling Cavity) has been estimated by using the interferometer.

The estimated MM efficiency is about 74 % when losses in the cavity are assumed to be zero.

 

(Motivation)

 There had been an issue that the recycling gain didn't go to the designed high value of about 42  (#5541).
One of the possibilities is a low efficiency in the MM to PRC (also see #5541).
Although the MM efficiency had been measured using a beam scan ( see a summary on the wiki) a long time ago, it haven't been verified.
Therefore the MM has to be reviewed by using the real interferometer.

(Measurement)

 The concept of this measurement is observe the amount of the reflected light from a power-recycled cavity and estimate the MM efficiency based on the measured reflectivities.
 Since using the real PRC (consisting of BS, ITMs and PRM) could be a too complicated system for this measurement,
simpler cavities, namely Power-Recycled ITMX and ITMY (PRX and PRY), were used to examine the MM efficiency.
 The measurement goes in the following order :
    (1) Measurement of the amount of the single-bounce reflection from PRM with BS and ITMs misaligned.
    (2) Lock PRX or PRY to carrier resonance.
    (3) Alignment of PRX/Y to maximize the intracavity power. This time ASDC was used as a monitor of the intracavity power.
    (4) Measurement of the amount of the reflected light when the cavity is in resonance. The value in REFLDC was averaged in 100 sec.
     => done by tdsavg 100 C1:LSC-REFLDC_OUT
The same measurement was performed for both PRX and PRY.
 
 - locking parameters -
  Sensor = REFL11_I
  Whitening gain = 10 (30 dB)
  PRCL_GAIN = 2
  UGF ~ 200 Hz

(Analysis)

In order to estimate the relation between the MM efficiency vs. the reflected light, two models are considered:
   (1) simple model => no loss and no sidebands
   (2) sideband-included model => no loss but sidebands are taken into the account of the reflection.
 
(1) In the simple model the reflectivity Prefl / Pin is expressed by
         [Reflectivity]  = Prefl / Pin = Z * Rcav +  (1- Z) * Rprm
 
where Z is MM efficiency and Rprm is the reflectivity of PRM
and Rcav is the reflectivity of PRX/Y when it's resonance and it is defined by
         Rcav = | rprm - ritm t2BS|2 / |1 -rprm ritm t2BS |2
 
Tprm = 5.75% and Titm = 1.4 % are assumed in all the calculations.
In the first equation the first term represents the mode matched light and hence it couples with the cavity through Rcav.
The second term is the non-mode-matched light and because they are not interacting with the cavity they will be simply reflected by PRM through Rprm.
 
(2) In reality two phase-modulated light (11 MHz and 55 MHz) will behave differently from the carrier.
For example when the carrier is in resonance the sidebands will be anti-resonance against the cavity.
So that the amount of REFLDC will be slightly bigger because of the reflection of the sidebands.
 
      Prefl = Z * Rcav * Pc + Z * Ranti * Ps +  (1- Z) * Rprm * (Pc + Ps)
 
where Pc and Ps are the power in the carrier light and the sidebands respectively.
And Ranti is the reflectivity of the anti-resonance PRX/Y, which can be obtained by replacing the minus sign by the plus sign in the equation of Rcav shown above.
It is assumed that the sum of the carrier power and sidebands power is the incident power Pin = Pc + Ps.
The power in the carrier and the sidebands were estimated based on the OSA measurement (#5519), so that
          Pc / Pin = |J0(0.14)|2 * |J0(0.17)|2 =  0.976
          Ps / Pin = 2 * |J1(0.14)|2 + 2 * |J1(0.17)|2 =  0.024
 

(Results)

Here are the measured values in REFLDC

 -- Measurement 1 : PRX
    Single bounce from PRM = 4802.27 counts
       ==> the incident power = 5095.25 counts
    Reflected light from PRX = 4433.88 counts
      ==> Reflectivity = 0.8702
 
-- Measurement 2 : PRY
    Single bounce from PRM = 4833.05 counts
       ==> the incident power = 5127.05 counts
    Reflected light from PRX = 4444.48 counts
      ==> Reflectivity = 0.86672
 
On average the reflectivity of power-recycled ITM cavity was 0.868 with a standard deviation of  0.001744.
Actually the standard deviation estimated here is not fair because the measurement was done by only twice,
but my intention was that I wanted to see how the error can affect the estimation of the MM efficiency.
Here is a plot comparing the model curves and the measured values with 5 sigma error box (5 times of measured standard deviation).

mm_reflection.png

It is shown that the mode matching efficiency is 73.7 % when the sideband-included model is considered.
With the 5 sigma deviation it can go from 65% to 82% but it is still low and lower than the beam scan measurement ( see a summary on the wiki). 

Anyways the estimated MM efficiency with the sidebands effect included and without loss effect is

        MM efficiency = 73.7 +/- 1.7 % (1 sigma error)  or +/- 8.7 % (5 sigma error)

  5640   Mon Oct 10 00:01:26 2011 KojiUpdateLSCFirst attempt to estimate mode matching efficiency using interferometer

"^2"s are missing in the second equation, but the calculation results seem correct.

PRX and PRY have different mode matching because of the Michelson asymmetry.
Are individually estimated mode matching indicates any sign of reasonable mode mismatch?
(The difference can be very small because the asymmetry is not so big.)

  13250   Thu Aug 24 18:02:16 2017 GabrieleSummaryLSCFirst cavity length reconstruction with a neural network

1) Introduction

In brief, I trained a deep neural network (DNN) to recosntuct the cavity length, using as input only the transmitted power and the reflection PDH signals. The training was performed with simulated data, computed along 0.25s long trajectories sampled at 8kHz, with random ending point in the [-lambda/4, lambda/4] unique region and with random velocity.

The goal of thsi work is to validate the whole approach of length reconstruction witn DNN in the Fabry-Perot case, by comparing the DNN reconstruction with the ALS caivity lenght measurement. The final target is to deploy a system to lock PRMI and DRMI. Actually, the Fabry-Perot cavity problem is harder for a DNN: the cavity linewidth is quite narrow, forcing me to use very high sampling frequency (8kHz) to be able to capture a few samples at each resonance crossing. I'm using a recurrent neural network (RNN), in the input layers of the DNN, and this is traine using truncated backpropagation in time (TBPT): during training each layer of RNN is unrolled into as many copies as there are input time samples (8192 * 0.25 = 2048). So in practice I'm training a DNN with >2000 layers! The limit here is computational, mostly the GPU memory. That's why I'm not able to use longer data stretches.

But in brief, the DNN reconstruction is performing well for the first attempt.

2) Training simulation

In the results shown below, I'm using a pre-trained network with parameters that do not match very well the actual data, in particular for the distribution of mirror velocity and the sensing noises. I'm working on improving the training.

I used the following parameters for the Fabry-Perot cavity:

The uncertaint is assumed to be the 90% confidence level of a gaussian distribution. The DNN is trained on 100000 examples, each one a 0.25/8kHz long trajectory with random velocity between 0.1 and 5 um/s, and ending point distributed as follow: 33% uniform on the [-lambda/4, lambda/4] region, plus 33% gaussian distribution peaked at the center with 5 nm width. In addition there are 33% more static examples, distributed near the center. 

For each point along the trajectory, the signals TRA, POX11_I and POX11_Q are computed and used as input to the DNN.

3) Experimental data

Gautam collected about 10 minutes of data with the free swinging cavity, with ALS locked on the arm. Some more data were collected with the cavity driven, to increase the motion. I used the driven dataset in the analysis below.

3.1) ALS calibration

The ALS signal is calibrated in green Hz. After converting it to meters, I checked the calibration by measuring the distance between carrier peaks. It turned out that the ALS signal is undercalibrated by about 26%. After correcting for this, I found that there is a small non-linearity in the ALS response over multiple FSR. So I binned the ALS signal over the entire range and averaged the TRA power in each bin, to get the transmission signals as a function of ALS (in nm) below:

I used a peak detection algorithm to extract the carrier and 11 MHz sideband peaks, and compared them with the nominal positions. The difference between the expected and measured peak positions as a function of the ALS signal is shown below, with a quadratic fit that I used to improve the ALS calibration

The result is

z_initial = 1e9 * L*lamba/c *1.26. * ALS
z_corrected = 2.1e-06 z^2  -1.9e-02 z  -6.91e+02

The ALS calibrated z error from the peak position is of the order of 3 nm (one sigma)

3.2) Mirror velocity

Using the calibrated ALS signal, I computed the cavity length velocity. The histogram below shows that this is well described by a gaussian with width of about 3 um/s. In my DNN training I used a different velocity distribution, but this shouldn't have a big impact. I'm retraining with a different distirbution.

4) DNN results

The plot below shows a stretch of time domain DNN reconstruction, compared with the ALS calibrated signal. The DNN output is limited in the [-lambda/4, lambda/4] region, so the ALS signal is also wrapped in the same region. In general the DNN reconstruction follows reasonably well the real motion, mostly failing when the velocity is small and the cavity is simultanously out of resonance. This is a limitation that i see also in simulation, and it is due to the short training time of 0.25s.

I did not hand-pick a good period, this is representative of the average performance. To get a better understanding of the performance, here's a histogram of the error for 100 seconds of data:

The central peak was fitted with a gaussian, just to give a rough idea of its width, although the tails are much wider. A more interesting plot is the hisrogram below of the reconstructed position as a function of the ALS position, Ideally one would expect a perfect diagonal. The result isn't too far from the expectation:

The largest off diagonal peak is at (-27, 125) and marked with the red cross. Its origin is more clear in the plot below, which shows the mean, RMS and maximum error as a function of the cavity length. The second peak corresponds to where the 55 MHz sideband resonate. In my training model, there were no 55 MHz sidebands nor higher order modes. 

5) Conclusions and next steps

The DNN reconstruction performance is already quite good, considering that the DNN couldn't be trained optimally because of computation power limitations. This is a validation of the whole idea of training the DNN offline on a simulation and then deploy the system online.

I'm working to improve the results by

  • training on a more realistic distribution of velocity
  • adding the 55 MHz sidebands
  • maybe adding HOMs
  • tune the DNN architecture

However I won't spend too much time on this, since I think the idea has been already validated.

 

  13251   Thu Aug 24 18:51:57 2017 KojiSummaryLSCFirst cavity length reconstruction with a neural network

Phenomenal!

  13256   Sat Aug 26 09:56:34 2017 GabrieleSummaryLSCFirst cavity length reconstruction with a neural network

Update

I included the 55 MHz sideband and higher order modes in my training examples. To keep things simple, I just assumed there are higher order modes up to n+m=4 in the input beam. The power in each HOM is randomly chosen from a random gaussian distribution with width determined from experimental cavity scans. I used a value of 0.913+-0.01 rad for the Gouy phase (again estimated from cavity scans, but in reasonable agreement with the nominal radius of curvature of ETMX)

Results are improved. The plot belows show the performance of the neural network on 100s of experimental data

For reference, the plots below show the performance of the same network on simulated data (that includes sensing noise but no higher order modes)

  13258   Mon Aug 28 08:47:32 2017 JamieSummaryLSCFirst cavity length reconstruction with a neural network
Quote:

Phenomenal!

truly.

  12065   Wed Apr 6 17:52:21 2016 gautamUpdateendtable upgradeFirst contact cleaning commenced

I've begun cleaning the optics that will eventually go back onto the newly installed X-endtable. We decided that First Contact was the way to go (as opposed to methanol drag wiping). Koji demonstrated the application of the (red) First Contact solution onto a 2" mirror - I then proceeded to work on the rest of the optics. We are broadly following the procedure in E1000079 - first one coat of First Contact solution is applied, then a small piece of PEEK is embedded by applying a second layer of solution over it (this will enable us to pull off the First Contact once we are ready - the plan is to do this after roughly placing the optic on the table. As of now, I've finished coating most of the optics that are part of the IR Transmon path - I will continue later in the evening.

The new endtable is almost ready for re-population. Steve just needs to shim the enclosure which will be done tomorrow morning. The game-plan as discussed at the meeting today is to first try and set up the IR Transmon path. This will allow us to verify that the endtable height is such that we can maintain a beam height of 4" everywhere on the table (I suspect we may have to compromise at some poing and do some fine adjustment of 1/4 to 1/2" somewhere though). It will also allow me to define the cavity axis relative to the table, which will be useful to place the green steering optics eventually. Doing this will be challenging though as right now, I can't see any of the arm flashes on the endtable using an IR card. Ideally, we want to somehow lock the X arm and then do the checks mentioned at the endtable, before beginning to put the endtable back together. 

  14613   Thu May 16 13:07:14 2019 gautamUpdateSUSFirst contact residue removal

I  used a pair of tweezers to remove the stray fiber of first contact. As Koji predicted, this was rather dry and so it didn't have the usual elasticity, so while I was able to pull most of it off, there is a small spot remaining on the HR surface of the ETM. We will remove this with a fresh application of a small patch of FC.

I the meantime, I'm curious if this has actually fixed the suspension woes, so yet another round of freeswinging data collection is ongoing. From the first 5 mins, looks positive, I see 4 peaks around 1Hz cool!

The following optics were kicked:
ETMY
Thu May 16 13:06:39 PDT 2019
1242072418

Update 730pm: There are now four well-defined peaks around 1 Hz. Together with the Bounce and Roll modes, that makes six. The peak at 0.92 Hz, which I believe corresponds to the Yaw eigenmode, is significantly lower than the other three. I want to get some info about the input matrix but there was some NDS dropout and large segments of data aren't available using the python nds fetch method, so I am trying again, kicked ETMY at 1828 PDT. It may be that we could benefit from some adjustment of the OSEM positions, the coupling of bounce mode to LL is high. Also the SIDE/POS resonances aren't obviously deconvolved. The stray first contact has to be removed too. But overall I think it was a successful removal, and the suspension characteristics are more in line with what is "expected". 

Attachment 1: etmy_sensors.pdf
etmy_sensors.pdf
Attachment 2: etmy_BRmode.pdf
etmy_BRmode.pdf
  43   Thu Nov 1 01:28:04 2007 waldmanOtherOMCFirst digital lock of OMC
[Pinkesh, Sam]

We locked a fiber based NPRO to the suspended OMC tonight using the TPT digital control system. To control the laser frequency, we took the PZT AI output and ran it on a BNC cable down the hallway to the Thorlabs HV box. The Thorlabs is a singled ended unit so we connected the AI positive terminal only and grounded the BNC to the AI shield. We could get a -6 to 1.5 V throw in this method which fed into the 10 k resisotr + 9 V battery at the input of the HV box. The HV out ran to the NPRO PZT fast input.

We derived our error signal from a PDA255 in reflection with a 29.5 MHz PDH lock. The signal feeds into one of the unused Tip/Tilt AA channels and is passed to the PZT LSC drive through the TPT_PDH1 filter bank. In the PZT_LSC filter we put a single pole at 1 Hz which, together with the phase we mentioned the other night (180 degrees at 3 kHz) should allow a 1 kHz-ish loop. In practice, as shown below, we got a 650 Hz UGF with 45 degrees of phase margin and about 6 dB of gain margin.

The Lower figure shows the error point spectrum with 3 settings. REF0 in blue shows lots of gain peaking at 1.5 kHz-ish, just where its expected - the gain was -40. The REF1 has gain of -20 and shows no gain peaking. The current trace in red shows some gain peaking cuz the alignment is better but it also has included a 1^2:20^2 boost which totally crushes the low frequency noise. We should do a better loop sweep after getting the alignment right so we can see how much boost it will really take.

Just for fun, we are leaving it locked overnight and recording the PZT_LSC data for posterity.
Attachment 1: 071101_PZT_firstLoopSweep.pdf
071101_PZT_firstLoopSweep.pdf
Attachment 2: 071101_PZT_firstLoopSweep.gif
071101_PZT_firstLoopSweep.gif
Attachment 3: 071101_OMC_FirstLock_spectra.pdf
071101_OMC_FirstLock_spectra.pdf
Attachment 4: 071101_OMC_FirstLock_spectra.gif
071101_OMC_FirstLock_spectra.gif
  13932   Fri Jun 8 01:08:22 2018 johannesUpdatePSLFirst light of AUX at YEND

Among the things that we hadn't taken care of yesterday before beginning to look for transmission signals were the polarization of the AUX beam on the AS table and optimizing the PLL feedback. The AUX beam is s-polarized on the PSL table (choice due to availablility of mirrors), and I added a half waveplate in front of the fiber to match it's axes. I placed another half-waveplate at the fiber output and send the reflection port of a PBS cube onto a PDA1CS photodetector. By alternatingly turning the waveplates I minimized the reflected light, giving strongly p-polarized light on the AS table for best results when interfering with the IFO beam. I wiggled the fiber and found no strong dependency of the output polarization on fiber bending. Attachment 2 shows the current layout.

The beat signal between AUX and PSL table is at -20dBm, and I adjusted the PLL gain and PI-corner to get reliable locking behavior. I think it's a good idea to keep the AUX beam on the AS table blocked while it's not in use, and only unblock it when it is phaselocked to avoid a rogue beam with no fixed phase relation to the PSL in the IFO.I blocked the beam after completing this work today.

I used the signal chain that Keerthana, Koji, and I set up yesterday to look for mode flashed of the AUX light in the YARM using the RF beat with the PSL carrier in transmission. To align the AUX beam to the arm the following steps were performed:

  1. Using a spectrum analyzer to look at the RF power at the target frequency between frequency-shifted AUX beam and PSL carrier on AS110, align the beam using the mirror pair closest to the fiber coupler for maximum signal.
  2. Initiate a sweep of the PLL LO frequency sourced by the Marconi using GPIB scripts over about 1 FSR. A strong peak was visible at ~31.76 MHz offset frequency
  3. Tune and hold LO frequency (in this case at 48.2526 MHz) such that AUX beam resonates in the arm. Optimize alignment by maximizing RF signal on PD in transmission.

This was followed by a sweep over two full FSRs. Attachment #1 shows the trace recorded by the AG4395 using the max data hold setting during the sweep. Essentially the beat between AUX and PSL carrier traced out the arm's transmission curve. At minimum transmission there was still a ~82dB beat on the transmission PD visible.

The YEND QPD is currently blocked and sees no light.

Attachment 1: AG4395A_07-06-2018_205019.pdf
AG4395A_07-06-2018_205019.pdf
Attachment 2: PSL_AUX_SETUP.pdf
PSL_AUX_SETUP.pdf
Attachment 3: AS_AUX_SETUP.pdf
AS_AUX_SETUP.pdf
  15534   Thu Aug 20 00:21:51 2020 gautamUpdateElectronicsFirst look at HV coil driver

Summary:

A single channel of this board was stuffed (and other channels partially populated). The basic tests passed, and nothing exploded! Even though this is a laughably simple circuit, it's nice that it works.

HV power supplies:

A pair of unused KEPCO BHK300-130 switching power supplies that I found in the lab were used for this test. I pulled the programmable cards out at the rear, and shorted the positive output of one unit to the negative of the other (with both shorted to the supply grounds as well), thereby creating a bipolar supply from these unipolar models. For the purposes of this test, I set the voltage and current limits to 100V DC, 10mA respectively. I didn't ramp up the supply voltage to the rated 300 V maximum. The setup is shown in Attachment #1.

Tests:

  1. With the input to the channel shorted to ground, I confirmed with a DMM that the output was (nearly) zero (there was an offset of ~40mV but I think this is okay).
  2. Used the calibrated voltage source, and applied +/- 3 V in steps of ~0.5 V, while monitoring the output with a DMM. Confirmed the output swing of ~ +/-90 V, which is what is expected, since the design voltage gain of this circuit is 31.
  3. Drove a 0.1 Hz, 500mVpp sine wave at the input while monitoring the output and the Vmon testpoints, see Attachment #2. Note the phasing between input and output, and also the fact that the gain is slightly lower than the expected gain of 31, because there are three poles at ~0.7 Hz, which already start showing some influence on the transfer function at 0.1 Hz.
  4. Noise measurement 
    • The whole point of this circuit is to realize sub 1pA/rtHz current noise to the coil, when it is connected.
    • For this test, no load was connected (i.e. voltage noise was measured at the output of the 25 kohm resistor), and the input was shorted to ground so that the DC value of the output was close to 0 (the idea was to not overload the SR560/SR785 with high voltage).
    • An SR560 preamp with gain x50 (DC coupled) was used to preamplify the signal. This was the maximum gain that could be used with the unit DC coupled, due to the small DC offset. I opted to keep the DC coupling to get a look at the low frequency noise as well, but in hindsight, maybe I should have used AC coupling as we only care about the current noise at ~100 Hz.
    • See Attachment #3 for results. The measurement is close to the model above ~100 Hz

Need to think more about how to better characterize this noise. An estimate of the required actuation can be found here.

Attachment 1: IMG_8724.JPG
IMG_8724.JPG
Attachment 2: timeDomain.pdf
timeDomain.pdf
Attachment 3: HVampNoise.pdf
HVampNoise.pdf
  15536   Sun Aug 23 23:36:58 2020 gautamUpdateElectronicsFirst look at HV coil driver

Summary:

A more careful analysis has revealed some stability problems. I see oscillations at frequencies ranging from ~600kHz to ~1.5 MHz, depending on the voltage output requested, of ~2 V pp at the high-voltage output in a variety of different conditions (see details). My best guess for why this is happening is insufficient phase margin in the open-loop gain of the PA95 high voltage amplification stage, which causes oscillations to show up in the closed loop. I think we can fix the problem by using a larger compensation capacitor, but if anyone has a better suggestion, I'm happy to consider it

Details:

The changes I wanted to make to the measurement posted earlier in this thread were: (i) to measure the noise with a load resistor of 20 ohms (~OSEM coil resistance) connected, instead of the unloaded config previously used, and (ii) measure the voltage noise on the circuit side (= TP5 on the schematic) with some high voltage output being requested. The point was to simulate conditions closer to what this board will eventually be used in, when it has to meet the requirement of <1pA/rtHz current noise at 100 Hz. The voltage divider formed by the 25 kohm series resistor and the 20 ohm OSEM coil simulated resistance makes it hopeless to measure this level of voltage noise using the SR785. On the other hand, the high voltage would destroy the SR785 (rated for 30 V max input). So I made a little Pomona box to alllow me to do this measurement, see Attachment #1. Its transfer function was measured, and I confirmed that the DC high voltage was indeed blocked (using a Fluke DMM) and that the output of this box never exceeded ~1V, as dictated by the pair of diodes - all seemed okay .

Next, I wanted to measure the voltage noise with ~10mA current flowing through the output path - I don't expect to require more than this amount of current for our test masses. However, I noticed some strange features in the spectrum (viewed continuously on the SR785 using exponential averaging setting). Closer investigation using an oscilloscope revealed:

  1. 600kHz to 1 MHz oscillations visible, depending on output voltage.
  2. The oscillations vanish if I drive output above +30 V DC (so input voltage > 1 V).
  3. The oscillations seem to be always present when the output voltage is negative.
  4. No evidence of this offset if circuit is unloaded and voltage across 25k resistor is monitored. But they do show up on scope if connected to circuit side even in this unloaded config.

Some literature review suggested that the capacitor in the feedback path, C4 on the schematic, could be causing problems. Specifically, I think that having that capacitor in the feeddback path necessitates the use of a larger compensation capacitor than the nominal 33pF value (which itself is higher than the 4.7pF recommended on the datasheet, based on experience of the ESD driver circuit which this is based on, oscillations were seen there too but the topology is a bit different). As a first test of this idea, I removed the feedback capacitor, C4 - this seemed to do the trick, the oscillations vanished and I was able to drive the output between the high voltage supply rails. However, we cannot operate in this configuration because we need to roll off the noise gain for the input voltage noise of the PA95 (~6 nV/rtHz at 100 Hz will become ~200 nV/rtHz, which I confirmed using the SR785). Using a passive RC filter at the output of the PA95 (a.k.a. a "snubber" network) is not an option because we need to sum in the fast actuation path voltage at the output of the 25 kohm resistor.

Some modeling confirms this hypothesis, see Attachment #2.  The quantity plotted is the open-loop gain of the PA95 portion of the circuit. If the phase is 0 degrees, then the system goes unstable.

So my plan is to get some 470pF capacitors and test this idea out, unless anyone has better suggestions? I guess usually the OpAmps are compensated to be unconditionally stable, but in this case maybe the power op-amp is more volatile?

Quote:

Need to think more about how to better characterize this noise. An estimate of the required actuation can be found here.

Attachment 1: IMG_5379.JPG
IMG_5379.JPG
Attachment 2: stabilityCriterion.pdf
stabilityCriterion.pdf
  13274   Wed Aug 30 11:04:08 2017 GabrieleSummaryLSCFirst look at neural network reconstruction of PRMI motion

Introduction

I trained a deep neural network (DNN) to reconstruct MICH and PRCL degrees of freedom in the PRMI configuration. For details on the DNN architecture please refer to G1701455 or G1701589. Or if you really want all the details you can look at the code. I used the following signals as input to the DNN: POPDC, POP22_Q, ASDC, REFL11_I/Q, REFL55_I/Q, AS55_I/Q.

Gautam took some PRMI data in free swinging and driven configuration:

  • 1187819331 + 10mins: Free swinging PRMI (after first locking PRMI on carrier and dither aligning).
  • 1187820070 + 5mins: PRM driven at low freq.
  • 1187820446 + 5mins: BS driven at low freq.

In contrast to the Fabry-Perot cavity case, we don't have a direct measurement of the real PRCL/MICH degrees of freedom, so it's more difficult to assess if the DNN is working well.

Results

All MICH and PRCL values are wrapped into the unique region [-lambda/4, lambda/4]^2. It's even a bit more complicated than simpling wrapping. Indeed, MICH is periodic over [-lambda/2, lambda/2]. However, the Michelson interferometer reflectivity (as seen from PRC) in the first half of the segment is  the same as in the second half, except for a change in sign. This change of sign in Michelson reflectivity can be compensated by moving PRCL by lambda/4, thus generating a pi phase shift in the PRC round trip propagation that compensate for the MICH sign change. Therefore, the unit cell of unique values for all signals can be taken as [-lambda/4, lambda/4] x [-lambda/4, lambda/4] for MICH x PRCL. But when we hit the border of the MICH region, PRCL is also affected by addtion of lambda/4. Graphically, the square regions A B C below are all equivalent, as well as more that are not highlighted:

This makes it a bit hard to un-wrap the resonstructed signal, especially when you add in the factor that in the reconstruction the wrapping is "soft".

The plot below shows an example of the time domain reconstruction of MICH/PRCL during the free swinging period.

It's hard to tell if the positions look reasonable, with all the wrapping going on.

Two-dimensional maps of signals

Here's an attempt at validating the DNN reconstruction. Using the reconstructed MICH/PRCL signal, I can create a 2d map of the values of the optical signals. I binned the reconstructed MICH/PRCL in a 51x51 grid, and computed the mean value of all optical signals for each bin. The result is shown in the plot below, directly compared with the expectation from a simulation.

The power signals (POP_DC, AS_DC, PO22_Q) looks reasonably good. REFL11_I/Q also looks good (please note that due to an early mistake in my code, I reversed the convention for I/Q, so PRCL signal is maximized in Q instead than in I). The 55MHz signals look a bit less clear...

Steps forward

  • I'm quite confident in the tuning of demodulation phase and signs for REFL11 and POP22, but less so for REFL55 and not sure at all for AS55. So it would be useful to measure a full sensing matrix of PRCL and MICH against those signals, to compare with my simulation
  • I'm working on an idea to fine tune the DNN using the real interferometer data, more to follow when the idea crystallizes in a clear form.
  13674   Thu Mar 8 23:50:27 2018 gautamUpdateALSFirst look at new ALS electronics
  • Locked single arms, dither aligned, and saved offsets to EPICS (slow) sliders in anticipation of having to reboot all vertex FEs.
  • Shutdown ETMY watchdog, stopped all models on c1iscey, and shutdown that frontend.
  • Walked down to Y-end, powered of c1iscey expansion chassis, and removed the ADC adaptor card.
  • Stopped all models on c1lsc. Shutdown watchdog on all optics in anticipation of c1sus model failing. Shutdown the c1lsc frontend.
  • Powered off the c1lsc expansion chassis. Installed the borrowed adapter card from c1iscey in c1lsc expansion chassis. Connected it to the "spare" ADC card Koji and I had installed in c1lsc expansion chassis last Wednesday.
  • Connected differential output of demod board to differential input of AA chassis. Connected SCSI connector from output of AA chassis to the newly installed adapter card.
  • Powered the c1lsc expansion chassis back on. Then powered c1lsc FE on.
  • Walking back out to the control room, I saw that all vertex FEs had crashed. I had to go back in and hard-reboot c1sus.
  • Before bringing back any models, I backed up the existing c1lsc model, and then modified c1x04 and c1lsc to use the newly acquired ALS signals for the X arm ALS signal chain.
  • Restarted all vertex FE models. Everything came back up smooth. DC light is still red on c1oaf but I didn't bother trying to rectify it tonight for these tests.
  • Reset appropriate LSC offsets with PSL shutter closed. Locked X arm on IR. Reset phase tracker servo gain for X arm ALS. Engaged slow temperature servo on EX laser.

Then I looked at  the spectrum, see Attachment #1. Disappointingly, it looks like the arm PDH servo is dominating the noise, and NOT unsuppressed EX laser frequency noise,. Not sure why this is so, and I'm feeling too tired to debug this tonight. But encouragingly, the performance of the new ALS signal chain looks very promising. Once I tune up the X arm loop, I'm confident that the ALS noise will be at least as good as the reference trace.

I am leaving c1iscey shutdown until this is fixed. So ETMY is not available for the moment.

Random factoid: Trying to print a DTT trace with LaTeX in the label text on pianosa causes the DTT window to completely crash - so if you dont save the .xml file, you lose your measurement.

Quote:

I made a LISO fit of the measured TF of the daughter board, so that I can digitally invert the daughter board whitening. Results attached. (Inverse) Filters have been uploaded to the ALS X Foton filter banks.

 

Attachment 1: BeatMouth_OOL.pdf
BeatMouth_OOL.pdf
  13675   Fri Mar 9 01:07:01 2018 gautamUpdateALSFirst look at new ALS electronics

[koji, gautam]

I was going to head out but then it occurred to me that I could do another simple test, which is to try and lock the X arm on ALS error signal (i.e. actuate on MC length to keep the beat between EX laser and PSL fixed, while the EX frequency is following the Xarm length). Comparing the in loop (i.e. ALS) error signal with the out-of-loop sensor (i.e. POX), it seems like POX is noisy. The curves were lined up by eye, by scaling the blue curve to match the red at the ~16Hz peaks. This supports my hypothesis in the previous elog. On the downside, could be anything. Electronics in the POX chain? The demod unit itself? Will look into it more tomorrow..

As an aside, controlling the arm with ALS error signal worked quite well, and the lock was maintained for ~1 hour.

Attachment 1: ALS_vs_POXnoise.pdf
ALS_vs_POXnoise.pdf
  1714   Thu Jul 2 06:31:35 2009 ClaraUpdatePEMFirst mic in place and connected

I clamped Bonnie (microphone) to the top of a chamber near the vertex of the arms and placed Clyde (pre-amp) on the table right below (see picture). The cable was laid and Bonnie and Clyde are plugged into port #13 on the ADC. The second cable was plugged into port #14, but it is not connected to anything. I placed the looped up cable on top of the cabinet holding the ADC.

Note: the angle in the photograph is such that we are looking along the y-arm.


Attachment 1: bonnieandclyde.JPG
bonnieandclyde.JPG
  8874   Thu Jul 18 20:20:52 2013 gautamConfigurationendtable upgradeFirst mirror glued to PZT and mounted in modified mounts

 

 Yesterday, I mounted the first PZT in one of the modified mounts, and then glued a 2-inch Y2 mirror on it using superglue.

Details:

-The mirror is a 2-inch, Y2 mirror with HR and AR coatings for 532 nm light.

-The AR side of the mirror had someone's fingerprint on it, which I removed (under Manasa's guidance) using tweezers wrapped in lens cleaning paper, and methanol.

-Before gluing the mirror, I had to assemble the modified mount. Manasa handed over the remaining parts of the mounts (which are now in my newly acquired tupperware box along with all the other Piezo-related hardware). I took the one labelled A, and assembled the holder part. I then used one of the new mounts (2.5 inches, these are with the clean mounts in a cardboard box in the cupboard holding the green optics along the Y-arm) and mounted the holder on it. 

-Having assembled the mount, I inserted the piezo tip-tilt into the holder. The wedge that the machine shop supplied is useful (indeed required) for this. 

-I then cleaned the AR surface of the mirror and the top-surface of the tip-tilt. 

-The gluing was done using superglue which Steve got from the bookstore (the remaining tube is in the small fridge). We may glue the other mirror using epoxy. I placed 4 small drops of superglue on the tip-tilt's top surface, placed the mirror with its AR face in contact with the piezo, and applied some pressure for a short while until the glue spread out fairly evenly. I then left the whole setup to dry for about half an hour.

-Steve suggested using a reference piece (I used two small bolts) to verify when the glue had dried.

-Finally, I attached the whole assembly to a base.

Here it is in action in my calibration setup (note that it has not been oriented yet. i.e. the two perpendicular axes of the piezo are for the time being arbitrarily oriented. And maybe the spreading of the glue wasn't that even after all...):

Piezo-mirror.jpg

 

Sidenote:

Yesterday, while setting stuff up, I tested the piezo with a 0.05 Hz, 10Vpp input from the SR function generator just to see if it works, and also to verify that I had set up all my electronics correctly. Though the QPD was at this point calibrated, I did observe periodic motion of both the X and Y outputs of my QPD amp! Next step- calibration... 

  11979   Fri Feb 5 16:50:24 2016 gautamUpdateGreen LockingFirst pass at mode-matching

I've done a first pass at trying to arrive at a mode-matching solution for the X-end table once we swtich the lasers out. For this rough calculation, I used a la mode to match my seed beam (with z = 0 being defined as the shutter housing on the current position of the Innolight laser head, and the waist of the beam from the NPRO being taken as the square-root of the X and Y waists as calculated here), to a target beam which has a waist of 35um at the center of the doubling oven (a number I got from this elog). I also ignored the optical path length changes introduced by the 3 half-wave plates between the NPRO and the doubling oven, and also the Faraday isolator. The best a la mode was able to give me, with the only degrees of freedom being the position of the two lenses, was a waist of 41um at the doubling oven. I suppose this number will change once we take into account the effects of the HWPs and the Faraday. Moreover, the optimized solution involves the first lens after the NPRO, L1, being rather close to the second steering mirror, SM2 (see labels in Attachment #2, in cyan), but I believe this arrangement is possible without clipping the beam. Moreover, we have a little room to play with as far as the absolute physical position of the z=0 coordinate is - i.e. the Lightwave NPRO head can be moved ~2cm forward relative to where the Innolight laser head is presently, giving a slightly better match to the target waist (see attachment #3). I will check the lenses we have available at the 40m to see if a more optimal solution can be found, but I'm not sure how much we want to be changing optics considering all this is going to have to be re-done for the new end table... Mode-matching code in Attachment #4...

Attachment 1: Modematch_AUXx.pdf
Modematch_AUXx.pdf
Attachment 2: NewSetUp.png
NewSetUp.png
Attachment 3: Modematch_AUXx_2.pdf
Modematch_AUXx_2.pdf
Attachment 4: XendModeMatch.m.zip
  3593   Tue Sep 21 16:05:21 2010 josephbUpdateCDSFirst pass at rack diagram

I've made a first pass at a rack diagram for the 1X1 and 1X2 racks, attached as png.

Gray is old existing boards, power supplies etc.  Blue is new CDS computers and IO chassis, and gold is for the Alberto's new RF electronics.  I still need to double check on whether some of these boards will be coming out (perhaps the 2U FSS ref board?).

Attachment 1: 1X1_1X2_racks.png
1X1_1X2_racks.png
  5157   Tue Aug 9 16:21:59 2011 Manuel, IshwitaUpdateWienerFilteringFirst results of our simulations

We did the simulation of the stacks by defining a transfer function for one stack (green plot) and another similar transfer function for the other stack.

We simulated the ground motion by filtering a white noise with a low pass filter with a cutoff frequency at 10Hz. (blue plot) (the ground motion for the 2 stacks are completely uncorrelated)

We simulated the electronic white noise for the seismic measurements. (black plot)

We filtered the ground motion (without the measurements electronic noise) with the stack's transfer function and subtracted them to find the mirror response (red plot), which is the target signal for the wiener filter.

We computed the static wiener filter with the target signal (distance between the mirrors) and the input data (seismic measurements = ground motion + electronic noise).

We filtered the input and plotted the output (light blue plot).

We subtracted the target and the output to find the residual (magenta plot).

We didn't figure out why the residual is above the electronic noise only under ~6hz. We tried to increase and decrease the electronic noise and the residual follows the noise still only under ~6Hz.

It also shows that the residues are above the target at frequencies over 20Hz. This means that we are injecting noise here.

simu1.png

We tried to whiten the target and the input (using an high pass filter) to make the wiener filter to care even of higher frequencies.

The residues are more omogeneously following the target.

We also plotted the Wiener filter transfer function without making whitening and with making whitening. It shows that if we do whitening we inject no noise at high frequency. But we loose efficency at low frequencies.

simu2w.png

TFwhitened2.png

We shouldn't care about high frequency, because the seismometers response is not good over 50Hz. So, instead of whitening, we should simply apply a low pass filter to the filter output to do not inject noise and keep a good reduction at low frequencies.

 

  2941   Mon May 17 19:42:11 2010 JenneUpdateIOOFirst steps toward MC mode measuring

[Jenne, Kevin, Steve]

We made some progress toward getting the MC's beam profile measured.  In the end, no changes were made to anything today, but we're more prepared to go for tomorrow.

What we did:

* Grabbed the scanning slit beam scan from the PSL lab.  It's the same kind as we had here at the 40m, so Kevin was able to hook it up to the computer, and confirmed that it works.

* Opened the IOO and OMC chamber doors, and locked the MC.  Unfortunately the MC mode was awful in Yaw.  Awful like TEM(0,10+). But it still locked.  

* Confirmed that the beam went through the Faraday.  I looked at the beam before and after the Faraday on a card, and it was the same nasty beam both before and after.  So it looks like Zach did a good job aligning the Faraday and everything else.  I was going to clamp the Faraday, but I didn't yet, since I wanted to see the nice happy TEM00 mode go through without clipping before risking moving the Faraday during clamping (I don't know how heavy it is, so I'm not sure how much it might potentially move during clamping.)

* Noticed that there is a whole lot of crap on both the OMC and BS tables that's going to have to move.  In particular, one of the weights leveling the OMC table is right where I need to put MMT2.  Steve suggested putting the optic there, in its approximate place, before doing too much other stuff, since it could potentially affect the leveling of the table, and thus the input pointing to the MC.  Unfortunately, to do that I'll need to move the weight, which is definitely going to change things.  Sad face.  Moving the weight will likely be one of the first things I do tomorrow, so that all 3 profile measurements have the same configuration. 

* Before closing up, I tried to align the MC, to get back to TEM00, to no avail.  I got as far as achieving TEM11 flashing, along with a bunch of other crappy modes, but didn't get 00.  That's also on the to-do list.

What we're going to do:

* Open the chambers, and align the MC to TEM00 (using the sliders on the MC align screen).

* Check with an IR card that the beam goes through the Faraday.

* Clamp the Faraday, reconfirm.

* Remove the weight on the OMC table.

* Place MMT2 on the OMC table in it's approximate final location.

* Realign the MC, and make sure the beam goes through the Faraday.  If this doesn't happen smoothly, I may need more instruction since I've never dealt with aligning the Faraday before.  What are the appropriate mirrors to adjust? 

* Move the PZT flat steering mirror from the BS table to the IOO table.  (Thoughts on this?  This will change the table leveling, and also includes the trickiness of needing to move the connectors for the PZT.)

* Place a flat mirror on the BS table to route the MC beam out to the BS/PRM/SRM oplev table. 

* Measure the mode using the beam scan: on the BS oplev table, on the POX table, and then perhaps by shooting the beam through the beamtube on the ETMY (new convention) table.

* Place MMT1 on the BS table, use flat mirrors to get it out of the chambers, repeat measurements.

* Place MMT2 in the correct position, use flat mirrors to get it out of the chambers, repeat measurements.

All of this may require some serious cleaning-up of the BS table, which is going to be ugly, but it has to happen sometime. Hopefully I can get away with only moving a minimal number of things, in order to get these measurements done.

 

Another note: Don't trust the PSL shutter and the switch on the MEDM screens! Always use a manual block in addition!!! We discovered upon closeup that hitting the "Closed" button, while it reads back as if the shutter is closed (with the red box around the buttons), does not in fact close the shutter.  The shutter is still wide open.  This must be fixed.

  2942   Tue May 18 01:40:56 2010 KojiUpdateIOOFirst steps toward MC mode measuring

OK. Don't worry. This is just an initial confusion which we also had for the suspensions a while ago.

The faraday must be clamped. It shakes the table terribly but it is fine. The leveling may change a bit but should be small enough. Otherwise, just tweak the weights. In fact, the faraday has enough large apertures and we hope we don't need to move it again, as far as the MC incident beam is not moved. But if necessary, we don't move the mirrors but move the faraday itself.

Usually the alignment of the MC is taken by MC2/MC3 such that we don't  move the refl. But if you think what have moved is the MC1/MC3 (i.e. activity in the IMC chamber), take the alignment of the MC1/MC3.

It is just a matter of time to get TEM00. If you get TEM11, it is already close. If you align for TEM11, it is enough aligned to lock TEM10 or TEM01. Once you got better mode, align for it again. Eventually you will get TEM00.

The leveling may change by moving the optics and the weight again. But once the leveling is recovered by arranging the weights somewhere else,
the pointing must be fine again. If necessary, You can remove two optics for squeezing injection (strange motorized rotating mirror and a mount sticking out from the table to south.)

Yes, we need to move the PZT mirror. For the connection, only Steve can give us the right way to do it. If it is too much hussle, just move only the mirror and ignore the wiring for now.

I will update how the mirrors should be migrated from the table to the table.

 

  2945   Tue May 18 12:04:13 2010 robUpdateIOOFirst steps toward MC mode measuring

Quote:

Another note: Don't trust the PSL shutter and the switch on the MEDM screens! Always use a manual block in addition!!! We discovered upon closeup that hitting the "Closed" button, while it reads back as if the shutter is closed (with the red box around the buttons), does not in fact close the shutter.  The shutter is still wide open.  This must be fixed.

 Has anyone tried pushing the "reset" button on the Uniblitz driver?

ELOG V3.1.3-