## Squeezed-light-enhanced atom interferometry below the standard quantum limit

Stuart S. Szigeti 1 2 , , ∗ , Behnam Tonekaboni , Wing Yung S. Lau , Samantha N. Hood , and Simon A. Haine 1 1 1 1

1 School of Mathematics and Physics, University of Queensland, Brisbane, QLD, 4072, Australia and

2 ARC Centre for Engineered Quantum Systems, The University of Queensland, Brisbane, QLD 4072, Australia ∗

Weinvestigate the prospect of enhancing the phase sensitivity of atom interferometers in the MachZehnder configuration with squeezed light. Ultimately, this enhancement is achieved by transferring the quantum state of squeezed light to one or more of the atomic input beams, thereby allowing operation below the standard quantum limit. We analyze in detail three specific schemes that utilize (1) single-mode squeezed optical vacuum (i.e. low frequency squeezing), (2) two-mode squeezed optical vacuum (i.e. high frequency squeezing) transferred to both atomic inputs, and (3) two-mode squeezed optical vacuum transferred to a single atomic input. Crucially, our analysis considers incomplete quantum state transfer (QST) between the optical and atomic modes, and the effects of depleting the initially-prepared atomic source. Unsurprisingly, incomplete QST degrades the sensitivity in all three schemes. We show that by measuring the transmitted photons and using information recycling [Phys. Rev. Lett. 110 , 053002 (2013)], the degrading effects of incomplete QST on the sensitivity can be substantially reduced. In particular, information recycling allows scheme (2) to operate at the Heisenberg limit irrespective of the QST efficiency, even when depletion is significant. Although we concentrate on Bose-condensed atomic systems, our scheme is equally applicable to ultracold thermal vapors.

PACS numbers: 03.75.Dg, 42.50.Gy, 42.50.Dv, 03.75.Be, 42.50.-p

## I. INTRODUCTION

Atom interferometry is a leading precision measurement technology, having demonstrated state-of-the-art measurements of accelerations and rotations [1-6], gravity gradients [7, 8], magnetic fields [9], the fine structure constant ( α ) [10, 11], and Newton's gravitational constant ( G ) [12-15]. Further increases to the sensitivity of atom interferometers would allow for some exciting science, such as improved tests of the weak equivalence principle [16-18], searches for quantum gravitational effects [19], and the measurement of gravitational waves [20, 21]. Current state-of-the-art atom interferometers utilize uncorrelated sources, which can operate no better than the standard quantum limit (SQL) - i.e. the sensitivity scales as 1 / √ N where N is the number of detected atoms. Unfortunately, current atomic sources have low fluxes (in comparison with photon sources), and there exist considerable technical barriers to increasing this flux [22]. Hence, developing atom interferometers that operate below the SQL, which have a better 'per atom sensitivity' than current devices, is of great interest.

schemes require large inter-atomic interactions, small atom number, and give little control over the motional atomic state. These are the opposite conditions required for precision atom interferometry. Alternatively, squeezed atomic states can be generated by mapping the quantum state of squeezed light to an atomic field [3843]. Given that squeezed light is known to give sub-SQL sensitivities in optical interferometers [44], transferring the entanglement to atomic degrees of freedom should similarly allow sub-SQL atom interferometry. Importantly, since the squeezing is generated independently of the atomic source, in principle this technique gives high flux (relative to state-of-the-art atomic sources), weaklyinteracting squeezed atomic states in targeted motional states - ideal for atom interferometry.

Sub-SQL atom interferometers necessarily exploit entanglement, and a number of proposals exist for generating the required entanglement between the atomic degrees of freedom. These are based on phenomena as diverse as molecular dissociation [23], spin-exchange collisions [24-26], atomic four-wave mixing [27-30], and atomic Kerr squeezing [31-37]. However, all these

∗ s.szigeti@uq.edu.au

In this paper, we present three squeezed-lightenhanced atom interferometry schemes, and show that these are all capable of sub-SQL phase sensitivities. We explicitly consider Bose-condensed atomic sources, as the narrow velocity distribution and large coherence length of a Bose-Einstein condensate (BEC) offers considerable advantages over thermal sources, including more precise manipulation of the motional state, increased visibility, and the prospect of feedback-stabilization under high flux outcoupling [45-51]. However, many of the results in this paper are equally applicable to ultracold thermal sources. Our analysis considers the effects of incomplete quantum state transfer (QST) between the optical and atomic modes, and the effect of depleting the initial condensate mode. We also incorporate the technique of information recycling [52] into our schemes, and demonstrate that this can be used to combat the negative effects of incomplete QST. Given the maturity of squeezed light generation

technology, and the high efficiency of photon detection, it would be relatively straightforward to incorporate our schemes into existing state-of-the-art atom interferometers. Consequently, squeezed-light enhancement and information recycling offer a promising path to improved sensitivity in atom-interferometer-based technologies.

The structure of our paper is as follows. In Sec. II we derive a simplified quantum model of atom-light coupling, based on two-photon Raman transitions, and show how this atom-light coupling can be used to achieve QST between the optical and atomic modes (Sec. II A), and coherent atomic beamsplitting (Sec. II B). Sec. III briefly reviews atom interferometry in the Mach-Zehnder configuration. Our first interferometry scheme, where enhancement is achieved with a single-mode squeezed optical vacuum (i.e. low frequency squeezing), is presented in Sec. IV. The sensitivity of this scheme is derived for complete QST and an undepleted initial atomic source (Sec. IVA), and compared to truncated Wigner simulations that include depletion. The effects of incomplete QST on this scheme, and how such effects can be reduced with information recycling, are considered in Sec. IV B. The effects of losses are briefly explored in Sec. IVC. In Sec. V, we present our second scheme, which utilizes a two-mode squeezed optical vacuum (i.e. high frequency squeezing) transferred to both atomic inputs. Once again, we quantitatively examine the phase sensitivity when QST is both complete (Sec. V A) and incomplete (Sec. V B), and numerically consider the effects of depletion (Sec. VC). Section VI presents and analyses our final atom interferometry scheme, where only a single atomic input is enhanced with a two-mode squeezed optical vacuum. These three schemes are then compared and summarized in Sec. VII.

## II. THEORETICAL MODEL FOR ATOM-LIGHT COUPLING

Our system has been previously described in [40-42]. We begin with a BEC consisting of atoms with two hyperfine states | 1 〉 and | 2 , 〉 separated in energy by an amount glyph[planckover2pi1] ω 2 . These two levels are coupled with a Raman transition via two optical fields, ˆ E 1 (probe beam, wavevector k 1 ) and ˆ E 2 (control beam, wavevector k 2 ), which are both detuned from an excited state | 3 〉 (see Fig. 1). We assume that the control field ˆ E 2 is much more intense than the probe field ˆ E 1 , allowing us to ignore depletion and quantum fluctuations and approximate the control field as a classical plane wave i.e. ˆ E 2 ( r , t ) ≈ E 2 exp[ ( i k 2 · r -ω t c )]. Furthermore, by design ∆ p ≡ ω 3 -ω p and ∆ c ≡ ω 3 -ω 2 -ω c are large compared with the Rabi frequencies of the | 1 〉 → | 〉 3 and | 2 〉 → | 〉 3 transitions. Therefore, the excited state can be adiabatically eliminated [53-55], giving an effective coupling between atomic states | 1 〉 and | 2 . 〉 Finally, we assume a dilute atomic sample, since it is generally optimal for atom interferometry to operate in a regime

FIG. 1. Energy level scheme for a three-level Raman transition comprising two non-degenerate hyperfine ground states, | 1 〉 and | 2 . 〉 The BEC is initially formed in the state | 1 , and 〉 population is transferred to | 2 〉 via the absorption of a photon from E 1 (the probe beam) and the emission of a photon into E 2 (the control beam). The probe and control beams are detuned from the excited state | 3 〉 by an amount ∆ p = ω 3 -ω p and ∆ c = ω 3 -ω 2 -ω c , respectively.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0067v3-figure-6d33650b-3af3-42cc-adef-1b903d142260.jpg)

where the inter-atomic interactions are negligible [45, 5658]. Under these approximations and the rotating-wave approximation [59], the Hamiltonian for the system becomes

$$\text{igh} \quad & \text{comes} \\ \text{is} \text{-} \\ \text{m-} \\ \text{of} \quad & \text{$\hat{\Psi}_{1}(\mathbf r)$} \hat{\Psi}_{2}(\mathbf r)$} \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{2}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{2}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbb { \Phi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{2}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbb { \Phi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r)\hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psis}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}\mathbf r \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Phi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_\text{$\hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Phis}_\text{$\hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathsf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \text{$\hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r) \hat{\Psi}_{1}(\mathbf r)$$

where ˆ H light is the Hamiltonian for the free photon field. Here, H 1 ( r ) = H ( r ) and H 2 ( r ) = H ( r ) + glyph[planckover2pi1] ω 2 , where H ( r ) is the single-atom Hamiltonian common to both hyperfine states. ˆ ( ψ 1 r ) and ˆ ( ψ 2 r ) are the usual bosonic field operators for atoms in hyperfine states | 1 〉 and | 2 , 〉 respectively, satisfying the commutation relation [ ˆ ( ψ i r ) , ψ ˆ ( † j r ′ )] = δ ij δ ( r -r ′ ). Similarly, ˆ E 1 ( r ) is the annihilation operator for the probe field satisfying [ ˆ E 1 ( r ) , E ˆ † 1 ( r ′ )] = δ ( r -r ′ ). We assume that the probe field ˆ E 1 ( r ) has a small spread of frequencies around ω p = c | k 1 | , in which case the effective coupling strength is given by

$$g = d _ { 1 2 } \sqrt { \frac { \omega _ { p } } { 2 \hbar { \epsilon } _ { 0 } } } \frac { \Omega } { \Delta _ { p } },$$

where d 12 is the dipole moment of the | 1 〉 → | 3 〉 transition, Ω is the Rabi frequency for the | 2 〉 → | 3 〉 transition, effected by the classical control field, and glyph[epsilon1] 0 is the permittivity of free space.

As is typical in atom interferometry, we assume that all atoms are initially in a single motional mode u 0 ( r ) of state | 1 . 〉 Furthermore, we will assume that our probe field is vacuum, except for occupation in a pulse characterized by a wave packet propagating in the z -direction:

$$\Gamma _ { 0 } ^ { \prime } \quad u _ { p } ( { \mathbf r }, t ) e ^ { i { \mathbf k } _ { 1 } \cdot { \mathbf r } } = u _ { \mathrm t r a n s } ( x, y ) u _ { \mathrm p r o p } ( z - c t ) e ^ { i { \mathbf k } _ { 1 } \cdot { \mathbf r } }, \quad ( 3 )$$

satisfying L | k 1 | glyph[greatermuch] 1, where L is the characteristic length scale of u prop ( z ). Provided the timescale for population transfer is fast compared with the timescale for atomic motion, the Hamiltonian (1) can be simplified considerably by expanding the field operators in the appropriate mode basis, and keeping only those modes that are highly occupied. Explicitly, we approximate

$$\hat { \psi } _ { 1 } ( { \mathbf r } ) \approx u _ { 0 } ( { \mathbf r } ) \hat { a } _ { 1 }, \quad \ \ \ ( 4 \mathbf a ) \quad \ \ \ n m r$$

$$\hat { \psi } _ { 2 } ( { \mathbf r } ) \approx u _ { 0 } ( { \mathbf r } ) e ^ { i ( { \mathbf k } _ { 1 } - { \mathbf k } _ { 2 } ) \cdot { \mathbf r } } \hat { a } _ { 2 }, \quad \quad ( 4 b ) \quad \text{ now}$$

$$\hat { E } _ { 1 } ( { \mathbf r } ) \approx u _ { p } ( { \mathbf r }, t ) e ^ { i { \mathbf k } _ { 1 } \cdot { \mathbf r } } \hat { b }. \quad \quad ( 4 c ) \quad \text{whi}$$

The simplified Hamiltonian is therefore ˆ H ≈ H ˆ 0 + ˆ H ′ int , where

$$\hat { \mathcal { H } } _ { 0 } = \left ( \hbar { \omega } _ { 2 } + \frac { \hbar { ^ { 2 } } } { 2 m } | { \mathbf k } _ { 1 } - { \mathbf k } _ { 2 } | ^ { 2 } \right ) \hat { a } _ { 2 } ^ { \dagger } \hat { a } _ { 2 } + \hbar { \omega } _ { p } \hat { b } ^ { \dagger } \hat { b } \quad ( 5 \mathbf a ) \quad \text{the d} \quad$$

$$\hat { \mathcal { H } } ^ { \prime } _ { \text{int} } = \hbar { g } f ( t ) \left ( \hat { a } _ { 1 } \hat { a } _ { 2 } ^ { \dagger } \hat { b } \, e ^ { i \omega _ { c } t } + h. c. \right ), \quad \quad ( 5 b ) \quad \hat { a } _ { 2 } ($$

and

$$f ( t ) = \int d { \mathbf r } \, | u _ { 0 } ( { \mathbf r } ) | ^ { 2 } u _ { p } ( { \mathbf r }, t ). \quad \quad ( 6 ) \quad \text{when}$$

Note that we have neglected the single-atom energy contribution due to H j ( r ), which is approximately a constant energy offset on the timescale of the population transfer. The time dependence in ˆ is due to the propH agation of the probe wave packet u p ( r , t ). Interestingly, an identical Hamiltonian is obtained if the probe beam is continuous-wave and the classical control field is shaped by the temporal function f ( ). t This would couple the mode of the probe field defined by the mode function u p ( r , t ) to the condensate.

Finally, after moving to the interaction picture ˆ H → ˆ H int = ˆ U † ˆ H ′ int ˆ U with ˆ U = exp( i ˆ H 0 t/ glyph[planckover2pi1] ), we obtain the following Heisenberg equations of motion:

$$i \dot { \hat { a } } _ { 1 } = g f ( t ) \hat { a } _ { 2 } \hat { b } ^ { \dagger } e ^ { i \delta t } \, \quad \, \quad \, ( 7 \tt a ) \quad \, \ s c r$$

$$i \dot { \hat { a } } _ { 2 } = g f ( t ) \hat { a } _ { 1 } \hat { b } e ^ { - i \delta t } \, \quad \ \ ( 7 b ) \quad \text{th}$$

$$i \dot { \hat { b } } = g f ( t ) \hat { a } _ { 2 } \hat { a } _ { 1 } ^ { \dagger } e ^ { i \delta t } \,, \quad \quad \quad ( 7 c )$$

where δ = ω p -ω c -ω 2 -glyph[planckover2pi1] | k 1 -k 2 | 2 / (2 m ) is the two-photon detuning, which is an experimental parameter freely adjustable by tuning the frequency offset between E 1 and E 2 . We only consider the optimal case δ = 0, when the system is on-resonance such that the energy transferred by the two-photon transition perfectly matches the change in electronic and kinetic energies of the atom. Equations (7) form the basis for the rest of this paper, as they allow us to describe the process of QST (of the light state to the atomic state) and our conventional coherent atomic beamsplitters. Each of these processes is discussed in more detail below.

## A. Quantum state transfer (QST)

The goal of this process is to take a pulse of light, described by annihilation operator ˆ , b and coherently map its quantum state onto the atomic mode ˆ . a 2 Thus, if ˆ b was initially in some interesting quantum state, such as a squeezed state, or was entangled with another mode, after QST ˆ a 2 will also be in this state and/or be entangled with this other mode. To achieve perfect QST, the number of photons in mode ˆ must be much less than the b initial number of condensate atoms in mode ˆ . a 1 Then, for a sufficiently short atom-light coupling time, a small number of atoms can be transferred to mode ˆ , which a 2 now have the initial quantum state of the light in mode ˆ , b while leaving the number of atoms in mode ˆ a 1 essentially unchanged. We can therefore make the undepleted reservoir approximation ˆ a 1 → √ N a 1 , where N a 1 is the mean number of atoms in mode ˆ . Under this approximation, a 1 the dynamics of QST are described wholly by Eqs (7b) and (7c), which can be solved exactly:

$$\hat { a } _ { 2 } ( t ) = \hat { a } _ { 2 } ( t _ { 0 } ) \cos ( \theta _ { Q S T } / 2 ) - i \hat { b } ( t _ { 0 } ) \sin ( \theta _ { Q S T } / 2 ), \ \ ( 8 a )$$

$$\hat { b } ( t ) = \hat { b } ( t _ { 0 } ) \cos ( \theta _ { Q S T } / 2 ) - i \hat { a } _ { 2 } ( t _ { 0 } ) \sin ( \theta _ { Q S T } / 2 ) \,, \ \ ( 8 b )$$

where θ QST ( ) t ≡ 2 g √ N a 1 ∫ t 0 dt ′ f ( t ′ ). As was shown in Jing et al. [38], when θ QST = π we have complete QST , such that the quantum state of ˆ ( b t 0 ) is perfectly mapped to the quantum state of ˆ a 2 ( t 1 ), up to a phase factor. However, achieving the required coupling strength for complete QST is likely to be challenging in practice. Thus incomplete QST , with 0 &lt; θ QST &lt; π , is the likely experimental regime, and consequently is considered in detail throughout this paper. Of course, incomplete QST also occurs when the undepleted reservoir approximation breaks down, which occurs when the initial number of photons in ˆ becomes comparable to the number of atoms b in ˆ a 1 . In this regime the solution (8) is invalid, and we require a numeric solution to Eqs (7). Moreover, as θ QST is no longer well-defined, the QST process must be described by an alternative metric. The natural choice is the QST efficiency :

$$\mathcal { Q } ( t ) \equiv \frac { \langle \hat { a } _ { 2 } ^ { \dagger } ( t ) \hat { a } _ { 2 } ( t ) \rangle } { \langle \hat { b } ^ { \dagger } ( t _ { 0 } ) \hat { b } ( t _ { 0 } ) \rangle },$$

which is a measure of the percentage of atoms outcoupled compared with the total number of input photons. Although this is a somewhat cruder metric than the fidelity, it is operationally more convenient, and certainly more than adequate for our purposes. When ˆ a 1 is treated as an undepleted reservoir, and Eqs (8) apply, then Q = sin ( 2 θ QST / 2).

It is instructive to conceptualize the process of QST as an atom-light beamsplitter . That is, a type of beamsplitter that distributes an initial quantum state of light ˆ ( b t 0 ) amongst an atomic mode ˆ ( a 2 t 1 ) and an outgoing mode of light ˆ ( b t 1 ), in much the same way as a conventional beamsplitter distributes a quantum state of light amongst two outgoing modes of light. This helpful analogy is illustrated diagrammatically in Fig. 2. As a concrete example, note that when θ QST = π/ 2, we have a

FIG. 2. Analogy between atom-light QST process, and a beamsplitter. A photon from mode ˆ can either implement b a Raman transition, resulting in an atom being outcoupled into mode ˆ , or the photon can be transmitted, remaining in a 2 mode ˆ . b In the regime where the number of photons in mode ˆ is b much less than the number of atoms in the condensate ˆ , we can treat ˆ a 1 a 1 as an undepletable reservoir.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0067v3-figure-6c98695a-f893-4857-91c1-9077bc4eb215.jpg)

50 50 atom-light beamsplitter, corresponding to a QST / efficiency of 50%. We invoke the atom-light beamsplitter analogy throughout this paper, as it allows quite complicated atom-interferometric schemes to be conceptualized as simple linear-optical setups.

## B. Coherent atomic beamsplitters

The coherent beamsplitting and reflection of two atomic modes via two-photon Raman transitions is a mature experimental technique that has been used with much success in atom interferometry [60]. In these experiments, both light pulses have a mean photon number much larger than the number of atoms in modes ˆ a 1 and ˆ , and are therefore well-approximated as undepletable a 2 coherent states. In this regime, the atom-light coupling is a conventional Raman transition, in the sense that the coupling coherently transfers population between the two atomic modes without significantly affecting the state of the optical modes. This can be explicitly seen by making the replacement ˆ b → √ N b , for mean photon number N b , and subsequently solving Eqs (7a) and (7b), yielding

$$\hat { a } _ { 1 } ( t ) = \hat { a } _ { 1 } ( t _ { 0 } ) \cos ( \theta _ { \text{BS} } / 2 ) - i \hat { a } _ { 2 } ( t _ { 0 } ) \sin ( \theta _ { \text{BS} } / 2 ), \ \ ( 1 0 \text{a} )$$

$$\hat { a } _ { 2 } ( t ) = \hat { a } _ { 2 } ( t _ { 0 } ) \cos ( \theta _ { \text{BS} } / 2 ) - i \hat { a } _ { 1 } ( t _ { 0 } ) \sin ( \theta _ { \text{BS} } / 2 ), \ \ ( 1 0 \text{b} )$$

with θ BS ( ) t ≡ 2 g √ N b ∫ t 0 dt ′ f ( t ′ ). When θ BS = π/ 2 and θ BS = π , we have a 50 / 50 atomic beamsplitter and mirror, respectively.

It is worth noting that Eqs (10) apply to other coherent beamsplitting techniques, such as Bragg pulses [61, 62] and Bloch oscillations [63]. Indeed, given that inertial sensors based on Bose-condensed sources and large momentum transfer beamsplitters offer a promising alternative route to improved sensitivity [46, 64, 65], incorporating such interferometers within the squeezed-light enhanced schemes outlined below is a most attractive prospect.

FIG. 3. Analogy between a conventional two-photon Raman transition and a beamsplitter. An atom in mode ˆ a 1 can absorb a photon from mode ˆ , and emit it into b E 2 via stimulated emission, producing one atom in mode ˆ a 2 . An atom can be transferred from ˆ a 2 to ˆ a 1 via the reverse process. In the regime where the number of photons in both ˆ and b E 2 is much greater than the number of photons in mode ˆ a 1 , we can treat the optical modes as undepletable reservoirs, and the process behaves as an atomic beamsplitter.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0067v3-figure-db247d3f-f3d2-41e5-ba5f-6e34c4cb3ecd.jpg)

## III. REVIEW: THE MACH-ZEHNDER ATOM INTERFEROMETER

A standard atom interferometer in the Mach-Zehnder (MZ) configuration operates by first splitting a singlemode population of atoms (in mode ˆ a 1 , say) with a 50 / 50 beamsplitter, and then letting a phase difference between the two modes accumulate due to the physical process of interest (e.g. linear acceleration). The two modes are redirected together via an atomic mirror, and subsequently mixed with a second 50 / 50 beamsplitter. This converts the phase shift φ into a population difference between the two modes, ˆ N a 1 -ˆ N a 2 , where ˆ N a i = ˆ ˆ . a a † i i The population difference at the output can be measured, allowing for an accurate estimate of the phase shift and hence the relevant physical quantity of interest. This output signal can be expressed in terms of the two input modes (i.e. the atomic modes prior to the first beamsplitter) by repeated application of Eqs (10):

$$\hat { S } _ { a } & = \hat { N } _ { a _ { 1 } } ( t _ { f } ) - \hat { N } _ { a _ { 2 } } ( t _ { f } ) \\ & = \cos \phi \left ( \hat { N } _ { a _ { 2 } } ( t _ { 1 } ) - \hat { N } _ { a _ { 1 } } ( t _ { 1 } ) \right ) \\ & + \sin \phi \left ( \hat { a } _ { 1 } ^ { \dagger } ( t _ { 1 } ) \hat { a } _ { 2 } ( t _ { 1 } ) + \hat { a } _ { 2 } ^ { \dagger } ( t _ { 1 } ) \hat { a } _ { 1 } ( t _ { 1 } ) \right ) \,, \quad ( 1 1 )$$

where t 1 and t f are the times immediately before the first beamsplitter and immediately after the second beamsplitter, respectively. The phase sensitivity of the interferometer can be determined from this output signal via the expression [66]

$$\Delta \phi = \sqrt { \frac { V ( S _ { a } ) } { ( d \langle \hat { S } _ { a } \rangle / d \phi ) ^ { 2 } } } \,, \quad \quad ( 1 2 )$$

where V ( Q ) ≡ 〈 ˆ Q 2 〉 - 〈 ˆ Q 〉 2 is the variance. If one input is either a traditionally prepared Bose-condensed source

(modelled as a coherent state or Fock state) or a lasercooled thermal source, and the other input is vacuum, then the standard MZ interferometer can achieve a sensitivity no better than the SQL ∆ φ = 1 / √ N t , where N t is the total number of atoms measured at the output. SubSQL sensitivities require a more exotic initial state [67]. As shown below, QST is a neat and practical method of generating just such a state.

## IV. ENHANCED ATOM INTERFEROMETRY WITH SINGLE-MODE SQUEEZED LIGHT

We first consider using a single-mode squeezed optical vacuum to enhance the sensitivity of atom interferometry. Although generating the single-mode squeezed vacuum state considered here is feasible, it is likely to be a technically challenging procedure. Ultimately, the difficulty stems from the frequencies of the light field where squeezing can be observed. Conceptually, it is impossible to squeeze only a single frequency of light; naturally squeezing occurs across a range of frequencies. More precisely, an optically squeezed state has quantum correlations between sidebands symmetrically distributed around a central carrier frequency, ω p = c | k 2 | [see Fig. 11(a)]. Below some critical frequency ω crit , technical considerations usually ensure that these correlations are masked by uncorrelated classical noise [the red frequencies in Fig. 11(a)]. Hence, in order for the optical mode taking part in the QST process to display quantum correlations, we require ∆ ω glyph[greatermuch] 2 ω crit , where ∆ ω is the characteristic width of F ω ( ), the Fourier transform of f ( ). t Although this is technically challenging, there has recently been demonstrations of significant squeezing at frequencies below 100 Hz [68]. We consider the application of higher-frequency squeezed light sources, which are easier to generate, in later sections of this paper.

To begin, suppose that the squeezed light is generated via an optical parametric oscillator (OPO) [69]. Without loss of generality, we assume that k 1 points in the z -direction, i.e. ˆ ( E r , t ) ≈ u trans ( x, y ) ˆ ( E z, t ). Then, in the momentum basis, the photon fields that form the inputs and outputs of the OPO, at times t i and t 0 respectively, are related by:

$$\hat { \phi } ( q, t _ { 0 } ) = \hat { \phi } ( q, t _ { i } ) \cosh r - i e ^ { i \theta _ { s q } } \hat { \phi } ^ { \dagger } ( - q, t _ { i } ) \sinh r \,, \ \ ( 1 3 ) \quad \text{the} \, 1$$

where q = k -k 1 , and r and θ sq are the squeezing parameter and angle, respectively. Strictly, ˆ ( φ q, t ) and ˆ ( φ † q, t ) are defined in terms of the position-space photon field ˆ ( E z, t ):

$$\hat { \phi } ( q, t ) = \frac { 1 } { \frac { } { \sqrt { 2 \pi } } } \int _ { \substack { d z \, e ^ { - i ( q + k _ { 1 } ) z } \hat { E } ( z, t ) \quad ( 1 4 \text{a} ) } }$$

$$\hat { \phi } ^ { \dagger } ( q, t ) = \frac { \cdot } { \sqrt { 2 \pi } } \int d z \, e ^ { i ( q + k _ { 1 } ) z } \hat { E } ^ { \dagger } ( z, t ). \quad ( 1 4 b ) \quad \text{when}$$

Only the portion of the photon field under the pulse envelope u p ( r , t ) interacts with the atoms at a given time t . Furthermore, the physics of the atom-light interaction is determined by the temporal area of the probe pulse. Consequently, the relevant mode of the photon field is

$$\dashrightarrow \lambda \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \lambda \\ \hat { b } ( t ) & \equiv \int \, d { \mathbf r } \, u _ { p } ^ { * } ( { \mathbf r }, t ) e ^ { - i { \mathbf k } _ { 1 } \cdot { \mathbf r } } \hat { E } ( { \mathbf r }, t ) \\ & = \int d q \, e ^ { i q c t } U _ { \text{prop} } ^ { * } ( q ) \hat { \phi } ( q, t ) \\ & \approx \int d q \, U _ { \text{prop} } ^ { * } ( q ) \hat { \phi } ( q, t ), & ( 1 5 )$$

where

$$U _ { \text{prop} } ( q ) = \frac { 1 } { \sqrt { 2 \pi } } \int d z \, e ^ { - i q z } u _ { \text{prop} } ( z ). \quad ( 1 6 )$$

The second line of Eq. (15) follows from Eq. (14a) and ∫ dxdy u | trans ( x, y ) | 2 = 1, and the dominant contribution to the integral occurs close to the carrier frequency k 1 (i.e. q = 0), thereby justifying the approximate expression in the third line [70]. It then follows from Eq. (13) that

$$\hat { b } ( t _ { 0 } ) & = \left [ \int d q \, U _ { \text{prop} } ^ { * } ( q ) \hat { \phi } ( q, t _ { i } ) \right ] \cosh r \\ & \quad - i e ^ { i \theta _ { \text{sq} } } \left [ \int d q \, U _ { \text{prop} } ^ { * } ( q ) \hat { \phi } ^ { \dagger } ( - q, t _ { i } ) \right ] \sinh r. \quad ( 1 7 )$$

The first term in square brackets is clearly ˆ ( b t i ). The second term in square brackets equals [ ˆ ( b t i )] † = ˆ ( b † t i ) provided U prop ( k ) is real and symmetric, which in practice is easy to satisfy. Consequently, the photon mode output from the OPO is simply

$$\hat { b } ( t _ { 0 } ) = \hat { b } ( t _ { i } ) \cosh r - i e ^ { i \theta _ { s q } } \hat { b } ^ { \dagger } ( t _ { i } ) \sinh r, \quad ( 1 8 )$$

which for initial vacuum input is a single-mode squeezed state. Eq. (17) and Eq. (18) illustrate the relationship between the squeezing spectra, which is typically measured in optical squeezing experiments, and the temporal modes relevant for quantum state transfer.

Our scheme is summarized in Fig. 4. An initial squeezed vacuum state ˆ ( b t 0 ) is used to transfer a small number of atoms from mode ˆ a 1 to ˆ a 2 via the QST process (8), thereby transferring some or all of the quantum state from ˆ ( b t 0 ) to ˆ a 2 ( t 1 ). The modes ˆ ( a 1 t 1 ) and ˆ ( a 2 t 1 ) then form the two input modes for a MZ atom interferometer (i.e. are coherently split, reflected and recombined), yielding the two outputs ˆ a 1 ( t f ) and ˆ a 2 ( t f ), used to construct the difference signal ˆ S a [see Eq. (11)]. Expectations are calculated with respect to the initial state | Ψ(0) , defined such that 〉

$$\hat { \phi } ( q, t _ { i } ) | \Psi ( 0 ) \rangle = \hat { a } _ { 2 } ( t _ { 0 } ) | \Psi ( 0 ) \rangle = 0 \quad \ ( 1 9 \text{a} )$$

$$\hat { a } _ { 1 } ( t _ { 1 } ) | \Psi ( 0 ) \rangle = \sqrt { N _ { a _ { 1 } } ( t _ { 1 } ) } | \Psi ( 0 ) \rangle \,, \quad ( 1 9 b )$$

where conservation of total atom number N t implies that N a 1 ( t 1 ) = N t -〈 ˆ N a 2 ( t 1 ) 〉 . Of course, in writing Eq. (19b) we have assumed that the condensate is initially in a

FIG. 4. A scheme for enhancing an atom interferometer with a single-mode squeezed optical vacuum. The squeezed light is used to outcouple a small number of atoms from a BEC. These outcoupled atoms and the remaining condensate atoms form the two inputs to a MZ atom interferometer.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0067v3-figure-f7aa2f58-a8d9-4930-a416-a2a5c2c97b77.jpg)

coherent state, and remains in a coherent state during the QST process. As shown below, this assumption is only valid when the number of outcoupled atoms 〈 ˆ N a 2 ( t 1 ) 〉 is much less than N t .

## A. Complete quantum state transfer

We first consider the optimal regime of complete QST θ QST = π (i.e. Q = 1). In this case Eq. (8a) implies that ˆ a 2 ( t 1 ) = -ib ˆ ( t 0 ). The expectation of the difference signal (11) simplifies to

$$\langle \hat { S } _ { a } \rangle = \left ( N _ { t } - 2 \sinh ^ { 2 } r \right ) \cos \phi \,. \quad \ \ ( 2 0 ) \quad \hat { \mathfrak { n r } }$$

In order to achieve minimum phase sensitivity, the variance in the signal must attain a minimum when the slope of the output signal is a maximum [see Eq. (12)]. This occurs at phase φ = π/ 2 and squeezing angle θ sq = π/ 2. Figure 5 shows the signal and phase sensitivity as a function of φ for θ sq = π/ 2. For r = 3 8, we achieve an en-. hancement in sensitivity of approximately 30 times better than the SQL. For comparison, we have shown the case where the two inputs to the MZ interferometer are a coherent state and a vacuum state, respectively.

For these optimal values, the variance is simply

$$V ( S _ { a } ) = N _ { t } e ^ { - 2 r } + 2 e ^ { - r } \sinh ^ { 3 } r. \quad \ \ ( 2 1 ) \quad \text{ one}$$

The minimum phase sensitivity, as a function of r and

FIG. 5. (Top). The expectation value (black, solid line) and quantum uncertainty (light blue, shading) for the signal ˆ normalized by total number of atoms ( S N t = 10 ) for 6 θ sq = π/ 2, and r = r opt ≈ 3 8. . The inset shows that the variance is less than the SQL near φ = π/ 2. (Bottom) The phase sensitivity for complete QST (blue, solid line), compared to the SQL.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0067v3-figure-03d1aed0-03af-430b-8ed6-a4e54c250425.jpg)

## N t , is therefore

$$\Delta \phi _ { \min } & = \frac { \sqrt { N _ { t } e ^ { - 2 r } + 2 e ^ { - r } \sinh ^ { 3 } r } } { N _ { t } - 2 \sinh ^ { 2 } r } \quad \ ( 2 2 ) \\ & \approx \frac { e ^ { - r } } { \sqrt { N _ { t } } },$$

where the approximate expression in the second line is only true in the limit 〈 ˆ ( N t b 0 ) 〉 = sinh 2 r glyph[lessmuch] N t , where ˆ N b = ˆ ˆ . b b †

Figure 6 shows the minimum interferometer sensitivity [Eq. (22)] as a function of the squeezing parameter r , which determines the average number of input photons via 〈 ˆ ( N t b 0 ) 〉 = sinh 2 r , for a range of initial BEC atom numbers. When N t glyph[greatermuch] 1, our analytic model predicts that an optimal squeezing parameter of r opt ≈ ln(4 N / t ) 4 yields a minimum sensitivity of ∆ φ min ≈ 1 /N 3 / 4 t . This is significantly less than the SQL, and furthermore is the best sensitivity possible in this undepleted regime provided sinh 2 r glyph[lessmuch] N t [71, 72]. For N t = 10 , this gives 6 an enhancement in sensitivity of approximately 32 compared with the SQL, which is equivalent to increasing the total number of atoms by a factor of 10 3 at the SQL. For this value of N t , the number of atoms outcoupled at r = r opt is sinh 2 r opt ≈ 500, suggesting that the undepleted reservoir model is still reasonably valid in this

FIG. 6. (Top) Minimum phase sensitivity (i.e. at φ = θ sq = π/ 2) as a function of squeezing parameter r for initial atom numbers N t = 10 4 (blue, top), N t = 10 5 (green, middle) and N t = 10 6 (magenta, bottom). The solid curves are the analytic solution [Eq. (22)], while the points correspond to a TW numerical solution. The standard error in the TW solutions is no larger than the point width. There is good agreement between the analytics and numerics when 〈 ˆ ( N t b 0 ) 〉 glyph[lessmuch] N t , where the initial condensate is not significantly depleted. (Bottom) The maximum QST efficiency Q max = max t Q ( ) t as a function of r . Theoretically, complete QST ( Q = 1) is achievable provided less than ∼ 10% of the total condensate number is outcoupled. The analytics predict Q = 1 always, so any deviation from this is due to depletion from the condensate mode ˆ . Note, however, that there exist regimes where a 1 mode ˆ a 1 must be treated quantum mechanically even though Q max = 1.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0067v3-figure-f5b83c6d-0a6c-472c-8539-a3da03b6e48f.jpg)

regime.

To include the effects of depletion from the condensate, we need to treat mode ˆ a 1 quantum mechanically, and simulate the full quantum dynamics of the QST process, which are governed by Eqs (7). This can be done via the truncated Wigner (TW) phase space method [7376]. Following standard methods [69, 77], the Heisenberg equations of motion are converted into a partial differential equation (PDE) for the Wigner quasi-probability distribution. Once third and higher-order derivatives are truncated (an uncontrolled approximation, but one that is typically valid provided the occupation per mode is not too small for appreciable time periods [78, 79]), this PDE takes the form of a Fokker-Planck equation, which can be efficiently simulated by a set of stochastic differential equations (SDEs) for complex numbers α j ( ). t In our case, the set of SDEs corresponding to Eqs (7) (with

δ = 0) is

$$i \dot { \alpha _ { 1 } } = g f ( t ) \alpha _ { 2 } \beta ^ { * }, \text{ \quad \ \ } ( 2 3 \text{a} )$$

$$i \dot { \alpha _ { 2 } } = g f ( t ) \alpha _ { 1 } \beta, \text{ \quad \ \ } ( 2 3 b )$$

$$i \dot { \beta } = g f ( t ) \alpha _ { 2 } \alpha _ { 1 } ^ { * }, \text{ \quad \ \ } ( 2 3 c )$$

where we have made the correspondences ˆ ( ) a i t → α i ( ) t and ˆ ( ) b t → β t ( ). The initial conditions for these SDEs are randomly sampled from the Wigner distribution corresponding to the initial quantum state [80]. Specifically, just before the QST process, ˆ a 1 is in a coherent state of mean number N t , ˆ a 2 is in a vacuum state, and ˆ is in b a single-mode squeezed vacuum state. The initial conditions corresponding to these initial states are

$$\alpha _ { 1 } ( t _ { 0 } ) = \sqrt { N _ { t } } + \eta _ { \alpha _ { 1 } },$$

$$\alpha _ { 2 } ( t _ { 0 } ) = \eta _ { \alpha _ { 2 } },$$

$$\beta ( t _ { 0 } ) = \eta _ { \beta } \cosh r - i e ^ { i \theta _ { s q } } \eta _ { \beta } ^ { * } \sinh r. \quad ( 2 4 c )$$

The η i are complex, independent Gaussian noises satisfying η i = 0 and η η ∗ i j = δ ij . The expectation value of some arbitrary operator function h is then computed by averaging over solutions to Eqs (23) with initial conditions (24):

$$\ a n a ^ { \dots \dots } _ { \text{TW} } \quad & \langle \{ h ( \hat { a } _ { 1 } ^ { \dagger }, \hat { a } _ { 2 } ^ { \dagger }, \hat { b } ^ { \dagger }, \hat { a } _ { 1 }, \hat { a } _ { 2 }, \hat { b } ) \} _ { s y m } \rangle = \overline { h \left ( \alpha _ { 1 } ^ { * }, \alpha _ { 2 } ^ { * }, \beta ^ { * }, \alpha _ { 1 }, \alpha _ { 2 }, \beta \right ) } \,, \\ \text{ms is } \quad & \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \.$$

where 'sym' denotes symmetric ordering [69], and the overline denotes the average of simulated trajectories.

Since beamsplitting and mirror operations are implemented by strong coherent optical fields, we can approximate them as linear, and the complex amplitudes α 1 ( t f ) and α 2 ( t f ) can be directly evolved from α 1 ( t 1 ) and α 2 ( t 1 ) by repeated application of Eqs (10). The mean 〈 ˆ S a 〉 and variance V ( S a ) can then be computed using the relations

$$\langle \hat { N } _ { i } \rangle = \overline { | \alpha _ { i } | ^ { 2 } } - 1 / 2,$$

$$\langle \hat { N } _ { i } ^ { 2 } \rangle = \overline { | \alpha _ { i } | ^ { 4 } } - \overline { | \alpha _ { i } | ^ { 2 } },$$

$$\langle \hat { N } _ { 1 } \hat { N } _ { 2 } \rangle = \overline { \left ( | \alpha _ { 1 } | ^ { 2 } - 1 / 2 \right ) \left ( | \alpha _ { 2 } | ^ { 2 } - 1 / 2 \right ) }. \quad \ \ ( 2 6 c )$$

The effects of depletion on the minimum interferometer sensitivity, as numerically modelled by TW simulations, is shown in Fig. 6. Since the efficiency of the QST process does not uniquely depend upon the dynamics of the QST process (i.e. Q is not uniquely defined by the choice of f ( )), for simplicity simulations were performed t with a uniform f ( ). t For small r we find good agreement between the undepleted reservoir approximation and the full quantum dynamics. However, as r increases, the full quantum simulations actually predict a better sensitivity than the undepleted reservoir model, reaching a minimum at an r = r TW &gt; r opt . We can understand this feature by looking at the quadratures of the modes, for example ˆ X θ a 1 = exp( iθ )ˆ a 1 + exp( -iθ )ˆ a † 1 . As shown in Fig. 7, the discrepancy is due to changes in ˆ a 1 under

FIG. 7. (Top) TW simulations for the variance of atomic quadratures ˆ X θ c = exp( iθ )ˆ + exp( c -iθ )ˆ , c † where ˆ is an arc bitrary mode, for an initial condensate number of N t = 10 . 6 The quadratures for the initial squeezed optical input ˆ b 0 are analytic, and can easily be computed from Eq. (18). According to the undepleted reservoir analytics, which assumes ˆ a 1 is always coherent, the variance of any quadrature ˆ X θ a 1 ( t 1 ) is unity, and ˆ X 0 a 2 ( t 1 ) = ˆ X π/ 2 b t ( 0 ) and ˆ X π/ 2 a 2 ( t 1 ) = ˆ X 0 b t ( 0 ) when QST is complete. However, even though Q = 1 for all values of r shown here, TW simulations show that ˆ X 0 a 1 ( t 1 ) and ˆ X 0 a 2 ( t 1 ) diverge from the undepleted reservoir solution once r glyph[greaterorsimilar] 3 11. . (Bottom) Plot showing that the quantum noise (i.e. variance) on the signal ˆ S a also diverges from the undepleted reservoir model for r glyph[greaterorsimilar] 3 11. .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0067v3-figure-ca5ac2e5-0d0d-4486-a2dd-901399f78937.jpg)

depletion, evidenced by V ( ˆ X 0 a 1 ( t 1 ) ) and V ( ˆ X 0 a 2 ( t 1 ) ) deviating from the analytic solution. Although the state of ˆ a 1 is coherent under the undepleted reservoir approximation, quantum depletion creates a state with decreased variance in ˆ J x = [ˆ ( a 1 t 1 )ˆ a † 2 ( t 1 ) + ˆ a 2 ( t 1 )ˆ a † 1 ( t 1 )] / 2 (and increased variance in ˆ X 0 a 1 ( t 1 ) and ˆ X 0 a 2 ( t 1 ) ). This gives a reduction in the sensitivity at φ = π/ 2, since the noise in the signal is directly proportional to ˆ J x by virtue of Eq. (11).

For further increases in r , the effects of depletion become significant, and complete QST (i.e. Q = 1) is no longer possible. The maximum possible QST efficiency, as a function of r , is shown in Fig. 6. Unsurprisingly, this is contrary to the undepleted reservoir model, where Q = 1 always for θ QST = π .

## B. Incomplete quantum state transfer and information recycling

In practice, it may be difficult to achieve the required coupling strength g for complete QST. When θ QST &lt; π , and optimal values φ = π/ 2 and θ sq = π/ 2 are chosen, the undepleted reservoir model predicts a signal slope and variance of

$$\frac { d \langle \hat { S } _ { a } \rangle } { d \phi } = - \left ( N _ { t } - 2 \sin ^ { 2 } \left ( \theta _ { Q S T } / 2 \right ) \sinh ^ { 2 } r \right ), \quad ( 2 7 \mathbf a )$$

and

$$V ( S _ { a } ) & = N _ { t } e ^ { - r } \left ( \cosh r + \cos \theta _ { Q S T } \sinh r \right ) \\ & + 2 e ^ { - r } \sin ^ { 4 } \left ( \theta _ { Q S T } / 2 \right ) \sinh ^ { 3 } r, \text{ \quad \ \ } \text{(27b)}$$

respectively. This gives a minimum phase sensitivity of

$$\overset {. } {. } \underset { 6 } { \Delta } \phi _ { \min } = \frac { \sqrt { N _ { t } \left ( 1 - 2 Q e ^ { - r } \sinh r \right ) + 2 \mathcal { Q } ^ { 2 } e ^ { - r } \sinh ^ { 3 } r } } { N _ { t } - 2 \mathcal { Q } \sinh ^ { 2 } r }.$$

(28)

The blue curve in Fig. 8 shows ∆ φ min vs. Q at optimum squeezing parameter ( r = r opt when the undepleted reservoir approximation holds) for an initial condensate of N = 10 6 atoms. Although sub-SQL sensitivities are still possible for incomplete QST, the enhancement sharply reduces from the optimal 1 /N 3 / 4 t scaling as Q decreases from unity. For instance, at Q = 0 5, correspond-. ing to θ QST = π/ 2 in the undepleted reservoir regime, the enhancement to the sensitivity beyond the SQL has dropped by a factor of 20 to ∼ √ 2. The blue curve in Fig. 8 shows ∆ φ min vs. Q at optimum squeezing parameter ( r = r opt when the undepleted reservoir approximation holds) for an initial condensate of N = 10 6 atoms. Although sub-SQL sensitivities are still possible for incomplete QST, the enhancement sharply reduces from the optimal 1 /N 3 / 4 t scaling as Q decreases from unity. For instance, at Q = 0 5, corresponding to . θ QST = π/ 2 in the undepleted reservoir regime, the enhancement to the sensitivity beyond the SQL has dropped by a factor of 20 to ∼ √ 2.

Fortunately, this degradation to the sensitivity due to incomplete QST can be ameliorated with the technique of information recycling [52]. Specifically, after the atomlight beamsplitter, a quadrature of the transmitted field ˆ ( b t 1 ) can be measured via homodyne detection by mixing the field ˆ ( b t 1 ) with a bright local oscillator ˆ b LO ( t 1 ), which is assumed to be a large amplitude coherent state (see Fig. 9). In order to pick out the mode corresponding to ˆ , b the local oscillator would need to be temporally shaped such that it is mode-matched to u p ( r , t ). The noise on this homodyne signal is correlated with the noise on the atomic signal ˆ , S a and hence can be combined with the atomic signal to reduce the overall noise of the phase measurement. That is, from the apparatus depicted in Fig. 9 we construct the signal

$$\hat { S } = \hat { S } _ { a } - \mathcal { G } \hat { S } _ { b }, \text{ \quad \ \ } \text{(2)}$$

FIG. 8. Plots showing the QST efficiency dependence of the minimum phase sensitivity ∆ φ min for an atom interferometer enhanced by single-mode squeezing, and initial atom number N t = 10 . 6 For r = r opt ≈ 3 8, . the sensitivity of the purely atomic signal ˆ S a sharply degrades with decreasing Q , and the TW simulations (blue diamonds) agree with the analytic undepleted reservoir prediction (solid blue line) except near Q = 1. In contrast, the information-recycled signal ˆ = ˆ S S a -G ˆ S b displays considerably better sensitivities and a slower degradation as Q decreases. For r = r opt , the TW simulations (red circles) predict a slightly better sensitivity for Q glyph[greaterorsimilar] 0 4 than . the analytic undepleted reservoir prediction (red dashed line). For comparison, the upper and lower horizontal black dotted lines show the SQL and the theoretical limit reached by perfect QST (i.e. ∆ φ = 1 /N 3 / 4 t ), respectively. Interestingly, the TW simulations (red crosses) demonstrate that better sensitivities can be achieved for r = r TW ≈ 4 8 . &gt; r opt ; for ˆ S a this is only true near Q = 1 (blue squares).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0067v3-figure-34582c7b-7290-4e61-a2ce-592469b73b9b.jpg)

where

$$\hat { S } _ { b } = \hat { N } _ { b _ { L O } } ( t _ { f } ) - \hat { N } _ { b } ( t _ { f } ), \text{ \quad \ \ } ( 3 0 \text{a} )$$

$$\hat { b } ( t _ { f } ) = \frac { 1 } { \sqrt { 2 } } \left ( \hat { b } ( t _ { 1 } ) - i \hat { b } _ { \text{LO} } ( t _ { 1 } ) \right ), \quad \ \ ( 3 0 b ) \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

$$\hat { b } _ { \text{LO} } ( t _ { f } ) = \frac { \cdot } { \sqrt { 2 } } \left ( \hat { b } _ { \text{LO} } ( t _ { 1 } ) - i \hat { b } ( t _ { 1 } ) \right ), \quad ( 3 0 c ) \quad \text{$\substack{\partial \mu$\mu$}$}$$

and G is an adjustable gain parameter .

In order to see more clearly how information recycling improves the phase sensitivity, it is instructive to first consider the optical analogy of the incomplete QST process - a beamsplitter. We know from the quantum optics literature that a single-mode squeezed state ˆ b sq incident on a 50 / 50 beamsplitter leads to entanglement between the two output fields. Specifically, the two beamsplitter outputs are ˆ = ( ˆ a ϑ -ib ˆ sq ) / √ 2 and ˆ = ( ˆ b b sq -iϑ / ˆ ) √ 2, where ˆ is ϑ a vacuum input. The entanglement leads to a variance V ( ( ˆ X 0 a -ˆ X π/ 2 b ) / √ 2 ) = exp( -2 ), r which is less than for two uncorrelated coherent inputs, where we have defined the generalized quadrature of each output as ˆ X θ a b ( ) = ˆ( ˆ ) exp( a b iθ ) + ˆ ( ˆ a † b † ) exp( -iθ ). More generally,

FIG. 9. An atom interferometer enhanced by both singlemode squeezed light and information recycling. The transmitted component of ˆ ( b t 1 ) is interfered with a bright local oscillator, thereby allowing a homodyne measurement. The information from this measurement is then combined with the atomic signal. Information recycling gives an improvement to the sensitivity when the QST efficiency Q is less than unity.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0067v3-figure-9f924ad4-ee89-40a3-960c-e47cbbe8d13c.jpg)

for an asymmetric beamsplitter with beamsplitting ratio (i.e. reflection coefficient ) sin 2 ( θ QST / 2), the quantity

$$\underset { \text{er-} } { \text{ted} } \quad V \left ( \sin ( \theta _ { Q \text{ST} } / 2 ) \hat { X } _ { a } ^ { 0 } - \cos ( \theta _ { Q \text{ST} } / 2 ) \hat { X } _ { b } ^ { \pi / 2 } \right ) = e ^ { - 2 r } \ \ ( 3 1 )$$

is less than for two uncorrelated states. Essentially the same argument can be used to show that the informationrecycled signal (29) has a smaller variance, and therefore gives a smaller phase sensitivity, than the signal ˆ . Since S a ˆ b LO ( t 1 ) is a bright coherent state of complex amplitude √ N LO exp( iθ LO ), N LO glyph[greatermuch] 〈 ˆ ( N t b 1 ) 〉 and therefore

$$\hat { S } _ { b } \approx \sqrt { N _ { \text{LO} } } \hat { X } _ { b ( t _ { 1 } ) } ^ { \theta _ { \text{LO} } } \,. \quad \quad ( 3 2 )$$

Similarly, N a 1 ( t 1 ) glyph[greatermuch] 〈 ˆ N a 2 ( t 1 ) 〉 , and so at the most sensitive point of the atom interferometer, φ = π/ 2,

$$\hat { S } _ { a } \approx \sqrt { N _ { a _ { 1 } } ( t _ { 1 } ) } \hat { X } _ { a _ { 2 } ( t _ { 1 } ) } ^ { 0 } \,. \quad \quad ( 3 3 )$$

Thus, choosing θ LO = π/ 2 (i.e. we measure the phase quadrature of the transmitted photons), and

$$\det _ { \ v e e n } \quad \mathcal { G } = \sqrt { \frac { N _ { a _ { 1 } } ( t _ { 1 } ) } { N _ { \text{LO} } } \left ( \frac { 1 - \mathcal { Q } } { \mathcal { Q } } \right ) } = \sqrt { \frac { N _ { a _ { 1 } } ( t _ { 1 } ) } { N _ { \text{LO} } } } \cot ( \frac { \theta _ { \text{QST} } } { 2 } ), \ ( 3 4 )$$

yields the combined (information-recycled) signal

$$\underset { \substack { \i d s \ t o \\ \text{the is} \\ \text{re we} \\ \text{out as} \\ \text{and}... } } \quad \hat { S } \approx \frac { \sqrt { N _ { a _ { 1 } } ( t _ { 1 } ) } } { \sin ( \frac { \theta _ { Q S T } } { 2 } ) } \left ( \sin ( \frac { \theta _ { Q S T } } { 2 } ) \hat { X } _ { a _ { 2 } ( t _ { 1 } ) } ^ { 0 } - \cos ( \frac { \theta _ { Q S T } } { 2 } ) \hat { X } _ { b ( t _ { 1 } ) } ^ { \pi / 2 } \right ),$$

which can have smaller variance than ˆ . S a

Figure 8 plots the minimum phase sensitivity ∆ φ min as a function of QST efficiency Q corresponding to both the purely atomic signal ˆ S a (solid, blue curve) and the information-recycled signal ˆ = S ˆ S a - G ˆ S b (dashed, red curve). Although both instances suffer a degradation of sensitivity with poorer QST, this degradation is significantly arrested for the information-recycled signal. As a specific comparison, when Q = 0 5, . the analytic model predicts that information recycling gives ∆ φ min ≈ 0 035 . / √ N t , compared with ∆ φ min ≈ 0 71 . / √ N t in the absence of information recycling. Even at 10% QST efficiency, information recycling gives a sensitivity more than a factor of ten better than the SQL, whereas there is a negligible enhancement in the absence of information recycling. For very low levels of QST ( &lt; 0 05% for . r = r opt ), the information recycling scheme gives sensitivities above the SQL. This is because V ( S ) is very sensitive to slight imperfections in the estimates of the quadratures when QST is very low. Such imperfections arise due to the finite size of the condensate initially populating mode ˆ . For although the approximation (33) is a 1 exact in the limit of an infinitely large condensate, we are typically only working with 10 4 -10 6 atoms. Hence the deviation of ˆ S a from a perfect quadrature measurement is considerable, which is the cause of the discrepancy at low values of Q .

Figure 8 also compares the sensitivities predicted by the analytic undepleted reservoir solutions to TW simulations [see Eqs (23)]. Without information recycling, the agreement is excellent except near Q = 1. In contrast, TW simulations predict better sensitivities from the information-recycled signal ˆ S (by a factor between four and five), occurring at a squeezing parameter r TW larger than the analytic optimum r opt . As discussed in Sec. IV A, this improvement is due to ˆ a 1 deviating from a coherent state, leading to a reduced variance in ˆ J x at the output.

## C. Effect of losses on phase sensitivity

Quantifying the effects of losses on the sensitivity is an important experimental consideration. The simplest method of accounting for losses is by introducing virtual beamsplitters with a transmission coefficient of η that input vacuum noise ˆ at various points within the interϑ ferometry scheme. Specifically, this maps some mode ˆ c to √ ηc ˆ+ √ 1 -ηϑ ˆ .

Using this approach, we considered four types of losses:

- 1. Losses in the generation and transmission of the squeezed optical state before the QST process - i.e. ˆ = ˆ ( c b t 0 ).
- 2. Losses in the mode ˆ a 2 after the QST process, due to imperfections in the QST process such as spontaneous emission - i.e. ˆ = ˆ c a 2 ( t 1 ).
- 3. Losses in the transmitted optical state, including detection inefficiency - i.e. ˆ = ˆ ( c b t 1 ).
- 4. Symmetric losses within the atom interferometer, which also accounts for inefficient atom detection -i.e. ˆ = [ˆ c a 1 ( t 1 ) -ia ˆ ( 2 t 1 )] / √ 2 and ˆ = [ˆ ( c a 2 t 1 ) -ia ˆ ( 1 t 1 )] / √ 2.

FIG. 10. An illustration of the effects of losses ( η = 0 95) on . the phase sensitivity within the undepleted reservoir model for N t = 10 6 and r = 3. Although losses clearly degrade the sensitivity, in most cases information recycling ameliorates this degradation.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0067v3-figure-852c3d0a-4afb-4b71-a132-96e4632c127b.jpg)

Figure 10 illustrates the various effects of these losses on the sensitivity for an inefficiency of η = 0 95. . Unsurprisingly, losses degrade the phase sensitivity, both with and without the inclusion of information recycling. Nevertheless, it is interesting that any type of loss is never worse than losses affecting the initial squeezed optical state. Furthermore, information recycling still delivers sensitivities below the SQL, and for values of Q &gt; 5 -10% these are much better than what is possible without information recycling.

Since losses affect other squeezed-light enhanced atom interferometry schemes in a qualitatively similar fashion, we will not discuss losses further. We simply note that losses degrade the effects of squeezing, as they do in any optical squeezing experiment, and that if losses are not too great then information recycling can somewhat ameliorate the effects of this degradation.

## V. ENHANCED ATOM INTERFEROMETRY WITH TWO-MODE SQUEEZED LIGHT

We now consider the case of high frequency squeezed light, which is experimentally less challenging to generate than low frequency squeezed light. We model this as a state of light that is dominated by uncorrelated classical noise at sideband frequencies below some frequency

ω crit , while containing quantum correlations above this frequency. Specifically, the mode with frequency ω p + ω is correlated with the mode with frequency ω p -ω only when ω &gt; ω crit . Implementing the scheme from Sec. IV would result in atoms being outcoupled without quantum correlations [see Fig. 11(a)]. We can, however, exploit this source of high-frequency squeezing in order to generate an atomic state that displays two-mode squeezing. By adjusting the frequency of the classical control field from ω c to ω c + ω s , where ω s &gt; ω crit , the mode that is on two-photon resonance δ = 0 becomes far from the probe carrier frequency ω p , and outside the noisy low frequency region. Consequently, QST occurs within a spectral region of width ∆ ω centered around ω p + ω s . However, this does not yield an enhancement to the sensitivity of our atom interferometer, as a mode in this region displays no special quantum correlations in isolation; it only displays squeezing when considered in conjunction with the modes centered around ω p -ω s . We therefore include a second classical control beam of frequency ω c -ω s , which is on two-photon resonance and thereby causes QST at those frequencies between ω p -( ω s -∆ ω/ 2) and ω p -( ω s + ∆ ω/ 2) [see Fig. 11(b)]. By making these two control fields counter-propagating (i.e. wavevectors k 2 and -k 2 ), the atom-light interaction results in two correlated modes of outcoupled atoms with different momenta, which can be easily distinguished.

More precisely, modifying Hamiltonian (1) such that there are two counter-propagating classical control fields at frequencies ω c ± ω s and wavevectors ± k 2 , respectively, both shaped by the pulse envelope u c ( r , t ) = u trans ( x, y ) u prop ( z -ct ), gives

$$\hat { \mathcal { H } } = \sum _ { j = 1, 2 } \int d { \mathbf r } \, \hat { \psi } _ { j } ^ { \dagger } ( { \mathbf r } ) H _ { j } ( { \mathbf r } ) \hat { \psi } _ { j } ( { \mathbf r } ) + \hat { \mathcal { H } } _ { \text{light} }$$

$$& + \hbar { g } \overset { \prime - \cdot, \cdot } { \int d \mathbf r } \left [ \hat { \psi } _ { 1 } ( \mathbf r ) \hat { \psi } _ { 2 } ^ { \dagger } ( \mathbf r ) \hat { E } _ { 1 } ( \mathbf r ) u _ { c } ^ { * } ( \mathbf r, t ) \left ( e ^ { - i [ \mathbf k _ { 2 } \cdot \mathbf r - ( \omega _ { c } + \omega _ { s } ) t ] } \\ & + \ e ^ { - i [ - \mathbf k _ { 2 } \cdot \mathbf r - ( \omega _ { c } - \omega _ { s } ) t ] } \right ) + h. c. \right ]. \quad & ( 3 6 ) \quad \text{$\theta_{\theta}$} \\ & \quad. / N _ { \cdot }.$$

As in Sec. II, we assume that the entire atomic population is initially in hyperfine state | 1 〉 and in spatial mode u 0 ( r ). Conservation of energy and momentum implies that the modes resonant with the atom-light coupling are:

$$\hat { a } _ { 1 } = \int d { \mathbf r } \, u _ { 0 } ^ { * } ( { \mathbf r } ) \hat { \psi } _ { 1 } ( { \mathbf r } ), \quad \quad ( 3 7 \mathbf a ) \quad \text{$Not$} _ { \mathbf Q S }.$$

$$\hat { a } _ { \pm } = \int d { \mathbf r } \, u _ { 0 } ^ { * } ( { \mathbf r } ) e ^ { i ( { \mathbf k } _ { \pm } \mp { \mathbf k } _ { 2 } ) \cdot { \mathbf r } } \hat { \psi } _ { 2 } ( { \mathbf r } ), \quad ( 3 7 b ) \quad \overset { \text{two} } { \text{ref} } \quad \Phi$$

$$\hat { b } _ { \pm } ( t ) = \, \int \, d { \mathbf r } \, u _ { c } ^ { * } ( { \mathbf r }, t ) e ^ { - i { \mathbf k } _ { \pm } \cdot { \mathbf r } } \hat { E } _ { 1 } ( { \mathbf r } ), \quad ( 3 7 c ) \quad \text{corr}$$

where k ± = ( ω p ± ω s )ˆ z /c . Note that the mode expansion in Eq. (37b) is approximate, and only valid in the regime where the two mode functions u 0 ( r ) exp[ ( i k ± ∓ k 2 ) · r ] are approximately orthonormal. This is true provided the wavelength of the control fields is much smaller than the spatial extent of the condensate, or equivalently that the momentum kick imparted to the atoms is larger than the momentum width of the atomic cloud. Such a condition is easily satisfied in typical atom interferometers with Bose-condensed sources. Substituting Eqs (37) into Eq. (36) gives ˆ H ≈ H ˆ 0 + ˆ H ′ int , where

$$\text{ploit} \begin{array} { c } \kappa \cdots \ni \cdots & \cdots & \cdots \ni \cdots \\ \text{m} \end{array} \quad \text{.} \\ \text{gen-} \\ \text{zing.} \quad \hat { \mathcal { H } } _ { 0 } & = \left ( \hbar { \omega } _ { 2 } + \frac { \hbar { 2 } ^ { 2 } } { 2 m } | \mathbf k _ { 1 } - \mathbf k _ { 2 } | ^ { 2 } \right ) \left ( \hat { a } _ { + } ^ { \dagger } \hat { a } _ { + } + \hat { a } _ { - } ^ { \dagger } \hat { a } _ { - } \right ) \\ \text{field} \quad & + \hbar { ( \omega } _ { p } + \omega _ { s } ) \hat { b } _ { + } ^ { \dagger } \hat { b } _ { + } + \hbar { ( \omega } _ { p } - \omega _ { s } ) \hat { b } _ { - } ^ { \dagger } \hat { b } _ { - }, & ( 3 8 ) \\ \text{is on} \\ \text{robe} & \hat { \mathcal { H } } _ { \text{int} } ^ { \prime } & = \hbar { g } f ( t ) \hat { a } _ { 1 } ( \hat { a } _ { + } ^ { \dagger } e ^ { i ( \omega _ { c } + \omega _ { s } ) t } + \hat { a } _ { - } ^ { \dagger } e ^ { i ( \omega _ { c } - \omega _ { s } ) t } ) ( \hat { b } _ { + } + \hat { b } _ { - } ) \\ \text{ency} & + \hbar {. } c. \,, & ( 3 9 ) \\ \text{al re} \quad &. \quad. \quad. \quad. \quad. \quad. \quad. \quad.$$

where we have chosen the orientation of our control and probe fields such that | k + -k 2 | 2 = | k -+ k 2 | 2 ≡ | k 1 + k 2 | 2 , where k 1 = ( ω /c p ) ˆ z . Making the transformation ˆ H → ˆ H int = ˆ U † ˆ H ′ ˆ U with ˆ = exp( U - H i ˆ 0 t/ glyph[planckover2pi1] ) yields

$$\hat { \mathcal { H } } _ { \text{int} } & = \hbar { g } f ( t ) \left [ \hat { a } _ { + } ^ { \dagger } \hat { a } _ { 1 } \left ( \hat { b } _ { + } e ^ { - i \delta t } + \hat { b } _ { - } e ^ { i ( 2 \omega _ { s } - \delta ) t } \right ) \\ & + \hat { a } _ { - } ^ { \dagger } \hat { a } _ { 1 } \left ( \hat { b } _ { + } e ^ { - i ( 2 \omega _ { s } + \delta ) t } + \hat { b } _ { - } e ^ { - i \delta t } \right ) + h. c. \right ], ( 4 0 )$$

where δ = ω p -ω c -ω 2 - | glyph[planckover2pi1] k 1 -k 2 | 2 / (2 m ). If we assume ω s glyph[greatermuch] ∆ , and adjust ω ω c such that δ = 0, the contribution from the fast rotating terms can be neglected, and we obtain

$$\hat { \mathcal { H } } _ { \text{int} } \approx \hbar { g } f ( t ) \left [ \hat { a } _ { 1 } \left ( \hat { a } _ { + } ^ { \dagger } \hat { b } _ { + } + \hat { a } _ { - } ^ { \dagger } \hat { b } _ { - } \right ) + h. c. \right ]. \quad ( 4 1 )$$

The Heisenberg equations of motion for this Hamiltonian are [ c.f. Eqs (7)]:

$$i \dot { \hat { a } } _ { 1 } = g f ( t ) \left ( \hat { a } _ { + } \hat { b } _ { + } ^ { \dagger } + \hat { a } _ { - } \hat { b } _ { - } ^ { \dagger } \right ), \quad \ \ ( 4 2 \mathbf a )$$

$$i \dot { \hat { a } } _ { \pm } = g f ( t ) \hat { a } _ { 1 } \hat { b } _ { \pm }, & & ( 4 2 b )$$

$$i \dot { \hat { b } } _ { \pm } = g f ( t ) \hat { a } _ { 1 } ^ { \dagger } \hat { a } _ { \pm }.$$

In the regime of low depletion from the condensate, we make the undepleted reservoir approximation ˆ a 1 → √ N a 1 as in Sec. II A, and obtain

$$\hat { a } _ { \pm } ( t _ { 1 } ) = \hat { a } _ { \pm } ( t _ { 0 } ) \cos ( \frac { \theta _ { Q S T } } { 2 } ) - i \hat { b } _ { \pm } ( t _ { 0 } ) \sin ( \frac { \theta _ { Q S T } } { 2 } ) ( 4 3 \mathbf a )$$

$$\hat { b } _ { \pm } ( t _ { 1 } ) = \hat { b } _ { \pm } ( t _ { 0 } ) \cos ( \frac { \theta _ { Q S T } } { 2 } ) - i \hat { a } _ { \pm } ( t _ { 0 } ) \sin ( \frac { \theta _ { Q S T } } { 2 } ) ( 4 3 b )$$

$$\colon & \cos ( \frac { \theta _ { Q S T } } { 2 } ) - i \hat { b } _ { \pm } ( t _ { 0 } ) \sin ( \frac { \theta _ { Q S T } } { 2 } ) ( 4 3 \tt a ) \\ & \cos ( \frac { \theta _ { Q S T } } { 2 } ) - i \hat { a } _ { \pm } ( t _ { 0 } ) \sin ( \frac { \theta _ { Q S T } } { 2 } ) ( 4 3 b )$$

Note that the '+' and '-' modes decouple, and so the QST process depicted in Fig. 12 can be conceptualized as two independent atom-light beamsplitters with the same 'reflectivity'.

It is not difficult to define input states ˆ b ± ( t 0 ) that are correlated, and more specifically correspond to a twomode squeezed vacuum state. The argument follows that laid out in Sec. IV. First, we assume ˆ E 1 ( r , t ) ≈ u trans ( x, y ) ˆ E 1 ( z, t ), and rewrite Eq. (37c) in terms of the momentum modes of ˆ E 1 ( z, t ) [ c.f. Eq. (15)]:

$$\hat { b } _ { \pm } ( t ) = \int d q \, U _ { \text{prop} } ^ { * } ( q \mp q _ { s } ) \hat { \phi } ( q, t ), \quad \ ( 4 4 )$$

FIG. 11. (a): Single-mode outcoupling scheme. If 2 ω crit &gt; ∆ , then only the uncorrelated part of the spectrum of ω ˆ E 1 takes part in the outcoupling process. No enhancement to the sensitivity is possible in this case. (b): An outcoupling method for utilizing high frequency squeezing. Two control beams (frequencies ω c -ω s and ω c + ω s ) are used to outcouple from the region of the spectrum of ˆ E 1 centered around ω p -ω s and ω p + ω s , respectively. Since ∆ p glyph[greatermuch] ω s , the coupling strengths of the two control beams are nearly identical. Here δ 0 = ω p -ω c -ω 2 , which equals the change in the atomic kinetic energy glyph[planckover2pi1] | k 1 -k 2 | 2 / (2 m ) on two-photon resonance.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0067v3-figure-8cd2e096-e9d0-4e8d-ba2f-5b6d9413a4cf.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0067v3-figure-9dd62bd5-ef4a-4a9e-8d43-df8ff329f03a.jpg)

where q s = | k ± -k 1 | = ω /c s . Therefore, since Eq. (13) is the optical state output from the OPO:

the modes ˆ b ± ( t 0 ) comprise a two-mode squeezed vacuum state.

$$\hat { b } _ { \pm } ( t _ { 0 } ) = \left [ \int d q \, U _ { \text{prop} } ^ { * } ( q \mp q _ { s } ) \hat { \phi } ( q, t _ { i } ) \right ] \cosh r & & \text{Ou} \\ - \, i e ^ { i \theta _ { \text{sq} } } \left [ \int d q \, U _ { \text{prop} } ^ { * } ( q \mp q _ { s } ) \hat { \phi } ^ { \dagger } ( - q, t _ { i } ) \right ] \sinh r. \ \ ( 4 5 ) & & \text{all or} \\ & \text{and } \hat { a }$$

The term in the first set of square brackets is ˆ b ± ( t i ). Since U prop ( q ) is assumed real and symmetric,

$$\int d q \, U _ { \text{prop} } ^ { * } ( q \mp q _ { s } ) \hat { \phi } ^ { \dagger } ( - q, t _ { i } ) & = \int d q \, U _ { \text{prop} } ( q \pm q _ { s } ) \hat { \phi } ^ { \dagger } ( q, t ) \\ & = \hat { b } _ { \mp } ^ { \dagger } ( t _ { i } ).$$

Hence, we arrive at

$$\hat { b } _ { \pm } ( t _ { 0 } ) = \hat { b } _ { \pm } ( t _ { i } ) \cosh r - i e ^ { i \theta _ { s q } } \hat { b } _ { \mp } ^ { \dagger } ( t _ { i } ) \sinh r. \quad ( 4 7 ) \quad \frac { 1 } { 2 } \mathfrak { a } ^ { \dagger } _ { \neq \infty }$$

Since the initial quantum state is chosen such that

$$\hat { \phi } ( q, t _ { i } ) | \Psi ( 0 ) \rangle = \hat { a } _ { \pm } ( t _ { 0 } ) | \Psi ( 0 ) \rangle = 0, \quad ( 4 8 \tt a ) \quad \text{cho}$$

$$\hat { a } _ { 1 } ( t _ { 0 } ) | \Psi ( 0 ) \rangle = \sqrt { N _ { a _ { 1 } } ( t _ { 0 } ) } | \Psi ( 0 ) \rangle, \quad ( 4 8 b ) \quad \overset { \text{on a} } { \langle \Delta ^ { l } }$$

Our scheme that utilizes this two-mode squeezed optical vacuum is shown in Fig. 13. The modes ˆ b + ( t 0 ) and ˆ b -( t 0 ) form the inputs for the QST process, transferring all or part of their quantum state to the modes ˆ a + ( t 1 ) and ˆ a -( t 1 ) respectively [see Eqs (43)]. These two atomic modes form the input for the atom interferometer, where after the usual Mach-Zehnder interferometry sequence (beam split, reflect, beam split), the atom number difference at the outputs is measured:

$$\Delta \hat { N } _ { a } ( t _ { f } ) & \equiv \hat { N } _ { a _ { + } } ( t _ { f } ) - \hat { N } _ { a _ { - } } ( t _ { f } ) \\ & = 2 \hat { J } _ { z } \cos \phi + 2 \hat { J } _ { x } \sin \phi. \quad \quad ( 4 9 )$$

Here we have used the pseudo-spin operators ˆ J k ≡ 1 2 a † σ k a , where a = (ˆ a + ( t 1 ) , a ˆ -( t 1 )) T and σ k are the set of Pauli spin matrices; this notation is convenient for some of the expressions below. However, ∆ ˆ N a is a poor choice for our signal, since the average number difference of a two-mode squeezed vacuum state is zero, and hence 〈 ∆ ˆ N a 〉 = 0 even for the case of incomplete QST. Conse-

FIG. 12. QST scheme for a two-mode squeezed optical input. Two counter-propagating classical control beams E + (frequency ω c + ω s ) and E -(frequency ω c -ω s ) are used to implement a Raman transition with ˆ E 1 . Due to energymomentum resonance, E ± is resonant with the region of ˆ E 1 centered at frequency ω p ± ω s . This outcouples atoms from mode ˆ , and results in QST between the optical modes ˆ a 1 b ± and atomic modes ˆ a ± . Note that ˆ a ± recoils with momentum glyph[planckover2pi1] k 1 ± glyph[planckover2pi1] k 2 , where ± k 2 is the wavevector of E ± .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0067v3-figure-6efb4497-d973-41f2-8216-05f70cb95b02.jpg)

quently, we choose a signal based on the fluctuations of the output number difference:

$$\hat { S } _ { a } = \left ( \hat { N } _ { a _ { + } } ( t _ { f } ) - \hat { N } _ { a _ { - } } ( t _ { f } ) \right ) ^ { 2 }. \quad \ \ ( 5 0 ) \quad \ \text{is c}$$

## A. Complete quantum state transfer

When QST is perfect (i.e. θ QST = π ) and ˆ a ± ( t 1 ) = -ib ˆ ± ( t 0 ), the average and variance of the signal (50) are

$$\langle \hat { S } _ { a } \rangle & = \sinh ^ { 2 } ( 2 r ) \sin ^ { 2 } \phi, \quad & ( 5 1 a ) \quad \text{than} \quad \text{order} \\ V ( \hat { S } _ { a } ) & = \frac { 1 } { 2 } \sinh ^ { 2 } ( 2 r ) \left [ 1 + \cosh ( 4 r ) \left ( 3 - 4 \cos ( 2 \phi ) \right ) \right ] \quad \text{the s} \quad \text{of } \, 1 0 \\ & + \sinh ^ { 4 } ( 2 r ) \cos ( 4 \phi ), \quad & ( 5 1 b ) \quad \text{with} \quad \text{with} \quad \text{$N_{t} \sim$}$$

respectively. Note that these expressions are independent of the squeezing angle θ sq . Equations (51) yield a phase sensitivity of

$$\Delta \phi = \frac { \sqrt { 1 + \cosh ( 4 r ) \tan ^ { 2 } \phi } } { \sinh ( 2 r ) }. \quad \quad ( 5 2 )$$

A minimum sensitivity of ∆ φ min = 1 / sinh(2 r ) occurs at φ = π . However, since the total number of atoms detected at the interferometer output is N t = 〈 ˆ N a + ( t f )+ ˆ N a -( t f ) 〉 = 2sinh 2 r , we can rewrite

$$\Delta \phi _ { \min } = \frac { 1 } { \sqrt { N _ { t } ( N _ { t } + 2 ) } }, \quad \ \ ( 5 3 ) \quad ( \text{wl}$$

FIG. 13. Analogous optical circuit for an atom interferometer enhanced with two-mode squeezed light. The atom-light coupling (partially) maps the quantum state of the optical modes ˆ b + ( t 0 ) and ˆ b -( t 0 ) onto the atomic modes ˆ a + ( t 1 ) and ˆ a -( t 1 ), respectively. These two atomic modes form the inputs to a Mach-Zehnder atom interferometer. The modes ˆ b + and ˆ b -are co-propagating, but shown here to be spatially separated for the purposes of visual clarity. Likewise, the blue atomlight beamsplitter represents the condensate mode ˆ a 1 , which is common to both QST processes (see Fig. 12).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0067v3-figure-7f16fef6-eb6c-460e-8c67-70a61abb7e1e.jpg)

which is approximately the Heisenberg limit 1 /N t in the limit of large N t .

Na¨ ıvely, this interferometer scheme compares favourably to the single-mode squeezed-light enhanced scheme considered in Sec. IV. For instance, in order to achieve a sensitivity of ∆ φ ∼ 10 -5 , which is the sensitivity obtained by an ideal atom interferometer of 10 6 atoms enhanced by single-mode squeezed light with r = r TW ≈ 4 8 (see Fig. 6), we need to outcouple . N t ∼ 10 5 atoms, which requires a squeezing parameter of r ∼ 6 1. . However, the story is not quite as simple once the effects of incomplete QST and depletion are considered.

## B. Incomplete quantum state transfer and information recycling

For incomplete QST, the minimum phase sensitivity (which occurs at φ = tan -1 [( V ( ˆ J 2 z ) /V ( ˆ J 2 x )) 1 / 4 ]) is [81]

$$( \Delta _ { \ell } )$$

$$( \Delta \phi _ { \min } ) ^ { 2 } & = \frac { 2 \sqrt { V ( \tilde { J } _ { z } ^ { 2 } ) V ( \tilde { J } _ { z } ^ { 2 } ) } + C ( \tilde { J } _ { z } ^ { 2 }, \tilde { J } _ { z } ^ { 2 } ) + \langle ( \tilde { J } _ { x } \tilde { J } _ { z } + \tilde { J } _ { x } \tilde { J } _ { x } ) ^ { 2 } \rangle } { 4 ( \langle \tilde { J } _ { z } ^ { 2 } ) - ( \tilde { J } _ { z } ^ { 2 } ) ) ^ { 2 } }, \\ & = \frac { \sqrt { ( 1 - \Omega ) ( 1 + 5 ( 1 - \Omega ) N _ { t } ) \left ( 4 + 2 ( N _ { t } - 2 \Omega ) + \frac { N _ { t } ( 9 - \Omega ( 4 2 - 3 7 \Omega ) ) - 6 4 \Omega ^ { 2 } ( 1 - \Omega ) + 7 \Omega + 1 } { 4 ( N _ { t } + 2 \Omega ) ^ { 2 } } \right ) } { N _ { t } ( N _ { t } + 2 \Omega ) } \\ & + \frac { \left ( 1 + \frac { 1 - \Omega } { 2 ( N _ { t } + 2 \Omega ) } \right ) + N _ { t } ( 1 - \Omega ) \left ( 3 + \frac { 2 N _ { t } + 5 - \Omega } { 2 ( N _ { t } + 2 \Omega ) } \right ) } { N _ { t } ( N _ { t } + 2 \Omega ) },$$

where V ( ˆ ) X = 〈 ˆ X 2 〉 - 〈 ˆ X 〉 2 is the variance of ˆ , X C X,Y ( ˆ ˆ ) = 〈 ˆ ˆ XY + ˆ ˆ Y X 〉 -2 〈 ˆ X 〉〈 ˆ Y 〉 is the symmetrized covariance of ˆ X and ˆ , Y Q = sin ( 2 θ QST / 2) and N t = 2 Q sinh 2 r . As shown in Fig. 14, even small decreases from complete QST result in a rapid degradation of sensitivity. Indeed, provided N t glyph[greatermuch] 1 and (1 -Q glyph[greatermuch] ) 1 /N t , both of which are easily satisfied in practice, then

$$\Delta \phi _ { \min } \approx \sqrt { \frac { ( 1 -$$

Although sensitivities below the SQL are obtainable for Q &gt; (2 + √ 10) / 6 ≈ 0 86, . Heisenberg scaling is lost for extremely small perturbations from complete QST.

As for the interferometer enhanced by single-mode squeezed light, we can arrest the loss of sensitivity at incomplete QST via information recycling. Specifically, we directly measure the number of photons output after the QST process [i.e. ˆ N b ± = ˆ b † ± ( t 1 ) ˆ b ± ( t 1 )] via photodetection, and subtract these counts from the number of atoms outcoupled during the QST process (i.e. ˆ N a ± ). Explicitly, we construct the signal

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0067v3-figure-8d75ea4c-192b-45ee-bf13-36f9c4459018.jpg)

$$\hat { S } = \left ( \Delta \hat { N } _ { a } ( t$$

where ∆ ˆ N b ≡ ˆ N b + -ˆ N b -Thus, at the optimum phase φ = π , the minimum sensitivity with information recycling takes the remarkably simple form:

$$\begin{array} {$$

which for large N t is approximately the Heisenberg limit, but most importantly is almost independent of Q and approximately equal to the sensitivity when Q = 1 [see Eq. (53)].

Unlike the atom interferometry scheme enhanced with single-mode squeezed light, information recycling almost completely removes any degradation due to incomplete QST. However, this does not imply that the QST efficiency does not matter, since N t ∝ Q . Nevertheless,

Q

FIG. 14. Plots demonstrating the QST efficiency dependence of the minimum phase sensitivity ∆ φ min for an atom interferometer enhanced by two-mode optical squeezing [Eq. (55)]. The horizontal dashed lines that intercept each curve at Q = (2 + √ 10) / 6 indicate the SQL for N t = 10 6 (black, top), 10 5 (red, middle) and 10 4 (blue, bottom). The Heisenberg scaling is quickly lost for small perturbations of Q from unity. In contract, when information recycling is applied ∆ φ min remains at the Heisenberg limit for all Q [see Eq. (58)].

Eq. (58) represents the minimum sensitivity attainable for a fixed N t , and certainly gives better sensitivities than for the purely atomic signal ˆ . Perhaps surprisingly, this S a remains true even under the effects of depletion, as shown in the next subsection.

## C. Effects of depletion

The analytics presented in this section thus far have assumed the undepleted reservoir approximation ˆ a 1 → √ N a 1 . However, this approximation breaks down for even moderate values of r . Thus, following Sec. IV we treated mode ˆ a 1 as a quantum dynamical degree of freedom by performing numerical TW simulations of the QST process. Under TW, the operators are mapped to

complex stochastic amplitudes via the correspondences ˆ a i → α i and ˆ b i → β i . These stochastic variables evolve according to the SDEs [ c.f. Eqs (42)]:

$$i \dot { \alpha } _ { 1 } = g f ( t ) \left ( \alpha$$

$$i \dot { \alpha } _ { \pm } = g f ( t ) \alpha _ {$$

$$i \dot { \beta } _ { \pm } = g f ( t ) \alpha _ { 1 } ^ { * } \alpha _ { \pm }, \quad \quad \quad ( 5 9 c ) \quad \Big | \check { \gtrsim } ^ { \cdot }$$

with initial conditions

$$\alpha _ { 1 } ( 0 ) & = \sqrt { N _ { a _ { 1 } } ( t _ { 0 } ) } + \eta _ { \alpha _ { 1 } }, & ( 6 0 \text{a} ) \quad \stackrel { \Xi } { \Xi }$$

$$\alpha _ { \pm } ( 0 ) & = \eta _ { \alpha _ { \pm } }, & ( 6 0 b ) & \quad \check { \triangleleft 1 }$$

$$\beta _ { \pm } ( 0 ) = \eta _ { \beta _ { \pm } } \cosh r - i e ^ { i \theta _ { s q } } \eta _ { \beta _ { \mp } } ^ { * } \sinh r. \quad ( 6 0 c )$$

The η i are complex, independent Gaussian noises with zero mean and variance 1 / 2. Without loss of generality, we assumed a uniform f ( ) t for the TW simulations presented below.

The phase sensitivity without information recycling was calculated using Eq. (54) and the following expressions:

$$\langle \hat { J } _ { x } ^ { 2 } \rangle = \overline { \mathcal { J } _ { x } ^ { 2 } } - \frac { 1 } { 8 },$$

$$\langle \hat { J } _ { z } ^ { 2 } \rangle = \overline { \mathcal { J } _ { z } ^ { 2 } } - \frac { \vec { 1 } } { \vec { 8 } }, \quad \quad \quad ( 6 1 b )$$

$$\langle \hat { J } _ { x } ^ { 4 } \rangle = \overline { \mathcal { J } _ { x } ^ { 4 } } - \frac { \bar { 5 } } { 4 } \overline { \mathcal { J } _ { x } ^ { 2 } } + \frac { 1 } { 1 6 }, \quad \quad \quad ( 6 1 c ) \quad \text{fig} \quad \text{fig}$$

$$\langle \hat { J } _ { z } ^ { 4 } \rangle = \overline { \mathcal { J } _ { z } ^ { 4 } } - \frac { \bar { 5 } } { 4 } \overline { \mathcal { J } _ { z } ^ { 2 } } + \frac { \bar { 1 } } { 1 6 }, \quad \quad ( 6 1 d ) \quad \begin{array} { c } \text{$\text{$visit$} \\ \text{Eq.} \end{array}$$

$$\langle \hat { J } _ { x } ^ { 2 } \hat { J } _ { z } ^ { 2 } + \hat { J } _ { z } ^ { 2 } \hat { J } _ { x } ^ { 2 } \rangle = 2 \overline { \mathcal { J } _ { x } ^ { 2 } \mathcal { J } _ { z } ^ { 2 } } + \frac { 5 } { 4 } \overline { \mathcal { J } _ { x } ^ { 2 } } + \frac { 1 } { 4 } \overline { \mathcal { J } _ { z } ^ { 2 } } - \overline { | \alpha _ { + } | ^ { 2 } | \alpha _ { - } | ^ { 2 } }, \quad \text{$\emph{atoms} l$} \\ ( 6 1 e ) \quad \text{error} \ \cdot \quad \text{$\vdots$} \ \end{matrix}$$

$$\langle ( \hat { J } _ { x } \hat { J } _ { z } + \hat { J } _ { z } \hat { J } _ { x } ) ^ { 2 } \rangle = 4 \overline { \mathcal { J } _ { x } ^ { 2 } \mathcal { J } _ { z } ^ { 2 } } - \frac { 5 } { 2 } \overline { \mathcal { J } _ { x } ^ { 2 } } - \frac { 3 } { 2 } \overline { \mathcal { J } _ { z } ^ { 2 } } + \overline { | \alpha _ { + } | ^ { 2 } | \alpha _ { - } | ^ { 2 } }, \quad \text{$the$so$} \\ ( 6 1 f ) \quad \text{$eval tl$} \\ \text{$hold $fc$}$$

where J z = | α + ( t 1 ) | 2 - | α -( t 1 ) | 2 and J x = α + ( t 1 ) α ∗ -( t 1 ) + α ∗ + ( t 1 ) α -( t 1 ). The phase sensitivity with information recycling was computed using ∆ φ min = 1 / √ 4 〈 ˆ J 2 x 〉 [81].

Figure 15 shows that the effects of depletion cause a serious degradation to the sensitivity, even when QST is complete. The cause of this degradation is shown in the bottom panel of Fig. 15, which plots the variance in ˆ J z (which is proportional to the number difference of the two input ports of the atom interferometer). A large, nonzero variance implies that there exist low (but not negligible) probability trajectories where the quantum state of a photon in ˆ b + is successfully mapped to an atom in mode ˆ a + , but ˆ b -is not mapped to ˆ a -(and similarly for an exchange of '+' and '-'). The inclusion of unequal atom numbers in the two paths of the atom interferometer inevitably degrades the phase sensitivity. Fortunately, the effect is entirely reversed with information recycling, as measurement of the transmitted photons provides information about the full quantum correlations within the system.

FIG. 15. (Top) TW calculations of the minimum phase sensitivity, normalized to the sensitivity for complete QST [see Eq. (53)], with (red asterisks) and without (blue circles) information recycling, assuming an initial BEC of N a 1 ( t 0 ) = 10 6 atoms. N t was varied by adjusting the number of input photons N t b ( 0 ) = 2sinh 2 r from zero to ∼ 10 . 8 The standard error of each point is on the order of the point width, and the solid red line is the analytic solution (53), which is approximately the Heisenberg limit. The TW simulations reveal that the undepleted reservoir approximation ceases to hold for relatively small values of N t ; however, the negative effect of depletion is (almost) completely removed when information recycling is used. (Bottom) TW calculations of the variance in ˆ J z and the maximum possible QST efficiency, Q max = max t Q ( ), t as a function of the total number of atoms detected at the outputs. V ( J z ) has been normalized to N / t 4, which is the variance in ˆ J z for two uncorrelated coherent states each with a mean number of N / t 2. In the undepleted reservoir approximation, the analytic solution predicts V ( J z ) = 0 always. However, as N t increases beyond ∼ 10 , a full quantum treatment of mode ˆ 5 a 1 shows that there are increasingly large fluctuations in the number difference of the two MZ input ports, which suppress correlations between ˆ a + ( t 1 ) and ˆ a -( t 1 ).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0067v3-figure-7b69f06d-e73a-47ad-8c7f-e085cd4bc200.jpg)

The bottom panel of Fig. 15 shows that the maximum possible QST efficiency rapidly drops below 100% once N t is glyph[greaterorsimilar] 10% of the initial number of atoms populating mode ˆ . a 1 This is what we intuitively expect; it is impossible to outcouple more than N t = N a 1 ( t 0 ) atoms from the condensate, and so increasing r only results in a lower Q . Furthermore, in this regime depletion causes

a slight increase to the minimum phase sensitivity of the information-recycled signal. Here, the finite size of the BEC necessarily truncates the atom number probability distribution for the MZ input state at N = N a 1 ( t 0 ). Consequently, at best the input atomic state will have an atom number probability distribution corresponding to a two-mode squeezed optical state with a truncated 'tail'. When N t ≈ N a 1 ( t 0 ), this distribution will be somewhere between a twin-Fock state and a two-mode squeezed vacuum state, and as discussed in [81] this yields a minimum phase sensitivity of ∆ φ glyph[lessorsimilar] √ 2 /N t . This is consistent with the TW simulations shown in the top panel of Fig. 15.

## VI. A SINGLE INPUT ENHANCEMENT WITH TWO-MODE SQUEEZED-LIGHT AND INFORMATION RECYCLING

Although the previously presented two-mode squeezed light enhanced atom interferometry scheme has many advantages (e.g. it utilizes high frequency squeezing and can attain Heisenberg scaling), one disadvantage is that the number of atoms detected at the output is proportional to the number of atoms outcoupled during the QST process. Consequently, a high atom number interferometer requires both a large squeezing parameter and good QST efficiency. In contrast, the number of detected atoms in the scheme enhanced by single-mode squeezed light is independent of r and Q ; it only depends on the atom number of the initial condensate.

In this section we present an alternative scheme based on high frequency squeezing [i.e. a two-mode squeezed state; see Eq. (47)] that also utilizes all the atoms in the condensate, independent of the squeezing parameter and QST efficiency. This scheme is summarized in Fig. 16. Instead of using two control beams, this scheme uses just one, detuned from two-photon resonance by an amount ω s &gt; ω crit . This ensures that the QST process only occurs for a region of probe frequencies centered about ω p + ω s . This corresponds to the optical mode ˆ b + ( t 0 ) interacting with the initial condensate mode ˆ a 1 , leading to outcoupled atoms in mode ˆ . Since there is only one a 2 control beam, the portion of the spectrum correlated with this region (i.e. the region centered around ω p -ω s ) does not couple to the atoms, and is therefore transmitted. This transmitted light is described by the mode ˆ b -. By monitoring this transmitted light with homodyne detection, the quantum correlations of the two-mode squeezed state can still be used to enhance the sensitivity via information recycling. Furthermore, any incomplete QST can similarly be ameliorated by monitoring any transmitted light in the region centered around ω p + ω s . The outcoupled atomic mode and the remaining condensate are then used as the two inputs to an atom interferometer. The relevant (information-recycled) signal is then a linear combination of the atomic number difference at the output of the atom interferometer and the photon number differences of the transmitted light centered around

FIG. 16. Analgous optical circuit for modified atom interferometry utilizing two-mode optical squeezing. Only one of the modes (which corresponds to photons in a particular frequency band) is coupled to the atoms, while the other mode is measured directly. ˆ b + and ˆ b -are co-propagaiting, but shown here to be spatially separated to highlight their distiniguishability.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0067v3-figure-5275ced1-6a40-49f0-8591-dcb11c2cb9a7.jpg)

ω p + ω s and ω p -ω s . We call this scheme a single input enhancement with two-mode squeezed light, since only one mode of the two-mode squeezed optical vacuum directly interacts with the atoms. This is in contrast to the scheme considered in Sec. V, which we could describe as a double input enhancement with two-mode squeezed light, since both modes of the two-mode squeezed optical vacuum couple to the BEC.

We compute the phase sensitivity of this interferometry scheme from the information-recycled signal

$$\hat { S } = \sqrt { \mathcal { Q } } \hat { S } _ { a } - \mathcal { G } _ { + } \sqrt { 1 - \mathcal { Q } } \hat { S } _ { b _ { + } } + \mathcal { G } _ { - } \hat { S } _ { b _ { - } } \,, \quad ( 6 2 )$$

where ˆ S a is given by Eq. (11),

$$\hat { S } _ { b _ { \pm } } = \hat { b } ^ { \dagger } _ { \text{LO} _ { \pm } } ( t _ { f } ) \hat { b } _ { \text{LO} _ { \pm } } ( t _ { f } ) - \hat { b } ^ { \dagger } _ { \pm } ( t _ { f } ) \hat { b } _ { \pm } ( t _ { f } ), \quad ( 6 3 )$$

and

$$\mathcal { G } _ { \pm } = \sqrt { \frac { \langle \hat { a } _ { 1 } ^ { \dagger } ( t _ { 1 } ) \hat { a } _ { 1 } ( t _ { 1 } ) \rangle } { \langle \hat { b } _ { \text{LO} _ { \pm } } ^ { \dagger } ( t _ { 1 } ) \rangle } } \,. \quad \quad ( 6 4 )$$

Note that the photon modes at time t f are related to those at t 1 by the simple beamsplitting relations:

$$\hat { b } _ { \pm } ( t _ { f } ) = \frac { 1 } { \sqrt { 2 } } \left ( \hat { b } _ { \pm } ( t _ { 1 } ) - i \hat { b } _ { \text{LO} _ { \pm } } ( t _ { 1 } ) \right ) \quad \ \ ( 6 5 \text{a} )$$

$$\hat { b } _ { \text{LO} _ { \pm } } ( t _ { f } ) = \frac { \cdot } { \sqrt { 2 } } \left ( \hat { b } _ { \text{LO} _ { \pm } } ( t _ { 1 } ) - i \hat { b } _ { \pm } ( t _ { 1 } ) \right ), \quad ( 6 5 b )$$

and of course the trivial relation ˆ b -( t 1 ) = ˆ b -( t 0 ) holds. Our justification for this choice of signal follows an argument similar to that used to justify the choice of ˆ in S Sec. IV B. In brief, in the regime where depletion from

FIG. 17. The minimum phase sensitivity ∆ φ min (i.e. for φ = θ sq = π/ 2) for a single input two-mode optical squeezed vacuum enhanced atom interferometer, corresponding to a purely atomic signal ˆ S a (undepleted reservoir approximation: solid blue curve; TW: blue squares), partial information-recycled signal ˆ S a -G -ˆ S b -(undepleted reservoir approximation: green dashed curve; TW: green circles), and complete informationrecycled signal ˆ = S √ Q ˆ S a -G + √ 1 -Q ˆ S b + + G -ˆ S b -(undepleted reservoir approximation: red dot-dashed curve; TW: red crosses). All these plots assumed N t = 10 6 and r = r opt ≈ 3 8. . For convenience, these curves have been normalized relative to the SQL 1 / √ N t . The upper and lower horizontal dotted lines show the SQL, and the theoretical limit reached by a single-mode squeezed optical vacuum enhancement with perfect QST (i.e. 1 /N 3 / 4 t ), respectively.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0067v3-figure-1171cb10-01cb-480c-b771-2337b4c04af5.jpg)

the condensate is minimal during the QST process, the undepleted reservoir approximation ˆ a 1 → √ N a 1 gives the familiar atom-light beamsplitter relations:

$$\hat { a } _ { 2 } ( t _ { 1 } ) & = \hat { a } _ { 2 } ( t _ { 0 } ) \cos ( \frac { \theta _ { Q S T } } { 2 } ) - i \hat { b } _ { + } ( t _ { 0 } ) \sin ( \frac { \theta _ { Q S T } } { 2 } ), ( 6 6 a ) \quad \text{choic} \\ \hat { b } _ { + } ( t _ { 1 } ) & = \hat { b } _ { + } ( t _ { 0 } ) \cos ( \frac { \theta _ { Q S T } } { 2 } ) - i \hat { a } _ { 2 } ( t _ { 0 } ) \sin ( \frac { \theta _ { Q S T } } { 2 } ), ( 6 6 b ) \quad \text{be a} \\ \quad \text{detur}$$

$$\nu _ { 2 } ^ { 2 } \nu _ { 2 } ^ { 2 } \nu _ { 0 / 2 } \nu _ { 0 / 2 } \nu _ { 0 / 2 } ^ { 2 } \nu _ { 0 / 2 } ^ { 2 }$$

and Q = sin ( 2 θ QST / 2). Since ˆ S a ≈ √ N a 1 ( t 1 ) ˆ X 0 a 2 ( t 1 ) and ˆ S b ± ≈ √ N LO ± ˆ X θ LO ± b ± ( t 1 ) , then for θ LO ± = π/ 2:

$$\hat { S } \approx \sqrt { N _ { a _ { 1 } } ( t _ { 1 } ) } \left ( \hat { X } _ { b _ { - } ( t _ { 0 } ) } ^ { \pi / 2 } - \hat { X } _ { b _ { + } ( t _ { 0 } ) } ^ { \pi / 2 } \right ), \quad ( 6 7 ) \quad \text{furl}$$

which for optimal squeezing angle θ sq = π/ 2 has a variance V ( S ) = 2 √ N a 1 ( t 1 ) exp( -2 ), which can be smaller r than V ( S a ).

The dependence of the minimum sensitivity on the QST efficiency is shown in Fig. 17. There are clear similarities to the plots shown in Fig. 8 corresponding to the scheme enhanced by single-mode squeezed optical vacuum. Nevertheless, we highlight some important differences. Firstly, the pure atomic signal ˆ S a always gives a sensitivity worse than the SQL. Secondly, when the QST efficiency is close to 100%, it is sufficient to use a partial information-recycled signal ˆ S a -G -ˆ S b -. Thirdly, even for the signal ˆ the sensitivity remains above the S analytic limit ∆ φ = 1 /N 3 / 4 t derived for the single-mode squeezed optical vacuum enhanced atom interferometer, since for finite levels of squeezing measurements of ˆ S b -are not perfectly correlated with √ Q ˆ S a -G + √ 1 -Q ˆ S b + . Finally, the TW simulations predict that the effects of depletion lead to a poorer sensitivity than that given by the analytic sensitivity for the information-recycled signal ˆ . S

## VII. SUMMARY AND DISCUSSION.

In this paper we have shown how squeezed light can be used to enhance the sensitivity of atom interferometers. We have specifically considered three schemes: (Sec. IV) enhancement with a single-mode squeezed optical vacuum (i.e. low frequency squeezing), (Sec. V) double input enhancement with two-mode squeezed optical vacuum (i.e. high frequency squeezing), and (Sec. VI) single input enhancement with two-mode squeezed optical vacuum. We have shown that all three schemes give sensitivities below the SQL - even when the effects of depletion from the initial condensate and incomplete QST are included. Furthermore, we have demonstrated that information recycling provides a further enhancement to the sensitivity when QST between the atoms and light is incomplete.

Table I provides a quantitative comparison of the sensitivities for the three schemes for complete QST, and for incomplete QST ( Q = 0 2) with and without information . recycling. This is a concise demonstration of the different sensitivities obtainable (including the effects of depletion), the degrading effects of incomplete QST, and how information recycling ameliorates this degradation. The choice of Q = 0 2 is not unrealistic given current technol-. ogy. For example, the 780 nm D2 transition of 87 Rb can be addressed with a Rabi frequency of Ω ∼ 20 MHz and detunings ∆ p ≈ ∆ c ∼ 70 GHz. Since the electric dipole moment of this transition is d 12 ∼ 2 × 10 -29 C/m, Eq. (2) implies that an effective coupling of g ∼ 65 × 10 -3 m 3 / 2 /s is obtainable. If we assume a condensate of Gaussian spatial profile u 0 ( r ), with transverse width R ⊥ ∼ 15 µ m, and further assume the pulse envelope u p (r , t ) is also Gaussian with the same transverse width R ⊥ and pulse duration T ∼ 1 ms, then

$$\theta _ { \text{QST} } = \frac { 4 } { 3 } ( 2 \pi ) ^ { 1 / 4 } g \sqrt { \frac { N T } { c A _ { \perp } } } \sim \frac { 3 \pi } { 1 0 }, \quad \quad ( 6 8 )$$

where A ⊥ = πR 2 ⊥ and we have used N = 10 . This gives 6 a QST efficiency of Q ∼ 20%.

Furthermore, losses due to spontaneous emission during the QST process are small in this regime. An estimate for the (time-dependent) rate at which atoms in

|                                                                                                                         | ∆ φ min            | ∆ φ min           | ∆ φ min         |
|-------------------------------------------------------------------------------------------------------------------------|--------------------|-------------------|-----------------|
| Interferometer scheme                                                                                                   | ˆ S a and Q = 100% | ˆ S a and Q = 20% | ˆ S and Q = 20% |
| Enhancement with single-mode squeezed state, r = 4 . 8, and N t = 10 6 (Sec. IV)                                        | 1 × 10 - 5         | 9 × 10 - 4        | 1 . 9 × 10 - 5  |
| Double input enhancement with two-mode squeezed state and                                                               | 1 . 6 × 10 - 6     | 2 . 1 × 10 - 3    | 2 × 10 - 6      |
| N t = 6 . 2 × 10 5 (Sec. V) Single input enhancement with two-mode squeezed state, r = 3 . 8 and N t = 10 6 (Sec. VI) a | 3 . 2 × 10 - 2     | 1 . 4 × 10 - 2    | 1 . 6 × 10 - 4  |

a Note that ˆ gives ∆ S φ = 3 8 . × 10 -5 for Q = 100%.

TABLE I. Numerical values of the minimum phase sensitivities for the three different squeezed-light enhanced atom interferometry schemes, in three distinct scenarios. These values were obtained via TW simulations. All three schemes assumed an initial condensate of N a 1 ( t 0 ) = 10 6 atoms; a large but achievable [82] number of atoms. All these atoms can be detected at the outputs for the two schemes utilizing a single port enhancement (i.e. the first and third rows of the table). In contrast, N t depends on the value of r for the double input enhancement (i.e. the middle row of the table); the choice N t ∼ 6 2 . × 10 5 corresponds to Q ≈ 0 2 and . r ≈ 7 8. . For comparison, an atom interferometer operating at the SQL with N a 1 ( t 0 ) = 10 6 atoms has a sensitivity of ∆ φ = 10 -3 .

the excited state | 3 〉 are lost due to spontaneous emission is Γ( ) t = Ω eff ( ) t γ/ ∆ , where p γ = 2 π × 6 07 MHz . is the natural linewidth for 87 Rb and Ω eff is the effective Rabi frequency between modes ˆ a 1 and ˆ a 2 during the QST process. Assuming a beamsplitter-like coupling [such as described by Eqs (10)], then N a 1 ( t 1 ) = N t cos 2 ( ∫ dt ′ Ω eff ( t ′ ) / 2) = N t - 〈 ˆ ( N t b 0 ) 〉 sin 2 ( θ QST / 2), where the final equality follows from conservation of atom number. Therefore, the total loss is

$$l = \int d t ^ { \prime } \, \Gamma ( t ^ { \prime } ) = \frac { \gamma } { \Delta _ { p } } \cos ^ { - 1 } \sqrt { 1 - \frac { \langle \hat { N } _ { b } ( t _ { 0 } ) \rangle } { N _ { t } } \sin ^ { 2 } ( \frac { \theta _ { Q S T } } { 2 } ) }. \quad \stackrel { \text{comp} } { \text{Morti} } \\ \pi _ { \dots } \, - \, - \, - \, - \, - \, - \, - \, - \, - \, - \, - \, - \, - \, - \, - \, - \, - \, - \, - \, - \, - \, - \, - \, - \, - \, - \, - \, - \, - \, - \, - \, - \, - \, - \, - \, - \, - \, - \, - \, - \, - \, - \, - \, - \, - \, - \, - \, - \, - \, - \, - \, - \, - \, - \, - \, - \, - \, - \, - \, - \, - \, - \, - \, - \, - \, - \, - \tau } \, \lambda \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nu \, \nolimits }$$

For a squeezing parameter of r = 3 5 (i.e. . ˆ ( N t b 0 ) ∼ 270) and the parameters specified above, l ∼ 6 5 . × 10 -7 .

[1] A. Peters, K. Y. Chung, and S. Chu, Nature 400 , 849 (1999).

[2] H. M¨ uller, S.-w. Chiow, S. Herrmann, S. Chu, and K.-Y. Chung, Phys. Rev. Lett. 100 , 031101 (2008).

However, l is an estimate of the fraction of the total atoms lost; what is more important is the number of atoms lost compared to the total number of atoms in ˆ , which is a 2 lN /N t a 2 ∼ 1%.

As mentioned throughout the paper, there are some important differences between the three schemes. The single-mode squeezed optical vacuum enhancement has the advantages of conceptual simplicity and full utilization of all the atoms in the condensate. However, the generation of a single-mode squeezed optical vacuum occurs at low frequencies which, while possible, is technically challenging. The double input enhancement with two-mode squeezed optical vacuum uses technically less challenging high-frequency squeezing, and its approximate Heisenberg scaling independent of the QST efficiency with information recycling - is clearly the best of the three schemes. However, achieving a large N t , and therefore a small absolute sensitivity, requires a relatively large squeezing parameter r and good QST efficiency Q . The single input enhancement with two-mode squeezed optical vacuum can be thought of as a compromise between the former two schemes. Although it attains a minimum sensitivity similar in magnitude to the singlemode squeezed optical vacuum enhancement, it utilizes both high-frequency squeezing and all the atoms in the initially prepared condensate.

Ultimately, however, all three schemes attain sensitivities that are substantially below the SQL, that are robust to imperfect QST once information recycling is incorporated, and are strongly compatible with the practical requirements of current state-of-the-art atom interferometry. This provides a compelling case for the further development of atom interferometers that are enhanced by squeezed optical states and/or information recycling.

## ACKNOWLEDGEMENTS

We would like to acknowledge useful discussions with Warwick Bowen, John Close, Joel Corney, Saleh RahimiKeshari and Nicholas Robins. Numerical simulations were performed using XMDS2 [83] on the University of Queensland (UQ) School of Mathematics and Physics computer 'Obelix', with thanks to Elliott Hilaire and Ian Mortimer for computing support. SSS acknowledges the support of Ian P. McCulloch and the Australian Research Council (ARC) Centre of Excellence for Engineered Quantum Systems (project no. CE110001013). SAH acknowledges the support of ARC Project DE130100575.

[3] P. A. Altin, M. T. Johnsson, V. Negnevitsky, G. R. Dennis, R. P. Anderson, J. E. Debs, S. S. Szigeti, K. S. Hardman, S. Bennetts, G. D. McDonald, L. D. Turner, J. D. Close, and N. P. Robins, New J. Phys. 15 , 023009 (2013).

- [4] B. Canuel, F. Leduc, D. Holleville, A. Gauguet, J. Fils, A. Virdis, A. Clairon, N. Dimarcq, C. J. Bord´, A. Lane dragin, and P. Bouyer, Phys. Rev. Lett. 97 , 010402 (2006).
- [5] A. Lenef, T. D. Hammond, E. T. Smith, M. S. Chapman, R. A. Rubenstein, and D. E. Pritchard, Phys. Rev. Lett. 78 , 760 (1997).
- [6] T. L. Gustavson, P. Bouyer, and M. A. Kasevich, Phys. Rev. Lett. 78 , 2046 (1997).
- [7] M. J. Snadden, J. M. McGuirk, P. Bouyer, K. G. Haritos, and M. A. Kasevich, Phys. Rev. Lett. 81 , 971 (1998).
- [8] J. M. McGuirk, G. T. Foster, J. B. Fixler, M. J. Snadden, and M. A. Kasevich, Phys. Rev. A 65 , 033608 (2002).
- [9] M. Vengalattore, J. M. Higbie, S. R. Leslie, J. Guzman, L. E. Sadler, and D. M. Stamper-Kurn, Phys. Rev. Lett. 98 , 200801 (2007).
- [10] S. Gupta, K. Dieckmann, Z. Hadzibabic, and D. E. Pritchard, Phys. Rev. Lett. 89 , 140401 (2002).
- [11] R. Bouchendira, P. Clad´ e, S. Guellati-Kh´lifa, e F. Nez, and F. Biraben, Phys. Rev. Lett. 106 , 080801 (2011).
- [12] J. B. Fixler, G. T. Foster, J. M. McGuirk, and M. A. Kasevich, Science 315 , 74 (2007).
- [13] G. Lamporesi, A. Bertoldi, L. Cacciapuoti, M. Prevedelli, and G. M. Tino, Phys. Rev. Lett. 100 , 050801 (2008).
- [14] M. Andia, R. Jannin, F. Nez, F. Biraben, S. GuellatiKh´ elifa, and P. Clad´ e, Phys. Rev. A 88 , 031605 (2013).
- [15] G. Rosi, F. Sorrentino, L. Cacciapuoti, M. Prevedelli, and G. M. Tino, Nature 510 , 518 (2014).
- [16] S. Fray, C. A. Diez, T. W. H¨nsch, a and M. Weitz, Phys. Rev. Lett. 93 , 240404 (2004).
- [17] S. Dimopoulos, P. W. Graham, J. M. Hogan, and M. A. Kasevich, Phys. Rev. Lett. 98 , 111102 (2007).
- [18] D. Schlippert, J. Hartwig, H. Albers, L. L. Richardson, C. Schubert, A. Roura, W. P. Schleich, W. Ertmer, and E. M. Rasel, Phys. Rev. Lett. 112 , 203002 (2014).
- [19] G. Amelino-Camelia, C. L¨mmerzahl, F. Mercati, a and G. M. Tino, Phys. Rev. Lett. 103 , 171302 (2009).
- [20] G. M. Tino and F. Vetrano, Class. Quantum Grav. 24 , 2167 (2007).
- [21] S. Dimopoulos, P. W. Graham, J. M. Hogan, M. A. Kasevich, and S. Rajendran, Phys. Rev. D 78 , 122002 (2008).
- [22] N. P. Robins, P. A. Altin, J. E. Debs, and J. D. Close, Physics Reports 529 , 265 (2013).
- [23] K. V. Kheruntsyan, M. K. Olsen, and P. D. Drummond, Phys. Rev. Lett. 95 , 150405 (2005).
- [24] H. Pu and P. Meystre, Phys. Rev. Lett. 85 , 3987 (2000).
- [25] B. L¨ ucke, M. Scherer, J. Kruse, L. Pezz´ e, F. Deuretzbacher, P. Hyllus, O. Topic, J. Peise, W. Ertmer, J. Arlt, L. Santos, A. Smerzi, and C. Klempt, Science 334 , 773 (2011).
- [26] C. Gross, H. Strobel, E. Nicklas, T. Zibold, N. Bar-Gill, G. Kurizki, and M. K. Oberthaler, Nature 480 , 219 (2011).
- [27] J.-C. Jaskula, M. Bonneau, G. B. Partridge, V. Krachmalnicoff, P. Deuar, K. V. Kheruntsyan, A. Aspect, D. Boiron, and C. I. Westbrook, Phys. Rev. Lett. 105 , 190402 (2010).
- [28] R. Bucker, J. Grond, S. Manz, T. Berrada, T. Betz, C. Koller, U. Hohenester, T. Schumm, A. Perrin, and J. Schmiedmayer, Nat Phys 7 , 608 (2011).
- [29] S. A. Haine and A. J. Ferris, Phys. Rev. A 84 , 043624 (2011).
- [30] R. J. Lewis-Swan and K. V. Kheruntsyan, Nat Commun 5 (2014).
- [31] M. Kitagawa and M. Ueda, Phys. Rev. A 47 , 5138 (1993).
- [32] A. Søndberg Sørensen, Phys. Rev. A 65 , 043610 (2002).
- [33] C. Gross, T. Zibold, E. Nicklas, J. Est`ve, e and M. K. Oberthaler, Nature 464 , 1165 (2010).
- [34] M. F. Riedel, P. B¨hi, Y. Li, T. W. H¨nsch, A. Sinatra, o a and P. Treutlein, Nature 464 , 1170 (2010).
- [35] M. T. Johnsson and S. A. Haine, Phys. Rev. Lett. 99 , 010401 (2007).
- [36] S. A. Haine and M. T. Johnsson, Phys. Rev. A 80 , 023611 (2009).
- [37] S. A. Haine, J. Lau, R. P. Anderson, and M. T. Johnsson, Phys. Rev. A 90 , 023613 (2014).
- [38] H. Jing, J.-L. Chen, and M.-L. Ge, Phys. Rev. A 63 , 015601 (2000).
- [39] M. Fleischhauer and S. Gong, Phys. Rev. Lett. 88 , 070404 (2002).
- [40] S. A. Haine and J. J. Hope, Phys. Rev. A 72 , 033601 (2005).
- [41] S. A. Haine and J. J. Hope, Laser Phys. Lett. 2 , 597 (2005).
- [42] S. A. Haine, M. K. Olsen, and J. J. Hope, Phys. Rev. Lett. 96 , 133601 (2006).
- [43] K. Hammerer, A. S. Sørensen, and E. S. Polzik, Rev. Mod. Phys. 82 , 1041 (2010).
- [44] C. M. Caves, Phys. Rev. D 23 , 1693 (1981).
- [45] J. E. Debs, P. A. Altin, T. H. Barter, D. D¨ring, G. R. o Dennis, G. McDonald, R. P. Anderson, J. D. Close, and N. P. Robins, Phys. Rev. A 84 , 033610 (2011).
- [46] S. S. Szigeti, J. E. Debs, J. J. Hope, N. P. Robins, and J. D. Close, New J. Phys. 14 , 023009 (2012).
- [47] S. A. Haine, A. J. Ferris, J. D. Close, and J. J. Hope, Phys. Rev. A 69 , 013605 (2004).
- [48] S. S. Szigeti, M. R. Hush, A. R. R. Carvalho, and J. J. Hope, Phys. Rev. A 80 , 013614 (2009).
- [49] S. S. Szigeti, M. R. Hush, A. R. R. Carvalho, and J. J. Hope, Phys. Rev. A 82 , 043632 (2010).
- [50] S. S. Szigeti, S. J. Adlong, M. R. Hush, A. R. R. Carvalho, and J. J. Hope, Phys. Rev. A 87 , 013626 (2013).
- [51] M. R. Hush, S. S. Szigeti, A. R. R. Carvalho, and J. J. Hope, New J. Phys. 15 , 113060 (2013).
- [52] S. A. Haine, Phys. Rev. Lett. 110 , 053002 (2013).
- [53] A. Sinatra, F. Castelli, L. A. Lugiato, P. Grangier, and J. P. Poizat, Quantum Semiclass. Opt. 7 , 405 (1995).
- [54] E. Brion, L. H. Pedersen, and K. Mølmer, J. Phys. A: Math. Theor. 40 , 1033 (2007).
- [55] V. Paulisch, H. Rui, H. Ng, and B.-G. Englert, Eur. Phys. J. Plus 129 , 12 (2014), 10.1140/epjp/i2014-140128.
- [56] A. O. Jamison, J. N. Kutz, and S. Gupta, Phys. Rev. A 84 , 043643 (2011).
- [57] P. A. Altin, G. McDonald, D. D¨ring, J. E. Debs, T. H. o Barter, J. D. Close, N. P. Robins, S. A. Haine, T. M. Hanna, and R. P. Anderson, New J. Phys. 13 , 065020 (2011).
- [58] P. A. Altin, G. McDonald, D. D¨ring, J. E. Debs, T. H. o Barter, J. D. Close, N. P. Robins, S. A. Haine, T. M. Hanna, and R. P. Anderson, New J. Phys. 13 , 119401 (2011).
- [59] M. P. Fewell, Opt. Commun. 253 , 125 (2005).
- [60] M. Kasevich and S. Chu, Appl. Phys. B 54 , 321 (1992), 10.1007/BF00325375.
- [61] D. M. Giltner, R. W. McGowan, and S. A. Lee, Phys. Rev. A 52 , 3966 (1995).

- [62] H. M¨ uller, S.-w. Chiow, and S. Chu, Phys. Rev. A 77 , 023609 (2008).
- [63] P. Clad´ e, T. Plisson, S. Guellati-Kh´lifa, e F. Nez, and F. Biraben, Eur. Phys. J. D 59 , 349 (2010-09-01).
- [64] S.-w. Chiow, T. Kovachy, H.-C. Chien, and M. A. Kasevich, Phys. Rev. Lett. 107 , 130403 (2011).
- [65] K. S. Hardman, C. C. N. Kuhn, G. D. McDonald, J. E. Debs, S. Bennetts, J. D. Close, and N. P. Robins, Phys. Rev. A 89 , 023626 (2014).
- [66] M. O. Scully and M. S. Zubairy, Quantum Optics , 1st ed. (Cambridge University Press, 1997).
- [67] L. Pezz´ e and A. Smerzi, Phys. Rev. Lett. 102 , 100401 (2009).
- [68] M. S. Stefszky, C. M. Mow-Lowry, S. S. Y. Chua, D. A. Shaddock, B. C. Buchler, H. Vahlbruch, A. Khalaidovski, R. Schnabel, P. K. Lam, and D. E. McClelland, Class. Quantum Grav. 29 , 145015 (2012).
- [69] D. F. Walls and G. J. Milburn, Quantum Optics , 2nd ed. (Springer-Verlag, Berlin and Heidelberg, 2008).
- [70] Actually, this phase factor exp( iqct ) is an artefact of writing the pulse shape u p ( r , t ) exp( i k 1 · r ) as a single frequency pulse with a slowly-varying envelope. Strictly, this pulse contains a range of frequencies about k 1 which exactly cancel exp( iqct ) in the above integral.
- [71] L. Pezz´ e and A. Smerzi, Phys. Rev. Lett. 100 , 073601 (2008).
- [72] M. D. Lang and C. M. Caves, Phys. Rev. Lett. 111 , 173601 (2013).
- [73] M. J. Steel, M. K. Olsen, L. I. Plimak, P. D. Drummond, S. M. Tan, M. J. Collett, D. F. Walls, and R. Graham, Phys. Rev. A 58 , 4824 (1998).
- [74] P. B. Blakie, A. S. Bradley, M. J. Davis, R. J. Ballagh, and C. W. Gardiner, Advances in Physics 57 , 363 (2008).
- [75] A. Polkovnikov, Ann. Phys. (N.Y.) 325 , 1790 (2010).
- [76] B. Opanchuk and P. D. Drummond, J. Math. Phys. 54 , 042107 (2013).
- [77] C. W. Gardiner and P. Zoller, Quantum Noise: A Handbook of Markovian and Non-Markovian Quantum Stochastic Methods with Applications to Quantum Optics , 3rd ed. (Springer, Berlin and Heidelberg, 2004).
- [78] A. Sinatra, C. Lobo, and Y. Castin, J. Phys. B: At. Mol. Opt. Phys. 35 , 3599 (2002).
- [79] M. T. Johnsson, G. R. Dennis, and J. J. Hope, New J. Phys. 15 , 123024 (2013).
- [80] M. Olsen and A. Bradley, Opt. Commun. 282 , 3924 (2009).
- [81] S. A. Haine, S. S. Szigeti, M. D. Lang, and C. M. Caves, arXiv:1411.5111 (2014).
- [82] K. M. R. van der Stam, E. D. van Ooijen, R. Meppelink, J. M. Vogels, and P. van der Straten, Rev. Sci. Instrum. 78 , 013102 (2007).
- [83] G. R. Dennis, J. J. Hope, and M. T. Johnsson, Comput. Phys. Commun. 184 , 201 (2013).
Stuart S. Szigeti Behnam Tonekaboni Wing Yung S. Lau Samantha N. Hood Simon A. Haine
Squeezed-light-enhanced atom interferometry below the standard quantum limit
We investigate the prospect of enhancing the phase sensitivity of atom
interferometers in the Mach-Zehnder configuration with squeezed light.
Ultimately, this enhancement is achieved by transferring the quantum state of
squeezed light to one or more of the atomic input beams, thereby allowing
operation below the standard quantum limit. We analyze in detail three specific
schemes that utilize (1) single-mode squeezed optical vacuum (i.e. low
frequency squeezing), (2) two-mode squeezed optical vacuum (i.e. high frequency
squeezing) transferred to both atomic inputs, and (3) two-mode squeezed optical
vacuum transferred to a single atomic input. Crucially, our analysis considers
incomplete quantum state transfer (QST) between the optical and atomic modes,
and the effects of depleting the initially-prepared atomic source.
Unsurprisingly, incomplete QST degrades the sensitivity in all three schemes.
We show that by measuring the transmitted photons and using information
recycling [Phys. Rev. Lett. 110, 053002 (2013)], the degrading effects of
incomplete QST on the sensitivity can be substantially reduced. In particular,
information recycling allows scheme (2) to operate at the Heisenberg limit
irrespective of the QST efficiency, even when depletion is significant.
Although we concentrate on Bose-condensed atomic systems, our scheme is equally
applicable to ultracold thermal vapors.