We study statically homogeneous Bose-Einstein condensates with spatially inhomogeneous interactions and outline an experimental realization of compensating linear and nonlinear potentials that can yield constant-density solutions. We illustrate how the presence of a step in the nonlinearity coefficient can only be revealed dynamically and examine how to reveal it by exploiting the inhomogeneity of the sound speed with a defect-dragging experiment. We conduct computational experiments and observe the spontaneous emergence of dark solitary waves. We use effective-potential theory to perform a detailed analytical investigation of the existence and stability of solitary waves in this setting, and we corroborate these results computationally using a Bogoliubov--de Gennes linear stability analysis. We find that dark solitary waves are unstable for all step widths, whereas bright solitary waves can become stable through a symmetry-breaking bifurcation as one varies the step width. Using phase-plane analysis, we illustrate the scenarios that permit this bifurcation and explore the dynamical outcomes of the interaction between the solitary wave and the step.
We study statically homogeneous Bose-Einstein condensates with spatially inhomogeneous interactions and outline an experimental realization of compensating linear and nonlinear potentials that can yield constant-density solutions. We illustrate how the presence of a step in the nonlinearity coefficient can only be revealed dynamically and examine how to reveal it by exploiting the inhomogeneity of the sound speed with a defect-dragging experiment. We conduct computational experiments and observe the spontaneous emergence of dark solitary waves. We use effective-potential theory to perform a detailed analytical investigation of the existence and stability of solitary waves in this setting, and we corroborate these results computationally using a Bogoliubov--de Gennes linear stability analysis. We find that dark solitary waves are unstable for all step widths, whereas bright solitary waves can become stable through a symmetry-breaking bifurcation as one varies the step width. Using phase-plane analysis, we illustrate the scenarios that permit this bifurcation and explore the dynamical outcomes of the interaction between the solitary wave and the step.
We study the structure, stability, and dynamics of dark solitary waves in parabolically trapped, collisionally inhomogeneous Bose-Einstein condensates (BECs) with spatially periodic variations of the scattering length. This collisional inhomogeneity yields a nonlinear lattice, which we tune from a small-amplitude, approximately sinusoidal structure to a periodic sequence of densely spaced spikes. We start by investigating time-independent inhomogeneities, and we subsequently examine the dynamical response when one starts with a collisionally homogeneous BEC and then switches on an inhomogeneity either adiabatically or nonadiabatically. Using Bogoliubov-de Gennes linearization as well as direct numerical simulations of the Gross-Pitaevskii equation, we observe dark solitary waves, which can become unstable through oscillatory or exponential instabilities. We find a critical wavelength of the nonlinear lattice that is comparable to the healing length. Near this value, the fundamental eigenmode responsible for the stability of the dark solitary wave changes its direction of movement as a function of the strength of the nonlinearity. When it increases, it collides with other eigenmodes, leading to oscillatory instabilities; when it decreases, it collides with the origin and becomes imaginary, illustrating that the instability mechanism is fundamentally different in wide-well versus narrow-well lattices. When starting from a collisionally homogeneous setup and switching on inhomogeneities, we find that dark solitary waves are preserved generically for aligned lattices. We briefly examine the time scales for the onset of solitary-wave oscillations in this scenario.
In the present work, we revisit two-component Bose-Einstein condensates in their fully three-dimensional (3D) form. Motivated by earlier studies of dark-bright solitons in the 1D case, we explore the stability of these structures in their fully 3D form in two variants. In one the dark soliton is planar and trapping a planar bright (disk) soliton. In the other case, a dark spherical shell soliton creates an effective potential in which a bright spherical shell of atoms is trapped in the second component. We identify these solutions as numerically exact states (up to a prescribed accuracy) and perform a Bogolyubov-de Gennes linearization analysis that illustrates that both structures can be dynamically stable in suitable intervals of sufficiently low chemical potentials. We corroborate this finding theoretically by analyzing the stability via degenerate perturbation theory near the linear limit of the system. When the solitary waves are found to be unstable, we explore their dynamical evolution via direct numerical simulations which, in turn, reveal wave forms that are more robust. Finally, using the SO(2) symmetry of the model, we produce multi-dark-bright planar or shell solitons involved in pairwise oscillatory motion.
In recent years, bright soliton-like structures composed of gaseous Bose-Einstein condensates have been generated at ultracold temperature. The experimental capacity to precisely engineer the nonlinearity and potential landscape experienced by these solitary waves offers an attractive platform for fundamental study of solitonic structures. The presence of three spatial dimensions and trapping implies that these are strictly distinct objects to the true soliton solutions. Working within the zero-temperature mean-field description, we explore the solutions and stability of bright solitary waves, as well as their interactions. Emphasis is placed on elucidating their similarities and differences to the true bright soliton. The rich behaviour introduced in the bright solitary waves includes the collapse instability and symmetry-breaking collisions. We review the experimental formation and observation of bright solitary matter waves to date, and compare to theoretical predictions. Finally we discuss the current state-of-the-art of this area, including beyond-mean-field descriptions, exotic bright solitary waves, and proposals to exploit bright solitary waves in interferometry and as surface probes.
Abstract We systematically construct vector solitary waves in harmonically trapped one-dimensional two-component Bose–Einstein condensates with unequal dispersion coefficients by a numerical continuation in chemical potentials from the respective analytic low-density linear limits to the high-density nonlinear Thomas-Fermi regime. The main feature of the linear states herein is that the component with the larger quantum number has instead a smaller linear eigenenergy, enabled by suitable unequal dispersion coefficients, leading to new series of solutions compared with the states similarly obtained in the equal dispersion setting. Particularly, the lowest-lying series gives the well-known dark-anti-dark waves, and the second series yields the dark-multi-dark states, and the following series become progressively more complex in their wave structures. The Bogoliubov-de Gennes spectra analysis shows that most of these states bear unstable modes, but they can be long-lived and remarkably all of them can be fully stabilized in suitable parameter regimes.
Dilute atomic Bose-Einstein condensates are inherently nonlinear systems and support solitary wave solutions. An important distinction from optical systems is the inhomogeneous background density, which results from the traps used to confine the atoms. As in optical systems, dark solitary waves in three-dimensional geometries are unstable to transverse excitations, which lead to a bending of the dark soliton plane and decay into vortex rings. Highly elon-gated geometries can now be achieved experimentally, in which the condensate dynamics are effectively one-dimensional, and the motion of the dark soliton is governed by the inhomogeneous longitudinal density. We show that a dark soliton is fundamentally unstable to such a changing background density, by means of numerical simulations of the soliton under various potentials (e.g. steps, ramps, harmonic traps, and optical lattices). This leads to the emission of radiation in the form of sound waves. The power emitted is found to be proportional to the square of the soliton acceleration. The latter quantity is shown to be proportional to the deformation of the apparent soliton profile, arising from the sound field in the region of the soliton. We demonstrate that the ensuing interactions between the soliton and sound field, and therefore the dynamics of the soliton, can be controlled experimentally via manipulation of the emitted sound, achieved by modifying the trap geometry. In this manner, it is possible to induce a rapid decay of the soliton, stabilise the soliton, or even pump energy into the soliton by means of parametric driving. 1
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.