# Suppression of molecular decay in ultracold gases without Fermi statistics

## Abstract

We study inelastic processes for ultracold three-body systems in which only one interaction is resonant. We have found that the decay rates for weakly bound molecules due to collisions with other atoms can be suppressed not only without fermionic statistics but also when bosonic statistics applies. In addition, we show that at ultracold temperatures three-body recombination involving a single resonant pair of atoms leads mainly to formation of weakly bound molecules which, in turn, are stable against decay. These results indicate that recombination in three-component atomic gases can be used as an efficient mechanism for molecular formation, allowing the achievement of high molecular densities.

###### pacs:

34.50.-s,34.50.Cx,67.85.-d,31.15.xIn recent years, the efficiency of diatomic molecule formation in ultracold quantum gases and their stability once formed have become key ingredients for many experiments. For example, the remarkable stability of weakly bound molecules in ultracold gases of fermions in two different spin states LongLivedMol () has greatly helped experimental studies in the BEC-BCS crossover regime BECBCS-Studies () when tuning the two-body -wave scattering length from to near a Feshbach resonance. The long lifetimes recently observed for ultracold heteronuclear molecules RbK () might also pave the way for future studies on ultracold polar molecules. The stability of weakly bound molecules where at least one atom is a fermion has been interpreted as a result of suppression due to the fermionic statistics governing both atom-molecule and molecule-molecule collisions Petrov (); OurPapers (). When bosonic statistics plays a role, however, weakly bound molecules tend to decay rapidly, limiting the experimental possibilities for studies near a Feshbach resonance OurPapers (); Paul (); BraatenVrel (). Obviously, besides stability, efficient production of molecules is also of crucial experimental importance.

In this Letter, we present a series of new predictions that opens up new possibilities for obtaining long molecular lifetimes and that also provide an alternative path for efficient molecular formation. Surprisingly, we have found that weakly bound molecules of two resonantly interacting atoms can still be stable against collisions with other atoms even when bosonic statistics applies. Specifically, inelastic collisions between these weakly bound molecules and any atom not resonant with either of the molecule’s atoms is suppressed as , therefore allowing long molecular lifetimes near a Feshbach resonance. This suppression in the absence of fermionic statistics has recently been observed in Ref. RbK () in an ultracold mixture of weakly bound RbK molecules and K atoms in an nonresonant state. Here, we show that the mechanism that determines the suppression can be traced to the Efimov physics governing the atom-molecule interactions OurPapers (). In contrast to the cases with fermionic suppression Petrov (); OurPapers (), our predicted suppression holds for any mass ratio between the collision partners. The sole requirement is that only the atoms bound in the molecule interact resonantly.

We will also show that three-body recombination in such systems displays interesting and potentially useful behavior. In the one- and two-component bosonic gases and in bose-fermi mixtures studied to date, recombination can lead to a substantial fraction of deeply bound molecules OurPapers (); OurK3Papers (); BraatenReview (), giving them in the process enough kinetic energy to escape from typical traps. Moreover, any weakly bound molecule formed can eventually decay rapidly due to collisions with other atoms and molecules K3RudiBosons (). The situation is better for fermionic gas mixtures of two spin species: molecules are stable Petrov (); OurPapers () and the recombination rate into weakly bound molecules () dominates recombination into deeply bound molecules () OurPapers (). These two ingredients ensured the achievement of highly efficient molecular formation and long molecular lifetimes K3Fermions (). However, because recombination for fermionic systems scales as , efficient molecule formation is restricted to relatively high temperatures or large scattering lengths K3Fermions (). In this regime, though, collisions can potentially lead to high molecular dissociation rates ( OurPapers ()), limiting the molecular lifetimes.

As we will show, a three-component gas of atoms with only a single resonant interaction retains the advantages of these fermionic systems — efficient formation of weakly bound molecules and long lifetimes — with the additional advantage of having a non-zero recombination rate at , allowing the formation of molecules with much lower temperature which is important to explore the quantum degenerate regime. We have found that for these systems three-body recombination into weakly bound molecules scales as at ultracold temperatures while recombination into deeply bound molecules scales as . This weaker -dependence ensures that near a Feshbach resonance recombination will lead mainly to formation of weakly bound molecules. As shown in Refs. K3RudiBosons (); K3Fermions (), these weakly bound molecules can stay trapped due to the small energy released, minimizing atomic and molecular losses. The most important aspect of molecular formation using this mechanism, however, is that even in the absence of fermionic statistics, the molecules formed are stable against decay since the decay rate scales as . This combination of higher formation rate and lower decay rate maximizes molecule production.

We study these ultracold three-body collisions by solving the Schrödinger equation in the adiabatic hyperspherical representation. The three-body inelastic collision rates are determined by solving the hyperradial Schrödinger equation given by (in atomic units)

(1) |

Here, is the three-body reduced mass, is the total energy, the hyperradial wave function, and a collective index that represents all quantum numbers necessary to label each channel. Since the hyperradius gives the overall size of the system, this equation describes the collective radial motion under the influence of the effective potential , with inelastic transitions driven by the nonadiabatic couplings OurPapers ().

In the regime where the two-body interactions are resonant, i.e., when (where is the characteristic range of two-body interactions), the effective potentials become universal, and several analytical properties of ultracold three-body systems can be derived. In fact, in our framework OurPapers () both the and the dependence of the three-body collision rates can be determined by simply identifying the attractive or repulsive character of as determined by the influence of Efimov physics OurPapers (). Previous work in ultracold three-body collisions Petrov (); OurPapers (); BraatenVrel (); BraatenReview () has analyzed the influence of Efimov physics for systems where at least two of the pairwise two-body interactions are resonant. In this paper, we will show that Efimov physics still has an important impact even when just one interaction is resonant.

For systems with single resonant interaction, the effective potentials are repulsive in the range and can be conveniently parametrized by the coefficients and OurPapers ():

(2) |

For , is associated with the atom-dimer channel and therefore plays an important role in atom-molecule processes. The effective potentials associated with describe three-body continuum channels. The coefficients and in Eq. (2) can be determined by assuming a zero-range model potential for the interatomic interactions (we have confirmed that finite-range two-body interactions give the same results) and writing the three-body wave function in terms of the Faddeev components OurPapers (). Doing this, we have determined that and are

(3) |

where is the total orbital angular momentum. Interestingly, in contrast to systems where two or three interactions are resonant, the strengths of the potentials and do not depend on either the mass ratio or the permutational symmetry of the resonant pair.

Our derivation of these rates’ scaling behavior rests on the observation that the rate limiting step is tunneling through the potential barriers represented by Eq. (2) OurPapers (). This, together with and from Eq. (3), gives the scaling laws shown in Table 1 for systems with one resonant pair of atoms. We show the results for systems (three distinguishable atoms) and for systems (two identical bosons and a third distinguishable atom). For systems, we have assumed that the - pair interacts resonantly and the atom is in an nonresonant state. For , we assumed that only the intraspecies interaction is resonant and, therefore, the atom is in an nonresonant state. Note that the scaling laws for both and systems are exactly the same. In the table, we show the first three partial wave contributions to vibrational relaxation of weakly bound molecules, (); to three-body recombination into weakly bound molecules, (); to recombination into deeply bound molecules, (); and to dissociation of weakly bound molecules, ().

(), | ||||||
---|---|---|---|---|---|---|

const | const | , | ||||

, | ||||||

, |

Figure 1(a) shows our numerical results for , obtained from a direct solution of Eq. (1) (see Ref. EsryHelium () for details) for three distinguishable atoms, , for ( is the total parity). We have assumed a finite-range model two-body interaction in which the resonant - interaction is adjusted to support one weakly bound -wave state whose position controls the scattering length, one deeply bound -wave state, and one deeply bound -wave state. Our model - interaction supports only one -wave state, and, for simplicity, we assumed the - interaction has no bound state. The figure shows that when , our predicted scaling law is verified by the numerical results. Thus, we have clear evidence that vibrational relaxation is indeed suppressed for large even in the absence of fermionic statistics. By extension, this result also support our arguments that molecular decay in collisions are also suppressed as even though bosonic statistics apply. In a two-component gas of and atoms, however, long lifetimes for molecules can only be achieved if there are no free atoms and with the additional condition of low molecular density, to prevent decay due to molecule-molecule collisions.

Recombination for systems involving nonresonant atoms also presents some peculiar properties. As mentioned previously, for , recombination into weakly bound molecules () scales as while recombination into deeply bound molecules () scales as (see Table 1). Therefore, recombination for large leads mostly to the formation of weakly bound molecules. This differential scaling with is in stark contrast to previously studied systems with two or three resonant interactions for which recombination is non-zero at , as is the case for one- and two-component bosonic systems and boson-fermion mixtures. In these systems, recombination into all two-body states scales as . In those systems, recombination into the weakly bound state is typically dominant as well, but only via the coefficient of . Systems with only a single resonant interaction share this property in addition to the differential scaling with . Figure 1(b) shows our numerical calculations of recombination rates for the same system as in Fig. 1(a). The numerical results support all of the above predictions, including not only the smaller magnetude of relative to but also their scaling with . Another consequence of this differential scaling is that for , i.e., in the absence of weakly bound molecules, recombination is proportional to . So, recombination at for systems with a single resonant pair could be suppressed compared to for large .

Gas Mixture | , | , | , | ||||||
---|---|---|---|---|---|---|---|---|---|

-- | , | , | |||||||

, | , | ||||||||

, | , | ||||||||

-- | , | , | , | ||||||

, | , | , | |||||||

, | , | , | |||||||

-- | , | , | , | ||||||

, | , | , | |||||||

, | , | , | |||||||

-- | , | ||||||||

, | |||||||||

, |

The results obtained here for three-body processes involving a single resonant pair — summarized in Table 1 — combined with previous results OurPapers () exhaust all possible three-body processes relevant to three-component atomic gases. The competition between the various collision processes in a multi-component gas is dictated fundamentally by their energy and scattering length dependence, although it is also crucial to take into account the atomic and molecular densities. In fact, manipulating the atomic and molecular densities can make favorable a particularly desired collisional behavior.

For completeness, we show in Table 2 the three-body collision rates for all possible ultracold three-component atomic gases. For each mixture, we have listed the three-body processes according to which interatomic interaction is resonant by indicating the resonant pair of atoms, and we have only included the dominant contribution at threshold. Three-body processes involving a single resonant pair are indicated by “” while all the others were obtained in previous work Petrov (); OurPapers (); BraatenVrel (); BraatenReview (). Futher, we have for and for , depending on the mass ratio between the collision partners OurPapers (), and we have omitted modulation factors that produce the interference minima and resonant peaks related to Efimov physics OurPapers ().

We can see from the table that in all possible three-component mixtures weakly bound molecules can be expected to be long-lived with relaxation scaling as , provided that all free atoms are nonresonant. In fact, for gas mixtures where distinguishable fermions are resonant, weakly bound molecules are expected to be long-lived irrespective of the identity of the free atoms.

Given the suppression possible in atom-molecule collisions, it might be that molecule-molecule collisions are the dominant relaxation mechanism in a three-component gas. For molecules, that are themselves fermions, molecule-molecule collisions are suppressed at ultracold temperatures due to their -wave character. For bosonic molecules, however, long lifetimes can only be reached for low molecular densities, in order to prevent molecule-molecule collisions. In contrast, for bosonic molecules, relaxation due to molecule-molecule collisions scales as with Petrov () for mass ratios (a condition satisfied for commonly used alkali atoms) and the main molecular collisional decay is due to atom-molecule collisions involving a single resonant interaction, scaling as .

From Table 2 we see that the optimal three-component systems utilize molecular formation close to an - Feshbach resonance. In this case, recombination involving two identical fermions is suppressed at ultracold temperatures while recombination in both -- and -- mixtures leads mainly to the formation of weakly bound molecules. The most important property of this system is that after the weakly bound molecules are formed, molecular loss due to collisions with the remaining atoms is suppressed at least as . Similar benefits can also be realized near a - resonance under certain circumstances. For instance, in -- mixtures, efficient - molecular formation could be achieved if the number of nonresonant atoms is higher than the and numbers. If in the end all bosons are bound in molecules, then collisions with the remaining fermions is suppressed, ensuring the desired stability of weakly bound molecules.

To summarize, we have shown that vibrational relaxation of weakly bound molecules can be suppressed near a Feshbach resonance even in the absence of Fermi statistics. Remarkably, the root cause of this suppression can be traced to the same kind of universality that yields the Efimov effect. The systems that show this behavior can also exhibit differential scaling with of the recombination rate into weakly bound molecules that is finite at zero temperature. These two results combine to make ultracold three-component gases with a single resonant interaction an attractive alternative for efficiently producing stable, ultracold molecules.

This work was supported by the National Science Foundation, the Air Force Office of Scientific Research, and by the W. M. Keck Foundation.

### References

- M. Greiner, C. A. Regal, and D. Jin, Nature (London) 426, 537 (2003); J. Cubizolles et al., Phys. Rev. Lett. 91, 240401 (2003); C. A. Regal, M. Greiner, and D. S. Jin Phys. Rev. Lett. 92, 083201 (2004); S. Jochim et al., Science 302 2101 (2003); M. W. Zwierlein et al., Phys. Rev. Lett. 92, 120403 (2004); K. E. Strecker, G. B. Partridge, and R. G. Hulet, ibid., 91 080406 (2003).
- C. A. Regal, M. Greiner, and D. S. Jin, Phys. Rev. Lett. 92, 040403 (2004); M. W. Zwierlein et al., ibid. 92, 120403 (2004); C. Chin et al., Science 305, 1128 (2004); M. W. Zwierlein, J. R. Abo-Shaeer, A. Schirotzek, C. H. Schunck, and W. Ketterle, Nature (London) 435, 1047 (2005).
- J. J. Zirbel, K.-K. Ni, S. Ospelkaus, J. P. D’Incao, C. E. Wieman, J. Ye, and D. S. Jin, cond-mat.other/0710.2479v1
- D. S. Petrov, C. Salomon, and G. V. Shlyapnikov, Phys. Rev. Lett. 93, 090404 (2004); J. Phys. B 38, 645 (2005).
- J. P. D’Incao and B. D. Esry, Phys. Rev. Lett. 94, 213201 (2005).
- V. A. Yurovsky, A. Ben-Reuven, P. S. Julienne, and C. J. Williams, Phys. Rev. A 60, R765 (1999).
- E. Braaten and H.-W. Hammer, Phys. Rev. A 70, 042706 (2004).
- J. P. D’Incao, H. Suno, and B. D. Esry, Phys. Rev. Lett. 93, 123201 (2004); J. P. D’Incao, and B. D. Esry, Phys. Rev. A 73, 030702(R) (2006); cond-mat/0703269.
- E. Braaten and H.-W. Hammer, Phys. Rep. 428, 259 (2006).
- T. Weber, et al., Phys. Rev. Lett. 91, 123201 (2003).
- J. Cubizolles, et al., Phys. Rev. Lett. 91, 240401 (2003); S. Jochim, et al., Phys. Rev. Lett. 91, 240402 (2003); M. W. Zwierlein, et al., Phys. Rev. Lett. 91, 250401 (2003).
- H. Suno, B. D. Esry, C. H. Greene and J. P. Burke, Jr., Phys. Rev. A 65, 042725 (2002).