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

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.