Advanced phasemeter for deep phase modulation interferometry

We present the development of an advanced phasemeter for the deep phase modulation interferometry technique. This technique aims for precise length measurements with a high dynamic range using little optical hardware. The advanced phasemeter uses fast ADCs and an FPGA to implement a design of multiple single-bin Fourier transforms running at high sampling rates. Non-linear noise sources in the design were analyzed and suppressed. A null measurement with an optical beatnote signal using λ = 1064nm was conducted. It showed a sensitivity of 0.8μrad/ √ Hz below 10Hz and 13.3μrad/ √ Hz above, with a large dynamic range. The shown performance could enable the measuring of optical path lengths with sensitivities down to 0.14pm/ √ Hz and 2.3pm/ √ Hz, respectively, over several fringes in an interferometric setup. © 2014 Optical Society of America OCIS codes: (120.0120) Instrumentation, measurement, and metrology; (120.3180) Interferometry; (120.4640) Optical instruments; (120.5050) Phase measurement; (120.5060) Phase modulation. References and links 1. O. Sasaki and K. Takahashi, “Sinusoidal phase modulating interferometer using optical fibers for displacement measurement,” Appl. Opt. 27(19), 4139–4142 (1988). 2. K. Danzmann and A. Rüdiger, “LISA technology concept, status, prospects,” Classical Quant. Grav. 20(10), 1-9 (2003). 3. S. Anza, M. Armano, E. Balaguer, M. Benedetti, C. Boatella, P. Bosetti, D. Bortoluzzi, N. Brandt, C. Braxmaier, M. Caldwell, L. Carbone, A. Cavalleri, A. Ciccolella, I. Cristofolini, M. Cruise, M. Da Lio, K. Danzmann, D. Desiderio, R. Dolesi, N. Dunbar, W. Fichter, C. Garcia, E. Garcia-Berro, A. Garcia, R. Gerndt, A. Gianolio, D. Giardini, R. Gruenagel, A. Hammesfahr, G. Heinzel, J. Hough, D. Hoyland, M. Hueller, O. Jennrich, U. Johann, S. Kemble, C. Killow, D. Kolbe, M. Landgraf, A. Lobo, V. Lorizzo, D. Mance, K. Middleton, F. Nappo, M. Nofrarias, G. Racca, J. Ramos, D. Robertson, M. Sallusti, M. Sandford, J. Sanjuan, P. Sarra, A. Selig, D. Shaul, D. Smart, M. Smit, L. Stagnaro, T. Sumner, C. Tirabassi, S. Tobin, S. Vitale, V. Wand, H. Ward, W. J. Weber, and P. Zweifel, “The LTP experiment on the LISA Pathfinder mission,” Classical Quant. Grav. 22(10), 125-138 (2005). 4. R. Pierce, J. Leitch, M. Stephens, P. Bender, and R. Nerem, “Intersatellite range monitoring using optical interferometry,” Appl. Opt. 47(27), 5007–5019 (2008). 5. J. H. Galeti, P. L. Berton, C. Kitano, R. T. Higuti, R. C. Carbonari, and E. C. N. Silva, “Wide dynamic range homodyne interferometry method and its application for piezoactuator displacement measurements,” Appl. Opt. 52(28), 6919–6930 (2013). 6. O. Sasaki and H. Okazaki, “Sinusoidal phase modulating interferometry for surface profile measurement,” Appl. Opt. 25(18), 3137–3140 (1986). #211944 $15.00 USD Received 12 May 2014; revised 16 Jun 2014; accepted 22 Jun 2014; published 21 Jul 2014 (C) 2014 OSA 28 July 2014 | Vol. 22, No. 15 | DOI:10.1364/OE.22.018214 | OPTICS EXPRESS 18214 7. O. Sasaki, H. Okazaki, and M. Sakai, “Sinusoidal phase modulating interferometer using the integrating-bucket method,” Appl. Opt. 26(6), 1089–1093 (1987). 8. G. Heinzel, C. Braxmaier, M. E. Caldwell, K. Danzmann, F. Draaisma, A. Garcia, J. Hough, O. Jennrich, U. Johann, C. J. Killow, K. Middleton, M. te Plate, D. Robertson, A. Rüdiger, R. Schilling, F. Steier, V. Wand, and H. Ward, “Successful testing of the LISA technology package (LTP) interferometer engineering model,” Classical Quant. Grav. 22(10), 149–154 (2005). 9. M. Dehne, M. Tröbs, G. Heinzel, and K. Danzmann, “Verification of polarising optics for the LISA optical bench,” Opt. Express 20(25), 27273–27287 (2012). 10. G. Heinzel, F. G. Cervantes, A. Garcia, J. Kullmann, W. Feng, and K. Danzmann, “Deep phase modulation interferometry,” Opt. Express 18(18), 19076–19086 (2010). 11. S. P. Lipshitz, R. A. Wannamaker, and J. Vanderkooy. “Quantization and dither: A theoretical survey,” in Audio Engineering Society Convention, (Audio Engineering Society, 1991), pp. 3141. 12. T. Xinguang, L. Xin, C. Hong, and D. Miyi. “DDFS spurious signals due to amplitude quantization in absence of phase-accumulator truncation,” J. Syst. Eng. Electron. 20(3), 485–492 (2009). 13. D. Shaddock, B. Ware, P. G. Halverson, R. E. Spero, and B. Klipstein, “Overview of the LISA Phasemeter,” AIP Conf. Proc. 873(1), 654–660 (2006).


Introduction
Optical interferometry systems for high-precision metrology exist in a wide variety of techniques and applications.In this paper, we present advances on the phase measurement unit (phasemeter) of the deep phase modulation interferometry (DPMI) technique.DPMI can be classified as a special case of sinusoidal phase modulation (SPM) homodyne interferometry [1].It aims for precise measurements of optical path length changes with high sensitivity over a large dynamic range (many fringes) while using a simple optical setup.In particular, these optical path length measurements can be applied for the continuous tracking of the displacement of free-floating test masses that act as interferometer end mirrors.This task is required, for example, in space-based gravitational wave detectors such as the Laser Interferometer Space Antenna (LISA) [2,3].In addition, a possible application is the tracking of test masses for other experiments and missions, in the context of satellite geodesy [4] for example.As DPMI is a special case of SPM interferometry, it could naturally be used to perform the tasks carried out by other implementations of SPM interferometry, like characterization of piezo actuators [5] or surface profiling [6,7].The conceptual difference to these methods is the use of a phase modulation with a high ("deep") modulation index to produce a large number of harmonics of this modulation frequency in the beatnote signal.These harmonic amplitudes can then be processed by a fit algorithm to extract the measurement phase.
A comparable technique used for test mass tracking is the heterodyne method developed for LISA pathfinder [8].It showed a sensitivity of better than 1.42 pm/ √ Hz in the frequency band from 3 mHz to 1 Hz in a displacement measurement using a laser wavelength λ = 1064 nm [9].An initial implementation of the DPMI technique was developed in [10].While operating with only one optical phase modulator, the initial implementation of the DPMI reached a sensitivity in the order of 20 pm/ √ Hz for optical path length measurements with λ = 1064 nm.In this work, steps towards further improvement of this sensitivity are presented.The initial experiment in the quoted work was conducted on a stable optical bench in vacuum using a frequency-stabilized laser, and it furthermore utilized corrections for distortions due to the transfer function of the photodetector and data acquisition.Therefore, this work concerns the next key component in the signal chain that shows potential to improve the overall performance: the phasemeter.Its new implementation is based on fast analog-digital converters (ADCs) connected to a field programmable gate array (FPGA), while its initial implementation used a standard PC with an ADC card.The new design's details, pitfalls and performance will be discussed in detail throughout this paper.

Theory
As mentioned before, DPMI can be classified as SPM homodyne interferometry.Therefore, the signal at one output port of the interferometer can be described as: where A is a common scale factor, c is the interferometer contrast, m is the modulation depth, f mod = ω mod /2π is the modulation frequency, ψ is the modulation phase and ϕ is the interferometer phase, which is the main measurand.This expression can be decomposed by Fourier analysis: with the following complex Fourier coefficients: with where A and c are merged to one single variable k = 2Ac.These complex Fourier coefficients are called harmonic amplitudes in the following.While other SPM interferometry techniques only use some of the first harmonic amplitudes, the DPMI technique exploits the fact that Eq. ( 4) forms an over-determined system of equations for four unknowns (k, m, ϕ, ψ).This can be solved by a numerical fit algorithm.The number of useful equations effectively depends on the modulation depth m.For a deeper modulation over several radian, the signal's power is distributed into more and higher harmonic bins.Therefore, an experimental setup is required that is able to apply a high phase modulation and extract the harmonic amplitudes for subsequent processing by the numerical fit.Earlier theoretical analysis [10] investigated optimal values for m and the number of utilized bins, N. Accordingly, the number of bins was chosen to be N = 10.The modulation depths used were adapted slightly depending on which value led to the best figure of merit of the numerical fit in the specific simulation or measurement.They are stated separately for each case.

Initial and advanced phasemeter
The initial setup for the DPMI phasemeter [10] implemented the extraction of the harmonic amplitudes with a commercial ADC card (sampling rate f s = 20 kHz, 16 bit resolution) connected to a standard PC, on which a fast Fourier transform (FFT) was performed.Only the desired bins at the harmonic frequencies were selected and processed by a Levenberg-Marquardt fit to solve the associated system of equations.The phase modulation of the signal in one interferometer arm was carried out by a function generator driving a ring piezo.The ring piezo stretches a fiber which is coiled around it.This increases the optical path length and therefore introduces a phase shift.The modulation had a frequency of f mod = 280 Hz and a depth of m ≈ 9.7.The frequency value was harmonically related to a FFT resolution of 20 Hz, as the frequency itself and higher harmonics should lie in the bins of the FFT.For this work, the hardware responsible for the generation of the modulation signal as well as the harmonic amplitude extraction was replaced.Both are now embedded in an FPGA receiving the beatnote signal from an FPGA Mezzanine card (FMC) ADC ( f s = 40 MHz, 12 bit resolution).Additionally, a digital-analog converter (DAC) card is attached to the FPGA to drive the ring piezo.The FPGA contains a direct digital synthesis (DDS) unit for the generation of the modulation signal, a digital signal processing (DSP) unit for the harmonic amplitude extraction as well as infrastructural components like DAC/ADC interfaces and an Ethernet interface for the data link to a PC.With this hardware configuration, a total amount of eight separate phasemeter channels can be implemented and used for simultaneous readout of eight interferometers, limited by the eight input channels of the ADC card.

Advanced phasemeter design details
The task of the DSP unit is to extract the harmonic amplitudes.This is done by performing multiple single-bin Fourier transforms (SBFTs).Instead of calculating a complete FFT, only the bins with the frequency values of interest -which are the harmonics of f mod -are computed.To achieve this, two pipelines for every chosen bin are implemented, one for the real part and one for the imaginary part of the (complex) harmonic amplitude.Each pipeline consists of a mixer and a cascaded integrator comb (CIC) filter, which averages over several consecutive samples and is therefore sampling down the mixer signal with a rate R. The mixer is fed by the beatnote signal converted by the ADC as well as by a sinusoidal wave with the frequency of the bin to be extracted or demodulated.The latter one is generated in the DDS unit, which consists of N numerically controlled oscillators (NCOs).Figure 1 shows a sketch of the design for one bin.The mixer effectively converts down just the desired harmonic amplitude to DC.
from ADC: from NCO: Fig. 1.Sketch of two pipelines performing the SBFT to calculate the harmonic amplitude's real and imaginary part of one bin at the n-th harmonic frequency.
All others harmonics are removed by the CIC filter, as illustrated in Fig. 2. It shows the photodetector signal and the transfer function of the filter.The frequency axis is scaled to multiples of the sampling frequency after downsampling, f s /R, which are the points of maximum filter suppression.R has to be set to a value that makes f mod an integer multiple of f s /R to filter out the undesired AC components.As mentioned earlier, the DDS furthermore transmits the fundamental frequency f mod to the DAC for driving the phase modulator.Figure 3 shows an overview over the whole digital design on the FPGA.The harmonic amplitudes sent to the PC are processed by a C program implementing the Levenberg-Marquardt fit.It can be classified as a least-square fit algorithm.Subsequently to the fit, a software CIC filter is applied to sample down ϕ and thus keep the amount of data sufficiently low.Improvements in comparison to the initial implementation are expected first of all due to the higher ADC sampling rate.This effectively reduces the level of white noise added by the signal quantization.The behavior was put into a formula [11], which can be altered to:  where ASD noise denotes the amplitude spectral density of the white noise introduced by the ADC, LSB is the voltage quantization of the least significant bit of the ADC given in Volts and f s is the ADC sampling rate.It should be noted that this describes only the additive noise to the beatnote signal amplitude.The calculation of the associated phase noise -and thus the expected improvement -is rather complicated due to the behavior of the DSP system and the fit.It will be subject to future studies.

Non-linear noise sources
Before the design was implemented, investigations of different non-linear noise sources had been conducted.These were done by simulating MATLAB/Simulink models of the system.In addition to the DDS and DSP components presented, a model for the optical part, i.e. the interferometer, was developed by evaluating Eq. ( 1).This provides the possibility of generating a reference input signal for the measurement phase ϕ in and comparing it to the output of the phasemeter ϕ out .For this analysis, sinusoidal signals for ϕ in were chosen to probe the design.
One of the critical elements turned out to be the fact that one has to take into account that the SBFT of the photodetector signal is computed in finite time windows, the lengths of which are controlled by the downsampling rate of the CIC filter.In case of constant values for the four parameters k, m, ψ and ϕ, this causes no problems.However, one has to consider that while k, m, ψ can actually be treated as constants in an ideal setup, ϕ changes in the experiment and is therefore time-dependent.A sinusoidal change of ϕ within the time window for example leads to sidebands besides the harmonic frequencies.While these sidebands carry the actual signal, the problem becomes visible when recalling the principle of the SBFT and in particular the CIC filter.The AC components remaining after mixing are supposed to be canceled out completely.That is not the case for the sidebands, which are then aliased and corrupt the desired signal.Figure 4 illustrates this behavior.The issue was mitigated by choosing a higher CIC filter order and therefore a better suppression of the unwanted sidebands.However, this method is limited in terms of hardware consumption and therefore a trade-off had to be made.As a second order filter proved sufficient in the simulations shown in the next section, it was chosen for the final implementation with which the measurements for this work were conducted.Nevertheless, a design with a higher CIC order or another filter might yield improvements for future measurements.
A second critical item to be taken care of were the NCOs.An NCO consists of a phase accumulator and a lookup table containing the period of a sine or cosine.The phase accumulator sums up a value called phase increment, which defines the NCO frequency, and drives the lookup table.The phase accumulator's bitwidth is usually higher than the one of the lookup table input.Therefore, its signal has to be truncated, which generates spurious signals.Those were dealt with using a dither applied to the signal before truncation.For high frequencies and thus high phase increments, spurious signal due to amplitude truncation in the lookup table have to be considered too [12].These can naturally be reduced by choosing higher amplitude resolutions.Furthermore, choosing optimized values for the frequency and the phase increment, respectively, will decrease the repetition rate of values in the phase accumulator signal.This in turn spreads the spurious signal over several spurs, decreasing their level and hence increasing the spurious free dynamic range of the NCOs.Additionally, higher modulation frequencies lead to a higher free dynamic range, as they require more bits in the phase increment.They are therefore desirable up to a point where the Nyquist theorem or the bandwidth of the used devices limit a further increment.
were caused by a non-flat transfer function of the data acquisition system and the photodetector.This leads to unequal damping and phase shifts of the different complex harmonic amplitudes, which in turn creates errors in a n not covered by the used model.It is assumed that this issue is also the reason for the shoulder visible here.Nevertheless, the measured phase noise would correspond to an optical path length noise of 2.3 pm/ √ Hz.Furthermore, a significantly higher dynamic range compared to the purely electrical setup is observable here.It reaches approximately up to seven orders of magnitude.This shows the possible applicability of DPMI for tracking over many fringes.
It should be noted, that the used null measurement cannot reveal nonlinearities that are common in both phasemeter channels.This would require a more sophisticated so-called three signal test, an analysis scheme that has already been developed for other phasemeters [13].Three different measurement phases ϕ a , ϕ b and ϕ c are combined to three new signals containing ϕ 1 = ϕ a − ϕ b , ϕ 2 = ϕ b − ϕ c and ϕ 3 = ϕ a − ϕ c .The latter ones are put into three phasemeter channels.Their output ϕ 1 + ϕ 2 + (−ϕ 3 ) should in an ideal case be zero.As different signals are fed into the channels, nonlinearities that are common in the null measurement can be revealed this way.Despite the lack of a three signal test, high confidence is put in the current version of the phasemeter due to fact that the absolute digital simulation showed no nonlinearities.

Conclusions
In this paper, advances in the development of a phasemeter for DPMI were discussed.This technique comes with a simple optical setup, using only one phase modulator, here implemented by a ring piezo.The new phasemeter is using a multi-channel ADC (40 MHz) attached to an FPGA implementing multiple SBFTs per channel for fast sampling and simultaneous processing of up to eight beatnote signals.The phasemeter performance in terms of phase noise was first estimated with a simulation of the digital design, which allowed an absolute comparison between an artificial input signal and the phasemeter output.A phase noise of 0.6 µrad/ √ Hz with a dynamic range of seven orders of magnitude was shown.Furthermore, a null measurement in an electrical testbed yielded a similar performance with a white noise floor of 0.8 µrad/ √ Hz, limited to lower frequencies by temperature effects.This measurement was repeated with an optical setup, i.e. performing a null measurement with an optical beatnote captured on a photodetector.It showed a phase noise of 13.3 µrad/ √ Hz for frequencies below 1 Hz.For higher frequencies, a convergence towards the white noise floor measured with the electrical testbed was observed.This measured white noise floor -which was also achieved in the simulations with a dynamic range of up to seven orders of magnitude -would correspond to 0.14 pm/ √ Hz in an optical path length measurement (λ = 1064 nm).At frequencies below 1 Hz, the measured noise would correspond to 2.3 pm/ √ Hz with a dynamic range of also up to seven orders of magnitude.Comparisons to earlier experiments suggest that the white noise floor can also be reached for lower frequencies by transfer function corrections of the data acquisition line and the photodetector.This could enable high measurement precision with high bandwidth and a large dynamic range over several orders of magnitude and therefore over many beatnote fringes.In particular if assuming the phasemeter of earlier experiments as the major noise source, improvement of DPMI to even higher sensitivities is possible.
To actually show these sensitivities, future work will include further testing of the phasemeter, for example with more sophisticated methods like a three signal test, and putting it into an interferometric setup for an optical displacement measurement.

Fig. 3 .
Fig. 3. Sketch of the digital system and its individual units.The green lines symbolize commands, black lines internal digital signals and blue lines digital signals connected to the analog

Fig. 4 .
Fig.4.Illustration of the origin of distortions occurring through the dynamics of ϕ: the sidebands on all harmonics can't be canceled out completely by the CIC filter and are aliased down to DC, disturbing the measurement of the desired harmonic amplitude in the specific pipeline.
Fig.2. of the DSP pipeline principle: after the desired harmonic amplitude is mixed down to DC, all other signal components are filtered out by the CIC filter, if its downsampling rate R is set correctly.