PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of ncommsLink to Publisher's site
 
Nat Commun. 2017; 8: 14854.
Published online 2017 March 23. doi:  10.1038/ncomms14854
PMCID: PMC5376650

Inelastic collisions of ultracold triplet Rb2 molecules in the rovibrational ground state

Abstract

Exploring and controlling inelastic and reactive collisions on the quantum level is a main goal of the developing field of ultracold chemistry. For this, the preparation of precisely defined initial atomic and molecular states in tailored environments is necessary. Here we present experimental studies of inelastic collisions of metastable ultracold Rb2 molecules in an array of quasi-1D potential tubes. In particular, we investigate collisions of molecules in the absolute lowest triplet energy level where any inelastic process requires a change of the electronic state. Remarkably, we find similar decay rates as for collisions between rotationally or vibrationally excited triplet molecules where other decay paths are also available. The decay rates are close to the ones for universal reactions but vary considerably when confinement and collision energy are changed. This might be exploited to control the collisional properties of molecules.

Recent advances in the preparation of ultracold molecular samples in well-defined quantum states1,2,3,4,5,6,7 sparked increasing interest in studying molecular collisions and chemical reactions on a pure and fundamental level8,9,10,11,12. Such experiments were first carried out with highly excited molecules13,14,15,16,17,18,19,20,21,22,23,24 (also in the context of Efimov physics25,26,27,28,29,30), which can vibrationally relax in a collision. For molecules in the vibrational ground state this decay is barred, but other interesting reaction paths remain. First investigations of reactive or inelastic loss of singlet molecules in the rovibronic ground state with polar KRb and RbCs have recently been carried out5,31,32,33,34. Besides the electronic ground state An external file that holds a picture, illustration, etc.
Object name is ncomms14854-m1.jpg, the metastable triplet state An external file that holds a picture, illustration, etc.
Object name is ncomms14854-m2.jpg is of special interest for collision experiments. If collisionally long-lived, such triplet molecules would allow for many interesting applications, such as tunable Feshbach resonances, due to their sizeable magnetic moment. To investigate chemical reactions between cold molecules, optical lattices are a convenient testbed as they allow for either isolating the molecules from each other or letting them collide. In addition optical lattices offer the possibility to control the dimensionality of the scattering process and to tune the interaction35,36,37,38. On the basis of this approach, it has been shown that strong inelastic collisions induce correlations and can inhibit particle loss in a molecular sample, a manifestation of the quantum Zeno effect15,34.

In this work, we present measurements on ultracold collisions of metastable triplet molecules that are internally in their lowest energy state. Specifically, we use 87Rb2 dimers of the An external file that holds a picture, illustration, etc.
Object name is ncomms14854-m3.jpg state which are in the lowest hyperfine level of the rovibrational ground state. As reference measurements, we carry out collision experiments with rotationally excited molecules (rotational quantum number R=2) and with vibrationally highly excited Feshbach molecules. Initially, the dimers are prepared in a cubic 3D optical lattice with at most a single molecule per lattice site (see Methods). By quickly ramping down one of its directions, the lattice is converted into an array of quasi-1D potential tubes (Fig. 1a). Subsequently, molecules within the same tube collide with tunable relative energies on the order of μK × kB, far above the Tonks gas regime of the work of Syassen et al.15. A single tube is typically filled with only a few molecules and can be considered as a closed few-body system since tunnelling between the tubes is negligible.

Figure 1
Initial spatial distribution of molecules.

Whenever a collision between two molecules is inelastic or reactive, enough energy is released to expel all products out of the lattice. This is because the lattice depths are comparatively shallow, being on the order of about 10 μK × kB, which corresponds to ~210 kHz × h. After a given interaction time t we measure the total number of remaining Rb2 molecules N(t) and the width of the whole cloud An external file that holds a picture, illustration, etc.
Object name is ncomms14854-m4.jpg(t) along the tubes. From the observed decay of N we conclude that a large part of the molecules is already lost in the first possible collision, quite independently of the internal vibrational excitation. Using a simple model we extract the decay rate constants and investigate how they depend on the confinement of the potential tubes. These results are then compared to predictions of a quantum defect model39,40.

Results

Collisions of Feshbach molecules

In the following, we first discuss the collision experiments with Rb2 Feshbach molecules as this will help us analysing the data for the v=0 states. Figure 2a shows three data sets of N(t) corresponding to different confinements of the potential tubes. We observe a strong loss of molecules within the first tens of milliseconds. It is striking that the decay takes place in a step-wise fashion which, after ~100 ms, gives way to a much slower exponential decay with a corresponding time constant of >1 s. This slow exponential decay is similar to the one we observe in a deep 3D optical lattice (Fig. 2a, inset), which is due to background gas collisions and spontaneous photon scattering41. Figure 2b shows that the width An external file that holds a picture, illustration, etc.
Object name is ncomms14854-m5.jpg oscillates synchronously to the steps.

Figure 2
Dynamics of Feshbach molecules for various confinements.

We interpret these dynamics as follows: as one direction of the 3D lattice is quickly ramped down (within 400 μs), the particles are suddenly released from their individual lattice sites into one-dimensional (1D) tubes. Along these tubes there is a harmonic confinement with trap frequency ωz due to the Gaussian intensity profile of the 2D lattice laser beams. The molecules will synchronously undergo an oscillatory motion along this direction with a period T=2π/ωz, while being strongly confined transversally. As a result, the width of the observed cloud An external file that holds a picture, illustration, etc.
Object name is ncomms14854-m6.jpg oscillates with 2ωz (Fig. 2b). Whenever the cloud is small and dense, the probability for molecular encounters and losses is increased. On the other hand, if the cloud is large the dimers are separated from each other and N(t) stays almost constant. Thus, the longitudinal oscillatory motion explains the step-like decay of the molecules. Once all inelastic collisions have taken place the fast losses stop and the remaining signal corresponds to single molecules in the tubes. This regime is typically reached after 50–100 ms.

Modelling of collision dynamics and analysis

To model the molecular decay in a quantitative way, we first reconstruct the distribution of the molecules in our 3D optical lattice. The molecules are initially produced from a Gaussian-shaped cloud of ultracold atoms in the optical lattice via magnetic Feshbach ramping (see Methods). We assume Poisson statistics for the atomic occupation of each individual lattice site. Only lattice sites that are occupied by exactly two atoms will finally be occupied by a single molecule. All other lattice sites will end up empty. Thus, a given total atom number and cloud size fixes the initial molecular distribution, see Fig. 1b. We have verified that the predicted cloud size and total number of the Feshbach molecules indeed agree with our measurements (see Methods). Taking the molecule distribution of Fig. 1b the histogram of the resulting filling of the 1D tubes is depicted in Fig. 1c. The average occupation is ~2.2. The histogram helps us to gain additional insights. If we assume inelastic two-body collisions to be the only source of particle loss, then the number of molecules remaining after a decay time of ≈100 ms equals the number of tubes with an odd initial filling. Evenly occupied ones, by contrast, end up empty. By comparing the experimentally measured fraction of remaining particles to this prediction we can check the model for consistency (see Supplementary Note 1 and associated Supplementary Fig. 1).

We now discuss the dynamics within a single 1D tube in more detail. Starting with a classical treatment, we consider each molecule as a point-like particle, which is initially localised and at rest in a single lattice site of the 3D lattice. The molecules are all released exactly at the same time and consequently meet in the centres of the tubes precisely at t=T/4=π/(2ωz). At that particular moment the molecular cloud size ideally vanishes along the z-direction, inelastic collisions take place and the total molecule number decreases abruptly. In our measurements, however, we do not observe such abrupt steps, but rather smoothed ones (Fig. 2a). We explain this fact by a non-vanishing initial velocity distribution of the molecules as a result of the Heisenberg uncertainty relation. Therefore, we leave the classical picture of point-like particles and rather describe each molecule i as a 1D quantum mechanical wave packet An external file that holds a picture, illustration, etc.
Object name is ncomms14854-m7.jpg, which is centred at χi(t) and has the width An external file that holds a picture, illustration, etc.
Object name is ncomms14854-m8.jpg (cf. Fig. 3a). This leads to a finite width in the molecule's velocity distribution, initially given by Δv(0)≈ħ/(mAn external file that holds a picture, illustration, etc.
Object name is ncomms14854-m9.jpg(0)), where m is the molecular mass. As t progresses, χi(t), An external file that holds a picture, illustration, etc.
Object name is ncomms14854-m10.jpg and Δv(t) oscillate with ωz, 2ωz and 2ωz, respectively (Supplementary Note 2). The dynamics for three particles in a tube is depicted in Fig. 3b (upper part). As expected, the collision times are now somewhat smeared out, however, every particle will still pass by every other one in the 1D tube within T/2 (assuming that no collisions occur in the meantime). Whenever two molecules collide, there is a certain likelihood for an inelastic process. A collision between two molecules i and j is possible as long as the spatial overlap An external file that holds a picture, illustration, etc.
Object name is ncomms14854-m11.jpg of their wave packets (cf. Fig. 3b, lower part) does not vanish. To describe the loss rate in a single 1D tube we choose the following ansatz (see Supplementary Note 3 for a derivation)

Figure 3
Dynamics of molecular wave packets in a quasi-1D potential tube.
An external file that holds a picture, illustration, etc.
Object name is ncomms14854-m12.jpg

where we sum twice over all molecular pairs (i, j). An external file that holds a picture, illustration, etc.
Object name is ncomms14854-m13.jpg is a normalization factor. An external file that holds a picture, illustration, etc.
Object name is ncomms14854-m14.jpg is the mean decay coefficient. We note that in the limit of high tube occupation numbers An external file that holds a picture, illustration, etc.
Object name is ncomms14854-m15.jpg describes the dynamics of a 1D thermal gas with density n via the rate equation An external file that holds a picture, illustration, etc.
Object name is ncomms14854-m16.jpg.

We use equation (1) to analyse the measured decay curves and to extract An external file that holds a picture, illustration, etc.
Object name is ncomms14854-m17.jpg. For this, we fit the solution of equation (1) (see Methods) to the measurements. Indeed, the only free fitting parameter in our model is An external file that holds a picture, illustration, etc.
Object name is ncomms14854-m18.jpg. This is because we can experimentally determine the initial spatial distribution of the molecules, and the initial width An external file that holds a picture, illustration, etc.
Object name is ncomms14854-m19.jpg(t=0) of the molecular wave packets is derived by calculating the release dynamics from the 3D optical lattice into the 1D tubes (see Supplementary Note 2 and associated Supplementary Fig. 2). For a given molecular sample we assume An external file that holds a picture, illustration, etc.
Object name is ncomms14854-m20.jpg(t=0) to be identical for all particles. The values of An external file that holds a picture, illustration, etc.
Object name is ncomms14854-m21.jpg(t=0) for the present work range between 0.17 to 0.26alat, where alat=532 nm is the lattice constant. Generally, the width An external file that holds a picture, illustration, etc.
Object name is ncomms14854-m22.jpg(t=0) determines how smoothed out the decay steps are, while An external file that holds a picture, illustration, etc.
Object name is ncomms14854-m23.jpg determines their relative heights. Therefore, a variation in An external file that holds a picture, illustration, etc.
Object name is ncomms14854-m24.jpg(t=0) does not have a strong influence on the extracted rate coefficient An external file that holds a picture, illustration, etc.
Object name is ncomms14854-m25.jpg. As can be seen in Fig. 2a the fitted curves agree well with the measurements. Our results for the decay rates are An external file that holds a picture, illustration, etc.
Object name is ncomms14854-m26.jpg mm s−1 for the corresponding radial trap frequencies of ωr=2π × (8.2, 11.6, 17.0) kHz, respectively.

We now get back to Fig. 2b. A strong damping of the oscillations of the cloud size An external file that holds a picture, illustration, etc.
Object name is ncomms14854-m27.jpg(t) along the tubes is observed, which at first might be unexpected for 1D systems where thermalization is generally suppressed42. We attribute the damping mainly to the fact that the reaction rate increases with collision energy (as will be shown below) and therefore particles with higher kinetic energy are lost faster. A further discussion also of other possible contributions to the damping is provided in Supplementary Note 4. In addition, we would like to note that the observed oscillating cloud size in Fig. 2b is generally larger than expected from our model calculations. This is due to limitations in the effective imaging resolution (see Methods).

Collisions of vibrational ground state molecules

Next, we study the inelastic collisions of the Rb2 triplet v=0 molecules that are produced in precisely defined internal quantum states via coherent optical transfer starting from the Feshbach state (see Methods). Figure 4b,c shows decay curves of these dimers with rotational quantum number R=0 and R=2, respectively. As for the Feshbach molecules, the observed loss is almost entirely due to collisions since the measured lifetime in the absence of molecular encounters is on the order of several seconds43. For direct comparison, we also present in Fig. 4a a data set obtained with Feshbach molecules. Within 5% the laser intensities of the optical lattice are the same for all three data sets. Remarkably, the measurements clearly reveal that the decay of molecules in state v=R=0 takes place on a similar timescale as compared to the v=0, R=2 molecules or the highly excited Feshbach molecules. This is not obvious because the relaxation paths are potentially different for these states. Specifically, while the Feshbach molecules can vibrationally relax within the triplet potential An external file that holds a picture, illustration, etc.
Object name is ncomms14854-m28.jpg, our v=R=0 molecules are already energetically in the absolute lowest level of the triplet manifold, also with respect to the hyperfine and Zeeman structure (see Methods). Thus, in an inelastic or reactive collision of two of our v=R=0 molecules either a Rb trimer must form or at least one of the two dimers must undergo a spin flip towards the singlet electronic ground state. Nevertheless, judging from our measurements presented here, there is no indication for a suppression of the molecular loss rate due to these restrictions, which is an important result of our experiments. These findings go along with theoretical predictions for collisions of polar triplet molecules where a spin flip of the electronic state was not suppressed either44.

Figure 4
Decay curves for various molecular quantum states.

The step-like loss discussed earlier for the Feshbach molecules is also visible in the data on the deeply bound states, albeit less pronounced. The softening of the steps is caused by smaller initial wave packet widths An external file that holds a picture, illustration, etc.
Object name is ncomms14854-m29.jpg(t=0) of the v=0 molecules since their polarizabilities at a wavelength of λ=1,064 nm are by factors between two and three higher as compared to the Feshbach molecules45. This leads to a stronger lattice confinement for the same laser intensities and, in addition, to an earlier non-adiabatic release of the wave packets when the lattice is ramped down.

Tuning of reaction rates

We now investigate the dependence of the reaction rate coefficients on the confinement and collision energy. For this purpose, we measure decay curves for various trap frequencies ωr of the tubes. Generally, each of these trap frequencies corresponds to a different collision energy Ecol, because in our current setup we cannot tune ωr and Ecol independently as both are controlled via the laser intensity I of the optical lattice. Specifically, ωr[proportional, variant]An external file that holds a picture, illustration, etc.
Object name is ncomms14854-m30.jpg and the collision energy scales as Ecol[proportional, variant]An external file that holds a picture, illustration, etc.
Object name is ncomms14854-m31.jpg[proportional, variant]I due to the initial potential energies of the particles. The decay rate constants An external file that holds a picture, illustration, etc.
Object name is ncomms14854-m32.jpg are shown as black circles in Fig. 5 where we use the average trap frequency An external file that holds a picture, illustration, etc.
Object name is ncomms14854-m33.jpg as a scale for radial confinement (Note that the v=0, R=2 state exhibits different polarizabilities in x- and y-direction45).

Figure 5
Reaction rate coefficients An external file that holds a picture, illustration, etc.
Object name is ncomms14854-m77.jpg and An external file that holds a picture, illustration, etc.
Object name is ncomms14854-m78.jpg.

We gain additional insights about the dependencies of the rate coefficient from theoretical considerations. For this, we use a quantum defect reaction model where an inelastic process takes place at short range with probability Pre (refs 39, 40). On the basis of the experimental observation that the majority of the molecules is lost in their first encounter with another molecule, quite independently of the initially prepared molecular state, we expect the reaction probability Pre to be close to unity. In the limit that Pre is unity, we are in the universal regime and the resulting s-wave scattering length can be written in a simple form An external file that holds a picture, illustration, etc.
Object name is ncomms14854-m34.jpg (see Methods for the general formula with Pre≠1). Here An external file that holds a picture, illustration, etc.
Object name is ncomms14854-m35.jpg=2πR62(1/4) is the mean scattering length of the van der Waals potential46, Γ is the gamma function and R6=(2μC6ħ−2)1/4 where μ is the reduced mass of the molecules. The C6 coefficient is C6≈17,550 a.u. (C6≈18,800 a.u.) for the v=0 molecules (Feshbach molecules), respectively47. In free space these parameter values would correspond to universal reaction rate constants39 An external file that holds a picture, illustration, etc.
Object name is ncomms14854-m36.jpg cm3 s−1 and 1.35 × 10−10 cm3 s−1, respectively, which roughly agree with measured reaction rate constants for 87Rb2 Feshbach molecules15,24.

The description changes when reducing the dimensionality of the scattering process35,48. Generally, a system enters the quasi-1D regime for large trap aspect ratios An external file that holds a picture, illustration, etc.
Object name is ncomms14854-m37.jpg and low enough collision energies Ecol<2ħωx,y. For our experiments, we estimate the average maximal energies to be a factor of two to three below this boundary (see Supplementary Note 4 and associated Supplementary Fig. 3) and ωx,y/ωz is at least >300 for the three investigated molecular states. Thus, the system can be described by an effective 1D model characterized by the interaction potential An external file that holds a picture, illustration, etc.
Object name is ncomms14854-m38.jpg with complex 1D scattering length a1D (refs 49, 50), which is a function of An external file that holds a picture, illustration, etc.
Object name is ncomms14854-m39.jpg and the trap confinement (see Methods). Parametrizing a1D as a1D=(α1D1D)−1 allows for writing the universal 1D reaction rate constant in the form35

An external file that holds a picture, illustration, etc.
Object name is ncomms14854-m40.jpg

where ħk is the relative momentum of the colliding molecules. The rate constant An external file that holds a picture, illustration, etc.
Object name is ncomms14854-m41.jpg can be considered a function of the collision energy Ecol, since Ecol=ħ2k2/(2μ). In addition, through the connection between 1D and 3D scattering lengths, An external file that holds a picture, illustration, etc.
Object name is ncomms14854-m42.jpg also depends on the transverse confinement ωr, that is, ωx and ωy (see equation (6) in Methods). In the low-energy limit, k→0, the 1D rate constant An external file that holds a picture, illustration, etc.
Object name is ncomms14854-m43.jpg vanishes as k2, in contrast to 3D where rate constants generally approach a constant value48. This is a manifestation of the change of Wigner threshold laws under confinement. To compare scaling of the 1D and 3D cases in more detail we examine the decay rates An external file that holds a picture, illustration, etc.
Object name is ncomms14854-m44.jpg, where i[set membership]{1D, 3D}, and n1D represents a 1D density of molecules. In the low-energy limit, the ratio Γ1D3D is proportional to (kd)2 × d2/An external file that holds a picture, illustration, etc.
Object name is ncomms14854-m45.jpg × (n1D/d2)/n3D (refs 35, 48), where d is the transverse confinement size, as given by the harmonic oscillator length (see Methods). We note that n1D/d2 can be thought of as an equivalent 3D density of the confined system. The decay rate ratio strongly depends on both confinement and collision energy and indicates that at low enough temperatures and strong confinement the 1D gas would be much more stable than the 3D one.

Figure 5 displays the full dependencies of the universal An external file that holds a picture, illustration, etc.
Object name is ncomms14854-m46.jpg on collision energy and confinement as described by equation (2) for both the Feshbach and the v=0 molecules. We note, however, that the calculated values of An external file that holds a picture, illustration, etc.
Object name is ncomms14854-m47.jpg cannot directly be compared to the experimentally determined values of An external file that holds a picture, illustration, etc.
Object name is ncomms14854-m48.jpg for two reasons. First, the measured An external file that holds a picture, illustration, etc.
Object name is ncomms14854-m49.jpg are only given as a function of confinement ωr and the corresponding average collision energies are not constant. Second, we have to take into account the dynamics of the molecules and their oscillating energy distribution within the lattice. Therefore, we calculate an approximate theoretical An external file that holds a picture, illustration, etc.
Object name is ncomms14854-m50.jpg by time-averaging over An external file that holds a picture, illustration, etc.
Object name is ncomms14854-m51.jpg according to

An external file that holds a picture, illustration, etc.
Object name is ncomms14854-m52.jpg

An alternative approach for comparison to the experimental data is described in Supplementary Note 5 (see also Supplementary Fig. 4). In equation (3), the summations cover all possible colliding pairs of molecules (i, j) in the whole sample with their respective collision energy Ecol. The theoretical calculations for An external file that holds a picture, illustration, etc.
Object name is ncomms14854-m53.jpg assuming the universal reaction are shown in Fig. 5 as white solid lines. The universal model reproduces the overall slopes of the data, but generally underestimates the measured decay rate constants. This can be explained by a slight non-universal character of the inelastic collisions for all three investigated states as will be shown in the following.

To account for non-universality in our calculations, we use the full expression of a3D as given by equation (5) of the Methods section. The theoretical values of An external file that holds a picture, illustration, etc.
Object name is ncomms14854-m54.jpg are calculated, again using equations (2) and (3). The non-universal model introduces two free fit parameters, the short-range parameter s and the reaction probability Pre (see Methods). The white dashed lines in Fig. 5 show the results of fits of this non-universal model to the data. For the Feshbach and the v=0, R=0 molecules good agreement with the measurements is achieved for a wide range of parameter values of s at reaction probabilities Pre between 0.4 and 0.9. The enhancement of the reaction rate relative to the universal regime is due to a shape resonance caused by a near-threshold bound state. For the v=0, R=2 state the best agreement with the data is obtained when maximizing the calculated rate which is achieved for Pre≈0.8 and a negative s. However, still some discrepancy remains. These rotationally excited molecules feature a partial spatial alignment of the molecular axes45. Specifically, they collide with their axes pointing dominantly perpendicular to the longitudinal direction of the tubes, which is in contrast to the R=0 molecules where the axis distribution is isotropic. This leads to additional quadrupole-quadrupole interactions. We checked, however, that these interactions are far too weak to significantly influence the reaction rates.

Discussion

In conclusion, we presented the experimental determination of reaction rates for deeply bound triplet molecular states. In general, we observe that the majority of the molecules are already lost in the first encounter. We also demonstrate that the decay rate constant in 1D depends on the confinement strength and the relative energy of the molecules. This offers the possibility to tune the inelastic interaction by adjusting the trap parameters. In this context, it could also be insightful to extend the collision experiments to the quasi-2D and 3D regimes since theory predicts that the dimensionality influences the energy dependence of the reaction rates35,36,37,38. Our measurements indicate that inelastic collisions of Rb2 molecules are characterized by rate constants close to universality. It would be interesting to pursue the investigation further, by mapping out how the rate constant depends individually on the confinement strength and the relative energy of the molecules. Furthermore, stereochemical aspects can be investigated by adjusting the relative alignment of the molecular axes via the preparation of specific rotational states45. The presented experimental results together with the provided theoretical model have to be considered as a step that might pave the way for a fundamental understanding of ultracold chemical reactions and their spatiotemporal control.

Methods

Preparation of cold molecules

The molecules are created in an optical lattice which consists of three perpendicular, retro-reflected laser beams of wavelength λ=1,064 nm. For this, a thermal sample of roughly 106 87Rb atoms at a temperature of ~1 μK is loaded into this lattice such that a significant number of sites is doubly occupied. The atoms reside in the lowest Bloch band of the 3D optical lattice and the optical potential is deep enough such that tunnelling is strongly suppressed. To first approximation the atoms are normally distributed in configuration space [proportional, variant] An external file that holds a picture, illustration, etc.
Object name is ncomms14854-m55.jpg with the widths of An external file that holds a picture, illustration, etc.
Object name is ncomms14854-m56.jpg(t=0)≈(26, 27, 26) μm. The atomic occupation of each individual lattice site is described by Poisson statistics. The atoms are in the electronic ground state with the total angular momentum quantum numbers f=1 and mf=1. By using magnetic Feshbach association at 1007.4 G we create weakly bound s-wave molecules at lattice sites occupied with exactly two atoms41,51. All remaining atoms are removed in a subsequent purification step41, such that a pure molecular sample is obtained. The molecular cloud has a size of An external file that holds a picture, illustration, etc.
Object name is ncomms14854-m57.jpg(t=0)≈(22, 24, 23) μm after production and contains ~4.5 × 104 Feshbach molecules. The uncertainties for the given values are ~10%. Afterwards, for a large part of our experiments, the molecules are transferred from the Feshbach state to the vibrational ground state v=0 of the energetically lowest triplet potential An external file that holds a picture, illustration, etc.
Object name is ncomms14854-m58.jpg. This is done via stimulated Raman adiabatic passage (STIRAP)1,45 with a transfer efficiency of ~80%. We prepare the molecules either in one of two quantum levels: The first one is described by the quantum numbers R=0, I=3, F=2, mF=2 where I, F, mF denote the total nuclear spin, the total angular momentum and its projection, respectively. This level is energetically the absolute lowest level of the triplet state, see52. The second level is rotationally excited by two units of angular momentum and has the quantum numbers R=2, I=3, F=4, mF=4. Compared to the R=0 level it has ~2 GHz × h higher energy, see52. The clouds of the v=0 molecules (both R=0 and R=2) are smaller than the cloud of Feshbach molecules, that is, An external file that holds a picture, illustration, etc.
Object name is ncomms14854-m59.jpg(t=0)≈(19, 20, 18) μm, and their particle numbers range from 2.5 × 104 to 3.2 × 104. For the given values the uncertainties are again ~10%.

Measuring the molecule number and cloud size

To measure the total number N and the cloud size An external file that holds a picture, illustration, etc.
Object name is ncomms14854-m60.jpg of the molecules in the 1D tubes at a particular point in time the 3D lattice is quickly switched on. This locks the molecules in their current positions. If they are in a v=0 state we transfer them back to the Feshbach state using STIRAP and subsequently dissociate them by magnetically ramping back over the Feshbach resonance. The resulting atoms are suddenly released from the optical lattice and after a short (200 μs) time of flight the atom cloud is imaged using standard absorption imaging for another 200 μs. We note that in our imaging procedure the number of molecules is in general underestimated and the measured cloud size is too large. In Supplementary Note 6 (see also Supplementary Fig. 5) we discuss the underestimation of the particle number in detail and derive a correction, which is applied to all our measurements. The overestimated cloud size is a result of the limited resolution of ~4 μm of the imaging optics and of the expansion of the atomic cloud during time of flight, and during absorption as well as due to tunneling in the 3D optical lattice shortly before detection. This tunneling is enhanced because a sizeable fraction of molecules will be excited to higher Bloch bands when they are rapidly reloaded back into the 3D lattice. This leads to an overestimation of the molecular cloud size An external file that holds a picture, illustration, etc.
Object name is ncomms14854-m61.jpg of ~3 μm. Therefore the relative influence is largest on smaller clouds. While we present the originally measured values for An external file that holds a picture, illustration, etc.
Object name is ncomms14854-m62.jpg in the text (see also Fig. 2b), we used corrected values for the simulations.

Determination of trap frequencies

Both, the radial (ωx,y) and longitudinal (ωz) trap frequencies are determined from measurements with Feshbach molecules. We use modulation spectroscopy43 to obtain ωx and ωy, respectively, while ωz is inferred from the periodicity of the steps in the molecular decay curves (cf. Fig. 2). The corresponding trap frequencies for the deeply bound v=0 states are derived by comparing their known polarizabilities to the ones of the Feshbach molecules43,45. As a consistency check, we find agreement between the so-predicted and the experimentally observed periodicity of the decay steps for both v=0 states (Fig. 4).

Numerical integration of the rate equation

To numerically integrate equation (1) we use the random number generated distribution of molecules (Fig. 1b), which assigns an initial location χi(t=0) to each particle i. Next, we propagate the wave packets of all molecules in small time steps Δt. The decay probability of each molecule pair (i, j) during Δt is given by An external file that holds a picture, illustration, etc.
Object name is ncomms14854-m63.jpg. If an inelastic collision takes place, both involved particles are removed from the sample.

Scattering in quasi-1D geometry

Since in our setup An external file that holds a picture, illustration, etc.
Object name is ncomms14854-m64.jpg, the collision in the presence of the trap can be described within the pseudopotential approximation with (in general) energy-dependent 3D scattering length50,53. In this treatment the interaction potential is replaced by the regularized Dirac delta function and in addition an effective 1D model can be derived49, using An external file that holds a picture, illustration, etc.
Object name is ncomms14854-m65.jpg. The resulting 1D scattering length a1D is connected to the 3D one via54,

An external file that holds a picture, illustration, etc.
Object name is ncomms14854-m66.jpg

Here C is a numerical factor depending only on the transverse trap anisotropy ωx/ωy. For the R=0 state (and the Feshbach state) C=−ζ(1/2)≈1.46, where ζ denotes the Riemann zeta function, while for the R=2 state C≈1.57. The complex valued 3D scattering length can be written as40,

An external file that holds a picture, illustration, etc.
Object name is ncomms14854-m67.jpg

Here y is defined via Pre=4y/(1+y)2, where Pre denotes the short-range inelastic process probability. The parameter s is the value of the scattering length in the absence of inelastic collisions in units of An external file that holds a picture, illustration, etc.
Object name is ncomms14854-m68.jpg.

The short-range reaction probability Pre approaches unity in the limit y→1. Inserting y=1 into equation (5), one observes that the scattering length approaches a3D=An external file that holds a picture, illustration, etc.
Object name is ncomms14854-m69.jpg(1−i) regardless of the value of s. This means that the reaction dynamics becomes universal in the sense that it is independent of the short-range details of the potential, which normally determine the scattering length. For y=1, equation (4) can be rewritten as

An external file that holds a picture, illustration, etc.
Object name is ncomms14854-m70.jpg

Data availability

The data that support the findings of this study are available from Björn Drews on reasonable request.

Additional information

How to cite this article: Drews, B. et al. Inelastic collisions of ultracold triplet Rb2 molecules in the rovibrational ground state. Nat. Commun. 8, 14854 doi: 10.1038/ncomms14854 (2017).

Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Material

Supplementary Information:

Supplementary Figures, Supplementary Notes and Supplementary References.

Acknowledgments

We thank Tommaso Calarco for fruitful discussions and Eberhard Tiemann for valuable information. K.J. is grateful for the hospitality of Tommaso Calarco. This work was supported by the German Research Foundation (DFG) within the project DE 510/2-1, by the Alexander von Humboldt Foundation, by the Foundation for Polish Science within the START programme and by the National Science Centre (Poland) project 2015/17/B/ST2/00592.

Footnotes

The authors declare no competing financial interests.

Author contributions B.D. and M.D. performed the experiments. All authors contributed to the data analysis and to the writing of the paper.

References

  • Lang F., Winkler K., Strauss C., Grimm R. & Hecker Denschlag J. Ultracold triplet molecules in the rovibrational ground state. Phys. Rev. Lett. 101, 133005 (2008). [PubMed]
  • Ni K.-K. et al. . A high phase-space-density gas of polar molecules. Science 322, 231–235 (2008). [PubMed]
  • Danzl J. G. et al. . An ultracold high-density sample of rovibronic ground-state molecules in an optical lattice. Nat. Phys. 6, 265–270 (2010).
  • Molony P. K. et al. . Creation of ultracold 87Rb133Cs molecules in the rovibrational ground state. Phys. Rev. Lett. 113, 255301 (2014). [PubMed]
  • Takekoshi T. et al. . Ultracold dense samples of dipolar RbCs molecules in the rovibrational and hyperfine ground state. Phys. Rev. Lett. 113, 205301 (2014). [PubMed]
  • Park J. W., Will S. A. & Zwierlein M. W. Ultracold dipolar gas of fermionic 23Na40K molecules in their absolute ground state. Phys. Rev. Lett. 114, 205302 (2015). [PubMed]
  • Guo M. et al. . Creation of an ultracold gas of ground-state dipolar 23Na87Rb molecules. Phys. Rev. Lett. 116, 205303 (2016). [PubMed]
  • Krems R. V. Cold controlled chemistry. Phys. Chem. Chem. Phys. 10, 4079–4092 (2008). [PubMed]
  • Carr L. D., DeMille D., Krems R. V. & Ye J. Cold and ultracold molecules: science, technology and applications. New J. Phys. 11, 055049 (2009).
  • Żuchowski P. S. & Hutson J. M. Reactions of ultracold alkali-metal dimers. Phys. Rev. A 81, 060703 (R) (2010).
  • Quéméner G. & Julienne P. S. Ultracold molecules under control!. Chem. Rev. 112, 4949–5011 (2012). [PubMed]
  • Richter F. et al. . Ultracold chemistry and its reaction kinetics. New J. Phys. 17, 055005 (2015).
  • Staanum P., Kraft S. D., Lange J., Wester R. & Weidemüller M. Experimental investigation of ultracold atom-molecule collisions. Phys. Rev. Lett. 96, 023201 (2006). [PubMed]
  • Zahzam N., Vogt T., Mudrich M., Comparat D. & Pillet P. Atom-molecule collisions in an optically trapped gas. Phys. Rev. Lett. 96, 023202 (2006). [PubMed]
  • Syassen N. et al. . Strong dissipation inhibits losses and induces correlations in cold molecular gases. Science 320, 1329–1331 (2008). [PubMed]
  • Wang T. T., Heo M.-S., Rvachov T. M., Cotta D. A. & Ketterle W. Deviation from universality in collisions of ultracold 6Li2 molecules. Phys. Rev. Lett. 110, 173203 (2013). [PubMed]
  • Mukaiyama T., Abo-Shaeer J. R., Xu K., Chin J. K. & Ketterle W. Dissociation and decay of ultracold sodium molecules. Phys. Rev. Lett. 92, 180402 (2004). [PubMed]
  • Zirbel J. J. et al. . Collisional stability of fermionic Feshbach molecules. Phys. Rev. Lett. 100, 143201 (2008). [PubMed]
  • Chin C. et al. . Observation of Feshbach-like resonances in collisions between ultracold molecules. Phys. Rev. Lett. 94, 123201 (2005). [PubMed]
  • Regal C. A., Greiner M. & Jin D. S. Lifetime of molecule-atom mixtures near a Feshbach resonance in 40K. Phys. Rev. Lett. 92, 083201 (2004). [PubMed]
  • Jochim S. et al. . Pure gas of optically trapped molecules created from fermionic atoms. Phys. Rev. Lett. 91, 240402 (2003). [PubMed]
  • Ferlaino F. et al. . Collisions between tunable halo dimers: exploring an elementary four-body process with identical bosons. Phys. Rev. Lett. 101, 023201 (2008). [PubMed]
  • Cubizolles J., Bourdel T., Kokkelmans S. J. J. M. F., Shlyapnikov G. V. & Salomon C. Production of long-lived ultracold Li2 molecules from a Fermi gas. Phys. Rev. Lett. 91, 240401 (2003). [PubMed]
  • Syassen N., Volz T., Teichmann S., Dürr S. & Rempe G. Collisional decay of 87Rb Feshbach molecules at 1005.8 G. Phys. Rev. A. 74, 062706 (2006).
  • Knoop S. et al. . Observation of an Efimov-like trimer resonance in ultracold atom-dimer scattering. Nat. Phys. 5, 227–230 (2009).
  • Zaccanti M. et al. . Observation of an Efimov spectrum in an atomic system. Nat. Phys. 5, 586–591 (2009).
  • Pollack S. E., Dries D. & Hulet R. G. Universality in three- and four-body bound states of ultracold atoms. Science 326, 1683–1685 (2009). [PubMed]
  • Bloom R. S., Hu M.-G., Cumby T. D. & Jin D. S. Tests of universal three-body physics in an ultracold Bose-Fermi mixture. Phys. Rev. Lett. 111, 105301 (2013). [PubMed]
  • Nakajima S., Horikoshi M., Naidon P. & Ueda M. Nonuniversal Efimov atom-dimer resonances in a three-component mixture of 6Li. Phys. Rev. Lett. 105, 023201 (2010). [PubMed]
  • Lompe T. et al. . Atom-dimer scattering in a three-component Fermi gas. Phys. Rev. Lett. 105, 103201 (2010). [PubMed]
  • Ni K.-K. et al. . Dipolar collisions of polar molecules in the quantum regime. Nature 464, 1324–1328 (2010). [PubMed]
  • Ospelkaus S. et al. . Quantum-state controlled chemical reactions of ultracold potassium-rubidium molecules. Science 327, 853–857 (2010). [PubMed]
  • de Miranda M. H. G. et al. . Controlling the quantum stereodynamics of ultracold bimolecular reactions. Nat. Phys. 7, 502–507 (2011).
  • Yan B. et al. . Observation of dipolar spin-exchange interactions with lattice-confined polar molecules. Nature 501, 521–525 (2013). [PubMed]
  • Idziaszek Z., Jachymski K. & Julienne P. S. Reactive collisions in confined geometries. New J. Phys. 17, 035007 (2015).
  • Li Z. & Krems R. V. Inelastic collisions in an ultracold quasi-two-dimensional gas. Phys. Rev. A 79, 050701 (R) (2009).
  • Li Z., Alyabyshev S. V. & Krems R. V. Ultracold inelastic collisions in two dimensions. Phys. Rev. Lett. 100, 073202 (2008). [PubMed]
  • Petrov D. S. & Shlyapnikov G. V. Interatomic collisions in a tightly confined Bose gas. Phys. Rev. A 64, 012706 (2001).
  • Idziaszek Z. & Julienne P. S. Universal rate constants for reactive collisions of ultracold molecules. Phys. Rev. Lett. 104, 113202 (2010). [PubMed]
  • Jachymski K., Krych M., Julienne P. S. & Idziaszek Z. Quantum theory of reactive collisions for 1/rn potentials. Phys. Rev. Lett. 110, 213202 (2013). [PubMed]
  • Thalhammer G. et al. . Long-lived Feshbach molecules in a three-dimensional optical lattice. Phys. Rev. Lett. 96, 050402 (2006). [PubMed]
  • Kinoshita T., Wenger T. & Weiss D. S. A quantum Newton's cradle. Nature 440, 900–903 (2006). [PubMed]
  • Deiß M. et al. . Polarizability of ultracold Rb2 molecules in the rovibrational ground state of . New J. Phys. 17, 065019 (2015).
  • Tomza M., Madison K. W., Moszynski R. & Krems R. V. Chemical reactions of ultracold alkali-metal dimers in the lowest-energy 3Σ state. Phys. Rev. A 88, 050701 (R) (2013).
  • Deiß M., Drews B., Deissler B. & Hecker Denschlag J. Probing the axis alignment of an ultracold spin-polarized Rb2 molecule. Phys. Rev. Lett. 113, 233004 (2014). [PubMed]
  • Gribakin G. F. & Flambaum V. V. Calculation of the scattering length in atomic collisions using the semiclassical approximation. Phys. Rev. A 48, 546–553 (1993). [PubMed]
  • Tomza M. et al. . Interatomic potentials, electric properties and spectroscopy of the ground and excited states of the Rb2 molecule: ab initio calculations and effect of a non-resonant field. Mol. Phys. 111, 1781–1797 (2013).
  • Micheli A. et al. . Universal rates for reactive ultracold polar molecules in reduced dimensions. Phys. Rev. Lett. 105, 073202 (2010). [PubMed]
  • Olshanii M. Atomic scattering in the presence of an external confinement and a gas of impenetrable bosons. Phys. Rev. Lett. 81, 938–941 (1998).
  • Idziaszek Z. & Calarco T. Two atoms in an anisotropic harmonic trap. Phys. Rev. A 71, 050701 (R) (2005).
  • Köhler T., Góral K. & Julienne P. S. Production of cold molecules via magnetically tunable Feshbach resonances. Rev. Mod. Phys. 78, 1311–1361 (2006).
  • Strauss C. et al. . Hyperfine, rotational, and vibrational structure of the state of 87Rb2. Phys. Rev. A 82, 052514 (2010).
  • Bolda E. L., Tiesinga E. & Julienne P. S. Effective-scattering-length model of ultracold atomic collisions and Feshbach resonances in tight harmonic traps. Phys. Rev. A 66, 013403 (2002).
  • Peng S.-G., Bohloul S. S., Liu X.-J., Hu H. & Drummond P. D. Confinement-induced resonance in quasi-one-dimensional systems under transversely anisotropic confinement. Phys. Rev. A 82, 063633 (2010).

Articles from Nature Communications are provided here courtesy of Nature Publishing Group