It is commonly known that two-dimensional mean-field models of optical and matter waves with cubic self-attraction cannot produce stable solitons in free space because of the occurrence of collapse in the same setting. By means of numerical analysis and variational approximation, we demonstrate that the two-component model of the Bose-Einstein condensate with the spin-orbit Rashba coupling and cubic attractive interactions gives rise to solitary-vortex complexes of two types: semivortices (SVs, with a vortex in one component and a fundamental soliton in the other), and mixed modes (MMs, with topological charges 0 and ±1 mixed in both components). These two-dimensional composite modes can be created using the trapping harmonic-oscillator (HO) potential, but remain stable in free space, if the trap is gradually removed. The SVs and MMs realize the ground state of the system, provided that the self-attraction in the two components is, respectively, stronger or weaker than the cross attraction between them. The SVs and MMs which are not the ground states are subject to a drift instability. In free space (in the absence of the HO trap), modes of both types degenerate into unstable Townes solitons when their norms attain the respective critical values, while there is no lower existence threshold for the stable modes. Moving free-space stable solitons are also found in the present non-Galilean-invariant system, up to a critical velocity. Collisions between two moving solitons lead to their merger into a single one.
It is commonly known that two-dimensional mean-field models of optical and matter waves with cubic self-attraction cannot produce stable solitons in free space because of the occurrence of collapse in the same setting. By means of numerical analysis and variational approximation, we demonstrate that the two-component model of the Bose-Einstein condensate with the spin-orbit Rashba coupling and cubic attractive interactions gives rise to solitary-vortex complexes of two types: semivortices (SVs, with a vortex in one component and a fundamental soliton in the other), and mixed modes (MMs, with topological charges 0 and ±1 mixed in both components). These two-dimensional composite modes can be created using the trapping harmonic-oscillator (HO) potential, but remain stable in free space, if the trap is gradually removed. The SVs and MMs realize the ground state of the system, provided that the self-attraction in the two components is, respectively, stronger or weaker than the cross attraction between them. The SVs and MMs which are not the ground states are subject to a drift instability. In free space (in the absence of the HO trap), modes of both types degenerate into unstable Townes solitons when their norms attain the respective critical values, while there is no lower existence threshold for the stable modes. Moving free-space stable solitons are also found in the present non-Galilean-invariant system, up to a critical velocity. Collisions between two moving solitons lead to their merger into a single one.
By means of variational methods and systematic numerical analysis, we demonstrate the existence of metastable solitons in three dimensional (3D) free space, in the context of binary atomic condensates combining contact self-attraction and spin-orbit coupling, which can be engineered by available experimental techniques. Depending on the relative strength of the intra- and intercomponent attraction, the stable solitons feature a semivortex or mixed-mode structure. In spite of the fact that the local cubic self-attraction gives rise to the supercritical collapse in 3D, and hence the setting produces no true ground state, the solitons are stable against small perturbations, motion, and collisions.
We report families of two-dimensional (2D) composite solitons in spinor dipolar Bose-Einstein condensates, with two localized components linearly mixed by the spin-orbit coupling (SOC), and the intrinsic nonlinearity represented by the dipole-dipole interaction (DDI) between atomic magnetic moments polarized in plane by an external magnetic field. Recently, stable solitons were predicted in the form of semivortices (composites built of coupled fundamental and vortical components) in the 2D system combining the SOC and contact attractive interactions. Replacing the latter by the anisotropic long-range DDI, we demonstrate that, for a fixed norm of the soliton, the system supports a continuous family of stable spatially asymmetric vortex solitons (AVSs), parameterized by an offset of the pivot of the vortical component relative to its fundamental counterpart. The offset is limited by a certain maximum value, while the energy of the AVS practically does not depend on the offset. At small values of the norm, the vortex solitons are subject to a weak oscillatory instability. In the present system, with the Galilean invariance broken by the SOC, the composite solitons are set in motion by a kick the strength of which exceeds a certain depinning value. The kicked solitons feature a negative effective mass, drifting along a spiral trajectory opposite to the direction of the kick. A critical angular velocity, up to which the semivortices may follow rotation of the polarizing magnetic field, is found too.
It is known that two-dimensional two-component fundamental solitons of the semivortex (SV) type, with vorticities (s_{+},s_{-})=(0,1) in their components, are stable ground states (GSs) in the spin-orbit-coupled (SOC) binary Bose-Einstein condensate with the contact self-attraction acting in both components, in spite of the possibility of the critical collapse in the system. However, excited states (ESs) of the SV solitons, with the vorticity set (s_{+},s_{-})=(S_{+},S_{+}+1) and S_{+}=1,2,3,..., are unstable in the same system. We construct ESs of SV solitons in the SOC system with opposite signs of the self-interaction in the two components. The main finding is stability of the ES-SV solitons, with the extra vorticity (at least) up to S_{+}=6. The threshold value of the norm for the onset of the critical collapse, N_{thr}, in these excited states is higher than the commonly known critical value, N_{c}≈5.85, associated with the single-component Townes solitons, N_{thr} increasing with the growth of S_{+}. A velocity interval for stable motion of the GS-SV solitons is found too. The results suggest a solution for the challenging problem of the creation of stable vortex solitons with high topological charges.
We study two-dimensional (2D) matter-wave solitons in spinor Bose–Einstein condensates under the action of the spin–orbit coupling and opposite signs of the self- and cross-interactions. Stable 2D two-component solitons of the mixed-mode type are found if the cross-interaction between the components is attractive, while the self-interaction is repulsive in each component. Stable solitons of the semi-vortex type are formed in the opposite case, under the action of competing self-attraction and cross-repulsion. The solitons exist with the total norm taking values below a collapse threshold. Further, in the case of the repulsive self-interaction and inter-component attraction, stable 2D self-trapped modes, which may be considered as quantum droplets (QDs), are created if the beyond-mean-field Lee–Huang–Yang terms are added to the self-repulsion in the underlying system of coupled Gross–Pitaevskii equations. Stable QDs of the mixed-mode type, of a large size with an anisotropic density profile, exist with arbitrarily large values of the norm, as the Lee–Huang–Yang terms eliminate the collapse. The effect of the spin–orbit coupling term on characteristics of the QDs is systematically studied. We also address the existence and stability of QDs in the case of SOC with mixed Rashba and Dresselhaus terms, which makes the density profile of the QD more isotropic. Thus, QDs in the spin-orbit-coupled binary Bose–Einstein condensate are for the first time studied in the present work.
We introduce one- and two-dimensional (1D and 2D) continuum and discrete models for the two-component BEC, with the spin-orbit (SO) coupling of the Rashba type between the components, and attractive cubic interactions, assuming that the condensate is fragmented into a quasidiscrete state by a deep optical-lattice potential. In 1D, it is demonstrated, in analytical and numerical forms, that the ground states of both the discrete system and its continuum counterpart switch from striped bright solitons, featuring deep short-wave modulations of its profile, to smooth solitons, as the strength ratio of the inter- and intracomponent attraction, γ, changes from γ<1 to γ>1. At the borderline, γ=1, there is a continuous branch of stable solitons, which share a common value of the energy and interpolate between the striped and smooth ones. Unlike the 2D system, the 1D solitons, which do not represent the ground state at given γ, are nevertheless stable against small perturbations, and they remain stable too in collisions with other solitons. In 2D, a transition between two different types of discrete solitons, which represent the ground state, viz., semivortices and mixed modes, also takes place at γ=1. A specific property of 2D discrete solitons of both types is their discontinuous transition into a delocalized state at a critical value of the SO-coupling strength. We also address the continuum 2D model in the borderline case of γ=1, which was not studied previously, and demonstrate the existence of an energy-degenerate branch of dynamically stable solitons connecting the semivortex and the mixed mode. Last, it is demonstrated that 1D and 2D discrete solitons are mobile, in a limited interval of velocities.
We present an analysis of two-dimensional (2D) matter-wave solitons, governed by the pseudospinor system of Gross-Pitaevskii equations with self- and cross attraction, which includes the spin-orbit coupling (SOC) in the general Rashba-Dresselhaus form, and, separately, the Rashba coupling and the Zeeman splitting. Families of semivortex (SV) and mixed-mode (MM) solitons are constructed, which exist and are stable in free space, as the SOC terms prevent the onset of the critical collapse and create the otherwise missing ground states in the form of the solitons. The Dresselhaus SOC produces a destructive effect on the vortex solitons, while the Zeeman term tends to convert the MM states into the SV ones, which eventually suffer delocalization. Existence domains and stability boundaries are identified for the soliton families. For physically relevant parameters of the SOC system, the number of atoms in the 2D solitons is limited by ∼1.5×10^{4}. The results are obtained by means of combined analytical and numerical methods.