We derive and study the hyperbolic reduced Maxwell-Bloch (HRMB) equations, a simplified model for the dynamics of strongly correlated Bose-Einstein condensates (BECs), and in particular for the interaction between the BEC atoms and its evaporated atoms under the strong interactions. This equation is one among four, which are proven to be integrable via the existence of a Lax pair, and thus the method of inverse scattering transform. Another equation is the reduced Maxwell-Bloch equation of quantum optics and the two others do not have physical applications yet. By studying the linear stability of the constant solutions of these four equations we observe various regimes, from stable, to modulational unstable, and unstable at all frequencies. The finite-dimensional reduction of the RMB equations is also used to give more insight into the constant solutions of these equations. From this study, we find that the HRMB equation arising from strongly correlated BECs is stable under the particular condition that the transition rate of evaporation is not too large compared to the number of evaporated atoms. We then derive explicit soliton solutions of the RMB equations and use numerical simulations to show collisions of solitons and kink solitons.
We derive and study the hyperbolic reduced Maxwell-Bloch (HRMB) equations, a simplified model for the dynamics of strongly correlated Bose-Einstein condensates (BECs), and in particular for the interaction between the BEC atoms and its evaporated atoms under the strong interactions. This equation is one among four, which are proven to be integrable via the existence of a Lax pair, and thus the method of inverse scattering transform. Another equation is the reduced Maxwell-Bloch equation of quantum optics and the two others do not have physical applications yet. By studying the linear stability of the constant solutions of these four equations we observe various regimes, from stable, to modulational unstable, and unstable at all frequencies. The finite-dimensional reduction of the RMB equations is also used to give more insight into the constant solutions of these equations. From this study, we find that the HRMB equation arising from strongly correlated BECs is stable under the particular condition that the transition rate of evaporation is not too large compared to the number of evaporated atoms. We then derive explicit soliton solutions of the RMB equations and use numerical simulations to show collisions of solitons and kink solitons.
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.
We present the theoretical and mathematical foundations of stability analysis for a Bose–Einstein condensate (BEC) at Dirac points of a honeycomb optical lattice. The combination of s -wave scattering for bosons and lattice interaction places constraints on the mean-field description, and hence on vortex configurations in the Bloch-envelope function near the Dirac point. A full derivation of the relativistic linear stability equations (RLSE) is presented by two independent methods to ensure veracity of our results. Solutions of the RLSE are used to compute fluctuations and lifetimes of vortex solutions of the nonlinear Dirac equation, which include Anderson–Toulouse skyrmions with lifetime s. Beyond vortex stabilities the RLSE provide insight into the character of collective superfluid excitations, which we find to encode several established theories of physics. In particular, the RLSE reduce to the Andreev equations, in the nonrelativistic and semiclassical limits, the Majorana equation, inside vortex cores, and the Dirac–Bogoliubov–de Gennes equations, when nearest-neighbor interactions are included. Furthermore, by tuning a mass gap, relative strengths of various spinor couplings, for the small and large quasiparticle momentum regimes, we obtain weak-strong Bardeen–Cooper–Schrieffer superconductivity, as well as fundamental wave equations such as Schrödinger, Dirac, Klein–Gordon, and Bogoliubov–de Gennes equations. Our results apply equally to a strongly spin–orbit coupled BEC in which the Laplacian contribution can be neglected.
We review our recent theoretical advances in the dynamics of Bose-Einstein condensates with tunable interactions using Feshbach resonance and external potential. A set of analytic and numerical methods for Gross-Pitaevskii equations are developed to study the nonlinear dynamics of Bose-Einstein condensates. Analytically, we present the integrable conditions for the Gross-Pitaevskii equations with tunable interactions and external potential, and obtain a family of exact analytical solutions for one- and two-component Bose-Einstein condensates in one and two-dimensional cases. Then we apply these models to investigate the dynamics of solitons and collisions between two solitons. Numerically, the stability of the analytic exact solutions are checked and the phenomena, such as the dynamics and modulation of the ring dark soliton and vector-soliton, soliton conversion via Feshbach resonance, quantized soliton and vortex in quasi-two-dimensional are also investigated. Both the exact and numerical solutions show that the dynamics of Bose-Einstein condensates can be effectively controlled by the Feshbach resonance and external potential, which offer a good opportunity for manipulation of atomic matter waves and nonlinear excitations in Bose-Einstein condensates.
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
The time-dependent extended Gross-Pitaevskii equation for Bose-Einstein condensates with attractive 1/r interaction is investigated with both a variational approach and numerically exact calculations. We show that these condensates exhibit signatures known from the nonlinear dynamics of autonomous Hamiltonian systems. The two stationary solutions created in a tangent bifurcation at a critical value of the scattering length are identified as elliptical and hyperbolical fixed points, corresponding to stable and unstable stationary states of the condensate. The stable stationary state is surrounded by elliptical islands, corresponding to condensates periodically oscillating in time, whereas condensate wave functions in the unstable region undergo a collapse within finite time. For negative scattering lengths below the tangent bifurcation no stationary solutions exist, i.e., the condensate is always unstable and collapses.