No Arabic abstract
The exactly solvable Lieb-Liniger model of interacting bosons in one-dimension has attracted renewed interest as current experiments with ultra-cold atoms begin to probe this regime. Here we numerically solve the equations arising from the Bethe ansatz solution for the exact many-body wave function in a finite-size system of up to twenty particles for attractive interactions. We discuss the novel features of the solutions, and how they deviate from the well-known string solutions [H. B. Thacker, Rev. Mod. Phys. textbf{53}, 253 (1981)] at finite densities. We present excited state string solutions in the limit of strong interactions and discuss their physical interpretation, as well as the characteristics of the quantum phase transition that occurs as a function of interaction strength in the mean-field limit. Finally we compare our results to those of exact diagonalization of the many-body Hamiltonian in a truncated basis. We also present excited state solutions and the excitation spectrum for the repulsive 1D Bose gas on a ring.
We propose a variational approximation to the ground state energy of a one-dimensional gas of interacting bosons on the continuum based on the Bethe Ansatz ground state wavefunction of the Lieb-Liniger model. We apply our variational approximation to a gas of dipolar bosons in the single mode approximation and obtain its ground state energy per unit length. This allows for the calculation of the Tomonaga-Luttinger exponent as a function of density and the determination of the structure factor at small momenta. Moreover, in the case of attractive dipolar interaction, an instability is predicted at a critical density, which could be accessed in lanthanide atoms.
In this work we perform an ab-initio study of an ideal two-dimensional sample of 4He atoms, a model for 4He films adsorbed on several kinds of substrates. Starting from a realistic hamiltonian we face the microscopic study of the excitation phonon-roton spectrum of the system at zero temperature. Our approach relies on Path Integral Ground State Monte Carlo projection methods, allowing to evaluate exactly the dynamical density correlation functions in imaginary time, and this gives access to the dynamical structure factor of the system S(q,omega), containing information about the excitation spectrum E(q), resulting in sharp peaks in S(q,omega). The actual evaluation of S(q,omega) requires the inversion of the Laplace transform in ill-posed conditions, which we face via the Genetic Inversion via Falsification of Theories technique. We explore the full density range from the region of spinodal decomposition to the freezing density, i.e. 0.0321 A^-2 - 0.0658 A^-2. In particular we follow the density dependence of the excitation spectrum, focusing on the low wave--vector behavior of E(q), the roton dispersion, the strength of single quasi--particle peak, Z(q), and the static density response function, chi(q). As the density increases, the dispersion E(q) at low wave--vector changes from a super-linear (anomalous dispersion) trend to a sub-linear (normal dispersion) one, anticipating the crystallization of the system; at the same time the maxon-roton structure, which is barely visible at low density, becomes well developed at high densities and the roton wave vector has a strong density dependence. Connection is made with recent inelastic neutron scattering results from highly ordered silica nanopores partially filled with 4He.
We use the coordinate Bethe ansatz to exactly calculate matrix elements between eigenstates of the Lieb-Liniger model of one-dimensional bosons interacting via a two-body delta-potential. We investigate the static correlation functions of the zero-temperature ground state and their dependence on interaction strength, and analyze the effects of system size in the crossover from few-body to mesoscopic regimes for up to seven particles. We also obtain time-dependent nonequilibrium correlation functions for five particles following quenches of the interaction strength from two distinct initial states. One quench is from the non-interacting ground state and the other from a correlated ground state near the strongly interacting Tonks-Girardeau regime. The final interaction strength and conserved energy are chosen to be the same for both quenches. The integrability of the model highly constrains its dynamics, and we demonstrate that the time-averaged correlation functions following quenches from these two distinct initial conditions are both nonthermal and moreover distinct from one another.
Based on the inhomogeneous T-Q relation constructed via the off-diagonal Bethe Ansatz, a systematic method for retrieving the Bethe-type eigenstates of integrable models without obvious reference state is developed by employing certain orthogonal basis of the Hilbert space. With the XXZ spin torus model and the open XXX spin-1/2 chain as examples, we show that for a given inhomogeneous T-Q relation and the associated Bethe Ansatz equations, the constructed Bethe-type eigenstate has a well-defined homogeneous limit.
We use the coordinate Bethe ansatz to study the Lieb-Liniger model of a one-dimensional gas of bosons on a finite-sized ring interacting via an attractive delta-function potential. We calculate zero-temperature correlation functions for seven particles in the vicinity of the crossover to a localized solitonic state and study the dynamics of a system of four particles quenched to attractive interactions from the ideal-gas ground state. We determine the time evolution of correlation functions, as well as their temporal averages, and discuss the role of bound states in shaping the postquench correlations and relaxation dynamics.