Two-Dimensional Vortex Solitons in Spin-Orbit-Coupled Dipolar Bose – Einstein Condensates

Solitons are self-trapped modes existing in various nonlinear systems. Creating stable solitons in twoand three-dimensional settings is a challenging goal in various branches of physics. Several methods have been developed theoretically and experimentally to achieve this, but few of them can support stable multi-dimensional solitons in free space. Recently, a new scheme using spin-orbit-coupling (SOC) has been proposed to create stable 2D solitons in Bose–Einstein condensates (BECs). This paper reviews recent theoretical progress on creating stable 2D solitons in spinor dipolar BEC with SOC, combined with long-range dipole-dipole interaction (DDI), Zeeman splitting (ZS) and contact nonlinearity, in free space. The continuous family of stable symmetric vortex solitons (SVS), asymmetric vortex solitons (AVS), as well as gap solitons (GS) is found via different settings. Their existence and stability conditions are summarized and discussed in detail. The mobility properties of these types of solitons are also addressed. For SVS, a potential method to manipulate its shape and mobility is investigated. These results are supposed to enrich our understanding of 2D solitons and help create multi-dimensional solitons in experiments.


Introduction
Creating stable solitons or self-trapped modes in systems with a dimension higher than one remains a challenging and attractive problem in physics, both in theory and experiments [1,2].The main obstacle lies in the ubiquitous cubic nonlinearity, which could be the Kerr effect in optics [3] or the intrinsic ones induced by two-body scattering in Bose-Einstein condensates (BECs) [4][5][6].Although these nonlinearities can support stable soliton-like modes in 1D geometry, in high dimensions like 2D or 3D systems, the collapse driven by cubic self-focusing nonlinearities will destroy the stable modes [7][8][9].Several methods have been developed, theoretically and experimentally, to solve this problem.One is using alternative nonlinearities other than cubic, such as quadratic self-focusing nonlinearity, which can support stable 2D spatial solitons [10] and spatiotemporal soliton chains [11] in optics, mixed cubic-quintic or quadratic-cubic ones, which can support stable fundamental spatial solitons [12][13][14][15][16].The other is introducing spatially periodic potentials, which is known to be able to help stabilize 2D and even 3D fundamental and vortical solitons in plentiful settings [17][18][19][20], and the spectral band these periodical potentials host can create gap solitons, as well [17,[21][22][23][24].Another method is to use nonlocal nonlinearities, which can came from the long-range interaction between the molecules in liquid crystals [25], the diffusive behavior in nonlinear medium [26,27], the van der Waals interactions between Rydberg atoms formed in BEC [28], the dipole-dipole interaction (DDI) and quadrupolar interaction between the permanent [29,30] or the external fields' induced [31][32][33][34][35][36] magnetic or electric dipole moments of the BEC atoms.
Besides the above methods, there is a new approach recently found theoretically by using a two-component BEC with linear spin-orbit coupling (SOC) between the two components.Combined with repulsive nonlinearity, such a system is known to support 1D and 2D bright solitons [37,38], 2D dark vortices [39][40][41][42] and bright gap solitons (with an optical lattice) [43][44][45][46].With intrinsic attractive nonlinearity, BEC with SOC was thought to be able to support 1D solitons only [37,47,48], but new work revealed, very surprisingly, that with Rashba-type SOC (see Equation (2)), BEC with attractive cubic nonlinearity can support stable 2D semi-and mixed-mode vortex solitons in free space [49][50][51][52][53].A similar mechanism can also be established with optical settings [54,55], and the excited state soliton can be created by introducing the inverted nonlinear potential [56,57].Combined with long-range interactions, dipolar BEC with SOC is found to support 2D bright solitons, as well, which could be a plane-wave soliton with a spatially varying phase or a stripe soliton with a spatially oscillating density for each component [38].These findings not only open new possibilities to create stable multi-dimensional solitons, but can be used to investigate 2D or 3D solitons with vortices and other topological properties, as well.In this article, we reviewed recent theoretical results of creating stable 2D solitons in free space by using BEC with SOC and DDI nonlinearity.In the setting of SOC BEC with pure DDI interactions, families of stable symmetric semivortices (SVS) are found, and the properties of these SVS families are quantified [58].Families of stable asymmetric semivortices (AVS) are also found in this setting, which are absent in systems with contact interactions [58].Control parameters and stable conditions of these AVS families are discussed and analyzed [58].Both families of soliton are anisotropic in space due to the anisotropic nature of the DDI.By introducing a Zeeman splitting (ZS) into this setting, it is found that the shape of the SVS can be controlled.By increasing the strength of the Zeeman splitting, a transition from horizontal SVS to vertical SVS occurs [59].The relation of the transition with the controlled parameters is analyzed numerically, and a physical explanation is given.When the ZS switched to the opposite direction, families of inverted SVS can also stably exist in the system [59].The mobility of these SVS was analyzed numerically, as well [58,59].Settings of SOC dipolar BEC combined with cubic contact interactions are also considered [60], and stable isotropic and anisotropic gap solitons are found and their properties and stable conditions analyzed, as well as the mobility of these solitons.These results will help to investigate 2D stable vortices in an SOC BEC system experimentally in the future and provide potential manipulation methods [58][59][60].
This review is outlined as follows.Section 2 introduces basic models and parameter settings.Section 3 discusses two types of stable vortex solitons (SVS and AVS) in a spinor SOC BEC with pure DDI, and their properties and mobilities are also investigated.Section 4 discusses the effect when Zeeman splitting is introduced, and a new type of vortex soliton (inverted SVS) is found stable in this setting, which serves as an excited state of the system.By dropping the kinetic term, an energy gap can be opened by the Zeeman splitting, the gap solitons and their properties are discussed in Section 5.The review is concluded in Section 6.

The Model of SOC BEC with Dipole-Dipole Interactions
In a two-component (spinor) BEC, a linear SOC between the two components can be induced by the combination of optical and magnetic fields [61][62][63].Under mean filed approximation, the evolution of the spinor BEC wavefunction, ψ = (ψ + , ψ − ), in a 2D setting is governed by the coupled Gross-Pitaevskii equation, which can be written in a normalized full form: In the equation, D[±] is the anisotropic Rashba SO-coupling term: where λ x,y ∈ [0, 1] is the strength of SO-coupling in the x and y direction, respectively.Ω represents the strength of Zeeman splitting between the two components, which may be replaced by the Stark-Lo Surdo splitting induced by the DC electric field.g and γ are coefficients of self-and cross-interactions of the components.β is the strength of the DDI interactions; for simplicity, only density-density interactions are considered; thus, the exchange and relaxation term between different spin components has been neglected [38] (the full DDI interaction term can be seen in [62]).What is more, as the following results will show that contact nonlinearity is not a dominant factor for the physics studied here, for convenience, the intra-component contact interactions are taken to be equal.R(r − r ) is the DDI's kernel, which reads: where cutoff is determined by the confinement in the transverse (third) dimension [30,60,64,65], all the dipoles are polarized to the same direction and θ is the angle between the polarization direction and (r − r ).Note that the sign of Ω, g, γ, β can all be switched between positive and negative experimentally [58][59][60].
Depending on the specific settings addressed, Equation (1) reduces to different forms, which are listed briefly as follows and investigated in detail in the following sections.

Semi-Vortices in Spinor SOC BEC with Pure DDI
In Equation (1), if Ω, g, γ is set to be zero and β is chosen to be > 0 (Setting (i)), the system then corresponds to a spinor SOC BEC with pure DDI, and Equation ( 1) is reduced to: The DDI kernel follows Equation (3) and cos 2 θ ≡ (x − x ) 2 /|r − r | 2 , which means the dipoles are all polarized in the 2D plane along the positive direction of the x axis.In order to obtain 2D solitons, the size of the soliton needs to be much larger than the cutoff , which requires proper setting of and λ.In this section, = 0.05, λ = 1 is fixed [58].
The SOC term in Equation (4) breaks the Galilean invariance of the system, while the SOC and DDI terms both break the conservation of the angular momentum of the system.Thus, for spinor SOC BEC with pure DDI, the system is invariant with respect to specular reflections in the x-y plane: where * denotes the complex conjugate, which provides a combined symmetry of the system, which reads: The localized solution to Equation ( 4) is said to be symmetric if it keep this symmetry, otherwise it is said to be asymmetric [58].Equation ( 4) can be numerically solved by the imaginary-time method [66,67].

Symmetric Vortex Solitons
By using input trial wavefunction: stationary symmetric vortex solitons (SVS) can be obtained from Equation ( 4), and their stability can be checked via real-time evolution.Typical results of the stable SVS are shown in Figure 1.Traditionally, only anisotropic fundamental solitons are found to be stable in various settings, and no stable anisotropic vortex soliton is reported in literature.However, numerical analysis reveals that the anisotropic SVSs are stable if the total norm N ≡ N + + N − = |φ + (x, y)| 2 dr + |φ − (x, y)| 2 dr exceeds a threshold N th ∼ 0.12.Under this threshold, both components of the SVS will become a breather mode showing oscillatory instability [58].This instability is different from the system with contact attractive nonlinearity [49] where the nonlinearity itself will lead to the collapse of the soliton, and the long-range DDI will not cause the collapse, such that the stability of the SVS under DDI is easier to archive.Stable families of SVSs are found, and they fulfill the Vakhitov-Kolokolov (VK) criterion dµ/dN < 0, which further sustain their stability.

Asymmetric Vortex Soliton
By using input trial wavefunction: the asymmetric vortex soliton can also be produced in this setting; where R ps = X in ps î + Y in ps ĵ is the initial offset of the vortex component with respect to the fundamental component and θ ps is the angular coordinate for the vortex's pivot placed at r = R ps .Typical results are shown in Figure 2. The vortical component of these asymmetric vortex solitons (AVS) apparently breaks the symmetry of Equations ( 5) and ( 6), which is not seen in the settings with contact nonlinearity [49].The AVS holds the same threshold N ∼ 0.12 as the SVS, below which the AVS will convert to robust breather mode, as well.For any fixed value of N > 0.12, there exists a continuous family of AVSs, parametrized by the offsets of the vortex's pivot, (X ps , Y ps ).The energy of the AVSs in each family is highly degenerated, depending only on N (see Figure 2(b1)).The pivot of the vortex component in both directions exhibits saturation behavior, and the max value of the pivot decreases with the total norm N (see Figure 2(b2,b3)) [58].These interesting results may originate from the anisotropic nature of DDI and open new possibilities for creating stable vortex 2D solitons in experiments, but no simple and intuitive explanation has yet been found.

Mobility of the SVS and AVS
Although the SOC breaks the Galilean invariance of the system, analysis of the SVS and AVS reveals that they can be both kicked to move by applying a initial kick to the solution: By tracking the shape change and trajectory of the center of mass (c.m.) of the solitons (mainly determined by the fundamental component) via real-time evolution, three kinds of motion of the SVS can be identified, as shown in Figure 3.

1.
For k smaller than a depinning threshold k d , the SVS will be pinned to move in a circle trajectory (see Figure 3a) and is isotropic in both directions x and y. 2.
For k > k d , the SVS is depinned and move with a negative effective mass (i.e., opposite the direction of the kick) in a spiral trajectory (see Figure 3b), but still isotropic in both directions.

3.
For k > k upper where k upper is another threshold larger than k d , the SVS will undergo a more complex "lacy" trajectory and become anisotropic (see Figure 3c).The AVS can also be kicked to move, but due to the anisotropy of the DDI and the asymmetric shape of the vortex component in that case, the motion pattern is far more complex than that of the SVS.The motion is highly anisotropic to the horizontal and vertical kick, and a kick in one direction will cause not only motion in this direction, but also in another direction.Under different kick parameters, the AVS can either exhibit pinning, spiral or "lacy" motions, but no easy pattern, as SVS is yet to be clarified [58].Typical examples of the motion of AVS are shown in Figure 4.These behaviors of mobility are important in the control and manipulation of solitons both in theory and experiments.

Semi-Vortices in Spinor SOC BEC with DDI and ZS
By adding a Zeeman splitting (ZS) term into Equation (4), the system now changes to spinor SOC BEC with DDI and ZS (Setting (ii)), with ZS strength Ω, the sign of which can be altered between positive and negative.Using the same trial input Equation (9) and the same numerical methods, stable SVS similar to that in Figure 1 can also be produced in this setting [59] (However, in [59], AVS is used to describe the soliton to emphasize that they are the anisotropic vortex soliton.Actually, they are the SVS we discussed in the previous section with no pivot of the vortex component.To avoid conflict, we stick to using SVS here.).Typical results of the SVSs are displayed in Figure 5.

Shape Control of the SVS by ZS
It is found that the strength of the ZS can change the shape of the SVS.As shown in Figure 5, when the strength of ZS Ω = 0, the system goes back to Equation (4), and similar SVS are produced (Panels (a1), (a2)).Just like the SVS in Figure 1, the anisotropic nature of DDI deforms the fundamental and vortex component in the horizontal direction (longitudinal directions compared to the dipole orientation), creating two density maxima separated in this direction.For this reason, these SVSs are called horizontal SVS [59].When Ω increases and exceeds a certain critical value Ω c , the shape of the SVS changes, especially for the vortex component (Panels (b1), (b2)), which now has two density maxima along the vertical directions (transverse directions compared to the dipole orientation; compare Panels (b2) and (a2)); therefore, it is called the vertical SVS.This change of shape is continuous with Ω increasing from zero, so at Ω = Ω c , the SVS features an elliptic ring with no density maxima for any specific direction (Panel c2).The value of Ω c can be determined numerically by defining a shape parameter: and the results are displayed in Figure 6.

FIG. 2. (a)
The shape parameter , defined as per Eq. ( 10), vs Figure 6 shows clearly that the strength of ZS, Ω, can be used as a shape control parameter to transform the SVS between horizontal and vertical mode.The critical value of this transition Ω c has an almost linear relation with the total norm N, which can be explained qualitatively by noticing that at large Ω and a large N limit, the energy of DDI and ZS scales as N 2 and ΩN, respectively [59].Since this transition is the result of competition between DDI and ZS, the transition then must occur at Ω = const.× N, as indicated in Panel (b) of Figure 6.

Inverted SVS with Ω < 0
The SVSs discussed before this subsection all have Ω > 0, with their fundamental (with vorticity S = 0) and vortical component (with vorticity S = 1) taking the lower and higher energy position, respectively, which can be regarded as the ground state of the system.If Ω is switched to < 0, then the fundamental and vortex component will switch their energy position, to be in the higher and lower energy state, respectively.This inverted state can be considered as the excited state of the system [59].Such inverted SVS is confirmed to exist and be stable via methods used in the previous sections, setting Ω < 0. Typical results are displayed in Figure 7.The value of M=∞ depends on the total norm, as well as of the size η of the initial kick.This dependence is illustrated in Fig. 5(b) by means of curves of M=∞(η) for two fixed values of N, which show that M=∞ increases with the growth of η.It is worthy to note that the essential dependence of the effective mass on the ZS strength occurs at quite small values of , therefore all the relevant AVSs are of the horizontal type [the transition to the vertical shape happens at much higher values of , see Fig. 2(b)].It follows from here that all the AVSs of the vertical type have the positive effective mass.This conclusion explains why the horizontal AVSs in Ref. [47] and vertical ones in Ref. [48] feature opposite signs of their dynamical mass, as mentioned above.Finally, the numerical results demonstrate that the effective mass approaches the norm of the soliton, i.e., M → N in the limit of → ∞.This conclusion is easily explained by the fact that, in the limit of → ∞, the AVS turns into a usual single-component soliton [see Eq. ( 14)], whose dynamical mass is always identical to N .

V. INVERTED ANISOTROPIC SEMIVORTICES, WITH < 0
As stressed above, we considered the soliton modes of semivortex type, with their fundamental and vortex constituents sitting, respectively, in lower-and higher-energy components, in terms of the ZS with > 0. Here, we aim to consider the reverse setting, with < 0, when the fundamental and vortex parts find themselves in the higherand lower-energy components, respectively.Because the fundamental constituents tends to carry a much larger share of the total norm than its vortical counterpart, one may expect that the inverted state will be an excited one, while the stable states considered above were ground states in the system.For this reason, it is quite interesting to explore the shape and, especially, stability of these excited states.Surprisingly, it is found below that they are stable in a considerable part of their parameter space.
Generic examples of the density and phase structure of the inverted AVSs are displayed in Fig. 6.As well as these examples, all the inverted AVSs feature the horizontal structure, as it is defined above (with two density maxima separated along the x axis, i.e., along the direction of the polarization of dipoles in the underlying condensate).Unlike the case of large positive considered above, vertical inverted AVSs have not been found.
The dependence of the AVS's chemical potential and energy on at fixed norm N is displayed in Fig. 7, both for the mode under the consideration and its flipped counterpart with vorticities −1 and 0 in components ψ+ and ψ−, respectively.In the absence of the ZS, = 0, both semivortex modes are fully equivalent to each other [40].Furthermore, the present mode at given and the flipped one at the ZS strength − remain equivalent, too.Thus, Fig. 7(b) makes it obvious that the present mode is an excited state at < 0, as its energy is higher than the energy of the ground state, which is realized, at < 0, by the flipped AVS mode.
Furthermore, Fig. 8(a) represents families of the AVS modes by displaying the share of the norm in the vortex component, F−, defined according to Eq. ( 5), vs at several fixed values of N. Naturally, F− increases with the growth of − .A nontrivial finding is that, although the inverted AVS is an excited state, as shown above, it has a vast stability region, > s(N ), with the dependence of the stability boundary on N, i.e., s(N ), shown in Fig. 8(c).In particular, Fig. 8(a) demonstrates that, at N > 2, the AVS stabilizes itself by shifting more than half of the total norm Families of stable inverted SVS can be identified through numerical results, which are all horizontal type (compare Figure 7 with Figures 1 and 5), and no vertical type inverted SVS is found to exist.Despite all be excited states of the system, as indicated in Figure 7c,d, they surprisingly occupy a vast stable region, given Ω > Ω s .Below, this Ωs inverted SVS still exists, but suffers spontaneous fragmentation instability (Panel (b)) [59].Previously, only unstable excited states of SVS were found [49,68], and no stable excited SVS has been reported in the literature, especially stable SVS with a dominant vortex component.These findings provide new possibilities to create and manipulate 2D excited vortex solitons experimentally.

Mobility and Effective Mass Controlled by ZS
Besides controlling the shape of the SVS, enabling transitions between the horizontal and vertical type, Zeeman splitting also exhibits crucial effects on the mobility of the SVS.By defining the effective mass M x,y = P x,y /V x,y where: and V x,y is the velocity of the SVS, the motion of the SVS can be studied quantitatively.Numerical results show that as in Section 3.3, the SVS can be kicked to move and undergoes different types of motion depending on the kick strength.The effective mass defined above only applies to the spiral motion (see Figure 3b), and it is found that for spiral motion, changing Ω can switch the effective mass between positive and negative, via infinite divergence at a specific value Ω = Ω M=∞ .Typical results are displayed in Figure 8.General conclusions can be drawn from these results: 1.
For Ω < Ω M=∞ , the SVS possesses a negative effective mass, and it will move in the opposite directions of the kick.This is in accordance with results in Section 3.3, where Ω = 0.

2.
For Ω = Ω M=∞ , the effective mass of SVS diverges to infinity, and the soliton will be pinned to move in a circular trajectory, similar to the motion shown in Panel (a) of Figure 3, which can be regarded as the immobility of the soliton.

3.
For Ω > Ω M=∞ , the SVS possesses a positive effective mass, and it will move along the directions of the kick.4.
For Ω → ∞, the effective mass approaches the total norm M → N, since it becomes a usual single-component soliton.

Gap Solitons in Spinor SOC BEC
Gap solitons (GS) are self-trapped modes existing in the spectral band gap of the system.The band gap generally is generated by certain types of spatially periodic potentials, such as Bragg gratings [69], wave guide arrays [21,70] and photonic crystals [17] in optics, or the optical lattice in BEC [23,24] and plasma medium [20,71].Combined with SOC, 1D and 2D BEC gap solitons can also be produced in the optical lattice [43,44,72,73].In these settings, the external potential, which creates the band structure, also exerts a spatial confinement on the solitons.
On the other hand, theoretical and numerical analysis reveals that by adding a Zeeman splitting, the SOC spinor BEC can support gap solitons in free space, with the collaboration of a DDI.Without spatial modulation, the energy gap is opened solely by the Zeeman splitting, which is a crucial ingredient to the SOC realization in experiments [74].Families of stable isotropic and anisotropic gap solitons exist at the edge of the band gap, depending on the orientation of the dipolar BEC [60].

Isotropic and Anisotropic Families of the Gap Soliton
By setting g = γ in Equation (1), the system now returns to a more general case with DDI, ZS and contact nonlinearity (Setting (iii)).Note that the sign of β can be altered by rotating the magnetic field [75], and the sign of g can be switched through Feshbach resonance.This equation actually describes a 2D spinor BEC tightly trapped in the transverse (i.e., third) direction with a trapping size a ⊥ , which requires the size of any soliton results to fulfill l a ⊥ .Meanwhile, the SOC strength needs to be strong enough to manifest effects experimentally, which requires β ≥ 1/a ⊥ .Together, these requirements give: and Equation (1) shows the ratio of kinetic energy and the SOC term is ∼ 1/(lβ), which implies under these settings that the kinetic term can be neglected safely.Equation (1) then reduces to: Linearization to Equation ( 14) with ψ ± ∼ exp(ipx + iqy − iwt) leads to the dispersion relation: which opens a band gap of width 2Ω in free space, as shown in Figure 9, where p, q are the quasi-momentum of the soliton in the x, y direction, respectively.Note that this energy gap will close if the kinetic terms are included [60].Analytic and numerical solutions to Equation ( 14) can be carried our [60], and the results show that without DDI (β = 0), the equation gives no soliton solutions when the contact nonlinearity is attractive (g < 0) and gives the unstable Townes' solitons when contact nonlinearity is repulsive (g > 0); thus, no stable solitons can be obtained under this setting without DDI.It needs to be noted that the above conclusions are drawn under conditions that the strengths of self-and cross-interactions are equal, g = γ in Equation (1).Actually, the relative ratio between g and γ can be adjusted by the Feshbach resonance technique experimentally, by applying a DC magnetic field to the spinor BEC [76][77][78], and very recent research revealed [79] that for unequal g and γ, even without DDI, a continuous family of stable gap solitons exists in the upper band provided that γ/g is less than a certain critical value (γ/g < 0.77 in [79]).
On the contrary, without contact nonlinearity (g = 0), but with DDI present (β > 0), the system can produce stable isotropic solitons near the bottom edge of the gap when the DDI takes the isotropic kernel (i.e., θ = 0 in Equation ( 3) where the dipoles are oriented perpendicular to the system's plane) and stable anisotropic solitons near the top edge of the gap when the DDI takes the anisotropic kernel (i.e., θ = 0 in Equation ( 3) where the dipoles are oriented parallel to the system's plane), as shown in Figure 9.With increasing total norm N, the isotropic gap solitons will extend as a continuous stable family across the full band into the upper Bloch band to become an embedded soliton (ES) [80,81], and it will terminate its stability as an ES when N exceeds a certain critical value N > N max (see Figure 9(b1)).The anisotropic gap solitons will also extend to a continuous family with increasing N, but it always terminates stability in the band (see Figure 9(b2)).To compare with experimental settings, an estimate gives the number of atoms constituting the gap soliton ranging from N ∼ 10 3 to N ∼ 10 4 , by using values of the magnetic moment of 52 Cr or 164 Dy atoms and the strength of transverse trapping potential ω ⊥ ∼ 100 Hz, which is feasible by current experimental techniques [60].

Mobility of the Gap Soliton
Both isotropic and anisotropic gap solitons can be kicked to move.For isotropic GS, it effective mass is always negative, in accordance with the generic features of GS [20,71].However, for anisotropic GS, its effective mass is positive.Typical results are shown in Figure 10.For both types of GS, the effective mass varies almost linearly with the total norm N, and at small N, the effective mass is isotropic along the x and y directions, but exhibits anisotropy with increasing of N.

Effects of Contact Nonlinearity
The conclusions in Section 5.1 are drawn without contact nonlinearity, i.e., for g = 0 in Equation (14).The effects of contact nonlinearity can be investigated numerically.The system is found to have an intrinsic stable boundary g cr ±0.25 for isotropic and anisotropic GS, respectively.Thus, the isotropic GSs are stable if g < 0.25, and the anisotropic GSs are stable if g > −0.25 (recall that g < 0 and g > 0 correspond to repulsive contact linearity for isotropic GS and anisotropic GS, respectively), which means the GSs are completely stable when arbitrarily strong contact nonlinearity is repulsive.On the other hand, due to the critical collapse effect of the contact nonlinearity, the norm of isotropic and anisotropic solitons cannot exceed certain limits, which is determined by the norm of the Townes' soliton N T 5.85, resulting in an exist condition for isotropic and anisotropic GS: g aniso > − N T /(ΩN) (16) Between these two limits, that is 0.25 < g iso < N T /(ΩN) and −N T /(ΩN) < g aniso < −0.25, the GSs are unstable and spontaneously transform into breathers in propagation [60].

Conclusions
We discussed in this review several settings of creating stable vortex 2D solitons in free space.By SOC combined with attractive anisotropic DDI, the continuous family of symmetric (SVS) and asymmetric (AVS) vortex soliton can be produced in spinor BEC.By adding a Zeeman splitting (ZS), the shape of the SVS can be switched between the horizontal and vertical type, the effective mass of the SVS can be altered from positive to negative and cross a divergence point where the effective mass goes to infinity and the SVS can be effectively pinned to be immobile.Inverted SVS is found to exist stably with the sign of ZS altered, despite the fact that it is actually the excited state of the system, which is new to the understanding of excited vortex solitons.The ZS also generates a spectral band structure in the system, which supports the continuous family of gap solitons (GS) in free space.Extensions to these research works is interesting to consider, especially on how to embed these findings into 3D settings, where the collapse effects are far stronger than in 2D, and even with the help of long-range DDI, d-wave collapse can easily destroy the soliton.Recent studies point out that combining SOC, DDI with the Lee-Huang-Yang (LHY) correction effect can suppress the collapsing or expanding of atoms, creating a stable quantum droplet in 3D free space [82].Thus, extending the research works in this review to 3D settings under LHY corrections is an intriguing topic to investigate.

Figure 1 .
Figure 1.Stable symmetric vortex solitons (SVS) numerically obtained from Equation (4).(a) Stable SVS with total norm N = 0.15.(a1,a2) present the density distribution of the fundamental and vortical component, |φ + (x, y)| 2 and |φ − (x, y)| 2 , respectively.(a3) is the phase structure of the vortical component.(a4) is the direct run of this SVS.(b) The chemical potential µ versus N. The yellow area indicates where the SVS is subjected to oscillatory instability and transforms into breather mode [58].As shown in Figure 1, all the semivortices (half vortices) have one component (φ + ) carrying vorticity S = 0 (called the fundamental component), and the other component (φ − ) carries vorticity S = +1 (called the vortical component).What is more, all the stable SVs feature anisotropic patterns due to the anisotropic nature of the DDI, and at the same time, keep the symmetry in Equation (7).Traditionally, only anisotropic fundamental solitons are found to be stable in various settings, and no stable anisotropic vortex soliton is reported in literature.However, numerical analysis reveals that the anisotropic SVSs are stable if the total norm N ≡ N + + N − = |φ + (x, y)| 2 dr + |φ − (x, y)| 2 dr exceeds a threshold N th ∼ 0.12.Under this threshold, both components of the SVS will become a breather mode showing oscillatory instability [58].This instability is different from the system with contact attractive nonlinearity[49] where the nonlinearity itself will lead to the collapse of the soliton, and the long-range DDI will not cause the collapse, such that the stability of the SVS under DDI is

Figure 2 .
Figure 2. Stable asymmetric vortex solitons (AVS) numerically obtained from Equation (4).(a) Stable AVS with total norm N = 0.15.The first column presents the density distribution of the fundamental and vortical component, |φ + )| 2 and |φ − | 2 , and the phase structure of φ − , with (X ps , Y ps ) = (0.04, 0).The second column is the AVS with (X ps , Y ps ) = (0, 0.065).(b) Various relations between the input and output X ps with energy and the total norm of the AVS [58].

Figure 6 .
Figure 6.Shape change of the SVS controlled by the Zeeman splitting Ω.(a) Shape parameter ∆ versus Ω with different N; the dashed line determines Ω c .(b) (Almost) linear relationship between Ω c and N [59].

FIG. 6 .
FIG. 6.The same as in Fig. 1, but for the inverted AVSs, with < 0. The first row represents a stable AVS with (N, ) = (6,−1.82),with equal norms in the fundamental and vortex components.The second row shows the soliton with (N, ) = (10,−10.8),which is a dipole mode for φ− and is located at the stability boundary [see Fig. 8(a) below].

Figure 9 .
Figure 9. Band gap and gap solitons of spin-orbit-coupling (SOC) spinor BEC with dipole-dipole interaction (DDI), Zeeman splitting (ZS) and contact nonlinearity.(a) The band gap structure described by Equation (15) with (β, Ω) = (1, 10).(b) Isotropic and anisotropic 2D gap solitons.The yellow area is the band; the white area is the gap.The blue curve is the branch of gap solitons found numerically, divided by the blue dot into stable (solid) and unstable (dashed) segments.The red dashed curve is the analytical prediction based on Equation (14)[60].GS, gap solitons.

Figure 10 .
Figure 10.Effective mass for isotropic (left) and anisotropic (right) 2D GS under kick in the x and y directions [60].