Precision inertial sensing with quantum gases

Quantum sensors based on light-pulse atom interferometers allow for high-precision measurements of inertial and electromagnetic forces such as the accurate determination of fundamental constants as the fine structure constant or testing foundational laws of modern physics as the equivalence principle. These schemes unfold their full performance when large interrogation times and/or large momentum transfer can be implemented. In this article, we demonstrate how precision interferometry can benefit from the use of Bose-Einstein condensed sources when the state of the art is challenged. We contrast systematic and statistical effects induced by Bose-Einstein condensed sources with thermal sources in three exemplary science cases of Earth- and space-based sensors.


Introduction
Atom interferometers are mainly used for inertially sensitive measurements [1] and a variety of tests of fundamental physics [2,3]. Key levers to increase the sensitivity are the transfer of a large number of photons during the beamsplitting processes, extending the time of free fall while maintaining contrast and atomic flux. At the same time, the characterization of errors requires an increased level of control over the manipulation and preparation of atoms. Limitations of interferometers operating with molassescooled atoms and their mitigation by reducing the residual expansion rates were studied theoretically [4] and experimentally [5,6].
In this paper, together with the trade-off between flux and expansion rate, we contrast the appropriateness of the two regimes of atomic ensembles to perform precision tests by evaluating the respective shot noise, cycle times, excitation rates and most prominent systematics such as gravity gradients (GGs), Coriolis force, wave-front aberrations (WFA) and mean-field interactions.
Proposals for space missions, in particular, rely on Delta-Kick collimation (DKC) via optical or magnetic potentials to exploit extended times of free fall in microgravity [7,8,9,10] and achieve extremely low wave packet ex-Send offprint requests to: gaaloul@iqo.uni-hannover.de pansion rates, corresponding to pK temperatures in thermal ensembles. Bose-Einstein-condensed (BEC) ensembles are better suited for DKC aiming at long interrogation times [4], but seem to suffer from a reduced atomic flux due to the evaporation despite recent promising studies [11]. Molasses-cooled atoms feature a higher number of atoms, but are typically velocity-filtered in 1D [12], which ultimately implies a lower flux of atoms as we will detail in our study.
To illustrate our comparative study between condensed and thermal sources, we consider three prominent cases for free-fall atom interferometers: a gravimeter, a gravitygradiometer, and a test of the universality of free fall (also known as Weak Equivalence Principle (WEP) test). Thermal sources are defined, in this study, as atomic ensembles with a vanishing condensed fraction. Conversely, BEC sources possess a condensed fraction of 100%. For each case, we limit the maximum allowed diameter of the atomic ensemble at the recombination pulse to preserve contrast. Subsequently, this enables the determination of the shot noise and other error terms for the trade-off.
This article is structured as follows. Starting with a brief overview of the state of the art of light-pulse atom interferometry, we continue with relevant error contributions to the read-out phases of atom interferometers, quantitatively evaluating them in the three study cases and -finally -discussing the limits of the condensed or thermal regime of the respective interferometry source.
The most prominent examples of atomic inertial sensors with thermal ensembles are gravimeters [32,33] that reach for example a sensitivity of 1.4×10 −8 g at 1 s [34] and gyroscopes, which measure rotations below 10 −10 rad/s in a few 100 s [18]. A compact prototype of a gravimeter using BECs reaches a sensitivity of δg/g = 3.7 × 10 −6 per cycle [35]. This suggests the possibility of a targeted stability on the order of 7.8 × 10 −10 after 100 s of integration time with a state of the art BEC source [11]. In a proofof-principle experiment [36], a point-like BEC source for atom interferometers is implemented in a large fountain experiment to achieve sensitivities of 6.7 × 10 −12 g.
The gravitational constant G is determined to a value G = 6.67191(99)×10 −11 m 3 kg −1 s −2 with a relative uncertainty of 150 ppm [22] limited by the initial velocity spread of the atoms in the interferometer.
Most recently [21,37,38], there has been extensive work on the determination of the fine-structure constant α via determination of the ratioh/m with matter-wave interferometry, where m is the atomic mass andh is Planck's reduced constant. In a fountain with thermal cesium atoms, the fine-structure constant is determined with an expected statistical error of 0.008 ppb [37]. The systematics are at the 0.12 ppb level, mainly stemming from spurious accelerations. With thermal rubidium,h/m is measured at the 5 × 10 −9 level [38]. An ytterbium contrast interferometer with BECs [39] is used to demonstrate anh/mmeasurement using large momentum transfer, controlling diffraction effects and atomic interactions with suppression of vibrational effects allowing sub-ppb precision.
In [26], a dual-species WEP test with 85 Rb and 87 Rb reaches a statistical uncertainty of η = 0.8 × 10 −8 and is limited by systematic effects, e.g. the Coriolis effect to η = (2.8 ± 3.0) × 10 −8 . Most recently, this limit has been pushed further down to η = 1.6 ± 5.2 × 10 −12 [3]. The STE-QUEST mission [40] aims at testing the WEP at the 10 −15 level and beyond by measuring the differential acceleration of a 87 Rb BEC and a 41 K BEC over a total mission time of 5 years [41]. The concept Quantum Test of the Equivalence principle and Space Time (QTEST) [42] intends to determine η at the 10 −15 level with thermal ensembles over four integration periods of three months each aboard the International Space Station.

Performance indicators
In the previous subsection, the state-of-the-art sensitivities for measurements of rotations, accelerations, the finestructure constant and the Eötvös ratio have been stated. The phase that is to be determined depends on several experimental parameters like the effective wave vector k eff , the interrogation time 2T, the velocity v of the atomic ensemble perpendicular to the sensitive axis, the length of the detector baseline L or D and the frequency f of the gravitational wave. For the commonly-proposed interferometry schemes discussed above, one can summarize the performance-defining scaling factors to be: k eff T 2 for gravimetry, WEP tests and G measurements, -k 2 eff T forh/m measurements, -k eff L cos(f T ) for gravitational wave detection, -k eff DT 2 for gravity gradiometry, -k eff T 2 v for rotations.
Increasing these scaling factors allows an improvement in sensitivity.
In this paper we consider shot-noise-limited measurements, where the single-shot phase uncertainty is given by defined by the number of interfering atoms N at and interferometric contrast C. The experiments are repeated n cycle times to average ('integrate') the noise of a singleshot phase.
With the assumption of a shot-noise-limited measurement, a fixed atom number and no reduction in contrast C, an increase in the scale factor by enhancing the free evolution time T or the effective wave number k eff can increase the single-shot phase sensitivity. Interrogation times 2T of several seconds were realized [36] and an extension to 10 s was proposed on space platforms [40]. The effective momentum transfer ranges from few 10 k eff for a single multiphoton pulse up to a few 100 s of k eff [43] for benchmark experiments. The integration time to reach the desired performance may range from typically 10 4 s up to several months. Generally, the large number of atoms in thermal ensembles is an advantage over BECs to reduce shot noise. Their increased spatial extension is expected to suppress mean-field effects efficiently compared to BECs. On the other hand, the exact same position and velocity distributions might limit the scaling factors that one could achieve due to large systematic uncertainties and atom losses. In the next section, these systematics and other potentially performance-limiting effects are quantified.

Performance limiting effects in atom interferometry
Apart from shot noise considerations, a variety of physical phenomena constitute limiting factors for precision experiments by coupling to the velocity spread or the spatial extension of the ensemble, as is the case for the Coriolis effect, GGs, WFA or mean-field effects. In the following, we characterize different systematic and statistical effects that might limit near-future experiments beyond state of the art such as long-fountain atomic gravimeters, spacebased atom interferometers and atom interferometers operated in ground-based laboratories in micro-gravity environments. Starting with intrinsic loss mechanisms due to matter-light interaction, we go on to discuss DKC as a technique to suppress systematic effects. An analysis of the effects of imperfect detection of the atomic sample concludes this section.

Coherent manipulation
The fidelity of the interferometric beam splitters and mirrors realized by the coherent manipulation of the atoms using light is closely connected to the phase-space properties of the atomic ensemble.
First, homogeneous excitation of the atomic ensemble requires a constant Rabi frequency over the spatial extent of the atoms, which in turn implies a laser beam size much larger than the ensemble size. In cases of optical power constraints, e.g. typical for space missions, the ensemble size is hence restricted in order to maintain contrast by achieving reasonable rates for coherent manipulation. For large free-fall times such a requirement can be translated into a maximum expansion rate of the ensemble. Figure 1 shows the significant difference in expansion rate between thermal and condensed ensembles which indicates a clear advantage of the latter especially for large interferometry time scales. Second, for all applications outlined in the previous section, typically Doppler-sensitive two-photon couplings are employed, such that the longitudinal atomic velocity introduces a detuning. The temporal profile of light pulses determines their velocity acceptance and hence defines criteria for the atoms' velocity dispersion. In general, smaller velocity widths, corresponding to a smaller distribution of detunings, are desirable for an efficient addressing of the atoms.
An effective, simplified model [44] allows to assess the effects of spatial and velocity selectivity quantitatively. The convolution of the atoms' radial density distribution n(r, t) and longitudinal dispersion f (v) with the positionand velocity-dependent excitation rate of the pulse determines the total excitation rate Note, that the transverse velocity distribution modifies the time-dependent spatial distribution n(r, t). This treatment assumes a box pulse of duration t and an effective Rabi frequency where the spatial dependence of the Rabi frequency Ω 0 (r) is given by the radial beam profile. A similar treatment was applied in the case of a recent gravitational wave detection proposal [4].

Wave-front aberrations
Matter-light interactions in the atom interferometric cycle are typically subject to the beam's natural wave-front curvature (e.g. of a gaussian beam) and additional imperfections of the laser beam profile [6,5] caused by non-ideal optics. While errors due to the initial collimation of large beams (>2 cm) are negligible, retro-reflection still introduces WFA that lead to a considerable systematic uncertainty. We employ a second order approximation to the deviation from flat wave-fronts for the combined effects of beam and optics. The resulting spatial dependence of the laser phase fronts imprints a position-dependent phase on the atoms. Depending on the amplitude and wavelength of the distortion relative to the size of the atomic ensemble, the resulting phase shift may average out, reduce contrast, lead to phase patterns that can be resolved during detection or result in an average phase shift. In state-of-the-art cold atom gravimeters, WFA induce the limiting systematic uncertainty of 30 nm/s 2 to 40 nm/s 2 [32,45]. A more recent analysis of the device in [45] evaluated the systematic error to 55 nm/s 2 with an uncertainty of 13 nm/s 2 [46]. We limit our study case to long-scale WFA, assuming a quadratic dependency of the wave-fronts on the transverse position of the atoms, as introduced by a curvature of the retro-reflecting mirror. In this case, the resulting wavefront curvature with radius R couples to the finite velocity spread σ v and induces the phase shift for a Mach-Zehnder (MZ) geometry, a spatial Gaussian density distribution of the ensemble and a Gaussian velocity spread σ v [5]. Here, k eff denotes the effective wave vector, T the interrogation time, k B is the Boltzmann constant and T at , m at refer to the ensemble temperature and atomic mass, respectively.

Mean-field effects
Mean-field effects arise due to atom-atom interactions in atomic ensembles, scale with growing densities and are an additional source for statistical errors. The mean-field energy reads and depends on the local density n(r) of the ensemble and the interaction strength g BEC int = 4πh 2 a sc /m at , where a sc is the s-wave scattering length. For a thermal ensemble, g thermal int = 2 g BEC int [47]. Following [48], the average mean-field energy for a spherical ensemble of volume V (t) = 4πr(t) 3 /3 with N at atoms is consequently given by E MF = g int N at /V . This assumes, however, a uniform density distribution with radius r while thermal and BEC ensembles in fact follow a Gaussian or parabolic distribution, respectively. Hence, we take the average of the mean-field energy by weighting it with the respective density distribution: In case of an equal g int , but unequal atomic density on the two arms of the atom interferometer, a spurious phase shift arises. Following [48], we model the contribution by linking the imbalance in density to the initial beam splitter and neglect effects due to overlap of the two arms or losses. If the initial beam splitter creates a superposition that deviates by σ N from equal probability in both states, the phase shift occurs, corresponding to the integral of the differential frequency shift between the arms. In our assessment, we assume the initial superposition to have equal probabilities of both states on average, but to be affected by white noise with a standard deviation of σ N per cycle. Without relying on quantum correlations [49], characterization of the beam splitter is limited by quantum projection noise, implying an upper limit of σ N = 1/ √ N at per cycle, which we adopt for our analysis. Consequently, the mean-fieldinduced phase uncertainty in our model depends on the atom number and implicitly on the time-dependent ensemble size, which allows a trade-off between maximal atom number fluctuation and minimum ensemble size at the first beam splitter.

Coriolis effect and gravity gradients
Two of the most relevant systematic effects are related to the Coriolis force and GGs. The first arises due to the transverse motion of the atoms with respect to the incident beam in combination with Earth's rotation, which forms an effective Sagnac interferometer [50,51]. The second is the acceleration uncertainty due to the mass distribution of Earth and the apparatus surrounding the experiment. Both give rise to systematic effects as they couple to the initial kinematic conditions of the ensemble [52].
In the case of a MZ geometry, the uncertainty in the atom's mean velocity δ v and mean position δ r couple to GGs γ parallel and perpendicular to the sensitive axis. The uncertainty in phase related to gravity gradients is given by Due to low order of magnitude of terms related to the coupling of rotations to transverse GGs, we restrict ourselves to phase uncertainties due to gradients parallel to the sensitive axis. Residual rotations Ω perpendicular to the sensitive axis couple to the atom's velocity. The phase uncertainty due to the Coriolis force (subscript C), is given by As the density distribution of thermal ensembles is governed by the statistics of a Gaussian density distribution, the uncertainties of the mean position δ r and mean velocity δ v are related to the spatial σ r and velocity σ v spread, respectively, via The number of atoms N at and the number of pre-requisite measurements ν 0 (i.e. the number of times the ensemble is imaged before the actual experiment) equally contribute to the statistical repetition. BECs follow a parabolic density distribution. However, this can be approximated by a Gaussian distribution via Equation 23, as explained in the next section.

Expansion rate and collimation
Since expansion rates have sizable effects on the precision and accuracy of atom interferometers, we discuss the possibilities offered by DKC to reduce them in this section. DKC is an established tool to further reduce the effective expansion rate of an ensemble (see [53] and references therein). This phase-space-manipulation technique exploits that the free evolution of a gas released from a trap leads to a linear correlation of momentum and radial position of atoms within the ensemble. Recapturing the ensemble in a quadratic potential for a well-chosen duration leads to its collimation. Preservation of phase space density in the collisionless case requires that this reduction in momentum spread is accompanied by an increased ensemble size and hence necessitates a trade-off between desired expansion rate reduction and required growth in ensemble size. For non-interacting gases, i.e. thermal or fermionic gases in all regimes, this relation is captured by Liouville's theorem, which states that the ratio of initial (σ v0 ) to final (σ v f ) velocity width is inversely proportional to the relative increase in ensemble size (σ r f /σ r0 ). A similar relation can be found for interacting degenerate gases -for which the asymptotic final expansion rate is determined by the initial localization through Heisenberg's uncertainty principle and a mean-field contribution [4,54,55] -and for interacting, non-degenerate gases [47,56]. Figure 1 illustrates a DKC sequence collimating a BEC and a thermal ensemble with typical parameters. The free expansion of the thermal ensemble is governed by the expansion law whereas the BEC dynamics are captured by corresponding scaling laws [54]. As an exemplary case, we take a thermal ensemble of 10 9 87 Rb atoms at 2 µK with a diameter of 2σ r = 0.2 mm and collimate it down to 80 nK, such that the required size at lens is 2σ r (t DKC ) = 1 mm. For a BEC, we assume an ensemble of 10 6 87 Rb atoms and collimate it to 50 pK in a trap with frequencies of 50 × 2π Hz. These are regime-typical parameters for experiments with either thermal ensembles or BECs [4,9,10]. Following [4], the expansion energies for the thermal atoms and the chemical potential µ BEC of the BEC are obtained from For a better comparison of the two fundamentally different density distributions, the Thomas-Fermi radius R TF of the isotropic BEC with parabolic density distribution can be related to a Gaussian spatial width σ r via [53] R TF (t) = σ r (t) √ 7.
The resulting characteristics of the collimation sequence for 87 Rb and 41 K ensembles are given in Table 1.  Figure 1: Size evolution of thermal ensembles (red) and BECs (blue) after release from a trap. A DKC stage is reducing the expansion energies down to 80 nK and 50 pK for the thermal and BEC ensembles, respectively (see Table 1 for the exact parameters). The dashed lines illustrate the expansion in the freely expanding case without collimation.
The illustrated collimation sequence assumes similar free expansion time t DKC prior to application of the lens for both regimes. In order to achieve a final expansion behaviour of the thermal ensemble similar to the BEC case, t DKC would need to be significantly increased (about two orders of magnitude, depending on the initial temperature), corresponding to a substantially larger ensemble size at the time of the lens [4]. This is a distinct disadvantage for thermal ensembles, since the DKC technique crucially depends on the harmonicity of the applied lens potential, which has to be verified over the entire spatial extent of the ensemble [11]. Application of velocityselective pulses for temperature reduction in 1D has the disadvantage of atom loss [12]. It is therefore not a promising pathway to reach expansion rates for thermal ensembles comparable to those of BECs. Raman sideband cooling [57,58,59] might be a better alternative for thermal ensembles even if it is limited to about an order of magnitude larger temperatures than what the BEC ensembles could reach.

Contrast and detection
Output states of atom interferometers can be detected by absorption or fluorescence imaging methods. Which method is appropriate depends on the experimental parameters and the information one wants to acquire. One main distinction is whether the relative atom numbers in the output port are counted or if atomic density distributions have to be spatially resolved. Atom number counting is commonly done with fluorescence imaging of the ensembles at the output ports, which relies on the excitation and successive emission of photons by the atoms that are then detected by a simple photo diode or CCD camera. The ensemble has to be excited by a laser beam, which means that it has to have a reasonably compact size to be illuminated. This can usually be achieved with thermal ensembles as well as BECs. The number of atoms contributing constructively to the signal at the output port is given by the product of the excitation probabilities at all interaction times t i : The contrast C is given by Equation 24 as the convoluted excitation probability for a given phase-space distribution of the ensemble. Inhomogeneous excitation efficiency or phase gradients e.g. caused by GGs may wash out the contrast [40].
In experiments employing spatial mapping of the output port wave function for the determination of the phase or analysis of features within the atomic ensemble [60], good spatial resolution along with a high signal-to-noise ratio and minimal systematic effects during detection are required. BECs with small spatial spread and expansion rates are thus favored over thermal ensembles to increase the spatial resolution of the CCD picture. Indeed, the high expansion rates of thermal ensembles may at long times lead to densities challenging for absorption imaging due to decreasing signal per volume.

Comparative performance study
Based on the discussion of phase shifts and performancelimiting effects in the previous section, we now elaborate on three study cases in which we compare quantum degenerate ensembles to thermal sources. In highly dynamical environments, such as inertial sensing and for navigation purposes, thermal sources may be beneficial since they typically feature more atoms and shorter cycle times, which decreases both, shot noise and integration time. However, a trade-off has to be found for every particular situation due to their relatively high expansion rates and spatial extension. The three specific examples selected for the comparison consist of inertial sensors that could operate beyond state-of-the-art in the near future: (i) A ground-based gravimeter with a relative uncertainty of ∆g/g = 10 −9 , (ii) a space gradiometer with a 2.5 mE resolution [61] and (iii) a WEP-test with an uncertainty of 2 × 10 −15 in the determination of the Eötvös ratio [40]. The interferometer geometries are illustrated in Figure 3.
In order to evaluate the performance of every regime, we will take typical parameters for thermal ensembles and BECs and assess their performance in each of the three experiments. Details of the ensemble sizes and velocity distributions for both regimes are given in Table 1 at key times of the interferometric sequences. We assume a Gaussian beam of 3 cm 1/e 2 -radius with Rabi-frequencies of π/(25 × 10 −6 ) Hz (gravimeter and WEP-test) and π/(55 × 10 −6 ) Hz (gradiometer) for second order beam splitting processes, GGs of 10 −6 s −2 and spurious rotations on the order of 1 µrad/s due to imperfect rotation compensation of the mirror and limited attitude control of the satellite for the two space-borne missions. The wave-front curvature is assumed to be 2.3 km (gravimeter), 5.6 km (gradiometer) and 250 km (WEP-test).
As illustrated in the previous sections, systematic effects linked to GGs and the Coriolis force are connected to the uncertainty of the mean position δ r and velocity δ v at the beginning of the interferometry sequence. The number of measurements ν 0 required for their characterization (see Equation 16) is determined for each application such that the largest systematic phase uncertainty related to GGs or rotations (Equation 11 and Equation 15) is below the target uncertainty chosen for the respective measurement.
The number of prerequisite experiments sufficient to suppress the systematic effects below the target uncertainty may differ for the BEC and the thermal case. In our cases, thermal ensembles have three orders of magnitude more atoms and are about 15-20 times larger than BECs after the DKC. Hence, the minimal number of characterization measurements ν 0 is reduced by a factor of 2 to 5 compared to the BEC case. For the sake of comparability, we choose to compute all systematic effects with ν BEC 0 . This enables an evaluation of the performance with a fixed set of parameters.
As the systematic phase uncertainty due to WFA depends on the velocity width and not the mean velocity (see Equation 5), its magnitude does not depend on the number verification measurements. Thus, it can neither be integrated down nor reduced by prerequisite measurements, conversely to the GGs. Its value is completely predetermined by the ensemble's expansion rate.
In order to adapt statistical error contributions such as shot noise and mean-field effects to the desired precision of every type of measurement, we calculate the minimum number of iterations n cycle until the target uncertainty is reached. The integrated (denoted by subscript i) shot noise is given by where the contrast C is given by Equation 24, i.e. the convolved excitation probability. Non-perfect contrast (C < 1) increases the shot noise as the number of atoms constituting the statistical sample is reduced. Mean-field effects contribute a statistical phase uncertainty expressed by Equation 10 where the ensemble expansion over time is taken into account. This effect integrates down with the number of experiments n cycle following Once the number of cycles is determined by the desired performance, the associated integration time t int depends on the preparation time t prep and the interrogation time 2T: t int = (t prep + 2T ) n cycle = t cycle n cycle .
For a straightforward comparison of the performance differences between BEC and thermal ensemble, the integration time is also chosen to be the same for both regimes, initially determined by the number of cycles needed to suppress the statistical effects of the BEC below the target uncertainty. Since thermal sources can be generated within a shorter preparation time, more cycles can be performed during the same integration time according to In order to estimate the various uncertainties, ensemble properties as spatial and velocity spreads are computed at each atom-light interaction pulse, to take into account the spatial and velocity selectivity of the pulses applying Equation 3 and Equation 24. The modified spatial and velocity spreads are the evaluation input for the mean-field effects, the WFA and the estimation of the contrast according to the formulae given in Section 2. The results of this study are summarized in Table 2 where the phase uncertainties are normalized as fractional phases φ/k eff gT 2 for the gravimeter and the WEP test. In the gradiometer case, the orders of magnitude are given in units of φ/k eff DT 2 = Γ for a baseline D. Here, φ target is the upper limit for any systematic or statistical phase uncertainty.
Detailed results for the three science cases are presented in the consecutive sections.

Gravimeter
We start by comparing two ground-based 87 Rb gravimeters operated with thermal atoms or BECs, the source characteristics of which are similar to the ones reported in [5] and [27,11]. In both cases, the interrogation time 2T equals 150 ms.
The first column in Table 2 shows the magnitudes of the different performance limiting effects discussed in the previous section in units of ∆g/g. The scenario targets an uncertainty of 1 µGal=10 −8 m/s 2 , corresponding to a fractional phase uncertainty of δ φtarget /k eff gT 2 = 10 −9 .
A thermal ensemble with 10 9 87 Rb atoms would have a diameter of 2σ r (t DKC ) = 1 mm after the DKC pulse. The velocity spread at the lens is 2.7 mm/s. The convolved excitation efficiency at the last beam splitter is 57%.
Adopting a BEC-source as in [11], it is reasonable to assume a collimation of the ensemble to an effective temperature of 50 pK for 10 6 atoms. The convolved excitation efficiency at the last beam splitter is 99% for an ensemble diameter of 2R TF (t DKC ) = 0.19 mm and a velocity spread of 183 µ m/s. Although the order of magnitude of the initial atom number is three times larger for the thermal atoms, the shot noise is only one order of magnitude below the one of the BEC due to the reduced contrast.
Theoretically, the target uncertainty of ∆g/g = 10 −9 is reached in both cases after only one verification shot and integration over seven (thermal, full cycle time 0.48 s) or three (BEC, full cycle time 1.15 s) experimental cycles corresponding to a few seconds of integration time. All statistical and systematic effects are below the target uncertainty hinting towards the possibility to use either of the source concepts.
However, if a better performance of the gravimeter is sought for, the first limit to tackle would be the meanfield effects at 9.2 × 10 −10 in the BEC case and wavefront distortions at 3.4 × 10 −10 in the thermal one. Meanfield effects are a statistical phenomenon, hence one can integrate the phase uncertainty for the BEC-case down by increasing the number of experiments.
As mentioned, WFA are not a statistical but a systematic phenomenon and cannot be integrated down by adding verification shots. Therefore thermal sources are limited by WFA to the 10 −10 level, whereas the BECs could improve on the accuracy up to the 10 −13 level.

Gradiometer
Here, we are address a satellite gradiometer as proposed in [61]. It features a baseline of D = 0.5 m separating the two interferometers, an interrogation time of 2T=10 s (full cycle time of 20 s) and a targeted uncertainty of 2.5 mE, clearly beyond the current state of the art. We adapt the interferometry time from 2T=10 s to 2T=0.5 s to constrain the ensemble size at a detectable level in the thermal case.
The center column of Table 2 displays the order of magnitude of uncertainties in the gradient determination related to the different effects. Here, the numbers are given as spurious gradients in units of Γ by dividing each phase uncertainty by k eff DT 2 .
The phase uncertainties due to GGs, the Coriolis force and mean-field effects receive contributions from both individual interferometers. Through Equation 16, the initial spatial and velocity spreads σ r,v;1,2 of interferometer 1,2 enter the systematic uncertainties given in Equation 11 and Equation 15 as supposing uncorrelated source noise. For initial spreads with the same width σ r,v;1 = σ r,v;2 , this yields a factor of √ 2, which results in an integration behaviour during the verification measurements of in the case of GGs and Coriolis effect. The same holds true for the integrated shot noise, which is increased by a factor of √ 2 compared to Equation 25: and for the mean-field effects, which are uncorrelated between the two branches of the interferometer: Interestingly, following our treatment in Subsection 2.3, we find that due to the different expansion behaviour and a higher atom number, the mean-field effects of the thermal ensemble are comparable, in magnitude, to that of the BEC on the time scales we are investigating (see Table 2 and Figure 2). Table 2: Estimation of statistical and systematic uncertainties for three scenarios: a lab-based 87 Rb gravimeter [5], a space-borne 87 Rb gradiometer [61] and a satellite 87 Rb/ 41 K WEP-test analogous to the STE-QUEST mission [40].
The phase uncertainties are given as fractions ∆a/g = δ φ /k eff gT 2 (gravimeter and WEP-test) and ∆Γ = δ φ /k eff DT 2 (space gradiometer). The expansion sequence over the course of the atom interferometer is calculated in Table 1. Systematic effects are denoted by δ, while statistical effects are denoted by σ and calculated after integrating over a number n cycle of experiments with N at atoms in each cycle. Gravity gradients are abbreviated with GG, the Coriolis effect with C, wave-front aberrations by WFA, shot noise with SN and mean-field effects by MF.  Assuming the same velocity spread for the two ensembles in the two interferometers and -as for the gravimeter case -the simplification of a constant curvature, the differential phase shift vanishes. Here, we drop this simplification and consider the propagation of a Gaussian beam which leads to a local dependency of the curvature, and thus to a non-vanishing phase shift in the differential signal. Assuming a residual radius of curvature of the retroreflection mirror of 4 km and a propagating the laser beam as outlined in [61] introduces the differential phase shift as reported in Table 2.
To reach the uncertainty goal of 2.5 mE, an integration time on the order of 1 day is required in both cases. Albeit the high expansion rate of the thermal ensemble is accounted for with a shorter interrogation time of 0.5 s (full cycle time of 1.2 s), the shot noise uncertainty in the thermal case exceeds the one of the BEC case as the contrast drops to 33% (97% for the BEC after 10 s). Again, the WFA would hinder any further attempts for a better performance below the 10 −12 s −2 level. All other systematic effects are complying with the performance required from this sensor and can be reduced by an increased number of verification shots.
Being limited by WFA at the level of 4.4 × 10 −15 , the BEC clearly leaves more room for improvement compared to thermal sources, which are bound three orders of magnitude above. Therefore, using a BEC source in this scenario is advantageous since it does neither suffer from a worse integration time nor a larger mean-field effect, yet significantly extends the desired sensitivity range compared to the thermal source.

WEP-test
The WEP-test example is based on the parameters of the STE-QUEST satellite mission proposed in [40]. We here assume the test pair 41 K and 87 Rb. It aims at an Eötvös ratio η determined with an uncertainty of 2 × 10 −15 . As for the gradiometer, we compare a thermal ensemble at 0.5 s of interrogation time (full cycle time of 0.83 s) with a BEC at 10 s (full cycle time of 20 s) of interrogation time to ensure non-vanishing contrast and technically manageable ensemble sizes in both cases. In the last column of Table 2, the results of the comparison are displayed.
For the determination of the systematic effects one has to take into account the differences of the two atomic species such as the different expansion rates, atomic masses and initial conditions like spatial spread. The two different species-specific excitation rates are also evaluated and the minimum given in Table 2.
The effects of GGs and Coriolis force do not scale with √ 2 in this case, but rather with the mean square of two uncorrelated spreads (see Equation 29). With two different species propagating in the arms, the mean-field effects are calculated as the mean-square sum of the individual mean-field effects of 87 Rb and 41 K as in Equation 10. As for the gradiometer, one benefits from the differential suppression of WFA when matching the expansion rates of the ensembles [40]. As two different sources and two different species are used for the production of the ensembles, the systematic differential phase uncertainty is given by  Table 2. (a) Time-dependent density ρ of the ensembles during the interferometer time based on the DKC sequence described in Table 1. (b) Fractional statistical phase uncertainty O(σ φ MF ) due to mean-field effects according to equation Equation 10 integrated over a number of n cycle experiments.
analogous to Equation 5. By experimentally matching the velocity spreads of the two ensembles to the 20% level in both arms, the WFA are suppressed by a factor of 3.
In the BEC case and in order for the statistical effects to be consistent with the mission performance goal, 10 6 experimental cycles n cycle are needed. This leads to a full mission time on the order of five years in case of a highly elliptical orbit operation. To reach the targeted uncertainty, an additional 10 6 verification measurements ν 0 are necessary. These measurements are included in the total mission time as they are performed on the transition between perigee and apogee parts of the orbit not dedicated to science measurements [40]. The contrast at the end of the sequence is at 99% for both, 87 Rb and 41 K, and for the chosen set of parameters it is feasible to achieve the missions goals with condensed ensembles.
With a thermal ensemble, one would need an integration time on the order of 10 8 s to suppress mean-field effects, and 10 6 verification shots to estimate the phase uncertainty due to GGs and the Coriolis effect to a sufficient level. Moreover, even for the reduced interrogation time of 0.5 s, the contrast is at 43% and the shot noise is almost 6 times larger than in the BEC case.
We conclude that -for the chosen set of parametersit is possible to achieve the mission goals with condensed π/2 π π/2 Figure 3: The geometry for a gradiometer consists of two MZ sequences that are operated simultaneously (black and blue lines). Each MZ sequence is formed by three consecutive light pulses that constitute a beam splitter (first π/2-pulse), a mirror (π-pulse) and a merging beam splitter after 2T . The dashed lines indicate a change of the internal state due to a momentum transfer by the atomlight interaction. A single MZ sequence corresponds to a gravimeter sequence to measure the coupling constant g of the gravitational field to one species. Two MZ sequences operated with two different species allow for a determination of η. In this case the differential velocity at t = t start of the MZ sequences is set such that they spatially overlap to eliminate the gravity gradient, but remain sensitive to a potentially species specific differential acceleration.
ensembles. The thermal case is, however, limited to the 10 −12 level due to the WFA effects stemming from the large expansion rates of the thermal ensembles [46]. In order to reach a reduced velocity spread with a thermal ensemble comparable to that of a BEC, the DKC would be required to handle ensembles with a diameter on the order of 0.5 m after a pre-lens free expansion time of several seconds. As for the space-gradiometer, this is unpractical for obvious manipulation and excitation reasons. In the condensed regime, a simultaneous collimation of the dualsource was recently considered in reference [62] and shown to be feasible.

Discussion and conclusion
In this paper, the current limits for state-of-the-art precision experiments with atom interferometry were analyzed. A particular emphasis was put on the comparison of the statistical and systematic uncertainties between condensed and thermal ensembles. Three detailed study cases of a lab-based gravimeter, a space gradiometer and a satellite WEP-test were chosen to illustrate the limits of each regime.
Thermal sources benefit from a shorter cycle time and larger atom numbers compatible with experiments where moderate scale factors suffice or rapid readouts are required. This is, however, beneficial at short interferometry times only. When going beyond state-of-the-art, i.e. from drift times of a fraction of a second to a few seconds, this advantage is lost. In our study cases, the scenarios utilizing BECs show -at equal integration time -the same shotnoise level and the magnitude of the mean-field effects is comparable to that of thermal sources. Moreover, the condensed sources benefit from a very large contrast (close to 1) when compared to their thermal counterparts. More dramatically, the WFA set an ultimate limit for thermal ensembles that would not be compatible with long interrogation times, which are required for advanced scenarios. For BEC ensembles, their compact sizes make this limit at least three orders of magnitude lower, highlighting their potential in the field of metrology.
Small scale distortions (few µm) of the optical beams [63], not considered in this article, can hint to a disadvantage for the BEC samples by means of averaging effects for WFA. However, the flexibility in tuning their initial size [62] mitigates this effect and could bring them to starting sizes similar to thermal ensembles if necessary. This engineering of the BEC size allows for a distinct analysis of WFA with long and short periodicity [5,46]. This might be especially relevant on short time scales, i.e. for very small ensemble sizes. Their subsequent expansion could still be limited to a few mm thanks to the DKC technique. In consequence, a trade-off between the size-stretch-induced phase uncertainties, e.g. to balance the level of GGs or Coriolis systematics versus WFA effects is required. This appears to be feasible, especially if one considers gravity gradient compensation schemes as the one in [64,61,65].
Other considerations that are not reflected by our study would further consolidate the BEC choice. Indeed, we optimistically anticipate here that thermal ensembles could be collimated to the 80 nK level and that the same level of efficiency in preparing, transporting and engineering of their the quantum states can be achieved as for BECs. As a conclusion, thermal and BEC sources could equally be employed in relatively short interferometry times (a few hundred ms) for the same performance. With respect to longer times, BEC sources are clearly more advantageous since size-related systematic effects are several orders of magnitude smaller than those of thermal ensembles.

Authors contributions
All the authors were involved in the preparation of the manuscript. All the authors have read and approved the final manuscript.