Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:cond-mat/0602062 · Physical Review A

Three-body recombination of ultracold Bose gases using the truncated Wigner method

A. A. Norrie, R. J. Ballagh, Crispin W. Gardiner, Ashton S. Bradley

Published Apr 27, 2006• 46 citations
arXivHTMLPDF

Abstract

We apply the truncated Wigner method to the process of three-body recombination in ultracold Bose gases. We find that within the validity regime of the Wigner truncation for two-body scattering, three-body recombination can be treated using a set of coupled stochastic differential equations that include diffusion terms, and can be simulated using known numerical methods. As an example we investigate the behavior of a simple homogeneous Bose gas, finding a very slight increase of the loss rate compared to that obtained by using the standard method.

DOI
10.1103/physreva.73.043618

Related papers

Variational truncated Wigner approximation for weakly interacting Bose fields: Dynamics of coupled condensates

Christopher D. Mink, Axel Pelster, Jens Benary +2 more

The truncated Wigner approximation is an established approach that describes the dynamics of weakly interacting Bose gases beyond the mean-field level. Although it allows a quantum field to be expressed by a stochastic c-number field, the simulation of the time evolution is still very demanding for most applications. Here, we develop a numerically inexpensive scheme by approximating the c-number field with a variational ansatz. The dynamics of the ansatz function is described by a tractable set of coupled ordinary stochastic differential equations for the respective variational parameters. We investigate the non-equilibrium dynamics of a three-dimensional Bose gas in a one-dimensional optical lattice with a transverse isotropic harmonic confinement. The accuracy and computational inexpensiveness of our method are demonstrated by comparing its predictions to experimental data.

arXiv:2106.05354•SciPost Physics•Feb 2, 2022•13 citations

Truncated Wigner dynamics and conservation laws

P. D. Drummond, Bogdan Opanchuk

Ultracold Bose gases can be used to experimentally test many-body theory predictions. Here we point out that both exact conservation laws and dynamical invariants exist in the topical case of the one-dimensional Bose gas, and these provide an important validation of methods. We show that the first four quantum conservation laws are exactly conserved in the approximate truncated Wigner approach to many-body quantum dynamics. Center-of-mass position variance is also exactly calculable. This is nearly exact in the truncated Wigner approximation, apart from small terms that vanish as ${N}^{\ensuremath{-}3/2}$ as $N\ensuremath{\rightarrow}\ensuremath{\infty}$ with fixed momentum cutoff. Examples of this are calculated in experimentally relevant, mesoscopic cases.

Physical Review A•Oct 18, 2017•34 citations

The truncated Wigner method for Bose-condensed gases: limits of validity and applications

Alice Sinatra, Carlos Lobo, Yvan Castin

We study the truncated Wigner method applied to a weakly interacting spinless Bose-condensed gas which is perturbed away from thermal equilibrium by a time-dependent external potential. The principle of the method is to generate an ensemble of classical fields ψ( r ) which samples the Wigner quasi-distribution function of the initial thermal equilibrium density operator of the gas, and then to evolve each classical field with the Gross–Pitaevskii equation. In the first part of the paper we improve the sampling technique over our previous work (Sinatra et al 2000 J. Mod. Opt. 47 2629–44) and we test its accuracy against the exactly solvable model of the ideal Bose gas. In the second part of the paper we investigate the conditions of validity of the truncated Wigner method. For short evolution times it is known that the time-dependent Bogoliubov approximation is valid for almost pure condensates. The requirement that the truncated Wigner method reproduces the Bogoliubov prediction leads to the constraint that the number of field modes in the Wigner simulation must be smaller than the number of particles in the gas. For longer evolution times the nonlinear dynamics of the noncondensed modes of the field plays an important role. To demonstrate this we analyse the case of a three-dimensional spatially homogeneous Bose-condensed gas and we test the ability of the truncated Wigner method to correctly reproduce the Beliaev–Landau damping of an excitation of the condensate. We have identified the mechanism which limits the validity of the truncated Wigner method: the initial ensemble of classical fields, driven by the time-dependent Gross–Pitaevskii equation, thermalizes to a classical field distribution at a temperature T class which is larger than the initial temperature T of the quantum gas. When T class significantly exceeds T a spurious damping is observed in the Wigner simulation. This leads to the second validity condition for the truncated Wigner method, T class − T ≪ T , which requires that the maximum energy max of the Bogoliubov modes in the simulation does not exceed a few k B T .

arXiv:cond-mat/0201217•Journal of Physics B Atomic Molecular and Optical Physics•Aug 21, 2002•352 citations

Ultracold three-body recombination of fermionic atoms

Brett D. Esry, Hiroya Suno, Chris H. Greene

Three-body recombination can play an important role in experiments on ultracold atomic gases because it can produce significant atom loss. Such losses have been seen in Bose-Einstein condensates, but have been thought to play little role for degenerate Fermi gases. We show that under some circumstances, three-body recombination can be just as large in fermi systems as in bose systems. In particular, when there is a two-body resonance near threshold — as is often the case experimentally — the threshold suppression due to the Pauli principle no longer applies. 1

International Conference on Atomic Physics (ICAP)•Apr 1, 2003•2 citations

Can Three-Body Recombination Purify a Quantum Gas?

Lena Hannah Dogra, Jake Glidden, Timon Hilker +4 more

Three-body recombination in quantum gases is traditionally associated with heating, but it was recently found that it can also cool the gas. We show that in a partially condensed three-dimensional homogeneous Bose gas three-body loss could even purify the sample, that is, reduce the entropy per particle and increase the condensed fraction η. We predict that the evolution of η under continuous three-body loss can, depending on small changes in the initial conditions, exhibit two qualitatively different behaviors-if it is initially above a certain critical value, η increases further, whereas clouds with lower initial η evolve towards a thermal gas. These dynamical effects should be observable under realistic experimental conditions.

arXiv:1905.09811•Physical Review Letters•Jul 12, 2019•21 citations

Five-body recombination of identical bosons

Michael D. Higgins, Chris H. Greene

This work treats resonant collisions between five identical ultracold bosons in the framework of the adiabatic hyperspherical representation. The five-body recombination rate coefficient is quantified using a semiclassical description in conjunction with an analysis of the lowest five-body hyperspherical adiabatic potential curves in a scattering length regime with no universal weakly bound tetramers, trimers, or dimers. A comparison is made between these results and the only existing experimental measurement of five-body loss in an ultracold gas of bosonic cesium atoms and with the lone theoretical estimation of the loss rate. The recombination rate for the process B + B + B + B + B → B 4 + B is also computed in a different regime of scattering lengths where there is one universal bound tetramer by implementing a few-channel quantum scattering calculation based on five-body hyperspherical potential curves and nonadiabatic couplings. Our calculations predict regions where five-body recombination can cause decay of the atom cloud in an ultracold gas that is even faster than 3-body and 4-body recombination, which can ideally be tested by using the current generation of box traps having nearly uniform density.

Proceedings of the National Academy of Sciences•Apr 25, 2025•2 citations
We apply the truncated Wigner method to the process of three-body recombination in ultracold Bose gases. We find that within the validity regime of the Wigner truncation for two-body scattering, three-body recombination can be treated using a set of coupled stochastic differential equations that include diffusion terms, and can be simulated using known numerical methods. As an example we investigate the behavior of a simple homogeneous Bose gas, finding a very slight increase of the loss rate compared to that obtained by using the standard method.
The truncated Wigner approximation is an established approach that describes the dynamics of weakly interacting Bose gases beyond the mean-field level. Although it allows a quantum field to be expressed by a stochastic c-number field, the simulation of the time evolution is still very demanding for most applications. Here, we develop a numerically inexpensive scheme by approximating the c-number field with a variational ansatz. The dynamics of the ansatz function is described by a tractable set of coupled ordinary stochastic differential equations for the respective variational parameters. We investigate the non-equilibrium dynamics of a three-dimensional Bose gas in a one-dimensional optical lattice with a transverse isotropic harmonic confinement. The accuracy and computational inexpensiveness of our method are demonstrated by comparing its predictions to experimental data.
Ultracold Bose gases can be used to experimentally test many-body theory predictions. Here we point out that both exact conservation laws and dynamical invariants exist in the topical case of the one-dimensional Bose gas, and these provide an important validation of methods. We show that the first four quantum conservation laws are exactly conserved in the approximate truncated Wigner approach to many-body quantum dynamics. Center-of-mass position variance is also exactly calculable. This is nearly exact in the truncated Wigner approximation, apart from small terms that vanish as ${N}^{\ensuremath{-}3/2}$ as $N\ensuremath{\rightarrow}\ensuremath{\infty}$ with fixed momentum cutoff. Examples of this are calculated in experimentally relevant, mesoscopic cases.
We study the truncated Wigner method applied to a weakly interacting spinless Bose-condensed gas which is perturbed away from thermal equilibrium by a time-dependent external potential. The principle of the method is to generate an ensemble of classical fields ψ( r ) which samples the Wigner quasi-distribution function of the initial thermal equilibrium density operator of the gas, and then to evolve each classical field with the Gross–Pitaevskii equation. In the first part of the paper we improve the sampling technique over our previous work (Sinatra et al 2000 J. Mod. Opt. 47 2629–44) and we test its accuracy against the exactly solvable model of the ideal Bose gas. In the second part of the paper we investigate the conditions of validity of the truncated Wigner method. For short evolution times it is known that the time-dependent Bogoliubov approximation is valid for almost pure condensates. The requirement that the truncated Wigner method reproduces the Bogoliubov prediction leads to the constraint that the number of field modes in the Wigner simulation must be smaller than the number of particles in the gas. For longer evolution times the nonlinear dynamics of the noncondensed modes of the field plays an important role. To demonstrate this we analyse the case of a three-dimensional spatially homogeneous Bose-condensed gas and we test the ability of the truncated Wigner method to correctly reproduce the Beliaev–Landau damping of an excitation of the condensate. We have identified the mechanism which limits the validity of the truncated Wigner method: the initial ensemble of classical fields, driven by the time-dependent Gross–Pitaevskii equation, thermalizes to a classical field distribution at a temperature T class which is larger than the initial temperature T of the quantum gas. When T class significantly exceeds T a spurious damping is observed in the Wigner simulation. This leads to the second validity condition for the truncated Wigner method, T class − T ≪ T , which requires that the maximum energy max of the Bogoliubov modes in the simulation does not exceed a few k B T .
Three-body recombination can play an important role in experiments on ultracold atomic gases because it can produce significant atom loss. Such losses have been seen in Bose-Einstein condensates, but have been thought to play little role for degenerate Fermi gases. We show that under some circumstances, three-body recombination can be just as large in fermi systems as in bose systems. In particular, when there is a two-body resonance near threshold — as is often the case experimentally — the threshold suppression due to the Pauli principle no longer applies. 1
Three-body recombination in quantum gases is traditionally associated with heating, but it was recently found that it can also cool the gas. We show that in a partially condensed three-dimensional homogeneous Bose gas three-body loss could even purify the sample, that is, reduce the entropy per particle and increase the condensed fraction η. We predict that the evolution of η under continuous three-body loss can, depending on small changes in the initial conditions, exhibit two qualitatively different behaviors-if it is initially above a certain critical value, η increases further, whereas clouds with lower initial η evolve towards a thermal gas. These dynamical effects should be observable under realistic experimental conditions.
This work treats resonant collisions between five identical ultracold bosons in the framework of the adiabatic hyperspherical representation. The five-body recombination rate coefficient is quantified using a semiclassical description in conjunction with an analysis of the lowest five-body hyperspherical adiabatic potential curves in a scattering length regime with no universal weakly bound tetramers, trimers, or dimers. A comparison is made between these results and the only existing experimental measurement of five-body loss in an ultracold gas of bosonic cesium atoms and with the lone theoretical estimation of the loss rate. The recombination rate for the process B + B + B + B + B → B 4 + B is also computed in a different regime of scattering lengths where there is one universal bound tetramer by implementing a few-channel quantum scattering calculation based on five-body hyperspherical potential curves and nonadiabatic couplings. Our calculations predict regions where five-body recombination can cause decay of the atom cloud in an ultracold gas that is even faster than 3-body and 4-body recombination, which can ideally be tested by using the current generation of box traps having nearly uniform density.