We study the relaxation of a highly collisional, ultracold but nondegenerate gas of polar molecules. Confined within a harmonic trap, the gas is subject to fluid-gaseous coupled dynamics that lead to a breakdown of first-order hydrodynamics. An attempt to treat these higher-order hydrodynamic effects was previously made with a Gaussian ansatz and coarse-graining model parameter [R. R. W. Wang and J. L. Bohn, Phys. Rev. A 108, 013322 (2023)], leading to an approximate set of equations for a few collective observables accessible to experiments. Here we present substantially improved reduced-order models for these same observables, admissible beyond previous parameter regimes, discovered directly from particle simulations using the WSINDy algorithm (Weak-form Sparse Identification of Nonlinear Dynamics). The interpretable nature of the learning algorithm enables estimation of previously unknown physical quantities and discovery of model terms with candidate physical mechanisms, revealing new physics in mixed collisional regimes. Our approach constitutes a general framework for data-driven model identification leveraging known physics.
We study the relaxation of a highly collisional, ultracold but nondegenerate gas of polar molecules. Confined within a harmonic trap, the gas is subject to fluid-gaseous coupled dynamics that lead to a breakdown of first-order hydrodynamics. An attempt to treat these higher-order hydrodynamic effects was previously made with a Gaussian ansatz and coarse-graining model parameter [R. R. W. Wang and J. L. Bohn, Phys. Rev. A 108, 013322 (2023)], leading to an approximate set of equations for a few collective observables accessible to experiments. Here we present substantially improved reduced-order models for these same observables, admissible beyond previous parameter regimes, discovered directly from particle simulations using the WSINDy algorithm (Weak-form Sparse Identification of Nonlinear Dynamics). The interpretable nature of the learning algorithm enables estimation of previously unknown physical quantities and discovery of model terms with candidate physical mechanisms, revealing new physics in mixed collisional regimes. Our approach constitutes a general framework for data-driven model identification leveraging known physics.
Simulating and predicting dynamics of quantum many-body systems is extremely challenging, even for state-of-the-art computational methods, due to the spread of entanglement across the system. However, in the long-wavelength limit, quantum systems often admit a simplified description, which involves a small set of physical observables and requires only a few parameters such as sound velocity or viscosity. Unveiling the relationship between these hydrodynamic equations and the underlying microscopic theory usually requires a great effort by condensed matter theorists. In the present paper, we develop a new machine-learning framework for automated discovery of effective equations from a limited set of available data, thus bypassing complicated analytical derivations. The data can be generated from numerical simulations or come from experimental quantum simulator platforms. Using integrable models, where direct comparisons can be made, we reproduce previously known hydrodynamic equations, strikingly discover novel equations and provide their derivation whenever possible. We discover new hydrodynamic equations describing dynamics of interacting systems, for which the derivation remains an outstanding challenge. Our approach provides a new interpretable method to study properties of quantum materials and quantum simulators in non-perturbative regimes.
In Chapters 11–13, we gave a detailed discussion of the dynamics of a trapped Bose gas at finite temperatures in a region where the collisions described by the C 12 and C 22 terms in the kinetic equation (3.42) do not play the central role. In this “collisionless” region, the dominant interaction effects are associated with the self-consistent fields which both the condensate and noncondensate atoms feel. Thus the dynamics can be understood to a first approximation by neglecting the C 12 and C 22 collision integrals in the kinetic equation and, at the next stage, treating them as a weak perturbation on the collisionless dynamics. In the rest of this book (Chapters 15–19), we turn to the study of the coupled ZNG equations in the opposite limit, where the C 12 and C 22 collision integrals completely determine the dynamics of the thermal cloud. Specifically, the collisions lead to the thermal cloud being in local hydrodynamic equilibrium, and hence this regime is described by the equations of collisional hydrodynamics. Its characteristic feature is that the nonequilibrium behaviour of the thermal cloud atoms can be completely described in terms of a few differential equations involving coarse-grained variables that are dependent on position and time, analogous to the condensate variables n c ( r , t ) and v c ( r , t ). In the present chapter, assuming that the thermal cloud distribution function f ( p, r , t ) is given by the Bose distribution describing partial local equilibrium, (15.16), we show how the ZNG coupled equations lead precisely to Landau's two-fluid equations, reviewed in Chapter 14. This equivalence is not obvious, mainly because Landau's equations are expressed in terms of thermodynamic variables, which are not used in a more microscopic analysis such as that used in the ZNG approach.
Describing and understanding the motion of quantum gases out of equilibrium is one of the most important modern challenges for theorists. In the groundbreaking Quantum Newton Cradle experiment [Kinoshita, Wenger and Weiss, Nature 440, 900 (2006)], quasi-one-dimensional cold atom gases were observed with unprecedented accuracy, providing impetus for many developments on the effects of low dimensionality in out-of-equilibrium physics. But it is only recently that the theory of generalized hydrodynamics has provided the adequate tools for a numerically efficient description. Using it, we give a complete numerical study of the time evolution of an ultracold atomic gas in this setup, in an interacting parameter regime close to that of the original experiment. We evaluate the full evolving phase-space distribution of particles. We simulate oscillations due to the harmonic trap, the collision of clouds without thermalization, and observe a small elongation of the actual oscillation period and cloud deformations due to many-body dephasing. We also analyze the effects of weak anharmonicity. In the experiment, measurements are made after release from the one-dimensional trap. We evaluate the gas density curves after such a release, characterizing the actual time necessary for reaching the asymptotic state where the integrable quasi-particle momentum distribution function emerges.
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.
Ultracold bosonic and fermionic quantum gases confined to quasi-one-dimensional (1D) geometry are promising candidates for probing fundamental concepts of Luttinger liquid (LL) physics. They can also be exploited for devising applications in quantum information processing and precision measurements. Here, we focus on 1D dipolar Bose gases, where evidence of super-strong coupling behavior has been demonstrated by analyzing the low-energy static and dynamical structures of the fluid at zero temperature by a combined reptation quantum Monte Carlo (RQMC) and bosonization approach. Fingerprints of LL behavior emerge in the whole crossover from the already strongly interacting Tonks–Girardeau at low density to a dipolar density wave regime at high density. We have also shown that a LL framework can be effectively set up and utilized to describe this strongly correlated crossover physics in the case of confined 1D geometries after using the results for the homogeneous system in LL hydrodynamic equations within a local density approximation. This leads to the prediction of observable quantities such as the frequencies of the collective modes of the trapped dipolar gas under the more realistic conditions that could be found in ongoing experiments. The present paper provides a description of the theoretical framework in which the above results have been worked out, making available all the detailed derivations of the hydrodynamic Luttinger equations for the inhomogeneous trapped gas and of the correlation functions for the homogeneous system.
In the collisional region at finite temperatures, the collective modes of superfluids are described by the Landau two-fluid hydrodynamic equations. This region can now be probed over the entire BCS--Bose-Einstein-condensate crossover in trapped Fermi superfluids with a Feshbach resonance, including the unitarity region. Building on the approach initiated by Zaremba, Nikuni, and Griffin in 1999 for trapped atomic Bose gases, we present a variational formulation of two-fluid hydrodynamic collective modes based on the work of Zilsel in 1950 developed for superfluid helium. Assuming a simple variational Ansatz for the superfluid and normal fluid velocities, the frequencies of the hydrodynamic modes are given by solutions of coupled algebraic equations, with constants only involving spatial integrals over various equilibrium thermodynamic derivatives. This variational approach is both simpler and more physical than a direct attempt to solve the Landau two-fluid differential equations. Our two-fluid results are shown to reduce to those of Pitaevskii and Stringari for a pure superfluid at T=0.