Ji-Guo Wang^,1,2, Yue-Qing Li^1,21Department of Mathematics and Physics, Shijiazhuang TieDao University, Shijiazhuang 050043, China
²Institute of Applied Physics, Shijiazhuang TieDao University, Shijiazhuang 050043, China

Received:2020-01-15Revised:2020-05-2Accepted:2020-05-15Online:2020-08-04


Abstract
We study the ground-state phases of two-dimensional rotating spin–orbit coupled spin-1/2 Bose–Einstein condensates (BECs) in a gradient magnetic field. The competition between gradient magnetic field, spin–orbit coupling and rotation leads to a variety of ground-state phase structures. In the weakly rotation regime, as the increase of gradient magnetic field strength, the BECs experiences a phase transition from the unstable phase to the single vortex-line phase. The unstable phase presents the vortex lines structures along the off-diagonal direction. With magnetic field gradient strength increasing, the number of vortex lines changes accordingly. As the magnetic field gradient strength increases further, the single vortex-line phase with a single vortex line along the diagonal direction is formed. The phase diagram shows that the boundary between the two phases is linear with the relative repulsion λ≥1 and is nonlinear with λ<1. In the relatively strong rotation regime, in addition to the unstable phase and the single vortex-line phase, the vortex-ring phase is formed for the strong magnetic field gradient and rapid rotation. The vortex-ring phase shows the giant and hidden vortex structures at the center of ring. The strong magnetic field gradient makes the number of the vortices around the ring unchanged.
Keywords： Bose–Einstein condensate;spin–orbit coupling;gradient magnetic field


PDF (1371KB)MetadataMetricsRelated articlesExportEndNote|Ris|BibtexFavorite
Cite this article
Ji-Guo Wang, Yue-Qing Li. Rotating spin-1/2 Bose–Einstein condensates in a gradient magnetic field with spin–orbit coupling. Communications in Theoretical Physics, 2020, 72(9): 095701- doi:10.1088/1572-9494/ab95fb

1. Introduction

With the first experimental observation in 1995 [1, 2], the Bose–Einstein condensate (BEC) has been widely studied in both experiment and theory. The collective excitation of matter waves in BECs, such as the bright and dark solitons [3–8], vortex lattices [9–11] and rogue waves [12–14] draw a great deal of interest. The experimental realization of spin–orbit coupling (SOC) in BECs by the NIST group [15] arises a hot research on spin–orbit coupled spin-1/2 BECs. In spin–orbit coupled spin-1/2 BECs, the competition between the intraspecies repulsion and the interspecies repulsion leads to a variety of ground-state phases, such as the plane wave phase, stripe phase and meron with dipole–dipole repulsion [16–20]. The dynamic properties of spin–orbit coupled spin-1/2 BECs are also studied, such as solitons, quantized vortices dynamics and Josephson dynamics [21–29].

In the rotating spin–orbit coupled spin-1/2 BECs, a variety of vortex structures [30–33] are found, such as half-quantum vortex, ring-like structure with domain, giant vortex and circular-hyperbolic skyrmion. The Mermin-Ho vortex, quantum entanglement and quantum knot are observed in the spinor BECs with gradient magnetic field [34–36]. In [37], a tunable SOC is realized in spin-1 BECs by modulating gradient magnetic field in experiment. The SOC strength can be adjusted by changing the momentum impulse from the gradient magnetic field. Recently, the related researches of spin–orbit coupled spin-1/2 and spin-1 BECs in a harmonic trap with gradient magnetic field have been studied [38–41]. In spin-1/2 BECs, the gradient magnetic field and SOC can control the transformation between the single plane-wave phase and half-skyrmion phase. In spin-1 BECs, the influence of SOC and gradient magnetic field on the properties of the monopoles with polar-core vortex. Compared with the previous works [42–45], the monopoles have a long-lived time, therefore, the experimental observation of the monopoles becomes easier.

In this paper, we investigate the ground-state phases of two-dimensional (2D) rotating spin–orbit coupled spin-1/2 BECs loaded in a harmonic plus quartic trap with gradient magnetic field. The competition between gradient magnetic field, SOC and rotation leads to a variety of ground-state phase structures. In the weakly rotation regime, as the increase of gradient magnetic field strength, two types of ground-state phases are found with a fixed SOC, i.e. the unstable phase and the single vortex-line phase. The unstable phase presents the vortex lines structures along the off-diagonal direction. With magnetic field gradient strength increases, the number of vortex lines changes accordingly. As the magnetic field gradient strength increases further, a single vortex line along the diagonal direction, the single vortex-line phase is found. The phase diagram shows that the boundary between the two phases is linear with the relative repulsion λ≥1 and is nonlinear with λ<1. In the relatively strong rotation regime, in addition to the unstable phase and the single vortex-line phase, the vortex-ring phase is formed for the strong magnetic field gradient and rapid rotation. The vortex-ring phase shows the giant and hidden vortex structures at the center of ring, the strong magnetic field gradient makes the number of the vortices around the ring unchanged (n=6).

2. Model and Hamiltonian

In our paper, we study the ground-state phases of 2D rotating spin–orbit coupled spin-1/2 BECs loaded in a harmonic plus quartic trap with gradient magnetic field, the expectation value of the Hamiltonian are given as(1)$ \begin{eqnarray}\begin{array}{rcl}E[{\rm{\Psi }}] & \equiv & \langle \widehat{H}\rangle =\displaystyle \int {{\rm{d}}}^{2}{\boldsymbol{r}}{{\rm{\Psi }}}^{\dagger }\left[-\displaystyle \frac{{{\hslash }}^{2}}{2m}{{\rm{\nabla }}}^{2}+V({\boldsymbol{r}})+{\upsilon }_{\mathrm{soc}}-{\rm{\Omega }}{L}_{z}\right.\\ & & \left.+{g}_{F}{\mu }_{B}{\boldsymbol{B}}({\boldsymbol{r}})\cdot {\boldsymbol{\sigma }}\Space{0ex}{3.08ex}{0ex}\right]{\rm{\Psi }}+\displaystyle \int {{\rm{d}}}^{2}{\boldsymbol{r}}({c}_{1}{N}_{\uparrow }^{2}+{c}_{2}{N}_{\downarrow }^{2}\\ & & +{c}_{12}{N}_{\uparrow }{N}_{\downarrow }),\end{array}\end{eqnarray}$where ${\rm{\Psi }}=({{\rm{\Psi }}}_{\uparrow },{{\rm{\Psi }}}_{\downarrow })$ is the wave function of spin-1/2 BECs, m is the atomic mass and ${\boldsymbol{r}}=(x,y)$. The densities of spin-$\uparrow $ component ${N}_{\uparrow }=| {{\rm{\Psi }}}_{\uparrow }{| }^{2}$ and spin-$\downarrow $ component ${N}_{\downarrow }=| {{\rm{\Psi }}}_{\downarrow }{| }^{2}$. The repulsions ${c}_{1}=\tfrac{4\pi {{\hslash }}^{2}{a}_{1}}{m}$ and ${c}_{2}=\tfrac{4\pi {{\hslash }}^{2}{a}_{2}}{m}$ are the intraspecies repulsions, ${c}_{12}=\tfrac{4\pi {{\hslash }}^{2}{a}_{12}}{m}$ is the interspecies repulsion. ${\boldsymbol{\sigma }}$ is the spin density ${\boldsymbol{\sigma }}=({\sigma }_{x},{\sigma }_{y},{\sigma }_{z})=({{\rm{\Psi }}}^{\dagger }{\hat{\sigma }}_{x}{\rm{\Psi }},{{\rm{\Psi }}}^{\dagger }{\hat{\sigma }}_{y}{\rm{\Psi }},{{\rm{\Psi }}}^{\dagger }{\hat{\sigma }}_{z}{\rm{\Psi }})$ with the 2×2 spin-1/2 matrices $\widehat{\sigma }$. The isotropic Rashba SOC is ${\upsilon }_{\mathrm{soc}}\,=-{\rm{i}}{\hslash }\gamma ({\sigma }_{x}{\partial }_{y}-{\sigma }_{y}{\partial }_{x})$ with SOC strength γ. $V({\boldsymbol{r}})=\tfrac{1}{2}m{\omega }_{\perp }^{2}({{\boldsymbol{r}}}^{2}+\beta {{\boldsymbol{r}}}^{4})$ is a harmonic plus quartic trap, here we set β=0.5. The rotational frequency ω along the z direction with the orbit angular momentum ${L}_{z}=-{\rm{i}}{\hslash }(x{\partial }_{y}-y{\partial }_{x})$. The parameters g_F is Lande factor and μ_B is Bohr magnetic moment. The gradient magnetic field [46] ${\boldsymbol{B}}({\boldsymbol{r}})=B(x{{\rm{e}}}_{x}-y{{\rm{e}}}_{y})$, where B is the strength. g_F=−1/2 is Lander factor and μ_B is Bohr magnetic moment.

The time evolution of the mean field is governed by(2)$ \begin{eqnarray}{\rm{i}}{\hslash }\displaystyle \frac{\partial {\rm{\Psi }}}{\partial t}=\displaystyle \frac{\delta E}{\delta {{\rm{\Psi }}}^{* }}.\end{eqnarray}$The time-dependent Gross–Pitaevskii equations of spin–orbit coupled BECs can be obtained by substituting equation (1) into (2),(3)$ \begin{eqnarray}\begin{array}{rcl}{\rm{i}}\displaystyle \frac{\partial {{\rm{\Psi }}}_{\uparrow }}{\partial t} & = & \left[-\displaystyle \frac{{{\rm{\nabla }}}^{2}}{2}+V({\boldsymbol{r}})+{\rm{i}}{\rm{\Omega }}(x{\partial }_{y}-y{\partial }_{x})+{c}_{1}| {{\rm{\Psi }}}_{\uparrow }{| }^{2}\right.\\ & & \left.+{c}_{12}| {{\rm{\Psi }}}_{\downarrow }{| }^{2}]{{\rm{\Psi }}}_{\uparrow }+[B(x+{\rm{i}}y)+\gamma (-{\rm{i}}{\partial }_{y}+{\partial }_{x})\Space{0ex}{3.15ex}{0ex}\right]{{\rm{\Psi }}}_{\downarrow },\\ {\rm{i}}\displaystyle \frac{\partial {{\rm{\Psi }}}_{\downarrow }}{\partial t} & = & \left[\Space{0ex}{3.08ex}{0ex}(-\displaystyle \frac{{{\rm{\nabla }}}^{2}}{2}+V({\boldsymbol{r}})+{\rm{i}}{\rm{\Omega }}(x{\partial }_{y}-y{\partial }_{x})\Space{0ex}{3.08ex}{0ex})+{c}_{2}| {{\rm{\Psi }}}_{\downarrow }{| }^{2}\right.\\ & & \left.+{c}_{12}| {{\rm{\Psi }}}_{\uparrow }{| }^{2}]{{\rm{\Psi }}}_{\downarrow }+[B(x-{\rm{i}}y)+\gamma (-{\rm{i}}{\partial }_{y}-{\partial }_{x})\Space{0ex}{3.08ex}{0ex}\right]{{\rm{\Psi }}}_{\uparrow }.\end{array}\end{eqnarray}$

The length (x and y), energy (SOC, rotation and interaction) and magnetic field gradient are measured in units of $\sqrt{{\hslash }/m{\omega }_{\perp }}$, ℏω_⊥ and ${\hslash }{\omega }_{\perp }/({g}_{F}{\mu }_{B}{a}_{h})$, respectively. Where ${a}_{h}=\sqrt{{\hslash }/m{\omega }_{\perp }}$ is characteristic length of harmonic plus quartic trap. We assume that the intraspecies repulsions c₁=c₂=c and the relative repulsion λ=c₁₂/c in the present paper.

3. Ground-state phases

The ground-state phases of spin-1/2 BECs with isotropic SOC loaded in a harmonic plus quartic trap with gradient magnetic field can be obtained by using the time-splitting Fourier spectral method with the imaginary-time propagation [47, 48], it makes the substitution $t\to {\rm{i}}t$ in equation (3), which leads to an exponential decay of the wave function, and a corresponding decay of the eigenstates via. The eigenenergy governs the decay rate, and so the eigenstate with the lowest energy, i.e. the ground state of the system. The spatial and time steps employed in this paper are Δx(y)=0.05 and ${\rm{\Delta }}t=0.0001$. The interspecies repulsion c=1000 and the relative repulsion λ=1.2 are fixed.

We first study on the ground-state phases of the relatively simple case without SOC (γ=0). Figure 1 shows the ground-state densities of spin–orbit coupled spin-1/2 BECs with the rotational frequencies ω=1.0 in (a)–(b) and ω=1.4 in (c)–(d), respectively. In figure 1(a), the rotation frequency is too small to generate vortex in the harmonic plus quartic trap with magnetic field gradient strength B=0. The ground-state phase is a plane wave phase, which is similar with the case without rotation. However, the ground-state phase becomes a stable phase with gradient magnetic field B=10 in figure 1(b). The strength of the center of the gradient magnetic field is zero, which arises the singularity of the topological defect. The densities show a discrete soliton in spin-$\uparrow $ component and a Mermin-Ho vortex in spin-$\downarrow $ component at the center position of trap. As the rotational frequency is increased to ω=1.4, the vortex lines are shown and they are symmetric about the diagonal direction without magnetic field gradient strength B=0 in figure 1(c). We investigate the effects of magnetic field gradient on the vortices. The vortices realize the vortex rings with magnetic field gradient strength B=10 in figure 1(d). A variety of shapes can be formed by the vortices in the inner ring as the magnetic field gradient increases gradually, such as triangular, square and so on.

Figure 1.

New window|Download| PPT slide
Figure 1.The ground-state density profiles of rotating spin-1/2 BECs loaded in a harmonic plus quartic trap with gradient magnetic field. The columns in every panel from left to right are the densities of spin-$\uparrow $ component, spin-$\downarrow $ component, total and difference. The upper row with rotation strength ω = 1, the magnetic field gradient strengths B=0 in (a) and B=10 in (b). The lower row with rotation strength ω = 1.4, the magnetic field gradient strengths B=0 in (c) and B=10 in (d).


When considering SOC, the competition between SOC, rotation and gradient magnetic field leads to the bountiful ground-state phases. Now we study the effects of gradient magnetic field B on the ground-state phases and fix the SOC strength γ=2. The ground-state phases are shown in figures 2 and 3 with rotational frequencies ω=1.0 and ω=1.4, respectively. We find that the vortices arranged radially around the center of the circle without magnetic field gradient B=0 in figure 2(a), which realize the vortex lattices. As the magnetic field gradient strength increases, the vortex lattices structure evolved into the vortex lines. The vortex lines along the oblique direction, as shown in figures 2(b)–(c). The spin-flip of SOC can generate the vortices, however, the spin is not in the plane. The gradient magnetic field can flip the atomic magnetic moment in the non-plane, it will drive the spin along direction of the gradient magnetic field. The vortices are distributed uniformly along the off-diagonal and diagonal directions. However, the gradient magnetic field inhibits the generation of vortices, which results in the number of vortex lines is unstable, we take the ground-state phase as the unstable phase. As the magnetic field gradient strength increases further, a single vortex line along the diagonal direction is shown in figure 2(d). Though the number of vortices of vortex line decreases, the shape of the ground-state phase keeps unchanged as the gradient magnetic field increases. We take the ground-state phase as the single vortex-line phase.

Figure 2.

New window|Download| PPT slide
Figure 2.The ground-state density profiles of rotating spin–orbit coupled spin-1/2 BECs loaded in a harmonic plus quartic trap with gradient magnetic field. The SOC strength γ=2 and rotation strength ω=1. The columns in every panel from the left to right are the densities of spin-$\uparrow $ component, spin-$\downarrow $ component, total and difference. The magnetic field gradient strengths in (a)–(d) are: B=0, 1, 2, and 4, respectively.

Figure 3.

New window|Download| PPT slide
Figure 3.The ground-state density profiles of spin–orbit coupled spin-1/2 BECs in a harmonic plus quartic trap with gradient magnetic field. The SOC strength γ=2 and rotation strength ω = 1.4. The columns in every panel from the left to right are the densities of spin-$\uparrow $ component, spin-$\downarrow $ component, total and difference. The magnetic field gradient strengths in (a)–(f) are: B=0, 4, 7, 7.5, 11 and 30, respectively.


For the large rotation strength ω=1.4, the vortices of the ground-state phase arranged radially around the center of the circle without magnetic field gradient B=0 in figure 3(a). The unstable phases are also shown in figures 3(b)–(e). In the unstable phases regime, with the magnetic field gradient increasing, we find that the number of the vortices along the off-diagonal direction is a constant value in figures 3(c)–(e), which is different from the shown in figure 2. As the magnetic field gradient strength increases further, the strong magnetism of the gradient magnetic field leads to the atoms hard to condensed in the center of the harmonic plus quartic trap. The vortex line evolves into the vortex ring structure. The giant and hidden vortex structures are found at the center of ring. The hidden vortex is only reported in the single-component BEC in [49, 50]. We take the ground-state phase with the vortex ring structure as the vortex-ring phase, as shown in figure 3(f). The strong magnetic field gradient makes the number of the vortices around the ring is a constant value (n=6).

The phase distributions of the ground-state phase are studied in figure 4. For the small rotation, we can see that the phase distributions from the general vortex structure to the vortex alignment with increasing the magnetic field gradient in figures 4(a)–(b). In figures 4(c)–(d), it is found that the hidden vortex structure at the center of ring with the large rotation. Figure 5 plots the corresponding spin textures of the ground-state phases in figure 4. The spin texture is ${S}_{x}=({{\rm{\Psi }}}_{\uparrow }^{* }{{\rm{\Psi }}}_{\downarrow }+{{\rm{\Psi }}}_{\downarrow }^{* }{{\rm{\Psi }}}_{\uparrow })/| {\rm{\Psi }}{| }^{2}$, ${S}_{y}=-{\rm{i}}({{\rm{\Psi }}}_{\uparrow }^{* }{{\rm{\Psi }}}_{\downarrow }-{{\rm{\Psi }}}_{\downarrow }^{* }{{\rm{\Psi }}}_{\uparrow })/| {\rm{\Psi }}{| }^{2}$, ${S}_{z}\,=(| {{\rm{\Psi }}}_{\uparrow }{| }^{2}-| {{\rm{\Psi }}}_{\downarrow }{| }^{2})/| {\rm{\Psi }}{| }^{2}$. The Skymiron line structures are shown in the spin textures of figures 5(a)–(d).

Figure 4.

New window|Download| PPT slide
Figure 4.Panels (a)–(b) are the phase distributions of the ground-state phases in figures 2(b), (d) and (c)–(d) are the phase distributions of the ground-state phases in figures 3(c), (f), respectively.

Figure 5.

New window|Download| PPT slide
Figure 5.Panels (a)–(d) are the corresponding spin texture of the ground-state phases in figure 4. The color of the arrows indicates the value of S_z.

4. Phase diagrams

In order to observe the ground-state phases in experiment easily, we study the phase diagrams of rotating spin–orbit coupled spin-1/2 BECs in a gradient magnetic field. The ground-state phase diagram is spanned by the SOC strength γ and magnetic field gradient strength B with the different relative repulsions λ=1.2 and λ=0.6, one can be seen in figure 6. The rotation strength ω=1. As the increase of gradient magnetic field strength, the BECs experiences a phase transition from the unstable phase to the single vortex-line phase. The unstable phase presents the vortex lines structures along the off-diagonal direction. With magnetic field gradient strength increases, the number of vortex lines changes accordingly (see figures 1(b)–(c)). The single vortex-line phase has a single vortex line along the diagonal direction. The shape of the single vortex-line phase keeps unchanged as the gradient magnetic field increases (see figure 1(d)). The phase diagram shows that the boundaries between the two phases is linear with the relative repulsion λ≥1 (see black dotted line λ=1.2) and is nonlinear with λ<1 (see red dotted line λ=0.6). For the fixed gradient magnetic field, the time reversal symmetry of spin–orbit coupled spin-1/2 BECs is preserved when the relative repulsion λ≥1. The relative repulsion inhibits the formation of the single vortex-line phase.

Figure 6.

New window|Download| PPT slide
Figure 6.The ground-state phase diagram is spanned by the SOC strength γ and magnetic field gradient strength B with the different relative repulsions λ=1.2 and λ=0.6. The rotation strength ω=1. Two types of ground-state phases are shown.


We also discuss the ground-state phase diagram is spanned by rotation strength ω and magnetic field gradient strength B in figure 7. The SOC strength γ=2 and the relative repulsion λ=1.2. Three types of the ground-state phases are found. In addition to the unstable phase and the single vortex-line phase, the vortex-ring phase is formed for the strong magnetic field gradient and rapid rotation (ω≥1.3). The vortex-ring phase shows the giant and hidden vortex structures at the center of ring. The strong magnetic field gradient makes the number of the vortices around the ring unchanged (see figure 3(f)).

Figure 7.

New window|Download| PPT slide
Figure 7.The ground-state phase diagram is spanned by the rotation strength ω and magnetic field gradient strength B with the SOC strength γ=2. Three types of ground-state phases are shown.

5. Summary

In summary, we have investigated the ground-state phases of 2D rotating spin–orbit coupled spin-1/2 BECs loaded in a harmonic plus quartic trap with gradient magnetic field. The competition between gradient magnetic field, SOC and rotation leads to a variety of ground-state phase structures. In the weakly rotation regime, as the increase of gradient magnetic field strength, two types of ground-state phases are found with a fixed SOC, i.e. the unstable phase and the single vortex-line phase. The unstable phase presents the vortex lines structures along the off-diagonal direction. With magnetic field gradient strength increases, the number of vortex lines changes accordingly. As the magnetic field gradient strength increases further, a single vortex line along the diagonal direction, the single vortex-line phase is found. The phase diagram shows that the boundary between the two phases is linear with the relative repulsion λ≥1 and is nonlinear with λ<1. In the relatively strong rotation regime, in addition to the unstable phase and the single vortex-line phase, the vortex-ring phase is formed for the strong magnetic field gradient and rapid rotation. The vortex-ring phase shows the giant and hidden vortex structures at the center of ring. The strong magnetic field gradient makes the number of the vortices around the ring unchanged (n=6).

Acknowledgments

This work is supported by the NSF of China under Grant No. 11 904 242, and the NSF of Hebei province under Grant No. A2019210280.

Reference View Option By original order
By published year
By cited within times
By Impact factor

[1]Anderson M H Ensher J R Matthews M R Wieman C E Cornell E A 1995 Science 269 198
DOI:10.1126/science.269.5221.198 [Cited within: 1]

[2]Davis K B Mewes M O Andrews M R van Druten N J Durfee D S Kurn D M Ketterle W 1995 Phys. Rev. Lett. 75 3969
DOI:10.1103/PhysRevLett.75.3969 [Cited within: 1]

[3]Khaykovich L Schreck F Ferrari G Bourdel T Cubizolles J Carr L D Castin Y Salomon C 2002 Science 296 1290
DOI:10.1126/science.1071021 [Cited within: 1]

[4]Strecker K E Partridge G B Truscott A G Hulet R G 2002 Nature 417 150
DOI:10.1038/nature747

[5]Li L Li Z D Malomed B A Mihalache D Liu W M 2005 Phys. Rev. A 72 033611
DOI:10.1103/PhysRevA.72.033611

[6]Liang Z X Zhang Z D Liu W M 2005 Phys. Rev. Lett. 94 050402
DOI:10.1103/PhysRevLett.94.050402

[7]Wang D S Hu X H Hu J P Liu W M 2010 Phys. Rev. A 81 025604
DOI:10.1103/PhysRevA.81.025604

[8]Hamner C Zhang Y Chang J J Zhang C Engels P 2013 Phys. Rev. Lett. 111 264101
DOI:10.1103/PhysRevLett.111.264101 [Cited within: 1]

[9]Ji A C Liu W M Song J L Zhou F 2008 Phys. Rev. Lett. 101 010402
DOI:10.1103/PhysRevLett.101.010402 [Cited within: 1]

[10]Kuopanportti P Huhtamäki J A M Möttönen M 2012 Phys. Rev. A 85 043613
DOI:10.1103/PhysRevA.85.043613

[11]Kasamatsu K Eto M Nitta M 2016 Phys. Rev. A 93 013615
DOI:10.1103/PhysRevA.93.013615 [Cited within: 1]

[12]Wen L Li L Li Z D Song S W Zhang X F Liu W M 2011 Eur. Phys. J. D 64 473 478
DOI:10.1140/epjd/e2011-20485-4 [Cited within: 1]

[13]Zhao L C Liu J 2013 Phys. Rev. E 87 013201
DOI:10.1103/PhysRevE.87.013201

[14]Zhong W P Belic M Malomed B A 2015 Phys. Rev. E 92 053201
DOI:10.1103/PhysRevE.92.053201 [Cited within: 1]

[15]Lin Y J Jiménez-García K Spielman I B 2011 Nature 471 83
DOI:10.1038/nature09887 [Cited within: 1]

[16]Wang C Gao C Jian C M Zhai H 2010 Phys. Rev. Lett. 105 160403
DOI:10.1103/PhysRevLett.105.160403 [Cited within: 1]

[17]Wu C Mondragon-Shem I 2011 Chin. Phys. Lett. 28 097102
DOI:10.1088/0256-307X/28/9/097102

[18]Li Y Pitaevskii L P Stringari S 2012 Phys. Rev. Lett. 108 225301
DOI:10.1103/PhysRevLett.108.225301

[19]Chen X Rabinovic M Anderson B M Santos L 2014 Phys. Rev. A 90 043632
DOI:10.1103/PhysRevA.90.043632

[20]Wilson R M Anderson B M Clark C W 2013 Phys. Rev. Lett. 111 185303
DOI:10.1103/PhysRevLett.111.185303 [Cited within: 1]

[21]Achilleos V Frantzeskakis D J Kevrekidis P G Pelinovsky D E 2013 Phys. Rev. Lett. 110 264101
DOI:10.1103/PhysRevLett.110.264101 [Cited within: 1]

[22]Chiu T L Liu T Y Chan H N Chow K W 2017 Commun. Theor. Phys. 68 290
DOI:10.1088/0253-6102/68/3/290

[23]Kasamatsu K 2015 Phys. Rev. A 92 063608
DOI:10.1103/PhysRevA.92.063608

[24]Zhang D W Fu L B Wang Z D Zhu S L 2012 Phys. Rev. A 85 043609
DOI:10.1103/PhysRevA.85.043609

[25]Chai X D Yu Z F Zhang A X Xue J K 2017 Chin. Phys. Lett. 34 090301
DOI:10.1088/0256-307X/34/9/090301

[26]Xu Y Zhang Y Wu B 2013 Phys. Rev. A 87 013614
DOI:10.1103/PhysRevA.87.013614

[27]Sakaguchi H Li B Malomed B A 2014 Phys. Rev. E 89 032920
DOI:10.1103/PhysRevE.89.032920

[28]Wang X B Han B 2019 Commun. Theor. Phys. 71 152
DOI:10.1088/0253-6102/71/2/152

[29]Wu C J Ian M S Zhou X F 2011 Chin. Phys. Lett. 28 097102
DOI:10.1088/0256-307X/28/9/097102 [Cited within: 1]

[30]Xu X Q Han J H 2011 Phys. Rev. Lett. 107 200401
DOI:10.1103/PhysRevLett.107.200401 [Cited within: 1]

[31]Zhou X F Zhou J Wu C 2011 Phys. Rev. A 84 063624
DOI:10.1103/PhysRevA.84.063624

[32]Liu C F Fan H Zhang Y C Wang D S Liu W M 2012 Phys. Rev. A 86 053616
DOI:10.1103/PhysRevA.86.053616

[33]Ramachandhran B Opanchuk B Liu X J Pu H Drummond P D Hu H 2012 Phys. Rev. A 85 023606
DOI:10.1103/PhysRevA.85.023606 [Cited within: 1]

[34]Ray M W Ruokokoski E Tiurev K Mottonen M Hall D S 2015 Science 348 544
DOI:10.1126/science.1258289 [Cited within: 1]

[35]Hall D S Ray M W Tiurev K Ruokokoski E Gheorghe A H Mottonen M 2016 Nat. Phys. 12 478
DOI:10.1038/nphys3624

[36]Liu J S Li J Liu W M 2017 Acta Phys. Sin. 66 130305(in Chinese)
DOI:10.7498/aps.66.130305 [Cited within: 1]

[37]Luo X Y Wu L N Chen J Y Guan Q Gao K Y Xu Z F You L Wang R Q 2016 Sci. Rep. 6 18983
DOI:10.1038/srep18983 [Cited within: 1]

[38]Li J Yu Y M Zhuang L Liu W M 2017 Phys. Rev. A 95 043633
DOI:10.1103/PhysRevA.95.043633 [Cited within: 1]

[39]Li J Liu W M 2018 Acta Phys. Sin. 67 110302(in Chinese)
DOI:10.7498/aps.67.20180539

[40]Wang J G Yang S J 2018 J. Phys.: Condens. Matter 30 295404
DOI:10.1088/1361-648X/aacc42

[41]Jin J J Guo H Zhang S Y Yan S B 2019 Ann. Phys. 411 167953
DOI:10.1016/j.aop.2019.167953 [Cited within: 1]

[42]Ruokokoski E Pietilä V Möttönen M 2011 Phys. Rev. A 84 063627
DOI:10.1103/PhysRevA.84.063627 [Cited within: 1]

[43]Ray M W Ruokokoski E Kandel S Möttönen M Hall1 D S 2014 Nature 505 657
DOI:10.1038/nature12954

[44]Cho Y M 2001 Phys. Rev. Lett. 87 252001
DOI:10.1103/PhysRevLett.87.252001

[45]Kiffner M Li W H Jaksch D 2013 Phys. Rev. Lett. 110 170402
DOI:10.1103/PhysRevLett.110.170402 [Cited within: 1]

[46]Pritchard D E 1983 Phys. Rev. Lett. 51 1336
DOI:10.1103/PhysRevLett.51.1336 [Cited within: 1]

[47]Bao W Jin S Markowich P A 2002 J. Comput. Phys. 175 487
DOI:10.1006/jcph.2001.6956 [Cited within: 1]

Bao W Jaksch D Markowich P A 2003 J. Comput. Phys. 187 318
DOI:10.1006/jcph.2001.6956 [Cited within: 1]

Bao W Jaksch D Markowich P A 2004 Multiscale Model. Simul. 2 210
DOI:10.1006/jcph.2001.6956 [Cited within: 1]

Bao W Chern I L Zhang Y Z 2013 J. Comput. Phys. 253 189
DOI:10.1006/jcph.2001.6956 [Cited within: 1]

[48]Wang H 2007 Int. J. Comput. Math. 84 925
DOI:10.1080/00207160701458369 [Cited within: 1]

[49]Wen L H Xiong H W Wu B 2010 Phys. Rev. A 82 053627
DOI:10.1103/PhysRevA.82.053627 [Cited within: 1]

[50]Mithun T Porsezian K Dey B 2014 Phys. Rev. A 89 053625
DOI:10.1103/PhysRevA.89.053625 [Cited within: 1]