Next Article in Journal
Neutrons and Gamma-Ray Dose Calculations in Subcritical Reactor Facility Using MCNP
Next Article in Special Issue
Analysis of Polarizability Measurements Made with Atom Interferometry
Previous Article in Journal / Special Issue
A Wigner Function Approach to Coherence in a Talbot-Lau Interferometer
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Prospects for Precise Measurements with Echo Atom Interferometry

1
Department of Physics & Astronomy, York University, 4700 Keele, Toronto ON M3J 1P3, Canada
2
Institut d’Optique d’Aquitaine, rue François Mitterand, Talence 33400, France
*
Author to whom correspondence should be addressed.
Atoms 2016, 4(3), 19; https://doi.org/10.3390/atoms4030019
Submission received: 30 March 2016 / Revised: 20 June 2016 / Accepted: 21 June 2016 / Published: 27 June 2016
(This article belongs to the Special Issue Atom Interferometry)

Abstract

:
Echo atom interferometers have emerged as interesting alternatives to Raman interferometers for the realization of precise measurements of the gravitational acceleration g and the determination of the atomic fine structure through measurements of the atomic recoil frequency ω q . Here we review the development of different configurations of echo interferometers that are best suited to achieve these goals. We describe experiments that utilize near-resonant excitation of laser-cooled rubidium atoms by a sequence of standing wave pulses to measure ω q with a statistical uncertainty of 37 parts per billion (ppb) on a time scale of ∼50 ms and g with a statistical precision of 75 ppb. Related coherent transient techniques that have achieved the most statistically precise measurements of atomic g-factor ratios are also outlined. We discuss the reduction of prominent systematic effects in these experiments using off-resonant excitation by low-cost, high-power lasers.
PACS:
37.25.+k, 06.20.-f, 37.10.Jk

Graphical Abstract

1. Introduction

Over the past few decades, there has been sustained interest in using the exquisite sensitivity of atom interferometric techniques to gain a precise knowledge of the fundamental forces that govern our universe through an improved understanding of light-matter interactions. Among pioneering advances in this field include the diffraction of atomic beams using standing wave light fields [1,2] and micro-scale material beam splitters [3,4], and sensitive measurements of the index of refraction of an atomic gas [5]. The development of atom interferometers (AIs) to measure fundamental constants [6] and inertial effects [7,8] using laser-cooled atoms showed the potential of AIs for realizing precise studies of fundamental physics, and for industrial applications such as oil and mineral prospecting or inertial navigation. During the last 25 years, there has been steady progress toward developing AIs and coherent transient techniques [9,10,11] for measurements of fundamental constants such as α (or the ratio of Planck’s constant to the mass of the test atom h / M ) [12,13,14,15,16], the gravitational constant G [17,18,19,20], studies of inertial effects such as gravitational acceleration [8,21,22], gravity gradients [23,24,25] and rotations [7,26,27,28,29], sensing magnetic gradients [30,31,32,33] and for more sensitive tests of the equivalence principle [34,35,36,37,38,39,40,41,42] and general relativity [43,44].
Most of the progress in both short-term and long-term sensitivity has been achieved using Raman interferometers [8,22,45,46,47,48]. This AI relies on optical velocity-sensitive two-photon Raman transitions between two long-lived hyperfine ground states in alkali atoms, such as the F = 1 and F = 2 states in 87 Rb. Raman AIs were the first implementation of state-labeled interferometer [45], where the coherent exchange of photon momentum between the atoms and the optical fields is associated with a change in the internal atomic state. Hence, all the information regarding the interference between atoms is stored in the relative population between states F = 1 and F = 2 .
Despite the well-developed nature of Raman AIs, a number of alternate interferometer configurations have been developed [49,50,51,52,53], particularly for measurements of gravity. In this article, we report on recent developments and techniques for precision measurements using a unique, single-state echo interferometer [54,55]. This AI requires only a single excitation frequency, and does not require velocity selection. Recently, this configuration has achieved several milestones, including an extension of the timescale to the transit time limit (∼250 ms for experimental conditions) [32], as well as significant improvements in statistical precision relating to measurements of the atomic recoil frequency (related to h / M ) [56,57] and the acceleration due to gravity [58]. In what follows, we review recent progress on both atomic recoil (Section 2) and gravity (Section 3) measurements. In these sections, we also discuss various methods of reducing or eliminating the dominant systematic effects which are currently limiting the measurements. In Section 4, we review related coherent transient techniques [59,60] that have demonstrated precise measurements of atomic g-factor ratios. Finally, in Section 5 we describe the development of a new class of auto-locked semiconductor diode lasers operating at 780 nm and 633 nm [61]. These low-cost, high-power laser sources exhibit impressive long-term lock stability that will be implemented in future generations of the experiments described in Section 2 and Section 3. Finally, we conclude with some perspectives in Section 6.
In the following subsections, we compare the operating principles of a Raman interferometer with those of a single-state grating-echo AI, including a brief theoretical description of the two types.

1.1. Description of Raman-Type Interferometers

As previously mentioned, Raman-transition-based interferometers rely on coherently transferring atoms between internal states of the atom. To make these transitions, two counter-propagating optical fields are used, one at frequency ω 1 with wavevector k 1 , and the other at ω 2 with k 2 . These Raman fields satisfy the resonance condition ω 1 ω 2 = ω HF + k eff · v + ω k eff for making two-photon transitions between ground states F = 1 , p and F = 2 , p + k eff . Here, ω HF is the hyperfine splitting between F = 1 and F = 2 , p = M v is the initial momentum of the atom, k eff = k 1 + k 2 is the effective wavevector of the counter-propagating Raman fields, and ω k eff = k eff 2 / 2 M is the atomic recoil frequency associated with the Raman transition.
The most basic and widely used type of Raman interferometer is the Mach-Zehnder configuration [8] which consists of a π / 2 π π / 2 sequence of pulses separated by a time T, as shown in Figure 1a. The first π / 2 -pulse acts as a beam-splitter that creates a 50/50 superposition of the states 1 , p and 2 , p + k eff . The atoms then travel along two spatially separated pathways. The π-pulse at t = T acts as a mirror, exchanging the population between the two states and redirecting the wavepackets associated with each trajectory back toward one another. The final π / 2 -pulse at t = 2 T recombines the wavepackets by “closing” the interferometer pathways and producing the interference. The two output ports of the interferometer correspond to the relative populations in each state, for instance N 1 / ( N 1 + N 2 ) , where N F represents the number of atoms in state F . These populations are usually measured via resonant fluorescence, where many photons can be scattered per atom. Since the Raman interferometer excites only two pathways, the fringe pattern follows a simple sinusoidal function
N 1 , 2 N 1 + N 2 P 1 , 2 = 1 2 1 ± C cos Δ Φ tot ,
where P 1 , 2 is the probability of finding the atom in either state at the output of the interferometer, C is the contrast of the interference fringes, and Δ Φ tot is the total interferometer phase difference. The key idea of a Raman interferometer is that the population between internal states oscillates as a function of the phase difference between interfering pathways. In general, this phase consists of three main contributions [62]
Δ Φ tot = Δ ϕ prop + Δ ϕ sep + Δ ϕ las ,
where Δ ϕ prop = L ( r , v ) d t / is the propagation phase corresponding to the difference of classical action (integral of the lagrangian L ) along the upper and lower pathways, Δ ϕ sep = p · Δ r / is a phase associated with a spatial separation Δ r between the wavepackets during the final π / 2 -pulse, and Δ ϕ las is due to the Raman laser phase imprinted on the atoms during each pulse, which is given by
Δ ϕ las = k eff · r c ( 0 ) 2 r c ( T ) + r c ( 2 T ) + φ 1 2 φ 2 + φ 3 .
Here, r c ( t ) represents the center-of-mass trajectory of the atom and φ j is the phase difference between the Raman fields at the j th pulse.
Due to the limited bandwidth of two-photon Raman transitions given by the Rabi frequency Ω eff / 2 π 10 100 kHz, and the desire to use only atoms in magnetically “insensitive” m F = 0 states, a large percentage of atoms are lost during the sample preparation process. To give a typical example, a sample of N 10 8 laser-cooled 87 Rb atoms at a temperature of ∼5 μK has an initial velocity spread of σ v 3 cm/s. After preparing the atoms in the F = 1 , m F = 0 state using a sequence of near-resonant push beams and microwave π-pulses, typically 1 / 3 of the atoms remain. Using a velocity-selective Raman pulse with Ω eff / 2 π 10 kHz transfers a narrow velocity class ( σ v = Ω eff / k eff 0 . 4 cm/s) of atoms to the F = 2 , m F = 0 state, while the remaining atoms are removed—resulting in an 8-fold loss in atom number. Thus, in total, Raman AIs exhibit atom number loss factors on the order of ∼25. In this example, approximately N∼ 4 × 10 6 contribute to the interferometer. However, each atom can scatter several thousand photons during the resonant detection pulses, thereby ensuring an adequate signal-to-noise ratio. A measure of a Raman interferometer’s sensitivity is the so-called shot-noise or quantum-projection-noise limit [63], where the Poissonian fluctuations of atomic state measurements limit the minimum uncertainty of individual phase measurements to δ ϕ shot = 1 / N . Shot-noise limits of <1 mrad are considered to be state-of-the-art [64,65].
For the simple example of the Earth’s gravitational potential, where L = M v 2 / 2 M g z and r c ( t ) = ( z 0 + ( v 0 + k eff 2 M ) t 1 2 g t 2 ) z ^ , it is straightforward to show that the first two phase terms in Equation (2) vanish—leaving only the laser phase. Hence, if the Raman beams are aligned along z ^ such that k eff = k eff z ^ , the total phase shift is
Δ Φ tot = k eff g T 2 + Δ φ ,
with Δ φ = φ 1 2 φ 2 + φ 3 , which is usually used as a control parameter in the experiment to scan the interference fringes—enabling a direct measurement of the gravitational acceleration. Equation (4) illustrates the strong sensitivity of atom interferometers to inertial effects such as gravity. Since the phase shift scales as k eff T 2 , with a modest interrogation time of T ∼ 50 ms and k eff ∼ 1.6 × 10 7 rad/m for light at wavelength λ = 780 nm, the acceleration due to gravity induces a phase shift of Δ Φ tot 4 × 10 5 rad. Hence, with a phase uncertainty of 1 mrad, the single-shot sensitivity of the interferometer is ∼ 3 × 10 9 g . State-of-the-art cold-atom gravimeters have demonstrated precisions of 0.2 ppb after 1000 s of integration [47].
For measurements of the atomic recoil frequency with a Raman interferometer, typically the Ramsey-Bordé configuration is used [14,15,16,45]. In this case, the central π-pulse is replaced with two π / 2 -pulses separated by a time T , which has the effect of spatially separating the two output ports of the interferometer. The phase difference between interfering pathways is then Δ Φ tot ± = 2 ω k eff T ± k eff g T ( T + T ) , where the ± corresponds to the upper and lower output ports, respectively [66]. The phase shift due to gravity can be rejected by operating the interferometer in a conjugate mode where the upper and lower ports are detected simultaneously—yielding the sum of the two phases Δ Φ tot = 4 ω k eff T [67]. To increase the sensitivity to the recoil frequency ω k eff , large momentum transfer beam-splitters such as high-order Bragg transitions, have been used in place of two-photon Raman transitions [68,69,70]. This has the effect of replacing the effective wave with n k eff in the equations above, thus the recoil phase of the conjugate Ramsey-Bordé interferometers becomes 4 n 2 ω k eff T . Bloch oscillations have also been used to increase the common momentum transfer to the atoms between the two pairs of π / 2 -pulses [14,15]. In this case, the recoil frequency is measured from the phase shift between the two Ramsey fringe patterns created with the first and last pairs of π / 2 -pulses [6]. This phase shift is proportional to the number of photon momenta transferred to the atoms by the Bloch oscillations, where transfers as large as 1600 photon momenta have been demonstrated [14,71]. Presently, the state-of-the-art in terms of precision for a measurement of h / M is 1 . 2 × 10 9 [15].

1.2. Description of Grating-Echo-Type Interferometers

The grating-echo AI is a single-state Talbot-Lau interferometer [72,73,74], the principles of which can be understood on the basis of a plane-wave description of the two-pulse scheme shown in Figure 1b [32,54,55,75,76,77]. This AI relies on matter-wave interference produced by Kapitza-Dirac scattering of atoms by short, off-resonant standing wave (SW) pulses. Typically, the interferometer uses a sub-Doppler-cooled atomic sample with a momentum spread σ p k , initially prepared in a single hyperfine ground state F . Two SW pulses are applied to the sample separated by a time T. The first pulse excites a superposition of momentum states separated by integer multiples of q . The second excitation pulse further diffracts the momentum states, causing certain trajectories to interfere in the vicinity of t = 2 T , henceforth referred to as the “echo” time. This interference creates a spatial modulation in the atomic probability density with a phase that is proportional to inertial effects (i.e., ϕ a = q · a T 2 due to an acceleration a), and a contrast that is temporally modulated at a harmonic of the two-photon recoil frequency ω q = q 2 / 2 M . To measure the properties of this interference, a unique optical detection scheme is used—a traveling wave read-out pulse is applied the sample in the vicinity of the echo time when the density modulation is strongest. A certain spatial harmonic of this modulation satisfies the Bragg condition for scattering light in the backward direction (i.e., the harmonic with period λ / 2 ) as shown in Figure 1b. This back-scattered “echo signal” carries both the phase and contrast information about the atomic interference between certain classes of trajectories—namely those whose momenta differ by q at t = 2 T .
A common feature of the echo AI experiments described in this article is that the contrast of the density modulation (grating) is small, due to the relatively small atom-field coupling strength and the short pulse durations of the SW pulses. Consequently, the experiments are limited by the strength of the signal, which is defined by the reflectivity of the grating (≈ 0.2%). So although echo experiments do not experience the appreciable atom loss characteristic of Raman AIs, they require large atom numbers and high-contrast gratings to achieve appreciable signal strengths.
To understand how the atomic density grating comes about, we consider two overlapping momentum states F , n q and F , n q labelled by integers n and n . In comparison to Raman and Bragg interferometers [8,78], no atom optical “combiner” pulse is required to produce interference between two wavepacket trajectories since the momenta are in the same hyperfine ground state F . Spatial overlap is the only condition required to create an interference pattern, which can be described by
F , n q | F , n q e i ( n n ) q · r e i Δ ϕ n , n .
Here, Δ ϕ n , n is the phase difference between the wavepackets, which has contributions from the Doppler shift, atomic recoil and the SW laser phase, as we discuss below. If the integers n and n satisfy n = n ± 1 , then the real part of this interference is ∼ cos ( q · r + Δ ϕ n , n ± 1 ) . The density distribution follows this simple sinusoidal pattern, which exhibits a period of λ / 2 and hence satisfies the Bragg scattering condition for detection. In reality, the pair of SW pulses excite multiple interfering trajectories—each contributing its own spatial harmonic to the density distribution. To account for this multi-path interference, one must sum over all possible trajectories to arrive at the correct interference pattern. Although this can lead to extremely complex periodic structures in the atomic density [55,75], a simplification that can always be made is the fact that the read-out light will only scatter from the q-Fourier component of this structure.
We now give a detailed description of the plane-wave theory of grating-echo formation. Initially, the atomic wavefunction is in a hyperfine ground state labeled by total angular momentum F with momentum p, thus the wave function before the first SW pulse can be written as
ψ 0 ( p , t ) = F , p e i ( ω 0 + ω p ) t ,
where ω 0 is the internal energy of the ground state, and ω p = p 2 / 2 M is the frequency associated with the initial kinetic energy of the atom. Since the phase term e i ( ω 0 + ω p ) t is common to all diffracted momentum states, it is unimportant for interference and we henceforth ignore it. In the short pulse duration (i.e., Raman-Nath) regime, the interaction with the off-resonant SW pulse simply modulates the phase of the wavefunction as e i Ω eff τ j cos ( q · r + φ j ) ψ 0 , where Ω eff is the effective Rabi frequency of the light, τ j is the duration of the j th pulse and φ j is the phase of the standing wave (usually defined by the location of the node created by a retro-reflecting mirror). We can use the Jacobi-Anger expansion [79] to describe this modulation in a more convenient form
e i Ω eff τ j cos ( q · r + φ j ) = n = ( i ) n J n ( Ω eff τ j ) e i n ( q · r + φ j ) ,
where J n ( x ) is the n th order Bessel function of the first kind. Thus, after the first standing wave pulse applied at time t = t 1 , the wavefunction can be shown to be
ψ 1 ( p , t ) = n A n F , p + n q e i n φ 1 e i n q · v ( t t 1 ) e i n 2 ω q ( t t 1 ) ,
where A n = ( i ) n J n ( Ω eff τ 1 ) and v = p / M is the initial center-of-mass velocity of the atom. Here, we have used the fact that F , p e i n q · r = F , p + n q . After the second standing wave pulse at time t = t 2 = t 1 + T 21 , the wavefunction is given by
ψ 2 ( p , t ) = n , m A n B m F , p + ( n + m ) q e i ( n φ 1 + m φ 2 ) × e i q · v [ n T 21 + ( n + m ) ( t t 2 ) ] e i ω q [ n 2 T 21 + ( n + m ) 2 ( t t 2 ) ] ,
where T 21 = t 2 t 1 and B m = ( i ) m J m ( Ω eff τ 2 ) . To find the interference pattern at time t = t 1 + 2 T 21 , we compute the atomic density distribution ρ 2 = ψ 2 | ψ 2
ρ 2 ( p , t ) = n , m , n , m A n A n * B m B m * F , p + ( n + m ) q | F , p + ( n + m ) q e i [ ( n n ) φ 1 + ( m m ) φ 2 ] × e i q · v [ ( n n ) T 21 + ( n + m n m ) ( t t 2 ) ] e i ω q { ( n 2 n 2 ) T 21 + [ ( n + m ) 2 ( n + m ) 2 ] ( t t 2 ) } .
This time-dependent expression, although complex, has a simple interpretation. Each term in the sum is composed of three factors: (i) the complex amplitude factor A n A n * B m B m * which determines the relative strength of different interfering trajectories, (ii) the interference term F , p + ( n + m ) q | F , p + ( n + m ) q e i ( n + m n m ) q · r , which produces a modulation in the atomic density with spatial harmonic ( n + m n m ) q , and (iii) a series of phase factors that modify the phase of the density modulation due to the laser interaction ( e i [ ( n n ) φ 1 + ( m m ) φ 2 ] ), the Doppler shift ( e i ϕ D ( t ) ), and the atomic recoil ( e i ϕ q ( t ) ), where
ϕ D ( t ) = q · v [ ( n n ) T 21 + ( n + m n m ) ( t t 2 ) ] ,
ϕ q ( t ) = ω q { ( n 2 n 2 ) T 21 + [ ( n + m ) 2 ( n + m ) 2 ] ( t t 2 ) } .
The set of integers { n , m , n , m } label the momentum (in units of q ) transferred to the atom by the SW pulses, and represent a particular pair of interfering trajectories in Figure 1b. For instance, the integer labels corresponding to the trapezoidal trajectories of the Mach-Zehnder geometry shown in Figure 1a are { n , m , n , m } = { 1 , 1 , 0 , 1 } (Here, we interpret the unprimed integers { n , m } as the momenta transferred along the upper pathway, while the primed integers { n , m } correspond to the lower pathway of any two trajectories.). It follows that the interference between these trajectories produces a density modulation with a period of 2 π / | n + m n m | q = 2 π / q = λ / 2 , which is the ideal period for back-scattering the electric field of the read-out pulse at wavelength λ.
Since the velocity distribution of the sample is assumed to be broad compared to the scale of momentum transfer ( σ p k ), the macroscopic density grating produced in the sample is found by averaging the single-atom probability density (10) over this distribution of velocities. However, the velocity-dependent Doppler phase causes a strong dephasing effect on the grating at all times except certain “echo” times where this phase is zero [55,75,80]. One can show that these times must satisfy t = t 1 + ( 1 δ 1 / δ 2 ) T 21 , where δ 1 = n n and δ 2 = n + m n m . Here, we are concerned with only the first echo time at t 1 + 2 T 21 , which implies a ratio of δ 1 / δ 2 = 1 . The back-scattering detection method constrains δ 2 = ± 1 , thus we require δ 1 = 1 (i.e., n = n ± 1 ). Inserting this constraint into δ 2 yields m = m 2 . These two constraints define the class of trajectories that produce a macroscopic interference pattern at t = t 1 + 2 T 21 which can back-scatter light at wavelength λ. This interference pattern, which represents only a subset of the total density modulation given by Equation (10), can be shown to be
ρ ˜ 2 ( r , t ) s = 1 , 1 e i s q · r e i s ( φ 1 2 φ 2 ) e i s q · v ( t t 1 2 T 21 ) e i s 2 ω q ( t t 1 ) × n A n A n + s * e i 2 n s ω q ( t t 1 2 T 21 ) m B m B m 2 s * e i 2 m s ω q ( t t 1 T 21 ) .
We now average over the velocity distribution of the sample, which is assumed to be a Maxwell-Boltzmann distribution N ( v ) = 1 π 3 / 2 σ v 3 e ( v / σ v ) 2 with e 1 radius σ v
ρ ˜ 2 ( r , t ) e [ q σ v ( t t 1 2 T 21 ) / 2 ] 2 cos ( q · r + φ 1 2 φ 2 ) × J 1 2 Ω eff τ 1 sin ω q ( t t 1 2 T 21 ) J 2 2 Ω eff τ 2 sin ω q ( t t 1 T 21 ) .
Here, we have made use of the Bessel function identity α J α ( x ) J α + β ( x ) e i α ϕ = ( i ) β e i β ϕ J β ( 2 x sin ϕ ) [79] to simplify the sums over n and m in Equation (12). Two important features of the interference pattern are now clear. First, as a result of velocity dephasing, the grating only exhibits non-vanishing contrast for a timescale of 2 / q σ v 1 μs in the vicinity of the echo time. Second, the atomic recoil frequency, which initially appeared in the phase of the wavefunction, now affects only the contrast of the interference pattern. This feature of echo AIs alleviates the need for phase sensitivity in a recoil-sensitive experiment, since the effect can be measured in the back-scattered signal intensity. This type of AI has also been referred to as a “contrast” interferometer in the context of recoil measurements with ultra-cold atoms [81,82].
The final step is to compute from this macroscopic density the signal that is detected in the experiment by back-scattering the traveling wave read-out light. The physical mechanism that generates this light is elastic Rayleigh scattering from a spatial modulation of the sample’s refractive index that satisfies the Bragg condition [83,84,85]. This coherent scattering process results from a phase-matching condition along the Bragg angle. Whereas the intensity of diffuse atomic scattering scales as the number of scatters N, here the intensity scales as N 2 —a well-known feature of coherent Bragg scattering [84]. The drawback of this process is that, since it depends on a coherent superposition of momentum states, each atom scatters at most one photon before being projected into one of the two states. In comparison, the incoherent process of near-resonant fluorescence used in Raman interferometers allows one to scatter many photons per atom to increase the signal-to-noise ratio. This emphasizes the need for large atom numbers, low-sample temperatures, and high-contrast gratings to reach signal-to-noise ratios comparable with Raman AIs.
The macroscopic density grating described by Equation (13) acts as a linear reflector for light of wavelength λ [85,86]. Thus, a traveling-wave read-out field of E RO e i ( k · r ω t + φ 3 ) couples to the atomic grating and produces a back-scattered field given by
E BS ( t ) r ( t ) E RO e i ( k · r ω t + φ 3 ) ,
where r ( t ) is a time-dependent reflection coefficient [85], which depends on the detuning of the read-out light, and the contrast of the density modulation at spatial frequency q. Hence, for a fixed detuning, the back-scattered field is proportional to the probability density given by Equation (13)
E BS ( t ) E RO e i ( k · r + ω t ) e i ( φ 1 2 φ 2 + φ 3 ) e [ q σ v ( t t 1 2 T 21 ) / 2 ] 2 × J 1 2 Ω eff τ 1 sin ω q ( t t 1 2 T 21 ) J 2 2 Ω eff τ 2 sin ω q ( t t 1 T 21 ) .
This back-scattered field contains all the information about the atomic interference between momentum states differing by q . For instance, the time-integrated power of the back-scattered light (referred to as the echo energy) is a measure of the contrast produced by this interference. Experiments utilizing the two-pulse AI, where the echo energy is measured as a function of T 21 , are described in Refs. [32,54,55,56,57,76,77].
Similarly, the effect of gravity on the echo AI is to shift the phase of the atomic grating, which in turn causes a phase shift on the back-scattered electric field. In the same spirit as described in Section 1.1, the phase shift due to gravity can be computed solely by considering the interaction with the lasers. From Equation (15), the laser phase has the same form as for the Raman interferometer: Δ ϕ las = Δ φ φ 1 2 φ 2 + φ 3 . By replacing the each φ j with φ j + q · r c ( t j ) , where r c ( t ) = r 0 + v 0 t + 1 2 g t 2 is the center-of-mass trajectory of the atom under gravity, we find at t = t 1 + 2 T 21
Δ ϕ las = q · g T 21 2 + Δ φ .
Echo AI experiments that have demonstrated sensitivity to this gravitational phase shift are described in Refs. [30,54,55,58,87].

2. Measurements of Atomic Recoil Frequency

2.1. Introduction

There is an ongoing, international effort to develop precise, independent techniques for measuring the atomic fine structure constant, α—a dimensionless parameter that quantifies the strength of the electromagnetic force which lies at the heart of light-matter interactions. These measurements can be used to stringently test the theory of quantum electrodynamics (QED). Historically, two types of determinations of α have been carried out: (i) those that use other precisely measured quantities to determine α through challenging QED calculations [88,89], and (ii) those that are independent of QED, and depend on only the quantities appearing in the definition α e 2 / 2 ϵ 0 h c , where e is the elementary charge, ϵ 0 is the vacuum permittivity, h is Planck’s constant and c is the speed of light. Some examples of α determinations that require QED are the measurements of the anomalous magnetic moment of the electron (precise to 0.37 ppb) [90], and the fine structure intervals of helium (precise to 5 ppb) [91]. The most precise examples of QED-independent determinations are those based on measurements of the von Klitzing constant, R K = h / e 2 , using the quantum-Hall effect [92,93], and the ratio h / M using (i) Bloch oscillations in cold atoms [13,94] and (ii) atom interferometric techniques [6,12,14,15]. Within these examples, atom interferometry has emerged as a powerful tool because of its inherently high sensitivity to h / M , which can be related to α according to
α 2 = 2 R c h m e = 2 R c M m e h M .
Here, R is the Rydberg constant, m e is the electron mass, and M is the mass of the test atom. Since R is known to 6 parts in 10 12 , and the mass ratio M / m e is typically known to a few parts in 10 10 [95], the quantity that limits the precision of a determination of α using Equation (17) is the ratio h / M . The most precise measurement of this ratio was recently carried out with 87 Rb, where h / M was determined to 1 . 2 ppb after 15 hours of data acquisition [15]. The corresponding determination of α was precise to 0.66 ppb. Other interferometric techniques that have demonstrated high sensitivity to h / M include Refs. [16,81,82,96,97,98,99].
In this section, we describe recent improvements in measurements of the atomic recoil frequency ω q = q 2 / 2 M using echo AIs [54,55,57]. The appeal of this type of AI lies in it’s reduced sensitivity to common systematic effects, such as phase shifts due to the AC Stark or Zeeman effects, since it involves only a single internal state. In addition, since echo AIs rely on short-duration standing-wave pulses, only a single laser is required, and the large bandwidth of these pulses alleviates the need for velocity selection. Finally, as mentioned in Section 1.2, the signature of atomic recoil affects only the contrast of the interference pattern and is insensitive to low-frequency phase noise of the standing wave. Thus, in comparison to Raman interferometers, a phase-stable apparatus is not required to make high-precision measurements of ω q .
We have developed a “modified” three-pulse echo AI (This configuration is distinct from the “stimulated” three-pulse AI used for measurements of g that we present in Section 3.5.3.) which exhibits increased sensitivity to atomic recoil compared to the aforementioned two-pulse configuration [57]. This modified geometry has been described in previous work using the formalism of coherence functions [100] and a full quantum-mechanical treatment that accurately describes the fringe shape [54,56,100]. In the same articles, we also discussed connections to δ-kicked rotors and quantum chaos, as well as scaling laws that apply to excitation with multiple SW pulses that have been used in other work [101,102].

2.2. Description of the Modified Three-Pulse AI

For the modified three-pulse AI described in this section, an additional SW pulse is applied between the first two pulses at t = δ T , as shown in Figure 2b [54,56,100,101,102]. This pulse has the effect of diffracting the atom into higher-order momentum states that contribute additional harmonics of ω q to the temporal modulation of the macroscopic grating contrast. Intuitively, the third SW pulse acts as a phase mask analogous to the function of a multi-slit pattern in classical optics. More specifically, this pulse shifts the phase of the momentum states by η ω q δ T , where η is an integer that depends on the particular pathways that lead to interference at t = 2 T . Hence, varying the time of this pulse δ T is analogous to moving the slit pattern along the propagation axis of light—yielding periodic revivals of the contrast of the interference pattern. An example of a pair of low-order interfering trajectories created by this interferometer is shown in Figure 2b (We emphasize that only a small subset of the trajectories excited by the SW pulses will interfere at t = 2 T for an arbitrary third pulse time, δ T (i.e., trajectories which, when combined, exhibit a Doppler phase that is independent of δ T ). Specifically, the only momentum states contributing to the signal are those that differ by q after SW3 and after SW2.). When one accounts for all possible trajectories, the resulting signal consists of a series of narrow fringes separated by the recoil period, τ q = π / ω q (∼32 μs for rubidium), as a result of the interference between all excited momentum states differing by q .
When all relevant trajectories are summed over, it can be shown [56,100] that the resulting echo energy is modulated by J 0 [ 2 u 3 sin ( ω q δ T ) ] 2 , provided the third pulse area u 3 = Ω eff τ 3 1 . Here, J 0 ( x ) is the zeroth-order Bessel function of the first kind, Ω eff is the effective two-photon Rabi frequency, and τ 3 is the third SW pulse duration. Figure 2c illustrates the predicted dependence of the echo energy—that is, the energy in the back-scattered electric field—as a function of δ T . The sensitivity of this AI to ω q scales inversely with the time scale T over which the signal can be measured, and it scales in proportion to the width of the fringes. The advantage of using this AI over the two-pulse configuration is the ability to narrow the fringe width using the third pulse. Additionally, since T is fixed, the same number of atoms remain in the excitation beams at the time of detection—avoiding a loss of signal with increasing δ T due to effects like the thermal expansion of the sample. The fringe width is effectively determined by the width of the excited momentum distribution. By increasing the proportion of high-order momentum states (and thus the proportion of high-order harmonics of ω q ) that contribute to the signal, the fringes become more sharply defined. The excitation is controlled by the interaction strength and duration of the third SW pulse. It can be shown that for small pulse durations (i.e., τ 3 ( Ω eff ω q ) 1 / 2 ) the full-width at half-maximum (FWHM) of the fringe scales as 1 / u 3 , as illustrated in Figure 2c [54,56,100].

2.3. Experimental Setup

As described in Refs. [32,56,57], two major improvements to our experiment have enabled us to reach time scales of T 60 ms: (i) utilizing a non-magnetizable glass vacuum system, which reduced decoherence effects related to inhomogeneous B-fields [58] and improved the molasses cooling of the sample, and (ii) using large-diameter, chirped excitation beams, which increased the transit time of the atoms in the beam and compensated for the Doppler shift due to gravity (The non-uniform magnetic field produced by a stainless-steel vacuum chamber, and the gravity-induced Doppler shift limited previous experiments to T 10 ms [76,100].).
The experiment utilizes a laser-cooled sample of rubidium typically containing ∼5 × 10 9 atoms at temperatures of T 5 μK. Either 85 Rb or 87 Rb atoms are loaded into a six-beam, vapor-loaded magneto-optical trap (MOT) in 250 ms. Prior to the AI experiment, the sample is prepared in the upper hyperfine atomic ground state ( 5 S 1 / 2 F = 3 for 85 Rb or F = 2 for 87 Rb). The light for the AI is derived from a Ti:sapphire laser (linewidth ∼1 MHz) that is locked above the D2 cycling transition using Doppler-free saturated absorption spectroscopy. A network of acousto-optic modulators (AOMs) is used to generate the frequencies necessary for the AI excitation and the read-out beams. The read-out light is detuned to the blue of the cycling transition by Δ RO = 40 MHz, which optimizes the back-scattered signal intensity for our sample size and density [56]. The AI beams are detuned by Δ AI = 220 MHz, and a frequency chirp of δ ( t ) = g t / λ is added to (subtracted from) the downward- (upward-) traveling component of the SW pulses. This compensates the Doppler shift induced by the falling atoms, and ensures that the beams remain resonant with the two-photon transition [32,56,57]. A “gate” AOM is used upstream of the AI AOMs as both a frequency shifter and a high-speed shutter to reduce the amount of stray light in the experiment. All RF sources and digital-delay generators used to define the pulse timing for the AI are externally referenced to a 10 MHz rubidium clock.
The AI beams are coupled into two AR-coated, single-mode optical fibers and aligned through the sample, as shown in Figure 3a. At the output of the fibers, the beams are expanded to a e 2 diameter of d 1 . 7 cm and are circularly polarized in the same sense by a pair of λ / 4 wave plates. The timing sequence for the experiment is illustrated in Figure 3b. A mechanical shutter on the upper platform closes before the read-out pulse in order to block the back-scatter of stray read-out light produced by various optical elements. This light would otherwise interfere with the coherent signal from the atoms. A gated photo-multiplier tube (PMT, 8 × 10 5 W/V at 780 nm, noise equivalent power 100 nW) is used to detect the power in the back-scattered field. Figure 3c shows an example of the echo signal recorded by the PMT averaged over 16 repetitions of the two-pulse AI. This signal is converted to units of optical power and numerically integrated to obtain a quantity which we term the echo energy. This quantity is proportional to the contrast of the atomic density modulation and the intensity of read-out light incident on the atoms. As a result, the signal is sensitive to both atom number fluctuations and photon number shot noise. This is a drawback compared to fluorescence detection techniques, where the optical transition can be saturated and is therefore less sensitive to photon shot noise [65]. In these experiments, we typically observe a noise floor of 0.1 pJ per shot, or ∼0.025 pJ after averaging over 16 repetitions, which was dominated primarily by the NEP of the PMT.

2.4. Results

Measurements of ω q were obtained using the modified three-pulse AI by measuring the echo energy as a function of the third pulse time, δ T , as illustrated in Figure 4a. This figure shows a measurement of ω q in 85 Rb on a time scale of T 36 . 7 ms, which was acquired in 15 minutes. Clearly, the shape of the fringes does not resemble that predicted by the theory shown in Figure 2(c). This is due to the contribution from each of the magnetic sub-levels in the F = 3 ground state of 85 Rb, which tend to smear out the higher harmonics in the signal as a result of their different optical coupling strengths. Furthermore, the presence of additional, nearby excited states ( F = 2 and 3 in the case of 85 Rb) produces an asymmetry in the fringe lineshape [56]. This effect is reduced in 87 Rb because the frequency difference between neighboring excited states is larger. To measure ω q , the data are fit to a phenomenological model that consists of a periodic sum of exponentially-modified Gaussian functions
F ( δ T ; τ q ) = l A l exp 1 2 σ l υ 2 + δ T l τ q υ erfc 1 2 σ l υ + δ T l τ q σ l ,
and the recoil frequency, ω q = π / τ q , is extracted from the fit. In this model, erfc ( x ) is the complementary error function, and the parameter υ, which determines the amount of asymmetry in the lineshape, is the same for all fringes. The fit to the data shown in Figure 4a yielded a reduced chi-squared of χ 2 / dof = 0 . 51 for dof = 300 degrees of freedom. This corresponds to a relative statistical precision of ∼180 ppb in ω q —representing a factor of ∼9 improvement over previous work [100].
To demonstrate the long-term statistical sensitivity of the interferometer, 82 independent measurements of ω q in 87 Rb were recorded (see Figure 4b) while holding all experimental parameters fixed to the extent possible. Here, ω q is determined from a weighted average over all individual measurements, where the points are weighted inversely proportional to the square of their statistical errors. The mean value shown in Figure 4b, which has not been corrected for systematic effects, is found with a statistical uncertainty of 37 ppb as determined by the standard deviation of the mean. An autocorrelation analysis of these data indicate that they are correlated at the 20% level with measurements taken at a previous time. This is attributed to slowly-varying environmental conditions, such as temperature and magnetic field, over the 14 hours of data acquisition.

2.5. Discussion of Systematic Effects

We have investigated systematic effects on the measurement of ω q related to the angle between excitation beams, the refractive indices of the sample and the background Rb vapor, light shifts, Zeeman shifts, B-field curvature and the SW pulse durations [56]. The total systematic uncertainty in this measurement is estimated to be ∼5.7 parts per million (ppm), and is dominated by two effects: (i) the refractive index of the sample, and (ii) the curvature of the B-field that the atoms experience as they fall under gravity. We now discuss these two effects in detail.
The refractive index of the atomic sample affects the wave vector of the excitation beams, since a photon in a dispersive medium acts as if it has momentum n k , where n is the index of refraction [104]. For near-resonant light, the index becomes a function of both the density of the medium, ρ, and the detuning of the applied light from the atomic resonance, Δ AI . The systematic effect on the recoil frequency due to the refractive index can be expressed as ω q ( ρ , Δ AI ) = ω q ( 0 ) n 2 ( ρ , Δ AI ) , where ω q ( 0 ) is the recoil frequency in the absence of systematic effects. The index of refraction can be computed from the electric susceptibility and the light-induced polarization of the medium [104]. Taking into account the level structure of the atom, it can be shown that [56]:
n ( ρ , Δ H G ) = 1 ρ ϵ 0 Γ H μ H G 2 Δ H G / Γ 1 + ( Δ H G / Γ ) 2 .
Here, Δ H G ω ( ω H ω G ) is the atom-field detuning between the ground and excited manifolds, g , G and e , H , for laser frequency ω, G (H) is a quantum number representing the total angular momentum of a particular ground (excited) manifold, and μ H G is the reduced dipole matrix element for transitions between those manifolds [74]. The root-mean-squared density of the cold sample at the time of trap release was estimated to be ( 4 . 1 ± 1 . 2 ) × 10 10 atom/cm 3 based on time-of-flight images [56]. Hence, we estimate a refractive-index-induced shift in ω q of 10 . 5 ± 3 . 0 ppm at a detuning of Δ AI = 220 MHz.
The other dominating systematic effect is due to the inhomogeneity of the magnetic field sampled by the atoms during the total interrogation time ( 2 T ∼ 91 ms) of the interferometer. This field primarily originates from nearby ferromagnetic material, such as an ion pump magnet and a glass-to-metal adaptor, and from the set of quadrupole coils we use to cancel the residual field in the vicinity of the MOT [56]. At the end of the optical molasses cooling phase, the atoms are distributed roughly equally in population among the magnetic sub-levels of the upper ground state F = 2 . A spatially-varying B-field with a non-zero curvature (i.e., β 2 2 B / z 2 0 ) has a parasitic impact on the interferometer as a result of a position-dependent force on the m F 0 states similar to a harmonic oscillator (We have also considered the effect of a constant background B-field, and found that it produces a small systematic effect on ω q (∼7.5 ppb per Gauss of residual B-field) as a result of the distribution of sub-level populations. Similarly, linear magnetic gradients do not affect the measurement of ω q using the echo AI [32].). For each momentum state trajectory, the atom samples a different region of space and experiences a different acceleration than that of a neighbouring trajectory. Since the momentum of each trajectory is differentially modified between excitation pulses, the interference for each class of trajectories occurs at a slightly different time—causing both a systematic shift of ω q and a loss of interference contrast. For a given state F , m F , the corresponding systematic correction to ω q is
ω q ( β 2 , T ) = ω q ( 0 ) 1 + 2 3 m F g F μ B β 2 M T 2 .
where g F is a Landé g-factor, and μ B is the Bohr magneton. In the worst case scenario, the systematic shift in the recoil frequency is dominated by the state F = 2 , m F = 2 . Based on measurements of the B-field in the vicinity of the atoms, we estimate β 2 0 . 1 mG/cm 2 = 10 4 T/m 2 . Thus, for T = 50 ms we estimate a relative shift of ∼11 ppm. A more detailed analysis [56,57], which accounts of the distribution of magnetic sub-level populations, yields a more realistic estimate of 6 . 3 ± 4 . 4 ppm.
Including other minor systematic effects, such as the relative beam angle ( 15 ± 8 ppb), the refractive index of the background vapor ( 52 ± 31 ppb), light shifts due to the interferometer beams and the saturated absorption setup ( 55 ± 2 ppb), and the finite SW pulse durations ( 0 ± 2 ppm), we estimate the total systematic shift on our measurement of ω q to be 4 . 3 ± 5 . 7 ppm [56]. Correcting for this shift, we find that our measurement ω q = 94 . 77400 ( 54 ) rad/ms is within 1.6 ppm of the expected value of ω q ( 0 ) = 94 . 77384783 ( 51 ) rad/ms, as derived from the most precise measurement of h / M [15]. The combined statistical (37 ppb) and systematic (5.7 ppm) uncertainties of our measurement are enough to account for this discrepancy.
We now discuss techniques for reducing the aforementioned systematic effects. Equation (19) for the refractive index suggests that the relative correction to ω q can be reduced by decreasing the sample density, ρ, or by increasing the excitation beam detuning, Δ AI . However, the current configuration of the AI relies on a large number of atoms to achieve a sufficient signal-to-noise ratio. Thus, a decrease in the sample density leads to a reduction in the signal size. Furthermore, the sensitivity of the three-pulse AI relies on a relatively strong atom-field coupling in order to excite many orders of momentum states. An increase in the excitation beam detuning without a corresponding increase in the field intensity leads to a reduction in the sensitivity to ω q . The refractive index systematic could be reduced by a factor of 10 3 by a 10-fold reduction in the rms density of the sample, accompanied by a 100-fold increase in the detuning. This would require an increase in the excitation field intensity by a factor of 100 (corresponding to 10 W/cm 2 ) in order to retain the same sensitivity to ω q . A more promising way forward is to utilize the frequency-dependence of the refractive index, which exhibits “magic” frequencies where the systematic shift cancels, as shown in Figure 5. These frequency are located between two excited state manifolds, where the dispersive corrections to the refractive index due to each state have the same magnitude but opposite sign. Since these frequencies are independent of the density ρ, they are ideal for cancelling the refractive index shift due to both the cold-atom sample and the background vapor [56]. Using this feature of the refractive index, one can avoid both reducing the sample density and using high-intensity excitation beams.
The systematic shift due to the B-field curvature can be significantly reduced by preparing the sample in the magnetically “insensitive” state F = 2 , m F = 0 (This can be achieved by employing a combination of optical pumping, resonant push beams and microwave π-pulses resonant with the F = 1 , m F = 0 F = 2 , m F = 0 transition, which the standard technique in Raman interferometers and atomic clocks.). Then, any systematics due to the B-field would originate from the second-order Zeeman effect which shifts all sub-levels the F = 2 hyperfine manifold by a frequency 1 2 K B 2 , where K 575 Hz/G 2 is the Zeeman shift of the clock transition in 87 Rb [103]. Compared to the first-order Zeeman shift of m F g F μ B 1 . 4 MHz/G for the state F = 2 , m F = 2 , the corresponding force on the atom is reduced by many orders of magnitude for the same B-field strength. Utilizing only m F = 0 atoms in the interferometer will have the added benefit of significantly reducing decoherence due to the B-field curvature—enabling an increase in T and a corresponding reduction in the statistical error of each measurement. Under current conditions, we have achieved a maximum time scale of T ∼ 65 ms. However, previous studies indicate that the transit time of the atoms in the excitation beams is ∼270 ms [32], suggesting that T can be as large as ∼135 ms before the temperature of the sample becomes the limiting factor.
Another avenue for improvement to increase the signal-to-noise ratio in the experiment, which directly affects the sensitivity to ω q . As mentioned in Section 1.2, the macroscopic grating behaves as a linear reflector for an optical field of wavelength λ with a complex reflectivity r [85,86]. Using measurements of the energy in the back-scattered echo signal and the optical power in the read-out pulse, we estimate a reflection coefficient of R = | r | 2 0 . 001 under typical experimental conditions. The reflectivity can potentially be increased by pre-loading the sample in an optical lattice such that the initial spatial distribution has a significant λ / 2 -periodic component [106]. Experimental studies of MOTs loaded into an intense, off-resonant optical lattice have shown that the reflection coefficient of the light that is Bragg-scattered off the resulting atomic grating can be as large as R 0 . 8 [107]. This motivates the pursuit of high-contrast grating production using a far-detuned lattice pulse that precedes the AI excitations. Numerical simulations of the reflection coefficient from the grating produced by a lattice-loaded sample indicate that a 100-fold increase in the back-scattered signal is feasible. Such an endeavor would require an apparatus with good stability and control of the phase of the SW fields (i) to effectively channel atoms into the nodes of the lattice potential without significant heating, and (ii) to match the phases of the excitation and lattice fields.
By implementing these improvements to the echo AI, we anticipate that a future round of recoil measurements will yield results with both statistical and systematic uncertainties at competitive levels.

3. Measurements of Gravitational Acceleration

3.1. Overview of Gravity Measurements

Interest in precise measurements of the gravitational acceleration g have been stimulated in part by the connection of such measurements to the determination of the universal gravitational constant G [18,19,20,108] and the possibility of the variation of the gravitational force on small-length scales [109]. Since these measurements can be designed to measure the absolute value of g or relative changes due to temporal effects such as tides and positional variations due to changes in density, gravimeters have played a ubiquitous role in the exploration of natural resources by detecting characteristic density profiles associated with minerals, petroleum, and natural gas. An important practical consideration is the ability of these sensors to provide a non-invasive technique for exploration in wide area (air, sea, or submersible) mineral assays involving environmentally-sensitive areas. Other applications include borehole mapping for verifying properties of rocks, determination of bulk density for the detection of cavities, and tidal forecasts. The most precise relative measurements of g are derived from superconducting quantum interference-based devices (SQUIDs) [110,111], whereas absolute measurements of g based on an optical Mach-Zehnder interferometer [112] can achieve an absolute accuracy of 1 ppb in an integration time of 20 minutes. This sensor relies on recording the chirped accumulation of fringes when a corner-cube retro-reflector on one arm of the interferometer falls through a height 0.3 m.
Interest in cold-atom-based interferometers began with the path breaking experiments in Refs. [8,113], which relied on a Raman interferometer to achieve a statistical precision of 3 ppm in a measurement time of 1000 s. Raman AIs have also obtained the most sensitive atom-based measurements of g. To select some examples, Ref. [114] achieved a statistical precision of 1.3 ppb in 75 s of data acquisition, whereas Ref. [22] included a detailed study of systematic effects and reported a statistical precision of 3 ppb over 1 min of integration. A key feature of both experiments was the active vibration stabilization of the inertial reference frame (i.e., the surface of the retro-reflection mirror) with respect to which the measurements were carried out. Additionally, in these examples atoms were launched in a 50 cm atomic fountain to obtain a free-fall time of over 300 ms. More recently, a Raman AI with a 6.5 m drop zone achieved an inferred single-shot sensitivity of 7 × 10 12 g [36,37]. The Raman AI has also been developed to realize the best atom-based measurements of gravity gradients [23,24] and rotations [27,28,29,36,115]. As a consequence, this AI has been the preferred configuration for remote sensing applications [11,47,48,52,64,116,117,118,119,120,121].
Alternative techniques have also demonstrated competitive measurements of g. Experiments relying on Bloch oscillations report statistical precisions of 50–200 ppb after a few minutes of integration time [51,53], and 220 ppb of total systematic uncertainty after a few hours in the case of Ref. [122] (This latter measurement relies on the precision of Planck’s constant h, which is presently known to 12 ppb [95].). Additionally, a single-state Mach-Zehnder interferometer involving Bragg transitions reported a sensitivity of 2.7 ppb after 1000 seconds of data acquisition using a drop height of 20 cm and passive vibration stabilization [52]. More recently, the same group demonstrated a Bragg-pulse gravimeter using a Bose-condensed sample which yielded an asymptotic uncertainty of 2.1 ppb [33]. The echo AI described in Ref. [58], which we review in this article, achieved a statistical precision of 75 ppb in one hour using a drop height of 1 cm and an apparatus in which only key components were passively vibration stabilized.
We now review measurements of g using echo AIs that rely on samples of laser-cooled rubidium atoms released from a MOT. The theoretical background and earlier results are described in Refs. [30,32,54,55,58,87,123].

3.2. Description of Echo AI Techniques for Measuring g

We first provide a discussion of the physical principles of AI configurations used for measurements of g that are based on the earlier theoretical description. Figure 6a represents the recoil diagram which shows displacements of centre-of-mass trajectories of wavepackets for momentum states for the two-pulse configuration of the AI based on a billiard ball model [124,125,126]. A sample of laser-cooled 85 Rb or 87 Rb atoms is excited along the vertical by two blue-detuned standing wave (SW) pulses separated by a time T 21 . Each SW pulse is composed of two traveling wave components, each carrying a wave vector with wavenumber k = 2 π / λ . Atoms in each of the magnetic sublevels of the F = 3 ground state in 85 Rb or the F = 2 ground state in 87 Rb are diffracted into a superposition of momentum states separated by q at t = 0 , where q = 2 k . This process involves the absorption of a photon from one travelling wave component of the standing wave and stimulated emission along the counter-propagating traveling wave component. The durations of the SW pulses are sufficiently short that they meet the Raman-Nath criterion [127,128], where the displacement of the atoms due to the momentum transfer from standing wave pulses is small compared to the spacing of the quasi-sinusoidal standing wave potential. For counter-propagating traveling wave components, the wavelength of the potential is λ / 2 . The classical description of the effect of the standing wave interaction is that the atoms are focused toward the nodes of the standing wave potential. The focusing of atoms into the nodes produces a one-dimensional density grating with a period of λ / 2 . In the quantum mechanical description, it is the interference between momentum states that produces this one-dimensional density grating. In this latter description, the atomic wave function develops a recoil modulation on a time scale τ q = π / ω q 32 μs, where ω q = q 2 / 2 M is frequency associated with the two-photon recoil energy of the atom. The velocity distribution of the cold sample along the SW axis causes the grating to dephase on a much shorter time scale τ coh = 1 / k u , where u = 2 k B T / M is the 1 / e width of the velocity distribution and T is the temperature of the laser-cooled sample. For typical sample temperatures of ∼20 μK, the dephasing time scale is ∼2 μs. A long time T 21 after the grating has dephased, a second SW pulse is applied to diffract the momentum states. The effect of this SW pulse is to cause the momentum states separated by q to rephase at the echo time 2 T 21 . Momentum state interference produces a maximum contrast in the density grating just before and just after the echo time. The rephasing is reminiscent of a two-pulse photon echo experiment [129] that involves a superposition of ground and excited states. The echo technique is a general method of cancelling Doppler dephasing in an atomic gas. In echo atom interferometry, this technique has been extended to ground states so that velocity selection is not required. The effect of cooling the sample is simply to ensure that the time scale of the experiment is suitably long. Under ideal conditions, the experimental time scale is limited by the transit time of atoms across the excitation beam.
Since the atoms are in the ground state, it is necessary to apply a near-resonant, travelling-wave read-out pulse in the vicinity of the echo time to detect the contrast and phase of the re-phased density grating. The periodic array of atoms formed at the echo time coherently back-scatters the read-out pulse. This process is known as Bragg scattering. The grating spacing of λ / 2 causes a total path difference change of λ for light reflecting from adjacent planes of the grating. This effect produces constructive interference since the phase difference of reflections from adjacent planes is 2 π . In this case, the wavelength of the back-scattered light is matched with the Fourier component of the density grating with spacing λ / 2 . The back-scattered signal due to the read-out pulse is called the echo signal. To determine g, the phase of this coherent signal, which scales as q g T 21 2 , is measured as a function of the pulse separation, T 21 .
The read-out light is back-scattered as a consequence of the law of conservation of momentum. If an incoming read-out photon is backscattered, then the total momentum delivered to the sample is q . This momentum transfer allows the two arms of the interferometer that differ by q to recombine. This action of the read-out pulse also creates a coherent superposition of ground and excited states throughout the sample. The total radiation pattern from this system of dipole radiators is phase-matched only along the backward direction. The experiment measures the phase of the echo signal with respect to an inertial frame of reference defined by an optical local oscillator (LO) with a frequency ω LO . A convenient reference frame is defined by the nodal point of a standing wave, such as a reflecting surface.
The back-scattered light at frequency ω AI is detected as a beat note at frequency | ω LO ω AI | using a heterodyne technique [55]. Although the quadratic accumulation of phase as a function of T 21 is appealing for a precision measurement of g (among various interferometer geometries, this configuration also encloses the largest space-time area), the signal amplitude exhibits recoil modulation as well as a chirped frequency, resulting in the need for a complicated fit function to extract g. The signal from the two pulse AI is analogous to the interference fringes recorded by the falling corner-cube optical interferometer discussed in Ref. [112].
An alternate configuration involves a three-pulse stimulated echo AI [32,75,123,130], as shown in Figure 6b. Here, the first SW pulse creates a superposition of momentum states separated by q . A second SW pulse applied at t = T 21 produces momentum states that are co-propagating at a fixed spatial separation with the same momentum during a central time window of duration T 32 . A third SW pulse applied at t = T 21 + T 32 causes the co-propagating states to interfere at the echo time t = 2 T 21 + T 32 , resulting in a density grating. Just like the two-pulse AI, the grating formation is associated with interference of momentum states separated by q . The signal amplitude as a function of pulse separation T 32 shows no recoil modulation and exhibits a fixed angular frequency q g T 21 due to the velocity g T 21 acquired by the atoms during the time interval T 21 in the presence of gravity. The period of the signal amplitude is given by τ v = λ / 2 g T 21 .
The constant modulation period improves the quality of the fits to the data, thereby resulting in improved statistical precision in measurements of g. For our experimental conditions, in which the setup was partially shielded from vibrations, the stimulated echo AI proved particularly useful. This is because the time window T 21 can be made small compared to the time scale over which vibrations cause mirror positions to become uncorrelated, while T 32 can be relatively large, thereby realizing a larger total time scale than the two-pulse AI. This feature is due to the co-propagating momentum states that have a constant spatial separation during the time window T 32 . The disadvantages include the reduction of the signal amplitude due to the additional SW pulse and the inherent sensitivity to any initial velocity along the SW axis.

3.3. Experimental Setup

Figure 7a shows a block diagram of the experimental setup. The details are described in [58]. A Ti:Sapphire laser is used to generate light for atom trapping and interferometry using a chain of acousto-optic modulators (AOMs) that serve as frequency shifters and amplitude modulators. All these elements are placed on a pneumatically-supported optical table. Light from these AOMs is transported to the atom trap using angle-cleaved, anti-reflection (AR) coated optical fibers.
The experimental setup for atom trapping and atom interferometry is shown in Figure 7b. The vacuum chamber used for atom trapping is made of 316 L stainless steel and it is anchored to an optical table mounted on pneumatic vibration isolators. The chamber is maintained at 5 × 10 9 Torr by an ion pump with a pumping speed of 270 L/s located 1 m away to reduce ambient magnetic fields. The chamber is surrounded by three pairs of magnetic field and gradient cancelling coils. A separate set of tapered coils wound on the chamber provides the magnetic gradient for atom trapping. The trapping optics, vacuum chamber, anti-Helmholtz and cancellation coils, and ion pump are supported by the optical table. The MOT is loaded from background vapor, with approximately 5 × 10 8 atoms loaded in 1 second. Time-of-flight charge-coupled device (CCD) camera images of atoms released after molasses cooling [131] show that the typical sample temperature is 20 μK.
The fiber-coupled beam used for atom interferometry identified in Figure 7a is aligned through a single-pass AOM operating at 250 MHz, as shown in Figure 7b. The circularly-polarized diffracted beam from this AOM, which is directed along the vertical and used for excitation of atoms, is detuned by Δ∼ 55 MHz above resonance [76]. This beam is retro-reflected through the atom cloud by a corner-cube reflector to produce standing wave excitation. The undiffracted beam, with a frequency of ω 0 + 305 MHz is spatially separated from the excitation beam by 2.5 cm. It is aligned through the same optical elements as the excitation beam to minimize the impact of relative phase changes due to vibrations and serves as an LO. The LO is physically displaced upon reflection by the corner-cube. The background light entering the apparatus during the AI pulse sequence is minimized by pulsing the gate AOM in Figure 7a only when the AI AOM is turned on. The excitation and LO beams, combined on a beam splitter and a balanced heterodyne detector with two oppositely-biased Si photodiodes with rise-times of 1 ns, are used to record a beat signal at a frequency ω RF = | ω AI ω LO | = 250 MHz. During the read-out pulse, the retro-reflection of the excitation beam is blocked by a mechanical shutter with an open/close time of 1 ms [132].
The corner-cube reflector, AI AOM, balanced detector, and related optics are anchored to a vibration isolation platform with a resonance frequency of 1 Hz, which rests on the pneumatically-supported optical table, as shown in Figure 7b. The optical table is effective in suppressing vibration frequencies above 100 Hz, whereas the vibration isolation platform suppresses frequencies in the range of 1–100 Hz. The mechanical shutter is separately anchored to the ceiling of the laboratory to reduce vibrational coupling. In this setup, only critical components are passively isolated with the vibration isolation platform.
Digital delay generators with time bases controlled by a 10 MHz signal from a rubidium clock (Allan variance of 10 12 in 100 seconds) are used to produce RF pulses with an on/off contrast of 90 dB to drive the AOMs. The time delays of optical pulses are controlled with a precision of 50 ps. The read-out pulse intensity is comparable to the saturation intensity of Rb atoms so that the entire echo signal envelope can be recorded without appreciably decohering the signal. This signal, which is measured as a 250 MHz beat note, is recorded on an oscilloscope with an analog bandwidth of 3.5 GHz and mixed down to DC using the RF oscillator driving the AI AOM to produce the in-phase ( E 0 cos ( ϕ g ) ), and in-quadrature ( E 0 sin ( ϕ g ) ) components of the back-scattered electric field. While the atom trap is loaded, an attenuated excitation beam is turned on to record a 250 MHz beat note. This measurement re-initializes the RF phase used to mix the signal down to DC at the beginning of each repetition of the experiment and ensures that the relative phases between the excitation beam and the LO are the same at the start of the experiment. Although the LO and AI beams are strongly correlated at the beginning of the experiment, the phase uncertainty progressively increases with the time scale of the experiment and it cannot be corrected mainly because the motion of the corner-cube reflector is not measured. The typical repetition rate of the experiment varies between 0.8–3 Hz.

3.4. Theory

We now review the theoretical description of the signal shapes and characteristics for both two- and three-pulse stimulated AI configurations using simplified equations that apply to an atomic system with a single magnetic ground state sublevel as in Refs. [32,58,87]. Here, g represents a constant gravitational acceleration along the axis of SW excitation.

3.4.1. Two-Pulse AI

The backscattered electric field due to the readout pulse for the two-pulse AI can be written as
E g ( 2 ) = E 0 ( 2 ) e i ϕ g ( 2 ) ,
where E 0 ( 2 ) is the electric field amplitude and ϕ g ( 2 ) is the gravitational phase. The electric field amplitude for the two-pulse AI can be shown to be
E 0 ( 2 ) E RO e ( Δ t / τ coh ) 2 J 1 2 θ 1 sin ( ω q Δ t ) J 2 2 θ 2 sin ω q ( T 21 + Δ t ) e t echo / τ decay ,
where E RO is the electric field amplitude of the read-out pulse, J n ( x ) is the n t h -order Bessel function of the first kind, θ 1 and θ 2 are the pulse areas of the first and second SW pulses respectively, Δ t = t 2 T 21 is the time relative to the echo time t echo = 2 T 21 , and ω q = q 2 / ( 2 M ) is the two-photon recoil frequency. Here, τ coh = 1 / k u is the coherence time due to Doppler dephasing that defines the temporal width of the signal shown in Figure 8a, where u = 2 k B T / M is the 1 / e width of the one-dimensional velocity distribution along the excitation beams and T is the sample temperature. The last term in Equation (22) represents a phenomenological decay, with a time constant τ decay that models signal loss due to decoherence mechanisms in the experiment (e.g., from spontaneous emission and the spatial curvature of the ambient magnetic field), as well as the transit time of cold atoms through the interaction zone defined by the excitation beams.
In the presence of gravity, the space-time area enclosed by the interferometer determines the phase accumulation and for the two-pulse AI [30,32,54,58,87], the phase is given by
ϕ g ( 2 ) = q · g ( T 21 2 + 2 T 21 Δ t + Δ t 2 / 2 ) .
We decouple the expression for the echo signal into two parts that are dependent on the time scales T 21 and Δ t to explain its characteristics. Accordingly, the two-pulse signal becomes
E g ( 2 ) = E D ( 2 ) ( Δ t ) E AI ( 2 ) ( T 21 ) e i ϕ D ( 2 ) ( Δ t ) e i ϕ AI ( 2 ) ( T 21 ) ,
where
ϕ D ( 2 ) ( Δ t ) = q · g ( 2 T 21 Δ t + Δ t 2 / 2 )
is the Doppler phase which is dependent on Δ t , and
ϕ AI ( 2 ) ( T 21 ) = q · g T 21 2
is the AI phase which is dependent only on T 21 . The Doppler component of the electric field amplitude is given by
E D ( 2 ) ( Δ t ) E RO e ( Δ t / τ coh ) 2 J 1 2 θ 1 sin ( ω q Δ t ) ,
and in the limit Δ t T 21 , the AI electric field amplitude is given by
E AI ( 2 ) ( T 21 ) = J 2 2 θ 2 sin ( ω q T 21 ) e t echo / τ decay .
The measurement of gravity is based on detecting the amplitude and phase of the back-scattered light from the atomic grating (which has a frequency ω AI and phase ϕ AI ) with reference to an optical LO, which has a fixed frequency ω LO and phase ϕ LO . The signal is recorded as a beat note at the frequency ω AI ω LO and with a phase difference ϕ signal = ϕ AI ϕ LO . The phase shifts associated with the atoms are sensitive to optical phase shifts of the SW pulses and the LO due to the environment. The total signal amplitude can be expressed in terms of the in-phase and in-quadrature components E 0 ( 2 ) cos ( ϕ g ( 2 ) ) and E 0 ( 2 ) sin ( ϕ g ( 2 ) ) as
E 0 ( 2 ) = 1 2 E 0 ( 2 ) cos ( ϕ g ( 2 ) ) 2 + E 0 ( 2 ) sin ( ϕ g ( 2 ) ) 2 1 / 2 .
The recoil modulation and signal decay terms can be removed from the in-phase and in-quadrature components of the back-scattered field amplitude by normalizing with respect to E 0 ( 2 ) . We are then left with cos ( ϕ g ( 2 ) ) and sin ( ϕ g ( 2 ) ) as the two components of the signal.
The dashed lines in Figure 8a show the Doppler electric field amplitude as predicted by Equation (27). Here, a convenient T 21 was chosen to maximize the recoil modulated signal, modeled by Equation (28). The solid red line shows gravity-induced oscillations within the echo envelope, as predicted by E D ( 2 ) e i ϕ D ( 2 ) . The oscillations are attributed to the free fall of atoms through a grating spacing of λ / 2 , which results in a phase increment of 2 π . This effect can also be described as a Doppler shift of the backscattered field due to the falling grating.
The solid red line in Figure 8c shows the predicted in-phase component of the signal amplitude for the two-pulse AI as a function of T 21 given by E AI ( 2 ) e i ϕ AI ( 2 ) (see Equations (26) and (28)). The recoil modulation its readily apparent and the frequency-chirped oscillations due to gravity are illustrated by setting g = 980 m/s 2 . The dashed black line shows the recoil-modulated total signal amplitude E AI ( 2 ) .

3.4.2. Three-Pulse AI

Based on Refs. [32,58,87], the backscattered electric field for the three-pulse stimulated echo AI shown in Figure 6b can be written as
E g ( 3 ) = E 0 ( 3 ) e i ϕ g ( 3 ) ,
where E 0 ( 3 ) is the electric field amplitude and ϕ g ( 3 ) is the gravitational phase. The electric field amplitude can be shown to be:
E 0 ( 3 ) E RO e ( Δ t / τ coh ) 2 J 1 2 θ 1 sin ( ω q Δ t ) × J 1 2 θ 2 sin ω q ( T 21 + Δ t ) J 1 2 θ 3 sin ω q ( T 21 + Δ t ) e t echo / τ decay .
Here, θ 3 is the pulse area of the third SW pulse and the time relative to the echo time is Δ t = t 2 T 21 T 32 . This signal exhibits recoil modulation as a function of T 21 but not as a function of T 32 .
In the presence of gravity, the phase of the three-pulse stimulated echo signal can be shown to be
ϕ g ( 3 ) = q · g ( T 21 2 + T 21 T 32 + T 32 Δ t + 2 T 21 Δ t + Δ t 2 / 2 ) .
As in the two-pulse case, this phase is proportional to the space-time area enclosed by the interferometer. Setting T 32 = 0 reduces ϕ g ( 3 ) to the earlier result for ϕ g ( 2 ) .
To explain the characteristics of the echo signal, we once again decouple the prediction into a part that is dependent on the time scales T 21 , T 32 , and a second part that is dependent on Δ t . Therefore, the three-pulse echo signal is written as
E g ( 3 ) = E D ( 3 ) ( Δ t ) E AI ( 3 ) ( T 21 , T 32 ) e i ϕ D ( 3 ) ( Δ t ) e i ϕ AI ( 3 ) ( T 21 , T 32 ) ,
where
ϕ D ( 3 ) ( Δ t ) = q · g T 32 + 2 T 21 Δ t + Δ t 2 / 2
is the Doppler phase, and
ϕ AI ( 3 ) ( T 21 , T 32 ) = q · g ( T 21 2 + T 21 T 32 )
is the AI phase. The Doppler component of the electric field amplitude is given by
E D ( 3 ) ( Δ t ) E RO e ( Δ t / τ coh ) 2 J 1 2 θ 1 sin ( ω q Δ t ) ,
and the AI electric field amplitude is given by
E AI ( 3 ) ( T 21 , T 32 ) = J 1 2 θ 2 sin ω q ( T 21 + Δ t ) J 1 2 θ 3 sin ω q ( T 21 + Δ t ) e t echo / τ decay .
We note that ϕ D can be varied by changing either T 32 or T 21 . As in the two pulse AI, this term produces a modulation of the echo envelope due to gravitational acceleration. Since t echo = 2 T 21 for the two-pulse AI and t echo = 2 T 21 + T 32 for the three-pulse AI, the functional forms of ϕ D ( 2 ) and ϕ D ( 3 ) are in fact identical.
The functional forms of the electric field amplitudes in Equation (22) and Equation (31) are similar. However, since the three-pulse stimulated AI amplitude involves the additional experimental parameter T 32 , the envelope of this echo signal can be recorded as a function of T 32 for an optimized value of T 21 . For non-zero T 32 , ϕ AI ( 3 ) is maximized if T 32 = 2 T 21 .
The dashed lines in Figure 8b show the Doppler electric field amplitude predicted by Equation (36). This shape is plotted by choosing a convenient T 21 to maximize the recoil-modulated signal, modelled by Equation (37). The solid red line shows oscillations within the echo envelope due to gravity, as predicted by E D ( 3 ) e i ϕ D ( 3 ) .
The solid red line in Figure 8d shows the shape of the in-phase component of the signal amplitude for the three-pulse stimulated echo AI as a function of T 32 , as predicted by E AI ( 3 ) e i ϕ AI ( 3 ) (see Equations (35) and (37)). This signal exhibits a characteristic period τ v determined by T 21 and shows no recoil modulation. The total signal amplitude E ( 3 ) ( T 21 , T 32 ) predicted by Equation (37) as a function of T 32 is shown in Figure 8d as a grey line. This curve exhibits a smooth decay due to signal loss arising from transit time and decoherence effects.

3.5. Results

3.5.1. Doppler Phase Measurements

The characteristics of the echo envelope can be used to measure g along the axis of SW excitation. Equation (25) for the two-pulse AI and Equation (34) for the three-pulse stimulated echo AI show that the Doppler phase ϕ D produces a similar modulation of the echo envelope for the two configurations.
The echo envelope has a simple dispersion shape shown in Figure 8a, and predicted by Equation (27) if T 21 and T 32 are small. As T 21 and T 32 are increased, the signal envelope exhibits oscillations due to gravity as shown by the single sequence acquisitions in Figure 9a,b for the two-pulse and three-pulse stimulated echo configurations, respectively. The oscillations due to g are evident for the echo time 2 T 21 = 9 . 3 ms in Figure 9a. In Figure 9b, the echo time 2 T 21 + T 32 = 45 . 1 ms with T 21 = 1 . 5 ms. The increase in modulation frequency within the echo envelope as a function of T 21 for the two-pulse AI (see Equations (25) and (27)) and as function of T 32 for the three-pulse AI (see Equations (34) and (36)) were used in Refs. [58,87] to measure g with a precision ranging from 0.6% to 0.8%. For these measurements, the Doppler modulation frequency across the echo envelope was assumed to be a constant since T 21 , T 32 Δ t and the frequency was determined from eight repetitions. Although the relatively short measurement time scale is appealing, the sensitivity to the fit functions used to model the echo envelope and the temporal duration of the signal (a few microseconds) limited the statistical uncertainty. The utility of this technique can be re-examined in an actively vibration-stabilized apparatus that is discussed later in this section.

3.5.2. Two-Pulse AI Measurement

The best method of obtaining amplitudes of the in-phase and in-quadrature components of the signal requires fitting to the signal shape as in Figure 9, but the consistency of the results is affected by the complicated fit functions. As a result, the signal from the two pulse AI is usually background subtracted, and the points are squared and summed over the signal duration to extract the two components of the signal amplitude as a function of either T 21 . The quadrature sum of the component amplitudes gives the total signal amplitude E 0 ( 2 ) , and each of the component amplitudes is normalized with respect to E 0 ( 2 ) to obtain cos ( ϕ AI ( 2 ) ) and sin ( ϕ AI ( 2 ) ) . Although this procedure is suitable for extracting ϕ AI ( 2 ) predicted by Equation (26), it is particularly sensitive to background subtraction and the signal strength, and ignores the frequency variation across the echo envelope predicted by Equation (23). For each value of T 21 , the signal amplitude E 0 ( 2 ) is calculated by averaging over three successive points on either side. This procedure ensures that the sinusoidal fits are not skewed by the scatter in the signal strength.
Figure 10 shows a measurement of g using the quadratic dependence of ϕ AI ( 2 ) on T 21 as predicted by Equation (26). The amplitude of the in-phase component is recorded as a function of T 21 for one hour with four observational windows, with each window consisting of 200–325 points acquired in a randomized sequence. Each data point represents the average of 16 repetitions. The error bars are determined on the basis of a probability density function (PDF) analysis [120]. The weights of the error bars are assigned as the product of the error bars from the PDF analysis and the error bars based on the signal amplitude E 0 ( 2 ) . The overall time scale was limited to T 21 = 12 . 8 ms due to the progressive breakdown of the periodically initialized RF phase. These data show the expected chirped frequency dependence of cos ( ϕ AI ( 2 ) ) on T 21 . The data are fit to a multi-parameter function of the form cos ( q g T 21 2 + q v 0 T 21 + ϕ 0 ) , which yielded a measurement of g = 9 . 79123 ( 9 ) m/s 2 . Similarly, we obtain g = 9 . 79130 ( 9 ) m/s 2 from the in-quadrature component with a similar uncertainty. A weighted average of the in-phase and in-quadrature components allowed the acceleration to be determined with a statistical precision of 7 ppm. In the fit function, v 0 models a velocity parameter for the atoms, and ϕ 0 is the initial phase of the grating with respect to the nodal point on the corner-cube reflector. The typical value of v 0 from the fit was 0.107(1) mm/s, which was much smaller than results obtained by tracking the centroid of the falling cloud with a CCD camera. Since cloud launch does not affect the phase of the two-pulse AI, we speculate that intensity imbalances in the two SW components produce this effect.
We note that the size of the residuals in all four windows of Figure 10 is smaller than the standard deviation for a random distribution of points [58]. Here, the standard deviation of the histogram of the residuals in phase units for the entire data sets is 0.7 rad out of a total phase accumulation of | ϕ AI ( 2 ) | 2 . 6 × 10 4 rad for T 21 = 12 . 85 ms. The increasing size of the residuals for T 21 > 10 ms in Ref. [58] indicates the sensitivity of the two-pulse AI to vibrations and decoherence effects such as magnetic field curvature.

3.5.3. Three-Pulse Measurement

We now discuss the improvements obtained using the three-pulse stimulated echo AI. Due to the relative insensitivity of the three-pulse AI to vibrations compared to the two-pulse AI, two previously-mentioned analysis techniques that have distinct disadvantages, namely fitting to the echo envelope, as well as the faster square-sum method can be avoided. Instead, the instantaneous amplitude of the background-subtracted signal is found from a single time slice of the echo envelope, as shown in Figure 8a. The best statistical precision was obtained with a temporal slice duration of 10 ns in which there is effectively no change in the signal amplitude. The average amplitude of each slice was determined by averaging 16 repetitions. This method ensures that the signal amplitude can be measured for 200 slices over the the echo envelope for each value of T 32 .
Figure 11a,b show cos ( ϕ g ( 3 ) ) for a single slice as a function of T 32 with T 21 fixed at 7 . 431900 ms. The in-phase and in-quadrature components cos ( ϕ g ( 3 ) ) and sin ( ϕ g ( 3 ) ) were obtained by normalizing with respect to the total signal amplitude E 0 ( 2 ) as for the two-pulse AI. These data were recorded in one hour with 100 points in each window acquired in randomized sequence. The data shows that the signal exhibits a single frequency as predicted by Equation (35). Two widely spaced observational windows allow this frequency to be determined precisely. The frequency for a single slice can be written as
ω AI ( 3 ) = ϕ AI ( 3 ) T 32 = q · g T 21 .
Here, the frequency change across the echo envelope predicted by Equation (32) is ignored. We use Equation (38) to fit data for the in phase component with T 21 = 7 . 431900 ms and q = 16105651 . 65 rad/m [103], to obtain g = 9 . 833245 ( 4 ) m/s 2 , which represents a statistical precision of 0.4 ppm. This value can be compared to the 7 ppm statistical uncertainty for the two-pulse AI.
The enhancement in precision can be attributed to several factors. Firstly, the determination of a single frequency in the absence of recoil modulation increases the robustness of fits. Secondly, the measurement time scale is significantly extended in comparison to the two-pulse AI since these data represent a total time scale 2 T 21 + T 32 = 45 ms, while limiting T 21 to ∼7.5 ms. Therefore, there is reduced sensitivity to the effects of magnetic curvature and vibrations, which also leads to a more gradual decay of the signal amplitude. As a result, the standard deviation of the residuals (∼0.11) is similar in each observational window. In phase units, the standard deviation of the residuals for the entire data set is 0.2 rad. In comparison, the overall standard deviation is 0.7 rad for the two-pulse AI. The results show that the slicing technique works well only because of good phase stability.
The slicing technique also allows the frequency across the echo envelope predicted Equation (32) to be observed. Additionally, a further improvement in statistical uncertainty is achieved by determining the frequencies of all time slices across the echo envelope. Figure 12a shows the frequency of each time slice obtained using two widely spaced observational windows as in Figure 11 as a function of Δ t . Here, the echo envelope is divided into 200 slices, each with a duration of 10 ns. The data confirm the linear dependence of the angular frequency on Δ t predicted by
ω a ( 3 ) = ϕ a ( 3 ) T 32 = q · g ( T 21 + Δ t ) ,
with ϕ a ( 3 ) given by Equation (32). Each data point in Figure 12a has a typical error bar of 1 ppm. The reduction in error in comparison to the two-pulse AI is a result of the extended (30 ms) time scale. The improved sensitivity makes it possible to observe the change in frequency with Δ t across the echo envelope. In contrast, for the data associated with measurements of g using the Doppler phase modulation inside the echo envelope, we assume a constant frequency across the echo envelope since the observational window is only a few microseconds long and the measurement does not have the desired sensitivity.
For the data in Figure 12a, the limited time scale of the echo envelope and the scatter lead to an overall statistical error in the slope that is appreciable (600 ppm) despite the relatively small statistical error in each of the points (1 ppm). The scatter is attributed to magnetic field effects described in Refs. [58,87]. However, the uncertainty in the frequency intercept is much more tightly constrained since the data are closely clustered near Δ t = 0 . Based on Equation (39), the frequency intercept is q · g T 21 . From the linear fit in Figure 12a, we find the intercept at Δ t = 0 to be 188192.095(17) Hz. Using the values of q and T 21 , we obtain g = 9 . 8774458 ( 9 ) m/s 2 , which represents a statistical precision of 90 ppb. A weighted average from the in-phase and in-quadrature components gives a combined statistical precision of 75 ppb. Figure 12b shows the residuals to the straight line fit in Figure 12a. The residuals increase in size in the regions where the echo envelope is small such as in the extremities and in the vicinity of Δ t = 0 .

3.6. Improvements and Future Work

To summarize the gravity measurements described here, we note that an analysis of the Doppler phase oscillations of the echo envelope resulted in measurements statistically precise to 0.6%. Experiments with the amplitude and phase of the two-pulse AI yielded a statistical uncertainty of 7 ppm. Experiments with the three-pulse AI with a drop height of 1.2 cm in a passively vibration-stabilized chamber demonstrated the best statistical precision of 75 ppb, but it is currently not competitive with the precision of interferometers based on two-photon Raman or Bragg transitions [22,36,46,47,48,52,114].
The time scale of the experiments described here was principally limited by the magnetized vacuum chamber. The magnetization of the chamber and a correction due to the index of refraction were the dominant sources of systematic errors [58]. Other constraints were associated with the signal-to-noise ratio due to the power available from the laser source. The index correction, which impacts q, is dependent on both the sample density and the detuning of the excitation [104], and it is prominent because of the near-resonant nature of the experiment. Based on the improved signal-to-noise ratio obtained in recent echo experiments that utilized a non-magnetic apparatus [32,57], we have projected a measurement sensitivity of 0.6 ppb for the three-pulse AI and 0.3 ppb for the two-pulse AI. Such an experiment can be realized with drop heights of 30 cm in a non-magnetic apparatus with active stabilization of the inertial reference frame. Such an experiment will also require excitation beams with a power output of several Watts so that far-off resonant excitation is feasible. This seemingly challenging requirement can be addressed by using low-cost laser systems consisting of tapered amplifier waveguides seeded by auto-locked diode laser systems [61] that are discussed in Section 5. It is expected that the new laser source will produce a ten-fold increase in signal size due to excitation of higher-order momentum states. Active vibration stabilization will allow cold atoms to be preloaded into a one-dimensional optical lattice so that the initial spatial distribution of atoms has a significant λ / 2 -periodic component [106]. This approach will produce gratings with reflectivities approaching unity [107], which would represent an appreciable increase in the signal-to-noise ratio compared to current experiments in which the reflectivity is 0 . 001 . This increase will permit the experiment to be carried out at a lower density to further reduce the refractive index correction. Atoms will also be selected in the magnetically “insensitive” m F = 0 state to avoid systematic effects due to magnetic interactions. A reduction in measurement time will be realized by using techniques for under-sampling fringe patterns [133] developed for commercial corner-cube gravimeters. We now discuss the layout of the improved apparatus.
The main limitation in the passively-stabilized experiment discussed here is the lack of reference to a proper inertial frame. Since the interferometer’s SW excitation beam is generated by retro-reflecting a traveling-wave beam off of a corner-cube retro-reflector, the phase of the SW excitation was linked to the position of the corner-cube with respect to Earth’s gravity field. This means that during the course of the experiment, there would be an uncontrolled phase accumulation in the signal proportional to the drift of the position of the apparatus. This puts an upper limit on the possible time scale of experiments, which in turn limits the achievable precision. Depending on the vibration spectrum of the lab environment, there might also be aliasing effects for certain experimental repetition rates, as suggested by the data in Refs. [58,87].
The new apparatus shown in Figure 13 will include a high-precision accelerometer, attached to the corner-cube reflector to provide a more reliable reference of the experiment’s frame of reference to the Earth’s frame. The position of the reflector as a whole can be measured at the beginning of each experiment, and this data can be used to either post-correct the data, or to actively correct the frame of reference using a mechanical actuator such as a voice-coil [22,24]. The sensitivity of the interferometer to motion of the SW potential during the experiment can be easily calculated [40,134], making the post-correction process relatively straightforward.
The apparatus will be passively-isolated from ambient vibrations by a two-stage system. The first stage consists of a laminar flow pneumatic isolator, on which the optics table is mounted. This isolator has a vertical resonance around 1 Hz. To further decrease the system’s susceptibility to low-frequency vibrations, the atom trapping optics, AI optics, and vacuum chamber are mounted on a passively-tuned spring damping system that has a sub-1 Hz vertical resonance. This isolator is, in turn, mounted on top of the previously-mentioned pneumatic isolation system. The only component of the apparatus that is not mounted to the spring-based isolator is the field cancellation coils, due to the mass limit of the spring-based isolator. The coils are instead attached directly to the pneumatically-isolated table. Since the zero field/gradient region created by the cancelling coils is quite large compared to the expected motion of the apparatus due to vibrations, the mounting scheme for the coils is not expected to affect the stability of the measurements.
Because the AI relies on two traveling wave components to produce SW excitation, it is necessary to monitor the drift of the retro-reflecting mirror with respect to the other optics in the system using a Michelson-type interferometer and a far-off resonance probe beam. One arm of this interferometer will include the retro-reflecting mirror and the polarizing-cube beam splitter along the path of the vertical excitation beam path and the output of the interferometer will be recorded on a sensitive photodetector. This signal from the probe will be used to postcorrect the AI phase measurement. Such a correction will account for the movement of the lower optics with respect to the retro-mirror. Since the retro-mirror’s position with respect to Earth will be monitored by a sensitive accelerometer, the phase of the SW excitation with respect to Earth will be fully determined. In an alternative scheme, the signal from the interferometer could be used to actively control the position of the retro-mirror with respect to the lower optics through the aforementioned voice coil and a PID control loop [8,22].
Another important shortcoming of the apparatus used in this work related to the magnetic field canceling coils. In the absence of state selection, the presence of magnetic field gradients and magnetic field curvature in the vicinity of the atoms during the experiment results in an additional force, which cannot be separated from the gravitational force. As such, the available time scale for experiments was also limited by the incomplete magnetic field cancellation and the magnetization of the vacuum chamber. The cancellation was incomplete due to the non-optimal geometry of the field-canceling coils and the magnetization of the vacuum chamber. In brief, the B-field at the position of the atoms is canceled by 3 sets of mutually-perpendicular coils. Optimal cancellation is obtained when each set of coils is arranged in a “Helmholtz” configuration, where the (circular) coils are separated by a distance equal to their radius. In our experiment, the coils were set up in a pseudo-Helmholtz configuration, with square coils separated by their side length. This resulted in a well-controlled field only in the immediate vicinity of the atom trap, with no extended volume of zero field. The new apparatus will make use of three sets of square cross-section Helmholtz coils, whose ideal distance of separation was found using numerical modeling to be 0.55 times the side length. These coils produce an extended volume of zero magnetic field, allowing the atoms to experience free-fall for several hundred milliseconds. If the atoms are launched in a fountain, the experimental time scale can be further extended to 1 s. Extension of time scales beyond 50 ms used in this work is also expected to require chirped standing wave excitation to avoid loss of signal amplitude due to the differential Doppler shift of falling atoms with respect to the traveling wave components of the standing wave. This limitation can be overcome using chirped, counter-propagating traveling-wave beams to cancel the Doppler shift due to falling atoms [32,57].

4. Magnetic Coherences

Properties of atomic coherences have been exploited for interesting applications such as quantum state preparation [135] and control [136], nondemolition measurements [137,138], entanglement [139], and precision magnetometry [140,141]. Other interesting applications have related to studies of velocity-changing collisions using both radio frequency (RF) [142] and optical excitation [143] and experiments with coherent transient effects that involve coupling between Zeeman sublevels [144,145]. In this section, we review applications of coherent transient techniques developed for atom interferometry for precise measurements of the strength of magnetic interactions such as atomic g-factor ratios [59,60].
For the experiments considered in Ref. [60], the sample of atoms is excited using two simultaneous traveling-wave laser pulses applied at t = 0 with wave vectors k 1 and k 2 at a small angle θ 10 mrad, as shown in Figure 14a. The atomic sample is a cloud of laser-cooled atoms loaded into a MOT. The individual traveling waves pulses have orthogonal linear or circular polarizations and are detuned from the excited state. However, the pulses are resonant with the two-photon transition that couples two magnetic sublevels of the ground state, as shown in Figure 14b. The timing diagram is shown in Figure 14c. The excitation creates a spatially-periodic superposition between the magnetic sublevels of the ground state coupled by the laser fields. The superposition has a period of λ / θ , where λ is the wavelength of the excitation. This coherence grating is probed by a read-out pulse along k 2 and the resulting signal, referred to as magnetic-grating free induction decay (MGFID) [146], is coherently scattered along k 1 due to conservation of momentum. The grating dephases due to the thermal motion of the atoms, and the time scale of the signal decay is determined by the time taken by a typical atom to move a distance on the order of a grating spacing ( λ / θ u ), where u is the most probable speed associated with the Maxwell-Boltzmann velocity distribution. A magnetic grating echo (MGE) is observed using a second set of excitation pulses at t = T to rephase the coherence grating, as in Figure 14c [146,147,148]. The second pulse modifies the time-dependent coefficients that describe the coherent superposition of magnetic sublevels so that the grating reforms at t = 2 T . This is analogous to the reversal of the Doppler phases of individual atoms in a traditional two-pulse photon echo experiment [149] that can be used for measuring the excited state lifetime [150]. In the absence of decoherence due to collisions and background light, the MGE amplitude should decay on a time scale determined by the transit time of atoms through the laser beams [148].
The focus of initial experimental work using the MGE was related to observing effects due to atomic recoil [147]. The MGFID and MGE were used to verify the expected dependence of the dephasing time of the coherences on the velocity distribution of the sample [147,148], observing the effect of magnetic fields for particular experimental configurations [147], and studying the effects of collisions [148]. Other applications of the MGFID include measurements of the diffusion constants [151] and phase space imaging [152]. Echo techniques were also used to investigate applications related to detecting nanostructures and depositing periodic arrays of atoms on substrates [153].
In a static magnetic field, the Zeeman shift between magnetic sublevels causes temporal oscillations within the envelopes of the MGFID and MGE signals at multiples of the Larmor frequency ω L . In Ref. [60], an analytical calculation was developed to predict the functional form of the Larmor oscillations in the MGFID for arbitrarily-directed weak static magnetic fields in situations in which the excitation pulses consisted of traveling waves with orthogonal linear and circular polarizations. This theoretical treatment is based on a rotation matrix approach [154,155,156], in which the effect of the magnetic field can be described as a time-dependent rotation of the atomic system about the quantization axis. This treatment accurately described signals observed from room temperature vapor and laser-cooled atoms. In the absence of magnetic fields, the velocity distribution of a cold sample measured using the MGFID agreed with the results of an independent technique used to measure the sample temperature [131]. Additionally, a rate equation treatment [157,158] was used to understand the properties of the MGE in a magnetic field. Thus, Ref. [60] renewed interest in precise measurements of atomic g-factor ratios using the MGFID and MGE signals and compact laser-cooled samples that could be excited in uniform magnetic fields.
The best measurements of Zeeman shifts [159] and atomic g-factor ratios [160] were carried out with uniform ∼50 G B fields in centimeter-length, paraffin-coated atomic vapor cells to avoid decoherence due to wall collisions. These experiments relied on RF spectroscopy and lamp-based optical pumping to determine the centers of transitions between atomic ground states of Rb isotopes. The typical accuracies of a few ppm were sufficient for independent tests of the non-relativistic Zeeman Hamiltonian. However, a measurement of g-factor ratios precise to 500 ppb or better can be used to test corrections arising from the self-energy of the electron and vacuum polarization. These effects have been incorporated into atomic theory describing the electron and nuclear g-factors g J and g I , respectively [161], but have not been verified by experiments. The predicted corrections to g J and g I depend on the nuclear mass and spin, which of course differ among isotopes. Measurements of isotopic atomic g-factor ratios therefore constitute a sensitive test of specific aspects of QED that relate to magnetic interactions. The importance of such experiments for matter-antimatter comparisons has been highlighted by the recent measurements of the magnetic moment of the antiproton [162].
Although measurements of atomic g-factor ratios were undertaken using the MGFID [54], the statistical precision was overwhelmed by systematic effects due to AC Stark shifts. Therefore, an alternative technique, which is the basis of atomic magnetometers [163] was adopted in Ref. [59] to demonstrate the most sensitive measurement of g-factor ratios. Our statistical precision of 690 ppb, obtained on a 10-ms time scale, exceeded the 2 ppm sensitivity of a long-standing measurement [160]. In our experiment, Larmor oscillations of ground state coherences from a dual isotope MOT containing 85 Rb and 87 Rb atoms were simultaneously recorded. The measurement exploited the compact dimensions of the sample (few millimeters) over which a uniform, actively-stabilized magnetic field was applied. Coherences involving F = 3 ( F = 2 ) hyperfine ground states in 85 Rb ( 87 Rb) were optically excited with a circularly-polarized pulse with a wave vector perpendicular to a constant, weak ( B < 1 G) magnetic field. Coherences between adjacent magnetic sublevels in each isotope can be modeled as a magnetization that precesses in the B field at the Larmor frequency ω L = g F μ B B / . The precession can be measured with a linearly-polarized probe pulse directed orthogonal to the B field. A rotation in the polarization of the probe pulse due to the differential absorption of its σ + and σ components can be recorded on a balanced detector. This signal exhibits sinusoidal oscillations at ω L . Therefore, simultaneous measurements of ω L 85 and ω L 87 as a beat note (see Figure 15) can be used to measure the ratio of Lande g F -factors.
The precision of 690 ppb and data acquisition time (3 minutes) allowed systematic effects due to Breit-Rabi corrections and AC Stark shifts to be investigated [59]. Results showed that the nonlinear variation of ω L on B due to the Breit-Rabi effect did not cancel when the ratio of g F -factors was measured, resulting in a systematic shift that depended on B. The shift also depended on the magnetic quantum number m F . Therefore, parameters of the probe pulse and B fields affected the weighted sum (of pairs) of coherences across the ground state manifolds that contribute to the signal. These systematic effects were also qualitatively confirmed using numerical simulations.
Subsequently, atom interferometric experiments that shared the same setup [32,57] achieved even better control of B-fields and B-gradients resulting in observation times of ∼250 ms—limited only by the transit time of cold atoms. Such an experiment with a 100-ms time scale will be ideally suited for this measurement. The experiment will resolve Larmor oscillations associated with all pairs of coherences in the ground state manifolds of the two isotopes by using a larger, uniform, pulsed B field of about 5 G. Additionally, Larmor oscillations from coherences in the F = 2 ( F = 1 ) ground states of 85 Rb ( 87 Rb) will be measured in addition to signals from the F = 3 ( F = 2 ) ground states in 85 Rb ( 87 Rb). This approach will eliminate Breit-Rabi corrections and allow the measurement of g J and g I . The estimated measurement precision of ∼100 ppb will permit tests of predicted relativistic corrections at the desired level [161]. The possibility of realizing similar measurements with RF excitation and optical pumping techniques can also be investigated.

5. Auto-Locked Lasers

We have presented a review of recent results pertaining to measurements of g, atomic recoil frequency, and atomic g-factor ratios. Although all these experiments have recorded significant improvements in precision, the measurements are dominated by systematic effects such as the index of refraction due to the near-resonant excitation using laser beams with power outputs of ∼100 mW. We envision control of systematic effects at the ppb level by using far-detuned excitation beams with relatively high-power outputs of a few Watts. Ideally, several high-power lasers will be required for atom trapping, atom interferometry, and lattice loading—a crucial step for increasing the contrast (and signal-to-noise ratio) of density gratings in future echo interferometer experiments. Commercial laser sources such as Ti:Sapphire lasers with output powers of ∼1 W are inadequate apart from being expensive to maintain over the data collection times of several thousand hours. A possible option is to use frequency-doubled light from erbium-doped fiber-amplified laser sources operating in the 1560 nm telecommunications band [164,165], with output powers of >10 W. However, the cost is comparable to Ti:Sapphire lasers and the lifetime of these sources has not been tested beyond a few years.
To fulfill experimental requirements, we have relied on elegant and practical designs for external cavity diode lasers (ECDLs) and controllers [166,167,168,169,170,171] that produce power outputs of 1-100 mW and linewidths of ∼1 MHz, to develop unique, low-cost, vacuum-sealed, auto-locked external cavity diode laser systems (ALDLS) that can be integrated using components from original equipment manufacturers (OEM), specially machined parts, and powerful central processors [61]. The laser source depends on optical feedback from a narrow-band interference filter to realize a narrow laser linewidth and improved spatial mode quality. The thermally-stabilized laser cavity can be evacuated within minutes and vacuum-sealed for several months, so that laser system exhibits reduced sensitivity to environmental temperature and pressure fluctuations. These lasers can be locked or scanned with respect to atomic or molecular spectral lines without the need for human intervention using a digital controller. The laser cavity relies on an interchangeable optics kit consisting of a laser diode and optical feedback elements to operate in the desired wavelength range. The digital signal processor of the controller is capable of storing a variety of algorithms in its memory for laser frequency stabilization using techniques such as pattern matching and first or third derivative feedback. Additional features include power amplification to several Watts using semiconductor waveguides [172], frequency locking to external cavities, and rapid amplitude modulation for wide-ranging applications. These capabilities are not commercially available as a package despite the widespread availability of the constituent OEM components.
In recent trials, ALDLS systems, operating at both 780 nm and 633 nm, have proven their linewidth (∼300 kHz) and lock stability (∼500 kHz) through accurate measurements of gravitational acceleration g using a state-of-the-art industrial gravimeter (Scintrex FG5X) with an absolute accuracy of 1 ppb. The value of g determined using our prototype lasers was in agreement with the baseline established using an iodine-stabilized He-Ne laser (which is an industrial standard), and exhibited lower scatter. These results suggest that systematic effects in cold-atom measurements of g [58] can be characterized using an industrial gravimeter and a common ALDLS-based light source operating at 780 nm. Three improvements that were targeted have also been realized in recent tests. Firstly, the vacuum-sealed cavity reduced frequency drifts by a factor of fifty, as shown in Figure 16, so that the lock stability could be extended to 24-hour time scales. Figure 17 shows the Allan deviation of a typical OEM current controller used to operate the ALDLS. Our results suggest that significant reduction in current noise is possible, since the best current noise density specification achieved in commercial controllers (which use a Libbrecht-Hall design [171]) is ∼200 pA/ Hz . Secondly, interchangeable optics kits used in the same laser head have allowed operation at two widely-separated wavelengths accessible to laser diodes. Thirdly, the controller demonstrated the use of different algorithms for locking to atomic and molecular resonances. In the 780 nm band, the control system relied on atomic rubidium spectroscopy. To lock the laser, the controller implemented pattern matching between Doppler-free spectra obtained in real time by scanning the laser and the reference peaks stored in the processor’s memory by using a sliding correlation algorithm. This technique determines the frequency offset between the scanned pattern and the reference pattern. The offset is fed back as a control voltage to the piezo transducer that controls the length of the laser cavity. As a result, the laser frequency is iteratively brought closer to the desired frequency with variable scan amplitude. In the 633 nm band, the controller achieved tighter locks by generating the third derivatives of narrower molecular resonances in iodine. Additionally, the controller was automatically able to relock even in the event of a mode hop, because it could be programmed to vary the diode current and temperature and compare spectral peaks over a scan range of up to 10 GHz. The availability of tapered amplifier waveguides with power outputs of several Watts, and techniques to transiently increase the drive current [173] to achieve a several-fold enhancement in output power, suggest that suitable low-cost laser systems can be developed for pursuing echo experiments in a far-off resonance configuration.

6. Conclusion

We have presented an overview of the best-suited echo interferometric techniques to realize precision measurements of g and ω q , and we have reviewed coherent transient methods for obtaining the most sensitive measurement of atomic g-factor ratios. The proposed apparatus outlined at the end of Section 3 will combine features of long-lived echo experiments in a non-magnetic chamber in which atoms can be routinely cooled to 1 μK using polarization-gradient cooling with the advantages of vibration stabilization. As a result, it should be possible to load the laser-cooled sample into a deep optical lattice to achieve a several hundred-fold increase in grating contrast and signal size. Echo AIs operating in this setup in an off-resonant configuration using low-cost, high-power lasers appear to be capable of reducing both the statistical and systematic uncertainties to competitive levels. As a result, the relative simplicity of echo AI experiments may be well suited for several precision measurements.

Acknowledgments

This work was supported by the Canada Foundation for Innovation, Ontario Innovation Trust, the Natural Sciences and Engineering Research Council of Canada, Ontario Centres of Excellence, the US Army Research Office, and York University. We acknowledge the contributions of Carson Mok, who carried out the gravity experiments, and of Iain Chan, who carried out the magnetic coherence experiments. We thank Tycho Sleator of New York University for the generous loan of Ti:Sapphire laser components.

Author Contributions

Brynle Barrett wrote the introductory Section 1.1 and Section 1.2, carried out the recoil frequency measurements presented in Section 2, and contributed to the gravity and magnetic coherence experiments presented in Section 3 and Section 4, respectively. Adam Carew designed and constructed a modified apparatus for gravity measurements, and he contributed to the measurements of gravitational acceleration, atomic recoil, magnetic coherences, and the design of the auto-locked lasers presented in Section 5. Hermina C. Beica prepared the manuscript, and contributed to the design and characterization of the auto-locked lasers. Andrejs Vorozcovs led the design and construction of the auto-lock controller. Alexander Pouliot constructed and characterized diode laser controllers and carried out an analysis of gravity measurements from the two-pulse AI.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Moskowitz, P.E.; Gould, P.L.; Atlas, S.R.; Pritchard, D.E. Diffraction of an Atomic Beam by Standing-Wave Radiation. Phys. Rev. Lett. 1983, 51, 370–373. [Google Scholar] [CrossRef]
  2. Gould, P.L.; Ruff, G.A.; Pritchard, D.E. Diffraction of Atoms by Light: The Near-Resonant Kapitza-Dirac Effect. Phys. Rev. Lett. 1986, 56, 827–830. [Google Scholar] [CrossRef] [PubMed]
  3. Keith, D.W.; Ekstrom, C.R.; Turchette, Q.R.; Pritchard, D.E. An Interferometer for Atoms. Phys. Rev. Lett. 1991, 66, 2693–2696. [Google Scholar] [CrossRef] [PubMed]
  4. Carnal, O.; Mlynek, J. Young’s Double-Slit Experiment with Atoms: A Simple Atom Interferometer. Phys. Rev. Lett. 1991, 66, 2689–2692. [Google Scholar] [CrossRef] [PubMed]
  5. Schmiedmayer, J.; Chapman, M.S.; Ekstrom, C.R.; Hammond, T.D.; Wehinger, S.; Pritchard, D.E. Index of Refraction of Various Gases for Sodium Matter Waves. Phys. Rev. Lett. 1995, 74, 1043–1047. [Google Scholar] [CrossRef] [PubMed]
  6. Weiss, D.S.; Young, B.C.; Chu, S. Precision Measurement of the Photon Recoil of an Atom using Atomic Interferometry. Phys. Rev. Lett. 1993, 70, 2706–2709. [Google Scholar] [CrossRef] [PubMed]
  7. Riehle, F.; Kisters, T.; Witte, A.; Helmcke, J.; Bordé, C.J. Optical Ramsey spectroscopy in a rotating frame: Sagnac effect in a matter-wave interferometer. Phys. Rev. Lett. 1991, 67, 177–180. [Google Scholar] [CrossRef] [PubMed]
  8. Kasevich, M.; Chu, S. Atomic interferometry using stimulated Raman transitions. Phys. Rev. Lett. 1991, 67, 181–184. [Google Scholar] [CrossRef] [PubMed]
  9. Berman, P.R. (Ed.) Atom Interferometry; Academic Press: San Diego, CA, USA, 1997.
  10. Cronin, A.; Schmiedmayer, J.; Pritchard, D.E. Optics and interferometry with atoms and molecules. Rev. Mod. Phys. 2009, 81, 1051–1129. [Google Scholar] [CrossRef]
  11. Barrett, B.; Gominet, P.A.; Cantin, E.; Antoni-Micollier, L.; Bertoldi, A.; Battelier, B.; Bouyer, P.; Lautier, J.; Landragin, A. Mobile and remote inertial sensing with atom interferometers. In Proceedings of the International School of Physicc “Enrico Fermi”, Vol. 188 “Atom Interferometry”; Tino, G.M., Kasevich, M.A., Eds.; IOS: Amsterdam, The Netherland, 2014; pp. 493–555. [Google Scholar]
  12. Wicht, A.; Hensley, J.M.; Sarajlic, E.; Chu, S. A Preliminary Measurement of the Fine Structure Constant Based on Atom Interferometry. Phys. Scr. 2002, 2002, 82–88. [Google Scholar] [CrossRef]
  13. Cladé, P.; de Mirandes, E.; Cadoret, M.; Guellati-Khélifa, S.; Schwob, C.; Nez, F.; Julien, L.; Biraben, F. Determination of the Fine Structure Constant Based on Bloch Oscillations of Ultracold Atoms in a Vertical Optical Lattice. Phys. Rev. Lett. 2006, 96, 033001. [Google Scholar] [CrossRef] [PubMed]
  14. Cadoret, M.; de Mirandes, E.; Clade, P.; Guellati-Khélifa, S.; Schwob, C.; Nez, F.; Julien, L.; Biraben, F. Combination of Bloch Oscillations with a Ramsey-Bordé Interferometer: New Determination of the Fine Structure Constant. Phys. Rev. Lett. 2008, 101, 230801. [Google Scholar] [CrossRef] [PubMed]
  15. Bouchendira, R.; Clade, P.; Guellati-Khélifa, S.; Nez, F.; Biraben, F. New determination of the fine structure constant and test of the quantum electrodynamics. Phys. Rev. Lett. 2011, 106, 080801. [Google Scholar] [CrossRef] [PubMed]
  16. Lan, S.Y.; Kuan, P.C.; Estey, B.; English, D.; Brown, J.M.; Hohensee, M.A.; Müller, H. A Clock Directly Linking Time to a Particle’s Mass. Science 2013, 339, 554–557. [Google Scholar] [CrossRef] [PubMed]
  17. Bertoldi, A.; Lamporesi, G.; Cacciapuoti, L.; de Angelis, M.; Fattori, M.; Petelski, T.; Peters, A.; Prevedelli, M.; Stuhler, J.; Tino, G.M. Atom interferometry gravity-gradiometer for the determination of the Newtonian gravitational constant G. Eur. Phys. J. D 2006, 40, 271–279. [Google Scholar] [CrossRef]
  18. Fixler, J.B.; Foster, G.T.; McGuirk, J.M.; Kasevich, M.A. Atom Interferometer Measurement of the Newtonian Constant of Gravity. Science 2007, 315, 74–77. [Google Scholar] [CrossRef] [PubMed]
  19. Lamporesi, G.; Bertoldi, A.; Cacciapuoti, L.; Prevedelli, M.; Tino, G.M. Determination of the Newtonian Gravitational Constant Using Atom Interferometry. Phys. Rev. Lett. 2008, 100, 050801. [Google Scholar] [CrossRef] [PubMed]
  20. Rosi, G.; Sorrentino, F.; Cacciapuoti, L.; Prevedelli, M.; Tino, G.M. Precision measurement of the Newtonian gravitational constant using cold atoms. Nature 2014, 510, 518–521. [Google Scholar] [CrossRef] [PubMed]
  21. Peters, A.; Chung, K.Y.; Young, B.; Hensley, J.; Chu, S. Precision atom interferometry. Phil. Trans. R. Soc. Lond. A 1997, 355, 2223–2233. [Google Scholar] [CrossRef]
  22. Peters, A.; Chung, K.Y.; Chu, S. Measurement of gravitational acceleration by dropping atoms. Nature 1999, 400, 849–852. [Google Scholar]
  23. Snadden, M.J.; McGuirk, J.M.; Bouyer, P.; Haritos, K.G.; Kasevich, M.A. Measurement of the Earth’s Gravity Gradient with an Atom Interferometer-Based Gravity Gradiometer. Phys. Rev. Lett. 1998, 81, 971–974. [Google Scholar] [CrossRef]
  24. McGuirk, J.M.; Foster, G.T.; Fixler, J.B.; Snadden, M.J.; Kasevich, M.A. Sensitive absolute-gravity gradiometry using atom interferometry. Phys. Rev. A 2002, 65, 033608. [Google Scholar] [CrossRef]
  25. Sorrentino, F.; Bertoldi, A.; Bodart, Q.; Cacciapuoti, L.; de Angelis, M.; Lien, Y.H.; Prevedelli, M.; Rosi, G.; Tino, G.M. Simultaneous measurement of gravity acceleration and gravity gradient with an atom interferometer. Appl. Phys. Lett. 2012, 101, 114106. [Google Scholar] [CrossRef]
  26. Lenef, A.; Hammond, T.; Smith, E.; Chapman, M.; Rubenstein, R.; Pritchard, D.E. Rotation Sensing with an Atom Interferometer. Phys. Rev. Lett. 1997, 78, 760–763. [Google Scholar] [CrossRef]
  27. Gustavson, T.L.; Bouyer, P.; Kasevich, M.A. Precision Rotation Measurements with an Atom Interferometer Gyroscope. Phys. Rev. Lett. 1997, 78, 2046–2049. [Google Scholar] [CrossRef]
  28. Barrett, B.; Geiger, R.; Dutta, I.; Meunier, M.; Canuel, B.; Gauguet, A.; Bouyer, P.; Landragin, A. The Sagnac effect: 20 years of development in matter-wave interferometry. C. R. Phys. 2014, 15, 875–883. [Google Scholar] [CrossRef]
  29. Dutta, I.; Savoie, D.; Fang, B.; Venon, B.; Alzar, C.L.G.; Geiger, R.; Landragin, A. Continuous Cold-Atom Inertial Sensor with 1 nrad/s Rotation Stability. Phys. Rev. Lett. 2016, 116, 183003. [Google Scholar] [CrossRef] [PubMed]
  30. Weel, M.; Chan, I.; Beattie, S.; Kumarakrishnan, A.; Gosset, D.; Yavin, I. Effect of a magnetic field gradient and gravitational acceleration on a time-domain grating-echo interferometer. Phys. Rev. A 2006, 73, 063624. [Google Scholar] [CrossRef]
  31. Davis, J.P.; Narducci, F.A. A proposal for a gradient magnetometer atom interferometer. J. Mod. Opt. 2008, 55, 3173–3185. [Google Scholar] [CrossRef]
  32. Barrett, B.; Chan, I.; Kumarakrishnan, A. Atom-interferometric techniques for measuring uniform magnetic field gradients and gravitational acceleration. Phys. Rev. A 2011, 84, 063623. [Google Scholar] [CrossRef]
  33. Hardman, K.S.; Everitt, P.J.; McDonald, G.D.; Manju, P.; Wigley, P.B.; Sooriyabadara, M.A.; Kuhn, C.C.N.; Debs, J.E.; Close, J.D.; Robins, N.P. A quantum sensor: Simultaneous precision gravimetry and magnetic gradiometry with a Bose-Einstein condensate. 2016; arXiv:1603.01967. [Google Scholar]
  34. Bonnin, A.; Zahzam, N.; Bidel, Y.; Bresson, A. Simultaneous dual-species matter-wave accelerometer. Phys. Rev. A 2013, 88, 043615. [Google Scholar] [CrossRef]
  35. Müntinga, H.; Ahlers, H.; Krutzik, M.; Wenzlawski, A.; Arnold, S.; Becker, D.; Bongs, K.; Dittus, H.; Duncker, H.; Gaaloul, N.; et al. Interferometry with Bose-Einstein Condensates in Microgravity. Phys. Rev. Lett. 2013, 110, 093602. [Google Scholar] [CrossRef] [PubMed]
  36. Dickerson, S.M.; Hogan, J.M.; Sugarbaker, A.; Johnson, D.M.S.; Kasevich, M.A. Multiaxis Inertial Sensing with Long-Time Point Source Atom Interferometry. Phys. Rev. Lett. 2013, 111, 083001. [Google Scholar] [CrossRef] [PubMed]
  37. Sugarbaker, A.; Dickerson, S.M.; Hogan, J.M.; Johnson, D.M.S.; Kasevich, M.A. Enhanced Atom Interferometer Readout through the Application of Phase Shear. Phys. Rev. Lett. 2013, 111, 113002. [Google Scholar] [CrossRef] [PubMed]
  38. Aguilera, D.; Ahlers, H.; Battelier, B.; Bawamia, A.; Bertoldi, A.; Bondarescu, R.; Bongs, K.; Bouyer, P.; Braxmaier, C.; Cacciapuoti, L.; et al. STE-QUEST—test of the universality of free fall using cold atom interferometry. Class. Quan. Gravity 2014, 31, 115010. [Google Scholar] [CrossRef]
  39. Schlippert, D.; Hartwig, J.; Albers, H.; Richardson, L.L.; Schubert, C.; Roura, A.; Schleich, W.P.; Ertmer, W.; Rasel, E.M. Quantum Test of the Universality of Free Fall. Phys. Rev. Lett. 2014, 112, 203002. [Google Scholar] [CrossRef]
  40. Barrett, B.; Antoni-Micollier, L.; Chichet, L.; Battelier, B.; Gominet, P.A.; Bertoldi, A.; Bouyer, P.; Landragin, A. Correlative methods for dual-species quantum tests of the weak equivalence principle. New J. Phys. 2015, 17, 085010. [Google Scholar] [CrossRef]
  41. Bonnin, A.; Zahzam, N.; Bidel, Y.; Bresson, A. Characterization of a simultaneous dual-species atom interferometer for a quantum test of the weak equivalence principle. Phys. Rev. A 2015, 92, 023626. [Google Scholar] [CrossRef]
  42. Hartwig, J.; Abend, S.; Schubert, C.; Schlippert, D.; Ahlers, H.; Posso-Trujillo, K.; Gaaloul, N.; Ertmer, W.; Rasel, E.M. Testing the universality of free fall with rubidium and ytterbium in a very large baseline atom interferometer. New J. Phys. 2015, 17, 035011. [Google Scholar] [CrossRef]
  43. Dimopoulos, S.; Graham, P.W.; Hogan, J.M.; Kasevich, M.A. Testing general relativity with atom interferometry. Phys. Rev. Lett. 2007, 98, 111102. [Google Scholar] [CrossRef] [PubMed]
  44. Geiger, R.; Amand, L.; Bertoldi, A.; Canuel, B.; Chaibi, W.; Danquigny, C.; Dutta, I.; Fang, B.; Gaffet, S.; Gillot, J.; et al. Matter-wave laser Interferometric Gravitation Antenna (MIGA): New perspectives for fundamental physics and geosciences. 2015; arXiv:1505.07137. [Google Scholar]
  45. Bordé, C.J. Atom Interferometry With Internal State Labelling. Phys. Lett. A 1989, 140, 10–12. [Google Scholar] [CrossRef]
  46. Hu, Z.K.; Sun, B.L.; Duan, X.C.; Zhou, M.K.; Chen, L.L.; Zhan, S.; Zhang, Q.Z.; Luo, J. Demonstration of an ultrahigh-sensitivity atom-interferometry absolute gravimeter. Phys. Rev. A 2013, 88, 043610. [Google Scholar] [CrossRef]
  47. Gillot, P.; Francis, O.; Landragin, A.; Pereira Dos Santos, F.; Merlet, S. Stability comparison of two absolute gravimeters: optical versus atomic interferometers. Metrologia 2014, 51, L15. [Google Scholar] [CrossRef]
  48. Freier, C.; Hauth, M.; Schkolnik, V.; Leykauf, B.; Schilling, M.; Wziontek, H.; Scherneck, H.G.; Müller, J.; Peters, A. Mobile quantum gravity sensor with unprecedented stability. 2015; arXiv:1512.05660. [Google Scholar]
  49. Hughes, K.J.; Burke, J.H.T.; Sackett, C.A. Suspension of Atoms Using Optical Pulses, and Application to Gravimetry. Phys. Rev. Lett. 2009, 102, 150403. [Google Scholar] [CrossRef] [PubMed]
  50. Poli, N.; Wang, F.Y.; Tarallo, M.G.; Alberti, A.; Prevedelli, M.; Tino, G.M. Precision Measurement of Gravity with Cold Atoms in an Optical Lattice and Comparison with a Classical Gravimeter. Phys. Rev. Lett. 2011, 106, 038501. [Google Scholar] [CrossRef] [PubMed]
  51. Charrière, R.; Cadoret, M.; Zahzam, N.; Bidel, Y.; Bresson, A. Local gravity measurement with the combination of atom interferometry and Bloch oscillations. Phys. Rev. A 2012, 85, 013639. [Google Scholar] [CrossRef]
  52. Altin, P.A.; Johnsson, M.T.; Dennis, G.R.; Anderson, R.P.; Debs, J.E.; Szigeti, S.S.; Hardman, K.S.; Bennetts, S.; McDonald, G.D.; Turner, L.D.; et al. Precision atomic gravimeter based on Bragg diffraction. New J. Phys. 2013, 15, 023009. [Google Scholar] [CrossRef]
  53. Andia, M.; Jannin, R.; Nez, F.; Biraben, F.; Guellati-Khélifa, S.; Cladé, P. Compact atomic gravimeter based on a pulsed and accelerated optical lattice. Phys. Rev. A 2013, 88, 031605. [Google Scholar] [CrossRef]
  54. Barrett, B.; Chan, I.; Mok, C.; Carew, A.; Yavin, I.; Kumarakrishnan, A.; Cahn, S.B.; Sleator, T. Time-Domain Interferometry with Laser-Cooled Atoms. In Advances in Atomic, Molecular and Optical Physics; Arimondo, E., Berman, P.R., Lin, C.C., Eds.; Elsevier Inc.: Amsterdam, The Netherlands, 2011; Volume 60, Chapter 3; pp. 119–199. [Google Scholar]
  55. Cahn, S.B.; Kumarakrishnan, A.; Shim, U.; Sleator, T.; Berman, P.R.; Dubetsky, B. Time-Domain de Broglie Wave Interferometry. Phys. Rev. Lett. 1997, 79, 784–787. [Google Scholar] [CrossRef]
  56. Barrett, B. Techniques for Measuring the Atomic Recoil Frequency Using a Grating-Echo Atom Interferometer. Ph.D. Thesis, York University, Toronto, ON, Canada, 2012. [Google Scholar]
  57. Barrett, B.; Carew, A.; Beattie, S.; Kumarakrishnan, A. Measuring the atomic recoil frequency using a modified grating-echo atom interferometer. Phys. Rev. A 2013, 87, 033626. [Google Scholar] [CrossRef]
  58. Mok, C.; Barrett, B.; Carew, A.; Berthiaume, R.; Beattie, S.; Kumarakrishnan, A. Demonstration of improved sensitivity of echo interferometers to gravitational acceleration. Phys. Rev. A 2013, 88, 023614. [Google Scholar] [CrossRef]
  59. Chan, I.; Barrett, B.; Kumarakrishnan, A. Precise determination of atomic g-factor ratios from a dual isotope magneto-optical trap. Phys. Rev. A 2011, 84, 032509. [Google Scholar] [CrossRef]
  60. Chan, I.; Andreyuk, A.; Beattie, S.; Barrett, B.; Mok, C.; Weel, M.; Kumarakrishnan, A. Properties of magnetic sublevel coherences for precision measurements. Phys. Rev. A 2008, 78, 033418. [Google Scholar] [CrossRef]
  61. Beica, H.C.; Carew, A.; Vorozcovs, A.; Dowling, P.; Pouliot, A.; Singh, G.; Kumarakrishnan, A. An Auto-Locked Diode Laser System for Precision Metrology. 2016; in preparation. [Google Scholar]
  62. Bongs, K.; Launay, R.; Kasevich, M.A. High-order inertial phase shifts for time-domain atom interferometers. Appl. Phys. B 2006, 84, 599–602. [Google Scholar] [CrossRef]
  63. Itano, W.M.; Bergquist, J.C.; Bollinger, J.J.; Gilligan, J.M.; Heinzen, D.J.; Moore, F.L.; Raizen, M.G.; Wineland, D.J. Quantum projection noise: Population fluctuations in two-level systems. Phys. Rev. A 1993, 47, 3554–3570. [Google Scholar] [CrossRef] [PubMed]
  64. Le Gouët, J.; Mehlstäubler, T.E.; Kim, J.; Merlet, S.; Clairon, A.; Landragin, A.; Pereira Dos Santos, F. Limits to the sensitivity of a low noise compact atomic gravimeter. Appl. Phys. B 2008, 92, 133–144. [Google Scholar] [CrossRef] [Green Version]
  65. Rocco, E.; Palmer, R.N.; Valenzuela, T.; Boyer, V.; Freise, A.; Bongs, K. Fluorescence detection at the atom shot noise limit for atom interferometry. New J. Phys. 2014, 16, 093046. [Google Scholar] [CrossRef]
  66. Lan, S.Y.; Kuan, P.C.; Estey, B.; Haslinger, P.; Müller, H. Influence of the Coriolis Force in Atom Interferometry. Phys. Rev. Lett. 2012, 108, 090402. [Google Scholar] [CrossRef] [PubMed]
  67. Chiow, S.w.; Herrmann, S.; Chu, S.; Müller, H. Noise-Immune Conjugate Large-Area Atom Interferometers. Phys. Rev. Lett. 2009, 103, 050402. [Google Scholar] [CrossRef] [PubMed]
  68. Müller, H.; Chiow, S.; Long, Q.; Herrmann, S.; Chu, S. Atom Interferometry with up to 24-Photon-Momentum-Transfer Beam Splitters. Phys. Rev. Lett. 2008, 100, 180405. [Google Scholar] [CrossRef] [PubMed]
  69. Chiow, S.w.; Kovachy, T.; Chien, H.C.; Kasevich, M.A. 102 hk Large Area Atom Interferometers. Phys. Rev. Lett. 2011, 107, 130403. [Google Scholar] [CrossRef] [PubMed]
  70. McDonald, G.D.; Kuhn, C.C.N.; Bennetts, S.; Debs, J.E.; Hardman, K.S.; Johnsson, M.; Close, J.D.; Robins, N.P. 80hk momentum separation with Bloch oscillations in an optically guided atom interferometer. Phys. Rev. A 2013, 88, 053620. [Google Scholar] [CrossRef]
  71. Clade, P.; Guellati-Khélifa, S.; Nez, F.; Biraben, F. Large Momentum Beam Splitter Using Bloch Oscillations. Phys. Rev. Lett. 2009, 102, 240402. [Google Scholar] [CrossRef] [PubMed]
  72. Clauser, J.F.; Li, S. Talbot-vonLau atom interferometer with cold potassium. Phys. Rev. A 1994, 49, R2213–R2216. [Google Scholar] [CrossRef] [PubMed]
  73. Chapman, M.S.; Ekstrom, C.R.; Hammond, T.D.; Schmiedmayer, J.; Tannian, B.E.; Wehinger, S.; Pritchard, D.E. Near-field imaging of atomic diffraction gratings: the atomic Talbot effect. Phys. Rev. A 1995, 51, R14–R17. [Google Scholar] [CrossRef] [PubMed]
  74. Berman, P.R.; Malinovsky, V.S. Principles of Laser Spectroscopy and Quantum Optics; Princeton University Press: Princeton, NJ, USA, 2011. [Google Scholar]
  75. Strekalov, D.V.; Turlapov, A.; Kumarakrishnan, A.; Sleator, T. Periodic structures generated in a cloud of cold atoms. Phys. Rev. A 2002, 66, 023601. [Google Scholar] [CrossRef]
  76. Beattie, S.; Barrett, B.; Weel, M.; Chan, I.; Mok, C.; Cahn, S.B.; Kumarakrishnan, A. Influence of spontaneous emission on a single-state atom interferometer. Phys. Rev. A 2008, 77, 013610. [Google Scholar] [CrossRef]
  77. Barrett, B.; Yavin, I.; Beattie, S.; Kumarakrishnan, A. Numerical simulation of a multilevel atom interferometer. Phys. Rev. A 2010, 82, 023625. [Google Scholar] [CrossRef]
  78. Giltner, D.M.; McGowan, R.W.; Lee, S.A. Atom Interferometer Based on Bragg Scattering from Standing Light Waves. Phys. Rev. Lett. 1995, 75, 2638–2641. [Google Scholar] [CrossRef] [PubMed]
  79. Gradshteyn, I.S.; Ryzhik, I.M. Tables of Integrals, Series, and Products, 7th ed.; Elsevier: Amsterdam, Netherlands, 2007; p. 940. [Google Scholar]
  80. Dubetsky, B.; Berman, P.R.; Sleator, T. Grating stimulated echo. Phys. Rev. A 1992, 46, R2213–R2216. [Google Scholar] [CrossRef] [PubMed]
  81. Gupta, S.; Dieckmann, K.; Hadzibabic, Z.; Pritchard, D.E. Contrast Interferometry using Bose-Einstein Condensates to Measure h/m and α. Phys. Rev. Lett. 2002, 89, 140401. [Google Scholar] [CrossRef] [PubMed]
  82. Jamison, A.O.; Plotkin-Swing, B.; Gupta, S. Advances in precision contrast interferometry with Yb Bose-Einstein condensates. Phys. Rev. A 2014, 90, 063606. [Google Scholar] [CrossRef]
  83. Weidmuller, M.; Hemmerich, A.; Gorlitz, A.; Esslinger, T.; Hansch, T.W. Bragg Diffraction in an Atomic Lattice Bound by Light. Phys. Rev. Lett. 1995, 75, 4583–4586. [Google Scholar] [CrossRef] [PubMed]
  84. Weidemüller, M.; Görlitz, A.; Hänsch, T.; Hemmerich, A. Local and global properties of light-bound atomic lattices investigated by Bragg diffraction. Phys. Rev. A 1998, 58, 4647–4661. [Google Scholar] [CrossRef]
  85. Slama, S.; von Cube, C.; Deh, B.; Ludewig, A.; Zimmermann, C.; Courteille, P.W. Phase-Sensitive Detection of Bragg Scattering at 1D Optical Lattices. Phys. Rev. Lett. 2005, 94, 193901. [Google Scholar] [CrossRef] [PubMed]
  86. Slama, S.; von Cube, C.; Kohler, M.; Zimmermann, C.; Courteille, P.W. Multiple reflections and diffuse scattering in Bragg scattering at optical lattices. Phys. Rev. A 2006, 73, 023424. [Google Scholar] [CrossRef]
  87. Mok, C. Demonstration of Improved Sensitivity of Echo Atom Interferometers to Gravitational Acceleration. Ph.D. Thesis, York University, Toronto, ON, Canada, 2013. [Google Scholar]
  88. Aoyama, T.; Hayakawa, M.; Kinoshita, T.; Nio, M. Revised Value of the Eighth-Order Contribution to the Electron g−2. Phys. Rev. Lett. 2007, 99, 110406. [Google Scholar] [CrossRef] [PubMed]
  89. Pachucki, K.; Yerokhin, V.A.; Cancio Pastor, P. Quantum electrodynamic calculation of the hyperfine structure of 3He. Phys. Rev. A 2012, 85, 042517. [Google Scholar] [CrossRef]
  90. Hanneke, D.; Fogwell, S.; Gabrielse, G. New Measurement of the Electron Magnetic Moment and the Fine Structure Constant. Phys. Rev. Lett. 2008, 100, 120801. [Google Scholar] [CrossRef] [PubMed]
  91. Smiciklas, M.; Shiner, D. Determination of the Fine Structure Constant Using Helium Fine Structure. Phys. Rev. Lett. 2010, 105, 123001. [Google Scholar] [CrossRef] [PubMed]
  92. Jeffery, A.M.; Elmquist, R.E.; Lee, L.H.; Shields, J.Q.; Dziuba, R.F. NIST comparison of the quantized Hall resistance and the realization of the SI OHM through the calculable capacitor. IEEE Trans. Instrum. Meas. 1997, 46, 264. [Google Scholar] [CrossRef]
  93. Small, G.W.; Ricketts, B.W.; Coogan, P.C.; Pritchard, B.J.; Sovierzoski, M.M.R. A new determination of the quantized Hall resistance in terms of the NML calculable cross capacitor. Metrologia 1997, 34, 241. [Google Scholar] [CrossRef]
  94. Battesti, R.; Clade, P.; Guellati-Khélifa, S.; Schwob, C.; Grémaud, B.; Nez, F.; Julien, L.; Biraben, F. Bloch Oscillations of Ultracold Atoms: A Tool for a Metrological Determination of h/mRb. Phys. Rev. Lett. 2004, 92, 253001. [Google Scholar] [CrossRef] [PubMed]
  95. Mohr, P.J.; Newell, D.B.; Taylor, B.N. CODATA Recommended Values of the Fundamental Physical Constants: 2014. 2015; arXiv:1507.07956. [Google Scholar]
  96. Weitz, M.; Young, B.C.; Chu, S. Atomic Interferometer Based on Adiabatic Population Transfer. Phys. Rev. Lett. 1994, 73, 2563–2566. [Google Scholar] [CrossRef] [PubMed]
  97. Müller, H.; Chiow, S.W.; Long, Q.; Vo, C.; Chu, S. A new photon recoil experiment: Towards a determination of the fine structure constant. Appl. Phys. B 2006, 84, 633–642. [Google Scholar] [CrossRef]
  98. Chiow, S.w.; Herrmann, S.; Chu, S.; Müller, H. Noise-Immune Conjugate Large-Area Atom Interferometers. Phys. Rev. Lett. 2009, 103, 050402. [Google Scholar] [CrossRef] [PubMed]
  99. Chiow, S.W.; Kovachy, T.; Chien, H.C.; Kasevich, M.A. 102 k Large Area Atom Interferometers. Phys. Rev. Lett. 2011, 107, 130403. [Google Scholar] [CrossRef] [PubMed]
  100. Beattie, S.; Barrett, B.; Chan, I.; Mok, C.; Yavin, I.; Kumarakrishnan, A. Technique for measuring atomic recoil frequency using coherence functions. Phys. Rev. A 2009, 79, 021605(R). [Google Scholar] [CrossRef]
  101. Wu, S.; Tonyushkin, A.; Prentiss, M.G. Observation of Saturation of Fidelity Decay with an Atom Interferometer. Phys. Rev. Lett. 2009, 103, 034101. [Google Scholar] [CrossRef] [PubMed]
  102. Tonyushkin, A.; Wu, S.; Prentiss, M.G. Demonstration of a multipulse interferometer for quantum kicked-rotor studies. Phys. Rev. A 2009, 79, 051402(R). [Google Scholar] [CrossRef]
  103. Steck, D.A. Rubidium 87 D Line Data, Version 2.1.4, September 2012. Available online: http://steck.us/alkalidata/rubidium87numbers.pdf (accessed on 18 June 2016).
  104. Campbell, G.K.; Leanhardt, A.E.; Mun, J.; Boyd, M.; Streed, E.W.; Ketterle, W.; Pritchard, D.E. Photon Recoil Momentum in Dispersive Media. Phys. Rev. Lett. 2005, 94, 170403. [Google Scholar] [CrossRef] [PubMed]
  105. Steck, D.A. Rubidium 85 D Line Data, version 2.1.4, December 2010. Available online: http://steck.us/alkalidata/rubidium85numbers.pdf (accessed on 18 June 2016).
  106. Andersen, M.F.; Sleator, T. Lattice Interferometer for Laser-Cooled Atoms. Phys. Rev. Lett. 2009, 103, 070402. [Google Scholar] [CrossRef] [PubMed]
  107. Schilke, A.; Zimmermann, C.; Courteille, P.W.; Guerin, W. Photonic Band Gaps in One-Dimensionally Ordered Cold Atomic Vapors. Phys. Rev. Lett. 2011, 106, 223903. [Google Scholar] [CrossRef] [PubMed]
  108. Gundlach, J.H.; Merkowitz, S.M. Measurement of Newton’s Constant Using a Torsion Balance with Angular Acceleration Feedback. Phys. Rev. Lett. 2000, 85, 2869–2872. [Google Scholar] [CrossRef] [PubMed]
  109. Chiaverini, J.; Smullin, S.J.; Geraci, A.A.; Weld, D.M.; Kapitulnik, A. New Experimental Constraints on Non-Newtonian Forces below 100 μm. Phys. Rev. Lett. 2003, 90, 151101. [Google Scholar] [CrossRef] [PubMed]
  110. Prothero, W.A.; Goodkind, J.M. A Superconducting Gravimeter. Rev. Sci. Instrum. 1968, 39, 1257. [Google Scholar] [CrossRef]
  111. Goodkind, J.M. The superconducting gravimeter. Rev. Sci. Instrum. 1999, 70, 4131–4152. [Google Scholar] [CrossRef]
  112. Niebauer, T.M.; Sasagawa, G.S.; Faller, J.E.; Hilt, R.; Klopping, F. A new generation of absolute gravimeters. Metrologia 1995, 32, 159. [Google Scholar] [CrossRef]
  113. Kasevich, M.; Chu, S. Measurement of the gravitational acceleration of an atom with a light-pulse atom interferometer. Appl. Phys. B 1992, 54, 321–332. [Google Scholar] [CrossRef]
  114. Müller, H.; Chiow, S.; Herrmann, S.; Chu, S.; Chung, K. Atom-Interferometry Tests of the Isotropy of Post-Newtonian Gravity. Phys. Rev. Lett. 2008, 100, 031101. [Google Scholar] [CrossRef] [PubMed]
  115. Stockton, J.K.; Takase, K.; Kasevich, M.A. Absolute Geodetic Rotation Measurement Using Atom Interferometry. Phys. Rev. Lett. 2011, 107, 133001. [Google Scholar] [CrossRef] [PubMed]
  116. Yu, N.; Kohel, J.M.; Kellogg, J.R.; Maleki, L. Development of an atom-interferometer gravity gradiometer for gravity measurement from space. Appl. Phys. B 2006, 84, 647–652. [Google Scholar] [CrossRef]
  117. Merlet, S.; Bodart, Q.; Malossi, N.; Landragin, A.; Pereira Dos Santos, F.; Gitlein, O.; Timmen, L. Comparison between two mobile absolute gravimeters: optical versus atomic interferometers. Metrologia 2010, 47, L9. [Google Scholar] [CrossRef]
  118. Schmidt, M.; Senger, A.; Hauth, M.; Freier, C.; Schkolnik, V.; Peters, A. A mobile high-precision absolute gravimeter based on atom interferometry. Gyro. Nav. 2011, 2, 170–177. [Google Scholar] [CrossRef]
  119. Bodart, Q.; Merlet, S.; Malossi, N.; Pereira Dos Santos, F.; Bouyer, P.; Landragin, A. A cold atom pyramidal gravimeter with a single laser beam. Appl. Phys. Lett. 2010, 96, 134101. [Google Scholar] [CrossRef] [Green Version]
  120. Geiger, R.; Menoret, V.; Stern, G.; Zahzam, N.; Cheinet, P.; Battelier, B.; Villing, A.; Moron, F.; Lours, M.; Bidel, Y.; Bresson, A.; Landragin, A.; Bouyer, P. Detecting inertial effects with airborne matter-wave interferometry. Nat. Commun. 2011, 2, 474. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  121. Bidel, Y.; Carraz, O.; Charriere, R.; Cadoret, M.; Zahzam, N.; Bresson, A. Compact cold atom gravimeter for field applications. Appl. Phys. Lett. 2013, 102, 144107. [Google Scholar] [CrossRef]
  122. Poli, N.; Wang, F.Y.; Tarallo, M.G.; Alberti, A.; Prevedelli, M.; Tino, G.M. Precision Measurement of Gravity with Cold Atoms in an Optical Lattice and Comparison with a Classical Gravimeter. Phys. Rev. Lett. 2011, 106, 038501. [Google Scholar] [CrossRef] [PubMed]
  123. Su, E.J.; Wu, S.; Prentiss, M.G. Atom interferometry using wave packets with constant spatial displacements. Phys. Rev. A 2010, 81, 043631. [Google Scholar] [CrossRef]
  124. Beach, R.; Hartmann, S.R.; Friedberg, R. Billiard-ball echo model. Phys. Rev. A 1982, 25, 2658–2666. [Google Scholar] [CrossRef]
  125. Beach, R.; Brody, B.; Hartmann, S.R. Elliptical billiard-ball echo model. Phys. Rev. A 1983, 27, 2537–2547. [Google Scholar] [CrossRef]
  126. Friedberg, R.; Hartmann, S.R. Billiard balls and matter-wave interferometry. Phys. Rev. A 1993, 48, 1446–1472. [Google Scholar] [CrossRef] [PubMed]
  127. Raman, C.V.; Nagendra Nathe, N.S. The diffraction of light by high frequency sound waves: Part I. Proc. Indian Acad. Sci. (Math. Sci.) 1935, 2, 406–412. [Google Scholar]
  128. Raman, C.V.; Nagendra Nathe, N.S. The diffraction of light by high frequency sound waves: Part II. Proc. Indian Acad. Sci. (Math. Sci.) 1935, 2, 413–420. [Google Scholar]
  129. Kurnit, N.A.; Abella, I.D.; Hartmann, S.R. Observation of a Photon Echo. Phys. Rev. Lett. 1964, 13, 567–568. [Google Scholar] [CrossRef]
  130. Mossberg, T.W.; Kachru, R.; Hartmann, S.R.; Flusberg, A.M. Echoes in gaseous media: A generalized theory of rephasing phenomena. Phys. Rev. A 1979, 20, 1976–1996. [Google Scholar] [CrossRef]
  131. Vorozcovs, A.; Weel, M.; Beattie, S.; Cauchi, S.; Kumarakrishnan, A. Measurements of temperature scaling laws in an optically dense magneto-optical trap. J. Opt. Soc. Am. B 2005, 22, 943–950. [Google Scholar] [CrossRef]
  132. Scholten, R.E. Enhanced laser shutter using a hard disk drive rotary voice-coil actuator. Rev. Sci. Instrum. 2007, 78, 026101. [Google Scholar] [CrossRef] [PubMed]
  133. Niebauer, T.M.; Schiel, A.; van Westrum, D. Complex heterodyne for undersampled chirped sinusoidal signals. Appl. Opt. 2006, 45, 8322–8330. [Google Scholar] [CrossRef] [PubMed]
  134. Cheinet, P.; Canuel, B.; Pereira Dos Santos, F.; Gauguet, A.; Yver-Leduc, F.; Landragin, A. Measurement of the Sensitivity Function in a Time-Domain Atomic Interferometer. IEEE Trans. 2008, 57, 1141–1148. [Google Scholar] [CrossRef]
  135. Smith, G.A.; Silberfarb, A.; Deutsch, I.H.; Jessen, P.S. Efficient Quantum-State Estimation by Continuous Weak Measurement and Dynamical Control. Phys. Rev. Lett. 2006, 97, 180403. [Google Scholar] [CrossRef] [PubMed]
  136. Deutsch, I.H.; Jessen, P.S. Quantum-state control in optical lattices. Phys. Rev. A 1998, 57, 1972–1986. [Google Scholar] [CrossRef]
  137. Geremia, J.; Stockton, J.K.; Mabuchi, H. Suppression of Spin Projection Noise in Broadband Atomic Magnetometry. Phys. Rev. Lett. 2005, 94, 203002. [Google Scholar] [CrossRef] [PubMed]
  138. Chaudhury, S.; Smith, G.A.; Schulz, K.; Jessen, P.S. Continuous Nondemolition Measurement of the Cs Clock Transition Pseudospin. Phys. Rev. Lett. 2006, 96, 043001. [Google Scholar] [CrossRef] [PubMed]
  139. van Enk, S.J.; Lütkenhaus, N.; Kimble, H.J. Experimental procedures for entanglement verification. Phys. Rev. A 2007, 75, 052318. [Google Scholar] [CrossRef]
  140. Allred, J.C.; Lyman, R.N.; Kornack, T.W.; Romalis, M.V. High-Sensitivity Atomic Magnetometer Unaffected by Spin-Exchange Relaxation. Phys. Rev. Lett. 2002, 89, 130801. [Google Scholar] [CrossRef] [PubMed]
  141. Vengalattore, M.; Higbie, J.M.; Leslie, S.R.; Guzman, J.; Sadler, L.E.; Stamper-Kurn, D.M. High-Resolution Magnetometry with a Spinor Bose-Einstein Condensate. Phys. Rev. Lett. 2007, 98, 200801. [Google Scholar] [CrossRef] [PubMed]
  142. Tamm, C.; Buhr, E.; Mlynek, J. Raman heterodyne studies of velocity diffusion effects in radio-frequency-laser double resonance. Phys. Rev. A 1986, 34, 1977–1994. [Google Scholar] [CrossRef]
  143. Buhr, E.; Mlynek, J. Collision-induced Ramsey resonances in Sm vapor. Phys. Rev. A 1987, 36, 2684–2705. [Google Scholar] [CrossRef]
  144. Suter, D.; Rosatzin, M.; Mlynek, J. Optically induced coherence transfer echoes between Zeeman substates. Phys. Rev. Lett. 1991, 67, 34–37. [Google Scholar] [CrossRef] [PubMed]
  145. Suter, D.; Klepel, H.; Mlynek, J. Time-resolved two-dimensional spectroscopy of optically driven atomic sublevel coherences. Phys. Rev. Lett. 1991, 67, 2001–2004. [Google Scholar] [CrossRef] [PubMed]
  146. Dubetsky, B.; Berman, P.R. Magnetic Grating Free Induction Decay and Magnetic Grating Echo. Laser Phys. 1994, 4, 1017–1029. [Google Scholar]
  147. Kumarakrishnan, A.; Cahn, S.B.; Shim, U.; Sleator, T. Magnetic grating echoes from laser-cooled atoms. Phys. Rev. A 1998, 58, R3387–R3390. [Google Scholar] [CrossRef]
  148. Kumarakrishnan, A.; Shim, U.; Cahn, S.B.; Sleator, T. Ground-state grating echoes from Rb vapor at room temperature. Phys. Rev. A 1998, 58, 3868–3872. [Google Scholar] [CrossRef]
  149. Abella, I.D.; Kurnit, N.A.; Hartmann, S.R. Photon Echoes. Phys. Rev. 1966, 141, 391–406. [Google Scholar] [CrossRef]
  150. Rotberg, E.A.; Barrett, B.; Beattie, S.; Chudasama, S.; Weel, M.; Chan, I.; Kumarakrishnan, A. Measurement of excited-state lifetime using two-pulse photon echoes in rubidium vapor. J. Opt. Soc. Am. B 2007, 24, 671–680. [Google Scholar] [CrossRef]
  151. Shim, U.; Kumarakrishnan, A.; Cahn, S.B.; Sleator, T. Collisional revival of magnetic grating free-induction decay. In Proceedings of the 1998 International Quantum Electronics Conference; 1998 Osa Technical Digest Series; Institute of Electrical & Electronics Enginee: Piscataway, NJ, USA, 1998; p. 201. [Google Scholar]
  152. Cahn, S.B.; Kumarakrishnan, A.; Shim, U.; Sleator, T. Phase Space Diagnostics of Trapped Atoms By Magnetic Ground-State Manipulation. Bull. Am. Phys. Soc. 1997, 42, 1057. [Google Scholar]
  153. Tonyushkin, A.; Kumarakrishnan, A.; Turlapov, A.; Sleator, T. Magnetic coherence gratings in a high-flux atomic beam. Eur. Phys. J. D 2010, 58, 39–46. [Google Scholar] [CrossRef]
  154. Shore, B.W. The Theory of Coherent Atomic Excitation; Wiley: New York, NY, USA, 1990. [Google Scholar]
  155. Edmonds, A.R. Angular Momentum in Quantum Mechanics; Princeton University Press: Princeton, NJ, USA, 1996. [Google Scholar]
  156. Rochester, S.M.; Budker, D. Atomic polarization visualized. Am. J. Phys. 2001, 69, 450. [Google Scholar] [CrossRef]
  157. Berman, P.R.; Rogers, G.; Dubetsky, B. Rate equations between electronic-state manifolds. Phys. Rev. A 1993, 48, 1506–1513. [Google Scholar] [CrossRef] [PubMed]
  158. Berman, P.R. Nonlinear spectroscopy and laser cooling. Phys. Rev. A 1991, 43, 1470–1483. [Google Scholar] [CrossRef] [PubMed]
  159. Cohen-Tannoudji, C.; Kastler, A. Optical pumping. Prog. Opt. 1966, 5, 169–249. [Google Scholar]
  160. White, C.W.; Hughes, W.M.; Hayne, G.S.; Robinson, H.G. Determination of g-Factor Ratios for Free Rb85 and Rb87 Atoms. Phys. Rev. 1968, 174, 23–32. [Google Scholar] [CrossRef]
  161. Anthony, J.M.; Sebastian, K.J. Relativistic corrections to the Zeeman effect in hydrogen-like atoms and positronium. Phys. Rev. A 1994, 49, 192–206. [Google Scholar] [CrossRef] [PubMed]
  162. DiSciacca, J.; Marshall, M.; Marable, K.; Gabrielse, G.; Ettenauer, S.; Tardiff, E.; Kalra, R.; Fitzakerley, D.W.; George, M.C.; Hessels, E.A.; et al. One-Particle Measurement of the Antiproton Magnetic Moment. Phys. Rev. Lett. 2013, 110, 130801. [Google Scholar] [CrossRef] [PubMed]
  163. Suter, D.; Mlynek, J. Dynamics of atomic sublevel coherences during modulated optical pumping. Phys. Rev. A 1991, 43, 6124–6134. [Google Scholar] [CrossRef] [PubMed]
  164. Sané, S.S.; Bennetts, S.; Debs, J.E.; Kuhn, C.C.N.; McDonald, G.D.; Altin, P.A.; Close, J.D.; Robins, N.P. 11 W narrow linewidth laser source at 780nm for laser cooling and manipulation of Rubidium. Opt. Express 2012, 20, 8915–8919. [Google Scholar] [CrossRef] [PubMed]
  165. Chiow, S.W.; Kovachy, T.; Hogan, J.M.; Kasevich, M.A. Generation of 43 W of quasi-continuous 780 nm laser light via high-efficiency, single-pass frequency doubling in periodically poled lithium niobate crystals. Opt. Lett. 2012, 37, 3861–3863. [Google Scholar] [CrossRef] [PubMed]
  166. Ricci, L.; Weidemuller, M.; Esslinger, T.; Hemmerich, A.; Zimmermann, C.; Vuletic, V.; Konig, W.; Hansch, T.W. A compact grating-stabilized diode laser system for atomic physics. Opt. Comm. 1995, 117, 541–549. [Google Scholar] [CrossRef]
  167. Baillard, X.; Gauguet, A.; Bize, S.; Lemonde, P.; Laurent, P.; Clairon, A.; Rosenbusch, P. Interference-filter-stabilized external-cavity diode lasers. Opt. Comm. 2006, 266, 609–613. [Google Scholar] [CrossRef]
  168. Gilowski, M.; Schubert, C.; Zaiser, M.; Herr, W.; Wubbena, T.; Wendrich, T.; Muller, T.; Rasel, E.M.; Ertmer, W. Narrow bandwidth interference filter-stabilized diode laser systems for the manipulation of neutral atoms. Opt. Comm. 2007, 280, 443–447. [Google Scholar] [CrossRef]
  169. Carr, A.V.; Sechrest, Y.H.; Waitukaitis, S.R.; Perreault, J.D.; Lonij, V.P.A.; Cronin, A.D. Cover slip external cavity diode laser. Rev. Sci. Instrum. 2007, 78, 106108. [Google Scholar] [CrossRef] [PubMed]
  170. Saliba, S.D.; Junker, M.; Turner, L.D.; Scholten, R.E. Mode stability of external cavity diode lasers. Appl. Opt. 2009, 48, 6692–6700. [Google Scholar] [CrossRef] [PubMed]
  171. Libbrecht, K.G.; Hall, J.L. A low-noise high-speed diode laser current controller. Rev. Sci. Instrum. 1993, 64, 2133–2135. [Google Scholar] [CrossRef]
  172. Nyman, R.A.; Varoquaux, G.; Villier, B.; Sacchet, D.; Moron, F.; Coq, Y.L.; Aspect, A.; Bouyer, P. Tapered-amplified antireflection-coated laser diodes for potassium and rubidium atomic-physics experiments. Rev. Sci. Instrum. 2006, 77, 033105. [Google Scholar] [CrossRef]
  173. Takase, K.; Stockton, J.K.; Kasevich, M.A. High-power pulsed-current-mode operation of an overdriven tapered amplifier. Opt. Lett. 2007, 32, 2617–2619. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Raman and echo-type atom interferometer schemes. (a) Three-pulse Mach-Zehnder configuration of a Raman interferometer for measuring inertial effects such as gravity. Each light pulse is composed of two counter-propagating beams of frequencies ω 1 and ω 2 that induced velocity-sensitive Raman transitions between states 1 , p (blue lines) and 2 , p + k eff (green lines). The phase difference between these two beams during each pulse ( φ j ) is imprinted on the wavepackets wherever a change of momentum takes place. In the absence of any accelerations, the final atomic populations at the output of the interferometer oscillate sinusoidally as a function of the total laser phase difference Δ φ = φ 1 2 φ 2 + φ 3 . (b) Two-pulse configuration of the grating-echo AI, which is sensitive both to inertial effects and to the atomic recoil frequency. Each light pulse is composed of counter-propagating beams of the same frequency ω and polarization, which create a standing wave with phase φ j . The two SW pulses, separated by a time T, diffract atoms in the same internal state 1 , p into a superposition of momentum states 1 , p + n q separated by integer multiples of the two-photon momentum q = 2 k . At time t = 2 T , a subset of these momentum states interfere and create a spatially modulated density grating with period λ / 2 . A traveling-wave read-out pulse is applied at this time, and the back-scattered “echo” signal from the atoms is measured. The phase of this back-scattered light is proportional to the phase shift of the interference pattern due to gravity. Similarly, the back-scattered light intensity is modulated as a function of T at the atomic recoil frequency.
Figure 1. Raman and echo-type atom interferometer schemes. (a) Three-pulse Mach-Zehnder configuration of a Raman interferometer for measuring inertial effects such as gravity. Each light pulse is composed of two counter-propagating beams of frequencies ω 1 and ω 2 that induced velocity-sensitive Raman transitions between states 1 , p (blue lines) and 2 , p + k eff (green lines). The phase difference between these two beams during each pulse ( φ j ) is imprinted on the wavepackets wherever a change of momentum takes place. In the absence of any accelerations, the final atomic populations at the output of the interferometer oscillate sinusoidally as a function of the total laser phase difference Δ φ = φ 1 2 φ 2 + φ 3 . (b) Two-pulse configuration of the grating-echo AI, which is sensitive both to inertial effects and to the atomic recoil frequency. Each light pulse is composed of counter-propagating beams of the same frequency ω and polarization, which create a standing wave with phase φ j . The two SW pulses, separated by a time T, diffract atoms in the same internal state 1 , p into a superposition of momentum states 1 , p + n q separated by integer multiples of the two-photon momentum q = 2 k . At time t = 2 T , a subset of these momentum states interfere and create a spatially modulated density grating with period λ / 2 . A traveling-wave read-out pulse is applied at this time, and the back-scattered “echo” signal from the atoms is measured. The phase of this back-scattered light is proportional to the phase shift of the interference pattern due to gravity. Similarly, the back-scattered light intensity is modulated as a function of T at the atomic recoil frequency.
Atoms 04 00019 g001
Figure 2. (a) An example of two, low-order trajectories that contribute to the two-pulse atom interferometer (AI) (SWj = j th SW pulse, RO = read-out pulse). Three momentum states are shown ( p 0 , p 0 + q , and p 0 + 2 q ) corresponding to the solid, dashed, and dotted lines, respectively. The contrast of the interference fringe created at the echo time t = 2 T is modulated sinusoidally by a phase 2 ω q T . (b) Low-order trajectories for the modified three-pulse AI. Here, a third pulse is applied at t = δ T between SW1 and SW2 which further modulates the phase of the interference by 2 ω q δ T . (c) Echo energy as a function of δ T / τ q for the modified AI. Line shapes are shown for three different pulse areas, u 3 = Ω eff τ 3 , to illustrate the effect of fringe narrowing that occurs for increasing interaction strength.
Figure 2. (a) An example of two, low-order trajectories that contribute to the two-pulse atom interferometer (AI) (SWj = j th SW pulse, RO = read-out pulse). Three momentum states are shown ( p 0 , p 0 + q , and p 0 + 2 q ) corresponding to the solid, dashed, and dotted lines, respectively. The contrast of the interference fringe created at the echo time t = 2 T is modulated sinusoidally by a phase 2 ω q T . (b) Low-order trajectories for the modified three-pulse AI. Here, a third pulse is applied at t = δ T between SW1 and SW2 which further modulates the phase of the interference by 2 ω q δ T . (c) Echo energy as a function of δ T / τ q for the modified AI. Line shapes are shown for three different pulse areas, u 3 = Ω eff τ 3 , to illustrate the effect of fringe narrowing that occurs for increasing interaction strength.
Atoms 04 00019 g002
Figure 3. (a) Optical setup for the interferometer. The glass cell has dimensions 7 . 6 × 7 . 6 × 84 cm and is oriented along the vertical. (b) Timing diagram for the AI. The gate acousto-optic modulator (AOM) is pulsed on to allow light for each excitation pulse produced by the k 1 and k 2 AOMs. The pulse occurring at t = T 1 + δ T corresponds to the third standing wave (SW) pulse. The read-out pulse (which is independent of the gate AOM) and the photo-multiplier tube (PMT) gate are turned on for ∼9 μs in the vicinity of the echo time, t = T 1 + 2 T . (c) Example of a two-pulse grating-echo signal (from a 10 μK 87 Rb sample) recorded by the PMT, which corresponds to an echo energy of 130 pJ. AI pulse spacing: T = 1 . 06338 ms; pulse durations: τ 1 = 3 . 8 μs, τ 2 = 1 . 2 μs; AI and read-out beam detunings: Δ AI = 220 MHz, Δ RO = 40 MHz; AI and read-out beam intensity: I∼ 40 mW/cm 2 .
Figure 3. (a) Optical setup for the interferometer. The glass cell has dimensions 7 . 6 × 7 . 6 × 84 cm and is oriented along the vertical. (b) Timing diagram for the AI. The gate acousto-optic modulator (AOM) is pulsed on to allow light for each excitation pulse produced by the k 1 and k 2 AOMs. The pulse occurring at t = T 1 + δ T corresponds to the third standing wave (SW) pulse. The read-out pulse (which is independent of the gate AOM) and the photo-multiplier tube (PMT) gate are turned on for ∼9 μs in the vicinity of the echo time, t = T 1 + 2 T . (c) Example of a two-pulse grating-echo signal (from a 10 μK 87 Rb sample) recorded by the PMT, which corresponds to an echo energy of 130 pJ. AI pulse spacing: T = 1 . 06338 ms; pulse durations: τ 1 = 3 . 8 μs, τ 2 = 1 . 2 μs; AI and read-out beam detunings: Δ AI = 220 MHz, Δ RO = 40 MHz; AI and read-out beam intensity: I∼ 40 mW/cm 2 .
Atoms 04 00019 g003
Figure 4. (a) Demonstration of an individual recoil measurement in 85 Rb using the modified three-pulse AI. Data were recorded in two temporal windows separated by T = 1132 τ q 36 . 67 ms. A least-squares fit to the data yields a statistical uncertainty in ω q of 180 ppb. Inset: expanded view of the fringe near δ T = 64 μs. (b) 82 chronological measurements of ω q in 87 Rb using T 45 . 5 ms. Each measurement was acquired in 10 minutes and produced a typical statistical error of ∼380 ppb. We have scaled these results by the expected value of the recoil frequency, ω q ( 0 ) = 94 . 77384783 ( 12 ) rad/ms, which is based on the value of h / M ( 87 Rb) from Ref. [15] and the F = 2 F = 3 transition frequency in 87 Rb from Ref. [103]. The dashed grid lines indicate the weighted standard deviation of 339 ppb, and the standard deviation of the mean is 37 ppb. The corresponding reduced chi-squared is χ 2 / dof = 0 . 93 for dof = 81 degrees of freedom. The mean value, shown by the solid grid line, is ∼2.8 ppm below the expected value, which is due to systematic effects. AI pulse parameters: T = 45 . 4837 ms, τ 1 = 2 . 2 μs, τ 2 = 1 . 4 μs, τ 3 = 3 μs, Δ AI = 219 . 8 MHz, Δ RO ∼ 40 MHz, I∼ 95 mW/cm 2 .
Figure 4. (a) Demonstration of an individual recoil measurement in 85 Rb using the modified three-pulse AI. Data were recorded in two temporal windows separated by T = 1132 τ q 36 . 67 ms. A least-squares fit to the data yields a statistical uncertainty in ω q of 180 ppb. Inset: expanded view of the fringe near δ T = 64 μs. (b) 82 chronological measurements of ω q in 87 Rb using T 45 . 5 ms. Each measurement was acquired in 10 minutes and produced a typical statistical error of ∼380 ppb. We have scaled these results by the expected value of the recoil frequency, ω q ( 0 ) = 94 . 77384783 ( 12 ) rad/ms, which is based on the value of h / M ( 87 Rb) from Ref. [15] and the F = 2 F = 3 transition frequency in 87 Rb from Ref. [103]. The dashed grid lines indicate the weighted standard deviation of 339 ppb, and the standard deviation of the mean is 37 ppb. The corresponding reduced chi-squared is χ 2 / dof = 0 . 93 for dof = 81 degrees of freedom. The mean value, shown by the solid grid line, is ∼2.8 ppm below the expected value, which is due to systematic effects. AI pulse parameters: T = 45 . 4837 ms, τ 1 = 2 . 2 μs, τ 2 = 1 . 4 μs, τ 3 = 3 μs, Δ AI = 219 . 8 MHz, Δ RO ∼ 40 MHz, I∼ 95 mW/cm 2 .
Atoms 04 00019 g004
Figure 5. Relative correction to the recoil frequency due to the refractive index as a function of the detuning of the excitation field, Δ AI . These curves are based on Equation (19) with a density of ρ = 10 10 atoms/cm 3 . Predictions for both 85 Rb and 87 Rb are shown. The detuning is plotted with respect to the F = 3 F = 4 transition in 85 Rb, and the F = 2 F = 3 transition in 87 Rb. The dashed grid lines label the location of excited states [103,105]. The “magic” frequencies, where the relative correction crosses zero, are indicated with arrows at Δ AI 66 . 4 MHz for 85 Rb and at Δ AI 162 . 6 MHz for 87 Rb.
Figure 5. Relative correction to the recoil frequency due to the refractive index as a function of the detuning of the excitation field, Δ AI . These curves are based on Equation (19) with a density of ρ = 10 10 atoms/cm 3 . Predictions for both 85 Rb and 87 Rb are shown. The detuning is plotted with respect to the F = 3 F = 4 transition in 85 Rb, and the F = 2 F = 3 transition in 87 Rb. The dashed grid lines label the location of excited states [103,105]. The “magic” frequencies, where the relative correction crosses zero, are indicated with arrows at Δ AI 66 . 4 MHz for 85 Rb and at Δ AI 162 . 6 MHz for 87 Rb.
Atoms 04 00019 g005
Figure 6. Recoil diagrams for (a) two-pulse and (b) three-pulse AIs in the absence of gravity. Only a subset of all trajectories are shown. SW refers to standing wave pulses and RO is a traveling wave read-out pulse. The SW pulses, composed of two counter-propagating traveling wave components, each with wave vector | k | , diffract atoms into a superposition of momentum states separated by q . For both two-pulse and three-pulse AIs, the backscattered signal arises from interferences between states differing by q at the echo time.
Figure 6. Recoil diagrams for (a) two-pulse and (b) three-pulse AIs in the absence of gravity. Only a subset of all trajectories are shown. SW refers to standing wave pulses and RO is a traveling wave read-out pulse. The SW pulses, composed of two counter-propagating traveling wave components, each with wave vector | k | , diffract atoms into a superposition of momentum states separated by q . For both two-pulse and three-pulse AIs, the backscattered signal arises from interferences between states differing by q at the echo time.
Atoms 04 00019 g006
Figure 7. Schematic of experimental setup. (a) Block diagram of laser sources and frequency control. The frequency of the laser light is defined with respect to ω 0 , the resonant frequency of the 87 Rb F = 2 F = 3 transition or the 85 Rb F = 3 F = 4 transition. (b) Schematic of the atom interferometry setup. The lower-detection optics and the upper corner-cube reflector are anchored together and placed on a vibration isolation platform. The vacuum chamber and vibration isolation platform rest on an optical table supported by pneumatic legs. The photodiodes detect a 250 MHz beat note, which is the frequency difference between the AI and LO beams. The forms of the ion pump, cancelling coils, and anti-Helmholtz trapping coils are not shown.
Figure 7. Schematic of experimental setup. (a) Block diagram of laser sources and frequency control. The frequency of the laser light is defined with respect to ω 0 , the resonant frequency of the 87 Rb F = 2 F = 3 transition or the 85 Rb F = 3 F = 4 transition. (b) Schematic of the atom interferometry setup. The lower-detection optics and the upper corner-cube reflector are anchored together and placed on a vibration isolation platform. The vacuum chamber and vibration isolation platform rest on an optical table supported by pneumatic legs. The photodiodes detect a 250 MHz beat note, which is the frequency difference between the AI and LO beams. The forms of the ion pump, cancelling coils, and anti-Helmholtz trapping coils are not shown.
Atoms 04 00019 g007
Figure 8. Predicted shape of the in-phase component of the echo envelope at (a) 2 T 21 = 0 . 1 ms for the two-pulse AI and (b) 2 T 21 + T 32 = 100 ms for the three-pulse AI are shown as solid red lines. Here, g = 9 . 8 m/s 2 , and Δ t is the time measured with respect to the echo time. In (a), a time slice used in the analysis is shown as a gray rectangle. The two echo envelopes shown in red are described by E D ( 2 ) e i ϕ D ( 2 ) for the two-pulse AI (see Equations (25) and (27)), and E D ( 3 ) e i ϕ D ( 3 ) for the three-pulse AI (see Equations (34) and (36)). The black dashed lines show the envelopes in the absence of gravity. The echo envelope exhibits an increase in the oscillation frequency as the free-fall time increases. (c) Predicted shape of the in-phase component of the signal amplitude for the two-pulse AI as a function of T 21 as predicted by E AI ( 2 ) e i ϕ AI ( 2 ) (see Equations (26) and (28)), and shown as a solid red line. The signal exhibits chirped sinusoidal behaviour due to the quadratic dependence of ϕ AI ( 2 ) on T 21 , and also shows recoil modulation with a period τ q = 33 . 151 μ s. Here, we set g = 980 m/s 2 for illustrative purposes. The black dashed lines show the total signal amplitude. (d) Predicted shape of the in-phase component of the signal amplitude for the three-pulse AI as a function of T 32 as predicted by E AI ( 3 ) e i ϕ AI ( 3 ) (see Equations (35) and (37)), shown as a solid red line. The signal exhibits a constant modulation frequency with a period τ v and shows no recoil modulation. Here, T 21 = 15 ms. Again, we set g = 980 m/s 2 for illustrative purposes. The total signal amplitudes for both two-pulse and three-pulse AIs show a phenomenological exponential decay (gray line) to illustrate signal loss due to decoherence and transit time effects.
Figure 8. Predicted shape of the in-phase component of the echo envelope at (a) 2 T 21 = 0 . 1 ms for the two-pulse AI and (b) 2 T 21 + T 32 = 100 ms for the three-pulse AI are shown as solid red lines. Here, g = 9 . 8 m/s 2 , and Δ t is the time measured with respect to the echo time. In (a), a time slice used in the analysis is shown as a gray rectangle. The two echo envelopes shown in red are described by E D ( 2 ) e i ϕ D ( 2 ) for the two-pulse AI (see Equations (25) and (27)), and E D ( 3 ) e i ϕ D ( 3 ) for the three-pulse AI (see Equations (34) and (36)). The black dashed lines show the envelopes in the absence of gravity. The echo envelope exhibits an increase in the oscillation frequency as the free-fall time increases. (c) Predicted shape of the in-phase component of the signal amplitude for the two-pulse AI as a function of T 21 as predicted by E AI ( 2 ) e i ϕ AI ( 2 ) (see Equations (26) and (28)), and shown as a solid red line. The signal exhibits chirped sinusoidal behaviour due to the quadratic dependence of ϕ AI ( 2 ) on T 21 , and also shows recoil modulation with a period τ q = 33 . 151 μ s. Here, we set g = 980 m/s 2 for illustrative purposes. The black dashed lines show the total signal amplitude. (d) Predicted shape of the in-phase component of the signal amplitude for the three-pulse AI as a function of T 32 as predicted by E AI ( 3 ) e i ϕ AI ( 3 ) (see Equations (35) and (37)), shown as a solid red line. The signal exhibits a constant modulation frequency with a period τ v and shows no recoil modulation. Here, T 21 = 15 ms. Again, we set g = 980 m/s 2 for illustrative purposes. The total signal amplitudes for both two-pulse and three-pulse AIs show a phenomenological exponential decay (gray line) to illustrate signal loss due to decoherence and transit time effects.
Atoms 04 00019 g008
Figure 9. (a) Example of fit to the in-phase component of the two-pulse echo signal obtained on a single acquisition for 2 T 21 = 9 . 3 ms. (b) Example of fit to the in-phase component of the three-pulse echo signal obtained on a single acquisition for 2 T 21 + T 32 = 45 . 1 ms.
Figure 9. (a) Example of fit to the in-phase component of the two-pulse echo signal obtained on a single acquisition for 2 T 21 = 9 . 3 ms. (b) Example of fit to the in-phase component of the three-pulse echo signal obtained on a single acquisition for 2 T 21 + T 32 = 45 . 1 ms.
Atoms 04 00019 g009
Figure 10. (a–d) Four observational windows for the two-pulse AI experiment showing the amplitude of the in-phase component as a function of T 21 . The frequency of the signal increases linearly as a function of T 21 ( ϕ AI ( 2 ) / T 21 = 2 q g T 21 ). The solid red line indicates a simultaneous least-squares fit to all four windows. These results were obtained using 85 Rb, with SW pulse lengths τ 1 = 800 ns and τ 2 = 90 ns, optical intensity of 50 mW/cm 2 , at a detuning of Δ = 55 MHz relative to the F = 3 4 cycling transition.
Figure 10. (a–d) Four observational windows for the two-pulse AI experiment showing the amplitude of the in-phase component as a function of T 21 . The frequency of the signal increases linearly as a function of T 21 ( ϕ AI ( 2 ) / T 21 = 2 q g T 21 ). The solid red line indicates a simultaneous least-squares fit to all four windows. These results were obtained using 85 Rb, with SW pulse lengths τ 1 = 800 ns and τ 2 = 90 ns, optical intensity of 50 mW/cm 2 , at a detuning of Δ = 55 MHz relative to the F = 3 4 cycling transition.
Atoms 04 00019 g010
Figure 11. Two observational windows of the amplitude of the in-phase component for the three-pulse AI with T 21 = 7 . 4319 ms as a function of T 32 in the vicinity of (a) T 32 = 0 . 2 ms and (b) T 32 = 30 ms. These data correspond to a single time slice and exhibit a constant frequency as a function of T 32 ( ϕ AI ( 3 ) / T 32 = q g T 21 ). The solid red line indicates a simultaneous least-squares fit to both data windows, which yields a frequency of 187324.75(8) Hz. An analysis of the fit residuals indicates the standard deviation of phase noise for the entire data set is 0.2 rad. These results were obtained using 87 Rb.
Figure 11. Two observational windows of the amplitude of the in-phase component for the three-pulse AI with T 21 = 7 . 4319 ms as a function of T 32 in the vicinity of (a) T 32 = 0 . 2 ms and (b) T 32 = 30 ms. These data correspond to a single time slice and exhibit a constant frequency as a function of T 32 ( ϕ AI ( 3 ) / T 32 = q g T 21 ). The solid red line indicates a simultaneous least-squares fit to both data windows, which yields a frequency of 187324.75(8) Hz. An analysis of the fit residuals indicates the standard deviation of phase noise for the entire data set is 0.2 rad. These results were obtained using 87 Rb.
Atoms 04 00019 g011
Figure 12. (a) Measurement of ω a ( 3 ) across the echo envelope using 200 time slices. Here, we use T 21 = 7 . 431900 ms. The intercept of the line is 188192.095(17) Hz, giving a value of g = 9 . 8774458 ( 9 ) m/s 2 and a corresponding statistical error of 90 ppb (not corrected for systematic effects). Similar analysis of the in-quadrature component reduces the combined statistical error to 75 ppb. (b) Residuals of the linear fit. These results were obtained using 87 Rb.
Figure 12. (a) Measurement of ω a ( 3 ) across the echo envelope using 200 time slices. Here, we use T 21 = 7 . 431900 ms. The intercept of the line is 188192.095(17) Hz, giving a value of g = 9 . 8774458 ( 9 ) m/s 2 and a corresponding statistical error of 90 ppb (not corrected for systematic effects). Similar analysis of the in-quadrature component reduces the combined statistical error to 75 ppb. (b) Residuals of the linear fit. These results were obtained using 87 Rb.
Atoms 04 00019 g012
Figure 13. Improved apparatus to measure of g using an echo AI.
Figure 13. Improved apparatus to measure of g using an echo AI.
Atoms 04 00019 g013
Figure 14. (a) Laser pulses along k 1 and k 2 excite a sample of laser-cooled atoms; θ∼ 10 mrad. (b) Level diagram for the experiment with k 1 and k 2 having orthogonal linear polarizations; detuning Δ∼ 40 MHz. (c) Timing diagram for the magnetic grating echo (MGE) signal.
Figure 14. (a) Laser pulses along k 1 and k 2 excite a sample of laser-cooled atoms; θ∼ 10 mrad. (b) Level diagram for the experiment with k 1 and k 2 having orthogonal linear polarizations; detuning Δ∼ 40 MHz. (c) Timing diagram for the magnetic grating echo (MGE) signal.
Atoms 04 00019 g014
Figure 15. Ground-state coherence signal from dual-isotope magneto-optical trap (MOT) showing multiple-frequency components in the signal. Only one eighth of the points are displayed so that the fit (red line) and the data (black dots) can be distinguished. From a single shot, the ratio of effective atomic g-factors r = g F 87 / g F 85 is typically determined to a precision of ∼1.5 ppm.
Figure 15. Ground-state coherence signal from dual-isotope magneto-optical trap (MOT) showing multiple-frequency components in the signal. Only one eighth of the points are displayed so that the fit (red line) and the data (black dots) can be distinguished. From a single shot, the ratio of effective atomic g-factors r = g F 87 / g F 85 is typically determined to a precision of ∼1.5 ppm.
Atoms 04 00019 g015
Figure 16. (a) Correction signal (red) and lock signal (blue) obtained with laser auto-locked to an iodine transition. The laser head was maintained at atmospheric pressure so that pressure-induced drifts of the correction signal extend over nearly 100 MHz. (b) Same signals as in part (a) with the laser head pumped down to 1 mTorr. The correction signal is uncorrelated with pressure changes and the range of frequency drifts shows a tenfold reduction. (c) Typical histogram of lock signals (blue curves in parts (a) and (b)) exhibits a full width at half maximum of 500 kHz, which is a measure of lock stability.
Figure 16. (a) Correction signal (red) and lock signal (blue) obtained with laser auto-locked to an iodine transition. The laser head was maintained at atmospheric pressure so that pressure-induced drifts of the correction signal extend over nearly 100 MHz. (b) Same signals as in part (a) with the laser head pumped down to 1 mTorr. The correction signal is uncorrelated with pressure changes and the range of frequency drifts shows a tenfold reduction. (c) Typical histogram of lock signals (blue curves in parts (a) and (b)) exhibits a full width at half maximum of 500 kHz, which is a measure of lock stability.
Atoms 04 00019 g016
Figure 17. Allan deviation of an original equipment manufacturer (OEM) current controller showing shot-noise-limited behavior for τ < 0 . 3 s. The minimum Allan deviation is 4 . 2 × 10 6 at τ = 0 . 34 s. The long-term Allan deviation is 10 . 7 × 10 6 . This data was acquired with bandwidth filter cut-offs at 0.01 Hz and 10 kHz. The typical rms deviation in the DC current is 0.32 μA. The current noise density measured with a 10 kHz bandwidth is 0.35 μA/ Hz .
Figure 17. Allan deviation of an original equipment manufacturer (OEM) current controller showing shot-noise-limited behavior for τ < 0 . 3 s. The minimum Allan deviation is 4 . 2 × 10 6 at τ = 0 . 34 s. The long-term Allan deviation is 10 . 7 × 10 6 . This data was acquired with bandwidth filter cut-offs at 0.01 Hz and 10 kHz. The typical rms deviation in the DC current is 0.32 μA. The current noise density measured with a 10 kHz bandwidth is 0.35 μA/ Hz .
Atoms 04 00019 g017

Share and Cite

MDPI and ACS Style

Barrett, B.; Carew, A.; Beica, H.C.; Vorozcovs, A.; Pouliot, A.; Kumarakrishnan, A. Prospects for Precise Measurements with Echo Atom Interferometry. Atoms 2016, 4, 19. https://doi.org/10.3390/atoms4030019

AMA Style

Barrett B, Carew A, Beica HC, Vorozcovs A, Pouliot A, Kumarakrishnan A. Prospects for Precise Measurements with Echo Atom Interferometry. Atoms. 2016; 4(3):19. https://doi.org/10.3390/atoms4030019

Chicago/Turabian Style

Barrett, Brynle, Adam Carew, Hermina C. Beica, Andrejs Vorozcovs, Alexander Pouliot, and A. Kumarakrishnan. 2016. "Prospects for Precise Measurements with Echo Atom Interferometry" Atoms 4, no. 3: 19. https://doi.org/10.3390/atoms4030019

APA Style

Barrett, B., Carew, A., Beica, H. C., Vorozcovs, A., Pouliot, A., & Kumarakrishnan, A. (2016). Prospects for Precise Measurements with Echo Atom Interferometry. Atoms, 4(3), 19. https://doi.org/10.3390/atoms4030019

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