In the present work, we explore the potential of spin-orbit (SO) coupled Bose-Einstein condensates to support multi-component solitonic states in the form of dark-bright (DB) solitons. In the case where Raman linear coupling between components is absent, we use a multiscale expansion method to reduce the model to the integrable Mel'nikov system. The soliton solutions of the latter allow us to reconstruct approximate traveling DB solitons for the reduced SO coupled system. For small values of the formal perturbation parameter, the resulting waveforms propagate undistorted, while for large values thereof, they shed some dispersive radiation, and subsequently distill into a robust propagating structure. After quantifying the relevant radiation effect, we also study the dynamics of DB solitons in a parabolic trap, exploring how their oscillation frequency varies as a function of the bright component mass and the Raman laser wavenumber.
In the present work, we explore the potential of spin-orbit (SO) coupled Bose-Einstein condensates to support multi-component solitonic states in the form of dark-bright (DB) solitons. In the case where Raman linear coupling between components is absent, we use a multiscale expansion method to reduce the model to the integrable Mel'nikov system. The soliton solutions of the latter allow us to reconstruct approximate traveling DB solitons for the reduced SO coupled system. For small values of the formal perturbation parameter, the resulting waveforms propagate undistorted, while for large values thereof, they shed some dispersive radiation, and subsequently distill into a robust propagating structure. After quantifying the relevant radiation effect, we also study the dynamics of DB solitons in a parabolic trap, exploring how their oscillation frequency varies as a function of the bright component mass and the Raman laser wavenumber.
We investigate the different types of matter-wave solitons in spin-orbit-coupled spin-2 spinor Bose-Einstein condensates. Using mean-field theory and adopting the multiscale perturbation method, the original five-component Gross-Pitaevskii spin-orbit-coupled spin-2 spinor Bose-Einstein condensate model can be reduced to a single effective nonlinear Schrödinger equation, which allows us to find analytical soliton solutions of this system. In this way, for different regimes of the spin-orbit coupling, Raman coupling, and interatomic interactions, we find approximate bright and dark soliton solutions. Particularly, the type of solitons depends on the dispersion properties of the system. When the lowest-energy band has a single-well structure, we find there only exist positive mass bright or dark solitons due to the dispersion coefficient of effective nonlinear Shrödinger equation always positive. However, when the lowest-energy band has a double-well structure, there will appear positive (negative) mass bright or dark solitons because the sign of the dispersion coefficient can be positive (negative) under different momentum. We employ direct numerical simulation of the original five-component Gross-Pitaevskii equations to confirm the analytical results.
We consider vector solitons of mixed bright-dark types in quasi-one-dimensional spinor (F=1) Bose-Einstein condensates. Using a multiscale expansion technique, we reduce the corresponding nonintegrable system of three coupled Gross-Pitaevskii equations (GPEs) to an integrable Yajima-Oikawa system. In this way, we obtain approximate solutions for small-amplitude vector solitons of dark-dark-bright and bright-bright-dark types, in terms of the ${m}_{F}=+1,\ensuremath{-}1,0$ spinor components, respectively. By means of numerical simulations of the full GPE system, we demonstrate that these states indeed feature soliton properties, i.e., they propagate undistorted and undergo quasielastic collisions. It is also shown that in the presence of a parabolic trap the bright component(s) is (are) guided by the dark one(s) and, as a result, the small-amplitude vector soliton as a whole performs quasiharmonic oscillations. The oscillation frequency is found as a function of the spin-dependent interaction strength for both small-amplitude and large-amplitude solitons.
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.
We present families of beating dark-dark solitons in spin-orbit (SO) -coupled Bose-Einstein condensates. These families consist of solitons residing simultaneously in the two bands of the energy spectrum. The soliton components are characterized by two different spatial and temporal scales, which are identified by a multiscale expansion method. The solitons are ``beating'' ones, as they perform density oscillations. The characteristic frequency of the latter is relevant to Zitterbewegung (ZB) oscillations, which were recently observed in experiments with SO-coupled condensates [C. Qu et al., Phys. Rev. A 88, 021604(R) (2013); L. J. LeBlanc et al., New J. Phys. 15, 073011 (2013)]. We find that spin oscillations may occur, depending on the parity of each soliton branch, which consequently lead to ZB oscillations of the beating dark solitons. Analytical results are corroborated by numerical simulations, illustrating the robustness of the 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.
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.