Within the past five years Albert Einstein's concept of a dilute atomic Bose Condensate has been realized in many experimental laboratories. Temperatures in the nano-Kelvin regime have been achieved using magnetic and optical trapping of laser and evaporatively cooled atoms. At such temperatures the relative de Broglie wavelengths of the gaseous trapped atoms can become long compared to their mean spacing, and through a process of bosonic amplification a "quantum phase transition" takes place involving 10 3 to 10 11 atoms, most of which end up in identical single particle quantum states, whose length scales are determined by the external trap. Rb, Na, Li and atomic H have been trapped at variable densities of the order of 10 13 /cm 3 , and in traps of varying geometries. Of these all but Li have an effective repulsive atomic pair interaction, but utilization of molecular Feshbach resonances allows the interactions of other species to be tuned over wide ranges of strengths, including control of the sign of the effective atomic pair interactions. Fully quantum and macroscopic systems at such low densities are a theorist's dream: simple Hartree type mean field theory provides a startlingly accurate description of density profiles, low energy excitation frequencies, and such a description, commonly called the Non-linear Schrödinger Equation (NLSE) or the Gross-Pitaevskii (GP) Equation, will be explored here. The NLSE appropriate for attractive atomic interactions is known to lead to chaotically unstable dynamics, and eventual implosion of the condensate should the local number density exceed a critical value. In this work we illustrate this type of chaotic collapse for attractive condensates, and then explore the types of chaotic dynamics of solitons and vortices, which are the signatures of dynamical non-linearity in the repulsive case. Finally the implications of the symmetry-breaking associated with phase rigidity are explored in model simulations of repulsive condensates: condensates with repulsive atomic interactions break into phase domains when subject to weak shocks and, perhaps surprisingly, break into chaotic "laser-speckle" type patterns as the shock level increases. The fully quantum mechanical NLSE thus displays a full range of chaotic types of motion: from particle-like chaotic collisions of solitons and vortices to fully developed time dependent wave chaos.
Within the past five years Albert Einstein's concept of a dilute atomic Bose Condensate has been realized in many experimental laboratories. Temperatures in the nano-Kelvin regime have been achieved using magnetic and optical trapping of laser and evaporatively cooled atoms. At such temperatures the relative de Broglie wavelengths of the gaseous trapped atoms can become long compared to their mean spacing, and through a process of bosonic amplification a "quantum phase transition" takes place involving 10 3 to 10 11 atoms, most of which end up in identical single particle quantum states, whose length scales are determined by the external trap. Rb, Na, Li and atomic H have been trapped at variable densities of the order of 10 13 /cm 3 , and in traps of varying geometries. Of these all but Li have an effective repulsive atomic pair interaction, but utilization of molecular Feshbach resonances allows the interactions of other species to be tuned over wide ranges of strengths, including control of the sign of the effective atomic pair interactions. Fully quantum and macroscopic systems at such low densities are a theorist's dream: simple Hartree type mean field theory provides a startlingly accurate description of density profiles, low energy excitation frequencies, and such a description, commonly called the Non-linear Schrödinger Equation (NLSE) or the Gross-Pitaevskii (GP) Equation, will be explored here. The NLSE appropriate for attractive atomic interactions is known to lead to chaotically unstable dynamics, and eventual implosion of the condensate should the local number density exceed a critical value. In this work we illustrate this type of chaotic collapse for attractive condensates, and then explore the types of chaotic dynamics of solitons and vortices, which are the signatures of dynamical non-linearity in the repulsive case. Finally the implications of the symmetry-breaking associated with phase rigidity are explored in model simulations of repulsive condensates: condensates with repulsive atomic interactions break into phase domains when subject to weak shocks and, perhaps surprisingly, break into chaotic "laser-speckle" type patterns as the shock level increases. The fully quantum mechanical NLSE thus displays a full range of chaotic types of motion: from particle-like chaotic collisions of solitons and vortices to fully developed time dependent wave chaos.
Bose-Einstein condensate (BEC) is a new state of matter, experimentally realized in a dilute alkali gas of weakly interacting bosons. The gas is cooled to very low temperatures, of the order of nano-Kelvin, at which point a large fraction of the atoms occupy the lowest quantum state. Under such conditions, the quantum effects become apparent on a macroscopic scale. Many fundamental ideas
of quantum mechanics have been realized and proven to be true in experiments,
by studying the Bose-Einstein condensates. Amongst them, the observation of
macroscopic matter wave interference, the creation of vortices, dark and bright
solitons etc, find special mention. The dynamics of this system in presence of an
optical lattice reveals a superfluid to Mott insulator quantum phase transition.
This has opened a new avenue for research on BEC and its applications. From
a theoretical standpoint, the dynamics of BEC can be described by means of an
effective mean-field theory. This formalism is much simpler than the many-body
Schr¨odinger equation and based on a classical nonlinear evolution equation, known
as Gross-Pitaevskii equation. The nonlinear nature of the GPE leads to several interesting phenomena like soliton propagation, vortices, four wave mixing etc.,
which are closely related to nonlinear fiber optics.
The present thesis deals with nonlinear excitations and different phases of BEC
and is organized as follows: after a brief introduction, in chapter two, we start
with BEC in a one-dimensional shallow optical lattice, experiencing both cubic
and quintic nonlinearity. A superfluid density wave, of twice the periodicity of
lattice potential, is identified. This indicates matter redistribution, where atoms
from neighbouring sites have been depleted leading to a density wave behaviour.
Interestingly, it ceases to exist when only one of these interactions is operative.
We predict the loss of superfluidity through a classical dynamical phase transition,
where modulational instability leads to the loss of phase coherence. In presence
of pure two-body interaction, the superfluid density shows a periodicity commensurate
with the lattice potential. Apart from the superfluid state, a density wave
insulating phase is also found to exist, possessing two frequency modulations commensurate
with the lattice potential. A dynamical phase transition connects the
superfluid phase with an insulating phase, which in turn, is connected with the
density wave insulating phase. The solutions are found to be marginally stable, as
per the Vakhitov-Kolokolov (VK) criterion.
Chapter three describes the dynamics of the atomic-molecular Bose-Einstein
condensate (AMBEC) by mean-field Gross-Pitaevskii equation (GPE), taking into
account the atom-molecular interaction, as well as atom-molecule inter-conversion.
The two-body atom-atom, atom-molecule and molecule-molecule interactions are
analogous to the cubic nonlinearity of the standard GPE, whereas, the quadratic
nonlinear term arising from the inter-conversion between atoms and molecules is
identical to the scenario in strongly coupled BEC in one dimension. The presence
of both weak and strong interactions adds further richness to this system. We
find the molecular BEC possesses a complex dark (grey) soliton, while the atomic condensate is a bright one. The energy and momentum of this system is computed
in order to see the effect of the both weak and strong interactions on Lieb mode
profile and compare it with the pure atomic condensate profile. Another class
of solutions exist in this system, where both atomic and molecular condensates
possess bright solitons.
In chapter four, we investigate the effect of a harmonic trap on the dynamics
of grey solitons in Bose-Einstein condensates. It is observed that, these nonlinear
excitations of the cigar-shaped Bose-Einstein condensate have strong coupling
with the trap at low momenta and hence can be effectively isolated from the
Bogoliubov sound modes, which respond weakly to harmonic confinement. This
strong coupling with the trap also makes the grey soliton amenable for control and
manipulation through trap modulation and temporal variation of the two-body interaction.
We have shown that in presence of an expulsive oscillator trap, the grey
soliton gets amplified, while having a compression in its width. Exact solutions are
also obtained in case of the two component Bose-Einstein condensates (TBECs),
where one component is a grey soliton and the other a bright one. The presence
of the trap and inter-species interaction significantly affects the grey soliton dynamics.
Similar to the single component, in the present case, we have shown that
both grey and bright soliton can be compressed simultaneously. A detailed study
of the energy and momentum of the soliton is reported, in order to find the effect
of the trap on the Lieb mode.
In chapter five, we study a BEC-defect complex in a quasi-one dimensional
scenario, where the defect atom moves in the condensate with constant velocity.
Explicit solutions of the coupled mean-field equations, describing defect-grey soliton
dynamics are obtained, demonstrating the coexistence of grey soliton and a
localized defect. Unlike the case of dark soliton, where the defect trapping center
has vanishing superfluid density, the moving grey soliton necessarily possesses a finite superfluid component at the defect location. The presence of the defect can
lower the grey soliton’s maximum velocity, as compared to the defect free case,
where the maximum velocity is the sound velocity. We have also carried out the
stability analysis of the obtained solutions using renormalized momentum formalism
and Vakhitov-Kolokolov criteria and find that the obtained solutions are stable
through both the methods. The grey soliton’s energy gets substantially modified
through its interaction with the defect, opening up the possibility of its control
through defect dynamics.
Chapter six describes a parametrically forced Bose-Einstein condensate in the
presence of a general moving optical lattice. The interaction between the atoms in
the condensate and the time dependent lattice potential leads to a novel propagating
superfluid matter wave, which can be controlled through chirp management.
This system, when placed in a trap, accelerates and undergoes rapid nonlinear
compression, controlled by the chirp. The density achieves its maximum, when
the matter wave changes direction. It is also shown that this system exhibits dynamical
superfluid insulator transition, where superfluidity breaks down and the
condensate transits to an insulating phase. The stability has been investigated
using Vakhitov-Kolokolov criterion and found that it depends on the chirp, as well
as on the nature of the atom-atom interaction. The exact expression for energy
is obtained and analyzed in detail to gain physical understanding of the chirp
management of the sinusoidal excitations and also the dynamical phase transition.
Bose–Einstein Condensation (BEC) of atomic gases has attracted a renewed theoretical and experimental interest in quantum many-body systems at extremely low temperatures (Pethick & Smith; 2002). This excitement stems from two favorable features: (1) by applying magnetic fields and lasers, most of the system parameters, such as the shape, dimensionality, internal states of the condensates, and even the strength of the interatomic interactions, are controllable; (2) due to the diluteness, the mean-field theory explains experiments quite well. In particular, the Gross–Pitaevskii (GP) equation demonstrates its validity as a basic equation for the condensate dynamics. The GP equation is the counterpart of the nonlinear Schrodinger (NLS) equation in nonlinear optics. Thus, a study based on nonlinear analysis is possible and important. In nonlinear physics, a soliton is remarkable object not only for the fact that exact solutions can be obtained but also for its usefulness as a communications tool due to its robustness. In general, solitons are formed under the balance between nonlinearity and dispersion. For atomic condensates, the former is attributed to the interatomic interactions, while the latter comes from the kinetic energy. Either dark or bright solitons are allowable depending on the positive or negative sign of the interatomic coupling constants g, respectively, and indeed have been observed in a quasione-dimensional (q1D) optically constructed waveguide (Strecker et al.; 2002) (Khaykovich et al.; 2002). Such matter-wave solitons are expected in atom optics for applications in atom laser, atom interferometry, and coherent atom transport (Meystre; 2002). In this chapter, we extend the analysis of the matter-wave solitons to a multicomponent case by considering the so-called spinor condensate (Stenger et al.; 1998) whose spin degrees of freedom are liberated under optical traps. Based on theoretical and experimental results, we introduce a new integrable model which describes the dynamical properties of the matter-wave soliton of spinor condensates (Ieda et al.; 2004a). We employ the inverse scattering method to solve this model exactly. As a result, we predict the occurrence of undiscovered physical phenomena such as macroscopic spin precession and spin switching. The chapter is organized as follows. In Sec. 2, the mean field theory of condensate is briefly reviewed. Section 3 introduces an effective interatomic coupling in a q1D condensate. Using these results, we consider a spinor condensate in q1D regime in Sec. 4. Then, in Sec. 5, we
Bose-Einstein condensation is a quantum statistical phase transition which was theoretically predicted almost a hundred years ago. After years of seminal research, physicists realized the first almost ideal Bose-Einstein condensates in ultracold dilute atomic gases in 1995. Since then, the theoretical and experimental methods concerning such systems have been developing rapidly, and many fascinating phenomena have been found in these novel quantum systems.\n\nBose-Einstein condensation occurs in a system consisting of massive bosons when a single quantum state becomes macroscopically occupied as the temperature is lowered below the transition temperature. In general, condensates consisting of repulsively interacting bosons exhibit superfluidity: Particle currents can flow in the system without dissipation and viscosity. Moreover, the velocity fields of condensates have to be irrotational, which severely restricts the rotational characteristics of these systems. Apart from the center of mass motion, the system may carry angular momentum in the form of elementary excitations or so-called quantized vortices.\n\nThis Thesis is a theoretical study of subjects related to stability and dynamics of quantized vortices in dilute atomic Bose-Einstein condensates. The precession and instability of off-centered vortices in trapped condensates is investigated both in the zero-temperature limit and at finite temperatures. Dynamical stability of multiply quantized vortices and vortex clusters is studied in axisymmetric trap geometries. Splitting of energetically and dynamically unstable multiply quantized vortices into singly quantized vortices is also studied. Finally, as a separate subject, tunneling of a condensate through a potential barrier is investigated. Majority of this work relies on numerical methods for solving the Gross-Pitaevskii and Bogoliubov equations, which are of central importance in the study of dilute atomic Bose-Einstein condensates.
When a gas of bosonic atoms is cooled sufficiently near to absolute zero, it undergoes a startling transition from a collection of randomly-behaved individual atoms to a phase in which thousands, even millions, of atoms act in collective unison. In this new phase, known as a Bose-Einstein condensate (BEC), the atoms all occupy the exact same quantum mechanical ground state and can be collectively described by a single phase-coherent wavefunction of macroscopic size. But how exactly does this harmonious behaviour arise within a collection of incoherent atoms as their temperature decreases? Does the condensate wavefunction grow from a single central nucleus, or does it emerge simultaneously throughout the system? And how does the situation differ when the gas is cooled suddenly, away from thermodynamic equilibrium? Can non-equilibrium behaviour be inherited by the condensate that forms? In this thesis we aim to answer these and other questions by studying the formation dynamics of Bose-Einstein condensates and the development of phase coherence. The research described herein makes use of several theoretical tools, which we lay out on the workbench in Chapter 2: Semi-classical Hartree-Fock theory provides a mean-field description of a condensate and thermal cloud at equilibrium; Quantum kinetic theory is used to calculate the dynamics of condensate growth into an equilibrium state, from a non-equilibrium thermal cloud; The stochastic projected Gross-Pitaevskii equation (SPGPE) employs a classical-field description to encapsulate the non-equilibrium dynamics of the condensate and low-energy region of the atom cloud, coupled to an equilibrium thermal cloud. This thesis is composed of four principal studies, each constituting a separate chapter.In Chapter 3 we study the formation of a Bose-Einstein condensate in a cigar-shaped 3D harmonic trap, induced by the controlled addition of an attractive ``dimple'' potential along the weak axis in a recent experiment [Garrett et al., Phys. Rev. A 83, 013630 (2011)]. In this manner condensation is induced without cooling due to a localized increase in the phase space density. We perform a quantitative analysis of the thermodynamic transformation in both the sudden and adiabatic regimes for a range of dimple widths and depths. We find good agreement with equilibrium calculations based on self-consistent semi-classical Hartree-Fock theory describing the condensate and thermal cloud. We observe there is an optimal dimple depth that results in a maximum in the condensate fraction. We also study the non-equilibrium dynamics of condensate formation in the sudden turn-on regime, finding good agreement for the observed time dependence of the condensate fraction with calculations based on quantum kinetic theory. In Chapter 4 we study the onset of phase coherence during the formation of a BEC following a sudden quench to below the transition temperature. Using the SPGPE, we simulate a recent experiment performed at ANU in which a sudden RF quench is applied to a dilute gas of metastable Helium atoms confined in an elongated 3D harmonic magnetic trap. At regularly-spaced intervals following the quench, broadband RF pulses are applied to homogeneously outcouple small (<5%) pulses of atoms, which expand and fall to a microchannel plate detector. The second-order correlation function at the detector, postulated to provide a measure of phase coherence of the in-situ cloud from which they were outcoupled, is measured versus time and compared with the condensate growth curve. A delay was observed in the relaxation of the second-order correlation function, relative to the peak in the number of condensate atoms. In our simulations we incorporate the effects of periodic outcoupling both the condensate and thermal cloud, and monitor the time-dependence of first- and second-order correlation functions in both position and momentum space. We observe long-range phase coherence in the final equilibrium state, though mild persistent phase fluctuations are present. We also observe the curious phenomenon of ``superbunching'' - an enhancement of density fluctuations - at specific values of the momentum, which we relate to anomalous correlations of fluctuations. Despite obtaining a good quantitative fit to the experimentally-measured condensate growth curve in our SPGPE simulations, we do not observe any evidence of the delay in coherence observed in the experiment, which remains unexplained. In Chapter 5 we use the SPGPE to study the crossover from a true BEC with long-range phase coherence to a quasicondensate dominated by phase fluctuations, for the case of a Bose gas in an highly-elongated 3D harmonic trap. We use parameters corresponding to an experiment in which a sudden quench resulted in the formation of a non-equilibrium condensate, the dynamics of which we study in the chapter that follows. Sweeping through a range of temperatures and chemical potentials within the crossover region, we reveal the behaviour of phase fluctuations through analysis of correlation functions. We compare the ideal-gas approximation for the transition temperature with the more accurate value obtained via the Binder cumulant. By comparing the phase coherence length with the spatial extent across which density fluctuations are suppressed, we identify the BEC/quasicondensate crossover point. We compare directly with predictions of mean-field theory [Petrov et al., Phys. Rev. Lett. 87, 050404 (2001)] and find that phase fluctuations persist to lower temperatures in our SPGPE simulations. We calculate anomalous correlations of the fluctuation field and show that the resulting correction to the mean-field condensate self-interaction energy is significant within the quasicondensate regime. In Chapter 6 we study the formation of a BEC into a non-equilibrium state by simulating an experiment in which a dilute gas of Rubidium-87 atoms in a highly-elongated 3D harmonic magnetic trap is rapidly quenched to below the transition temperature. Following the quench, a condensate is observed to form with an excessively elongated shape, and subsequently undergo damped quadrupole shape oscillations [Shvarchuck et al., Phys. Rev. Lett. 89, 270404 (2002)]. We simulate this experiment by modifying the SPGPE to incorporate a hydrodynamic thermal cloud, parametrized by axially and temporally-dependent temperature and chemical potential. Comparison against simulations without this axial dependence suggests the hydrodynamic behaviour of the thermal cloud plays a crucial role in the formation of a non-equilibrium condensate. Analyses of the first-order correlation function and the condensate wavefunction reveal the spontaneous formation of quantum vortices, as well as strong phase fluctuations at equilibrium indicating quasicondensate behaviour.
Bright solitary-waves may be created in dilute Bose-Einstein condensates of at tractively interacting atoms in one dimensional regimes. In integrable systems, such solitary waves are particle-like objects called solitons. We investigate the consequences of non-integrability on the solitary waves in trapped Bose-Einstein condensates caused by an axial harmonic trap, and non-integrability caused by three dimensional effects. To analyse the soliton-like nature of the solitary-waves in an axial harmonic trap, a particle analogy for the solitary-waves is formulated. Exact Soliton solutions exist in the absence of an external trapping potential, which behave in a particle-like manner, and we find the particle analogy we employ to be a good model also when a harmonic trapping potential is present up to a gradual shift in the trajectories when the harmonic trap period is short compared with the, collision time of the solitons. We find that the collision time of the solions is dependent on the relative phase of the solitons as they collide. In the case of two solitons, the particle model is integrable, and the dynamics are completely regular. In the case of a system of two solitary waves of equal norm, the solitons are shown to retain their phase difference for repeated collisions. The extension to three particles supports both regular and chaotic regimes. The trajectory shift observed for two solitons carrier over to the case of three solitons. This shift aside, the agreement between the particle model and the wave dynamics remains good, even in chaotic regimes. We predict that these chaotic regimes will be an indicator of rapid depletion of the condensate due to quantum transitions of the condensate particles into non-condensate modes. To analyse the residual effects of the three dimensional nature of the solitary waves, we use a nonlinear Schrödinger equation with an additional quintic term. We perform variational calculations, and confirm the collapse of a soliton when the number of particles contained therein is increased past a critical number. We investigate the effects of varying the axial trap frequency and scattering length on the critical number. We propose a method to model particle exchange between solitons by extending the variational treatment to two solitons.
Trapped Bose-condensed atomic gases are remarkable because, in spite of the fact that these are very dilute systems, they exhibit robust coherent dynamic behaviour when perturbed. These quantum “wisps of matter” are a new phase of highly coherent matter. While binary collisions are very infrequent, the large coherent mean field associated with the Bose condensate ensures that interactions play a crucial role in determining the collective response of these trapped superfluid gases. In our discussion of the theory of collective oscillations of atomic condensates, the macroscopic Bose wavefunction Φ( r , t ) plays a central role. This wavefunction is the BEC order parameter. As discussed in Chapter 1, the initial attempts at defining this order parameter began with the pioneering work of London (1938a), were further developed by Bogoliubov (1947) and finally extended to deal with any Bose superfluid using the general quantum field theoretic formalism developed by Beliaev (1958a). Almost all this early theoretical work was limited to T = 0 where, in a dilute weakly interacting Bose gas, all the atoms are in the condensate. The first extension of these ideas to nonuniform Bose condensates was by Pitaevskii (1961) and, independently, by Gross (1961), which led to the now famous Gross–Pitaevskii (GP) equation of motion for Φ( r , t ). Before the discovery of BEC in trapped gases, the time-dependent GP equation was mainly used to study vortices in Bose superfluids, which involve a spatially nonuniform ground state. Apart from this application, the GP equation was largely unknown. The situation changed overnight in 1995 with the creation of trapped nonuniform Bose condensates in atomic gases.