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:2002.01423 · SciPost Physics

Exact out-of-equilibrium steady states in the semiclassical limit of the interacting Bose gas

Giuseppe Del Vecchio Del Vecchio, Alvise Bastianello, Andrea De Luca, Giuseppe Mussardo

Published Jul 6, 2020• 27 citations
arXivHTMLPDF

Abstract

We study the out-of-equilibrium properties of a classical integrable non-relativistic theory, with a time evolution initially prepared with a finite energy density in the thermodynamic limit. The theory considered here is the Non-Linear Schrödinger equation which describes the dynamics of the one-dimensional interacting Bose gas in the regime of high occupation numbers. The main emphasis is on the determination of the late-time Generalised Gibbs Ensemble (GGE), which can be efficiently semi-numerically computed on arbitrary initial states, completely solving the famous quench problem in the classical regime. We take advantage of known results in the quantum model and the semiclassical limit to achieve new exact results for the momenta of the density operator on arbitrary GGEs, which we successfully compare with ab-initio numerical simulations. Furthermore, we determine the whole probability distribution of the density operator (full counting statistics), whose exact expression is still out of reach in the quantum model.

DOI
10.21468/scipostphys.9.1.002

Related papers

Interaction quenches in the one-dimensional Bose gas

Márton Kormos, Aditya Shashi, Yang-Zhi Chou +2 more

The nonequilibrium dynamics of integrable systems are highly constrained by the conservation of certain charges. There is substantial evidence that after a quantum quench they do not thermalize but their asymptotic steady state can be described by a generalized Gibbs ensemble (GGE) built from the conserved charges. Most of the studies on the GGE so far have focused on models that can be mapped to quadratic systems, while analytic treatment in nonquadratic systems remained elusive. We obtain results on interaction quenches in a nonquadratic continuum system, the one-dimensional (1D) Bose gas described by the integrable Lieb-Liniger model. The direct implementation of the GGE prescription is prohibited by the divergence of the conserved charges, which we conjecture to be endemic to any continuum integrable systems with contact interactions undergoing a sudden quench. We compute local correlators for a noninteracting initial state and arbitrary final interactions as well as two-point functions for quenches to the Tonks-Girardeau regime. We show that in the long time limit integrability leads to significant deviations from the predictions of the grand canonical ensemble, allowing for an experimental verification in cold-atom systems.

arXiv:1305.7202•Physical Review B•Nov 22, 2013•136 citations

Interaction quenches in the Lieb-Liniger model

Márton Kormos, Aditya Shashi, Yang-Zhi Chou +1 more

We obtain exact results on interaction quenches in the 1D Bose gas described by the integrable Lieb-Liniger model. We show that in the long time limit integrability leads to significant deviations from the predictions of the grand canonical ensemble and a description within the generalized Gibbs ensemble (GGE) is needed. For a non-interacting initial state and arbitrary final interactions, we find that the presence of infinitely many conserved charges generates a non-analytic behavior in the equilibrated density of quasimomenta. This manifests itself in a dynamically generated Friedel-like oscillation of the non-local correlation functions with interaction dependent oscillation momenta. We also exactly evaluate local correlations and the generalized chemical potentials within GGE.

arXiv:1204.3889•Apr 17, 2012•4 citations

An exact stochastic field method for the interacting Bose gas at thermal equilibrium

Iacopo Carusotto, Yvan Castin

We present a new exact method to numerically compute the thermodynamical properties of an interacting Bose gas in the canonical ensemble. As in our previous paper (Carusotto I, Castin Y and Dalibard J 2001 Phys. Rev. A 63 023606), we write the density operator ρ as an average of Hartree dyadics | N :ϕ 1 ⟩⟨ N :ϕ 2 | and we find stochastic evolution equations for the wavefunctions ϕ 1,2 such that the exact imaginary-time evolution of ρ is recovered after averaging over noise. In this way, the thermal equilibrium density operator can be obtained for any temperature T . The method is then applied to study the thermodynamical properties of a homogeneous one-dimensional N -boson system: although Bose-Einstein condensation cannot occur in the thermodynamical limit, a macroscopic occupation of the lowest mode of a finite system is observed at sufficiently low temperatures. If k B T >>µ, the main effect of interactions is to suppress density fluctuations and to reduce their correlation length. Different effects such as a spatial antibunching of the atoms are predicted for the opposite k B T ⩽µ regime. Our exact stochastic calculations have been compared with existing approximate theories.

arXiv:cond-mat/0108042•Journal of Physics B Atomic Molecular and Optical Physics•Nov 23, 2001•37 citations

Analytic results for a quantum quench from free to hard-core one-dimensional bosons

Márton Kormos, Mario Collura, Pasquale Calabrese

It is widely believed that the stationary properties after a quantum quench in integrable systems can be described by a generalized Gibbs ensemble (GGE), even if all of the analytical evidence is based on free theories in which the pre- and postquench modes are linearly related. In contrast, we consider the experimentally relevant quench of the one-dimensional Bose gas from zero to infinite interaction, in which the relation between modes is nonlinear, and consequently Wick's theorem does not hold. We provide exact analytical results for the time evolution of the dynamical density-density correlation function at any time after the quench and we prove that its stationary value is described by a GGE in which Wick's theorem is restored.

arXiv:1307.2142•Physical Review A•Jan 13, 2014•173 citations

Nonequilibrium dynamics of a one-dimensional Bose gas via the coordinate Bethe ansatz

Jan C. Zill

This thesis is concerned with non-equilibrium phenomena in interacting quantum many-body systems. Specifically, we investigate the time-evolution and relaxation dynamics of Bose gases in a highly restricted geometry, which constrains the dynamics to one spatial dimension. This leads to a description of the system in terms of a simple model Hamiltonian which permits exact many-body quantum mechanical solutions due to its integrability. In the first part of this thesis, a computational method is developed to obtain experimentally relevant correlation functions in the framework of the coordinate Bethe ansatz. We employ this method to compute exact ground-state correlation functions of the Lieb-Liniger gas for up to seven particles covering the whole regime of repulsive interactions. We also investigate the dynamics of the system after an instantaneous change of the interaction strength. This quantum quench deposits large amounts of energy that cannot be dissipated due to the system being closed, and so the dynamics far from equilibrium are probed. We prepare the system in two different initial states, and quench to the same final interaction strength. The latter is determined in such a way that the added energy due to the quench is the same for both scenarios. Conventional statistical mechanics predicts the same relaxed state, but due to the integrability of the system all considered correlation functions of the relaxed states differ from each other and also from the thermal ones. We then investigate the dynamics and relaxed state for a quench from zero to repulsive interactions in more detail, focussing on the mechanism of relaxation and the involved time-scales. We find that local correlation functions relax on time-scales determined by the interaction strength, in contrast to non-local correlation functions, whose relaxation time-scale is proportional to the system size. Next, we employ the same methodology to study the one-dimensional Bose gas with attractive interactions. In this case many-body bound states are permissible solutions of the Lieb-Liniger model. We compare exact ground-state correlation functions of up to seven particles to their corresponding mean-field solution. The latter displays a quantum phase transition at a critical interaction strength, marking the transition from a uniform-density state to a localized bright soliton. Our exact results agree remarkably well with the corresponding mean-field solution past the critical point. We also investigate the dynamics following an interaction strength quench, starting again from the non-interacting ground state. Bound states strongly influence correlation functions for all post-quench interaction strengths, and local correlation functions are largely increased compared to their initial value. In the last part of this thesis, we investigate the behavior of the one-dimensional Bose gas under periodic driving of the interaction strength. To this end, we extend the coordinate Bethe-ansatz formalism by employing Floquet theory to obtain solutions for the full time- dependent Hamiltonian. This realizes an extension to non-integrable systems in a way that allows controlled breaking of integrability. We compare the dynamics of the system following an interaction quench to a fixed final interaction strength to that of a system with periodically modulated post-quench interaction strength. For fast driving, the system evolves according to the time-averaged Hamiltonian and correlation functions are nearly identical to those of the undriven system. The response of the system qualitatively changes at resonances of the time-averaged Hamiltonian, where we observe energy absorption and a marked change in correlation functions. However, the system does not absorb energy indefinitely.

•Jan 30, 2017•0 citations

Exact time evolution of space- and time-dependent correlation functions after an interaction quench in the one-dimensional Bose gas

Jorn Mossel, Jean-Sébastien Caux

We consider the non-equilibrium dynamics of the interacting Lieb–Liniger gas after instantaneously switching the interactions off. The subsequent time evolution of the space- and time-dependent correlation functions is computed exactly. Different relaxation behavior is observed for different correlation functions. The long time average is compared with the predictions of several statistical ensembles. The generalized Gibbs ensemble restricted to a fixed number of particles is shown to give correct results at large times for all length scales.

arXiv:1201.1885•New Journal of Physics•Jul 3, 2012•104 citations
We study the out-of-equilibrium properties of a classical integrable non-relativistic theory, with a time evolution initially prepared with a finite energy density in the thermodynamic limit. The theory considered here is the Non-Linear Schrödinger equation which describes the dynamics of the one-dimensional interacting Bose gas in the regime of high occupation numbers. The main emphasis is on the determination of the late-time Generalised Gibbs Ensemble (GGE), which can be efficiently semi-numerically computed on arbitrary initial states, completely solving the famous quench problem in the classical regime. We take advantage of known results in the quantum model and the semiclassical limit to achieve new exact results for the momenta of the density operator on arbitrary GGEs, which we successfully compare with ab-initio numerical simulations. Furthermore, we determine the whole probability distribution of the density operator (full counting statistics), whose exact expression is still out of reach in the quantum model.
The nonequilibrium dynamics of integrable systems are highly constrained by the conservation of certain charges. There is substantial evidence that after a quantum quench they do not thermalize but their asymptotic steady state can be described by a generalized Gibbs ensemble (GGE) built from the conserved charges. Most of the studies on the GGE so far have focused on models that can be mapped to quadratic systems, while analytic treatment in nonquadratic systems remained elusive. We obtain results on interaction quenches in a nonquadratic continuum system, the one-dimensional (1D) Bose gas described by the integrable Lieb-Liniger model. The direct implementation of the GGE prescription is prohibited by the divergence of the conserved charges, which we conjecture to be endemic to any continuum integrable systems with contact interactions undergoing a sudden quench. We compute local correlators for a noninteracting initial state and arbitrary final interactions as well as two-point functions for quenches to the Tonks-Girardeau regime. We show that in the long time limit integrability leads to significant deviations from the predictions of the grand canonical ensemble, allowing for an experimental verification in cold-atom systems.
We obtain exact results on interaction quenches in the 1D Bose gas described by the integrable Lieb-Liniger model. We show that in the long time limit integrability leads to significant deviations from the predictions of the grand canonical ensemble and a description within the generalized Gibbs ensemble (GGE) is needed. For a non-interacting initial state and arbitrary final interactions, we find that the presence of infinitely many conserved charges generates a non-analytic behavior in the equilibrated density of quasimomenta. This manifests itself in a dynamically generated Friedel-like oscillation of the non-local correlation functions with interaction dependent oscillation momenta. We also exactly evaluate local correlations and the generalized chemical potentials within GGE.
We present a new exact method to numerically compute the thermodynamical properties of an interacting Bose gas in the canonical ensemble. As in our previous paper (Carusotto I, Castin Y and Dalibard J 2001 Phys. Rev. A 63 023606), we write the density operator ρ as an average of Hartree dyadics | N :ϕ 1 ⟩⟨ N :ϕ 2 | and we find stochastic evolution equations for the wavefunctions ϕ 1,2 such that the exact imaginary-time evolution of ρ is recovered after averaging over noise. In this way, the thermal equilibrium density operator can be obtained for any temperature T . The method is then applied to study the thermodynamical properties of a homogeneous one-dimensional N -boson system: although Bose-Einstein condensation cannot occur in the thermodynamical limit, a macroscopic occupation of the lowest mode of a finite system is observed at sufficiently low temperatures. If k B T >>µ, the main effect of interactions is to suppress density fluctuations and to reduce their correlation length. Different effects such as a spatial antibunching of the atoms are predicted for the opposite k B T ⩽µ regime. Our exact stochastic calculations have been compared with existing approximate theories.
It is widely believed that the stationary properties after a quantum quench in integrable systems can be described by a generalized Gibbs ensemble (GGE), even if all of the analytical evidence is based on free theories in which the pre- and postquench modes are linearly related. In contrast, we consider the experimentally relevant quench of the one-dimensional Bose gas from zero to infinite interaction, in which the relation between modes is nonlinear, and consequently Wick's theorem does not hold. We provide exact analytical results for the time evolution of the dynamical density-density correlation function at any time after the quench and we prove that its stationary value is described by a GGE in which Wick's theorem is restored.
This thesis is concerned with non-equilibrium phenomena in interacting quantum many-body systems. Specifically, we investigate the time-evolution and relaxation dynamics of Bose gases in a highly restricted geometry, which constrains the dynamics to one spatial dimension. This leads to a description of the system in terms of a simple model Hamiltonian which permits exact many-body quantum mechanical solutions due to its integrability. In the first part of this thesis, a computational method is developed to obtain experimentally relevant correlation functions in the framework of the coordinate Bethe ansatz. We employ this method to compute exact ground-state correlation functions of the Lieb-Liniger gas for up to seven particles covering the whole regime of repulsive interactions. We also investigate the dynamics of the system after an instantaneous change of the interaction strength. This quantum quench deposits large amounts of energy that cannot be dissipated due to the system being closed, and so the dynamics far from equilibrium are probed. We prepare the system in two different initial states, and quench to the same final interaction strength. The latter is determined in such a way that the added energy due to the quench is the same for both scenarios. Conventional statistical mechanics predicts the same relaxed state, but due to the integrability of the system all considered correlation functions of the relaxed states differ from each other and also from the thermal ones. We then investigate the dynamics and relaxed state for a quench from zero to repulsive interactions in more detail, focussing on the mechanism of relaxation and the involved time-scales. We find that local correlation functions relax on time-scales determined by the interaction strength, in contrast to non-local correlation functions, whose relaxation time-scale is proportional to the system size. Next, we employ the same methodology to study the one-dimensional Bose gas with attractive interactions. In this case many-body bound states are permissible solutions of the Lieb-Liniger model. We compare exact ground-state correlation functions of up to seven particles to their corresponding mean-field solution. The latter displays a quantum phase transition at a critical interaction strength, marking the transition from a uniform-density state to a localized bright soliton. Our exact results agree remarkably well with the corresponding mean-field solution past the critical point. We also investigate the dynamics following an interaction strength quench, starting again from the non-interacting ground state. Bound states strongly influence correlation functions for all post-quench interaction strengths, and local correlation functions are largely increased compared to their initial value. In the last part of this thesis, we investigate the behavior of the one-dimensional Bose gas under periodic driving of the interaction strength. To this end, we extend the coordinate Bethe-ansatz formalism by employing Floquet theory to obtain solutions for the full time- dependent Hamiltonian. This realizes an extension to non-integrable systems in a way that allows controlled breaking of integrability. We compare the dynamics of the system following an interaction quench to a fixed final interaction strength to that of a system with periodically modulated post-quench interaction strength. For fast driving, the system evolves according to the time-averaged Hamiltonian and correlation functions are nearly identical to those of the undriven system. The response of the system qualitatively changes at resonances of the time-averaged Hamiltonian, where we observe energy absorption and a marked change in correlation functions. However, the system does not absorb energy indefinitely.
We consider the non-equilibrium dynamics of the interacting Lieb–Liniger gas after instantaneously switching the interactions off. The subsequent time evolution of the space- and time-dependent correlation functions is computed exactly. Different relaxation behavior is observed for different correlation functions. The long time average is compared with the predictions of several statistical ensembles. The generalized Gibbs ensemble restricted to a fixed number of particles is shown to give correct results at large times for all length scales.