Exact methods in analysis of nonequilibrium dynamics of integrable models: application to the study of correlation functions in nonequilibrium 1D Bose gas

In this paper we study nonequilibrium dynamics of one dimensional Bose gas from the general perspective of dynamics of integrable systems. After outlining and critically reviewing methods based on inverse scattering transform, intertwining operators, q-deformed objects, and extended dynamical conformal symmetry, we focus on the form-factor based approach. Motivated by possible applications in nonlinear quantum optics and experiments with ultracold atoms, we concentrate on the regime of strong repulsive interactions. We consider dynamical evolution starting from two initial states: a condensate of particles in a state with zero momentum and a condensate of particles in a gaussian wavepacket in real space. Combining the form-factor approach with the method of intertwining operator we develop a numerical procedure which allows explicit summation over intermediate states and analysis of the time evolution of non-local density-density correlation functions. In both cases we observe a tendency toward formation of crystal-like correlations at intermediate time scales.

💡 Research Summary

The paper presents a comprehensive study of the nonequilibrium dynamics of a one‑dimensional Bose gas (the Lieb‑Liniger model) from the perspective of integrable‑system methods. After a concise introduction that motivates the problem through recent advances in ultracold‑atom experiments and nonlinear quantum optics, the authors review four major analytical frameworks that have been employed for integrable models: (i) the inverse scattering transform (IST), (ii) intertwining‑operator techniques, (iii) q‑deformed algebraic constructions, and (iv) extended dynamical conformal symmetry. For each approach the strengths and limitations are discussed, with particular emphasis on the difficulty of handling generic, highly excited initial states within the IST formalism.

The core of the work is the development and application of a form‑factor based method combined with an intertwining‑operator mapping. Form factors are exact matrix elements of local operators (here the density operator) between Bethe‑Ansatz eigenstates. By using an intertwining operator the authors map a simple, non‑interacting basis onto the interacting Bethe basis, thereby expressing any chosen initial wavefunction as a linear combination of exact eigenstates. This mapping is crucial because it allows the subsequent summation over intermediate states to be performed with fully known energies and form‑factor amplitudes.

Two distinct initial conditions are examined, both in the regime of strong repulsive interactions (the Tonks‑Girardeau limit). The first is a zero‑momentum condensate, i.e. all particles occupying the (k=0) mode. The second is a Gaussian wavepacket prepared in real space, which mimics the experimental situation of a localized cloud released into a one‑dimensional waveguide. For each case the authors construct a truncated Bethe‑Ansatz spectrum by enumerating all quantum‑number configurations within a finite momentum window. They then compute the corresponding energies, rapidities, and the density‑density form factors using known determinant representations. The time‑dependent non‑local density‑density correlation function \