In this thesis we investigate the properties of ultra-cold quantum gases in reduced dimension and the effects of harmonic confinement on soliton-like properties. We study regimes of agreement between mean-field and many-body theories the generation of entanglement between initially independent finite sized atomic systems. Classical solitons are non-dispersing waves which occur in integrable systems, such as atomic Bose-Einstein condensates in one dimension. Bright and dark solitons are possible, which exist as peaks or dips in density. Quantum solitons are the bound-state solutions to a system satisfying quantum integrability, given via the Bethe Ansatz. Such integrability is broken by the introduction of harmonic confinement. We investigate the equivalence of the classical field and many-body solutions in the limit of large numbers of atoms and derive numerical and variational approaches to examine the ground state energy in harmonic confinement and the fidelity between a Hartree-product solution and a quantum soliton solution. Soliton collisions produce no entanglement between either state and result only in an asymptotic position and phase shift, however external potentials break integrability and thus give the possibility of entangling solitons. We investigate the dynamical entanglement generation between two atomic dimers in harmonic confinement via exact diagonalisation in a basis of Harmonic oscillator functions, making use of the separability of the centre-of-mass component of the Hamiltonian. We show repulsive states show complex dynamics, but with an overall tendency towards states of larger invariant correlation entropy, whereas attractive states resist entanglement unless a phase matching condition is satisfied. This phase matching condition could in theory be used to generate states with highly non-Poissonian number superpositions in atomic systems with controlled number.
In this thesis we investigate the properties of ultra-cold quantum gases in reduced dimension and the effects of harmonic confinement on soliton-like properties. We study regimes of agreement between mean-field and many-body theories the generation of entanglement between initially independent finite sized atomic systems. Classical solitons are non-dispersing waves which occur in integrable systems, such as atomic Bose-Einstein condensates in one dimension. Bright and dark solitons are possible, which exist as peaks or dips in density. Quantum solitons are the bound-state solutions to a system satisfying quantum integrability, given via the Bethe Ansatz. Such integrability is broken by the introduction of harmonic confinement. We investigate the equivalence of the classical field and many-body solutions in the limit of large numbers of atoms and derive numerical and variational approaches to examine the ground state energy in harmonic confinement and the fidelity between a Hartree-product solution and a quantum soliton solution. Soliton collisions produce no entanglement between either state and result only in an asymptotic position and phase shift, however external potentials break integrability and thus give the possibility of entangling solitons. We investigate the dynamical entanglement generation between two atomic dimers in harmonic confinement via exact diagonalisation in a basis of Harmonic oscillator functions, making use of the separability of the centre-of-mass component of the Hamiltonian. We show repulsive states show complex dynamics, but with an overall tendency towards states of larger invariant correlation entropy, whereas attractive states resist entanglement unless a phase matching condition is satisfied. This phase matching condition could in theory be used to generate states with highly non-Poissonian number superpositions in atomic systems with controlled number.
This paper investigates bright quantum-matter-wave solitons beyond the Gross-Pitaevskii equation (GPE). As proposals for interferometry and creating nonlocal quantum superpositions have been formed, it has become necessary to investigate effects not present in mean-field models. We investigate the effect of harmonic confinement on the internal degrees of freedom, as the ratio of zero-point harmonic oscillator length to classical soliton length, for different numbers of atoms. We derive a first-order energy correction for the addition of a harmonic potential to the many-body wave function and use this to create a variational technique based on energy minimization of this wave function for an arbitrary number of atoms, and include numerics based on diagonalization of the Hamiltonian in a basis of harmonic oscillator Fock states. Finally we compare agreement between a Hartree product ground state and the Bethe ansatz solution with a Gaussian envelope localizing the center of mass and show a region of good agreement.
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.
We investigate finite number effects in collisions between two states of an initially well defined number of identical bosons with attractive contact interactions, oscillating in the presence of harmonic confinement in one dimension. We investigate two N/2 atom bound states, which are initially displaced (symmetrically) from the trap center, and then left to freely evolve. For sufficiently attractive interactions, these bound states are like those found through use of the Bethe ansatz (quantum solitons). However, unlike the free case, the integrability is lost due to confinement, and collisions can cause mixing into different bound-state configurations. We study the system numerically for the simplest case of N=4, via an exact diagonalization of the Hamiltonian within a finite basis, investigating left-right number uncertainty as our primary measure of entanglement. We find that for certain interaction strengths, a phase-matching condition leads to resonant transfer to different bound-state configurations with highly non-Poissonian relative number statistics.
We investigate finite-number effects in collisions between two states of an initially well-known number of identical bosons with contact interactions, oscillating in the presence of harmonic confinement in one dimension. We investigate two N/2 (interacting) ground states, which are initially displaced from the trap center, and the effects of varying interaction strength. The numerics focus on the simplest case of N=4. In the noninteracting case, such a system would display periodic oscillation with a half harmonic oscillator period (due to the left-right symmetry). With the addition of contact interactions between the bosons, collisions generate entanglement between each of the states and distribute energy into other modes of the oscillator. We study the system numerically via an exact diagonalization of the Hamiltonian with a finite basis, investigating left-right number uncertainty as our primary measure of entanglement. Additionally, we study the time evolution and equilibration of the single-body von Neumann entropy for both the attractive and repulsive cases. We identify parameter regimes for which attractive interactions create behavior qualitatively different from that of repulsive interactions, due to the presence of bound states (quantum solitons), and explain the processes behind this.
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.
This thesis examines two aspects of the non-equilibrium dynamics of Bose-Einstein condensates (BECs). Our particular interests are the relaxation of near-integrable systems, and the inter-particle correlations present in BECs, both in equilibrium and following an abrupt change in the Hamiltonian. In the first part of this thesis we investigate the stability of solitary wave solutions to a system of two coupled one-dimensional non-linear Schrodinger equations (CNLSE). This system has an integrable point known as the Manakov[1] model which has analytic dark-bright soliton solutions[2]. We break the integrability of this model by varying the inter-species interaction strength, and employ variational and numerical methods in order to find dark-bright solitary wave solutions, away from the integrable point. The time evolution of these states is calculated via numerical integration of the CNLSE, which allows us to assess the stability of these solutions in isolation and to quantify their robustness against collisions with a dark soliton. We find that there is a broad region of the parameter space in which solutions for black-bright (stationary) solitary waves can be found. In this domain, the integrability breaking is revealed only during collisions. Prior to a collision, the numerical solutions are indistinguishable from true solitons, however, the collisional stability of these solitary waves is significantly affected by the extent to which integrability is broken. We find that there is a smooth transition between the non-dispersive particle-like nature of soliton interactions for integrable systems, and the destructive collisions which take place when the system is perturbed far from the integrable point. In the case of moving (grey-bright) solitary waves, we find a more restricted region of the parameter space which admits solitary wave solutions. In this domain we are able to find approximately stable solitary wave solutions using a variational ansatz. The collisional stability of these grey-bright solitary waves is also found to depend upon the extent to which integrability is broken. We conclude that the observation of stable, long-lived dark-bright solitary waves is not sufficient to indicate near-integrability of the two-component system. However, near-integrability may inferred if dark-bright solitary waves are observed to survive multiple collisions. These results are not only of theoretical interest, but may also inform future experiments which explore the behaviour of solitary waves in two-component Bose-Einstein condensates. Following this we investigate the properties of magnetic solitary wave solutions to the CNLSE. We find that the analytic ‘magnetic soliton’ solutions derived in reference [3] are dynamically unstable and relax to non-stationary magnetic solitary wave states over a short timescale. We investigate the collisional dynamics of these magnetic solitary waves and find that these excitations are remarkably robust. In the second part of this thesis we calculate the fluctuations in the number difference between two sub-regions of a Bose-Einstein condensate. This study is motivated by recent experiments performed by the Truscott group at the Australian National University (ANU), which indicate sub- Poissonian fluctuations even for condensates with large thermal depletion. This result cannot easily be reconciled with existing literature, which predicts strongly super-Poissonian fluctuations below the critical temperature (Tc). We develop one-dimensional, and three-dimensional, models to describe these experiments and employ the Bogoliubov formalism to explore the quantum fluctuations of the condensate, which have not previously been studied in detail. We calculate the relative number fluctuations between two halves of the condensate, and also between multiple pairs of spatial bins. For harmonically confined condensates at zero temperature, quantum correlations result in a suppression of the relative number fluctuations between the two halves of the condensate below the shot-noise level. However, the presence of even a small thermal component results in strongly super-Poissonian fluctuations below Tc; this is in agreement with the existing theoretical and experimental results. In addition, we find that the fluctuations of the harmonically confined condensate are non-uniform, with the relative number fluctuations being peaked towards the edges of the condensate. The recent studies of the Truscott group measure the fluctuations after free expansion. We therefore calculate the propagation of correlations following a sudden removal of the harmonic confinement. We find that free expansion results in a strong decrease in the relative number fluctuations. Condensates which exhibit super-Poissonian fluctuations when harmonically confined may therefore display number- squeezing after expansion. A second set of experiments by the Truscott group measured the density fluctuations present in small ‘atom laser’ pulses of atoms out-coupled from their BEC. We also model this procedure and find that immediately after the out-coupling process, the number fluctuations of the out-coupled atomic pulse are approximately Poissonian, irrespective of the correlations present in the atomic reservoir. However, the competing effects of interactions between the atomic pulse and the condensate reservoir, and free expansion of the pulse after out-coupling, can result in the measured fluctuations being either number squeezed or super-Poissonian. Our results will guide future experiments planned by the Truscott group.