arXiv:2206.00112

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

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

Published May 31, 2022• 0 citations

Abstract

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.

DOI: 10.48550/arxiv.2206.00112

Related papers

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

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

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.

Physical Review LettersApr 4, 202312 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.4998Physical Review AFeb 12, 201543 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, 20170 citations

Non-Equilibrium Steady State of the Lieb-Liniger model: exact treatment of the Tonks Girardeau limit

Spyros Sotiriadis

Aiming at studying the emergence of Non-Equilibrium Steady States (NESS) in quantum integrable models by means of an exact analytical method, we focus on the Tonks-Girardeau or hard-core boson limit of the Lieb-Liniger model. We consider the abrupt expansion of a gas from one half to the entire confining box, a prototypical case of inhomogeneous quench, also known as "geometric quench". Based on the exact calculation of quench overlaps, we develop an analytical method for the derivation of the NESS by rigorously treating the thermodynamic and large time and distance limit. Our method is based on complex analysis tools for the derivation of the asymptotics of the many-body wavefunction, does not make essential use of the effectively non-interacting character of the hard-core boson gas and is sufficiently robust for generalisation to the genuinely interacting case.

arXiv:2007.12683Jul 24, 20202 citations

Non-Equilibrium Steady State of the Lieb-Liniger model: multiple-integral representation of the time evolved many-body wave-function

Spyros Sotiriadis

We continue our study of the emergence of Non-Equilibrium Steady States in quantum integrable models focusing on the expansion of a Lieb-Liniger gas for arbitrary repulsive interaction. As a first step towards the derivation of the asymptotics of observables in the thermodynamic and large distance and time limit, we derive an exact multiple integral representation of the time evolved many-body wave-function. Starting from the known but complicated expression for the overlaps of the initial state of a geometric quench, which are derived from the Slavnov formula for scalar products of Bethe states, we eliminate the awkward dependence on the system size and distinguish the Bethe states into convenient sectors. These steps allow us to express the rather impractical sum over Bethe states as a multiple rapidity integral in various alternative forms. Moreover, we examine the singularities of the obtained integrand and calculate the contribution of the multivariable kinematical poles, which is essential information for the derivation of the asymptotics of interest.

arXiv:2010.03553Oct 7, 20205 citations

On computing non-equilibrium dynamics following a quench

Albertus J. J. M. de Klerk

Computing the non-equilibrium dynamics that follows a quantum quench is difficult, even in exactly solvable models. Results are often predicated on the ability to compute overlaps between the initial state and eigenstates of the Hamiltonian that governs time evolution. Except for a handful of known cases, it is generically not possible to find these overlaps analytically. Here we develop a numerical approach to preferentially generate the states with high overlaps for a quantum quench starting from the ground state or an excited state of an initial Hamiltonian. We use these preferentially generated states, in combination with a "high overlap states truncation scheme" and a modification of the numerical renormalization group, to compute non-equilibrium dynamics following a quench in the Lieb-Liniger model. The method is non-perturbative, works for reasonable numbers of particles, and applies to both continuum and lattice systems. It can also be easily extended to more complicated scenarios, including those with integrability breaking.

Dec 1, 202120 citations