The exactly solvable Lieb-Liniger model of interacting bosons in one dimension has attracted renewed interest as current experiments with ultracold 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 20 particles for attractive interactions. We discuss the features of the solutions, and how they deviate from the well-known string solutions [Thacker, Rev. Mod. Phys. 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 one-dimensional Bose gas on a ring.
The exactly solvable Lieb-Liniger model of interacting bosons in one dimension has attracted renewed interest as current experiments with ultracold 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 20 particles for attractive interactions. We discuss the features of the solutions, and how they deviate from the well-known string solutions [Thacker, Rev. Mod. Phys. 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 one-dimensional Bose gas on a ring.
The exact ground state of the many-body Schr\"odinger equation for N bosons on a one-dimensional ring interacting via a pairwise $\ensuremath{\delta}$-function interaction is presented for up to 50 particles. The solutions are obtained by solving Lieb and Liniger's system of coupled transcendental equations numerically for finite N. The ground-state energies for repulsive and attractive interactions are shown to be smoothly connected at the point of zero interaction strength, implying that the Bethe ansatz can be used also for attractive interactions for all cases studied. For repulsive interactions the exact energies are compared to (i) Lieb and Liniger's thermodynamic limit solution and (ii) the Tonks-Girardeau gas limit. It is found that the energy of the thermodynamic limit solution can differ substantially from that of the exact solution for finite N when the interaction is weak or when N is small. A simple relation between the Tonks-Girardeau gas limit and the solution for finite interaction strength is revealed. For attractive interactions we find that the true ground-state energy is given to a good approximation by the energy of the system of N attractive bosons on an infinite line, provided the interaction is stronger than the critical interaction strength of mean-field theory.
We show that quantum solitons in the Lieb-Liniger Hamiltonian are precisely the yrast states. We identify such solutions with Lieb's type II excitations from weak to strong interactions, clarifying a long-standing question of the physical meaning of this excitation branch. We demonstrate that the metastable quantum phase transition previously found in mean-field analysis of the weakly interacting Lieb-Liniger Hamiltonian [Phys. Rev. A 79, 063616 (2009)] extends into the medium- to strongly interacting regime of a periodic one-dimensional Bose gas. Our methods are exact diagonalization, finite-size Bethe ansatz, and the boson-fermion mapping in the Tonks-Girardeau limit.
Abstract The coordinate Bethe ansatz can be extended to a model, the Lieb–Liniger model, of a one-dimensional gas of Bosons interacting with repulsive δ-function potentials. It has attracted attention due to its relevance for experimental developments in the fields of ultracold gases and optical lattices. This chapter provides an exposition of the related classical nonlinear Schrödinger equation, followed by its generalization to the quantum model. It explores a limiting case, the Tonks-Girardeau gas. The δ-function potentials supply a kind of boundary condition on the wave functions allowing us to analyze the eigenfunctions of the Bethe ansatz, which are examined on the infinite line and for periodic boundary conditions. The latter leads to the Bethe ansatz equations. The solution of these equations is achieved in the thermodynamic limit for the ground state and for low-lying excited states.
We study the ground-state properties and excitation spectrum of the Lieb-Liniger model, i.e. the one-dimensional Bose gas with repulsive contact interactions. We solve the Bethe-Ansatz equations in the thermodynamic limit by using an analytic method based on a series expansion on orthogonal polynomials developed in and push the expansion to an unprecedented order. By a careful analysis of the mathematical structure of the series expansion, we make a conjecture for the analytic exact result at zero temperature and show that the partially resummed expressions thereby obtained compete with accurate numerical calculations. This allows us to evaluate the density of quasi-momenta, the ground-state energy, the local two-body correlation function and Tan’s contact. Then, we study the two branches of the excitation spectrum. Using a general analysis of their properties and symmetries, we obtain novel analytical expressions at arbitrary interaction strength which are found to be extremely accurate in a wide range of intermediate to strong interactions.
In this thesis we study the physics of quantum many-body systems confined to one-dimensional geometries. The work was motivated by the recent success of experimentalists in developing atom traps, capable of restricting the motion of the individual atoms to a single spatial dimension. Specifically, we look at aspects of the one-dimensional Bose gas including; excitation spectrum, correlation functions, and dynamical behaviour. In Chapter \ref{ch:excitation1D} we consider the Lieb-Liniger model of interacting bosons in one-dimension. 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 novel features of the solutions, including deviations from the well-known string solutions due to finite size effects. 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. Our results are compared to those obtained via exact diagonalization of the Hamiltonian in a truncated basis. In Chapter \ref{ch:g2} we analytically calculate the spatial nonlocal pair correlation function for an interacting uniform one dimensional Bose gas at finite temperature and propose an experimental method to measure nonlocal correlations. Our results span six different physical realms, including the weakly and strongly interacting regimes. We show explicitly that the characteristic correlation lengths are given by one of four length scales: the thermal de Broglie wavelength, the mean interparticle separation, the healing length, or the phase coherence length. In all regimes, we identify the profound role of interactions and find that under certain conditions the pair correlation may develop a global maximum at a finite interparticle separation due to the competition between repulsive interactions and thermal effects. In Chapter \ref{ch:casimirdrag} we study the drag force below the critical velocity for obstacles moving in a superfluid. The absence of drag is well established in the context of the mean-field Gross-Pitaevskii theory. We calculate the next order correction due to quantum and thermal fluctuations and find a non-zero force acting on a delta-function impurity moving through a quasi-one-dimensional Bose-Einstein condensate at all subcritical velocities and at all temperatures. The force occurs due to an imbalance in the Doppler shifts of reflected quantum fluctuations from either side of the impurity. Our calculation is based on a consistent extension of Bogoliubov theory to second order in the interaction strength, and finds new analytic solutions to the Bogoliubov-de Gennes equations for a gray soliton. In Chapter \ref{ch:solitons} we study the effect of quantum noise on the stability of a soliton. We find the soliton solutions exactly define the reflectionless potentials of the Bogoliubov-de Gennes equations. This results in complete stability of the solitons in a purely one dimensional system. We look at the modifications to the density profile of a black soliton due to quantum fluctuations.
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.