We numerically study free expansion of a few Lieb-Liniger bosons, which are initially in the ground state of an infinitely deep hard-wall trap. Numerical calculation is carried out by employing a standard Fourier transform, as follows from the Fermi-Bose transformation for a time-dependent Lieb-Liniger gas. We study the evolution of the momentum distribution, the real-space single-particle density, and the occupancies of natural orbitals. Our numerical calculation allows us to explore the behavior of these observables in the transient regime of the expansion, where they are nontrivially affected by the particle interactions. We derive analytically (by using the stationary phase approximation) the formula which connects the asymptotic shape of the momentum distribution and the initial state. For sufficiently large times, the momentum distribution coincides (up to a simple scaling transformation) with the shape of the real-space single-particle density (the expansion is asymptotically ballistic). Our analytical and numerical results are in good agreement.
Exact solutions of the Schrödinger equation describing a freely expanding Lieb-Liniger gas of delta-interacting bosons in one spatial dimension are constructed. We demonstrate that for any interaction strength the system enters a strongly correlated regime during such expansion. The asymptotic form of the wave function is shown to have the form characteristic for "impenetrable-core" bosons. Exact solutions are obtained by transforming a fully antisymmetric (fermionic) time-dependent wave function that obeys the Schrödinger equation for a free gas. This transformation employs a differential Fermi-Bose mapping operator that depends on the strength of the interaction and the number of particles.
We use Gaudin's Fermi-Bose mapping operator to calculate exact solutions for the Lieb-Liniger model in a linear (constant-force) potential (the constructed exact stationary solutions are referred to as the Lieb-Liniger-Airy wave functions). The ground-state properties of the gas in the wedgelike trapping potential are calculated in the strongly interacting regime by using Girardeau's Fermi-Bose mapping and the pseudopotential approach in the
1/c approximation (
c denotes the strength of the interaction). We point out that quantum dynamics of Lieb-Liniger wave packets in the linear potential can be calculated by employing an
N-dimensional Fourier transform as in the case of free expansion.
One-dimensional Bose gases are considered, interacting either through the hard-core potentials or through the contact delta potentials. Interest in these gases gained momentum because of the recent experimental realization of quasi-one-dimensional Bose gases in traps with tightly confined radial motion, achieving the Tonks-Girardeau (TG) regime of strongly interacting atoms. For such gases the Fermi-Bose mapping of wavefunctions is applicable. The aim of the present communication is to give a brief survey of the problem and to demonstrate the generality of this mapping by emphasizing that: (i) It is valid for nonequilibrium wavefunctions, described by the time-dependent Schrödinger equation, not merely for stationary wavefunctions. (ii) It gives the whole spectrum of all excited states, not merely the ground state. (iii) It applies to the Lieb-Liniger gas with the contact interaction, not merely to the TG gas of impenetrable bosons.
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.
Recent theoretical and experimental results demonstrate a close connection between the super-Tonks-Girardeau (STG) gas and a one-dimensional (1D) hard-sphere Bose (HSB) gas with hard-sphere diameter nearly equal to the 1D scattering length
a1D of the STG gas, a highly excited gaslike state with nodes only at interparticle separations $|{x}_{j\ensuremath{\ell}}|={x}_{\mathrm{node}}\ensuremath{\approx}{a}_{1\mathrm{D}}$. It is shown herein that when the coupling constant
gB in the Lieb-Liniger interaction ${g}_{B}\ensuremath{\delta}({x}_{j\ensuremath{\ell}})$ is negative and $|{x}_{12}|\ensuremath{\geqslant}{x}_{\mathrm{node}}$, the STG and HSB wave functions for
N=2 particles are not merely similar, but identical; the only difference between the STG and HSB wave functions is that the STG wave function allows a small penetration into the region
∣x12∣<xnode, whereas for a HSB gas with hard-sphere diameter
ahs=xnode, the HSB wave function vanishes when
∣x12∣<ahs. Arguments are given suggesting that the same theorem holds also for
N>2. The STG and HSB wave functions for
N=2 are given exactly in terms of a parabolic cylinder function, and for $N\ensuremath{\geqslant}2$,
xnode is given accurately by a simple parabola. The metastability of the STG phase generated by a sudden change of the coupling constant from large positive to large negative values is explained in terms of the very small overlap between the ground state of the Tonks-Girardeau gas and collapsed cluster states.