R Chang

Q Bouton

H Cayla

C Qu

A Aspect

C I Westbrook

D Clément

Momentum-Resolved Observation of Thermal and Quantum Depletion in a Bose Gas

We report on the single-atom-resolved measurement of the distribution of momenta hk in a weaklyinteracting Bose gas after a 330 ms time-of-flight. We investigate it for various temperatures and clearly separate two contributions to the depletion of the condensate by their k-dependence. The first one is the thermal depletion. The second contribution falls off as k -4 , and its magnitude increases with the in-trap condensate density as predicted by the Bogoliubov theory at zero temperature. These observations suggest associating it with the quantum depletion. How this contribution can survive the expansion of the released interacting condensate is an intriguing open question.

In quantum systems, intriguing many-body phenomena emerge from the interplay between quantum fluctuations and interactions. Quantum depletion is an emblematic example of such an effect, occurring in one of the simplest many-body systems: a gas of interacting bosons at zero temperature. In the absence of interactions, the ground state corresponds to all particles occupying the same single-particle state. Taking into account inter-particle repulsive interactions at the meanfield level leads to a similar solution where all particles are condensed in the same one-particle state whose shape is determined by the trapping potential and interactions. In a beyond mean-field approach, which can be interpreted as taking into account quantum fluctuations and twobody interactions, the description is dramatically different. The many-body ground state consists of several components: a macroscopically occupied single-particle state, the condensate, and a population of single-particle states different from the condensate, the depletion. This many-body description applies to diverse bosonic systems such as superfluid Helium [START_REF] Sokol | Bose-Einstein condensation[END_REF], exciton-polaritons [START_REF] Utsunomiya | [END_REF] and degenerate Bose gases [3]; it has also found analogies in phenomena such as Hawking radiation from a black-hole [4] and spontaneous parametric down conversion in optics [5]. The fraction of atoms not in the condensate at zero temperature, the quantum depletion, increases with the strength of inter-particle interactions and with the density, rising up to 90% in liquid 4 He [START_REF] Sokol | Bose-Einstein condensation[END_REF]. In ultracold gases, where the density is significantly smaller, the quantum depletion usually represents a small fraction (less than 1%) of the total population. At non-zero temperature there is an additional contribution to the population of single-particle states above the condensate, originating from the presence of thermal fluctuations.

For weakly interacting systems, Bogoliubov theory describes quantum and thermal contributions to the condensate depletion [6,7]. This approach shows that the elementary, low-energy excitations are collective quasiparticle (phonon) modes, as has been verified in experimental studies with liquid 4 He [8], degenerate quantum gases [9] and exciton-polaritons [START_REF] Utsunomiya | [END_REF]. At zero temperature, the many-body ground state is defined as a vacuum of these quasi-particle modes. When projected onto a basis of single-particle states with momentum hk, this manybody ground state exhibits a distribution n(k), which scales as k -4 at large k. These power law tails do not exist in mean-field descriptions, for which the momentum distribution has a finite extent. At non zero temperature, the contribution to n(k) induced by thermal fluctuations decays exponentially for energies larger than the temperature. Previous experiments with atomic gases [10,11] have observed the total depletion of the condensate after a time-of-flight expansion, but could not distinguish between the thermal and quantum contributions.

In this letter, we report on the observation of momentum-space signatures of thermal and quantum depletion in a gas of interacting bosons. We investigate, for various temperatures and atomic densities, the three-dimensional atomic distribution after a long timeof-flight (see Fig. 1), i.e., the asymptotic momentum distribution. Three components can be identified (see Fig. 2): the condensate (I), the thermal depletion (II) and a tail decaying as k -4 and increasing with the intrap condensate density (III). This suggests associating region III with the quantum depletion, but with two caveats. Firstly, k -4 -tails originated from contact interaction were observed to vanish during the expansion of interacting fermions [12]. Recent theoretical work investigating interacting bosons predicts that the k -4 -tails adiabatically decrease with the condensate density during the expansion [13]. Secondly, the magnitude of the k -4 -tails we measure is larger than the in-trap prediction of the Bogoliubov theory, by a factor of about 6. Our identification of the k -4 -tail with the quantum depletion thus demands that there exists either a non-adiabatic process in the expansion, decoupling the in-trap k -4 component from the expanding condensate, or an interactioninduced effect beyond the mean-field description of the expansion, leading to 1/k 4 tails.

Our experiment is performed with a Bose-Einstein condensate of metastable Helium-4 atoms ( 4 He * ). Cigarshaped condensates, of typically N = 2 × 10 5 4 He * atoms in the polarized 2 3 S 1 , m J = +1 state are produced in an optical dipole trap with trapping frequencies ω/2π = (438, 419, 89) Hz [14,15]. After abruptly turning off the optical trap (in less than 2 µs) we detect the gas with a micro-channel plate (MCP) [16,17], after a 55 cm free-fall corresponding to a time-of-flight (TOF) of ∼ 330 ms. A radio-frequency (RF) ramp transfers a fraction of the polarized m J = +1 atoms to the nonmagnetic m J = 0 state after 2 ms of TOF (see [18] for details). The presence of a magnetic gradient after the RF ramp ensures that only m J = 0 atoms fall onto the MCP and can be detected. The MCP allows us to detect 4 He * atoms individually and to record the twodimensional (2D) position and the arrival time of each atom in the plane of the detector. The arrival time of each atom directly translates into a vertical coordinate, so that we reconstruct the complete 3D atoms distribution, in contrast with usual optical imaging which yields a 2D column-integrated density. Another advantage of a MCP operated in counting mode is its extremely low dark count rate. Here it allows us to investigate the atomic density over more than 5 decades (see below).

The detector is placed 55 cm below the trapped cloud (see Fig. 1(a)), so that after switching off the trap, the observation is made in the far-field regime where finite-size effects of the source can be safely neglected. In the freefalling frame of reference, we thus identify the position r of a detected atom (with respect to the cloud center) with a momentum component k = mr/h t, where t = 330 ms is the time-of-flight of the cloud center [18]. This yields the asymptotic momentum distribution n ∞ (k) obtained from the density distribution of the expanding cloud n(r, t), according to the ballistic relationship

n ∞ (k) = (h t/m) 3 n(r, t) (1)
The distribution n ∞ (k) should not be confused with the in-trap momentum distribution n(k), since the initial phase of the expansion is affected by inter-atomic interactions. Nevertheless, as we argue below, the highmomentum tails of n ∞ (k) provide interesting information on the in-trap momentum distribution n(k).

An example of n ∞ (k) is shown in Fig. 1(b). The main component is the pancake-shaped distribution expected for the cigar-shaped condensate according to the meanfield description of the expansion [19,20]. We also distinguish momentum components beyond those of the condensate, with a much lower density and an isotropic distribution [18]. From 3D distributions n ∞ (k), we obtain radial and longitudinal profiles as shown in Fig. 2 [18]. The resulting profiles exhibit three distinct regions, as illustrated in Fig. 2(b).

Firstly, the observed distributions are dominated by the high-density condensate (region I: k ≡ |k| < 1.7 µm -1 ). The initial expansion of the condensate is driven by inter-particle interactions, resulting in an asymptotic distribution different from the in-trap momentum distribution. This dynamics is fully captured by a meanfield treatment, the scaling solution [19,20] calculated in the Thomas Fermi approximation, as shown in Fig. 2

(b).

We have checked that a numerical solution of the timedependent Gross-Pitaevskii equation for the pure condensate, beyond the Thomas Fermi approximation, yields negligible corrections [18].

Beyond the anisotropic mean-field distribution we observe high-momentum tails (k > 1.7 µm -1 ) with spherical symmetry [18]: we attribute these components (regions II and III in Fig. 2(b)) to the depletion of the condensate. The isotropic character of the atomic distribution in the tails indicates that the anisotropic mean-field potential describing the interactions in the condensate does not play a significant role in the expansion of particles belonging to regions II and III. Thus we assume in the following that the high-momentum components of the tails, corresponding to a kinetic energy larger than the chemical potential of the condensate, quickly escape the condensate and are not affected by the mean-field potential during the expansion. The tails are visible over three decades in density, allowing us to perform a detailed study of their momentum dependence. We observe two distinct regions: a middle region without a well-defined power-law variation (region II), and a high-momentum region with density varying as k -α (region III). A powerlaw fit to the data in region III yields α = 4.2 [START_REF] Utsunomiya | [END_REF].

A quantitative description of the condensate depletion close to zero temperature is provided by Bogoliubov's microscopic theory, which yields a beyond mean-field model taking into account quantum and thermal excitations [6,7]. In the many-body Hamiltonian, this approximated ). Solid lines are a smoothed average of the density data, and the lightly shaded band is the running standard deviation. The dark count rate corresponds to a level less than ∼ 10 -2 µm 3 , which is one order of magnitude below the lowest data point. The plotted range of momenta is limited by the physical size of the detector.

approach retains only quadratic terms in the particle operators a k , where a k is the operator annihilating a particle with momentum hk. The simplified Hamiltonian can then be diagonalized by introducing the quasi-particle operators b k defined by the Bogoliubov transformation [6,7]. At zero temperature, the many-body ground state is the vacuum of quasi-particles, defined as b † k b k = 0 for any k = 0. It corresponds to a non-zero population of excited single-particle states, a † k a k = |v k | 2 for k = 0. This is the quantum depletion of the Bose condensate, which has no classical analog, and can be seen as arising from the interplay of Bose statistics and interactions. At non-zero temperature, the particle occupation number of non-zero momentum k can be expressed analytically in terms of the quasi-particle occupation number:

b k = u k a k + v -k a † -k
n(k) = a † k a k = (|u k | 2 + |v k | 2 ) b † k b k + |v k | 2 (2)
with the occupation number of quasi-particles given by a Bose-Einstein distribution,

b † k b k = 1 exp( (k)/k B T ) -1 (3)
where (k) is the Bogoliubov dispersion relation, k B is Boltzmann's constant and T is the temperature. The interpretation of Eq. 2 is clear: the first term, proportional to b † k b k , represents the thermal depletion; the second term is associated with the quantum depletion.

The particle occupation number |v k | 2 corresponding to the quantum depletion varies as k -1 for kξ 1, and as k -4 for kξ 1, where ξ is the healing length of the condensate. The small k behavior is related to the Heisenberg inequality associated with the particle and the density operators [21]. The large k behavior arises due to the two-body contact interaction and is related to Tan's contact constant, a universal quantity that connects contact interactions to the thermodynamics of a many-body system [22,23]. By contrast, the depletion associated with the thermal excitations varies differently with k due to the additional term b † k b k . In particular, it decays exponentially for kλ dB 1, where λ dB = h 2 /2πmk B T is the de Broglie wavelength. These differences provide a means to unambiguously distinguish the quantum depletion from the thermal depletion.

The presence of an inhomogeneous trap does not modify the prediction for the condensate depletion at momenta large compared to the inverse system size 1/R 0.08 µm -1 , where the results of Bogoliubov theory can be averaged using the local density approximation (LDA). For a harmonically trapped gas, a LDA calculation keeps the large-momentum k -4 scaling of the homogeneous model [18]. On the other hand, the thermal depletion distribution in a harmonic trap approaches a polylog function at high temperature [3].

To identify the contributions in regions II and III, we have investigated the tails of the measured distribution n ∞ (k) as a function of temperature. Fig. 3 presents the radial distributions k 4 × n ∞ (k) for clouds subjected to a controlled heating sequence [18]. Assuming the LDA average of Eq. 2, we fit the tails in Fig. 3 (k > 2µm -1 ) with the function k 4 × n fit (k) where

n fit (k) = N th g 3/2 [exp(-k 2 λ 2 dB /4π)] 1.202 (2π/λ dB ) 3 + C ∞ (2π) 3 k 4 (4)
with T a (via λ dB = h/ √ 2πmk B T a ), N th and C ∞ are fitting parameters. The first term in Eq. 4 is the polylog function describing a thermal component with an atom number N th and an apparent temperature T a [24]. The second term corresponds to a distribution decaying as 1/k 4 . The function n fit (k) is an excellent fit to the experimental profiles (see Fig. 3). As the gas is heated, the

T a = 149(2) nK k 4 n (k) [m -1 ] f th = 0.13 f th = 0.07 f th = 0.11 f th = 0.28 = 257(3) nK ∞ k 4 n (k) [m -1 ] ∞ k 4 n (k) [m -1 ] ∞ k 4 n (k) [m -1 ] ∞ FIG. 3. Measured k 4
plotted as a function of k and at various temperatures. Solid blue lines are a smoothed average of the density data, and the lightly shaded blue band shows the running standard deviation. The dashed line is a fit using Eq. 4 which involves two terms: the thermal depletion (also shown as a dotted line) and the quantum depletion revealed by the k -4 scaling at large k. The fitting procedure yields the temperature Ta, the thermal atom number N th and the asymptotic constant C∞. Noted in each subplot are Ta, the ratio of the thermal energy and the chemical potential kBTa/µ and f th = N th /N .

temperature T a and the thermal fraction f th = N th /N increase. The variation of f th with T a /T c (T c being the critical temperature of condensation) is in excellent agreement with the semi-classical prediction [3], confirming our identification of region II with the thermal depletion.

Although they represent less than ∼ 0.5 % of the total atom number, the k -4 tails are clearly visible beyond the thermal component (see Fig. 3), and thus associated with a zero-temperature effect [START_REF]Contact interactions between thermal atoms lead to a thermal contribution to the k -4 -tails[END_REF]. In the weakly interacting regime we investigate, condensate lifetimes are on the order of seconds. We have measured identical k -4 -tails when holding the atoms for an extra second in the trap, showing that the gas is at equilibrium before the release.

The presence of k -4 -tails in a cloud released from a trap was previously reported in strongly interacting Fermi gases [12] but was not found with bosons [START_REF] Makotyn | [END_REF]. The observation of the k -4 -tails in a Fermi gas required ramping the interaction strength to zero before the expansion [12], on a time scale shorter than that associated with many-body effects. Recent theoretical work concluded that the k -4 -tails should adiabatically vanish during the expansion of a Bose gas when the strength of interaction is kept constant [13]. These considerations indicate that in order to associate the observed k -4 -tails in n ∞ (k) with the quantum depletion in the trapped cloud, we must invoke a non-adiabatic process. Since the scattering length in the m J = 0 state is expected to be smaller than that in the m J = +1 state [18, 27], a non-adiabatic transfer between these two states at the optical trap turnoff might explain our observation, but we have not yet found any evidence of this possibility in the experiment. On the other hand, there is no many-body treatment of the expansion of an interacting Bose gas. We thus cannot exclude the possibility that the tails result from a modification of the beyond mean-field momentum distribution during the time-of-flight dynamics.

In order to further investigate the origin of the k -4tails, we have studied their dependence upon the condensate density. The fitting parameter C ∞ of Eq. 4 is equal to the Tan contact constant [22,23], which, for a harmonically trapped gas, is found equal to C LDA = (64π 2 /7)a 2 s N 0 n 0 in the LDA approximation [18]. The experimental results are plotted in Fig. 4 where the error bars reflect the uncertainty on C ∞ from the fit, as well as those on the calibration of N 0 and n 0 [18]. The fitted contact constant C ∞ /N 0 per condensed particle is found proportional to n 0 , as expected. The measured values of C ∞ , however, are about 6.5 times larger than the expected value C LDA . Note that in order to increase the density n 0 we increase the trapping frequency, which results in a decrease of the density of the central, dominant part of the distribution n ∞ (k) measured after TOF. The observed proportionality between C ∞ /N 0 and n 0 in spite of the variation of n ∞ (k 0) rules out several possible spurious effects in the response of the MCP.

In conclusion, the measurement of the momentum distribution of a weakly interacting Bose gas released from a trap has allowed us to observe two components in the high momentum tails beyond the mean-field distribution. The first one is due to thermal depletion, and although some questions remain open, there are several observations which suggest associating the second one with the quantum depletion. The single-atom detection method of metastable Helium gases is also able to provide signals of atom-atom correlations in momentum, a feature we intend to use in future investigations of momentum-space signatures of many-body effects. Our condensates are produced in the 2 3 S 1 , m J = +1 internal electronic state in a crossed dipole trap which is described in [15]. This metastable state has a lifetime of 8000 seconds. To measure the asymptotic momentum distribution n ∞ (k) of the gas, we switch off the optical trap in 2 µs, allowing the cloud to undergo expansion while falling under the influence of gravity. The Micro-Channel Plate (MCP) detector and electronics are described in detail in [16]. The estimated detection efficiency is 25%. The spatial resolution in the plane of the detector is measured to be equal to 100 µm and the time resolution along the vertical axis is 10 µs.

C /N [µm ] ∞ n [µm ]
The electronic detector is positioned 55 cm below the trapped gas giving a TOF value of 330 ms. To ensure that the gas expansion is un-perturbed by residual magnetic field gradients present during the TOF, atoms are transferred to the magnetically insensitive m J = 0 state. This transfer is achieved with a radio-frequency (RF) which causes transitions between the magnetic sublevels whose energies have been split by a magnetic bias field with ∆E = E m J =+1 -E m J =0 h × 10 MHz. To transfer the atoms independently of their velocity, we use a RF sweep with central frequency 10 MHz and span of ±500 kHz, applied on the atoms after a 2 ms TOF. The RF sweep is 1 ms long at constant RF power. After the RF sweep, we apply a magnetic gradient to push the atoms remaining in the m J = ±1 states away from the detector. The RF power, combined with the removal of the m J = ±1 states, allows us to control the flux of atoms (m J = 0) striking the detector. We typically operate between 15 and 45% RF transfer efficiency (see 1D density profiles). The scattering length in the m J = +1 is a m J =+1 s 142a 0 , where a 0 is the Borh radius. In the m J = 0 state, it has never been measured so we have inferred its value from the predictions of [27] with the knowledge of a m J =+1 s to extract the contribution of the quintet potential. We find a m J =0 s ∼ 100a 0 , much smaller than a m J =+1 s . Note that an experimental measurement might reveal corrections to this approximate value. The MCP detector provides a three-dimensional histogram of atom numbers (see Fig. 1(b)). The position of an atom labelled with integer j is given by two-dimensional spatial coordinates (Y j , Z j ) in the plane of the MCP and the time of arrival t j . Similarly we note (Y 0 , Z 0 , t = t 0 ) the coordinates of the center of the cloud whose time of arrival defines the TOF t used throughout this work. In the frame centered on the falling could, the position r of the atom is r = (g t/2t j × ( t2t 2 j ), Y j -Y 0 , Z j -Z 0 ), accounting for the acceleration g of gravity. The use of the ballistic relation yields the asymptotic momentum distribution n ∞ (k) = (h t/m) 3 n(r = hk t/m, t) (Eq. 1 in the main text). The momentum-space resolution is 0.03 µm -1 along directions y and z (in the MCP plane), and 0.01 µm -1 along x (orthogonal to the MCP plane). The distributions studied here consist of the the average of roughly 1500 experimental shots. Individual shots have been re-centered to account for slight fluctuations in cloud center-of-mass after the TOF.

The condensate atom number N 0 and density n 0 are calibrated by comparing the 3D density profiles on the MCP detector with the predictions of the scaling solution for our trap frequencies [19,20]. The uncertainty on these values is evaluated to be 20%.

1D density profiles

From the three-dimensional asymptotic momentum density n ∞ (k), one-dimensional profiles are generated for the radial and longitudinal cloud directions. Longitudinal profiles are cuts along the z direction (at k x k y 0) with a small integration along the two transverse x and y directions (see below). Radial profiles in the xy plane are generated from an angular average over azimuthal angles |φ| > 45 • , where φ = arctan(k x /k y ) (at k z 0). This choice of azimuthal angles avoids a defect located on the surface of the micro-channel plate.

Since the experimental signals contain both a high-density condensate and low-density tails with a 4 decade separation in scales, it is necessary to divide the measurement of the density profile into two steps. To measure the condensate momentum components we first use a low RF transfer efficiency (typically 15%). In this low particle flux regime the MCP is far from electronic saturation, but the low-density tails are hard to detect. Secondly to measure the high-momentum tails, we use a high RF transfer efficiency (typically 45%). At high flux of detected particles, the dense condensate locally saturates the MCP, while the low-density wings remain unperturbed. We have verified that the two runs provide identical density profiles at intermediate k where local saturation is not a problem for the high-flux runs and there is sufficient signal in the low-flux runs. The 1D profile for the condensate (low-flux data) has a transverse integration of ± 0.1 µm -1 . The profile for the tails (high-flux data) has transverse integration of ± 0.8 µm -1 . The two-step measurement and transverse integration ensures sufficient signal in the high-momentum tails, while accurately capturing the condensate profile at low-momenta.

In order to illustrate the observed symmetry in region II and III, we plot measured 1D profiles with a small angular average on φ (±10 • degrees) along different directions separated by 30 • degrees in Fig. S1. The radial 1D profile shown in Fig. 2 (obtained from a large angular average) is reported as a black line in each subplot. The anisotropy of the condensate distribution appears clearly while the momentum profiles in region II and III have a spherical symmetry.

Mean-field momentum distribution of an interacting gas from time-of-flight

The condensate and its expansion is known to be well modeled by a mean-field interaction [3]. To verify that the observed tails (regions II and III in Fig. 2(b)) are not an artifact of the condensate, we compare our data to the complete 3D mean-field solution. The mean-field Gross-Pitaevskii (GP) solution beyond the Thomas-Fermi approximation may lead to the appearance of additional momentum components after TOF. These would result from slight modifications of the Thomas-Fermi real-space density occurring on the length scale (a 4 ho /R) 1/3 , which fixes the characteristic thickness of the boundary [3] (here a ho and R are, respectively, the oscillator length and the Thomas-Fermi radius and, for simplicity, we have assumed isotropic trapping).

The ground state is obtained by numerical simulation of the Gross-Pitaevskii equation in imaginary time, and the expansion dynamics are performed through real time propagation. The system size used in the calculations limits the simulated TOF. However, the distribution after 2.45 ms is observed to converge, indicating the complete conversion of mean-field interaction energy into kinetic energy. The numerical results for the distribution after a TOF of 2.45 ms are presented in Fig. S2, using the ballistic relation as defined in Eq. 1. The simulations clearly show that the mean-field Gross-Pitaevskii ground state is unable to reproduce the tails observed in the experiment.

Quantum depletion with local density approximation at T = 0

In uniform systems Bogoliubov's approach yields a population of momentum state equal to

n(k) = |v k | 2 = h2 k 2 /2m + mc 2 2 (k) - 1 2 (S1)
where (k) = h2 k 2 c 2 + (h 2 k 2 /2m) 2 1/2 is the Bogoliubov excitation spectrum. For a harmoncially trapped condensate in the Thomas-Fermi limit, the local speed of sound is c(r) = c 0 1 -(x/R x ) -(y/R y ) 2 -(z/R z ) 2 with c 0 = g s n 0 /m the speed of sound at the trap center and g s = 4πh 2 a s /m. In the LDA one finds

n( k) = 1 (2π) 3 |v k ( r)| 2 d r (S2) = R x R y R z (
R x R y R z 105π 2 1 k 4 ξ 4 (S4)
where ξ = h/ √ 2mc 0 is the healing length. Within the Bogoliubov approach in the LDA approximation, the Tan contact constant, defined as C/(2π) 3 = lim k→∞ n(k)k 4 , is equal to C LDA = (64π 2 /7)a 2 s N 0 n 0 , with a s the s-wave scattering length and N 0 (resp. n 0 ) the condensate atom number (resp. the condensate density). For metastable Helium, the calculation gives (C LDA /N 0 )/n 0 5.08 × 10 -15 m -2 .

Controlled heating sequence

To increase the temperature of our ultracold gas, we perform a controlled heating sequence using a 3D optical lattice. This sequence involves the adiabatic transfer of the gas from the optical dipole trap to the lattice in 30 ms, followed by a series of non-adiabatic lattice pulses of duration 0.5 ms during which the amplitude of the lattice is set to zero. The gas is then held in the optical lattice for 100 ms during which time it rethermalizes. Finally, we transfer the atom cloud adiabatically back to the initial optical dipole trap in 30 ms. Increasing the lattice depth while keeping the same sequence increases the final temperature of the gas at constant atom number and trapping frequencies.

FIG. 1 .

 1 FIG. 1. (a) Sketch of cloud expansion and detection by a micro-channel plate detector, yielding the 3D asymptotic momentum distribution (far-field regime). The initially cigarshaped Helium condensate (black) undergoes anisotropic expansion, inverting its aspect ratio. Quantum depletion and/or thermally excited atoms (grays) populate momentum states beyond those associated with the condensate, and are expected to have a spherical symmetry. (b) Measured 3D distribution of atoms n∞(k) after a 330 ms time-of-flight. The central dense part corresponds to the condensate while the isolated dots are excited particles outside of the condensate wavefunction. Also shown are the 2D projections, highlighting the condensate anisotropy.

µm 3 ]k [µm - 1 ]FIG. 2 .

 312 FIG.2. 1D momentum profiles obtained from cuts of the 3D data n∞(k). (a) Profiles along the radial (blue) and longitudinal (purple) directions. In linear scale, the tails are not visible (b) Log-log scale plot of the radial profile. The solid line is the scaling solution for the condensate in the Thomas-Fermi approximation (region I: k < 1.7 µm -1 ); the dotted line is a Bose distribution fitting the thermal wings (region II: 1.7 < k < 4.3 µm -1 ); the dashed line is a k -4 power-law fitting the high-momentum tails (region III: k > 4.3 µm -1 ). Solid lines are a smoothed average of the density data, and the lightly shaded band is the running standard deviation. The dark count rate corresponds to a level less than ∼ 10 -2 µm 3 , which is one order of magnitude below the lowest data point. The plotted range of momenta is limited by the physical size of the detector.

kkkk

  [µm -1 ] k B T a /µ = 1.73 T a = 286(5[µm -1 ] k B T a /µ = 1.24 [µm -1 ] k B T a /µ = 0.90 T a = 198(3[µm -1 ] k B T a /µ = 0.75

k 1 FIG. 4 .

 14 FIG.4. Contact constant C∞/N0 per condensed particle plotted as a function of the condensate density n0. The geometric trapping frequency ω/2π and the ratio kBTa/µ are indicated. The dashed line is the Bogoliubov prediction in the LDA, CLDA (see text), and the solid line is 6.5 × CLDA.

  Observation of Thermal and Quantum Depletion in a Bose Gas R. Chang, Q. Bouton, H. Cayla, C. Qu, A. Aspect, C. I. Westbrook and D. Clément Production of gaseous condensates and detection after time-of-flight

  FIG. S1. Observed momentum distributions along various directions, showing the anisotropic condensate distribution and the isotropic components associated with the thermal and the quantum depletion. The vertical dot-dashed line delimits the condensate (low momenta k) and the condensate depletion. (a) Illustration of the anisotropic momentum distribution of the condensate in the far-field regime (pancake-shaped distribution). The arrows represent the directions along which we plot the 1D momentum profiles with a small angular average (±10 • degrees) in the panels (b)-(g). (b)-(d) Blue dots are the measured 1D momentum profiles in the kz = 0 plane. (e)-(g) Green dots are the measured 1D momentum profiles in the ky = 0 plane.In each subplot, the solid line is the radial 1D momentum profile plotted in Fig.2of the main text (obtained from a large angular average, see above).

  FIG. S2. Simulation of the expansion dynamics of a 3DBose-Einstein condensate with the mean-field Gross-Pitaevskii equation. The parameters describing the condensate in the numerics (including atom number and trapping frequencies) are those of the experiment. The GP simulation is compared to the experimental data along the longitudinal (a) and radial (b) direction. It shows an excellent agreement in the region I associated with the condensate but does not reproduce the experimental tails. Numerical results and experimental distributions have been normalized to n∞(0, 0, 0) = 1.

We thank E. Cornell, Z. Hadzibabic, K. Mölmer, L. P. Pitaevskii, S. Stringari for thorough discussions about the interpretation of our results. We thank F. Nogrette for providing support on the detector and M. Mancini for reading the manuscript. We acknowledge financial support from the Région Ile-de-France (DIM Daisy), the RTRA Triangle de la Physique, the European Research Council (Senior Grant Quantatop), the LabEx PALM (ANR-10-LABX-0039), the International Balzan Prize Foundation (2013 Prize for Quantum Information Processing and Communication awarded to A. Aspect), the Direction Générale de l'Armement, the french National Research Agency (ANR 15-CE30-0017-04) and the Institut Francilien de Recherche sur les Atomes Froids.