Growth dynamics of a Bose-Einstein condensate in a dimple trap without cooling
We study the formation of a Bose-Einstein condensate in a cigar-shaped three-dimensional harmonic trap, induced by the controlled addition of an attractive “dimple” potential along the weak axis. In this manner we are able to induce condensation without cooling due to a localized increase in the phase space density. We perform a quantitative analysis of the thermodynamic transformation in both the sudden and adiabatic regimes for a range of dimple widths and depths. We find good agreement with equilibrium calculations based on self-consistent semiclassical Hartree-Fock theory describing the condensate and thermal cloud. We observe there is an optimal dimple depth that results in a maximum in the condensate fraction. We also study the non-equilibrium dynamics of condensate formation in the sudden turn-on regime, finding good agreement for the observed time dependence of the condensate fraction with calculations based on quantum kinetic theory.
Current address: ]Centre for Cold Matter, Blackett Laboratory, Imperial College London, London SW7 2AZ, United Kingdom Current address: ]ARC Centre of Excellence for Quantum-Atom Optics and Centre for Atom Optics and Ultrafast Spectroscopy, Swinburne University of Technology, Melbourne, VIC 3122, Australia Current address: ]Turku Centre for Quantum Physics, Department of Physics and Astronomy, University of Turku, FI-20014 Turku, Finland
The formation of a Bose-Einstein condensate (BEC) and the growth of long-range coherence from a gas of thermal atoms is a problem of interest in the field of ultracold atoms Stoof (2007). Before the first observations of Bose-Einstein condensation in a dilute gas there was some disagreement about the expected time scale for condensate formation Griffin et al. (1995). The first quantitative predictions were made by Gardiner et al. Gardiner et al. (1997), who derived a rate equation for the growth of a single condensate mode from a super-critical thermal vapor. This was soon followed by the first experimental measurements of condensate formation by Miesner et al. Miesner et al. (1998). Starting from just above the critical temperature for a BEC, this experiment implemented a sudden evaporative cooling ramp to remove the high-energy tail of a near-degenerate Bose gas. The ensuing rethermalization led to the formation of a Bose-Einstein condensate. These experiments were subsequently analyzed using improved formalisms by Gardiner and co-workers Gardiner et al. (1998); Lee and Gardiner (2000); Davis et al. (2000) and Bijlsma et al. Bijlsma et al. (2000). They found that, while their numerical calculations were qualitatively in agreement with experimental observations, quantitatively no agreement could be found, and this has remained unresolved.
In 2002 Köhl et al. Köhl et al. (2002) performed an experiment similar to that of Miesner et al. Miesner et al. (1998) but implemented continuous rather than sudden evaporation from near quantum degeneracy. For this experiment the data were generally in good agreement with the results of quantum kinetic calculations incorporating the details of the evaporation and the effects of three-body loss Davis et al. (2002). The same formalism was applied to later experiments in a quasi-condensate geometry by Hugbart et al. Hugbart et al. (2007), where the calculated shape of the condensate growth curves agreed well with experiment apart from an unexplained time delay. Other evaporative cooling experiments leading to BECs worth noting are those of Schvarchuck et al. Shvarchuck et al. (2002), who performed shock cooling in an elongated geometry and observed nonequilibrium dynamics in the resulting quasicondensate, and Ritter et al. Ritter et al. (2007), who measured the dynamics of the onset of long-range coherence in a three-dimensional condensate.
Bose-Einstein condensation without evaporative cooling was first induced by Stamper-Kurn et al. Stamper-Kurn et al. (1998), motivated by the earlier work of Pinkse et al. Pinkse et al. (1997). Stamper-Kurn et al. began with a near-degenerate Bose gas in a cigar-shaped harmonic trap and slowly applied an additional attractive “dimple” trap formed by a red-detuned optical dipole potential to adiabatically increase the phase-space density by a factor of 50. It was shown that this was reversible within the limits of heating caused by their dipole trap. Condensation was also induced by distillation without cooling in a double-well potential, demonstrated in an experiment by Shin et al. Shin et al. (2004). The lowering of a second well in this system caused the condensate in the first well to evaporate and re-form in the second at a higher temperature. Erhard et al. observed the formation of an BEC through spin collisions in an spinor condensate from initial partically condensed components in the states, and they modeled their experiment using rate equations Erhard et al. (2004). Recently an alternate approach to reversible BEC formation was demonstrated in an experiment by Catani et al. Catani et al. (2009), where entropy was exchanged between two atomic species, instead of between atoms inside and outside a dimple potential.
In this paper we revisit the method of Stamper-Kurn et al. Stamper-Kurn et al. (1998) to quantitatively study the dynamics and the thermodynamics of Bose-Einstein condensation, and we make comparisons of our experimental results with theoretical calculations. We induce condensate formation by the controlled application of a tightly focused laser sheet to a near-degenerate Bose gas in a cigar-shaped magnetic trap (illustrated in Fig. 1). The addition of the resulting one-dimensional dimple potential to the weakly confined dimension of the harmonic trap induces condensation by locally increasing density while the temperature remains almost constant, hence increasing the local phase-space density.
We divide our results into two sections. First, we have measured the final equilibrium state of the Bose gas following both the quasistatic (i.e., slow) and sudden turn-on of the dimple potential for a range of laser intensities beginning from a well-controlled initial nondegenerate state. The thermodynamics for an ideal gas with a delta function dimple have previously been studied in Ref. Uncu et al. (2007). Using a self-consistent mean-field model involving semiclassical Hartree-Fock theory for the thermal cloud and the Thomas-Fermi approximation for the condensate, we can predict the final condensate fraction for a given dimple depth for both the quasistatic and sudden turn-on. We perform a quantitative comparison of experiment and theory for a thermodynamic transformation through the BEC phase transition in an interacting Bose gas. Second, we observe the dynamics of condensate formation following sudden turn-on of the dimple potential and compare with a quantum kinetic model of condensate growth. This configuration allows a quantitative comparison with theory for condensate formation without evaporative cooling. We note that this scenario has been studied previously using stochastic classical fields in one dimension Stoof and Bijlsma (2001); Proukakis et al. (2006).
This paper is organized as follows: In Sec. II we summarize our experimental setup and procedure. In Sec. III we present our study of equilibrium thermodynamics by comparing the theoretical predictions with our experimental data and discussing the results. In Sec. IV, we present our study of condensate formation dynamics, first providing details of our quantum kinetic theory, and then comparing the theoretical predictions with our experimental data and discussing results. We finish with conclusions in Sec. V, and we provide additional theoretical details in the appendices.
Ii Experimental procedure
Our experiments are performed on a gas of ultracold Rb atoms in the state, confined in a magnetic harmonic potential provided by a Z-wire configuration on an atom chip, as previously described in Ref. Vale et al. (2004). Briefly, we collect atoms in a reflection MOT before transferring them to a magnetic trap at a distance of 200 m below the chip surface with axial and transverse frequencies of our cylindrically symmetric trapping potential Hz. Using rf evaporation we cool to near quantum degeneracy before decreasing the bias magnetic field to move the trap to 430 m below the chip surface, with a resulting decrease in trapping frequencies to Hz such that the resulting system is more three dimensional and less susceptible to phase fluctuations that exist in elongated condensates in the tighter trap Petrov et al. (2001). Further evaporative cooling results in a cloud of – atoms. We go as close as we can to the BEC transition while remaining above the transition temperature. This ensures the largest possible condensate fraction when the dimple potential is subsequently turned on.
Starting from these initial conditions, we next apply a red-detuned nm optical dipole potential to the system, intersecting the weak () direction of the magnetic trap at the center as illustrated in Fig. 1. This dipole beam is known by direct charge-coupled device (CCD) imaging to be Gaussian and diffraction limited. By observing the effect of shifting the focus back and forth along the dimension, we ensured the magnetic trap intersects the beam waist. The beam can be focused to two different half-widths in the dimension, thereby creating either a wide (32 m) or a narrow (11 m) dimple potential. This should be compared to the typical thermal cloud extent of m near the critical temperature. A cylindrical lens is used to expand the beam to widths of 350 m (wide) or 220 m (narrow) in the perpendicular dimension compared to the cloud width of 6 m, such that the intensity is approximately constant in these dimensions. Thus the cylindrically symmetric trapping potential can be approximated by
where is the half-width and the optical potential depth
is proportional to the peak laser intensity , and can reach a maximum depth of nK (wide) or 1610 nK (narrow). The other relevant parameters are the scattering cross sections , linewidths , the resonant frequencies of the and lines in Rb, and the laser frequency Grimm et al. (2000).
Our measurements are performed using absorption imaging after a time of flight of 20.3 ms, after turning off all trapping potentials. The temperature is determined by fitting a thermal cloud distribution to the wings of the resulting image, and the condensate fraction is determined from a two-component fit to the density. The condensate is sufficiently dense that condensate fractions of less that 1% can be distinguished from the thermal background.
Iii Thermodynamic transformations across the critical point
In our first set of experiments, we apply dimples of various depths and measure the resulting final equilibrium temperatures and condensate fractions. We do this for both wide and narrow dimples, and for both quasistatic and sudden turn-on. We then compare our measured values with the predictions of semiclassical Hartree-Fock theory, incorporating full mean-field interactions of both the condensate and thermal cloud, as well as accounting for the effects of three-body loss. While it is widely assumed that this is appropriate for the three-dimensional Bose gas, there have been relatively few comprehensive comparisons with experimental data aside from Ref. Gerbier et al. (2004).
iii.1 Theoretical procedure
We determine the initial entropy and total energy prior to dimple turn-on , given the experimentally measured initial temperature and atom number, using the semiclassical theory as outlined in Appendix A. To predict the final state, in the case of quasistatic turn-on we assume the system evolves isentropically and solve for the final temperature at which . In the case of sudden turn-on, we use the initial density to calculate the sudden change in energy imparted by the dimple,
where and are the densities of the condensate and thermal cloud, respectively. We then solve for the final temperature at which , which determines the final condensate fraction. We also estimate the effects of three-body loss in our calculations, as detailed in Appendix B.
iii.2 Comparison with experiment — wide dimple
Our experimental procedure is as follows. For the wide dimple measurements, we begin with an atomic cloud of atoms at a temperature of nK. This corresponds to a phase-space density at the center of the trap of , indicating that the cloud is very close to the BEC transition (). To turn on the dimple quasistatically, we ramp up the dimple potential linearly at a rate of 70s, and then hold it constant for a 300-ms equilibration time before turning off all potentials and imaging. To ensure that we are in the quasistatic regime, we have repeated this process for various ramp rates: For faster rates we observe a decrease in condensate fraction at large dimple depths due to nonadiabatic heating. To turn on the dimple suddenly, we ramp up the potential in less than 0.1 ms, and then hold it constant for a 1000-ms equilibration time before imaging.
The resulting temperatures and condensate fractions are plotted versus final dimple depth and compared with the predictions of semiclassical theory in Fig. 2. We find good agreement between theory and experiment for both the quasistatic and sudden turn-on, with the best fits obtained using initial conditions and nK. We also find the results of the quasistatic and sudden turn-on to be nearly indistinguishable from each other. This is because only a small fraction of the atoms are drawn from the harmonic trap into the wide dimple, which has a maximum attainable depth of order .
We compare the predictions of the semiclassical theory with our experimental data for both wide (Fig. 2) and narrow (Fig. 3) dimples. We plot temperature and condensate fraction versus dimple depth for both quasistatic and sudden turn-on. In all cases, the only fitting parameters used to generate the theoretical curves are the initial temperature and atom number, constrained to lie within their respective measurement uncertainties.
iii.3 Comparison with experiment — narrow dimple
Using the more tightly focused narrow dimple, we are able to attain a maximum depth much greater than and thereby observe differences between quasistatic and sudden turn-on. For the quasistatic turn-on we begin with an atomic cloud of atoms at a temperature of nK (). We ramp up the dimple potential linearly over a time of 1500 ms, and then hold the potential constant for a 300-ms equilibration time before imaging. For the sudden turn-on we begin with an atomic cloud of atoms at a temperature of nK () and follow the same procedure as for the wide dimple.
The resulting temperatures and condensate fractions are plotted versus final dimple depth and compared with theory in Fig. 3. In this case we observe a significant difference between the results of quasistatic and sudden turn-on, though this is partly due to the difference in initial conditions. We find good agreement between theory and experiment for the quasistatic turn-on, with the best fits obtained using initial conditions and nK.
However, we do not find good agreement for the sudden turn-on, with theory predicting a much smaller condensate fraction than was experimentally observed. The best fit shown is for and nK, though better fits can be obtained by using values of and that lie outside their respective measurement uncertainties.
A potential explanation for this discrepancy is that the turn-on is not sufficiently quick, and the density of the gas is not frozen during the turn-on. This would result in a smaller increase in total energy of the gas compared to the prediction of our model, leading to a smaller increase in temperature and hence a larger condensate fraction — as is observed experimentally. Indeed, in the limit of quasistatic turn-on from the same initial conditions (dotted lines in Fig. 3), the predicted condensate fraction is considerably larger than for sudden turn-on. Since the measured values lie between the limits of quasistatic and sudden turn-on, an intermediate turn-on time would likely provide a good fit to the data.
To model this would require a fully dynamical treatment that is beyond the limitations of our semiclassical theory and is numerically impractical within the quantum kinetic model we introduce later in Sec. IV.1. Furthermore, departures from the sudden case would only be expected at turn-on times similar to the time scale for rethermalization. Since the mean free time between collisions at initial conditions is approximately 6 ms, whereas the time for sudden turn-on is less than 0.1 ms, this explanation seems unlikely. We are therefore unfortunately forced to leave this discrepancy unresolved.
An interesting feature of the narrow dimple data is that there exists an optimal dimple depth ( ) at which a maximum condensate fraction is obtained. As the dimple depth is increased beyond this value, the condensate fraction gradually decreases back toward zero. We can understand this feature in the context of the condensate formation process as follows. At shallow depths, the dimple potential acts merely as a perturbation to the broader harmonic trap. To a first approximation the chemical potential and temperature are unchanged, whereas the energy of the translational ground state is decreased relative to the bottom of the harmonic trap (see Fig. 4). When the ground-state energy approaches the chemical potential, a condensate forms, as observed for the wide dimple and at shallow depths of the narrow dimple. However, at larger depths a significant fraction of the thermal cloud is drawn into the narrow dimple, causing non-negligible changes in and . At sufficiently large depths, the entire thermal cloud falls into the dimple potential (which is approximately harmonic near its center) and the condensate evaporates, in agreement with the well-known result that compression of a gas cannot alter the phase-space density Houbiers et al. (1997).
Iv Condensate formation dynamics
In our second set of experiments, we suddenly turn on the dimple and allow the system to evolve for various times before measuring the condensate fraction. We do this for both wide and narrow dimples, in each case for two different fixed dimple depths: one shallow and the other deep. We then compare our measured values with the predictions of quantum kinetic theory. In particular, we use the ergodic quantum Boltzmann equation Luiten et al. (1996) with the additional inclusion of the effects of the condensate mean-field Davis et al. (2000); Bijlsma et al. (2000) and three-body loss. While similar comparisons have been made previously Bijlsma et al. (2000); Davis et al. (2000); Köhl et al. (2002); Davis et al. (2002); Hugbart et al. (2007) the condensation transitions in these experiments were induced by evaporative cooling. Here condensation is introduced without loss by modifying the density of states of the trap.
iv.1 Quantum kinetic theory
The initial equilibrium state of the Bose gas above becomes nonequilibrium following the sudden turn-on of a dimple potential. We begin by calculating this initial nonequilibrium state, and then evolve it in time to final equilibrium using an ergodic quantum Boltzmann equation (EQBE), as previously described in Luiten et al. (1996); Davis et al. (2000); Bijlsma et al. (2000).
The ergodic assumption is that all semiclassical states for a given single-particle energy have the same mean occupation and thus the nonequilibrium phase-space distribution function depends only on time and energy: . Hence the position and momentum dependence in the full quantum Boltzmann equation is projected out, yielding a partial differential equation in energy and time only. The ergodic approximation is necessary in order to reduce the dimensionality of the phase space to render the quantum Boltzmann equation computationally tractable. The EQBE then governs the evolution of the energy density
where is the energy distribution function and
is the density of states. The semiclassical excitation energy is as defined in Eq. (8) but with the mean field of the thermal cloud neglected [i.e., setting to zero]. This eliminates the need to determine the densities self-consistently, and our results suggest that this approximation is reasonable. The condensate mean field is calculated via the Thomas-Fermi approximation, as per Eq. (10), but again neglecting the mean field of the thermal cloud. We provide a more detailed description of the EQBE in Appendix C and describe how to incorporate the effects of three-body loss in Appendix B.
A typical atom will travel a distance of less than 1 during the sudden ramp-up of the dimple potential. As this is more than an order of magnitude smaller than the width of the narrow dimple, we approximate the sudden turn-on as instantaneous. Under this assumption, the phase-space distribution immediately after the sudden turn-on must be the same as before the turn-on: , where we henceforth use primed (unprimed) variables to denote quantities immediately after (before) dimple turn-on at . Although the phase-space distribution function is unchanged, the post-dimple semiclassical excitation energy differs from the pre-dimple expression via the inclusion of the dimple potential:
Clearly, the change in alters the density of states in Eq. (5), and therefore it alters the energy density, which we calculate immediately after the dimple turn-on as
In doing this we are ergodically projecting the pre-dimple phase-space distribution using the post-dimple density of states. It should be noted that the physical initial phase space will actually be non-ergodic. However, previous Monte Carlo calculations of the Boltzmann equation have shown that ergodicity is restored relatively quickly Wu and Foot (1996), and hence we expect this should be a reasonable approximation.
The pre- and post-dimple energy densities are shown in Fig. 5, for parameters typical of the narrow dimple: nK, , and nK. The post-dimple energy density closely resembles the pre-dimple energy density, except that the peak in post-dimple energy density is slightly smaller because a small fraction of atoms lying within the dimple now have energies below the bottom of the harmonic potential (). Also plotted in Fig. 5 is the energy density at final equilibrium, calculated via the EQBE. As the system evolves toward equilibrium, successively more atoms are drawn into the dimple, occupying the energies in the range , and ultimately resulting in the formation of a condensate.
iv.2 Comparison with experiment
We compare the predictions of the quantum kinetic calculations just described with our experimental data for both wide (Fig. 6) and narrow (Fig. 7) dimples. In both cases, we plot condensate fraction versus time for two different dimple depths. In all cases, the only fitting parameters used to generate the theoretical curves are the initial temperature and atom number, constrained to lie within their respective measurement uncertainties.
For the wide dimple measurements, we begin with an atomic cloud of atoms at a temperature of nK (). We ramp up the potential suddenly, and then hold it constant for various equilibration times up to 1000 ms before imaging. We use dimples of two different depths: and nK. The resulting condensate fractions are plotted versus time and compared with the predictions of quantum kinetic theory in Fig. 6. We find good agreement between theory and experiment, with the best fits obtained using initial conditions (, nK) and (, nK) for dimple depths and nK, respectively.
For the narrow dimple measurements, we begin with an atomic cloud of atoms at a temperature of nK (), and we follow the same experimental procedure described earlier for a final dimple depth of nK. The results are plotted in Fig. 7, with the best fit obtained using initial conditions and nK. We find good agreement between theory and experiment, particularly at final equilibrium. In contrast to the wide dimple case, we observe a delay of about 50 ms in the onset of condensate formation, which agrees with the EQBE calculation.
We have studied the thermodynamics of the transition across the BEC critical point and the formation dynamics of Bose-Einstein condensation resulting from the application of attractive Gaussian dimple potentials of various widths and depths to a Bose gas in a cigar-shaped harmonic trap. We have measured the equilibrium temperature and condensate fraction over a range of dimple depths, for both quasistatic and sudden turn-on of the dimple, and compared our results with the predictions of semiclassical Hartree-Fock theory. For narrow dimples, we found that there exists an optimal dimple depth at which a maximum condensate fraction is attained. Beyond this depth, the dimple acts merely as a tighter harmonic trap, and it therefore does not increase the phase-space density. We found good agreement between theory and experiment, except in the case of sudden turn-on of deep, narrow dimples. We also measured the (nonequilibrium) condensate fraction over a range of times after sudden turn-on of the dimple and compared our results with the predictions of quantum kinetic theory. We observed a short delay in the onset of condensate formation in the case of narrow dimples, but not in the case of wide dimples. In both cases we found good agreement between theory and experiment. This provides further validation to the quantum kinetic model of condensate formation, for which previous comparisons were based on sudden evaporative cooling.
Acknowledgements.This research was supported under the Australian Research Council’s Discovery Projects funding scheme (Project Nos. DP0343094 and DP0985142) and the ARC Centre of Excellence for Quantum-Atom Optics (CE0348178). MCG acknowledges financial support from NSERC, Endeavour IPRS, and the University of Queensland.
Appendix A semiclassical theory
We make use of the Thomas-Fermi approximation for the condensate density and the semiclassical Hartree-Fock approximation for the thermal cloud, as outlined in Goldman et al. (1981); Huse and Siggia (1982); Oliva (1989); Bagnato et al. (1987); Shi and Zheng (1997a, b); Giorgini et al. (1997) and compared with experiment by Gerbier et al. Gerbier et al. (2004).
At equilibrium the thermal cloud is well described by the Bose-Einstein distribution, , where is the temperature and is the chemical potential. In the semiclassical Hartree-Fock approximation, the excitation energy (the energy required to add an atom with momentum at position ) is given by the expression
where is the atomic mass and is the interaction strength, proportional to the -wave scattering length . The mean-field density combined with the external potential constitute the effective potential experienced by atoms in the thermal cloud. The density of the thermal cloud is calculated as
while the density of the condensate is calculated using the Thomas-Fermi approximation,
Thus, the densities of both the thermal cloud and condensate must be solved self-consistently for a given and to give the experimentally measured total atom number .
From here other thermodynamics quantities such as the condensate fraction , total entropy, and total energy can be determined. The last two are relevant to quasistatic and sudden turn-on of the dimple potential, respectively. By writing the total entropy is given by (cf. p. 15 of Pitaevskii and Stringari (2003))
while the total energy of the system is
Appendix B Three-body loss
Three-body loss occurring between condensed and non-condensed atoms can be calculated from the three-body correlation function as
where Burt et al. (1997). Using a broken symmetry approach, we write the Bose field operator as a mean field plus fluctuations,
and substitute this into the previous expression. Identifying the condensate density as yields
The noncondensate density is given by . Using Wick’s theorem on the higher order operator moments of the fluctuations gives
The loss rates for the condensate and thermal cloud atoms can then be written separately by noting that, for example, the second term in this above expression represents a three-body process in which two condensate atoms and one thermal cloud atom are lost. The coefficients can be divided accordingly to give Söding et al. (1999)
which ensures . To incorporate these loss rates in our semiclassical theory (Appendix A), we treat these expressions for and as coupled ordinary differential equations, and we solve them in time at each spatial gridpoint via Euler’s method. This is not entirely straightforward, because the semiclassical theory is strictly static: The final equilibrium is determined from the initial equilibrium in a single leap, without any stepwise time evolution. In the case of sudden dimple turn-on, we begin with the final equilibrium densities, and we calculate three-body loss over the 1-s equilibration time. In essence, we are assuming that the three-body loss rates on the final equilibrium densities are approximately equal to the average loss rates during the 1-s evolution from initial to final equilibrium. In the case of the quasistatic dimple turn-on, we can be more precise: Because the system is never out of equilibrium as it evolves from initial to final equilibrium, we can break up the process into arbitrarily many steps. For each step, we evolve the system isentropically as the dimple depth is incrementally increased, and we calculate three-body loss over the corresponding time interval.
The inclusion of three-body loss has a significant effect on the predictions of both semiclassical Hartree-Fock and quantum kinetic theory, as we show in Fig. 8. This is particularly evident in the case of adiabatic turn-on, where the inclusion of three-body loss not only drastically reduces the maximum condensate fraction but also reduces the optimal dimple depth at which the maximum fraction is attained. All of the curves shown correspond to the narrow dimple, with initial conditions and nK.
Appendix C Ergodic quantum Boltzmann equation
In our implementation of the EQBE, we use an energy grid with uniform spacing , with the lowest energy bin corresponding to the condensate. Because the condensate energy level (equal to the bottom of the effective potential) changes with time, we redefine our energy grid to be , where is the condensate chemical potential in the Thomas-Fermi approximation. With this transformation the EQBE is Bijlsma et al. (2000)
where is the weighted density of states, defined as
which contains the time derivative of the effective potential and hence depends on the rate of condensate growth . The contribution of binary elastic collisions is given by the term
where , , and . The first (second) term within the square brackets represents the forward (backward) collisions . Factors of the form are due to Bose enhancement and vanish in the classical limit . For the case of collisions involving a condensate atom, , we make the replacement . Detailed derivations of the EQBE can be found in Luiten et al. (1996); Davis et al. (2000); Bijlsma et al. (2000).
We numerically evolve the EQBE in time using an explicit fourth-order Runge-Kutta method. At each time step, we calculate by summing over all combinations of energy bins satisfying the delta function in Eq. (C) and updating for all four bins involved in each combination. We then add the contributions from three-body loss [see Eq. (29) and (30)], and use the resulting value of to calculate in the first term on the right-hand side of the EQBE.
To incorporate three-body loss in our quantum kinetic calculations, the loss rate for the thermal cloud must be modified to include energy dependence:
where the energy-position density of the thermal cloud is calculated from the energy density as
and is the position-dependent density of states, obtained by omitting the spatial integral from Eq. (5). At each time step in the numerical evolution of the EQBE, we determine the condensate density from the Thomas-Fermi approximation [Eq. (10), with the thermal cloud mean-field neglected] and determine the thermal cloud density by integrating this expression for over energy. We then calculate the loss rate for the condensate, Eq. (22), and the energy-dependent loss rate for the thermal cloud, Eq. (27). Lastly, we integrate out the spatial dependence to get
which are simply added to the EQBE calculation of in Eq. (24).
- Stoof (2007) H. T. C. Stoof, Nature 447, 390 (2007).
- Griffin et al. (1995) A. Griffin, D. W. Snoke, and S. Stringari, Bose-Einstein Condensation, 1st ed. (Cambridge University Press, Cambridge, UK, 1995).
- Gardiner et al. (1997) C. W. Gardiner, P. Zoller, R. J. Ballagh, and M. J. Davis, Phys. Rev. Lett. 79, 1793 (1997).
- Miesner et al. (1998) H.-J. Miesner, D. M. Stamper-Kurn, M. R. Andrews, D. S. Durfee, S. Inouye, and W. Ketterle, Science 279, 1005 (1998).
- Gardiner et al. (1998) C. W. Gardiner, M. D. Lee, R. J. Ballagh, M. J. Davis, and P. Zoller, Phys. Rev. Lett. 81, 5266 (1998).
- Lee and Gardiner (2000) M. D. Lee and C. W. Gardiner, Phys. Rev. A 62, 033606 (2000).
- Davis et al. (2000) M. J. Davis, C. W. Gardiner, and R. J. Ballagh, Phys. Rev. A 62, 063608 (2000).
- Bijlsma et al. (2000) M. J. Bijlsma, E. Zaremba, and H. T. C. Stoof, Phys. Rev. A 62, 063609 (2000).
- Köhl et al. (2002) M. Köhl, M. J. Davis, C. W. Gardiner, T. W. Hänsch, and T. Esslinger, Phys. Rev. Lett. 88, 080402 (2002).
- Davis et al. (2002) M. J. Davis, S. A. Morgan, and K. Burnett, Phys. Rev. A. 66, 053618 (2002).
- Hugbart et al. (2007) M. Hugbart, J. A. Retter, A. F. Varon, P. Bouyer, A. Aspect, and M. J. Davis, Phys. Rev. A 75, 011602 (2007).
- Shvarchuck et al. (2002) I. Shvarchuck, C. Buggle, D. S. Petrov, K. Dieckmann, M. Zielonkowski, M. Kemmann, T. G. Tiecke, W. von Klitzing, G. V. Shlyapnikov, and J. T. M. Walraven, Phys. Rev. Lett. 89, 270404 (2002).
- Ritter et al. (2007) S. Ritter, A. Öttl, T. Donner, T. Bourdel, M. Köhl, and T. Esslinger, Phys. Rev. Lett. 98, 090402 (2007).
- Stamper-Kurn et al. (1998) D. M. Stamper-Kurn, H.-J. Miesner, A. P. Chikkatur, S. Inouye, J. Stenger, and W. Ketterle, Phys. Rev. Lett. 81, 2194 (1998).
- Pinkse et al. (1997) P. W. H. Pinkse, A. Mosk, M. Weidemüller, M. W. Reynolds, T. W. Hijmans, and J. T. M. Walraven, Phys. Rev. Lett. 78, 990 (1997).
- Shin et al. (2004) Y. Shin, M. Saba, A. Schirotzek, T. A. Pasquini, A. E. Leanhardt, D. E. Pritchard, and W. Ketterle, Phys. Rev. Lett. 92, 150401 (2004).
- Erhard et al. (2004) M. Erhard, H. Schmaljohann, J. Kronjäger, K. Bongs, and K. Sengstock, Phys. Rev. A 70, 031602(R) (2004).
- Catani et al. (2009) J. Catani, G. Barontini, G. Lamporesi, F. Rabatti, G. Thalhammer, F. Minardi, S. Stringari, and M. Inguscio, Phys. Rev. Lett. 103, 140401 (2009).
- Uncu et al. (2007) H. Uncu, D. Tarhan, E. Demiralp, and O. E. Müstecaplıoğlu, Phys. Rev. A 76, 013618 (2007).
- Stoof and Bijlsma (2001) H. T. C. Stoof and M. J. Bijlsma, J. Low Temp. Phys. 124, 431 (2001).
- Proukakis et al. (2006) N. P. Proukakis, J. Schmiedmayer, and H. T. C. Stoof, Phys. Rev. A 73, 053603 (2006).
- Vale et al. (2004) C. J. Vale, B. Upcroft, M. J. Davis, N. R. Heckenberg, and H. Rubinsztein-Dunlop, J. Phys. B: At. Mol. Opt. Phys. 37, 2959 (2004).
- Petrov et al. (2001) D. S. Petrov, G. V. Shlyapnikov, and J. T. M. Walraven, Phys. Rev. Lett. 87, 050404 (2001).
- Grimm et al. (2000) R. Grimm, M. Weidemüller, and Y. B. Ovchinnikov, Adv. At. Mol. Opt. Phy. 45, 95 (2000).
- Gerbier et al. (2004) F. Gerbier, J. H. Thywissen, S. Richard, M. Hugbart, P. Bouyer, and A. Aspect, Phys. Rev. A 70, 013607 (2004).
- Houbiers et al. (1997) M. Houbiers, H. T. C. Stoof, and E. A. Cornell, Phys. Rev. A 56, 2041 (1997).
- Luiten et al. (1996) O. J. Luiten, M. W. Reynolds, and J. T. M. Walraven, Phys. Rev. A 53, 381 (1996).
- Wu and Foot (1996) H. Wu and C. J. Foot, J. Phys. B 29, L321 (1996).
- Goldman et al. (1981) V. V. Goldman, I. F. Silvera, and A. J. Leggett, Phys. Rev. B 24, 2870 (1981).
- Huse and Siggia (1982) D. A. Huse and E. D. Siggia, J. Low Temp. Phys. 46, 137 (1982).
- Oliva (1989) J. Oliva, Phys. Rev. B 39, 4197 (1989).
- Bagnato et al. (1987) V. Bagnato, D. E. Pritchard, and D. Kleppner, Phys. Rev. A 35, 4354 (1987).
- Shi and Zheng (1997a) H. Shi and W.-M. Zheng, Phys. Rev. A 56, 2984 (1997a).
- Shi and Zheng (1997b) H. Shi and W.-M. Zheng, Phys. Rev. A 56, 1046 (1997b).
- Giorgini et al. (1997) S. Giorgini, L. P. Pitaevskii, and S. Stringari, J. Low Temp. Phys. 109, 309 (1997).
- Pitaevskii and Stringari (2003) L. Pitaevskii and S. Stringari, Bose-Einstein Condensation (Oxford Science, Oxford, UK, 2003).
- Burt et al. (1997) E. A. Burt, R. W. Ghrist, C. J. Myatt, M. J. Holland, E. A. Cornell, and C. E. Wieman, Phys. Rev. Lett. 79, 337 (1997).
- Söding et al. (1999) J. Söding, D. Guéry-Odelin, P. Desbiolles, F. Chevy, H. Inamori, and J. Dalibard, Appl. Phys. B 69, 257 (1999).