We consider a binary Bose-Einstein condensate (BEC) described by a system of two-dimensional (2D) Gross-Pitaevskii equations with the harmonic-oscillator trapping potential. The intraspecies interactions are attractive, while the interaction between the species may have either sign. The same model applies to the copropagation of bimodal beams in photonic-crystal fibers. We consider a family of trapped hidden-vorticity (HV) modes in the form of bound states of two components with opposite vorticities ${S}_{1,2}=\ifmmode\pm\else\textpm\fi{}1$, the total angular momentum being zero. A challenging problem is the stability of the HV modes. By means of a linear-stability analysis and direct simulations, stability domains are identified in a relevant parameter plane. In direct simulations, stable HV modes feature robustness against large perturbations, while unstable ones split into fragments whose number is identical to the azimuthal index of the fastest growing perturbation eigenmode. Conditions allowing for the creation of the HV modes in the experiment are discussed too. For comparison, a similar but simpler problem is studied in an analytical form, viz., the modulational instability of an HV state in a one-dimensional (1D) system with periodic boundary conditions (this system models a counterflow in a binary BEC mixture loaded into a toroidal trap or a bimodal optical beam coupled into a cylindrical shell). We demonstrate that the stabilization of the 1D HV modes is impossible, which stresses the significance of the stabilization of the HV modes in the 2D setting.
We consider a binary Bose-Einstein condensate (BEC) described by a system of two-dimensional (2D) Gross-Pitaevskii equations with the harmonic-oscillator trapping potential. The intraspecies interactions are attractive, while the interaction between the species may have either sign. The same model applies to the copropagation of bimodal beams in photonic-crystal fibers. We consider a family of trapped hidden-vorticity (HV) modes in the form of bound states of two components with opposite vorticities S1,2=\ifmmode±\else\textpm\fi1, the total angular momentum being zero. A challenging problem is the stability of the HV modes. By means of a linear-stability analysis and direct simulations, stability domains are identified in a relevant parameter plane. In direct simulations, stable HV modes feature robustness against large perturbations, while unstable ones split into fragments whose number is identical to the azimuthal index of the fastest growing perturbation eigenmode. Conditions allowing for the creation of the HV modes in the experiment are discussed too. For comparison, a similar but simpler problem is studied in an analytical form, viz., the modulational instability of an HV state in a one-dimensional (1D) system with periodic boundary conditions (this system models a counterflow in a binary BEC mixture loaded into a toroidal trap or a bimodal optical beam coupled into a cylindrical shell). We demonstrate that the stabilization of the 1D HV modes is impossible, which stresses the significance of the stabilization of the HV modes in the 2D setting.
We consider two-dimensional (2D) localized vortical modes in the three-wave system with the quadratic ( χ ( 2 ) ) nonlinearity, alias nondegenerate second-harmonic (SH)-generating system, guided by the isotropic harmonic-oscillator (alias parabolic) confining potential. In addition to the straightforward realization in optics, the system models mixed atomic-molecular Bose–Einstein condensates. The main issue is stability of the vortex modes, which is investigated through computation of instability growth rates for eigenmodes of small perturbations, and by means of direct simulations. The threshold of parametric instability for single-color beams, represented solely by the SH with zero vorticity, is found in an analytical form with the help of the variational approximation. Trapped states with vorticities ( + 1 , − 1 , 0 ) in the two fundamental-frequency components and the SH one (the so-called hidden-vorticity modes) are completely unstable. Also unstable are semi-vortices , with component vorticities ( 1 , 0 , 1 ) . However, full vortices, with charges ( 1 , 1 , 2 ) , have a well-defined stability region. Unstable full vortices feature regions of robust dynamical behavior, where they periodically split and recombine, keeping their vortical content.
Summary form only given. I present results of the detailed asymptotic and numerical results on stability of unit and higher order vortices in coupled Gross-Pitaevskii equations, which model the experimental work. Low lying eigenfrequencies, including ones generating instability of the unit vortex placed in one component of a binary mixture of BEC, have been calculated applying perturbative theory in the limit of the weak interparticle interaction. Results of the perturbative theory have been verified and significantly extended using numerical analysis of the linearised equations for excitations and direct numerical simulation of the coupled Gross-Pitaevskii equations. Two main instability mechanisms have been identified and their dependencies on the relative number of the particles in the condensates and on the strength of the inter-species interaction have been studied. The first mechanism is the drift instability induced by the vortex free component of the condensate and resulting in the motion of the vortex out of the condensate center to its periphery. The second mechanism is the splitting of the higher order vortices into unit ones. Absolute and/or relative increase of the number of particles in the vortex free condensate component has been found to have stabilizing effect on both drift and splitting instabilities. Different scenarios of the splitting instabilities have been found.
We study the stability of vortices in a binary system of Bose-Einstein condensat es, with their wave functions modeled by a set of coupled, time-dependent Gross-Pitaevskii equations. Beginning with an effective two-dimensional system, we identify miscible and immiscible regimes characterized by the inter- and intra-atomic interactions and the initial configuration of the system. We then consider a binary system of Bose-Einstein condensates placed in a rotating harmonic trap and study the single united vortex state in this system. We derive an approximate form for the energy of a single united vortex in the binary system and the critical angular velocity for the global stability of a vortex at the center of the trap. We also compute the metastability onset angular velocity for the local stability of a vortex at the center of the trap. Numerical solutions to the Gross-Pitaevskii equations support these expressions. These rotational results inform us of a novel bimodal formation within the miscible regime of the binary condensate system. We thus demonstrate the nontrivial aspects of vortex stability in interacting binary Bose-Einstein condensates as a result of their nonlinear interactions.
We consider a two-dimensional (2D) counterpart of the experiment that led to the creation of quasi-1D bright solitons in Bose-Einstein condensates (BECs) [Nature 417, 150--153 (2002)]. We start by identifying the ground state of the 2D Gross-Pitaevskii equation for repulsive interactions, with a harmonic-oscillator (HO) trap, and with or without an optical lattice (OL). Subsequently, we switch the sign of the interaction to induce interatomic attraction and monitor the ensuing dynamics. Regions of the stable self-trapping and catastrophic collapse of 2D fundamental solitons are identified in the parameter plane of the OL strength and BEC norm. The increase of the OL strength expands the persistence domain for the solitons to larger norms. For single-charged solitary vortices, in addition to the survival and collapse regimes, an intermediate one is identified, where the vortex resists the collapse but loses its structure, transforming into a fundamental soliton. The same setting may also be implemented in the context of optical solitons and vortices, using photonic-crystal fibers.
We introduce two- and one-dimensional (1D) models of a binary BEC (Bose-Einstein condensate) in a periodic potential, with repulsive interactions. We chiefly consider the most fundamental case of the interspecies repulsion with zero intraspecies interactions. The same system may also model a mixture of two mutually repulsive fermionic species. Existence and stability regions for gap solitons (GSs) supported by the interplay of the interspecies repulsion and periodic potential are identified. Two-component GSs are constructed by means of the variational approximation (VA) and in a numerical form. The VA provides an accurate description for the GS which is a compound state of two tightly bound components, each essentially trapped in one cell of the periodic potential. GSs of this type dominate in the case of intragap solitons, with both components belonging to the first finite bandgap of the linear spectrum (only this type of solitons is possible in a weak lattice). Intergap solitons, with one component residing in the second bandgap, and intragap solitons which have both components in the second gap, are possible in a deeper periodic potential, with the strength essentially exceeding the recoil energy of the atoms. Intergap solitons are, typically, bound states of one tightly and one loosely bound component. In this case, results are obtained in a numerical form. The number of atoms in experimentally relevant situations is estimated to be $\ensuremath{\sim}5000$ in the 2D intragap soliton, and $\ensuremath{\sim}25\phantom{\rule{0.2em}{0ex}}000$ in its intergap counterpart; in 1D solitons, it may be up to 105. For 2D solitons, the stability is identified in direct simulations, while in the 1D case it is done via eigenfrequencies of small perturbations, and then verified by simulations. In the latter case, if the intragap soliton in the first bandgap is weakly unstable, it evolves into a stable breather, while unstable solitons of other types (in particular, intergap solitons) get completely destroyed. The intragap 2D solitons in the first bandgap are less robust, and in some cases they are completely destroyed by the instability. Addition of intraspecies repulsion to the repulsion between the components leads to further stabilization of the GSs.
Bose-Einstein condensation(BEC) has been a research hotspot in cold atom physics for the past few decades. Recently, experiments on BECs have reported the observation of liquid-like quantum droplets, attracting significant interest. The physics for forming quantum droplets can be explained by zero-point quantum fluctuation which effectively prohibit the collapse caused by the mean-field interactions between the atoms. The stabilization mechanism beyond the mean-field theory, induced by quantum fluctuations, can be described by the Lee–Huang–Yang(LHY) correction. Normally, forming stable 2D or 3D solitons in free space is a challenging problem. Additionally, an intriguing question is how to generate 2D or 3D solitons with embedded vorticity. The LHY stability mechanism offers a new approach to producing stable solitons and vortex solitons in BECs. In this thesis, two-dimensional vortex quantum droplets in Bose-Einstein condensates, including those with and without external potential, are examined. The study of solitons and vortex solitons in BECs under mean-field interactions and external potential has yielded significant results. This thesis expands this field by introducing the LHY correction to study stable quantum droplets and vortex droplets in BECs. The work chapters begin with studying binary quantum droplets in free space, considering the case where interspecies interactions are weakly attractive and intraspecies interactions are weakly repulsive. The zero-vorticity, semi-vortex, and binary vortex quantum droplets are studied in Chapter 3. Stable semi-vortex quantum droplets with embedded vorticity up to 8 are found in this system. This result gives a way to form stable vortex droplets with high topological charges in free space by controlling one component with zero vorticity and the other component with embedded vorticity. The stability and properties of these droplets are systematically studied. Chapter 4 explores the vortex quantum droplets and hidden vortex droplets confined within a radial lattice trap. It is found that the existence curves of the droplets confined within this optical trap could violate the Vakhitov-Kolokolov(VK) criterion, which is a necessary condition to form a stable soliton. A remarkable finding is that the vortex and hidden vortex quantum droplets in this radial lattice remain stable when they carry the embedded vorticity up to 11. The double-ring vortex and hidden vortex quantum droplets confined within this radial lattice are also discussed. The work expands from non-magnetic atoms to magnetic atoms in Chapter 5, which needs to consider the dipole-dipole interactions(DDIs). The anisotropic dipolar vortex quantum droplets trapped in an annular potential are studied. Under the effect of the ring shape potential, the double peak ring-shaped vortex droplets remain stable when they carry the embedded vorticity up to 11. It is the first time stable droplets with high topological charges have been formed in a dipolar system. Additionally, the shape transformation of dipolar vortex droplets induced by the strength of the mean-field interaction is examined. Finally, the behavior of dipolar vortex droplets confined within a weak-depth annular potential is also discussed. The results of this thesis may prove highly useful in creating stable vortex quantum droplets in Bose-Einstein condensates and advancing the broader field of ultracold atoms.