ID |
Date |
Author |
Type |
Category |
Subject |
742
|
Sat Jul 26 15:09:57 2008 |
Aidan | Update | Computers | Reboot of op440m |
I was reviewing the PSL Overview screen this afternoon and op440m completely froze when I center-clicked on the REF CAVITY TRANSMISSION indicator. It was unresponsive to any keyboard or mouse control. The moon button had no effect to shut the machine down.
Called Alberto in and we logged into op440m from rosalba. From there we logged in as 'root' and run a shutdown script '/usr/sbin/shutdown -i S -g 1'. The medm screens started disappearing from the op440m display and we were eventually asked to enter System Maintenance Mode. From here we selected RUN LEVEL 5: "state 5: Shut the machine down so that it is safe to remove the power". Following this the machine turned itself off.
We powered it back on, logged back in as controls and restarted the medm screens. Everything seems to be running fine now.
Aidan. |
1205
|
Mon Dec 29 18:01:07 2008 |
Aidan | Update | Auxiliary locking | Updated 40m Upgrade Document T080074-00-R |
Added a paragraph to the 40m Upgrade document describing the fiber stabilization and frequency doubling proposed for auxiliary locking.
Also added a complete diagram of the fiber stabilization and a draft sketch of the frequency doubling.
Uploaded to https://nodus.ligo.caltech.edu:30889/svn/trunk/docs/upgrade08/ via svn. |
1687
|
Fri Jun 19 13:39:29 2009 |
Aidan | Update | General | Upper limit measurement of the scatter from the eLIGO beam dumps. |
I measured the scatter from the eLIGO beam dumps as best I could. The experiment setup is shown in the attached diagram.
After familiarizing myself with the equipment in the morning I noticed three issues with the setup
1 - around the minimum scatter the back scatter from the beam dump is very susceptible to the incident angle (makes sense since the Si plate inside the beam dump at Brewster's angle when there is minimum scatter).
2 - The mirrored plug (Part 20 in D0900095) which is suppose to be used for alignment is not very effective. It moves around too much in its hole in the front face of the beam dump. Just by touching it I could make the reflected beam jump around by about 0.1 radians.
- I think to align these properly we'll have to partly assemble the dumps. If we leave off the front plate of the horn then we can measure the reflection off the Si. If we measure this with a power meter then alignment becomes a simple matter of rotating until this reflection is minimized.
3. - For this measurement the incident beam was a small (~ 1mm diameter) central beam with a small amount of spray of laser light beyond that central region. This spray was hitting the aluminium front face of the beam dump and was scattering back to the photodiode. This was clearly the limiting factor in the measurement. Most of this light was spread horizontally so I placed a couple of pieces of black glass on either side of the aperture, just blocking the edges a little. This reduce the background reading at the minimum scatter from 17.0uV to around 4.5uV with still a little bit of light hitting the top and bottom of beam dump face.
The incident power on the beam dump fluctuated a little but was in the range 20.5 to 22mW. The response of the PD is approximately 0.2 A/W and the transimpedance is 7.5E4 V/A.
The SR830 Sensitivity was set to 1x1 mV.
It was difficult to measure the actual angle of incidence. The dump pivoted about a point directly under the input aperture at the front. By measuring the displacement of a point on the back of the dump as I rotated it and knowing the distance between this point and the pivot point I was able to make a reasonably accurate measurement of a range of angles about the minimum.
The measured scatter (in V measured directly by the PD and as a fraction of the incident power) is shown in the attached plots.
I think I can do a better job cleaning up the incident beam - so these numbers only represent an upper limit on the scatter.
attachment 1: beam dump assembly
attachment 2: experimental layout
attachment 3: scatter measurement
attachment 4: BRDF - (scatter divided by the solid angle = 1.1 m steradians)
attachment 5: (slightly blurred )photo of dump - overhead view |
Attachment 1: D0900095_ELIGO_35_WATT_AIR_COOLED_BEAM_DUMP_3INCH_P.POL-3.PDF
|
|
Attachment 2: beam_dump_expt.png
|
|
Attachment 3: scatter_measurement1.pdf
|
|
Attachment 4: BRDF1.pdf
|
|
Attachment 5: IMG_0308.JPG
|
|
1782
|
Thu Jul 23 07:34:45 2009 |
Aidan | Update | CDS | Added C2 MEDM screens to 40m SVN. |
See Adhikari eLOG entry: http://nodus.ligo.caltech.edu:8080/AdhikariLab/194 |
2792
|
Mon Apr 12 17:48:32 2010 |
Aidan | Update | Computer Scripts / Programs | elog restarted |
The elog crashed when I was uploading a photo just now. I logged into nodus and restarted it. |
2793
|
Mon Apr 12 19:50:30 2010 |
Aidan | Summary | Green Locking | Temperature sweep of the Lightwave: df/dT = 2.8GHz/K |
The beams from the Innolight and Lightwave NPROs were both incident on a 1GHZ New Focus PD. Mott and I swept the temperature of the Lightwave and tracked the change in frequency of the beatnote between the two. The Innolight temperature was set to 39.61C although the actual temperature was reported to be 39.62C.
Freq. vs temperature is plotted below in the attached PDF. The slope is 2.8GHz/K.
The data is in the attached MATLAB file. |
Attachment 1: LightWave_temp_sweep.pdf
|
|
Attachment 2: LightWave_Temp.m
|
% plot the data from the Lightwave Temperature sweep
% Lightwave temperature
LWTemp = [0.2744
0.2753
.2767
.2780
.2794
.2808
... 67 more lines ...
|
2807
|
Mon Apr 19 11:31:04 2010 |
Aidan | Update | Green Locking | 1W NPRO output profile |
Quote: |
Koji asked me to take a profile of the output of the 1W NPRO that will be used for green locking. I used the razor-scan method, plotting the voltage output of a PD vs the position of the razor across the beam, both vertically and horizontally. This was done at 6 points along the beam path out of the laser box.
I determined the beam spot size at each point by doing a least-squares fit on the plots above in Matlab (using w as one of the fitting parameters) to the cumulative distribution functions (error functions) they should approximate.
I then did another least-squares fit, fitting the above "measured" beam profiles to the gaussian form for w vs z. Below is a summary.
It seems reasonable, though I know that M2 < 1 is fishy, as it implies less divergence than ideal for that waist size. Also, like Koji feared, the waist is inside the box and thus the scan is almost entirely in the linear regime.

|
There is a clearly a difference in the divergence angle of the x and y beams - maybe 10-20%. Since the measurements are outside the Rayleigh range and approximately in the linear regime, the slope of the divergence in this plot should be inversely proportional to the waists - meaning the x- and y- waist sizes should differ by about 10-20%. You should check your fitting program for the waist.
|
2809
|
Mon Apr 19 16:27:13 2010 |
Aidan | Update | Green Locking | Raicol crystals arrived and we investigated them |
Jenne, Koji and I opened up the package from Raicol and examined the crystals under the microscope. The results were mixed and are summarized below. There are quite a few scratches and there is residue on some of the polished sides. There is a large chip in one and there appear to be gaps or bands in the AR coatings on the sides.
There are two albums on Picassa
1. The package is opened ...
2. The crystals under the microscope.
Crystal |
Summary |
724 |
Chip in the corner of one end face, Otherwise end faces look clean. Large scratch on one polished side. |
725 |
End faces look good. Moderate scratch on one polished face. Residue on one polished face. |
726 |
Tiny dot on one end face, otherwise look okay. Large bands in one polished face. Moderate scratch on polished face |
727 |
Large, but shallow chip on one polished face. End faces look clean. Bands in one of the polished faces. |
|
2816
|
Tue Apr 20 11:14:31 2010 |
Aidan | Update | Green Locking | Raicol crystals arrived and we investigated them |
Here is Crystal 724 polished side 2 with all photos along the length stitched together |
2910
|
Tue May 11 14:39:17 2010 |
Aidan | Update | Green Locking | Green Laser Beam Profile |
Here's a photo of the set-up used. The beam profile is measured relative to the f=-100mm lens. |
Attachment 1: P5110057_beams.jpg
|
|
2951
|
Wed May 19 14:36:46 2010 |
Aidan | HowTo | Phase Camera | Phase Camera algorithm and stuff |
I had a think about the algorithm we might use for the phase camera measurement. MATLAB has an fft function that will allow us to extract the data that we need with a single command.
We record a series of images from a camera and put them into a 3D array or movie, image_arr, where the array parameters are [x-position, y-position, time], i.e. a 2D slice is a single frame from the camera. Then we can do an FFT on that object with the syntax, f3D = fft(image_arr, [ ], 3), which only does the FFT on the temporal components. The resulting object is a 3D array where each 2D slice is an 2D array of amplitude and phase information across the image for a single temporal frequency of the movie.
So if we recorded a movie for 1s where the sample rate is 58Hz, then the 1st frame of f3D is just a DC image of the movie, the 2nd frame are the complex 1Hz components of the movie, etc all the way up to 29Hz.
Suppose then that we have a image, part of which is being modulated, e.g. a chopper wheel rotating at 20 or 24Hz, or a laser beam profile which contains a 1kHz beat between a sideband and a reference beam. All we have to do is sample at at least twice that modulation frequency, run the command in MATLAB, and then we immediately get an image which contains the phase and magnitude information that we're interested in (in the appropriate 2D slice o the FFT).
As an example, I recorded 58 frames of data from a camera, sampling at 58Hz, which was looking at a spinning chopper wheel. There was a white sheet of paper behind the wheel which was illuminated from behind by a flashlight. The outer ring was chopping at 24Hz and the inner ring was chopping at 20Hz. I stuck all the images into the 3D array in MATLAB, did the transformation and picked out the DC, 20Hz and 24Hz signals. The results are shown in the attached PDFs which are:
- phase_camera_DC_comp.pdf - a single image from the camera and the DC component (zoomed in) of the FFT
- phase_camera_F1_comp.pdf - the magnitude and phase information of the 20Hz component of the FFT
- phase_camera_F2_comp.pdf - the magnitude and phase information of the 24Hz component of the FFT (this PDF contains a typo that says 25Hz).
- load_raw_data.m - the MATLAB routine that loads the saved data from the camera and does the FFT
You can, and I have, run the MATLAB engine from C directly. This will allow you to transfer the data from the camera to MATLAB directly in memory, rather than via the disk, but it does need proper memory allocation to avoid segmentation faults - that was too frustrating for me in the short term. In this case, the 58 frames were recorded to a file as a contiguous block of data which I then loaded into MATLAB, so it was slower than it might've otherwise been. Also the computer I was running this on was a bit of a clunker so it took a bit of time to do the FFT.
The data rate from the camera was 58fps x (1024 x 1024) pixels per frame x 2 bytes per pixel = 116MB per second. If we were to use this technique in a LIGO phase camera, where we want to measure a modulation which is around 1kHz, then we'd need a sample rate of at least 2kHz, so we're looking at at least a 30x reduction in the resolution. This is okay though - the original phase camera had only ~4000 spatial samples. So we could use, for instance, the Dalsa Falcon VGA300 HG which can give 2000 frames per second when the region of interest is limited to 64 pixels high. |
Attachment 1: phase_camera_DC_comp.pdf
|
|
Attachment 2: phase_camera_F1_comp.pdf
|
|
Attachment 3: phase_camera_F2_comp.pdf
|
|
Attachment 4: load_raw_data.m
|
% load a raw data file into MATLAB
fid = fopen('phase_camera_data.dat');
n1 = 750;
A3D = ones(n1, n1, 58);
for jj = 1:58
A = fread(fid, [1024, 1024], 'uint16');
A3D(:,:,jj) = A((512-floor(n1/2)):(512-floor(n1/2))+n1-1, ...
(512-floor(n1/2)):(512-floor(n1/2))+n1-1);
... 64 more lines ...
|
2955
|
Thu May 20 10:06:56 2010 |
Aidan | HowTo | Phase Camera | Phase Camera- raw data video |
|
2995
|
Wed May 26 18:54:55 2010 |
Aidan | Summary | Green Locking | Mounted Crystal 724 in the Doubling Oven |
Andri and I mounted the Raicol Crystal #724 in one of the new Covesion Ovens. The procedure was the same as before - see elog entry here.
There was one issue - the glass plate that goes on top of the crystal is coated on one side with ITO (Indium-Tin Oxide) and it's not 100% certain that this was mounted in the correct orientation. It is virtually impossible to tell which side of the glass is coated.
The base plate of the oven was tapped for an M3 hole. We retapped it for an 8-32 and bolted it to a post and that one of the New Focus 4-axis translation stage. The assembly is currently bolted to the PSL table, awaiting use. |
2998
|
Thu May 27 08:22:57 2010 |
Aidan | Update | Computers | Restarted the elog this morning |
|
3082
|
Wed Jun 16 18:14:13 2010 |
Aidan | Update | General | Glass cover from overhead light smashed on PSL table |
I was giving a tour of the 40m yesterday. We were looking at the PSL table. About 30 seconds after I turned the lights on a glass cover from one of the lights (NW corner) popped out of its holder and smashed on the table.
I've cleaned up all the broken glass I could see but there may be some small shards there. Please use caution in that area. |
Attachment 1: DSC_1769.JPG
|
|
Attachment 2: DSC_1768.JPG
|
|
3113
|
Thu Jun 24 06:49:29 2010 |
Aidan | Update | Green Locking | a channel for PPKTP temperature |
Is this a setpoint temperature that we can change by writing to the channel or is it a readout of the actual temperature of the oven?
kiwamu:
This is a readout channel just to monitor the actual temperature.
Quote: |
We added a channel on c1psl in order to monitor the temperature of the PPKTP sitting on the PSL table.
|
|
3421
|
Fri Aug 13 15:29:35 2010 |
Aidan | Frogs | Photos | Here's the 40m team |
|
Attachment 1: 40m_team.JPG
|
|
3756
|
Thu Oct 21 19:10:39 2010 |
Aidan | Update | PSL | Found the beat at 1064nm |
Quote: |
[Koji Suresh]
We found the beat at 1064nm. T(PSL)=26.59deg, T(X-end)=31.15deg.
The X-end laser was moved to the PSL table.
The beating setup was quickly constructed with mode matching based on beam profile measurements by the IR cards.
We used the 1GHz BW PD, Newfocus #1611-FS-AC.
As soon as we swept the Xtal temp of the X-end laser, we found the strong beating.
|
|
3766
|
Fri Oct 22 21:28:27 2010 |
Aidan | Update | Locking | Green era: found a green beat note finally |
Nice work!
Quote: |
finally we found it !

|
|
3907
|
Fri Nov 12 11:36:06 2010 |
Aidan | Configuration | Computers | Fixed the PERL PATH errors from the scripts directory change |
I've been trying to get a PID loop running for the green locking and I've discovered that some of the directories in the path for the Perl Modules are now obselete. Specifically, since the scripts directory has been moved from /cvs/cds/caltech/scripts to /cvs/cds/rtcds/caltech/c1/scripts the following locations in the Perl @INC path list need to be changed:
/cvs/cds/caltech/scripts/general
/cvs/cds/caltech/scripts/general/perlmodules
I've added the above directories to the PERL5LIB path in /cvs/cds/caltech/cshrc.40m.
setenv PERL5LIB /cvs/cds/caltech/scripts/general:
/cvs/cds/caltech/scripts/general/perlmodules:
/cvs/cds/caltech/libs/solaris9/usr_local_lib/perl5/5.8.0/:
/cvs/cds/rtcds/caltech/c1/scripts/general:
/cvs/cds/rtcds/caltech/c1/scripts/general/perlmodules
This seems to fix the problem .. at least, you no longer get an error if you try nodus:~>perl -e 'user EpicsTools'
|
3908
|
Fri Nov 12 12:06:11 2010 |
Aidan | Configuration | Green Locking | PID script working - now it needs to be tuned |
I've set up a PID script that senses the EX-PSL Green Beat note (from the frequency counter) and actuates on the temperature of the end laser to drive the beat note to a given setpoint.
- I've added the necessary EPICS channels to
c1iscaux and rebooted it so that the channels are live. They are listed in a new database file slow_grnx_pid.db
- This database was added to the list of those loaded by
startup.cmd .
- The PID script,
GRNXSlowServo , is just a modified version of FSSSlowServo .
- The version I've been running is currently in /cvs/cds/caltech/users/abrooks/.
- There's also an MEDM screen in this directoy,
C1LSC_EX_GRN_SLOW.adl , there that shows the PID settings.
Right now the script only passes the initial sanities checks, that is:
- It runs.
- You can enable/disable it without any errors and it starts actuating.
The settings all need to be tuned up - e.g. maximum_increment, hard_stops, time_step, PID constants.
Additionally, the units in the whole thing are pretty useless - some of the channels are in VOLTS, others in WATTS. I'd like to change all these to be in Hz.
EPICS channels added:
grecord(ao,"C1:LSC-EX_GRN_SLOWKD")
grecord(ao,"C1:LSC-EX_GRN_SLOWKP")
- grecord(ao,"C1:LSC-EX_GRN_SLOWKI")
grecord(ao,"C1:LSC-EX_GRN_PID_SETPT")
grecord(ao,"C1:LSC-EX_GRN_TIMEOUT")
grecord(stringin,"C1:LSC-EX_GRN_SLOWVERSION")
grecord(bi,"C1:LSC-EX_GRN_SLOWLOOP")
grecord(bi,"C1:LSC-EX_GRN_DEBUG")
grecord(bi,"C1:LSC-EX_GRN_SLOWBEAT")
|
Attachment 1: Screen_shot_2010-11-12_at_12.05.01_PM.png
|
|
3922
|
Mon Nov 15 14:42:01 2010 |
Aidan | Update | PSL | C1PSL rebooted? |
Yeah. Joe and I rebooted c1psl a couple of times this morning. I didn't realize the burtrestore wasn't automatic.
Quote: |
Has C1PSL rebooted? Has burtrestore been forgotten? Even without elog?
We found some settings are wrong and the PMC has pretty low gain.
|
|
3930
|
Tue Nov 16 09:02:54 2010 |
Aidan | Update | Green Locking | PID loop - calibration of SR620 output |
[Aidan, Kiwamu]
Kiwamu and I roughly calibrated the analogue output from the SR620 frequency counter yesterday. The input channel, intuitively named C1:PSL-126MOPA_126MON, now reads the measured frequency in MHz with an error of about 0.1MHz - this is, I think, due to the bit noise on the D/A conversion that Kiwamu discovered earlier. That is, the output range of the SR620 corresponds to around 100MHz and is digitized at 10-bit resolution, and ...
100MHz/(10^2) ~= 0.098MHz. [Sad Face]
Calibration:
We set the EPICS range to [-100, 100] (corresponding to [-5V, 5V]), connected a Marconi to the Freq Counter, input a variety of different frequencies and measured the counts on the EPICS channel.
The linear fit to the calibration data was F = 2.006*EPICScount - 0.2942. From this we worked out the maximum and the minimum for the range settings that give the channel in MHz: EGUF = -200.8942 and EGUL = 200.3058. The previous range was [-410, 410]
Calibration of SR620 analogue output
Input Frequency (MHz) |
Measured EPICS Value |
10 |
5.191 |
20 |
9.98 |
30 |
15.21 |
40 |
20.00 |
50 |
25.18 |
60 |
29.99 |
70 |
35.18 |
71 |
35.565 |
72 |
35.9894 |
73 |
36.3861 |
74 |
37.17 |
75 |
37.576 |
76 |
37.9669 |
77 |
38.3575 |
78 |
39.166 |
79 |
39.5691 |
80 |
39.978 |
|
|
|
3931
|
Tue Nov 16 10:47:45 2010 |
Aidan | Update | Green Locking | Rebooted c1psl - added new GRNBEAT_FREQ channel |
I restored C1:PSL-126MOPA_126MON to its original settings (EGUF = -410, EGUL = 410) and added a new calc channel called C1:LSC-EX_GRNBEAT_FREQ that is derived from C1:PSL-126MOPA_126MON. The calibration in the new channel converts the input to MHz.
grecord(calc, "C1:LSC-EX_GRNBEAT_FREQ")
{
field(DESC,"EX-PSL Green Beat Note Frequency")
field(SCAN, ".1 second")
field(INPA,"C1:PSL-126MOPA_126MON")
field(PREC,"4")
field(CALC,"0.4878*A")
}
I rebooted c1psl and burtrestored. |
3932
|
Tue Nov 16 12:47:30 2010 |
Aidan | Update | Green Locking | PID loop but no green |
The PID loop is ready to be run on the green beat note but, since the tanks are open, there is no green transmission from the end getting to the PSL table. Nevertheless, here's the screen for the PID loop. The loop script is still in my directory /cvs/cds/caltech/users/abrooks/GRNXSlowServo
The medm screen is attached. It shows the current beat note frequency in MHz ()
In c1auxey/ETMYaux.db I added a couple of channels. These are all displayed on the MEDM screen. I added them to autoBurt.req as well.
- C1:LSC-EX_GRNLSR_TEMP_NOM: the zero-volt setpoint temperature of the end laser (as set on the front panel of the Mephisto controller). This must be entered manually in EPICS as there is no way to read it remotely. [Sad Face]
- C1:LSC-EX_GRNLSR_TEMP_CALC: the sum of the zero-volt set point temperature and the offset temperature set by the input voltage from C1:LSC-EX_GREENLASER_TEMP
I rebooted c1auxey to get these to work.
Once we get the green beat back again, the PID loop should servo on the end laser temperature to drive the Beat Frequency to the Frequency Setpoint, C1:LSC-EX_GRN_PID_SETPT, which can be set by the pink slider.
RA: All MEDM screens must be in the proper MEDM directory!! Also, all perl scripts must have a .pl extension!!! Also, all scripts must be in the scripts directory even if they are in development!!! And all scripts should use 'env' rather than have absolute pathnames for the location of perl, csh, tcsh, python, etc. |
Attachment 1: Screenshot-C1LSC_EX_GRN_SLOW.adl.png
|
|
3934
|
Tue Nov 16 16:00:26 2010 |
Aidan | Update | Green Locking | PID loop but no green |
Quote: |
RA: All MEDM screens must be in the proper MEDM directory!! Also, all perl scripts must have a .pl extension!!! Also, all scripts must be in the scripts directory even if they are in development!!! And all scripts should use 'env' rather than have absolute pathnames for the location of perl, csh, tcsh, python, etc.
|
That's not unreasonable. But if we try
grep "perl" /cvs/cds/rtcds/caltech/c1/scripts/*/*
you can see that we've got a fair amount of housekeeping to attend to. We might want to think about tidying up the scripts directory as part of the cds upgrade.
|
4040
|
Fri Dec 10 09:58:57 2010 |
Aidan | Update | SUS | beam pointing has been done |
Good news. I feel multi-chromatic-locking success is just around the corner.
By the way, there's a new presentation on the DCC from the ANU group where they've locked a short single cavity with both colors - G1000735:
https://dcc.ligo.org/cgi-bin/private/DocDB/ShowDocument?docid=14040
Quote: |
[Koji, Osamu and Kiwamu]
We aligned the beam axis pointing down to both X and Y arm.
Now the beams are hitting the centers of both ETMX and ETMY.
Amazingly Osamu made X arm flashing by aligning the cavity.
|
|
4177
|
Thu Jan 20 15:39:59 2011 |
Aidan | Update | Locking | Upper limit on frequency noise of ADC |
[Aidan, Kiwamu]
Kiwamu and I plugged the output from a DS3456 function generator into the ADC and started recording the data. The func. generator output a 237.8Hz, 1Vpp sine wave. We chose this value because it corresponds to the FSR of a 38.5m cavity (=3.896MHz) divided by 2^14, the frequency divider amount we intend to use.
Since 1 FSR divided down is 237.8Hz and corresponds to a length change of the cavity of 532nm/2 = 266nm, then we can, roughly, say that a frequency change of 1Hz corresponds to a length change in the cavity of approximately 1nm. The width of the 237.8Hz peak in the spectra corresponds to an upper limit on the noise floor due to digitizing the signal (this could be limited by the ADC, or the function generator, or the windowing on the FFT).
The FWHM of the peak in the spectrum was approximately 5mHz, corresponding to an uncertainty in the length of the cavity of about 6pm (we used a Hanning Window, 50% overlap and a BW of 2.92mHz, 7 averages). Regardless of what is the dominant contribution to the width of the peak, this implies that the frequency noise associated with digitizing a signal in the ADC is much smaller than we require and will not limit our performance if we choose to use a frequency divider and digital PLL with the Green Locking.
RA: Here's the previous measurement |
Attachment 1: Sine-wave-width-test_of_ADC.pdf
|
|
4178
|
Thu Jan 20 17:00:39 2011 |
Aidan | Configuration | Locking | Ballpark figures for Green Locking PLLs (Digital vs Analogue) |
If we use a digital PLL for locking the frequency of the PSL and END green lasers then we can expect a UGF of around 1kHz (assuming a sampling rate of 16kHz). Let's assume a simple 1/f loop giving a loop gain of ~1000x at 1Hz. If the free-swinging ETM pendulum motion at 1Hz is of the order of 1 micron, then the residual motion at 1Hz, once we lock the digital PLL by actuating on the ETM position, will be of the order of 1nm. This is bordering on too high.
Alternatively, is we use an analogue PLL then we can expect a much higher UGF and many orders of magnitude more gain at 1Hz (see here). So we would expect the residual motion of the pendulum to be much smaller - probably limited by some other noise source somewhere in the system (I doubt it's going to be reduced by 12 orders of magnitude).
RA: I think ballpark's not good enough for this. To see what's good enough, we need to to an analysis similar to what Bram has for the ALS. Get the 40m seismic spectrum from the arm locking spectrum or the green laser feedback signal and then correct it for a realistic loop shape.
KA: For this purpose I have made the simulink model for the green locking more than a year ago, but the entire green team has consistently neglected its presence...
https://nodus.ligo.caltech.edu:30889/svn/trunk/docs/upgrade08/Green_Locking/Servo_modeling/091121/ |
4182
|
Fri Jan 21 11:45:01 2011 |
Aidan | Configuration | Locking | Ballpark figures for Green Locking PLLs (Digital vs Analogue) |
Quote: |
If we use a digital PLL for locking the frequency of the PSL and END green lasers then we can expect a UGF of around 1kHz (assuming a sampling rate of 16kHz). Let's assume a simple 1/f loop giving a loop gain of ~1000x at 1Hz. If the free-swinging ETM pendulum motion at 1Hz is of the order of 1 micron, then the residual motion at 1Hz, once we lock the digital PLL by actuating on the ETM position, will be of the order of 1nm. This is bordering on too high.
Alternatively, is we use an analogue PLL then we can expect a much higher UGF and many orders of magnitude more gain at 1Hz (see here). So we would expect the residual motion of the pendulum to be much smaller - probably limited by some other noise source somewhere in the system (I doubt it's going to be reduced by 12 orders of magnitude).
RA: I think ballpark's not good enough for this. To see what's good enough, we need to to an analysis similar to what Bram has for the ALS. Get the 40m seismic spectrum from the arm locking spectrum or the green laser feedback signal and then correct it for a realistic loop shape.
KA: For this purpose I have made the simulink model for the green locking more than a year ago, but the entire green team has consistently neglected its presence...
https://nodus.ligo.caltech.edu:30889/svn/trunk/docs/upgrade08/Green_Locking/Servo_modeling/091121/
|
Agreed. It doesn't completely rule out the digital PLL. I'll check out Kiwamu's model. |
4205
|
Wed Jan 26 10:11:47 2011 |
Aidan | Update | Green Locking | cavity scan |
Quote: |

|
Whether or not it's as clean as we'd like, it's really nice to see this result with real data. |
4209
|
Wed Jan 26 14:49:48 2011 |
Aidan | Update | Environment | Turned on Control room AC |
80 degrees is too uncomfortable in the control room so I turned on the AC. The set point is 74F. |
4217
|
Fri Jan 28 09:03:38 2011 |
Aidan | Summary | Green Locking | Digital Frequency Discriminator |
Quote: |
That's some pretty fast work! I thought we would be taking up to a week to get that happening. I wonder what's the right way to measure the inherent frequency noise of this thing?
Also, should the comparator part have some hysteresis (ala Schmidt trigger) or is it best to just let it twirl as is? Is it sensitive to DC offsets on the input or is there a high pass filter? What's the correct low pass filter to use here so that we can have a low phase lag feedback to the ETM?
|
We could try inputing a 4kHz carrier modulated width a depth of a few Hz at a modulation frequency of F1. Then we could take an FFT of the output of the discriminator and measure the width of the peak at F1 Hz. This seems like an arduous way to measure the frequency noise at a single frequency though.
It'll definitely be sensitive to DC offsets but there is already a filter bank on the INPUT filter so we can shape that as necessary. We could probably band-pass that from [4.5 - 5.3kHz] (which would correspond to a range of [73,87] MHz into a 2^14 frequency divider.
|
4227
|
Sun Jan 30 17:15:09 2011 |
Aidan | Summary | Green Locking | Digital Frequency discriminator - frequency noise |
I've had a go at trying to estimate the frequency noise of the digital frequency discriminator (DFD). I input a 234.5Hz (0.5Vpp) signal from a 30MHz function generator into the ADC. The LP output of the DFD measured 234.5Hz. However, this signal is clearly modulated by roughly +/- 0.2Hz at harmonics of 234.5Hz (as you can see in the top plot in the dataviewer screenshot below). So the frequency noise can be estimated as rms of approximately 0.2Hz.
This is supported by taking the spectra of the LP output and looking at the RMS. Most of the power in the RMS frequency noise (above the minimum frequency) comes from the harmonics of the input signal and the RMS is approximately 0.2Hz.
I believe this stems from the rather basic LP filter (three or four poles around 10Hz?) that is used in the LP filter to remove the higher frequency components that exist after the mixing stage. (The currently loaded LPF filter is not the same as the saved one in Foton - and that one won't load at the moment, so I'm forced to remember the shape of the current filter).
The attached screen capture from data viewer shows the LP_OUT hovering around 234.5Hz. |
Attachment 1: _Untitled_(modified).png
|
|
4229
|
Mon Jan 31 07:03:59 2011 |
Aidan | Summary | Green Locking | DFD - noise spectra |
Quote: |
I've had a go at trying to estimate the frequency noise of the digital frequency discriminator (DFD). I input a 234.5Hz (0.5Vpp) signal from a 30MHz function generator into the ADC. The LP output of the DFD measured 234.5Hz. However, this signal is clearly modulated by roughly +/- 0.2Hz at harmonics of 234.5Hz (as you can see in the top plot in the dataviewer screenshot below). So the frequency noise can be estimated as rms of approximately 0.2Hz.
This is supported by taking the spectra of the LP output and looking at the RMS. Most of the power in the RMS frequency noise (above the minimum frequency) comes from the harmonics of the input signal and the RMS is approximately 0.2Hz.
I believe this stems from the rather basic LP filter (three or four poles around 10Hz?) that is used in the LP filter to remove the higher frequency components that exist after the mixing stage. (The currently loaded LPF filter is not the same as the saved one in Foton - and that one won't load at the moment, so I'm forced to remember the shape of the current filter).
The attached screen capture from data viewer shows the LP_OUT hovering around 234.5Hz.
|
Here is the spectrum of the input into the DFD (a 234.5Hz sine wave, 0.5 Vpp) and the spectrum and RMS of the LP output. The linewidth of the input signal is clearly much less than 0.1Hz, where as the RMS noise (above 2mHz) is approximately 0.2Hz and the main contributions are clearly the harmonics of the 234.5Hz signal. |
Attachment 1: DFD-bandwidth_noise.pdf
|
|
4230
|
Mon Jan 31 07:41:23 2011 |
Aidan | Update | Green Locking | DFD - medm screen |
I added an MEDM screen for the DFD to the GREEN screen. It is displayed in the attached screen shot.
This screen is located in: /cvs/cds/rtcds/caltech/c1/medm/c1gfd/C1GFD_DFD.adl |
Attachment 1: Screenshot-3.png
|
|
4234
|
Mon Jan 31 18:25:25 2011 |
Aidan | Update | Green Locking | DFD - results from the new filters (and running with AWG) |
Quote: |
This is a plot showing the old filters and the new ones we added this morning.
The new ones have a Cheby for AC coupling below 10 Hz and then a 500 Hz LP after the mixer. The LP frequency has been increased so that we can use this signal in a feedback loop to the ETM with a ~100 Hz UGF.
|
Joe injected a 234.567 etc. Hz sine wave into the excitation channel in the DFD INPUT filter. The spectrum of the output of the LP filter with the new filter is shown below with the RMS calculated from 300Hz down to 1mHz - see first attachment. The RMS is equal to about 2.5Hz. (Incidentally, the RMS is very much higher (slightly less than 400Hz - see second attachment) if you calculate it from 7kHz down to 1mHz). |
Attachment 1: DFD-bandwidth_noise_CLP500_to_300Hz.pdf
|
|
Attachment 2: DFD-bandwidth_noise_CLP500_to_7000kHz.pdf
|
|
4245
|
Thu Feb 3 16:08:06 2011 |
Aidan | Update | Computer Scripts / Programs | RCG VCO frequency error |
Joe and I were looking at the RCG VCO algorithm to determine if we could adapt it to run at a faster rate (you can currently change its frequency at 1Hz). I noticed that the algorithm that is used to calculated the values of sine and cosine at time T1 is a truncated Taylor series which uses the values of sine and cosine calculated at time T1 - Delta t . I was concerned that there would be an accumulating phase error so I tested the algorithm in MATLAB and compared it to a proper calculation of sine and cosine. It turns out that at a given 'requested' frequency there is a constantly accumulating phase error - which means that the 'actual' frequency of the RCG VCO is incorrect. So I have plotted the frequency error vs requested VCO frequency. It gets pretty bad!
Here's the code I used:
dt = 1/16384;
diffList = [];
% set the frequencies
flist = 1:5:8192;
for f = flist;
% get the 'accurate' values of sine and cosine
tmax = 0.05;
time1 = dt:dt:tmax;
sineT = sin(2.0*pi*f*time1);
cosineT = cos(2.0*pi*f*time1);
% determine the phase change per cycle
dphi = f*dt*2*pi;
cosT1 = 1:numel(time1);
sinT1 = 0*(1:numel(time1));
% use the RCG VCO algorithm to determine the values of sine and cosine
for ii = 1:numel(time1) - 1;
cosNew = cosT1(ii)*(1 - 0.5*dphi^2) ...
- (dphi)*sinT1(ii);
sinNew = sinT1(ii)*(1 - 0.5*dphi^2) ...
+ (dphi)*cosT1(ii);
cosT1(ii+1) = [ cosNew];
sinT1(ii+1) = [ sinNew];
end
% extract the phase from the VCO values of sine and cosine
phaseT = unwrap(angle(cosineT + i* sineT));
phaseT1 = unwrap(angle((cosT1 + i*sinT1)));
% determine the phase error for 1 cycle
diff = phaseT1 - phaseT;
% determine the frequency error
slope = (diff(2) - diff(1))/(dt);
diffList = [diffList, slope];
disp(f)
pause(0.001)
end
% plot the results
close all
figure
orient landscape
loglog(flist, abs(diffList/(2.0*pi)))
xlabel('Requested VCO Frequency (Hz)')
ylabel('Frequency error (Hz)')
grid on
print('-dpdf', '/users/abrooks/VCO_error.pdf')
|
Attachment 1: VCO_error.pdf
|
|
4259
|
Tue Feb 8 10:23:02 2011 |
Aidan | Summary | Green Locking | Digital Frequency Discriminator - reference |
Here's the reference for the self-reference frequency detection idea. See Figure 2.
http://www.phys.hawaii.edu/~anita/new/papers/militaryHandbook/mixers.pdf |
4260
|
Tue Feb 8 13:26:11 2011 |
Aidan | Summary | Green Locking | Temperature dependence of phase change of green on reflection |
I did a quick back of the envelope calculation of the expected green phase change on reflection from the aLIGO ITM.
The phase change per nm, K1 = delta phi/delta Lambda, around 532nm is ~1.5 degrees/nm (from the LMA data) [this number is approximately 100x smaller at 1064nm]
I assumed that very small changes in the thickness of the coating appear equivalent to shifting the spectra for reflection/transmission/phase-change-on-reflection up or down by delta lambda, where
delta Lambda/Lambda = delta h/h
where h is the total thickness of the coating and delta h is the change in the thickness of the coating.
Assume that delta h/h = alpha deltaT, where alpha is the coefficient of thermal expansion and delta T is the change in temperature. (approximately 1K)
Then delta phi = K1* Lambda * alpha * delta T = 1.5 degrees/nm * 532nm * 10^-5 K^-1 * 1.0 K = 8 * 10^-3 degrees.
Assume that 360 degree phase change corresponds to one FSR.
Therefore, the frequency shift due to temperature change in the coating = 8*10^-3/360 * FSR = 2.2 *10^-5 * FSR.
Therefore, the expected frequency shift per degree temperature change = 2.2*10^-5 * FSR [Hz/K] |
4271
|
Thu Feb 10 14:52:36 2011 |
Aidan | HowTo | Computers | Adding filenames in MATLAB plots |
The following code is incredibly useful when creating MATLAB plots as it adds the filename of the script to the plot itself. I think it should be used for all MATLAB plots that go on the elog.
For example, I have no idea where the data/script is that was used to generate these plots.
orient landscape
xposn = 0.0;
yposn = -0.13; % you sometimes have to tweak this value depending on the page size and the number of subplots
text(xposn,yposn,[filename('fullpath'), '.m'], ...
'Units', 'normalized', ...
'Interpreter', 'none', ...
'FontSize', 6)
print('-dpdf', [filename('fullpath'), '.pdf'])
|
4351
|
Thu Feb 24 17:42:00 2011 |
Aidan | Update | Green Locking | 15% of end laser sideband power transmitted through cavity |
I did a quick calculation to determine the amount of sideband transmission through the FP cavity.
The modulation frequency of the end PDH is 216kHz. The FSR of the cavity is about 3.9MHz. So the sidebands pick up about 0.17 radians extra phase on one round trip in the cavity compared to the carrier.
The ITM reflectance is r_ITM^2 = 98.5% of power, the ETM reflection is r_ETM^2 = 95%.
So the percentage of sideband power reflected from the cavity is R_SB = r_ITM*r_ETM*(exp(i*0.17) - 1)^2 / (1 - r_ETM*r_ITM exp(i*0.17) )^2 = 0.85 = 85%
So about 15% of the sideband power is transmitted through the cavity. The ratio of the sideband and carrier amplitudes at the ETM is 0.05
So, on the vertex PD, the power of the 80MHz +/-200kHz sidebands should be around sqrt(0.15)*0.05 = 0.02 = 2% of the 80MHz beatnote.
Once we get the green and IR locked to the arm again, we're going to look for the sidebands around the beatnote.
|
4362
|
Sun Feb 27 09:43:59 2011 |
Aidan | Update | Green Locking | sidebands on beatnote |
Can we set up a fiber-PD on the end table to look at the beat between the "end laser IR beam" and the "PSL IR beam fiber-transmitted end beam"?
We should see the same thing on that PD that we see on the green PD (plus any fiber noise and I'm not really sure how much that'll be off the top of my head). If we unlock the lasers from the arm cavity then the free-running noise of the lasers wrt to each other will probably swamp the 50kHz and 150kHz signals. Maybe we could lock the end laser to the free-running PSL by demodulating the beat note signal from the fiber-PD and then we could look for the extra sidebands in the IN-LOOP signal. Then we could progressively lock the PSL to the MC and arm cavity and see if the sidebands appear on the fiber-PD at some point in that process.
It's possible that the 216kHz drive of the PZT on the Innolight is somehow driving up some sub-harmonics in the crystal. I think this is unlikely though: if you look at Mott's measurements of the Innolight PZT response, there are no significant PM resonances at 50 or 150kHz.
Quote: |
Other than the +/-216 kHz sidebands, we can see some funny peaks at +/- 50 kHz and +/-150 kHz.
Quote: #4351 by Aidan |
So, on the vertex PD, the power of the 80MHz +/-200kHz sidebands should be around sqrt(0.15)*0.05 = 0.02 = 2% of the 80MHz beatnote.
Once we get the green and IR locked to the arm again, we're going to look for the sidebands around the beatnote.
|
|
|
4365
|
Tue Mar 1 08:42:18 2011 |
Aidan | Update | elog | Restarted the elog this morning |
The elog was dead this morning. I reanimated it. It is now undead. |
Attachment 1: Zombie.gif
|
|
4368
|
Wed Mar 2 17:19:58 2011 |
Aidan | Configuration | Green Locking | Moved PDH PD on end table |
As previously noted, there are multiple beams coming back from the ETM surface onto the PDH PD on the end table. They are spread out in a vertical pattern. All the spots swing together (as the ETM moves?).
I moved the PDH Green PD on the end table so that it was further away from the Faraday and I added an iris in between the Faraday and the PD. Now only the principle reflection from the ETM is incident on the PD. See attached photos. In order to sneak past the neighbouring optics, I had to steer the beam down a bit, so the PD is now lower than it previously was.
Just FYI: the angle between the returning beams is about 5 or 6 mrad at the PD. Before the beams get to the PD they go through a telescope that de-magnifies the beam by about 5 or 6 times. This implies that the angle between adjacent returning beams from the ETM is around 1 mrad at the ETM.
This does make the position of the spot on the PD more susceptible to the alignment of the ETM - we should use a short focal length lens and image the ETM plane onto the PD.
First image - overview of table
Second image - the three returning beams immediately before the IRIS
Third image - a close up of the IRIS and PDH PD.
|
Attachment 1: P1000223-a.jpg
|
|
Attachment 2: P1000218.jpg
|
|
Attachment 3: P1000224.jpg
|
|
4369
|
Wed Mar 2 18:08:43 2011 |
Aidan | Update | Green Locking | Ghost beams on green |
Kiwamu and I noticed that there is a ghost beam on the green beam going into the ETM. What we see is some interference fringes on the edge of the transmission of the green beam through the dichroic beam splitter (DCBS). If we look at the reflection from the dichroic beam splitter these are much more pronounced.
The spacing of the fringes (about 2 per 10mm) indicates an angle between the fields of around 0.1 mrad.
We were able to cause significant motion of the fringes by pushing on the knobs of the steering mirrors that steer the beam into the DCBS. A rough calculation of the derivative of optical path difference between the ghost and the primary beam as a function of input angle gives about 15 microns per mrad. What filtering the effect the arm cavity will have on the ghost beam is not immediately clear, but the numbers shouldn't be too difficult to determine.
|
Attachment 1: Ghost_Beam_at_ETM_DCBS.pdf
|
|
4428
|
Wed Mar 23 08:50:36 2011 |
Aidan | Update | Green Locking | servo handig off |
Nicely done!
Quote: |
Succeeded in handing off the servo from the green to the red.

|
|
Attachment 1: green-to-red.jpg
|
|
4435
|
Wed Mar 23 19:16:17 2011 |
Aidan | Summary | Green Locking | Y-END green equipment is all available |
With the exception of a 2" mirror mount, I've confirmed that we have everything for the Y-end green production and mode-matching.
We need to calculate a mode-matching solution for the Lightwave laser so that it gives the correct beam size in the doubling crystal.
Additionally, Rana has suggested that we change the pedestals from the normal 1" diameter pedestal+fork combo to the 3/4" diameter posts and wider bases that are used on the PSL table (as shown in the attached image). |
Attachment 1: three-quarter_inch_pedestal.jpg
|
|
4443
|
Fri Mar 25 08:58:38 2011 |
Aidan | Update | Treasure | Cleared stuff off SP table |
I tidied up some of the stuff that was on the SP table. The ISS box that has been sitting on there for months is now underneath the X-arm on top of the spare Marconi which is stored there.
|
Attachment 1: ISS_box.JPG
|
|
4460
|
Wed Mar 30 16:32:29 2011 |
Aidan | Configuration | Computer Scripts / Programs | Added a sitemap alias |
I added an alias to the sitemap MEDM screen in /cvs/cds/caltech/target/cshrc.40m
Now you can enjoy launching sitemap from a terminal.
alias sitemap 'medm -x /cvs/cds/rtcds/caltech/c1/medm/sitemap.adl'
|