Rapid extraction of the phase shift of the cold-atom interferometer via phase demodulation*
Cheng Binga), Wang Zhao-Yinga), Xu Ao-Penga), Wang Qi-Yua), Lin Qianga),b)
Institute of Optics, Department of Physics, Zhejiang University, Hangzhou 310027, China
Center for Optics and Optoelectronics Research, College of Science, Zhejiang University of Technology, Hangzhou 310023, China

Corresponding author. E-mail: zhaoyingwang@zju.edu.cn

Corresponding author. E-mail: qlin@zjut.edu.cn

*Project supported by the National Natural Science Foundation of China (Grant Nos. 11174249 and 61475139), the Ministry of Science and Technology of China (Grant No. 2011AA060504), the National Basic Research Program of China (Grant No. 2013CB329501), and the Fundamental Research Funds for the Central Universities, China (Grant No. 2015FZA3002).

Abstract

Generally, the phase of the cold-atom interferometer is extracted from the atomic interference fringe, which can be obtained by scanning the chirp rate of the Raman lasers at a given interrogation time T. If mapping the phase shift for each T with a series of measurements, the extraction time is limited by the protocol of each T measurement, and therefore increases dramatically when doing fine mapping with a small step of T. Here we present a new method for rapid extraction of the phase shift via phase demodulation. By using this method, the systematic shifts can be mapped though the whole interference area. This method enables quick diagnostics of the potential cause of the phase shift in specific time. We demonstrate experimentally that this method is effective for the evaluation of the systematic errors of the cold atomic gravimeter. The systematic phase error induced by the quadratic Zeeman effect in the free-falling region is extracted by this method. The measured results correspond well with the theoretic prediction and also agree with the results obtained by the fringe fitting method for each T.

Keyword: 37.25.+k; 32.60.+i; atom interferometer; phase extraction; phase demodulation
1. Introduction

The cold-atom interferometer has become an effective tool for high precision measurements. It is useful in both the fundamental physics[16] and the practical application. Quantum inertial sensors based on the cold-atom interferometer have a rapid development in the last two decades. High performance atomic gravimeters, [79] gyroscopes, [10] and gravity gradiometers[11, 12] have been realized. For instance, a sensitivity of 4.2  μ Gal/Hz1/2 has been achieved for the atomic gravimeter in 2013, [13] which has exceeded the best available traditional gravimeter (FG-5). The mobile instruments based on these sensors have been developed worldwide in recent years. Up to now, several transportable atomic gravimeters have been reported.[8, 1417] They are aiming at the applications in geodesy, inertial navigation, Earth gravitational field mapping, engineering geological prospecting, etc.

The evaluation of the accuracy of the atomic gravimeter is complicated due to a numerous class of systematic shifts, which have been discussed in detail in Ref. [18]. The current accuracy is mainly limited by the shifts induced by the Coriolis force and the wave-front aberration of the Raman beams.[19] The systematic error could be divided into two parts: independent items and dependent items on the direction of the effective wavevector of the Raman beams. The method of phase extraction by inverting the effective wavevector of the Raman beams is widely used in analyzing the specific phase shift. It is well known that the phase of the interference fringe is proportional to the chirp rate and scales as the square of the interrogation time T. Typically, another phase introduced during the interferometer is extracted from the interference fringe. In the most simple case, this is achieved by scanning the chirp rate of the Raman beams at a given interrogation time T and doing a sinuous fitting of the fringe. Also, similar to the Ramsey clock, a feedback can be applied to an added phase using the differential signal at the fringe slope as an error signal to achieve null total phase.[20, 21] If more than one interference happens at the same time, the phase difference of the interferometer can be extracted by the Bayesian estimation, [22] the combined graphs of two interference output with elliptic curves[23] or Lissajous curves[24] to reject some common noises.

Nowadays, the interrogation time gets longer and longer to obtain a better precision. The systematic errors induced by the spatial variations have to be considered carefully. In this paper, we present a new method to rapidly extract the total phase of the atom interferometer. The phase modulation is carried out when scanning the T2 at a given chirp rate (around the dark fringe). The total phase as a function of T can be extracted continuously via the demodulation. Therefore the systematic error in the whole interference area is rapidly obtained. As a benefit from this new method, the systematic error induced by the quadratic Zeeman is investigated in our experiment.

2. Theory
2.1. The method of phase demodulation

A Mach– Zehnder type atomic interferometer using the stimulated Raman transition[7] utilizes the combination of π /2– π π /2 Raman pulses of durations τ – 2τ τ . The three pulses are separated by free evolving time T. Generally, the phase of the cold-atom interferometer is extracted from the atomic interference fringe, which is obtained by scanning the chirp rate α of the Raman lasers at a given interrogation time T. For the atoms free falling in the gravitational field, this chirp of the Raman frequency is also used to compensate the Doppler effect, if the direction of the Raman beams are in the vertical configuration. The interference fringe can be described as

where C is a contrast, and A is the ratio of offset.

Now, we introduce a new method of phase demodulation. We fix the chirp rate α , introduce a small positive resident Δ α 0 = keffgα , then change the interrogation time T to step linearly one combined parameter Δ α 0T2, so that we can achieve a linear step of added phase. This means that the measured probability Pn is a function of the parameter , and the parameter Tn changes non-uniformly, here n is the measurement number. The measured probability Pn has the phase modulated cosine form (assuming T2 as the time and Δ α 0 as the frequency), and the modulating phase ϕ (T) is what we want to know. To demodulate the phase ϕ (T), we use a technique similar to phase demodulating for radio frequency by multiplying the fringe (after removing an estimated offset A0) with the reference wave,

where refquphase and refinphase are the quadrature and the in-phase parts of the reference signal, respectively. The reference wave serves as the local oscillations. Then, the products of the multiplication are

After filtering by low-pass and band-stop filters (at the frequency of Δ α 0 and its second harmonic), the two demodulated signal components Pqu and Pin have the form of the first term in Eq.  (3). Finally, we can calculate the demodulating phase ϕ (T) as a function of filtered Pqu and Pin. In this way, it is possible to track the phase larger than 2π . The result of arctan(Pin/Pqu) will give the phase ϕ (T) in the − π /2 to π /2 range. According to the sign of Pqu, we can add π to reach the range of − π /2 to 3π /2. Finally, by examining the gaps of 2π , we can add accumulated counts of 2π to each segment of the phase.

The measurement error in this demodulation method arises from two aspects: the filtering residual of the oscillations and the low-pass effect of filtering. By choosing the proper filter parameters, one can achieve more than 30  dB rejection at the oscillating frequency, which corresponds to 1  mrad in the phase error. The low-pass filter effect will lead to an error if the phase ϕ (T) signal has high frequency components, such as a step function, and the error also scales with the magnitudes of the high frequency components. If we measure close to time T where the step of the phase happens, with higher modulation frequency (larger Δ α 0), we can increase the filter bandwidth and thus increase the time resolution.

Comparing with the method doing a localized cosine fit around any T with a chosen width and weight function, the demodulating method needs less calculation and is easy to treat a large phase shift.

Another utility of this method is to quickly obtain the absolute g value. When there is a second order phase drift ϕ (T) = Δ α T2, the absolute g can be obtained as (α + Δ α 0 + Δ α ′ )/keff. However, when α in Eq.  (1) is larger than keffg, the result of demodulating is indistinguishable with that of the smaller case, and the g should be (α − Δ α 0 − Δ α ′ )/keff. Whether it is larger or smaller, the case can be verified by doing another linear scan of T2 with the chirp rate equal to α + Δ α 0 + Δ α ′ or α − Δ α 0 − Δ α ′ . The right case should give the transition probability always at the minimal value of fringes (if as usually ϕ (T) is less than π /2).

2.2. Extraction of the phase shift of interferometer

During the process of the interference, the quantum axis is defined by the magnetic field produced by the coils. So there is a quadratic Zeeman phase shift in the total phase of the interference, which should be analyzed carefully. In this paper, we introduce three methods for the measurement of the quadratic Zeeman phase shift. The phase of the interferometer due to the quadratic Zeeman energy shift is written as

where Δ ω Zeeman(t) = 2π KB(t)2 is the quadratic Zeeman energy shift, K is the coefficient of the quadratic Zeeman effect, Tc is the circle time of each drop of atoms, and gs(t) is the sensibility function of frequency shift at time t, which has the form

The quadratic Zeeman phase can be calculated once the magnetic field B(t) along the atom falling trajectory is measured. In this paper, we will also extract the quadratic Zeeman phase from the interference fringe.

Usually, reversing the direction of the two Raman lasers can form the interferometer with opposite effective wavevector of the Raman beams. By changing the sign of the chirp rate, the probabilities of atom transition from the initial state to another state are given by

where Δ α 0 = keffgα , ϕ dep is the phase term which changes sign when toggling the direction of the wavevector, and ϕ indep is independent of the direction of the effective wavevector, which includes the quadratic Zeeman phase shift ϕ Zeeman.

We use two methods to extract the quadratic Zeeman phase. The first one is called cosine fitting. For a given T and changing Δ α 0, the phase shift can be obtained by fitting the fringes with a cosine function. The second one is called the phase demodulation as shown in Section 2. The phase demodulations of the probability of two effective Raman wavevector configurations are done independently. Lastly, averaging the phase shifts of the two opposite wavevectors, we can obtain the total independent terms for both fitting and demodulation methods.

To reach the same precision of phase shift estimation σ Δ ϕ , the fringe fitting method needs shots of measurement, where σ ϕ , noise is the equivalent phase noise in each shot. Considering the number of T of the final result Npoint, the total shot of measurement is . On the other hand, the phase demodulation method takes the assumption that the phase shift has no quick change along T, thus one shot of Tn is contributing to adjacent Nfilter points in the final result due to signal filtering. The benefit in the measurement time is 1/Nfilter times for the same amount of result points and goal precision. In the cases where there are spikes or other fast changes in the phase shift, the filtering will spoil the signal, and introduce errors to adjacent points. We can overcome this by increasing the modulation frequency and making an adjustment of Nfilter in the fast change regions to recover the detailed phase shifts.

3. Experiment
3.1. Experimental setup

We use the atomic interferometer in our lab[25] to test the new phase demodulation method. The apparatus includes a cold rubidium 87 atom source and two phase locked lasers. About 108 cold atoms are collected in a three-dimensional magnetic-optical trap (MOT), cooled by sub-Doppler cooling to a temperature of 5  μ K, and then released to free falling. The atoms are prepared in the first order Zeeman effect insensitive states, in this case, the mF= 0 states, and the vertical velocity is selected by one π pulse of velocity sensitive Raman transition.

The Raman lasers are from the D2 line, which couple the two hyperfine structures of the 52S1/2 states by stimulated Raman transition. The Raman laser beams are emitted from one polarization-maintaining fiber and collimated by a lens system. The diameter of the Raman beams is 24  mm, the effective Rabi frequency is 41.7  kHz, and the duration of the π /2 pulse is around 6  μ s.

After three pulses of Raman lasers separated by time T, the interferometer signal, which is the population in the final state, is detected by collecting the fluorescence signal of the atoms scattering the detecting lasers.

3.2. Magnetic field measurement and phase calculating

Firstly, the magnetic field in the interaction area can be measured with the spectral method. The magnetic field is generated by a pair of square coils around the system, which is used to maintain the axis of the quantum states in a residual magnetic background. The overall field is measured by the microwave spectrum[25] of the cold atoms at different falling times. The result is shown in Fig.  1. For the differential measurement, the fields at the coil currents of 0.5  A and 1.0  A are measured independently.

Fig.  1. The magnetic field strength as a function of the falling time: (a)  bias field current 0.5  A, (b) bias field current 1.0  A.

The magnetic field strength is interpolated before integrating to remove the effects of the limit measurement points. Then according to Eq.  (4), the corresponding quadratic Zeeman phase shift can be calculated by integration of the quadratic Zeeman frequency shift along the two interference trajectories. The numerical integration result shows that the overall peak phase error is about 700  mrad for the 1  A current case in Fig.  2.

Fig.  2. The calculated quadratic Zeeman shift versus the interrogation time T.

3.3. Phase measurement by the interference fringe

The interference fringes of pulse separation time T from 10  ms to 70  ms with 5  ms step are realized. The interference signal phase stepping is carried out by changing the chirping rate α of the Raman laser’ s frequency in a range where the change of α T2 is about 4π . The cosine fitting results are shown as squares and triangles in Fig.  3. Usually, the fixed T method takes longer time and gives less result points.

Fig.  3. The measured phase shift terms independent of the Raman wavevector versus the interrogation time T. The solid and dash lines are the scanning T2 phase demodulations with 0.5  A and 1.0  A bias fields, respectively; and the squares and triangles are the fixed T fringe fitting measurements with 0.5  A and 1.0  A bias fields, respectively.

To use the phase demodulation method, the atom transition probabilities are acquired for T2 from 16  ms2 to 4900  ms2 with a step of 4  ms2 (T = 4  ms to 70  ms with non-uniform steps). Using the reference wave of frequency Δ α 0= 2π × 10  kHz to demodulate the probability signal as in Eq.  (3), we need to roughly estimate Δ α 0, which needs to be sufficient for signal demodulating and filtering. In the demodulation, we use averages of 25 points and 13 points as filters to remove the residual harmonics at 10 kHz/s and 20 kHz/s. This method takes about one 25th of the time compared to the T fringe fitting method. The results are shown as lines in Fig.  3.

By comparing the results in Fig.  3 with the theoretical result (Fig.  2), it is obvious that there are other phase errors in the independent term. By only changing the current of the magnetic coils, these other effects may not change. This can be verified by comparing the differential phase shift. In Fig.  4, the differential quadratic Zeeman phase shift between 0.5 A bias field and 1 A bias field is matching the theoretically integrated result. The measurement errors are dominant by the vibration of ground and the response of the passive anti-vibration platform. For T = 70  ms, the fringe fitting method gives a standard error of 17  mrad and the phase demodulation method gives a standard error of 42  mrad. The results of the three methods are consistent within about 100  mrad, which proves that this new phase demodulation method is effective for rapidly extracting the phase of the atom interferometer.

Fig.  4. Comparison of the results of differential phase shift between 0.5  A and 1.0  A bias fields. Line: scanning T2 phase demodulation; dash-dot line: calculated quadratic Zeeman shift; square: fixed T with fringe fitting measurement.

Other phase shifts independent of the wavevector of the Raman beams can be obtained by removing the calculated phase shift from the experiment phase shift, which is shown in Fig.  5. These phase shifts include those caused by the light shift and the frequency reference of the Raman lasers.

Fig.  5. Independent phase shifts other than the quadratic Zeeman shift. Line: scanning T2 phase demodulation; square: fixed T fringe fitting measurement.

4. Conclusion

We present a method of quickly mapping the phase shift of the atom interferometer as a function of the square of interrogation time (T2). Using non-uniformed steps of T but linear steps of one combined parameter Δ α 0T2, we can achieve a linear step of an added phase. Then the demodulation on the atom transition probability is used to extract the modulation phase.

The scan of T2 and phase demodulation method has an advantage of more detailed mapping of the phase shift in smaller steps of T. But as the signal is obtained using demodulation, although with a notch filter, there are still some residual oscillations at the modulation frequency and its harmonics. The residual oscillations may limit the precision of this technique.

To demonstrate the techniques, the quadratic Zeeman shift is analyzed at different bias field currents. The results are compared with the fixed T fringe fitting measurements and the calculated quadratic Zeeman phase shift using direct magnetic field mapping. The results of the three methods are consistent.

Reference
1 Cladé P, Mirand es E, Cadoret M, Khélifa S G, Schwob C, Nez F, Julien L and Biraben F 2006 Phys. Rev. Lett. 96 033001 DOI:10.1103/PhysRevLett.96.033001 [Cited within:1]
2 Fixler J B, Foster G T, McGuirk J M and Kasevich M A 2007 Science 315 74 DOI:10.1126/science.1135459 [Cited within:1]
3 Dimopoulos S, Graham P W, Hogan J M and Kasevich M A 2007 Phys. Rev. Lett. 98 111102 DOI:10.1103/PhysRevLett.98.111102 [Cited within:1]
4 Lamporesi G, Bertoldi A, Cacciapuoti L, Prevedelli M and Tino G M 2008 Phys. Rev. Lett. 100 050801 DOI:10.1103/PhysRevLett.100.050801 [Cited within:1]
5 Dimopoulos S, Graham P W, Hogan J M, Kasevich M A and Rajendran S 2009 Phys. Rev. Lett. 678 37 DOI:10.1016/j.physletb.2009.06.011 [Cited within:1]
6 Hohensee M A, Chu S, Peters A and Müller H 2011 Phys. Rev. Lett. 106 151102 DOI:10.1103/PhysRevLett.106.151102 [Cited within:1]
7 Kasevich M A and Chu S 1992 Appl. Phys. B 54 321 DOI:10.1007/BF00325375 [Cited within:2]
8 Le Gouët J, Mehlstäubler T E, Kim J, Merlet S, Clairon A, Land ragin A and Pereira Dos Santos F 2008 Appl. Phys. B 92 133 DOI:10.1007/s00340-008-3088-1 [Cited within:1]
9 Müller H, Chiow S W, Herrmann S, Chu S and Chung K Y 2008 Phys. Rev. Lett. 100 031101 DOI:10.1103/PhysRevLett.100.031101 [Cited within:1]
10 Gustavson T L, Bouyer P and Kasevich M A 1997 Phys. Rev. Lett. 78 2046 DOI:10.1103/PhysRevLett.78.2046 [Cited within:1]
11 Snadden M J, McGuirk J M, Bouyer P, Haritos K G and Kasevich M A 1998 Phys. Rev. Lett. 81 971 DOI:10.1103/PhysRevLett.81.971 [Cited within:1]
12 McGuirk J, Foster G T, Fixler J B, Snadden M and Kasevich M A 2002 Phys. Rev. A 65 033608 DOI:10.1103/PhysRevA.65.033608 [Cited within:1]
13 Hu Z K, Sun B L, Duan X C, Zhou M K, Chen L L, Zhan S, Zhang Q Z and Luo J 2013 Phys. Rev. A 88 043610 DOI:10.1103/PhysRevA.88.043610 [Cited within:1]
14 Bodart Q, Merlet S, Malossi N, Pereira Dos Santos F, Bouyer P and Land ragin A 2010 Appl. Phys. Lett. 96 134101 DOI:10.1063/1.3373917 [Cited within:1]
15 Merlet S, Bodart Q, Malossi N, Land ragin A, Pereira Dos Santos F, Gitlein O and Timmen L 2010 Metrologia 47 L9 DOI:10.1088/0026-1394/47/4/L01 [Cited within:1]
16 Bidel Y, Carraz O, Charriere R, Cadoret M, Zahzam N and Bresson A 2013 Appl. Phys. Lett. 102 144107 DOI:10.1063/1.4801756 [Cited within:1]
17 Hauth M, Freier C, Schkolnik V, Senger A, Schmidt M and Peters A 2013 Appl. Phys. B 113 49 DOI:10.1007/s00340-013-5413-6 [Cited within:1]
18 Peters A, Chung K Y and Chu S 2001 Metrologia 38 25 DOI:10.1088/0026-1394/38/1/4 [Cited within:1]
19 Louchet-Chauvet A, Farah T, Bodart Q, Clairon A, Land ragin A, Merlet S and Pereira Dos Santos F 2011 New J. Phys. 13 065025 DOI:10.1088/1367-2630/13/6/065025 [Cited within:1]
20 Merlet S, Le Gouët J, Bodart Q, Clairon A, Land ragin A, Pereira Dos Santos F and Rouchon P 2009 Metrologia 46 87 DOI:10.1088/0026-1394/46/1/011 [Cited within:1]
21 Duan X C, Zhou M K, Mao D K, Yao H B, Deng X B, Luo J and Hu Z K 2014 Phys. Rev. A 90 023617 DOI:10.1103/PhysRevA.90.023617 [Cited within:1]
22 Stockton J K, Wu X and Kasevich M A 2007 Phys. Rev. A 76 033613 DOI:10.1103/PhysRevA.76.033613 [Cited within:1]
23 Foster G T, Fixler J B, McGuirk J M and Kasevich M A 2002 Opt. Lett. 27 951 DOI:10.1364/OL.27.000951 [Cited within:1]
24 Chen X, Zhong J Q, Song H W, Zhu L, Wang J and Zhan M S 2014 Phys. Rev. A 90 023609 DOI:10.1103/PhysRevA.90.023609 [Cited within:1]
25 Wu B, Wang Z Y, Cheng B, Wang Q Y, Xu A P and Lin Q 2013 J. Phys. B: At. Mol. Opt. Phys. 47 015001 DOI:10.1088/0953-4075/47/1/015001 [Cited within:2]