# Kinetic Model of Trapped Finite Temperature Binary Condensates

## Abstract

We construct a non-equilibrium theory for the dynamics of two interacting finite-temperature atomic Bose-Einstein condensates and use it to numerically estimate the relative rates of the arising collisional processes near equilbrium. The condensates are described by dissipative Gross-Pitaevskii equations, coupled to quantum Boltzmann equations for the thermal atoms. The density-density interactions between atoms in different components facilitate a number of transport processes of relevance to sympathetic cooling: in particular, considering realistic miscible and immiscible trapped atomic Rb-K and Rb-Rb condensate mixtures, we highlight the dominance of an inter-component scattering process associated with collisional “exchange” of condensed and thermal atoms between the components close to equilibrium.

###### pacs:

03.75.Mn,67.85.-dIntroduction. The unprecedented control of trapped neutral cold atom experiments enables the creation and study of degenerate multi-component systems, including Bose-Bose (1), Bose-Fermi (2) and Fermi-Fermi (3); (4) mixtures, and the related problems of spinor gases (5) and artificial gauge fields (6). In two-component systems, cooling to quantum degeneracy is typically performed through a combination of evaporative (7) and sympathetic (8) cooling techniques. For bosonic mixtures, this leads to condensation in differing hyperfine states of an atom (Rb (8); (1); (9); (10)), differing isotopes (e.g. Rb-Rb (11), Yb-Yb (12)) or differing elements (e.g. Rb-K (13); (14), Rb-Cs (15); (16), Rb-Sr and Rb-Sr (17), Rb-Na (18)). Interest in such systems has focused on understanding numerous properties, including equilibrium profiles (19); (20); (21); (22); (23), stability properties (24); (25), collective excitations (20); (26), vortices (27), solitary waves (28), and dissipative (29) and quenched dynamics (30); these effects depend critically on the relative inter-atomic interactions within and between the species, which determine whether the emerging condensates overlap spatially or phase-separate (19); (20).

The description of bosonic mixtures is typically either focused on the low temperature (Gross-Pitaevskii) (31), or high-temperature (Boltzmann) limit (32), or on treating the condensate in contact with a static heat bath (33); (34); (35). Alternative approaches are based on classical field methods (which ignore the dynamics of the high-lying thermal modes (36); (30)), or on number-conserving methods (which explicitly include only the back-action of the thermal cloud on the condensate (37); (38)). Moreover, although related kinetic models have been derived in the context of spinor gases (39); (40), there has been to date no critical assessment of the relative importance of the various collisional processes at finite temperatures.

In this work we (i) present a self-consistent model for the study of partly-condensed bosonic mixtures in the presence of dynamical condensates and thermal clouds (shown schematically in Fig. 1), and (ii) perform detailed numerical calculations at equilibrium to shed light into the relative importance of those processes for experimentally-accessible miscible and immiscible mixtures. A key conclusion is the estimation of the importance of a novel “cross-condensate-exchange” term found to dominate as soon as both species exhibit condensation.

Our approach is a generalization of the “Zaremba–Nikuni–Griffin” (ZNG) kinetic model (41); (42); (43) which, in the context of single-component condensates, has proven extremely useful in describing collective modes (44), condensate growth (45), dynamics of macroscopic excitations (46) and surface evaporative cooling (47). This approach is based on the decomposition of the Bose field operator into a symmetry-breaking part and its fluctuations, and a separation of timescales argument (see also (48); (49); (50)): in our analysis, we explicitly separate slowly-evolving variables for the condensate and thermal clouds, whose non-vanishing mean values are defined by appropriate dynamical equations, from those evolving on the more rapid collisional timescales. A careful self-consistent perturbative treatment of the latter in the context of an appropriate perturbing Hamiltonian (details to appear elsewhere) enables us to quantify for the first time collisional processes not accounted for in earlier models of Refs. (39); (40).

Coupled Kinetic Theory. The second-quantized Hamiltonian describing the interacting binary mixture is

(1) |

where is the annihilation operator for a species- atom, and , where is the scattering length between atoms in component and , is the reduced mass, is the mass of a boson in component ; and is the trapping potential.

The condensed and non-condensed components are separated by the Beliaev decomposition , with condensate evolution equations obtained from the Heisenberg equations of motion for as

(2) |

where denotes averages (with ), , and is the effective condensate potential. Following established techniques (41), we identify the condensate fields (corresponding condensate densities ), and diagonal non-condensate densities as the only relevant slowly-varying quantities of interest. Triplet anomalous averages , are treated perturbatively (via adiabatic elimination (48); (49)), giving rise to well-known source terms, and we only maintain energy-conserving contributions, thus also discarding pair anomalous averages , (51). An important aspect of our work (see also Refs. (36); (38)) is that we explicitly maintain perturbative corrections arising from the off-diagonal terms . As a result, the equation of motion for component takes the form

(3) |

It contains a number of source terms: appearing in Eq. (3) describes the intra-component scattering of condensate and non-condensate particles and is well-known from single-component ZNG theory (41); is a direct generalization of to binary condensates, describing the inter-component scattering of condensed and non-condensed atoms. , is qualitatively different from the other two (see Eq. (6) below). , and are defined through their relation to the collision integrals () and given below.

To construct a self-consistent theory including the evolution of the non-condensate, the multi-component single-particle Wigner distribution is introduced as . We restrict our analysis to the diagonal terms of , i.e. , valid in the absence of coherent couplings between states. The kinetic equation for component is

(4) |

The Hartree-Fock potential that non-condensed atoms in component feel is , and the non-condensate density of component is obtained via . The source terms appearing in Eq. (3) are related to the collision integrals by the relationships (, ) and (). The various collisional contributions to the kinetic scattering of particles for the component are summarized in Fig. 2. Diagrams - concern collision integrals describing the scattering of condensed and non-condensed particles, while diagrams and show scattering amongst non-condensed atoms.

The collision integral (encapsulating ) appearing in Eq. (4) is defined by (using )

(5) |

It describes the scattering of a condensate atom and a non-condensate atom into thermal states, and its inverse process (Fig. 2(i)-(ii)), with the Kronecker delta accounting for enhanced scattering of atoms within the same species ( subcase). Here is the Hartree-Fock energy, defines the condensate energy, and gives the condensate momentum.

The collision integral is defined from the average of an off-diagonal pair of fluctuation operators, (), (in contrast to all other “ZNG” source terms originating from triplet terms (39); (40); (41)) as

(6) |

Although physically intuitive, this term is qualitatively different to the collision integral of Eq. (5), as it describes a process (Fig. 2(iii)) whereby one condensate and one non-condensate atom from different components scatter into a thermal and condensed state respectively.

The final collision processes and in Eq. (4) describe the scattering between non-condensate atoms of same () or differing () species, given by

(7) |

Numerical Results. To gain insight into the relative importance of the collision rates we compute their temperature dependence for experimentally-relevant equilibrium Rb-K and Rb-Rb mixtures in isotropic harmonic traps (Hz), and a total atom number in each component. These mixtures were chosen as their tunable scattering lengths (13); (14); (11) enable the probing of both miscible and immiscible () regimes. By re-expressing (and analogously for and ), we explicitly identify “in” and “out” scattering rates, which are equal at equilibrium. Following Refs. (43); (41), we define collisional rates that give the number of atoms leaving a phase-space volume per unit time as a result of collisions, for a perturbation from equilibrium. By transforming to the centre of mass frames, the collision rates can be written as where is the relative velocity and the cross-section. The collision rates between non-condensate atoms are while the collisions scattering rate, , takes the form

(8) |

with the relative velocity and .

Figure 3 shows the equilibrium condensate/thermal density profiles (top panels) and collision rates (middle/bottom panels) for a mixture of Rb and K at a temperature of , when condensate fractions , for both miscible (, left column) and immiscible (, right column) cases. The two condensates (dashed lines) mix (top panel, left) or phase-separate (top panel, right) with thermal clouds (solid lines) displaying peaks at the condensate edges, thus also leading to a mean-field-induced double-peaked thermal structure. The spatially-resolved collisional rates between condensate and non-condensate (middle panels) reveal peaks close to the condensate edges, which become more pronounced for immiscible condensates. In particular, collisional rates feature large localized peaks in regions where both components exhibit an appreciable condensate, which can locally dominate all other collisional processes in the centre of the immiscibility region. Collision rates between non-condensate atoms (bottom panels) are found to closely follow the shape of the thermal profiles, with the inter-species rates, (black-dotted line), affected by both thermal cloud distributions.

The cumulative effect of such collisional processes is best characterized through the integrated collision rates of , and over the cloud’s volume. These are shown in Fig. 4(a) for Rb-K (left panels), and Rb-Rb mixtures (right panels), in their experimentally-accessible miscibility (circles) and immiscibility (triangles) regimes for nK (Fig. 3 data, filled symbols) and nK (open symbols),. We find: (i) the collisions (green symbols) are at least as large as the other collisional terms in all cases; (ii) increasing the temperature (hollow to filled symbols) in the presence of an appreciable condensate fraction, enhances all collisional rates involving thermal atoms (including and ); (iii) collisional rates are largely sensitive to the interation strength (through prefactors), controllable through Feshbach resonances: an increase in the inter-species scattering length (13); (14) (left image) increases the relative importance of all Rb-K thermal-thermal, or condensate-thermal collisional terms (black/green points). Likewise, an increase in the tunable intra-component scattering length (11) (right), enhances the Rb-Rb (blue) collisional terms.

The significance of the collisional process becomes evident in Fig. 4(b), showing the temperature dependence of its ratio to the corresponding integrated rates for Rb-K mixtures. Since the efficiency of sympathetic cooling depends on the energy exchange rate, the rate of interspecies collisions (32) demonstrates the relevant role that the collisional terms may play

Rb | Rb | Rb | Rb | Rb | Rb | |

Rb | Rb | Rb | Rb | Rb | Rb | |

1 | ||||||

1 | ||||||

1 | - | - | 1 | - | - |

during various sympathetic cooling stages.

The evaluation of integrated collision rates enables us to define a typical collisional timescale , as , where is the relevant number of available thermal atoms taking part in collisions for each process; this allows us to distinguish between the collisionless [] and the hydrodynamic [] regimes (42); (31). Fig. 4(c) shows the variation of the hydrodynamic parameter for Rb-K mixtures in the miscible (left) and immiscible (right) regimes, focussing only on interspecies collisions: this reveals the dominance of , which approaches 1 in the immiscible case, with all intra-species collisions satisfying in both cases. Doubling the trap frequency to Hz roughly doubles , increases by and only slightly increases , which however remains the largest. Consideration of realistic Rb-Rb mixtures [characteristic values shown in Table 1] reveals that while is always the largest interspecies contribution, tuning the intraspecies to 900 (11) provides a plausible candidate for multi-component hydrodynamic behaviour in the sense of .

Conclusions. We have presented a system of coupled kinetic equations for the evolution of two-component condensates in the presence of dynamical thermal clouds,. Analytical considerations led to the identification of a “cross-condensate-exchange” term, which was numerically shown to dominate, even at relatively high temperatures, over other collisional terms close to equilibrium, both for Rb-K and Rb-Rb mixtures. Such a term could be highly relevant during sympathetic cooling, particularly in later stages, when the sympathetically-cooled component acquires appreciable condensation. Consideration of collisional timescales indicated the potential of generating hydrodynamic multi-component condensates with , at least for some collisional processes. Our work, which is generalizable to multi-component Bose gases and Bose-Fermi mixtures, sets the scene for future dynamical studies of sympathetic cooling, multi-component condensate formation and coupled expansion dynamics.

Acknowledgements. We acknowledge discussions with S.L. Cornish, S.A. Gardiner, P. Mason, E. Zaremba and M. Edwards, and support from EPSRC (Doctoral Prize Fellowship (MJE), grant EP/K03250X/1 (KLL, NPP)).

### References

- D. S. Hall, M. R. Matthews, J. R. Ensher, C. E. Wieman, and E. A. Cornell, Phys. Rev. Lett. 81, 1539 (1998).
- Z. Hadzibabic, C. A. Stan, C. A. Dieckmann, S. Gupta, M. W. Zwierlein, A. Görlitz, and W. Ketterle, Phys. Rev. Lett. 88, 160401 (2002).
- B. DeMarco and D. S. Jin, Science 285, 1703 (1999).
- M. Taglieber, A.-C. Voigt, T. Aoki, T. W. Hänsch, and K. Dieckmann, Phys. Rev. Lett. 100, 010401 (2008).
- D. M. Stamper-Kurn and M. Ueda, Rev. Mod. Phys. 85, 1191 (2013).
- J. Dalibard, F. Gerbier, G. Juzeliūnas, and P. Öhberg, Rev. Mod. Phys. 83, 1523 (2011).
- W. Ketterle and N. J. Van Druten, Adv. Atom. Mol. Opt. Phys. 37, 181 (1996).
- C. J. Myatt, E. A. Burt, R. W. Ghrist, E. A. Cornell, and C. E. Wieman, Phys. Rev. Lett. 78, 586 (1997).
- P. Maddaloni, M. Modugno, C. Fort, F. Minardi, and M. Inguscio, Phys. Rev. Lett. 85, 2413 (2000).
- K. M. Mertes, J. W. Merrill, R. Carretero-González, D. J. Frantzeskakis, P. G. Kevrekidis, and D. S. Hall, Phys. Rev. Lett. 99, 190402 (2007).
- S. B. Papp, J. M. Pino, and C. E. Wieman, Phys. Rev. Lett. 101, 040402 (2008)
- S. Sugawa, R. Yamazaki, S. Taie, and Y. Takahashi, Phys. Rev. A 84, 011610(R) (2011)
- G. Modugno, M. Modugno, F. Riboli, G. Roati, and M. Inguscio, Phys. Rev. Lett. 89, 190404 (2002)
- G. Thalhammer, G. Barontini, L. De Sarlo, J. Catani, F. Minardi, and M. Inguscio, Phys. Rev. Lett. 100, 210402 (2008)
- D. J. McCarron, H. W. Cho, D. L. Jenkin, M. P. Köppinger, and S. L. Cornish, Phys. Rev. A 84, 011603(R) (2011).
- A. D. Lercher, T. Takekoshi, M. Debatin, B. Schuster, R. Rameshan, F. Ferlaino, R. Grimm, H. -C. Nägerl, Eur. Phys. J. D 65, 3 (2011).
- B. Pasquiou, A. Bayerle, S. M. Tzanova, S. Stellmer, J. Szczepkowski, M. Parigger, R. Grimm, and F. Schreck, Phys. Rev. A 88, 023601 (2013)
- D. Xiong, X. Li, F. Wang, D. Wang, e-print arXiv:1305.7091
- T.-L. Ho and V. B. Shenoy, Phys. Rev. Lett. 77, 3276 (1996).
- H. Pu and N. P. Bigelow, Phys. Rev. Lett. 80, 1130 (1998); 80, 1134 (1998).
- P. Öhberg and S. Stenholm, Phys. Rev. A 57, 1272 (1998).
- M. Trippenbach, K. Góral, K. Rzążewski, B. Malomed and Y. B. Band, J. Phys. B: At. Mol. Opt. Phys. 33, 4017 (2000).
- R. W. Pattinson, T. P. Billam, S. A. Gardiner, D. J. McCarron, H. W. Cho, S. L. Cornish, N. G. Parker, and N. P. Proukakis, Phys. Rev. A 87, 013625 (2013).
- P. Öhberg, Phys. Rev. A 59, 634 (1999).
- N. P. Robins, W. Zhang, E. A. Ostrovskaya, and Y. S. Kivshar, Phys. Rev. A 64, 021601 (2001).
- C.-H. Zhang and H. A. Fertig, Phys. Rev. A 75, 013601 (2007).
- E. J. Mueller and T.-L. Ho, Phys. Rev. Lett. 88, 180403 (2002).
- N. G. Berloff Phys. Rev. Lett. 94, 120401 (2005).
- R. W. Pattinson, N. G. Parker, and N. P. Proukakis, J. Phys.: Conf. Ser. 497, 012029 (2014)
- I.-K. Liu, R. W. Pattinson, T. P. Billam, S. A. Gardiner, S. L. Cornish, T.-M. Huang, W.-W. Lin, S.-C. Gou, N. G. Parker and N. P. Proukakis, arXiv:1408.0891
- C. J. Pethick and H. Smith, Bose-Einstein Condensation in Dilute Gases (Cambridge University Press, Cambridge, 2002).
- G. Delannoy, S. G. Murdoch, V. Boyer, V. Josse, P. Bouyer, and A. Aspect, Phys. Rev. A 63, 051602(R) (2001); A. Mosk, S. Kraft, M. Mudrich, K. Singer, W. Wohlleben, R. Grimm, M. Weidemüller, Appl. Phys. B 73, 791 (2001); M. Anderlini, D. Ciampini, D. Cossart, E. Courtade, M. Cristiani, C. Sias, O. Morsch, and E. Arimondo, Phys. Rev. A 72, 033408 (2005).
- M. Lewenstein, J. I. Cirac, and P. Zoller, Phys. Rev. A 51, 4617 (1995).
- E. Timmermans and R. Côté, Phys. Rev. Lett. 80, 3419 (1998).
- T. Papenbrock, A. N. Salgueiro, and H. A. Weidenmüller, Phys. Rev. A 66, 025603 (2002).
- A. S. Bradley and P. B. Blakie, Phys. Rev. A 90, 023631 (2014).
- T. P. Billam, P. Mason, and S. A. Gardiner, Phys. Rev. A 87, 033628 (2013).
- P. Mason and S. A. Gardiner, Phys. Rev. A 89, 043617 (2014).
- T. Nikuni and J. E. Williams, J. Low Temp. Phys. 133, 323 (2003).
- Y. Endo and T. Nikuni, J. Low Temp. Phys. 163, 92 (2011).
- E. Zaremba, T. Nikuni, and A. Griffin, J. Low. Temp. Phys. 116, 277 (1999); A. Griffin, T. Nikuni, and E. Zaremba, Bose-Condensed Gases at Finite Temperatures (Cambridge University Press, Cambridge, 2009).
- T. Nikuni, E. Zaremba, and A. Griffin, Phys. Rev. Lett. 83, 10 (1999); T. Nikuni and A. Griffin, Phys. Rev. A 63, 033608 (2001).
- B. Jackson and E. Zaremba, Phys. Rev. A 66, 033606 (2002).
- B. Jackson and E. Zaremba, Phys. Rev. Lett. 87 100404 (2001); 89 150402 (2002); New J. Phys. 5 88 (2003).
- M. J. Bijlsma, E. Zaremba, and H. T. C. Stoof, Phys. Rev. A 62, 063609 (2000).
- B. Jackson, N. P. Proukakis, and C. F. Barenghi, Phys. Rev. A 75, 051601(R) (2007).
- J. Märkle, A. J. Allen, P. Federsel, B. Jetter, A. Günther, J. Fortágh, N. P. Proukakis, and T. E. Judd, Phys. Rev. A 90, 023614 (2014).
- N. P. Proukakis and K. Burnett, J. Res. Natl. Inst. Stand. Technol. 101, 457 (1996); N. P. Proukakis, K. Burnett, and H. T. C. Stoof, Phys. Rev. A 57, 1230 (1998).
- N. P. Proukakis, J. Phys. B: At. Mol. Opt. Phys. 34, 4737 (2001); N. P. Proukakis and B. Jackson, ibid. 41, 203002 (2008).
- R. Walser, J. Williams, J. Cooper, and M. Holland, Phys. Rev. A 59, 3878 (1999); J. Wachter, R. Walser, J. Cooper, and M. Holland, ibid. 64 053612 (2001).
- A. Griffin Phys. Rev. B 53, 9341 (1996).
- F. Dalfovo, S. Giorgini, L. P. Pitaevskii, and S. Stringari Rev. Mod. Phys. 71, 463 (1999).