Next Article in Journal
The Rise of Wearable Devices during the COVID-19 Pandemic: A Systematic Review
Next Article in Special Issue
Feasibility of Novel Rear-Side Mirage Deflection Method for Thermal Conductivity Measurements
Previous Article in Journal
A SLIC-DBSCAN Based Algorithm for Extracting Effective Sky Region from a Single Star Image
Previous Article in Special Issue
On-Board Monitoring of SO2 Ship Emissions Using Resonant Photoacoustic Gas Detection in the UV Range
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Investigation and Mitigation of Noise Contributions in a Compact Heterodyne Interferometer

1
Department of Aerospace Engineering, Texas A&M University, 701 H.R. Bright Bldg., College Station, TX 77843, USA
2
Wyant College of Optical Sciences, The University of Arizona, 1630 E. University Blvd., Tucson, AZ 85721, USA
*
Author to whom correspondence should be addressed.
Sensors 2021, 21(17), 5788; https://doi.org/10.3390/s21175788
Submission received: 23 July 2021 / Revised: 19 August 2021 / Accepted: 23 August 2021 / Published: 28 August 2021

Abstract

:
We present a noise estimation and subtraction algorithm capable of increasing the sensitivity of heterodyne laser interferometers by one order of magnitude. The heterodyne interferometer is specially designed for dynamic measurements of a test mass in the application of sub-Hz inertial sensing. A noise floor of 3.31 × 10 11 m / Hz at 100 m Hz is achieved after applying our noise subtraction algorithm to a benchtop prototype interferometer that showed a noise level of 2.76 × 10 10 m / Hz at 100 m Hz when tested in vacuum at levels of 3 × 10 5 Torr. Based on the previous results, we investigated noise estimation and subtraction techniques of non-linear optical pathlength noise, laser frequency noise, and temperature fluctuations in heterodyne laser interferometers. For each noise source, we identified its contribution and removed it from the measurement by linear fitting or a spectral analysis algorithm. The noise correction algorithm we present in this article can be generally applied to heterodyne laser interferometers.

1. Introduction

In the past decades, displacement measuring interferometry (DMI) has extended its application to gravitational wave (GW) detection, including free-falling test mass measurements in space-based gravitational wave detection as the Laser Interferometer Space Antenna (LISA) and its technology demonstrator, LISA Pathfinder [1,2,3,4,5], intersatellite displacement measurements as in the mission GRACE Follow-On that utilizes a Laser Ranging Interferometer [6], and inertial sensor development for seismic activity monitoring in ground-based Laser Interferometer Gravitational-Wave Observatory (LIGO) [7,8]. The above applications require high sensitivity or, in other words, a low noise floor in the frequency regime between tens of micro-Hz to a few Hz. The test mass motion is usually expected to range from a few microns to several millimeters, where a DMI with a large dynamic range is crucial to avoid phase readout ambiguities. Among various DMI techniques, heterodyne laser interferometry has the advantages of multi-fringe measurement range, inherent directional sensitivity, and fast detection speed. The development of a common-mode rejection scheme in heterodyne interferometer [1,9] provides a promising solution to displacement metrology applications due to the high rejection ratio to common-mode noise sources. Further improvements on the interferometer performance towards picometer level sensitivity require the knowledge of the residual noise sources and the contribution of each noise source.
Efforts have been made to identify noise sources in heterodyne interferometry [10,11,12,13,14,15,16,17,18,19,20,21,22,23,24]. The inherent periodic error in heterodyne interferometry is mitigated by spatially separating the beams to reduce frequency and polarization mixing [10,11]. The effect of non-linear optical path difference (OPD) noise caused by electric sideband interference is investigated in [12,13], where its origin and mitigation method by active feedback control are introduced. The laser frequency, as the “ruler” of interferometric measurements, has been studied in terms of noise behaviors [14] and various stabilization methods [15,16,17,18]. Photoreceivers are characterized in [19], where the predominant noise sources are identified, and a photodiode with low current noise is designed. The effects of temperature fluctuations in a low frequency regime are analyzed in [20,21]. Other common noise sources in DMI, such as environmental noises and phasemeter noises, are also identified in previous work [22,23,24]. However, an inclusive post-processing algorithm to characterize and correct multiple noise sources in the system has not yet been published.
In this article, we present a noise identification and mitigation algorithm that can be applied to general heterodyne interferometers. We have tested this in a compact laboratory benchtop interferometer. In Section 2, we introduce the interferometer design, which is built to have a high common-mode rejection ratio and be periodic-error free. In Section 3, we analyze residual noise sources based on the differential measurement results. The contributions from distinct noise sources including the non-linear OPD noise, laser frequency noise, and temperature fluctuations, are identified and subtracted from the differential measurement. The noise correction results show a noise floor of 3.31 × 10 11   m / Hz at 100 m Hz , which is enhanced by an order of magnitude from the original differential measurement, reaching the detection system noise limit above 1 Hz .

2. Compact Heterodyne Laser Interferometer

2.1. Design and Benchtop Prototype

Figure 1a shows the layout of the compact interferometer design. After passing through two individual acousto-optical frequency shifters (AOFS), two laser beams are shifted by frequencies δ f 1 and δ f 2 , respectively. The heterodyne frequency f het is the difference between the two laser frequencies, where f het = δ f 1 δ f 2 . The two beams enter a 50/50 lateral beam splitter (LBS), generating four parallel beams. All four beams propagate through a polarizing beam splitter (PBS) and a quarter-wave plate (QWP) towards the mirrors on the measurement end. They are then reflected back into the PBS, exiting along an orthogonal direction. In the second LBS, two beam pairs interfere and are detected by photodetectors PD R and PD M , respectively. The electric fields of the four beams at the interferometer output are given by
E 1 = E 0 exp i 2 π ( f + δ f 1 ) t + ϕ n ( 0 , δ y ) + Δ ϕ M ,
E 2 = E 0 exp i 2 π ( f + δ f 1 ) t + ϕ n ( δ x , δ y ) + Δ ϕ R ,
E 3 = E 0 exp i 2 π ( f + δ f 2 ) t + ϕ n ( 0 , 0 ) ,
E 4 = E 0 exp i 2 π ( f + δ f 2 ) t + ϕ n ( δ x , 0 ) ,
where E 0 is the nominal amplitude of the electric field. The term ϕ n ( x , y ) denotes the structual noises that imprints on the beam optical pathlengths while propagating through the optical components from the entry plane to the plane defined by the mirror M. This noise term ϕ n is dependent on the beam path (optical axis) position, which is represented by a coordinate ( x , y ) in the plane perpendicular to the beam entry directions, as shown in Figure 1a. The phase terms Δ ϕ M and Δ ϕ R represent the additional contributions from the optical pathlength differences between mirrors M and M M , as well as between M and M R . As depicted in Figure 1b, two beams with the heterodyne frequency difference ( E 1 and E 3 , E 2 and E 4 ) are combined into one beam pair.
Figure 1b shows a top view of the interferometer. Optical components can be bonded or cemented together to reduce the overall footprint, increase the system stability, and simplify the alignment process. In this design, two interferometers are constructed: (a) the measurement interferometer (MIFO) that measures the OPD between the target mirror M M and the fixed mirror M, where ϕ M = ϕ n ( 0 , δ y ) ϕ n ( 0 , 0 ) + Δ ϕ M = ϕ n ( 0 , δ y ) + Δ ϕ M ; (b) the reference interferometer (RIFO) that measures the OPD between the reference mirror M R and the fixed mirror M, where ϕ R = ϕ n ( δ x , δ y ) ϕ n ( δ x , 0 ) + Δ ϕ R = ϕ n ( 0 , δ y ) + Δ ϕ R . The mirror M is affixed to the same base plate with the interferometer to provide a common reference for MIFO and RIFO. The mirror M R is mounted as closely as possible to M M to provide a local reference for target displacement. The irradiance signal detected by the photodetectors PD R and PD M are described by
I R = I R 0 + A cos 2 π f het t + ϕ n ( 0 , δ y ) + Δ ϕ R ,
I M = I M 0 + B cos 2 π f het t + ϕ n ( 0 , δ y ) + Δ ϕ M .
The phases ϕ M and ϕ R can be extracted by various methods such as phase-locked loops (PLL) [25] or discrete Fourier transforms [1]. By taking the difference between individual phase readouts, the target displacement d is calculated by
d = Δ ϕ 2 π · 2 λ = ϕ R ϕ M 2 π · 2 λ = Δ ϕ R Δ ϕ M 2 π · 2 λ ,
where Δ ϕ is the differential phase readout, and λ is the source wavelength. This is referred to as the common-mode rejection scheme, where the noises in the common-mode optical paths ϕ n ( 0 , δ y ) are canceled with the presence of a common reference provided by mirror M.
We built a benchtop prototype with commercial optical components and mechanical mounting parts, as shown in Figure 2. Optical components are clamped together mechanically in this assembly to demonstrate the design concept. The laser source is a tunable diode laser (NP Photonics) operating at a nominal wavelength of 1064 n m . Two AOFS (G&H T-M150-0.4C2G-3-F2P) shift the frequencies upwards by 150 M Hz and 145 M Hz , respectively, in two paths, generating a 5 M Hz heterodyne frequency. A commercial phasemeter (Liquid Instruments Moku:Lab) is used to extract the phase from the heterodyne signal using a PLL algorithm with a sampling frequency of 30.5   Hz . All fibers in the system are polarization-maintaining fibers (PMF) to preserve the interference visibility.

2.2. Operation Environments

The benchtop system is tested in the vacuum chamber with a pressure of 3 × 10 5 Torr to reduce the effects from refractive index fluctuation, acoustic noises, and air turbulence. The laser source is coupled into the system by a fiber feedthrough port, while the optical signal is detected after transmitting through an optical window port. Figure 3 shows the operation environment for the preliminary test. A pendulum platform consisting of a breadboard and four stainless steel wires is set inside the chamber for vibration isolation above its resonance. The first resonance mode of the pendulum platform has a frequency of 305 m Hz measured by a ring-down test. Viton strips are inserted between the frame of the pendulum platform and the base of the vacuum chamber to mitigate thermal conduction. In addition, a thermal insulation box is built outside the chamber with Styrofoam to reduce thermal coupling from the ambient environment.

2.3. Preliminary Test

To evaluate the noise floor of the interferometer, we utilized a single static mirror instead of using individual fixed, reference, and test mirrors. Figure 4 shows the measurement results for MIFO and RIFO, as well as their difference, which represents the sensitivity level of the overall interferometer. An eight-hour measurement is taken in the operation environment described in Section 2.2. Figure 4a, the logarithmic average of the linear spectral density (LSD) plot, shows that the individual interferometers MIFO and RIFO achieve a sensitivity level of 1.64 × 10 7   m / Hz at 100 m Hz . The traces of MIFO and RIFO overlap highly due to the common paths shared between the two interferometers. With the common-mode rejection scheme, the sensitivity is enhanced to the level of 2.76 × 10 10   m / Hz at 100 m Hz . Figure 4b shows the time series of the measured displacement for the individual interferometers and the differential measurement, which is measured to be 7.56 × 10 10   m over the duration of 12 min.

3. Noise Source Characterization and Suppression

The preliminary test demonstrated that the benchtop compact interferometer reduces the interferometer noise floor by rejecting common-mode pathlength noises, thus enhancing the sensitivity level by three orders of magnitude. The residual noise sources in the differential measurement are investigated in this section to determine their contribution to the overall noise floor. The interferometer performance is further improved by combining a linear fit method [26,27,28] with a spectral analysis method.

3.1. Non-Linear OPD Noise

The noise investigation of the LISA Pathfinder interferometer demonstrated that the electromagnetic coupling of the radio-frequency (RF) signals between two AOFS drivers, generates sidebands that imprint on the optical signals [12,13]. Figure 5 shows the measured RF driving signals applied on both AOFS. The sidebands interfere with the RF signals at nominal shifting frequencies, and lead to a heterodyne signal with the frequency equal to the interferometer’s main heterodyne frequency f het . Therefore, this “ghost” signal is detected by the photodetector (PD) and mixed with the actual displacement readout, resulting in a non-linear noise that is dependent on the individual phase readouts ϕ M and ϕ R .
Based on the analysis presented in references [12,13], the non-linear OPD noise caused by the sideband interference can be expressed as
δ ϕ OPD = C 1 cos ϕ M + ϕ R 2 + C 2 sin ϕ M + ϕ R 2 · sin ϕ M ϕ R 2 + C 3 cos ϕ M + ϕ R + C 4 sin ϕ M + ϕ R · sin ϕ M ϕ R .
Equation (8) can be reformatted as
δ ϕ OPD = C · N ( ϕ M , ϕ R ) = [ C 1 , C 2 , C 3 , C 4 ] · cos ϕ M + ϕ R 2 · sin ϕ M ϕ R 2 sin ϕ M + ϕ R 2 · sin ϕ M ϕ R 2 cos ϕ M + ϕ R · sin ϕ M ϕ R sin ϕ M + ϕ R · sin ϕ M ϕ R ,
where C is the coupling coefficient vector, and N is the non-linear OPD term vector in the phase measurement. A linear fit is then performed between the differential phase Δ ϕ and the non-linear OPD terms to estimate the coupling factor C ˜ . The contribution of the non-linear OPD noise is removed, following the equation
Δ ϕ OPDcorr = Δ ϕ δ ϕ ˜ OPD = Δ ϕ C ˜ · N ( ϕ M , ϕ R ) .
Figure 6 shows the LSD for the original differential measurement results and the results after subtracting the non-linear OPD noise. Applying this noise correction leads to a reduction in the noise floor from 2.76 × 10 10   m / Hz to 3.86 × 10 11   m / Hz at 100 m Hz . Table 1 lists the fitted coefficients C and estimated errors δ C during the linear fit process.

3.2. Laser Frequency Noise

The traceability of interferometric measurements is dependent on the accurate knowledge of the laser source frequency. Therefore, fluctuations of the laser frequency in an interferometer with unequal armlengths directly impact the interferometer phase following the equation
δ ϕ freq = 2 π Δ L c δ ν ,
where Δ L is the pathlength difference between two arms of the interferometer, and δ ν is the laser frequency fluctuation. For the proposed interferometer, the phase readout is the differential measurement of two individual interferometers MIFO and RIFO. The contribution of laser frequency noise to the measurement can be expressed by
δ ϕ freq = δ ϕ fM δ ϕ fR = 2 π δ ν c ( L M 1 L R 1 ) + ( L M 2 L R 2 ) ,
where L Mi is the length of the i th arm in MIFO, and L Ri is the length of the i th arm in RIFO. For the nominal design, L M 1 is equal to L R 1 due to the common reference surface provided by the fixed mirror M. Therefore, the magnitude of the laser frequency noise is determined by the axial distance between the reference mirror M R and the measurement mirror M M that is mounted on the test mass. In the preliminary test, only one static mirror is used, which means both terms ( L M 1 L R 1 ) and ( L M 2 L R 2 ) are nominally canceled. However, residual arm length differences remain due to the manufacturing tolerance as well as the imperfect alignment. To evaluate the magnitude of the residual laser frequency noise in our benchtop system, a fiber-based delay-line interferometer (DIFO) is built and integrated with the prototype interferometer as shown in Figure 7.
The DIFO is a heterodyne interferometer with an intentionally designed unequal arm length using a 2-meter fiber in one interferometer arm. The laser frequency fluctuation then gets amplified when coupled into the DIFO phase readout ϕ d according to Equation (12).
Figure 8 shows the linear spectra of the interferometers when injecting a 2 Hz modulation to the laser frequency. The spectrum of DIFO at 2 Hz is 4.82 × 10 7   m , while the spectrum of RIFO is 5.30 × 10 10   m . It can also be noted from Figure 8 that, aside from the injected laser frequency modulation peak and its harmonics, the noise floor of DIFO is almost identical to RIFO due to predominantly common environmental noises. This makes it difficult to extract any excess laser frequency noise from the DIFO readout. Therefore, we first subtract the RIFO phase ϕ R from DIFO phase ϕ D , to mitigate the common path noise effects. The residual phase, ϕ D = ϕ D ϕ M , is used to estimate the noise coupling factor. Moreover, the residual phase is band-passed in the frequency regime between 1 m Hz and 1 Hz to mitigate temperature fluctuation and mechanical vibration effects.
We perform a linear fit of the band-passed phase, ϕ D , to the differential phase, to estimate the coupling coefficient K ˜ . The noise correction procedure, which is applied to the original data, can be described as
Δ ϕ freqcorr = Δ ϕ δ ϕ ˜ freq = Δ ϕ K ˜ · ϕ D .
Figure 9 shows the results of identifying the laser frequency noise contribution based on the OPD-noise-corrected differential phase described in Section 3.1. The laser frequency contribution is identified to be 2.28 × 10 12   m / Hz at 100 m Hz in this case, and is not the dominant noise source in this measurement. We injected a 2 Hz modulation to the laser frequency in order to independently determine its coupling factor and the effective optical arm length mismatch between the interferometers. The results are shown in Figure 8, and we obtain an effective arm length mismatch, Δ L , of 2.31 × 10 5   m . From Equation (11), it is inferred that the effects of laser frequency fluctuations are suppressed by a small Δ L , and therefore is not a limiting factor for the current experimental setup. By substituting Δ L into Equation (11), the laser frequency fluctuations are estimated to be 3.47 × 10 8   Hz / Hz at 1 m Hz .

3.3. Temperature Fluctuation Noise

Temperature fluctuations lead to OPD variations due to changes in the refractive index of the air, and thermo-elastic distortions of the optical elements and mounts. Temperature fluctuations usually limit the sensitivity of an interferometer in the frequency regime below 1 mHz [20,21]. Efforts such as inserting Viton strips and building Styrofoam boxes have been made to mitigate their effects to the operational environment during the experiment. To evaluate the residual thermal contributions to the differential measurements, ambient temperature is monitored by a temperature sensor attached to the outer wall of the vacuum chamber.
In contrast to the two noise sources introduced above, thermal fluctuation effects are recorded with a delay in the optical readout. In other words, the group phase between the temperature measurement and the actual effect on the differential interferometer must be considered in the analysis. Therefore, when analyzing temperature fluctuations, a spectral analysis method [20] is adopted instead of a time-domain linear fit.
Before applying the spectral analysis method, we first evaluate the correlation between the displacement and the temperature measurement by calculating the coherence between them. Figure 10 shows the amplitude and phase of the calculated coherence, and we observe a strong correlation (amplitude larger than 0.5) at frequencies below 1 mHz, while a non-zero phase is observed in the same frequency bandwidth. In our case, this shows that the major limitation presented by temperature fluctuations to our interferometer performance occurs in the sub-mHz regime and has a delay effect on the displacement measurement.
In this spectral analysis method, the transfer function between the temperature and the differential measurement is estimated by the equation
H ( ω ) = S T ϕ ( ω ) S TT ( ω ) ,
where S T ϕ ( ω ) is the cross power spectral density of temperature to differential phase, and S TT ( ω ) is the power spectral density of the temperature measurement. Once the transfer function is determined, the contribution of the temperature fluctuation can be removed from the differential measurement in the spectral domain following
Δ ϕ tempcorr ( ω ) = Δ ϕ ( ω ) H ( ω ) · T ( ω ) .
Figure 11 shows the amplitude of the calculated transfer function between the displacement and the temperature measurement. From Figure 10, the correlation in the high frequency bandwidth above 1 mHz is negligible. Hence, directly applying the transfer function in Equation (15) will add noises in the uncorrelated frequency regime. A low-pass spectral filter is applied to the transfer function in order to remove its contribution, which is also plotted in Figure 11. The cut-off frequency of this low-pass filter is set to be 1 mHz based on the coherence results.
Figure 12 shows the logarithmic averaged LSD of the OPD-noise-corrected differential phase and the temperature correction results based on the spectral analysis method described in Equations (14) and (15). Temperature fluctuations clearly dominate in the low frequency regime below 0.5 mHz. After subtracting this noise contribution, the noise floor at 0.1 mHz is reduced from 6.5 × 10 9   m / Hz to 9.8 × 10 11   m / Hz .

3.4. Detection System Noise Limit

The photodetector (PD) technical noise determines the minimum noise floor that an interferometric system can achieve, aside from fundamental limits such as shot noise. In a PD system, typical noise sources include dark current noise from the photodiode, shot noise originating from potential barriers, Johnson noise from thermal fluctuation in the resistor of the sensor circuits, and capacitive noise from the photodiode and circuit [29]. While other noise sources can be reduced by optimizing the circuit design [19], shot noise depends exclusively on the optical power received by the PD and its responsivity. In a heterodyne interferometer with a PLL phase readout, the minimum detectable displacement d min can be derived when the signal-to-noise ratio (SNR) of the system is equal to 1 [30].
Assuming a small displacement that can be modeled as A d sin ( ω d t ) , the electric field of the interference signal in a general case, is expressed by
E D = E 1 + E 2 = A 1 exp [ i ( ω 1 t + ϕ 1 ) ] + A 2 exp [ i ( ω 2 t + ϕ 2 + A d sin ( ω d t ) ) ] ,
where A i are the amplitudes of the electric fields, ω i are the optical frequencies, and ϕ i are the phase terms of the individual beams. The irradiance on PD then can be represented as
P D = | E 1 + E 2 | 2 = | A 1 | 2 + | A 2 | 2 + 2 A 1 A 2 cos [ Δ ω t + Δ ϕ + A d sin ( ω d t ) ] .
By applying trignometric identities and Bessel expansion to Equation 17, the irradiance on PD is rewritten as
P D = P 1 + P 2 + 2 P 1 P 2 { J 0 ( A d ) cos ( Δ ω t + Δ ϕ ) + m = 1 J 2 m ( A d ) [ cos ( 2 m ω d t + Δ ω t + Δ ϕ ) + cos ( 2 m ω d t Δ ω t Δ ϕ ) ] n = 0 J 2 n + 1 ( A d ) cos [ ( 2 n + 1 ) ω d t + Δ ω t + Δ ϕ ] cos [ ( 2 n + 1 ) ω d t Δ ω t Δ ϕ ] } ,
where P i is the optical power in each inteferometer arm, J α ( x ) is the Bessel function of the first kind. In this case, only small displacement is considered where A d 1 . Based on the small argument approximation for Bessel function where
J 0 ( x ) = 1 , J 1 ( x ) = x 2 , J 2 ( x ) = 1 2 x 2 2 ,
Equation (19) can be simplified as
P D = P 1 + P 2 + 2 P 2 P 2 { cos ( Δ ω t + Δ ϕ ) A d 2 [ cos ( ω d t + Δ ω t + Δ ϕ ) cos ( ω d t Δ ω t Δ ϕ ) ] + O ( A d 2 ) } .
Therefore, the time average of the irradiance and the time-averaged incoherent square sum of the irradiance on the photodetector are
P D = P 1 + P 2 , P D 2 = ( P 1 + P 2 ) 2 + P 1 P 2 A d 2 ,
The actual displacement to be measured, d, is converted by the relation d = λ · A d / ( 2 π · 2 ) . The PD system noise in this case can be decomposed to shot noise based on analytical calculation, and other noises based on experimental measurements, following the expression
i PD 2 = i shot 2 + i other 2 = 2 e ρ P D B + V M 2 G 2 ,
where ρ is the detector responsivity, B is the detector bandwidth, V M is the measured output voltage from PD without incident optical signal, and G is the transimpedance gain and is usually specified by the manufacturer. Assuming the interferometer has a visibility of V, the SNR when PD noises are present can be expressed as
SNR = I D 2 / i PD 2 = k 0 d 4 ρ 2 P 1 P 2 2 e ρ P 0 B + V M 2 / G 2 ,
where k 0 is the wavenumber. When SNR is equal to 1, the minimum detectable displacement is calculated as
d min = 2 e ρ P 0 + V M 2 / ( G 2 · B ) k ρ P 0 V [ m / Hz ] · B .
For the benchtop system, the RMS of output voltage square without detecting signal is V M 2 = 1.29 × 10 7 [ V ] 2 , the optical power on the detector is 4.5 × 10 5   W , and the interferometer visibility is 0.45. By substituting the numerical values into Equation (24), the minimum detectable displacement is calculated to be d min =   1.39 × 10 12   m / Hz .

3.5. Discussions

We investigated three common noise sources for heterodyne interferometry, including non-linear OPD noise, laser frequency noise, and temperature noise. The correction method for each noise source is identified, and the overall post-processing correction algorithm is described in the flow chart in Figure 13.
According to this correction algorithm, Figure 14 shows the differential measurement results (blue trace), the overall noise correction results (red trace), and the contribution of each noise source. The residual noise floor after noise correction shown in Figure 14 is 3.31 × 10 11   m / Hz at 100 m Hz . It is worth mentioning that this algorithm applies to uncorrelated noise sources in general heterodyne interferometry. If there are correlated noise sources, for example as in the DIFO readout and the temperature readout, the data needs to be pre-processed, before applying the correction algorithm to determine the frequency range of their correlation.
This post-processing noise correction algorithm provides an effective alternative to improve the sensitivity of heterodyne interferometers without implementing active stabilizations and can achieve high sensitivity over long-term measurements. Moreover, the application of the post-processing algorithm allows us to maintain an overall compact system footprint, and minimal complexity to the control system. The proposed interferometer and noise correction algorithm achieve a performance that is comparable to a few other techniques in GW applications, and has its own relative advantages, such as low hardware complexity and simple interferometer layout.
In addition to the benchtop system we built in the laboratory, this noise correction algorithm can be applied to general heterodyne interferometers [7,9]. Furthermore, this algorithm makes it possible to improve the performance of existing heterodyne optical readout systems without changing their configurations.

4. Conclusions and Outlook

We have presented a noise identification and noise correction algorithm that applies to general heterodyne interferometry. A compact heterodyne interferometer has been developed and tested to determine its performance, when applying this algorithm to its measurement results. We have developed a mitigation strategy for typical noise sources that are present in heterodyne interferometers such as non-linear OPD noise, laser frequency noise, and temperature fluctuations. The contribution from each noise source is estimated and subtracted from the main phase measurement. After correction, the noise floor is reduced by an order of magnitude to a level of approximately 3.31 × 10 11   m / Hz at 100 m Hz .
There are a few modifications to be made in the future to further improve the interferometer performance over the bandwidth above 1 mHz. On the operation environment side, the natural frequency of the pendulum platform can be reduced below 100 mHz to provide noise attenuation to a broader frequency band. To mitigate the non-linear OPD noise, better management of the crosstalk between RF signals can be implemented, such as building electromagnetic shield for each individual AOFS RF amplifier to reduce the amplitude of the sideband interference. On the instrument side, the interferometer design can be optimized to minimize the common-path noise term, ϕ n in Equation (5), in the individual interferometers MIFO and RIFO. In addition, a laser source with a higher optical output power will reduce the shot noise level and, therefore, the overall interferometer noise floor at higher frequencies. Lastly, it will be interesting to investigate possible implementations of this algorithm using machine learning techniques and predicting filters to achieve a real-time operation and real-time sensitivity enhancement to laser interferometers without adding complex hardware for active stabilizations.

Author Contributions

Conceptualization, Y.Z. and F.G.; Funding acquisition, F.G.; Investigation, Y.Z.; Methodology, Y.Z.; Project administration, F.G.; Supervision, F.G.; Validation, Y.Z., A.S.H. and G.V.; Writing—original draft, Y.Z.; Writing—review and editing, A.S.H., G.V. and F.G. All authors have read and agreed to the published version of the manuscript.

Funding

The authors gratefully acknowledge support from the National Geospatial-Intelligence Agency (NGA) through Grant HMA04762010016, the National Science Foundation (NSF) through Grants PHY-2045579 and ECCS-1945832, and the Jet Propulsion Laboratory (JPL) subcontract 1659175.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors acknowledge Cristina Guzman for editing the contents of this manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Heinzel, G. The LTP interferometer and phasemeter. Class. Quant. Grav. 2004, 21, S581–S587. [Google Scholar] [CrossRef] [Green Version]
  2. Wand, V.; Guzmán, F.; Heinzel, G.; Danzmann, K. LISA Phasemeter development. In Laser Interferometer Space Antenna: 6th International LISA Symposium; American Institute of Physics Conference Series; Merkovitz, S.M., Livas, J.C., Eds.; American Institute of Physics: Melville, NY, USA, 2006; Volume 873, pp. 689–696. [Google Scholar]
  3. Aston, S.M.; Speake, C.C. An Interferometric Based Optical Read-Out Scheme For The LISA Proof-Mass. In Laser Interferometer Space Antenna: 6th International LISA Symposium; American Institute of Physics Conference Series; Merkovitz, S.M., Livas, J.C., Eds.; American Institute of Physics: Melville, NY, USA, 2006; Volume 873, pp. 326–333. [Google Scholar]
  4. Schuldt, T.; Kraus, H.J.; Weise, D.; Braxmaier, C.; Peters, A.; Johann, U. A high sensitivity heterodyne interferometer as optical readout for the LISA inertial sensor. In Proceedings of the Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, San Diego, CA, USA, 13 August 2006; Volume 6293, p. 62930Z. [Google Scholar]
  5. Armano, M. LISA Pathfinder: The experiment and the route to LISA. Class. Quant. Grav. 2009, 26, 094001. [Google Scholar] [CrossRef]
  6. Abich, K.; Abramovici, A.; Amparan, B.; Baatzsch, A.; Okihiro, B.B.; Barr, D.C.; Bize, M.P.; Bogan, C.; Braxmaier, C.; Burke, M.J.; et al. In-Orbit Performance of the GRACE Follow-on Laser Ranging Interferometer. Phys. Rev. Lett. 2019, 123, 031101. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Hines, A.; Richardson, L.; Wisniewski, H.; Guzman, F. Optomechanical Inertial Sensors. Appl. Opt. 2020, 59, G167–G174. [Google Scholar] [CrossRef] [PubMed]
  8. Robertson, N.; Drever, R.; Kerr, I.; Hough, J. Passive and active seismic isolation for gravitational radiation detectors and other instruments. J. Phys. E Sci. Instruments 1982, 15, 1101–1105. [Google Scholar] [CrossRef]
  9. Joo, K.N.; Clark, E.; Zhang, Y.; Ellis, J.D.; Guzmán, F. A compact high-precision periodic-error-free heterodyne interferometer. J. Opt. Soc. Am. A 2020, 37, B11–B18. [Google Scholar] [CrossRef]
  10. Wu, C.M.; Deslattes, R. Analytical Modeling of the Periodic Nonlinearity in Heterodyne Interferometry. Appl. Opt. 1998, 37, 6696–6700. [Google Scholar] [CrossRef] [PubMed]
  11. Joo, K.N.; Ellis, J.D.; Spronck, J.W.; van Kan, P.J.; Schmidt, R.H. Simple heterodyne laser interferometer with subnanometer periodic errors. Opt. Lett. 2009, 34, 386–388. [Google Scholar] [CrossRef] [Green Version]
  12. Wand, V.; Bogenstahl, J.; Braxmaier, C.; Danzmann, K.; Garcia, A.; Guzmán, F.; Heinzel, G.; Hough, J.; Jennrich, O.; Killow, C.; et al. Noise sources in the LTP heterodyne interferometer. Class. Quant. Grav. 2006, 23, S159–S167. [Google Scholar] [CrossRef] [Green Version]
  13. Heinzel, G.; Wand, V.; Garcia, A.; Guzman, F.; Steier, F.; Killow, C.; Robertson, D.L.; Ward, H. Investigation of Noise Sources in the LTP Interferometer S2-AEI-TN-3028. 2008. Available online: http://hdl.handle.net/11858/00-001M-0000-0013-4724-5 (accessed on 7 April 2021).
  14. Salvadé, Y.; Dändliker, R. Limitations of interferometry due to the flicker noise of laser diodes. J. Opt. Soc. Am. A 2000, 17, 927–932. [Google Scholar] [CrossRef]
  15. Drever, R.W.P.; Hall, J.L.; Kowalski, F.V.; Hough, J.; Ford, G.M.; Munley, A.J.; Ward, H. Laser phase and frequency stabilization using an optical resonator. Appl. Phys. B 1983, 31, 97–105. [Google Scholar] [CrossRef]
  16. Supplee, J.M.; Whittaker, E.A.; Lenth, W. Theoretical description of frequency modulation and wavelength modulation spectroscopy. Appl. Opt. 1994, 33, 6294–6302. [Google Scholar] [CrossRef]
  17. Numata, K.; Yu, A.W.; Jiao, H.; Merritt, S.A.; Micalizzi, F.; Fahey, M.E.; Camp, J.B.; Krainak, M.A. Laser system development for the LISA (Laser Interferometer Space Antenna) mission. In Solid State Lasers XXVIII: Technology and Devices; Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series; SPIE: Bellingham, WA, USA, 2019; Volume 10896, p. 108961H. [Google Scholar]
  18. Ellis, J.D.; Joo, K.N.; Buice, E.S.; Spronck, J.W. Frequency stabilized three mode HeNe laser using nonlinear optical phenomena. Opt. Express 2010, 18, 1373–1379. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Guzman Cervantes, F.; Livas, J.; Silverberg, R.; Buchanan, E.; Stebbins, R. Characterization of photoreceivers for LISA. Class. Quant. Grav. 2011, 28, 094010. [Google Scholar] [CrossRef]
  20. Nofrarias, M.; Gibert, F.; Karnesis, N.; Garcia, A.F.; Hewitson, M.; Heinzel, G.; Danzmann, K. Subtraction of temperature induced phase noise in the LISA frequency band. Phys. Rev. D 2013, 87, 102003. [Google Scholar] [CrossRef] [Green Version]
  21. Gibert, F.; Nofrarias, M.; Karnesis, N.; Gesa, L.; Martin, V.; Mateos, I.; Lobo, A.; Flatscher, R.; Gerardi, D.; Burkhardt, J. Thermo-elastic induced phase noise in the LISA Pathfinder spacecraft. Class. Quant. Grav. 2015, 32, 045014. [Google Scholar] [CrossRef] [Green Version]
  22. Badami, V.; de Groot, P. Displacement Measuring Interferometry. In Handbook of Optical Dimensional Metrology, 3rd ed.; Taylor & Francis: Boca Raton, FL, USA, 2013; pp. 178–198. [Google Scholar]
  23. Gerberding, O.; Diekmann, C.; Kullmann, J.; Tröbs, M.; Bykov, I.; Barke, S.; Brause, N.C.; Esteban Delgado, J.J.; Schwarze, T.S.; Reiche, J.; et al. Readout for intersatellite laser interferometry: Measuring low frequency phase fluctuations of high-frequency signals with microradian precision. Rev. Sci. Instruments 2015, 86, 074501. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Ellis, J.D. Field Guide to Displacement Measuring Interferometry; Greivenkamp, J.E., Ed.; SPIE Press: Bellingham, WA, USA, 2014; Volume FG30. [Google Scholar]
  25. Gardner, F.M. Phaselock Techniques, 3rd ed.; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 2005. [Google Scholar]
  26. Cervantes, F.G.; Steier, F.; Wanner, G.; Heinzel, G.; Danzmann, K. Subtraction of test mass angular noise in the LISA Technology Package interferometer. Appl. Phys. B 2008, 90, 395–400. [Google Scholar] [CrossRef] [Green Version]
  27. Marquardt, D.W. An Algorithm for Least-Squares Estimation of Nonlinear Parameters. J. Soc. Ind. Appl. Math. 1963, 11, 431–441. [Google Scholar] [CrossRef]
  28. Guzman, F. Gravitational Wave Observation from Space: Optical Measurement Techniques for LISA and LISA Pathfinder. Ph.D. Thesis, Max Planck Institute for Gravitational Physics & Gottfried Wilhelm Leibniz Universität Hannover, Hannover, Germany, 2009. Available online: https://edocs.tib.eu/files/e01dh09/606793232.pdf (accessed on 7 April 2021).
  29. Watchi, J.; Cooper, S.; Ding, B.; Mow-Lowry, C.M.; Collette, C. Contributed Review: A review of compact interferometers. Rev. Sci. Instruments 2018, 89, 121501. [Google Scholar] [CrossRef] [Green Version]
  30. Dehne, M.; Tröbs, M.; Heinzel, G.; Danzmann, K. Verification of polarising optics for the LISA optical bench. Opt. Express 2012, 20, 27273–27287. [Google Scholar] [CrossRef] [PubMed]
Figure 1. (a) The layout of the compact benchtop heterodyne laser interferometer design; (b) top view of the layout. Two laser beams with different frequencies enter the lateral beam splitter (LBS) and split into four beams, constructing two interferometers measuring interferometer (MIFO) and reference interferometer (RIFO). The MIFO is to measure the optical pathlength difference (OPD), Δ ϕ M and the RIFO is to measure Δ ϕ R . The actual displacement of test mass M M is calculated from the differential measurement between MIFO and RIFO.
Figure 1. (a) The layout of the compact benchtop heterodyne laser interferometer design; (b) top view of the layout. Two laser beams with different frequencies enter the lateral beam splitter (LBS) and split into four beams, constructing two interferometers measuring interferometer (MIFO) and reference interferometer (RIFO). The MIFO is to measure the optical pathlength difference (OPD), Δ ϕ M and the RIFO is to measure Δ ϕ R . The actual displacement of test mass M M is calculated from the differential measurement between MIFO and RIFO.
Sensors 21 05788 g001
Figure 2. Layout of the benchtop compact heterodyne laser interferometer based on the design depicted in Figure 1. Only one static mirror M is used to characterize the system’s noise floor. Commercial optical and mechanical components are used to construct the benchtop system.
Figure 2. Layout of the benchtop compact heterodyne laser interferometer based on the design depicted in Figure 1. Only one static mirror M is used to characterize the system’s noise floor. Commercial optical and mechanical components are used to construct the benchtop system.
Sensors 21 05788 g002
Figure 3. Operation environments. (a) Pendulum platform inside the chamber; (b) styrofoam thermal insulator outside the chamber. The compact interferometer is placed on a pendulum stage consisting of a breadboard and four steel wires. Viton strips are applied between the pendulum frame and the vacuum chamber to reduce ambient thermal coupling. The entire pendulum structure is fixed inside the vacuum chamber.
Figure 3. Operation environments. (a) Pendulum platform inside the chamber; (b) styrofoam thermal insulator outside the chamber. The compact interferometer is placed on a pendulum stage consisting of a breadboard and four steel wires. Viton strips are applied between the pendulum frame and the vacuum chamber to reduce ambient thermal coupling. The entire pendulum structure is fixed inside the vacuum chamber.
Sensors 21 05788 g003
Figure 4. Preliminary test results. (a) LSD and its logarithmic average of the MIFO, RIFO, and differential measurements, respectively; (b) time series of a 12-min section of the MIFO, RIFO, and differential measurements. The traces of measurement results from MIFO and RIFO highly overlap in (a). Experimental results show a noise floor of 1.64 × 10 7   m / Hz at 100 m Hz for individual interferometers, and a noise floor of 2.76 × 10 10   m / Hz at 100 m Hz for the differential measurement, which is enhanced by three orders of magnitude from individual interferometers. The drift for differential measurement is 7.56 × 10 10   m over a period of 12 min.
Figure 4. Preliminary test results. (a) LSD and its logarithmic average of the MIFO, RIFO, and differential measurements, respectively; (b) time series of a 12-min section of the MIFO, RIFO, and differential measurements. The traces of measurement results from MIFO and RIFO highly overlap in (a). Experimental results show a noise floor of 1.64 × 10 7   m / Hz at 100 m Hz for individual interferometers, and a noise floor of 2.76 × 10 10   m / Hz at 100 m Hz for the differential measurement, which is enhanced by three orders of magnitude from individual interferometers. The drift for differential measurement is 7.56 × 10 10   m over a period of 12 min.
Sensors 21 05788 g004
Figure 5. Spectra of RF driving signals for both AOFS. The frequency interval between this sideband and the main peak equals to the heterodyne frequency, leading to a ghost signal of unstable phase responsible for the non-linear OPD noise.
Figure 5. Spectra of RF driving signals for both AOFS. The frequency interval between this sideband and the main peak equals to the heterodyne frequency, leading to a ghost signal of unstable phase responsible for the non-linear OPD noise.
Sensors 21 05788 g005
Figure 6. LSD logarithmic average of original differential measurements and the results after noise correction, and the OPD noise contribution. The noise floor is reduced from 2.76 × 10 10   m / Hz to 3.86 × 10 11   m / Hz at 100 m Hz after applying the noise correction algorithm.
Figure 6. LSD logarithmic average of original differential measurements and the results after noise correction, and the OPD noise contribution. The noise floor is reduced from 2.76 × 10 10   m / Hz to 3.86 × 10 11   m / Hz at 100 m Hz after applying the noise correction algorithm.
Sensors 21 05788 g006
Figure 7. System layout with the delay-line interferometer integrated inside the vacuum chamber. The delay-line interferometer is constructed by inserting a 2-m fiber in one arm of the fiber Mach–Zehnder interferometer to amplify the effects of laser frequency noise.
Figure 7. System layout with the delay-line interferometer integrated inside the vacuum chamber. The delay-line interferometer is constructed by inserting a 2-m fiber in one arm of the fiber Mach–Zehnder interferometer to amplify the effects of laser frequency noise.
Sensors 21 05788 g007
Figure 8. Linear spectra of DIFO, RIFO, and the differential measurement to the injected laser frequency noise by intentionally modulating the laser frequency at a rate of 2 Hz . The spectrum response of DIFO to the laser frequency noise is 4.82 × 10 7   m while the RIFO response is 5.30 × 10 10   m at 2 Hz . After performing the differential operation, the spectrum response is reduced to 5.07 × 10 12   m at 2 Hz .
Figure 8. Linear spectra of DIFO, RIFO, and the differential measurement to the injected laser frequency noise by intentionally modulating the laser frequency at a rate of 2 Hz . The spectrum response of DIFO to the laser frequency noise is 4.82 × 10 7   m while the RIFO response is 5.30 × 10 10   m at 2 Hz . After performing the differential operation, the spectrum response is reduced to 5.07 × 10 12   m at 2 Hz .
Sensors 21 05788 g008
Figure 9. Laser frequency noise contribution estimated by a linear fit, compared to the OPD-noise-corrected differential measurement. The laser frequency noise is estimated to be 2.28 × 10 12   m / Hz at 100 m Hz , and is not the limiting noise source in our benchtop experiment.
Figure 9. Laser frequency noise contribution estimated by a linear fit, compared to the OPD-noise-corrected differential measurement. The laser frequency noise is estimated to be 2.28 × 10 12   m / Hz at 100 m Hz , and is not the limiting noise source in our benchtop experiment.
Sensors 21 05788 g009
Figure 10. The amplitude and phase of the coherence function between the displacement and the temperature measurements. In the low frequency regime below 1 mHz, it shows a strong correlation between the displacement and the temperature measurement as the amplitude is large than 0.5, and a time delay effect as the phase is nonzero.
Figure 10. The amplitude and phase of the coherence function between the displacement and the temperature measurements. In the low frequency regime below 1 mHz, it shows a strong correlation between the displacement and the temperature measurement as the amplitude is large than 0.5, and a time delay effect as the phase is nonzero.
Sensors 21 05788 g010
Figure 11. Transfer function between the displacement and temperature measurements before and after applying a 1 mHz low-pass filter. The amplitude of the transfer function is attenuated by a low-pass filter in the frequency regime where the correlation is negligible.
Figure 11. Transfer function between the displacement and temperature measurements before and after applying a 1 mHz low-pass filter. The amplitude of the transfer function is attenuated by a low-pass filter in the frequency regime where the correlation is negligible.
Sensors 21 05788 g011
Figure 12. Contributions of temperature fluctuations estimated by a spectral analysis method, and the measurement results before and after correction of this noise source. The noise floor is reduced from 6.5 × 10 9   m / Hz to 9.8 × 10 11   m / Hz at 0.1 mHz after applying the noise correction algorithm.
Figure 12. Contributions of temperature fluctuations estimated by a spectral analysis method, and the measurement results before and after correction of this noise source. The noise floor is reduced from 6.5 × 10 9   m / Hz to 9.8 × 10 11   m / Hz at 0.1 mHz after applying the noise correction algorithm.
Sensors 21 05788 g012
Figure 13. Flow chart of the overall noise correction algorithm. The non-linear OPD noise, laser frequency noise and the temperature fluctuation noise are characterized, and the contribution from each noise source is removed with this algorithm.
Figure 13. Flow chart of the overall noise correction algorithm. The non-linear OPD noise, laser frequency noise and the temperature fluctuation noise are characterized, and the contribution from each noise source is removed with this algorithm.
Sensors 21 05788 g013
Figure 14. Overall noise correction results based on the differential measurement and contribution of each noise source. The residual noise floor is 3.31 × 10 11   m / Hz at 100 m Hz , and 9.8 × 10 11   m / Hz at 0.1 mHz after correction. The detection system noise level is 1.39 × 10 12   m / Hz and also marked in this plot.
Figure 14. Overall noise correction results based on the differential measurement and contribution of each noise source. The residual noise floor is 3.31 × 10 11   m / Hz at 100 m Hz , and 9.8 × 10 11   m / Hz at 0.1 mHz after correction. The detection system noise level is 1.39 × 10 12   m / Hz and also marked in this plot.
Sensors 21 05788 g014
Table 1. Linear fitting coefficients and estimated fitting errors.
Table 1. Linear fitting coefficients and estimated fitting errors.
i1234
C i 3.62 × 10 3 1.11 × 10 3 1.89 × 10 6 2.51 × 10 6
δ C i 1.96 × 10 6 1.96 × 10 6 1.23 × 10 6 1.24 × 10 6
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhang, Y.; Hines, A.S.; Valdes, G.; Guzman, F. Investigation and Mitigation of Noise Contributions in a Compact Heterodyne Interferometer. Sensors 2021, 21, 5788. https://doi.org/10.3390/s21175788

AMA Style

Zhang Y, Hines AS, Valdes G, Guzman F. Investigation and Mitigation of Noise Contributions in a Compact Heterodyne Interferometer. Sensors. 2021; 21(17):5788. https://doi.org/10.3390/s21175788

Chicago/Turabian Style

Zhang, Yanqi, Adam S. Hines, Guillermo Valdes, and Felipe Guzman. 2021. "Investigation and Mitigation of Noise Contributions in a Compact Heterodyne Interferometer" Sensors 21, no. 17: 5788. https://doi.org/10.3390/s21175788

APA Style

Zhang, Y., Hines, A. S., Valdes, G., & Guzman, F. (2021). Investigation and Mitigation of Noise Contributions in a Compact Heterodyne Interferometer. Sensors, 21(17), 5788. https://doi.org/10.3390/s21175788

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop