Composite-fringe atom interferometry for high dynamic-range sensing Chen Avinadav,1,2, ∗ Dimitry Yankelev,1,2, ∗ Ofer Firstenberg,1 and Nir Davidson1 1Department of Physics of Complex Systems, Weizmann Institute of Science, Rehovot 7610001, Israel 2Rafael Ltd, Haifa 3102102, Israel Atom interferometers offer excellent sensitivity to gravitational and inertial signals but have lim- ited dynamic range. We introduce a scheme that improves this trade-off by a factor of 50 using composite, moiré-like, fringes, obtained from sets of measurements with slightly varying interrogation times. We analyze analytically the performance gain in this approach and the trade-offs it entails between sensitivity, dynamic range, and bandwidth, and we experimentally validate the analysis over a wide range of parameters. Combining composite-fringe measurements with a particle-filter estimation protocol, we demonstrate continuous tracking of a rapidly varying signal over a span two orders of magnitude larger than the dynamic range of a traditional atom interferometer. I. INTRODUCTION changes from shot to shot. However, imperfections such as non-linearity of the classical sensors, transfer func- Atom interferometry (AI) [1] enables highly sensitive tion errors and misalignment may limit the usefulness of and accurate sensing of gravitational [2–4] and inertial this technique in harsh conditions [32], necessitating pro- forces [5–10]. In addition to laboratory-based experi- longed operation at short T at the expense of sensitivity. ments in fundamental physics, such as tests of general In this work, we introduce an approach to AI which relativity [11–15] and precision measurement of physical increases its dynamic range with minor penalty on sen- constants [16–20], AIs for field-applications are being de- sitivity. We perform a set of measurements with slightly veloped worldwide [21]. Mobile atomic gravimeters and varying values of T , corresponding to slightly different in- gravity-gradiometers [22–31] have been demonstrated for terferometer scale factors. Together, as in a moiré effect, geophysical surveys on land [22, 30], at sea [32], and in the these measurements constitute a composite fringe whose air [33], while cold atom accelerometers and gyroscopes frequency, as well as its phase, encodes the measured in- are developed for inertial navigation [34–36]. Such appli- ertial signal, providing a non-ambiguous dynamic range cations motivate the development of advanced AI tech- larger than measurements with a fixed T . The increase niques for robust operation under conditions of large un- in dynamic range scales inversely with the span of scale certainty and large temporal variations in the measured factors and can reach orders of magnitude, limited only signal. by the experimental SNR. In addition to a static demon- As a phase-measuring instrument, a trade-off exists be- stration, we apply the scheme together with a particle- tween an AI sensitivity and its ambiguity-free dynamic filter estimator to successfully track rapidly-varying sig- range. The ratio of dynamic range to sensitivity is in nals, which change by more than 2π between consecutive general fixed by the signal-to-noise ratio (SNR), whereas measurements and span hundreds of radians altogether, the scale factor, which determines their absolute values, while maintaining high sensitivity. may be controlled by changing the interferometer inter- rogation time T . When prior knowledge of the measured II. EXPERIMENTAL SETUP signal is insufficient, a standard approach involves ini- tial measurements with low sensitivity and high dynamic range (short T ) and gradual progress to measurements We apply the composite fringe approach in a Mach- with high sensitivity and low dynamic range (long T ) Zehnder AI which measures gravitational acceleration [30]. However, the time-averaged sensitivity per √Hz of [40]. A freely-falling cold-atom ensemble interacts with arXiv:1912.12304v2 [physics.atom-ph] 7 Jul 2020 this sequence is greatly reduced. Alternatively, simulta- pulses of counter-propagating laser beams which stimu- neous measurements using two interrogation times was late two-photon transitions. A sequence of three pulses demonstrated in a dual-species interferometer [37] with spatially splits, redirects, and recombines the atomic an improved dynamic range of 5 at the cost of added wavepackets [Fig. 1(a)]. The interferometer phase is × 2 experimental complexity. Ambiguities of AIs may also be φ = (keffg α) T + φL, where keff is the effective atom − resolved by hybridization with classical sensors with large optics two-photon wavevector, g is the gravitational ac- dynamic range [38, 39], especially relevant when the mea- celeration, and T is the time between pulses. The relative sured signal varies continuously in time and substantially frequency of the counter-propagating beams is chirped at a rate α = keffg0, where g0 is an approximate value of g, to compensate the changing Doppler shift of the falling ∗ These authors contributed equally to this work. atoms. φL is a tunable laser phase applied during the [email protected] final π/2-pulse. [email protected] Our apparatus is described in detail in Ref. [41]. 2 scaling: φ varies rapidly with g and represents a (a) /2 /2 (b) T comp max high-resolution measurement, similar to standard mea- T T n surements with fixed T , whereas ωcomp varies slowly with Height g and represents a coarse measurement [Fig. 2(d)]. This T sç ;g min beat frequency, which is revealed by varying T , provides Time the extended dynamic range and stands in close analogy to moiré interferometry [43, 44], while the phase of the composite fringe is still accessible and provides the high Figure 1. (a) Schematic diagram of an atomic gravimeter, operating in a π/2-π-π/2 Mach-Zehnder geometry. Counter- sensitivity. propagating laser beams interact with the atoms and gen- The composite-fringe frequency can be esti- erate the interferometer sequence. (b) Composite-fringe AI mated unambiguously up to ωcomp = π, result- comprises a set a of measurements with interrogation times ing in the extended dynamic range ∆gcomp = T that vary quadratically within a small range. 2 2 π (N 1) / keff Tmax Tmin . With respect to a fixed-−T measurement,− the dynamic range increases by a factor Briefly, we trap and cool 87Rb atoms and launch them upwards using moving optical molasses. Verti- cal, retroreflected Raman beams, derived from a single ∆gcomp 1 N 1 laser modulated at 6.834 GHz, drive Doppler-sensitive = 2 − 2 . (4) ∆gstandard 2 1 Tmin/Tmax two-photon transitions [42] between F =1,m =0 and − | F i F =2,m =0 for state initialization and interferome- At the same time, the gravity sensitivity per shot of a | F i try sequence. The population fraction in F = 2, deter- composite fringe is approximately (see Appendix A) mined by state-dependent fluorescence, constitutes the interferometer output signal. σφ Tmax σg,comp = σg,standard. (5) ≈ keffTminTmax Tmin · III. COMPOSITE FRINGE ANALYSIS AND For Tmin approaching Tmax, the dynamic range enhance- PERFORMANCE ment is significant while the penalty in σg,comp is small. Indeed, Fig. 2(c) presents the gravity residuals δg ob- The standard fringe of an AI is S (φ)= A (C/2)cos φ, tained from fits of the measured composite fringes, ex- where A and C are the fringe offset and− contrast. By hibiting similar sensitivity to measurements at fixed T = measuring this fringe at N points (scanning either α Tmax with no observed systematic bias. 2 or φL) at a fixed T , the gravitational phase keffgT is The potential gain in dynamic range is limited by σφ. determined with uncertainty σφ/√N, with σφ the to- Primarily, when the uncertainty in estimating the com- tal phase uncertainty per shot (see Appendix A). There- posite fringe phase and frequency is comparable to the fore, gravity is determined with uncertainty per shot line separation in Fig. 2(d), a line jump may occur, re- 2 σg,standard = σφ/ keffT , over an ambiguity-free dy- sulting in a large error in estimating g. The criterion for namic range ∆g =2π/ k T 2 . The ratio of dy- avoiding large errors is approximately (see Appendix A) standard eff namic range to sensitivity-per-shot is (∆g/σg)standard = √N T 2 T 2 2π/σφ, depending only on the measurement phase uncer- max − min 1. (6) tainty. σφ TminTmax ≫ Here instead, we form a composite fringe from a set of 2 Notably, the potential gain in dynamic range increases if N measurements Sn, varying the scale factor keffT lin- early by choosing variable interrogation times [Fig. 1(b)], σφ decreases, as Tmin can approach Tmax without increas- ing the error probability. Conversely, the error probabil- T 2 T 2 ity may be reduced by increasing N, which also serves T 2 = T 2 + n max − min for n =0,...,N 1, (1) n min N 1 − to increase the dynamic range, at the cost of temporal − bandwidth. We note that in the limit of T = T , with T < T , keeping α and φ fixed. As min max min max L Eq. (6) cannot be satisfied for finite σ , however such Figs. 2(a),(b) show, the series S forms a composite φ { n} measurements can still be analyzed as in traditional AI fringe, Sn = A (C/2)cos(nωcomp + φcomp), with − within the original dynamic range ∆gstandard, recovering 2 the standard measurement scheme at fixed T . φcomp = (keffg α) Tmin + φL, (2) − 2 2 To understand the trade-offs in the composite fringe Tmax Tmin ωcomp = (keffg α) − . (3) approach, we present in Fig. 3(a) the projected dynamic − N 1 − range enhancement compared to fixed-T measurements Unlike standard fringes with fixed T , here both phase and the corresponding error probability. The latter was and frequency depend on g, albeit with very different evaluated for total phase uncertainty σφ = 400mrad, 3 Figure 2.