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

Physical Review Letters

Relaxation of Phonons in the Lieb-Liniger Gas by Dynamical Refermionization

Isabelle Bouchoule, Jérôme Dubail, Léa Dubois, Dimitri M. Gangardt

Published Apr 4, 2023• 12 citations
Publisher

Abstract

Motivated by recent experiments, we investigate the Lieb-Liniger gas initially prepared in an out-of-equilibrium state that is Gaussian in terms of the phonons, namely whose density matrix is the exponential of an operator quadratic in terms of phonon creation and annihilation operators. Because the phonons are not exact eigenstates of the Hamiltonian, the gas relaxes to a stationary state at very long times whose phonon population is a priori different from the initial one. Thanks to integrability, that stationary state needs not be a thermal state. Using the Bethe-ansatz mapping between the exact eigenstates of the Lieb-Liniger Hamiltonian and those of a noninteracting Fermi gas and bosonization techniques we completely characterize the stationary state of the gas after relaxation and compute its phonon population distribution. We apply our results to the case where the initial state is an excited coherent state for a single phonon mode, and we compare them to exact results obtained in the hard-core limit.

DOI
10.1103/physrevlett.130.140401

Related papers

Relaxation of phonons in the Lieb-Liniger gas by dynamical refermionization

Isabelle Bouchoule, Jérôme Dubail, Léa Dubois +1 more

We investigate the Lieb-Liniger gas initially prepared in an out-of-equilibrium state that is Gaussian in terms of the phonons. Because the phonons are not exact eigenstates of the Hamiltonian, the gas relaxes to a stationary state at very long times. Thanks to integrability, that stationary state needs not be a thermal state. We characterize the stationary state of the gas after relaxation and compute its phonon population distribution. Technically, this follows from the mapping between the exact eigenstates of the Lieb-Liniger Hamiltonian and those of a non-interacting Fermi gas -- a mapping provided by the Bethe equations -- , as well as on bosonization formulas valid in the low-energy sector of the Hilbert space. We apply our results to the case where the initial state is an excited coherent state for a single phonon mode, and we compare them to exact results obtained in the hard-core limit.

arXiv:2206.00112•May 31, 2022•0 citations

Relaxation dynamics of the Lieb-Liniger gas following an interaction quench: A coordinate Bethe-ansatz analysis

Jan C. Zill, Tod M. Wright, Karen V. Kheruntsyan +2 more

We investigate the relaxation dynamics of the integrable Lieb-Liniger model of contact-interacting bosons in one dimension following a sudden quench of the collisional interaction strength. The system is initially prepared in its noninteracting ground state and the interaction strength is then abruptly switched to a positive value, corresponding to repulsive interactions between the bosons. We calculate equal-time correlation functions of the nonequilibrium Bose field for small systems of up to five particles via symbolic evaluation of coordinate Bethe-ansatz expressions for operator matrix elements between Lieb-Liniger eigenstates. We characterize the relaxation of the system by comparing the time-evolving correlation functions following the quench to the equilibrium correlations predicted by the diagonal ensemble and relate the behavior of these correlations to that of the quantum fidelity between the many-body wave function and the initial state of the system. Our results for the asymptotic scaling of local second-order correlations with increasing interaction strength agree with the predictions of recent generalized thermodynamic Bethe-ansatz calculations. By contrast, third-order correlations obtained within our approach exhibit a markedly different power-law dependence on the interaction strength as the Tonks-Girardeau limit of infinitely strong interactions is approached.

arXiv:1407.4998•Physical Review A•Feb 12, 2015•43 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

Relaxation dynamics of an exactly solvable electron-phonon model

Dante Marvin Kennes, V. Meden

We address the question whether observables of an exactly solvable model of electrons coupled to (optical) phonons relax into large time stationary state values and investigate if the asymptotic expectation values can be computed using a stationary density matrix. Two initial nonequilibrium situations are considered. A sudden quench of the electron-phonon coupling, starting from the noninteracting canonical equilibrium at temperature $T$ in the electron as well as in the phonon subsystems, leads to a rather simple dynamics. A richer time evolution emerges if the initial state is taken as the product of the phonon vacuum and the filled Fermi sea supplemented by a highly excited additional electron. Our model has a natural set of constants of motion, with as many elements as degrees of freedom. In accordance with earlier studies of such type of models, we find that expectation values which become stationary can be described by the density matrix of a generalized Gibbs ensemble which differs from that of a canonical ensemble. For the model at hand, it appears to be evident that the eigenmode occupancy operators should be used in the construction of the stationary density matrix.

arXiv:1006.0927•Physical Review B•Aug 12, 2010•23 citations

Gaussian-state Ansatz for the non-equilibrium dynamics of quantum spin lattices

Raphaël Menu, Tommaso Roscilde

The study of non-equilibrium dynamics is one of the most important challenges of modern quantum many-body physics, in relationship with fundamental questions in quantum statistical mechanics, as well as with the fields of quantum simulation and computing. In this work we propose a Gaussian Ansatz for the study of the nonequilibrium dynamics of quantum spin systems. Within our approach, the quantum spins are mapped onto Holstein-Primakoff bosons, such that a coherent spin state - chosen as the initial state of the dynamics - represents the bosonic vacuum. The state of the system is then postulated to remain a bosonic Gaussian state at all times, an assumption which is exact when the bosonic Hamiltonian is quadratic; and which is justified in the case of a nonlinear Hamiltonian if the boson density remains moderate. We test the accuracy of such an Ansatz in the paradigmatic case of the S=1/2 S = 1 / 2 transverse-field Ising model, in one and two dimensions, initialized in a state aligned with the applied field. We show that the Gaussian Ansatz, when applied to the bosonic Hamiltonian with nonlinearities truncated to quartic order, is able to reproduce faithfully the evolution of the state, including its relaxation to the equilibrium regime, for fields larger than the critical field for the paramagnetic-ferromagnetic transition in the ground state. In particular the spatio-temporal pattern of correlations reconstructed via the Gaussian Ansatz reveals the dispersion relation of quasiparticle excitations, exhibiting the softening of the excitation gap upon approaching the critical field. Our results suggest that the Gaussian Ansatz correctly captures the essential effects of nonlinearities in quantum spin dynamics; and that it can be applied to the study of fundamental phenomena such as quantum thermalization and its breakdown.

arXiv:2301.01363•SciPost Physics•Jun 12, 2023•8 citations

Dynamics of the attractive 1D Bose gas: analytical treatment from integrability

Pasquale Calabrese, Jean-Sébastien Caux

The physics of the attractive one-dimensional Bose gas (Lieb–Liniger model) is investigated with techniques based on the integrability of the system. Combining a knowledge of particle quasi-momenta to exponential precision in the system size with determinant representations of matrix elements of local operators coming from the algebraic Bethe ansatz, we obtain rather general analytical results for the zero-temperature dynamical correlation functions of the density and field operators. Our results thus provide quantitative predictions for possible future experiments in atomic gases or optical waveguides.

arXiv:0707.4115•Journal of Statistical Mechanics Theory and Experiment•Aug 28, 2007•51 citations
Motivated by recent experiments, we investigate the Lieb-Liniger gas initially prepared in an out-of-equilibrium state that is Gaussian in terms of the phonons, namely whose density matrix is the exponential of an operator quadratic in terms of phonon creation and annihilation operators. Because the phonons are not exact eigenstates of the Hamiltonian, the gas relaxes to a stationary state at very long times whose phonon population is a priori different from the initial one. Thanks to integrability, that stationary state needs not be a thermal state. Using the Bethe-ansatz mapping between the exact eigenstates of the Lieb-Liniger Hamiltonian and those of a noninteracting Fermi gas and bosonization techniques we completely characterize the stationary state of the gas after relaxation and compute its phonon population distribution. We apply our results to the case where the initial state is an excited coherent state for a single phonon mode, and we compare them to exact results obtained in the hard-core limit.
We investigate the Lieb-Liniger gas initially prepared in an out-of-equilibrium state that is Gaussian in terms of the phonons. Because the phonons are not exact eigenstates of the Hamiltonian, the gas relaxes to a stationary state at very long times. Thanks to integrability, that stationary state needs not be a thermal state. We characterize the stationary state of the gas after relaxation and compute its phonon population distribution. Technically, this follows from the mapping between the exact eigenstates of the Lieb-Liniger Hamiltonian and those of a non-interacting Fermi gas -- a mapping provided by the Bethe equations -- , as well as on bosonization formulas valid in the low-energy sector of the Hilbert space. We apply our results to the case where the initial state is an excited coherent state for a single phonon mode, and we compare them to exact results obtained in the hard-core limit.
We investigate the relaxation dynamics of the integrable Lieb-Liniger model of contact-interacting bosons in one dimension following a sudden quench of the collisional interaction strength. The system is initially prepared in its noninteracting ground state and the interaction strength is then abruptly switched to a positive value, corresponding to repulsive interactions between the bosons. We calculate equal-time correlation functions of the nonequilibrium Bose field for small systems of up to five particles via symbolic evaluation of coordinate Bethe-ansatz expressions for operator matrix elements between Lieb-Liniger eigenstates. We characterize the relaxation of the system by comparing the time-evolving correlation functions following the quench to the equilibrium correlations predicted by the diagonal ensemble and relate the behavior of these correlations to that of the quantum fidelity between the many-body wave function and the initial state of the system. Our results for the asymptotic scaling of local second-order correlations with increasing interaction strength agree with the predictions of recent generalized thermodynamic Bethe-ansatz calculations. By contrast, third-order correlations obtained within our approach exhibit a markedly different power-law dependence on the interaction strength as the Tonks-Girardeau limit of infinitely strong interactions is approached.
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 address the question whether observables of an exactly solvable model of electrons coupled to (optical) phonons relax into large time stationary state values and investigate if the asymptotic expectation values can be computed using a stationary density matrix. Two initial nonequilibrium situations are considered. A sudden quench of the electron-phonon coupling, starting from the noninteracting canonical equilibrium at temperature TTT in the electron as well as in the phonon subsystems, leads to a rather simple dynamics. A richer time evolution emerges if the initial state is taken as the product of the phonon vacuum and the filled Fermi sea supplemented by a highly excited additional electron. Our model has a natural set of constants of motion, with as many elements as degrees of freedom. In accordance with earlier studies of such type of models, we find that expectation values which become stationary can be described by the density matrix of a generalized Gibbs ensemble which differs from that of a canonical ensemble. For the model at hand, it appears to be evident that the eigenmode occupancy operators should be used in the construction of the stationary density matrix.
The study of non-equilibrium dynamics is one of the most important challenges of modern quantum many-body physics, in relationship with fundamental questions in quantum statistical mechanics, as well as with the fields of quantum simulation and computing. In this work we propose a Gaussian Ansatz for the study of the nonequilibrium dynamics of quantum spin systems. Within our approach, the quantum spins are mapped onto Holstein-Primakoff bosons, such that a coherent spin state - chosen as the initial state of the dynamics - represents the bosonic vacuum. The state of the system is then postulated to remain a bosonic Gaussian state at all times, an assumption which is exact when the bosonic Hamiltonian is quadratic; and which is justified in the case of a nonlinear Hamiltonian if the boson density remains moderate. We test the accuracy of such an Ansatz in the paradigmatic case of the S=1/2 S = 1 / 2 transverse-field Ising model, in one and two dimensions, initialized in a state aligned with the applied field. We show that the Gaussian Ansatz, when applied to the bosonic Hamiltonian with nonlinearities truncated to quartic order, is able to reproduce faithfully the evolution of the state, including its relaxation to the equilibrium regime, for fields larger than the critical field for the paramagnetic-ferromagnetic transition in the ground state. In particular the spatio-temporal pattern of correlations reconstructed via the Gaussian Ansatz reveals the dispersion relation of quasiparticle excitations, exhibiting the softening of the excitation gap upon approaching the critical field. Our results suggest that the Gaussian Ansatz correctly captures the essential effects of nonlinearities in quantum spin dynamics; and that it can be applied to the study of fundamental phenomena such as quantum thermalization and its breakdown.
The physics of the attractive one-dimensional Bose gas (Lieb–Liniger model) is investigated with techniques based on the integrability of the system. Combining a knowledge of particle quasi-momenta to exponential precision in the system size with determinant representations of matrix elements of local operators coming from the algebraic Bethe ansatz, we obtain rather general analytical results for the zero-temperature dynamical correlation functions of the density and field operators. Our results thus provide quantitative predictions for possible future experiments in atomic gases or optical waveguides.