Observation and control of quantized scattering halos

We investigate the production of s-wave scattering halos from collisions between the momentum components of a Bose-Einstein condensate released from an optical lattice. The lattice periodicity translates in a momentum comb responsible for the quantization of the halos' radii. We report on the engineering of those halos through the precise control of the atom dynamics in the lattice: we are able to specifically enhance collision processes with given center-of-mass and relative momenta. In particular, we observe quantized collision halos between opposite momenta components of increasing magnitude, up to 6 times the characteristic momentum scale of the lattice.

Introduction -Scattering experiments act as a probe revealing at a macroscopic scale the properties of collisional processes that occur at a microscopic scale. Since the seminal works of H. Geiger, E. Marsden and E. Rutherford [1][2][3], such experiments have remained a method of choice to probe atoms, molecules and their interactions. In short, the description of scattering in quantum mechanics gives rise to two remarkable features : firstly the quantum description of the collision process leads to its decomposition in terms of partial scattering waves [4], secondly, each of the collisional partners can itself be described as a matter-wave, which can be in a superposition of several components with well defined momenta, leading to multiple elementary collisional processes happening all "at once". As a result, quantum scattering exemplifies a key feature of quantum mechanics: waveparticle duality.
With the advent of cold atom samples, this topic has been revisited with only a few partial waves involved. The characterization of the collisional properties, and in particular of the s-wave scattering length, was performed either using photoassociation measurements [5][6][7][8] or studying the kinetics towards equilibrium of an atomic sample [9,10]. In this latter type of measurements, the interplay between partial waves turns out to be subtle: the thermalization rate involves partial waves interferences while the collision rate does not [11]. Using a 1D collider geometry, the experiments of Refs. [12,13] have captured quantum scattering in its purest form: at low energy, the s-wave collisions create a spherical shell of pair-correlated atoms, and at slightly higher energy the volume occupied by the scattered atoms reflect the interference between partial waves. Such colliderlike experiments have also been carried out with different species [14,15]. More recently, second-and third-order correlations between momentum-correlated atoms in a collisional halo have been investigated [16,17], opening quantum-nonlocality tests to ensembles of massive particles [18].
In this article, we engineer the collisions between atomic ensembles in a multiple 1D collider, using an outof-equilibrium Bose-Einstein condensate (BEC) of 87 Rb released from an optical lattice. The collisions occur Figure 1: Scheme of the experimental protocol. A Bose-Einstein condensate (red clouds) is initially loaded in an optical lattice (dashed blue line). a The lattice is suddenly phaseshifted by an amount ϕ = ϕ0 (solid blue line); here ϕ0 = 180 • . The subsequent dynamics in the lattice modifies the momentum distribution. b After a sufficiently long time-of-flight, the momentum distribution exhibits diffraction peaks, separated by kL = h/d, characteristic of the wave nature of the BEC in the lattice (red disks) along with isotropic scattering halos due to s-wave scattering between atoms (light red disk). through s-wave scattering, between the discrete momentum components of the BEC, in the course of the timeof-flight, leading to the appearance of scattering halos [19]. Through accurate control of the lattice phase and amplitude, we show that it is possible to engineer the dynamics of the BEC in the lattice before the release, and tailor the wavefunction of the atoms after release and expansion, in order to selectively enhance the collision processes between specific momentum components of the BEC. We can thus produce quantized scattering halos from collisions with a chosen relative and/or centerof-mass (c-o-m) momentum.  Background -When a BEC is loaded in an optical lattice, the periodic structure of the lattice imprints on the wavefunction. In particular the resulting momentum distribution is made up of equally spaced peaks, separated by an interval h/d, where d is the lattice spacing and h the Planck constant. The shape of each individual peak is set by atomic interactions and the extra confinement superimposed to the lattice. This momentum distribution can be measured by releasing the atoms from the lattice and allowing the atomic cloud to ballistically expand for a sufficiently long time-offlight (TOF), t T OF : the spatial density n(r) then reproduces the initial momentum densityñ(p) up to a scaling: n(r, t T OF ) =ñ(p = mr/t T OF , t = 0), with m the mass of the atom (see Fig. 1).
However, this comb pattern resulting from the wave nature of the BEC is only a partial description of the final atomic density distribution: for interacting atoms, collisions originating from the particle nature of matter may occur during the ballistic expansion. In the lowenergy regime characteristic of ultracold atoms, s-wave elastic scattering dominates which, due to energy and momentum conservation in the collision process, results in spherical halos centered on the center of mass of the two colliding atomic wavepackets (see [20,21] and Appendix A).
Thus, spherical collision halos are expected to appear in the momentum distribution between the diffrac-tion orders after time-of-flight, as depicted in Figure 1. Each diffraction order, centered on a position x j = j × h/(md) × t T OF (j ∈ Z), contains a fraction π j of the number of atoms N initially in the BEC. These fractions determine the atomic density in the halos : in a perturbative approach, the number of collisions between orders j and k, and therefore the number of atoms scattered in the corresponding halo, is proportional to the product π j π k ( [22], Appendix A).
In this work, we engineer the (dominant) scattering halos, by precisely controlling the state of the BEC in the lattice before release. To this end, we apply a sudden displacement to the lattice on a scale smaller than the lattice spacing. This triggers out-of-equilibrium dynamics of the BEC inside the lattice. Controlling the duration of this evolution in the lattice prior to release (see Fig. 1), we effectively tailor the momentum distribution {π j }.
Methods -We perform our experiments in a hybrid trap [23] in which we obtain pure rubidium-87 Bose-Einstein condensates of 2 · 10 5 atoms in the lowest hyperfine state |F = 1, m F = −1 . These BECs are loaded in a onedimensional optical lattice produced by two counterpropagating laser beams of wavelength λ = 1064 nm superimposed to the optical dipole beam of the hybrid trap. In the optical lattice, the atoms experience the following potential : where E L = 2 k 2 L /(2m) (E L = 4E R , with E R the recoil energy) and k L = 2π/d are respectively the energy and the wavevector associated to the lattice. The dimensionless depth of the lattice s is independently and precisely calibrated [24] for each experiment presented here. The phase ϕ is set by the relative phase between the two phase-locked acousto-optic modulators controlling the lattice beams. In 87 Rb, higher-order d−wave collisions occur for energies of ∼ 200 µK [12,13], while in the experiments shown here we impart at most ∼ 25 µK of collisional energy. We are therefore in the purely swave scattering regime. Within this regime, in order for the changes in momentum distribution to directly relate to the number of atoms in the halos, the scattering cross section needs to be independent from the relative speed across all our experiments (see Appendix A), and not be affected by scattering suppression due to the superfluid onset [25]. We have checked that the minimum relative velocity we can impart in a collision is of h/(md) = 8.6 mm · s −1 , which is much larger than the peak sound velocity for the BEC in one of the colliding diffraction orders (c ∼ 1.5 mm · s −1 ). Therefore the scattering cross-section is the same for all collision halos presented here.
In order to tailor the momentum distribution after expansion, we use the following general procedure: we first load adiabatically the BEC for ϕ = 0; then a sudden phase shift ϕ 0 is applied [31], triggering out-ofequilibrium dynamics (t = 0 + ); in the following holding time, the momentum distribution {π j (t)} evolves in the lattice; finally, the lattice is released at t = t hold . The atom cloud then expands for a duration t T OF , and an absorption image of the resulting density n(r) is taken to measure the momentum distribution and the scattering halos.
As a guide to our engineering, we numerically compute the evolution of a one-body wavefunction, initially in the ground state of an infinite optical lattice of depth s. After the phase shift ϕ 0 , the wavefunction is decomposed on the lattice bands, and evolves during the holding time t hold . We then extract the final momentum distribution (see for example Results -In a first series of experiments, we demonstrate control of the c-o-m momentum of the scattering halo. For that purpose, we use the phase shift ϕ 0 = 140 • , which puts the atom clouds in the lattice on the side slope of each well (see Fig. 1). As the atoms reach the bottom of the well in the following dynamics, we expect a high co-m momentum, which is higher the deeper the lattice is, with a small dispersion of the distribution over diffraction orders.
As the dominant collision halo will occur between the two mostly populated diffraction orders, we define a visi-bility parameter V as the difference between the fractions π j of the second and third most populated diffraction orders in the momentum distribution. It is a measure of how well the two diffraction orders contributing to the main halo "stand out", and varies between 0 (the third highest order is as populated as the second) and 0.5 (the two orders contributing to the main halo are the only ones populated). In Fig. 2a, the values of the visibility V obtained from the wavefunction simulation are plotted over a large range of lattice depths s and holding times t hold . As we are here interested in controlling mainly the c-o-m momentum, the visibility V is only plotted if the two dominant orders are next to each other (separated by 1 × h/d in momentum), and is otherwise set to 0. Finally we also indicate the value of the c-o-m momentum by a color code. As intuitively expected, higher c-o-m momenta can be reached for deeper lattices, for a holding time that gets shorter the deeper the lattice gets.
At a given depth, our analysis in terms of visibility shows that "patches" of lower c-o-m momentum values can be obtained for longer holding times. However the visibility in these patches is less pronounced, meaning there are likely more than two significantly populated orders of diffraction for these parameters, and we can expect that several collisional processes will occur simultaneously in the experiment. Therefore we have not represented experimental data for these values.
In Fig. 2b-g, we represent a few snapshots of the density distributions obtained after time-of-flight for parameters that are indicated in Fig. 2a. For each of these measurements, the lattice depth was calibrated independently [24]. The disk-shaped halo due to collisions between the two most populated orders, separated by 1 × h/d, is clearly visible, and is the main feature besides the regular structure of the diffraction orders. Alongside the absorption pictures we represent the measured histograms {π j } of the populations in the diffraction order p j = j × h/d. These are compared to the calculated histograms from the wavefunction simulation, and the visibility parameter from the simulation is indicated by a red-shaded area. Note that the sign of the c-o-m momentum can be easily changed by changing the sign of the shift ϕ 0 .
Even though the fraction of atoms in the halo can get rather large (see below), the agreement between the simulated distributions, which do not include interactions, and the measured distributions, for which collisions have occurred, is very good. This can seem surprising, as atoms are removed from the diffraction peaks during the collisions. However in any collision between two orders of diffraction, an equal number of atoms is removed from the involved peaks, after which the measured distribution (as shown in blue in Fig. 2) is the new normalized distribution in the peaks. We have simulated that even with a fraction of colliding atoms as high as 30%, the expected change between the non-interacting theoretical distribution and the distribution obtained after accounting for the removal of the colliding atoms and normaliz-ing remains small (<5%). This means that the one-body wavefunction simulation is a surprisingly robust guide in our investigation of collision halos, as is further evidenced below (see Fig. 4).
In a second series of experiments, our aim was to control specifically the relative momentum of the collisions. To that end, we apply a shift ϕ 0 = 180 • : the atom clouds in the lattice are placed at the top of the lattice potential, and split in two identical clouds that fall on either side (see Fig. 1). In a classical picture, after a holding time that allows the split clouds to reach the bottom of the wells, momentum orders of opposite sign should be equally populated, with zero c-o-m momentum. Their typical magnitude should also increase with lattice depth s.
In Fig. 3a, we plot the values of the visibility V , as defined before, over a range of values of lattice depths s and holding times t hold , for ϕ 0 = 180 • , as calculated from the wavefunction simulation. We do not impose a condition anymore on the separation between the most populated orders (note that in this case, the theoretical distributions are symmetrical and the two highest peaks are equally populated). The diameter of the expected collision halo (with zero c-o-m momentum) is indicated by a color code. As expected, the relative collisional momentum increases with lattice depth. In the main feature of this figure, the relative momentum of the most populated orders increases in steps of 2 × h/d as the lattice depth is increased, for decreasing holding times (a behaviour that is expected intuitively).
In Fig. 3b-g we represent snapshots of the density distributions obtained after time-of-flight for parameters that are indicated in Fig. 3a. Disk-shaped collision halos of increasing diameters, quantized in units of 2h/d, are clearly visible. To the left of the absorption images, we represent the momentum distributions π j as extracted from absorption images. We compare the experimentally measured histograms with those from the wavefunction simulation for the same parameters s and t hold , and find once again very good agreement. Since large collision halos become very dilute at the long time-of-flight needed for the reconstruction of the histograms (a few tens of ms), we used a different time-of-flight for the visualization of the halos (typically a few ms).
The observed fraction of colliding atoms we measure is in general larger than predicted by a perturbative theory (see Appendix A). To measure it, we count the atoms in the collision halos when the diffraction orders are masked, as in Fig. 4a. We account for the portion of disk hidden, and calculate the ratio of the number of atoms in the halo to the sum of the numbers in the halo and in the diffraction orders. We measure typically a fraction of 30 ± 3% in the data presented Fig. 2b and Fig. 3b (from a statistical average of the fractions in 15 absorption images), where the theory predicts 20%. This is in contrast to the results in [22] which were well described with a similar prediction; but while these results involved typically a few percent of the atoms of the BEC in the collisions between two orders, we observe collisions between highly populated orders, and get typically five times larger atom fractions involved in the collisions, at which point a non-perturbative approach may be needed to get quantitative agreement. We also independently checked that this higher number is not due to collisions happening inside the lattice during the non-equilibrium evolution time t hold (see Appendix C).
We have also verified that the diameter of the dominant collision halo in momentum space is indeed quantized in terms of h/d (the momentum scale set by the lattice), and that larger diameters are reached for larger lattice depths. This analysis was performed independently from the known position of the diffraction orders, or any a priori visibility calculation. This is also performed for a phase shift of ϕ 0 = 180 • , and for the experimental parameter values shown in Fig. 3a. For each of the parameter values {s,t hold }, we record a sequence of images with increasing TOF. On each of these images, as in Fig. 4a, we mask the diffraction orders, to focus on the most visible disk-shaped halo. By performing an angular average, we obtain a radial distribution, clearly showing the edge of the collision halo. We fit the position of this edge with a sigmoid function (Fig. 4b), which gives us the radius of the scattering halo for a given TOF. We then extract the speed of expansion of the scattering halo by fitting the linear growth of the radius with TOF (Fig. 4c).
This procedure was repeated for multiple values of s, choosing for each depth a holding time t hold yielding a scattering halo with the strongest signal. The results are shown in Fig. 4d. We find indeed that the expansion velocity of the dominant halo, as measured from the halo only, shows sharp jumps between discrete values that are an integer multiple of v L = h/(md), the velocity scale set by the lattice. As the lattice depth increases, so does the collisional energy, which is converted from the lattice potential energy during t hold , and so does the expansion velocity of the dominating halo. Alongside the experimental data in Fig. 4d, we represent the expected collision velocity for the dominant scattering halo, as given by the point with highest visibility for a given depth s on the map Fig. 3a, obtained from the wavefunction simulation. We find that the experimentally observed steps for the relative collision velocity in the dominant halos are in remarkably good agreement with the characterization in terms of visibility of the diffraction orders.
Discussion and conclusion -In this work, we have demonstrated that we can engineer at will the collision halos from ultracold atoms released from an optical lattice. This is achieved through the control of the dynamics of a BEC in the optical lattice, which allows us to tailor the momentum distribution giving rise to the collisions after release from the trap. We have shown that we can selectively populate diffraction orders to impart a large momentum either to the center-of-mass motion of colliding atom clouds, or to their relative motion. The search for appropriate experimental parameters has been guided by the simulation of the dynamics of a non-interacting BEC in the lattice, and the identification of two highly populated diffraction orders with the desired characteristics. The energies involved in the collisions are quantized in terms of the lattice characteristic momentum but can be adjusted through a wide range of values.
The protocol used here to demonstrate separate control on center-of-mass and relative momentum can be extended to realize arbitrary combinations of the two. In Fig. 5, we provide an example where we control simultaneously the collision halo diameter (2 × h/d) and the center-of-mass momentum, (1 × h/d). This still relies on the simple control of the static lattice depth and a sudden phase shift, and a visbility search for the desired diffraction orders. This work could be extended further using more elaborate dedicated schemes for the control of the momentum distribution, by an appropriate timedependent shaping of both the depth and phase of the lattice before time-of-flight expansion. This could provide a fairly simple and general technique alongside more conventional atom-optics schemes [26]. Such an optimization could also be of interest for matter-wave interferometry where the engineering of momentum superpositions of atomic ensembles is a key technique [27,28]. In this section, we provide a quantitative theoretical description of the atom-atom scattering processes giving rise to the halo structures that are observed in the experiments. The perturbative framework that we develop is inspired from previous works that were investigating colliding Bose-Einstein condensates [20,29]. The results that it yields are consistent with the calculations undertaken in Ref. [22].
Our starting point is a perfect Bose-Einstein condensate that is prepared within a one-dimensional optical lattice. The atoms of the condensate are sharing the single-particle wavefunction with r ≡ (x, y, z) and the lattice being oriented along the z axis of the coordinate system. ψ l are the Fourier components of the periodic condensate wavefunction in the lattice. The function is a Thomas-Fermi envelope that accounts for the presence of a weak overall elliptic confinement with the longitudinal and transverse frequencies ω || and ω ⊥ , respectively, and withω = (ω || ω 2 ⊥ ) 1/3 . The associated Thomas-Fermi radius R is straightforwardly evaluated as where N is the number of atoms in the condensate and a s denotes the s-wave scattering length. At time t = 0 the lattice and confinement potentials are switched off and the atomic cloud is thereby allowed to freely expand. If atom-atom interactions could be completely neglected from that instant on, the momentum distribution of the atoms would be simply given by the Fourier transform of the condensate wavefunction, i.e., we would have n 0 (p, t) = N |Φ 0 (p)| 2 with is the Fourier transform of the Thomas Fermi profile (A2). As we generally have k L R 1 in the experiment, Eq. (A5) suggests the appearance of tighly localized momentum density peaks centered about integer multiples of the lattice momentum l k L , each one of those populated with N |ψ l | 2 atoms.
The key approximation that we undertake here is to assume that we can account for the presence of atomatom interaction during the free expansion process of the condensate in a perturbative manner, using first order quantum perturbation theory. Employing the interaction representation, the time evolution of the many-body state |Ψ t describing the atomic cloud is approximately written as in linear order in the s-wave scattering length a s . Here, with g = 4π 2 a s /m is the two-body interaction Hamiltonian expressed in terms of the atomic creation and annihilation operators in momentum space which evolve according toψ † t (p) =ψ † (p) exp[itp 2 /(2m )] as well aŝ ψ t (p) =ψ(p) exp[−itp 2 /(2m )], respectively, and which fulfill the bosonic commutation relation We can therefore rewrite Eq. (A8) aŝ and use the standard definition of the momentum-space field operators giving rise tô whereb 0 is the annihilation operator associated with the condensate orbital.
In the following, we focus on regions in momentum space located far away from the lattice density peaks |φ(p−l k L e z )| 2 which we obtained in zeroth order in the interaction strength. For a momentum p being in such a scarcely populated region we can safely setΦ 0 (p) = 0 and henceψ(p)|Ψ 0 = 0. The atom density detected at such a momentum is therefore determined as in lowest nonvanishing order in the s-wave scattering length. Evaluating and using Eq. (A11) in combination with Eq. (A5) , we obtain the expression The time integral appearing in Eq. (A15) is straightforwardly evaluated yielding t 0 dt e iω l,l (q,q )t = t sinc[ω l,l (q, q )t/2]e iω l,l (q,q )t/2 (A17) with sinc(ξ) ≡ sin(ξ)/ξ. It has a similar effect as Dirac's delta distribution for large t → ∞ insofar as it would yield vanishing contributions within Eq. (A15) for values of ω l,l (q, q ) that are different from zero. Hence, given the fact that the prefactorφ(p 1 )φ(p 1 ) restrains the effective integration domain of p 1 and p 1 to a very narrow region about the origin, we can infer that for a sufficiently long evolution time t a nonvanishing momentum density can be encountered only in the immediate vicinity of the surface of spheres in momentum space that have their north and south poles at the points l k L e z and l k L e z and are therefore characterized by the equation for any pair of integers l, l ∈ Z.
To quantitatively calculate the momentum density according to Eq. (A12), we can set by the choice that we made for p, since we effectively have |p 1 + p 1 | k L according to the above reasoning. Similarly, we can evaluate when this commutator is multiplied with the prefactor φ * (p 2 )φ * (p 2 )φ(p 1 )φ(p 1 ). The two spheres on which p has to be simultaneously located in order to yield a significant momentum density would therefore have the same center at 1 2 (l 1 + l 1 ) ke z = 1 2 (l 2 + l 2 ) ke z , which implies that they are actually identical and we have either l 1 = l 2 and l 1 = l 2 or l 1 = l 2 and l 1 = l 2 . Using (A21) we can, consequently, write the momentum density (A12) as represents some sort of phase-space wavefunction that describes the collision process between two momentum components l and l of the condensate wavefunction, with p 0 = 1 2 (l − l ) k L e z and with the origin in momentum space being set at 1 2 (l + l ) k L e z . In the context of scattering halos, we are specifically interested in the case l = l , which implies a finite relative momentum 2p 0 ∼ k L of the two wavefunction components. We can therefore neglect the p 1 · p 1 term arising in the exponent within Eq. (A24) as it will become negligibly small with respect to the other terms in this exponent, owing to the presence of the prefactorφ(p 1 )φ(p 1 ). This consequently yields Making use of the fact that the Thomas-Fermi profile (A3) characterizing the condensate wavefunction is of finite extent, we can safely take the limit t → ∞ in the above expression and obtain for the momentum density of scattered atoms in the longtime limit, where we use p 0 = ±p 0 e z with p 0 = |p 0 |. This expression can be further simplified by using the fact that for p = p 0 the term exp[iR(p 2 − p 2 0 )(τ − τ )/( p 0 )] has, in the limit Rp 0 / → ∞, a similar effect as a delta function in τ − τ within Eq. (A26). We are therefore entitled to rewrite this expression as which for Rp 0 / → ∞ yields a tight and isotropic concentration of the momentum density around the sphere with the radius p = p 0 . As is clearly seen from Eq. (A27), the width of these scattering halos in momentum space is of the order of ∆p ∼ /R, which is in agreement with the experimental findings. This is confirmed in Figure 6 which shows a numerical evaluation of Eq. (A27) for the experimental parameters under consideration and for various choices of the condensate population and the overall confinement frequency (assuming a spherical trap). The atomic density in momentum space is concentrated about a sphere of radius p 0 = |l − l | k L /2 = k L /2 in this example (|l − l | = 1). The width (FWHM) of this scattering halo is found to roughly equal π /R in all of the four studies cases. Note that these density profiles cannot be directly compared with the profile shown in Fig. 4b since the latter is the result of the integration of the scattering sphere along the imaging axis, followed by angular averaging.
The total number of atoms that are scattered as a consequence of the collision between the l and l components of the condensate wavefunction (with l = l ) is then straightforwardly calculated as N coll.at. = 2|ψ l | 2 |ψ l | 2 N with (A29) In this last expression we have isolated the scattering cross section σ, which is here independent of the colliding orders l and l . Eq. A28 is quantitatively rather similar to the analogous Eq. (2) for the number of collisions within Ref. [22], which was obtained there in the framework of a three-dimensional optical lattice. The scaling with (ω ⊥ /ω || ) 2/3 translates the fact that two cigarshaped condensates moving across each other along their symmetry axis give rise to more atom-atom collisions than two pancake-shaped condensates. For a mean angular frequencyω = 2π × 50 Hz, a number of atoms equal to 2 × 10 5 , and an ideal, equal-weights collision, |ψ l | 2 = |ψ l | 2 = 0.5, we find a fraction N coll.at. /N 20 % of atoms that collide. This is somewhat smaller than what is observed experimentally, which may indicate that a non-perturbative approach is required for a quantitative agreement. This equation can easily be solved in matrix form to obtain the eigenenergies of the system and the decomposition of the eigenstates on the plane wave basis. The BEC is initially loaded in the lowest energy eigenstate of the lattice, corresponding to a quasi-momentum q = 0 : |Ψ(t = 0 ) = |Φ 0,0 . After the sudden phase shift, the wavefunction in the reference frame of the lattice takes the form : Since the lattice preserves its periodicity in the shift, the quasi-momentum is conserved. The wavefunction after the shift can therefore be decomposed onto the eigenstates of the shifted lattice with quasi-momentum q = 0 to obtain the evolution during the holding time t hold : |Ψ(t hold ) = n∈N α n e −i E n,0 t hold |Φ n,0 , with α n = Φ n,0 |Ψ(t = 0 + ) = l∈Z C * (n,0) l C (0,0) l e −ilϕ0 . Finally the population in the diffraction order of momentum j × h/d is obtained by calculating the projection of the state |Ψ(t hold ) on the corresponding eigenstate of the momentum operator |χ jk L : π j = | χ jk L |Ψ(t hold ) | 2 = n∈N α n C (n,0) j e −i E n,0 t hold 2 .
Appendix C: Impact of in-lattice dynamics on the collisional halos In this appendix we show that the collisional halos observed in this article are determined only by the momentum distribution at the end of the out-of-equilibrium dynamics in the lattice, of duration t hold . This is illustrated in figure 7, where we have represented the collision halos recorded for a phase shift ϕ 0 = 180 o and a lattice depth of s = 1.92 ± 0.07, for various durations t hold . These are expressed as a fraction of the period T of the dipolar motion in the lattice (here T 140 µs).
In the first couple images, no collision halo between orders ±1 × k L can be seen, and it only appears for t hold = T /4, when these orders are significantly populated. As the dynamics is extended further, and the final population in orders ±1 × k L reduces again, the corresponding halo is not observed anymore. If collisions were occurring inside the lattice, the collision halo would remain visible for all holding times after its first appearance. This demonstrates that the halos are determined by the momentum distribution upon release from the lattice, and not affected by the prior dynamics. t hold = 0 t hold T/8 t hold T/4 t hold 3T/8 t hold T/2 Figure 7: Single-shot images of diffraction orders for varying holding times t hold expressed as a fraction of the dipolar motion period T (here T 140 µs), for a shift of ϕ0 = 180 • , and a lattice depth of s = 1.92 ± 0.07. The color scale is set to enhance the collision halos (clipping high values on some diffraction orders). The collisional halo between orders ±1 × kL is only observed when these orders are populated at the time of release from the lattice, for t hold T /4.