TNI mode cleaner/ laser frequency stabilization system

Similar documents
Configuration Study of Pre-Mode Cleaner and Reference Cavity in the 40m PSL System

Multiply Resonant EOM for the LIGO 40-meter Interferometer

레이저의주파수안정화방법및그응용 박상언 ( 한국표준과학연구원, 길이시간센터 )

The Pre Stabilized Laser for the LIGO Caltech 40m Interferometer: Stability Controls and Characterization.

LASER INTERFEROMETER GRAVITATIONAL WAVE OBSERVATORY -LIGO- CALIFORNIA INSTITUTE OF TECHNOLOGY MASSACHUSETTS INSTITUTE OF TECHNOLOGY

Installation and Characterization of the Advanced LIGO 200 Watt PSL

Optical Vernier Technique for Measuring the Lengths of LIGO Fabry-Perot Resonators

Development of Optical lever system of the 40 meter interferometer

The VIRGO injection system

ISC RF Photodetector Design: LSC & WFS

Arm Cavity Finesse for Advanced LIGO

Quantum States of Light and Giants

The VIRGO detection system

AM Stabilized RF Amplifier Driver

Homework Assignment 13

7th Edoardo Amaldi Conference on Gravitational Waves (Amaldi7)

The VIRGO suspensions

Broadband Photodetector

LIGO II Photon Drive Conceptual Design

EXPERIMENTAL SETUP AIMED TO STUDY THE ELECTRICAL IMPEDANCE VARIATIONS OF A PLASMA COLUMN IN A WIDE FREQUENCY RANGE

Optical design of shining light through wall experiments

Notes on the Pound-Drever-Hall technique

Commissioning of Advanced Virgo

Interferometer signal detection system for the VIRGO experiment. VIRGO collaboration

Timing Noise Measurement of High-Repetition-Rate Optical Pulses

Part 2: Second order systems: cantilever response

How to Build a Gravitational Wave Detector. Sean Leavey

A gravitational wave is a differential strain in spacetime. Equivalently, it is a differential tidal force that can be sensed by multiple test masses.

This is a brief report of the measurements I have done in these 2 months.

Calibration of the LIGO displacement actuators via laser frequency modulation

LIGO SURF Progress Report II: Squeezer

A review of Pound-Drever-Hall laser frequency locking

Optical generation of frequency stable mm-wave radiation using diode laser pumped Nd:YAG lasers

The Virgo detector. L. Rolland LAPP-Annecy GraSPA summer school L. Rolland GraSPA2013 Annecy le Vieux

Control Servo Design for Inverted Pendulum

Output Mode Cleaner Design

FFP-TF2 Fiber Fabry-Perot Tunable Filter Technical Reference

Homework Assignment 13

Spatial Uniformity of Silicon Photodiodes at Radio Frequencies

Testbed for prototypes of the LISA point-ahead angle mechanism

Constant Current Control for DC-DC Converters

Mode mismatch and sideband imbalance in LIGO I PRM

Advanced Operational Amplifiers

Virgo status and commissioning results

Readout and control of a power-recycled interferometric gravitational wave antenna

DESIGN OF COMPACT PULSED 4 MIRROR LASER WIRE SYSTEM FOR QUICK MEASUREMENT OF ELECTRON BEAM PROFILE

SIGNAL RECOVERY: Sensors, Signals, Noise and Information Recovery

Wavelength Control and Locking with Sub-MHz Precision

New Technique Accurately Measures Low-Frequency Distortion To <-130 dbc Levels by Xavier Ramus, Applications Engineer, Texas Instruments Incorporated

LASER INTERFEROMETER GRAVITATIONAL WAVE OBSERVATORY - LIGO - CALIFORNIA INSTITUTE OF TECHNOLOGY MASSACHUSETTS INSTITUTE OF TECHNOLOGY

High Current, High Power OPERATIONAL AMPLIFIER

TCS beam shaping: optimum and achievable beam profiles for correcting thermo-refractive lensing (not thermo-elastic surface deformation)

A simple high-sensitivity interferometric position sensor for test mass control on an advanced LIGO interferometer

Loop Compensation of Voltage-Mode Buck Converters

Last Name Girosco Given Name Pio ID Number

HIGH POWER LASERS FOR 3 RD GENERATION GRAVITATIONAL WAVE DETECTORS

Advanced Virgo commissioning challenges. Julia Casanueva on behalf of the Virgo collaboration

LASER INTERFEROMETER GRAVITATIONAL WAVE OBSERVATORY -LIGO- CALIFORNIA INSTITUTE OF TECHNOLOGY MASSACHUSSETTS INSTITUTE OF TECHNOLOGY

SUPPLEMENTARY INFORMATION DOI: /NPHOTON

R. J. Jones Optical Sciences OPTI 511L Fall 2017

Core Technology Group Application Note 6 AN-6

Testing and Stabilizing Feedback Loops in Today s Power Supplies

DRAFT Expected performance of type-bp SAS in bkagra

The AEI 10 m Prototype. June Sina Köhlenbeck for the 10m Prototype Team

Twelve voice signals, each band-limited to 3 khz, are frequency -multiplexed using 1 khz guard bands between channels and between the main carrier

R. J. Jones College of Optical Sciences OPTI 511L Fall 2017

Class D audio-power amplifiers: Interactive simulations assess device and filter performance

Residual Phase Noise Measurement Extracts DUT Noise from External Noise Sources By David Brandon and John Cavey

Glossary of VCO terms

visibility values: 1) V1=0.5 2) V2=0.9 3) V3=0.99 b) In the three cases considered, what are the values of FSR (Free Spectral Range) and

Nonlinear Macromodeling of Amplifiers and Applications to Filter Design.

Increasing the laser power incident on the recycling mirrors in the LIGO interferometers

Design of a Regenerative Receiver for the Short-Wave Bands A Tutorial and Design Guide for Experimental Work. Part I

SC5307A/SC5308A 100 khz to 6 GHz RF Downconverter. Datasheet SignalCore, Inc.

Alignment control of GEO 600

Results from the Stanford 10 m Sagnac interferometer

Background (What Do Line and Load Transients Tell Us about a Power Supply?)

Noise Budget Development for the LIGO 40 Meter Prototype

Using a Negative Impedance Converter to Dampen Motion in Test Masses

INTERFEROMETRIC SENSING AND CONTROL

First and second order systems. Part 1: First order systems: RC low pass filter and Thermopile. Goals: Department of Physics

PLL Synchronizer User s Manual / Version 1.0.6

LIGO Photodiode Development and Optical Platform for LIGO Photodetectors Testing

Texas Components - Data Sheet. The TX53G1 is an extremely rugged, low distortion, wide dynamic range sensor. suspending Fluid.

Andrea Zanchettin Automatic Control 1 AUTOMATIC CONTROL. Andrea M. Zanchettin, PhD Winter Semester, Linear control systems design Part 1

Circuit Design and Implementation of Micro-Displacement Measurement System of Laser Self-Mixing Interference

Stable Recycling Cavities for Advanced LIGO

Differential Amplifier : input. resistance. Differential amplifiers are widely used in engineering instrumentation

University of Pittsburgh

Final Report for IREU 2013

Nonlinear Macromodeling of Amplifiers and Applications to Filter Design.

Power supplies are one of the last holdouts of true. The Purpose of Loop Gain DESIGNER SERIES

Theory and Applications of Frequency Domain Laser Ultrasonics

First step in the industry-based development of an ultra-stable optical cavity for space applications

Universal and compact laser stabilization electronics

Parametric signal amplification

Specify Gain and Phase Margins on All Your Loops

MASSACHUSETTS INSTITUTE OF TECHNOLOGY HAYSTACK OBSERVATORY WESTFORD, MASSACHUSETTS

LFR: flexible, clip-around current probe for use in power measurements

DUAL ULTRA MICROPOWER RAIL-TO-RAIL CMOS OPERATIONAL AMPLIFIER

Transcription:

LASER INTERFEROMETER GRAVITATIONAL WAVE OBSERVATORY -LIGO- CALIFORNIA INSTITUTE OF TECHNOLOGY MASSACHUSETTS INSTITUTE OF TECHNOLOGY Technical Note LIGO-T000077-00- R 8/10/00 TNI mode cleaner/ laser frequency stabilization system Eric Black, Seiji Kawamura, Luca Matone, Shanti Rao This is an internal working note of the LIGO Project. California Institute of Technology Massachusetts Institute of Technology LIGO Project - MS 51-33 LIGO Project - MS 20B-145 Pasadena CA 91125 Cambridge, MA 01239 Phone (626) 395-2129 Phone (617) 253-4824 Fax (626) 304-9834 Fax (617) 253-7014 E-mail: info@ligo.caltech.edu E-mail: info@ligo.mit.edu WWW: http://www.ligo.caltech.edu/ file /home/lmatone/tni/mc/t000077.ps- printed September 5, 2000

Preliminary Version last revised September 5, 2000 Contents 1 Introduction 2 2 Requirements 2 3 Topology 3 3.1 Experimental setup................................... 3 3.2 The open loop transfer function............................. 7 4 Effective frequency fluctuations 7 5 Measurements 8 5.1 Noise measurements.................................. 9 5.1.1 Shot noise and electronic noise......................... 9 5.1.2 MC and laser noise............................... 10 5.2 Spurious path measurement............................... 10 5.3 Finesse measurement.................................. 11 5.4 Transfer function measurement............................. 12 5.4.1 The coil driver transfer function measurement................. 12 5.4.2 First internal mass resonance.......................... 12 5.4.3 NPRO slow actuator.............................. 12 5.4.4 Old NPRO frequency stabilization servo.................... 12 6 Further optimization 12 6.1 Modifications to the slow path............................. 13 6.2 Modifications to the common path........................... 13 A The UGF measurement 14 B The crossover frequency 14 1

Abstract This document describes work done on the Thermal Noise Interferometer s mode cleaner/laser frequency stabilization servo while Seiji Kawamura was visiting, July 17 through August 11, 2000. Keywords Thermal Noise Interferometer 1 Introduction This document describes the current frequency stabilization scheme developed for the Thermal Noise Interferometer (TNI). It consists of 1. an NPRO 126 laser 500mW and 2. a suspended triangular cavity, or mode-cleaner designed to meet the requirement on the residual frequency fluctuation of δ ν 30 mhz/ Hz (1) 2 Requirements The laser frequency noise is approximately ( 2 10 4 /f ) Hz/ Hz (2) from 10 Hz to 10 khz [1] If the desired displacement sensitivity for the test cavity is then the tolerated residual noise is δ x =10 19 m/ Hz (3) δ ν = ν 0 L 0 δ x =3mHz/ Hz (4) where ν 0 =3 10 14 Hz is the laser frequency and L 0 =1cm is the test cavity length. Assuming a conservative Common Mode Rejection Ratio (CMRR) of 10, the requirement can be relaxed to 30 mhz/ Hz. In order to achieve such frequency stabilization, the design of the servo amplifier must have a gain G of approximately G = 7000 @ 100 Hz (5) G = 700 @ 1 khz 2

Laser Pockel cell Mode-Cleaner Photodiode Photodiode Filter Oscillator Mixer Amplifier 3 Topology Filter Figure 1: The topology used for the stabilization in frequency of the laser. The control scheme consists of two paths: 1. at low frequencies, the cavity is locked to the laser, and 2. at high frequency, the laser follows the cavity. We have observed that the laser is more stable than the damped, suspended cavity at frequencies below a few Hertz, and we expect that the cavity will be more stable than the laser at higher frequencies. The topology used is shown in fig.(1). 3.1 Experimental setup The block diagram equivalent of fig.(1) is shown in fig.(2) with the following transfer functions: 1. describing IC =2.5 ω MC µv s ω MC Hz (6) (a) the mode-cleaner transfer function where ω MC /2π is the cavity pole; (b) the photodetector transfer function; (c) the modulation and demodulation process; 3

2. A = 100 V V (7) is a common stage to the two paths. As of this writing, it is composed of an SR560 and a high-voltage (±30 V) amplifier in series; 3. E pzt = ω 1 s ω 2 V s ω 1 ω 2 V (8) describes the passive circuit, shown in fig.(3), build to compensate the high frequency path by acting on the PZT, where ω 1 /2π = 100Hz (9) ω 2 /2π = 60kHz 4. s ω 3 ω 4 ω 5 V E MC = G MC ω 3 s ω 4 s ω 5 V (10) describes the passive circuit, whose current configuration is shown in fig.(4) but will be replaced by the filter shown in fig.(9), that acts on the coil driver of the MC test mass, where 5. P pzt =4MHz/V is the transfer function for the laser actuator and G MC = 2.5 10 3 (11) ω 3 /2π = 1Hz ω 4 /2π = 10Hz ω 5 /2π = 30Hz P mc =13 ω 2 p s 2 ω p Q s ω2 p µm V (12) is the transfer function for the coil driver and suspension together, with a resonance at ω p /2π = 1 Hz and a damped quality factor Q of approximately 10; 6. f 0 =3 10 14 Hz and L 0 =0.5m are the nominal laser frequency and mode cleaner cavity length; 7. ν is the uncontrolled laser frequency fluctuations and 8. x mc is the uncontrolled cavity length fluctuations. Once in closed loop, the path of interest is the one to the laser and the block diagram can be greatly simplified as shown in fig.(2) where M [V/Hz] describes the IC, A and the MC path transfer functions. 4

Laser frequency fluctuations ν Residual fluctuations δν - M Ppzt PZT actuator Epzt Filter PZT y Laser frequency fluctuations ν Residual fluctuations δν - correction signal to PZT - Mode-cleaner, photodetector modulation and demodulation process I C demodulated voltage MC length fluctuations x mc f0/l0 Residual fluctuations δx mc correction signal to MC mass - Pmc Coil driver and suspension Ppzt Filter MC mass Emc Epzt Amplifier (HV) A y PZT actuator Filter PZT Figure 2: The block diagram of the experimental setup. The green path is acts on the laser PZT while the red path acts on the coil driver actuator of the MC end mass. 5

Amplitude x1 1/f 100Hz 60kHz Frequency R1=2k Ω C1=1.4nF C3=700nF Vin R1 C1 C2 Vout pole @ 100Hz zero @ 60kHz Figure 3: Amplitude x2.5e-2 f 1/f x2.5e-3 1Hz 10Hz 30Hz Frequency R1=9.1kΩ R2=80kΩ R3=200Ω C1=2 µ F C3=24 µ F Vin R2 C1 R1 zero @ 1Hz pole @ 10Hz Vout R3 C3 pole @ 30Hz Figure 4: 6

10 4 Total Open Loop OL PZT path OL MC path Amplitude 10 2 10 0 10 1 10 0 10 1 10 2 10 3 10 4 10 5 f(hz) Phase (deg) 0 50 100 Total Open Loop OL PZT path OL MC path 150 200 10 1 10 0 10 1 10 2 10 3 10 4 10 5 Figure 5: The open loop transfer function. Green: slow path; red: fast path; blue: global. 3.2 The open loop transfer function A model of the open loop transfer function is shown in fig.(5), where the slow path is shown in green, the fast path is in red and the global transfer function is in blue. The model shows a DC gain of 5000 with a Unity Gain Frequency (UGF) of about 100 khz and a crossover at about 4 Hz. We have measured the UGF and found it to be 80 khz. We have ample phase margin (about 90 ) at the UGF, and we expect to be able to push the UGF to 100 khz without difficulty. We have measured the phase at UGF and found it to be 140 : the model does not take into account the PZT, amplifiers and photodetector phase delays. f(hz) 4 Effective frequency fluctuations The laser frequency fluctuation is effectively attenuated only as we lock the laser onto the suspended mode cleaner cavity. However, no stabilization is achieved once we lock the laser onto the mode cleaner (the slow path) and therefore we d like to investigate what is the effective frequency stabilization for the servo designed so far. Referring to the top drawing in fig.(2), the residual frequency fluctuation δν can be written as δν = 1 1G eff ν (13) 7

where and G eff = G pzt 1G mc (14) G pzt = IC A E pzt P pzt (15) G mc = IC A E mc P mc f0 L 0 Let s look at two different frequency regions: low frequency, where the mode cleaner follows the laser, and high frequency where the laser follows the mode cleaner. 1. For the low frequency region and the effective gain G eff reduces to G mc > 1 (16) G eff G pzt (17) G mc and the effective frequency stabilization is the ratio between the fast path gain G pzt and the slow path gain G mc ; 2. at high frequency G mc < 1 (18) and the effective frequency stabilization is then dominated by the fast path G eff G pzt (19) 5 Measurements We performed a series of measurements to characterize the system: 1. noise, both electronic and photon-shot, 2. a spurious path around the PZT actuator, 3. the finesse of the mode cleaner, 4. the transfer functions of various stages in the loop, and 5. a possible test mass resonance measurement. We will describe each of these measurements in more detail below. 8

Lightbulb Photodetector Oscillator 12.33MHz Mixer Low-Pass Filter SR560 high pass 100Hz low pass 10kHz G=100 Vout 5.1 Noise measurements Figure 6: Three noise sources have been considered: electronic noise of demodulated signal, shot-noise and mode-cleaner/laser noise. 5.1.1 Shot noise and electronic noise We measured the electronic noise of just the photodiode and demodulator and compared it with the shot noise of the system. We generated shot noise by shining a flashlight on the photodiode, producing the same dc power level as the laser. (Remember that the laser beam is heavily attenuated before it is sent into the photodiode.) The setup for this measurement is shown in fig.(6), where the output goes to a spectrum analyzer. We performed three kinds of measurements: the noise with the flashlight off, the noise with the flashlight on, and the noise with the SR560 disconnected from the photodiode and demodulator. In each case, the noise with the SR560 disconnected was substantially less than the noise with the photodiode and demodulator connected. We did this zero measurement both by grounding the SR560 s input and by disconnecting the photodiode and oscillator from the mixer, leaving the mixer and its low-pass filter connected to the amplifier. The noise was the same for both configurations and was typically more than 20 db less than the noise with the electronics connected. A variety of measurements of the electronic and shot noise were taken, with SR560 gains of 10 2, 10 3, and 10 4. All measurements were done at 1 khz. Shot noise and electronic noise combined was typically not more than 2 or 3 db greater than electronic noise alone, indicating that the electronic noise and the shot noise are comparable. If N elec is the electronic noise alone, and N shotelec is the combined noise, we may find the shot noise by N shot = Nshotelec 2 N elec 2. The electronic noise depends weakly on the gain of the SR560 and was found to be between 50 and 65 nv/ Hz for all gains used. Not surprisingly, the shot noise did not depend on the SR560 s gain and was found to be 43.2 ± 8.1nV/ Hz. We may conclude that our measurement system is no more than about a factor of 2 noisier than the shot noise limit. Furthermore, since IC =2.5µV/ Hz, the limit to frequency stabilization determined by shot 9

noise is N shot IC = 43nV/ Hz 2.5µV/ Hz =1.7mHz/ Hz << 30m Hz/ Hz (20) which is well below the required stability of 30 mhz/ Hz without taking in consideration a CMRR of 10. 5.1.2 MC and laser noise By locking the laser to the reference cavity and the mode cleaner to the laser, a low frequency measurement was possible of the noise level of both the mode cleaner and the laser. For the laser PZT V pzt and coil driver V mc monitors, we found that V pzt = 40 db V RMS @ 100mHz (21) V mc = 85 db V RMS @ 100mHz corresponding to 40 khz RMS for the laser and 400 khz RMS for the mode cleaner: the mode cleaner is 10 times noisier that the laser @ 100 mhz. For this reason, a slow path is used to control the mode cleaner to the laser for frequencies below 100 mhz. 5.2 Spurious path measurement The possible presence of spurius paths can limit the gain used in the servo system. For the TNI successful operation, it is mandatory to measure and resolve any possible spurious paths. One of them is the amplitude modulation generated when acting on the laser PZT and it can be described by the block diagram shown in fig.(7). Here, I describes the cavity transfer function, E the electronics, F - I X E Figure 7: F the PZT actuator and by X we indicate and unknown transfer function that describes the spurious link. The objective of this measurement is to quantify X and understand if it does indeed play a role. The open loop transfer function G for the block diagram shown in figure is E G = IF 1XE (22) 10

Oscillator 12.33MHz Laser Photodiode signal out Pockel cell Mixer signal in Figure 8: If X E > 1 then G I F (23) X and the gain G is limited by the X transfer function. From both the measurement of the open loop transfer function and the model shown in fig.(5), we induce that IF =10 (24) and we know that E = 10 @ 1 khz. Fig.(8) shows the experimental setup to measure the X transfer function. Laser light is again modulated and directed immediately toward the detection photodiode. A signal is sent to the PZT actuator and the output signal is the demodulated signal. A measurement of the ratio of the two signal, performed at 1kHz, has given According to eq.(23), the limit in gain is X 3 10 5 (25) G < IF X = 10 3 10 =3 5 105 (26) which is a limit high enough compared to the required gain. 5.3 Finesse measurement The finesse has been measured in 3 different ways, always in closed loop and always measuring the cavity pole in the open loop transfer function. The first measurement consisted in exciting the PZT actuator and recording the transfer function. This measurement dated back to the beginning of June 2000. The result has given a pole at 30 khz, corresponding to a finesse of 5000. A drawback in the use of the PZT is an apparent phase delay at frequencies above 10 khz. 11

Two other measurements were made in the month of July 2000. The first of the series was done by acting on the Pockels cell, therefore eliminating the PZT phase delay. This resulted in a cavity pole of 60 khz corresponding to a finesse of 2500. On the other hand, by locking on the first higher order mode, TEM 01, and repeating the measurement using again the Pockels cell, a cavity pole was found at 18 khz, corresponding to a finesse of 8000. The cavity may have been contaminated. For the frequency control, the electronics have been designed and created assuming a cavity pole @ 60 khz. 5.4 Transfer function measurement A series of transfer function measurements have been performed in open and closed loop. Here we briefly describe the most important measurements done so far. 5.4.1 The coil driver transfer function measurement We also measured the coil driver efficiency by locking the mode cleaner to the laser with a UGF at 2 khz. By injecting a sinusoid @ 100 Hz into the PZT actuator, and by monitoring the signal applied to the coil driver, we inferred an efficiency of 13 µm/v at dc assuming a pendulum resonance at 1 Hz. 5.4.2 First internal mass resonance While locking the mode cleaner to the laser, with a UGF at 500 Hz, we observed a resonance at 27.9 khz, which we suspect may be an internal mode of the mode cleaner mirrors. The system was able to acquire and hold lock despite this resonance. 5.4.3 NPRO slow actuator We were able to measure the temperature actuator (slow) of the laser in use. The amplitude of the transfer function rolled off as f 2 above 0.2 Hz, but its phase did not hold steady at 180, as one would expect from Bode s gain-phase relations. This actuator may exhibit some nonlinear behavior above 0.2 Hz. 5.4.4 Old NPRO frequency stabilization servo Transfer functions from INTI/PC OUT INTI/FAST OUT and gain settings have been measured and studied, but we are not using that unit in our present configuration. 6 Further optimization One last optimization was performed on the electronics. While the electronics described thus far has been built, tested and in use, the last optimization was only tested with the aid of the low-noise SR560 preamplifiers and temporary electronics. At the writing of this document, a board is being built with 12

Amplitude Frequency 3mHz Vout R1=100kΩ R2=5kΩ R3=250kΩ C1=220 µ F Vin R1 R2 R3 C1 pole @ 3mHz Figure 9: the new features. These new features would improve the effective frequency stability of the laser beam. 6.1 Modifications to the slow path In order to improve the current apparatus, the idea is to limit the slow path bandwidth to 3 mhz. As a consequence, higher gain is necessary but we are limited by the HV power supply of ± 30 V. By replacing the HV power supply with one delivering ± 50 V (since the pzt path can only accept up to ± 50 V), we can obtain the necessary gain to keep the mode cleaner locked to the frequency excursions. The passive filter replacing the slow path electronics is shown in fig.(9) (recall that the filter used previously is shown in fig.(4)). The switch shown in figure enables or disables the filter due to the high capacitance. 6.2 Modifications to the common path We are also planning, but not yet tested, the use of the boosts located at 10 khz. Both boosts would integrate up to 10 khz and then become flat. We can afford to have two boosts due to a fairly large phase margin. 13

n ν - ε A r m B x δν C Figure 10: Appendix In this appendix we describe clever methods to measure directly the UGF point and crossover frequencies of the servo s open loop transfer function. A The UGF measurement It is possible to measure the OL transfer function of a servo system, in closed loop, without measuring the CL transfer function first. This is shown in fig.(10) where ν is the free laser frequency fluctuations and ABC is the OL transfer function. By injecting a signal n anywhere in the chain, and reading the signals right before and after the injection point (in this case signal r and signal m), then it can be shown with some algebraic manipulation that r n = ABC 1ABC m 1 = n 1ABC r m = ABC (27) B The crossover frequency The crossover frequency can be found by measuring at what frequency References G pzt G mc =1 (28) [1] Rich Abbott, James Mason, and Rick Savage NPRO frequency stabilization LIGO-T970051-00- R (2-6-97). 14