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.
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.
This dissertation focuses on the properties of nonlinear waves in Bose-Einstein condensates (BECs). The fundamental model here is the nonlinear Schrodinger equation, the so-called Gross-Pitaevskii (GP) equation, which is a mean-field description of BECs. The systematic analysis begins by considering the dark-bright (DB)-soliton interactions and multiple-dark-bright-soliton complexes in atomic two-component BECs. The interaction between two DB solitons in a homogeneous condensate and at the presence of the trap are both considered. Our analytical approximation relies in a Hamiltonian perturbation theory, which leads to an equation of motion of the centers of DB-soliton interacting pairs. Employing this equation, we demonstrate the existence of robust DB-soliton molecules, in the form of stationary two- and three-DB-soliton states. Also the equilibrium distance of the constituent solitons and the corresponding oscillation frequencies are found semianalytically, where the latter corresponds to the characteristic anomalous modes' eigenfrequencies that we numerically computed via a so called Bogoliubov-de Gennes (BdG) linearization analysis. Those studies are discussed in Chapter 2. Then, we extend our studies to the dynamics of dark-bright (DB) solitons in binary BECs at finite temperature using a system of two-coupled dissipative GPs. We show that the effect of the bright soliton is to partially stabilize dark solitons against temperature-induced dissipation, thus providing longer lifetimes in Chapter 3. Furthermore, the dark-dark (DD) solitons as a prototypical coherent structure that emerges in two-component BECs are studied and are connected to dark-bright (DB) solitons via SO(2) rotation. We obtained their beating frequency and their frequency of oscillation inside a parabolic trap. They are identified as exact periodic orbits in the Manakov limit of equal inter- and intra- species nonlinearity strengths with and without the trap and we showcase the persistence of such states upon weak deviations from this limit. Also we investigated in detail the effect of the deviation from the Manakov case by considering different from unity scattering length ratios in Chapter 4. Next, we revisited Hamiltonian eigenvalue problems that typically arise in the linearization around a stationary state of a Hamiltonian nonlinear PDE. Also we presented a overview of the known facts for the eigenvalue counts of the corresponding unstable spectra. In particular, we focused on a straightforward plan to implement finite-dimensional techniques for locating this spectrum via the singular points of the meromorphic Krein Matrix and illustrated the value of the approach by considering realistic problems for recently observed experimentally multivortex and multisoliton solutions in atomic Bose-Einstein condensates in Chapter 5. In the two dimensional scenario, we also examine the stability and dynamics of vortices under the effect of dissipation used as a simplified model for the inclusion of the effect of finite temperatures in atomic BECs, which enables an analytical prediction that can be compared directly to numerical results in Chapter 6. In all the above studies, our analytical prediction from the equation of motion are in good agreement with the numerical results from the BdG analysis.
In this thesis we investigate the static properties and non-equilibrium dynamics of bright solitary waves in atomic Bose-Einstein condensates in the zero-temperature limit, and we investigate the non-equilibrium dynamics of a driven atomic Bose-Einstein condensate at finite temperature. Bright solitary waves in atomic Bose-Einstein condensates are non-dispersive and soliton-like matter-waves which could be used in future atom-interferometry exper-iments. Using the mean-field, Gross-Pitaevskii description, we propose an exper-imental scheme to generate pairs of bright solitary waves with controlled velocity and relative phase; this scheme could form an important part of a future atom in-terferometer, and we demonstrate that it can also be used to test the validity of the mean-field model of bright solitary waves. We also develop a method to quanti-tatively assess how soliton-like static, three-dimensional bright solitary waves are; this assessment is particularly relevant for the design of future experiments. In reality, the non-zero temperatures and highly non-equilibrium dynamics occur-ring in a bright solitary wave interferometer are likely to necessitate a theoretical description which explicitly accounts for the non-condensate fraction. We show that a second-order, number-conserving description offers a minimal self-consistent treatment of the relevant condensate – non-condensate interactions at low temper-atures and for moderate non-condensate fractions. We develop a method to obtain a fully-dynamical numerical solution to the integro-differential equations of motion of this description, and solve these equations for a driven, quasi-one-dimensional test system. We show that rapid non-condensate growth predicted by lower-order descriptions, and associated with linear dynamical instabilities, can be damped by the self-consistent treatment of interactions included in the second-order descrip-tion. 1 Bright solitary waves and non-equilibrium dynamics in atomic Bose-Einstein condensates
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.
In this thesis, we have studied the behavior of two-component dark-bright solitons in multicomponent Bose-Einstein condensates (BECs) analytically and numerically in different situations. We utilized various analytical methods including the variational method and perturbation theory. By imprinting a linear phase on the bright component only, we were able to impart a velocity relative to the dark component and thereby we obtain an internal oscillation between the two components. We find that there are two modes of the oscillation of the dark-bright soliton. The first one is the famous Goldstone mode. This mode represents a moving dark-bright soliton without internal oscillation and is related to the continuous translational symmetry of the underlying equations of motion in the uniform potential. The second mode is the oscillation of the two components relative to each other. We compared the results obtained from the variational method with numerical simulations and found that the oscillation frequency range is 90 to 405 Hz and therefore observable in multicomponent Bose-Einstein condensate experiments. Also, we studied the binding energy and found a critical value for the breakup of the dark-bright soliton. Building on these results, we have studied another situation where we have the dark-bright soliton oscillate in a harmonic potential. We found for weak trapping the internal modes are nearly independent of center of mass motion of the dark-bright soliton. In contrast, in tighter traps the internal modes couple strongly to the center of mass motion, showing that for dark-bright solitons in a harmonic potential the center of mass and relative degrees of freedom are not independent. We found this result is robust against noise in the initial condition and should, therefore, be experimentally observable. In addition, we have studied the interaction between a moving dark-bright soliton in a uniform background with internal oscillation and a fixed impurity, modeled by a delta function potential. The interaction excites different modes in the system. Our analytical model capture two of these modes: the relative oscillation between the two components, as well as the in-sync oscillation of the widths. The numerical simulations allow further internal modes like out-of-sync oscillations of the soliton widths and even shape deformations of various kinds. We identify regions in parameter space for the transmission, reflection and inelastic scattering of the dark-bright soliton by the potential barrier. We have studied the velocity of dark-bright solitons described with an ansatz that uses one center of mass variable to represent the position of the two components. We found for a dark-bright soliton the maximum velocity is limited by the relative number of atoms in the bright component as compared to the size of the hole or density notch created by the dark component. Above this critical velocity the dark-bright soliton develops internal oscillations, and eventually unbinds and breaks apart.
This dissertation focuses on the properties of nonlinear waves in Bose-Einstein condensates (BECs). The fundamental model here is the nonlinear Schrodinger equation, the so-called Gross-Pitaevskii (GP) equation, which is a mean-field description of BECs. The systematic analysis begins by considering the dark-bright (DB)-soliton interactions and multiple-dark-bright-soliton complexes in atomic two-component BECs. The interaction between two DB solitons in a homogeneous condensate and at the presence of the trap are both considered. Our analytical approximation relies in a Hamiltonian perturbation theory, which leads to an equation of motion of the centers of DB-soliton interacting pairs. Employing this equation, we demonstrate the existence of robust DB-soliton molecules, in the form of stationary two- and three-DB-soliton states. Also the equilibrium distance of the constituent solitons and the corresponding oscillation frequencies are found semianalytically, where the latter corresponds to the characteristic anomalous modes' eigenfrequencies that we numerically computed via a so called Bogoliubov-de Gennes (BdG) linearization analysis. Those studies are discussed in Chapter 2. Then, we extend our studies to the dynamics of dark-bright (DB) solitons in binary BECs at finite temperature using a system of two-coupled dissipative GPs. We show that the effect of the bright soliton is to partially stabilize dark solitons against temperature-induced dissipation, thus providing longer lifetimes in Chapter 3. Furthermore, the dark-dark (DD) solitons as a prototypical coherent structure that emerges in two-component BECs are studied and are connected to dark-bright (DB) solitons via SO(2) rotation. We obtained their beating frequency and their frequency of oscillation inside a parabolic trap. They are identified as exact periodic orbits in the Manakov limit of equal inter- and intra- species nonlinearity strengths with and without the trap and we showcase the persistence of such states upon weak deviations from this limit. Also we investigated in detail the effect of the deviation from the Manakov case by considering different from unity scattering length ratios in Chapter 4. Next, we revisited Hamiltonian eigenvalue problems that typically arise in the linearization around a stationary state of a Hamiltonian nonlinear PDE. Also we presented a overview of the known facts for the eigenvalue counts of the corresponding unstable spectra. In particular, we focused on a straightforward plan to implement finite-dimensional techniques for locating this spectrum via the singular points of the meromorphic Krein Matrix and illustrated the value of the approach by considering realistic problems for recently observed experimentally multivortex and multisoliton solutions in atomic Bose-Einstein condensates in Chapter 5. In the two dimensional scenario, we also examine the stability and dynamics of vortices under the effect of dissipation used as a simplified model for the inclusion of the effect of finite temperatures in atomic BECs, which enables an analytical prediction that can be compared directly to numerical results in Chapter 6. In all the above studies, our analytical prediction from the equation of motion are in good agreement with the numerical results from the BdG analysis.
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.