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:2410.12204 · Physical Review Letters

Fate of Thermalization of Ultracold Fermions with Two-Body Dissipation

Xin-Yuan Gao, Yangqian Yan

Published Apr 17, 2025• 2 citations
arXivHTMLPDF

Abstract

Two-body inelastic collisions arising from chemical reactions are prevalent in ultracold fermionic and bosonic molecular gases. Although recent advancements have achieved quantum degeneracy in these systems, loss dynamics are typically modeled phenomenologically using rate equations that often assume thermalization during chemical reactions. In this study, we employ the inelastic quantum Boltzmann equation to analyze particle loss, temperature evolution, and momentum distributions in single-component Fermi gases from first principles. Our results demonstrate that the conventional particle-number rate equation accurately describes the dynamics in trapped systems but fails to capture the behavior in homogeneous systems. Notably, under pure p-wave inelastic collisions and zero elastic collisions, we find that systems prepared near or above quantum degeneracy remain in a thermal state, whereas systems initialized deep within degeneracy exhibit nonequilibrium dynamics. Our theoretical predictions align well with recent experimental observations in trapped systems, and our claim can be further verified in atomic systems with induced two-body loss in box potentials.

DOI
10.1103/physrevlett.134.153402

Related papers

Inelastic quantum Boltzmann equation of single-component Fermi gases

Xin-Yuan Gao, Yangqian Yan

This work is a companion paper to X.-Y. Gao et al., Phys. Rev. Lett. 134, 153402 (2025), where we discuss the nonequilibrium two-body loss dynamics of a single-component ultracold Fermi gas and its possible thermalization. This paper provides detailed information on the derivation and analysis of the inelastic quantum Boltzmann equation (IQBE) used to describe the system. We demonstrate that the Mellin transform is a powerful tool for solving and approximating the IQBE for free-space systems. In this case, the particle-number dynamics are beyond the description of the widely used phenomenological two-body equation. For harmonically trapped systems, we propose a fast-flowing approximation to simplify the numerical evaluation of the IQBE. We verify the approximation in an analogous quasi-one-dimensional system and apply it to three-dimensional calculations, obtaining satisfactory agreement with recent experimental results. Furthermore, we compare the nonequilibrium results with those obtained using a thermal ansatz in both situations, providing a systematic understanding of the antievaporation phenomena observed in such systems.

Physical Review A•Apr 17, 2025•1 citation

Effects of collisions against thermal impurities in the dynamics of a trapped fermion gas

P. Capuzzi, Patrizia Vignolo, Federico Toschi +2 more

We present a theoretical study of the dynamical behavior of a gas made of ultracold fermionic atoms, which during their motions can collide with a much smaller number of thermal bosonic impurities. The atoms are confined inside harmonic traps and the interactions between the two species are treated as due to $s$-wave scattering with a negative scattering length modeling the $^{40}\mathrm{K}\penalty1000-\hskip0pt^{87}\mathrm{Rb}$ fermion-boson system. We set the fermions into motion by giving a small shift to their trap center and examine two alternative types of initial conditions, referring to (i) a close-to-equilibrium situation in which the two species are at the same temperature (well below the Fermi temperature and well above the Bose-Einstein condensation temperature); and (ii) a far-from-equilibrium case in which the impurities are given a Boltzmann distribution of momenta while the fermions are at very low temperatures. The dynamics of the gas is evaluated by the numerical solution of the Vlasov-Landau equations for the one-body distribution functions, supported by some analytical results on the collisional properties of a fermion gas. We find that the trapped gaseous mixture is close to the collisionless regime for values of the parameters corresponding to current experiments, but can be driven towards a collisional regime even without increasing the strength of the interactions, either by going over to heavier impurity masses or by matching the width of the momentum distribution of the impurities to the Fermi momentum of the fermion gas.

arXiv:cond-mat/0408528•Physical Review A•Oct 21, 2004•4 citations

Limits of sympathetic cooling of fermions by zero-temperature bosons due to particle losses

Lincoln D. Carr, Thomas Bourdel, Yvan Castin

It has been suggested by Timmermans [Phys. Rev. Lett. 87, 240403 (2001)] that loss of fermions in a degenerate system causes strong heating. We address the fundamental limit imposed by this loss on the temperature that may be obtained by sympathetic cooling of fermions by bosons. Both a quantum Boltzmann equation and a quantum Boltzmann master equation are used to study the evolution of the occupation number distribution. It is shown that, in the thermodynamic limit, the Fermi gas cools to a minimal temperature ${k}_{\mathrm{B}}T∕\ensuremath{\mu}\ensuremath{\propto}{({\ensuremath{\gamma}}_{\text{loss}}∕{\ensuremath{\gamma}}_{\text{coll}})}^{0.44}$, where ${\ensuremath{\gamma}}_{\text{loss}}$ is a constant loss rate, ${\ensuremath{\gamma}}_{\text{coll}}$ is the bare fermion-boson collision rate not including the reduction due to Fermi statistics, and $\ensuremath{\mu}\ensuremath{\sim}{k}_{\mathrm{B}}{T}_{\mathrm{F}}$ is the chemical potential. It is demonstrated that, beyond the thermodynamic limit, the discrete nature of the momentum spectrum of the system can block cooling. The unusual nonthermal nature of the number distribution is illustrated from several points of view: the Fermi surface is distorted, and in the region of zero momentum the number distribution can descend to values significantly less than unity. Our model explicitly depends on a constant evaporation rate, the value of which can strongly affect the minimum temperature.

arXiv:cond-mat/0305441•Physical Review A•Mar 9, 2004•18 citations

Exploring chemical reactions in a quantum degenerate gas of polar molecules via complex formation

Peiru He, Thomas Bilitewski, Chris H. Greene +1 more

A recent experiment [L. De Marco, G. Valtolina, K. Matsuda, W. G. Tobias, J. P. Covey, and J. Ye, A degenerate Fermi gas of polar molecules, Science 363, 853 (2019)] reported for the first time the preparation of a Fermi degenerate gas of polar molecules and observed a suppression of their chemical reaction rate compared to the one expected from a treatment assuming classical Maxwell-Boltzmann statistics. While it was hypothesized that the suppression in the ultracold regime had its roots in the Fermi statistics of the molecules, this argument is inconsistent with the fact that the Fermi pressure should set a lower bound for the chemical reaction rate. Here we develop a simple model of chemical reactions that occur via the formation and decay of molecular complexes. We indeed find that pure two-body molecule losses are unable to explain the observed suppression. Instead we extend our description beyond two-body physics by including multibody complex-molecule interactions. Although our model is able to quantitatively reproduce recent experimental observations, it requires parameters physically unlikely for direct microscopic interactions. The underlying processes, however, might emerge effectively from many-body or medium effects. A detailed understanding of the direct microscopic mechanisms responsible for these higher-order interaction processes is therefore still pending.

arXiv:2008.05904•Physical Review A•Dec 18, 2020•9 citations

Kinetics of quantum reaction-diffusion systems

Federico Gerbino, Igor Lesanovsky, Gabriele Perfetto

We discuss many-body fermionic and bosonic systems subject to dissipative particle losses in arbitrary spatial dimensions d d , within the Keldysh path-integral formulation of the quantum master equation. This open quantum dynamics represents a generalisation of classical reaction-diffusion dynamics to the quantum realm. We first show how initial conditions can be introduced in the Keldysh path integral via boundary terms. We then study binary annihilation reactions A+A\to \emptyset A + A → ∅ , for which we derive a Boltzmann-like kinetic equation. The ensuing algebraic decay in time for the particle density depends on the particle statistics. In order to model possible experimental implementations with cold atoms, for fermions in d=1 d = 1 we further discuss inhomogeneous cases involving the presence of a trapping potential. In this context, we quantify the irreversibility of the dynamics studying the time evolution of the system entropy for different quenches of the trapping potential. We find that the system entropy features algebraic decay for confining quenches, while it saturates in deconfined scenarios.

arXiv:2406.20028•SciPost Physics Core•Feb 4, 2025•4 citations

Strong Dissipation Inhibits Losses and Induces Correlations in Cold Molecular Gases

Niels Syassen, D. M. Bauer, Matthias Lettner +6 more

Atomic quantum gases in the strong-correlation regime offer unique possibilities to explore a variety of many-body quantum phenomena. Reaching this regime has usually required both strong elastic and weak inelastic interactions because the latter produce losses. We show that strong inelastic collisions can actually inhibit particle losses and drive a system into a strongly correlated regime. Studying the dynamics of ultracold molecules in an optical lattice confined to one dimension, we show that the particle loss rate is reduced by a factor of 10. Adding a lattice along the one dimension increases the reduction to a factor of 2000. Our results open the possibility to observe exotic quantum many-body phenomena with systems that suffer from strong inelastic collisions.

arXiv:0806.4310•Science•Jun 6, 2008•437 citations
Two-body inelastic collisions arising from chemical reactions are prevalent in ultracold fermionic and bosonic molecular gases. Although recent advancements have achieved quantum degeneracy in these systems, loss dynamics are typically modeled phenomenologically using rate equations that often assume thermalization during chemical reactions. In this study, we employ the inelastic quantum Boltzmann equation to analyze particle loss, temperature evolution, and momentum distributions in single-component Fermi gases from first principles. Our results demonstrate that the conventional particle-number rate equation accurately describes the dynamics in trapped systems but fails to capture the behavior in homogeneous systems. Notably, under pure p-wave inelastic collisions and zero elastic collisions, we find that systems prepared near or above quantum degeneracy remain in a thermal state, whereas systems initialized deep within degeneracy exhibit nonequilibrium dynamics. Our theoretical predictions align well with recent experimental observations in trapped systems, and our claim can be further verified in atomic systems with induced two-body loss in box potentials.
This work is a companion paper to X.-Y. Gao et al., Phys. Rev. Lett. 134, 153402 (2025), where we discuss the nonequilibrium two-body loss dynamics of a single-component ultracold Fermi gas and its possible thermalization. This paper provides detailed information on the derivation and analysis of the inelastic quantum Boltzmann equation (IQBE) used to describe the system. We demonstrate that the Mellin transform is a powerful tool for solving and approximating the IQBE for free-space systems. In this case, the particle-number dynamics are beyond the description of the widely used phenomenological two-body equation. For harmonically trapped systems, we propose a fast-flowing approximation to simplify the numerical evaluation of the IQBE. We verify the approximation in an analogous quasi-one-dimensional system and apply it to three-dimensional calculations, obtaining satisfactory agreement with recent experimental results. Furthermore, we compare the nonequilibrium results with those obtained using a thermal ansatz in both situations, providing a systematic understanding of the antievaporation phenomena observed in such systems.
We present a theoretical study of the dynamical behavior of a gas made of ultracold fermionic atoms, which during their motions can collide with a much smaller number of thermal bosonic impurities. The atoms are confined inside harmonic traps and the interactions between the two species are treated as due to sss-wave scattering with a negative scattering length modeling the $^{40}\mathrm{K}\penalty1000-\hskip0pt^{87}\mathrm{Rb}$ fermion-boson system. We set the fermions into motion by giving a small shift to their trap center and examine two alternative types of initial conditions, referring to (i) a close-to-equilibrium situation in which the two species are at the same temperature (well below the Fermi temperature and well above the Bose-Einstein condensation temperature); and (ii) a far-from-equilibrium case in which the impurities are given a Boltzmann distribution of momenta while the fermions are at very low temperatures. The dynamics of the gas is evaluated by the numerical solution of the Vlasov-Landau equations for the one-body distribution functions, supported by some analytical results on the collisional properties of a fermion gas. We find that the trapped gaseous mixture is close to the collisionless regime for values of the parameters corresponding to current experiments, but can be driven towards a collisional regime even without increasing the strength of the interactions, either by going over to heavier impurity masses or by matching the width of the momentum distribution of the impurities to the Fermi momentum of the fermion gas.
It has been suggested by Timmermans [Phys. Rev. Lett. 87, 240403 (2001)] that loss of fermions in a degenerate system causes strong heating. We address the fundamental limit imposed by this loss on the temperature that may be obtained by sympathetic cooling of fermions by bosons. Both a quantum Boltzmann equation and a quantum Boltzmann master equation are used to study the evolution of the occupation number distribution. It is shown that, in the thermodynamic limit, the Fermi gas cools to a minimal temperature ${k}_{\mathrm{B}}T∕\ensuremath{\mu}\ensuremath{\propto}{({\ensuremath{\gamma}}_{\text{loss}}∕{\ensuremath{\gamma}}_{\text{coll}})}^{0.44}$, where ${\ensuremath{\gamma}}_{\text{loss}}$ is a constant loss rate, ${\ensuremath{\gamma}}_{\text{coll}}$ is the bare fermion-boson collision rate not including the reduction due to Fermi statistics, and $\ensuremath{\mu}\ensuremath{\sim}{k}_{\mathrm{B}}{T}_{\mathrm{F}}$ is the chemical potential. It is demonstrated that, beyond the thermodynamic limit, the discrete nature of the momentum spectrum of the system can block cooling. The unusual nonthermal nature of the number distribution is illustrated from several points of view: the Fermi surface is distorted, and in the region of zero momentum the number distribution can descend to values significantly less than unity. Our model explicitly depends on a constant evaporation rate, the value of which can strongly affect the minimum temperature.
A recent experiment [L. De Marco, G. Valtolina, K. Matsuda, W. G. Tobias, J. P. Covey, and J. Ye, A degenerate Fermi gas of polar molecules, Science 363, 853 (2019)] reported for the first time the preparation of a Fermi degenerate gas of polar molecules and observed a suppression of their chemical reaction rate compared to the one expected from a treatment assuming classical Maxwell-Boltzmann statistics. While it was hypothesized that the suppression in the ultracold regime had its roots in the Fermi statistics of the molecules, this argument is inconsistent with the fact that the Fermi pressure should set a lower bound for the chemical reaction rate. Here we develop a simple model of chemical reactions that occur via the formation and decay of molecular complexes. We indeed find that pure two-body molecule losses are unable to explain the observed suppression. Instead we extend our description beyond two-body physics by including multibody complex-molecule interactions. Although our model is able to quantitatively reproduce recent experimental observations, it requires parameters physically unlikely for direct microscopic interactions. The underlying processes, however, might emerge effectively from many-body or medium effects. A detailed understanding of the direct microscopic mechanisms responsible for these higher-order interaction processes is therefore still pending.
We discuss many-body fermionic and bosonic systems subject to dissipative particle losses in arbitrary spatial dimensions d d , within the Keldysh path-integral formulation of the quantum master equation. This open quantum dynamics represents a generalisation of classical reaction-diffusion dynamics to the quantum realm. We first show how initial conditions can be introduced in the Keldysh path integral via boundary terms. We then study binary annihilation reactions A+A\to \emptyset A + A → ∅ , for which we derive a Boltzmann-like kinetic equation. The ensuing algebraic decay in time for the particle density depends on the particle statistics. In order to model possible experimental implementations with cold atoms, for fermions in d=1 d = 1 we further discuss inhomogeneous cases involving the presence of a trapping potential. In this context, we quantify the irreversibility of the dynamics studying the time evolution of the system entropy for different quenches of the trapping potential. We find that the system entropy features algebraic decay for confining quenches, while it saturates in deconfined scenarios.
Atomic quantum gases in the strong-correlation regime offer unique possibilities to explore a variety of many-body quantum phenomena. Reaching this regime has usually required both strong elastic and weak inelastic interactions because the latter produce losses. We show that strong inelastic collisions can actually inhibit particle losses and drive a system into a strongly correlated regime. Studying the dynamics of ultracold molecules in an optical lattice confined to one dimension, we show that the particle loss rate is reduced by a factor of 10. Adding a lattice along the one dimension increases the reduction to a factor of 2000. Our results open the possibility to observe exotic quantum many-body phenomena with systems that suffer from strong inelastic collisions.