ORIGINAL RESEARCH article
Sec. Optics and Photonics
Volume 9 - 2021 | https://doi.org/10.3389/fphy.2021.768799
SU(3) Spin–Orbit Coupled Rotating Bose–Einstein Condensate Subject to a Gradient Magnetic Field
- 1College of Intelligent Manufacturing, Sichuan University of Arts and Science, DaZhou, China
- 2Industry Technology Research Institute of Intelligent Manufacturing, Sichuan University of Arts and Science, DaZhou, China
- 3College of Chemistry and Chemical Engineering, Sichuan University of Arts and Science, DaZhou, China
- 4Department of Physics, Shaanxi University of Science and Technology, Xi’an, China
We consider a harmonically trapped rotating spin-1 Bose–Einstein condensate with SU(3) spin–orbit coupling subject to a gradient magnetic field. The effects of SU(3) spin–orbit coupling, rotation, and gradient magnetic field on the ground-state structure of the system are investigated in detail. Our results show that the interplay among SU(3) spin–orbit coupling, rotation, and gradient magnetic field can result in a variety of ground states, such as a vortex ring and clover-type structure. The numerical results agree well with our variational analysis results.
Spin–orbit coupling (SOC) plays an important role in a variety of physics branches. The realization of SOC in neutral atomic gases has attracted major attention both theoretically and experimentally [1–10]. Previous works have shown that the SOC can induce a variety of topological defects, such as skyrmions and half-quantum vortices and solitons, which are useful for the design and exploration of new functional materials [11–18]. Besides the ground-state structure, the dynamic properties of the system are significantly altered in the presence of SOC due to the close relationship between the spin and motional degrees of freedom in the topological excitations [19–24].
With the development of low-temperature technology, various types of SOC, such as Rashba, Dresselhaus, Weyl, and spin–orbit angular momentum, can be realized in experiments [1, 4, 25, 26]. However, most of previous studies of SOC have focused on the type of SU(2), where the internal states are coupled to their momenta via the SU(2) Pauli matrices. Only few works consider the SU(3) SOC, where the spin operators are spanned by the Gell-Mann matrices, which is more effective in describing the internal couplings among three-component condensates, such as a spin-1 Bose–Einstein condensate (BEC). Recently, Han and his co-authors have considered a homogenous SU(3) SO-coupled Bose gas and obtained the double-quantum spin vortices . On the base of their pioneering research work, Li and Chen have studied the SU(3) SO-coupled BEC confined in a harmonic plus quartic trap . Very recently, the ground states of a harmonically trapped spin-1 BEC with SU(3) BEC affected by the external rotation are investigated in , where a clover-type ground-state structure is discovered.
Nowadays, the gradient magnetic field has attracted more and more attention. More specifically, it is found that a variety of topological defects, such as a magnetic monopole and quantum knot, and even the artificial SOC can be realized by controlling the gradient magnetic field [30–34]. Li and his co-authors have investigated the ground state of three-component BEC in the gradient magnetic field and obtained the central Mermin–Ho vortex, magnetic monopole, and symmetry vortex lattices [35, 36]. They have also investigated the ground state of SU(2) SO-coupled BEC in the magnetic field and found that the skyrmion chain can induced by the isotropic SU(2) SOC . As far as we know, there is little study on an SU(3) SO-coupled BEC subject to both rotation and gradient magnetic field, which is we attempt to do. In this work, we consider an SU(3) SO-coupled rotating BEC subject to a gradient magnetic field and show that the system has a rich variety of ground states. Here, we want to note that, in real experiments, it is difficult to realize a rotating BEC with fixed SOC [38–40]. To do so, we can rotate the lasers creating SOC for an isotropic trap or rotate both the lasers creating SOC and an anisotropic trap, which leads to an effective time-independent Hamiltonian, Heff = H − ΩLz, describing the system in a rotating frame of Ref. .
The rest of this paper is organized as follows. In Section 2, we formulate the theoretical model describing an SU(3) SO-coupled BEC subject to both rotation and gradient magnetic field and briefly introduce the numerical method. The numerical results and its corresponding theoretical analysis through the variational approach are presented in Section 3. Finally, in Section 4, the main results of the work are summarized.
2 Model and Method
To begin with, we consider a quasi-two-dimensional (Q2D) spin-1 BEC with SU(3) SOC , which is confined in a controllable magnetic field . Under the mean-field approximation, the Hamiltonian of such a system can be written as [42–44]
The SU(3) SOC considered in the present work can be written as νso = κλ ⋅ p, where κ is the strength of SU(3) SOC and
Note that the SU(3) SOC in the Hamiltonian involves all the pairwise couplings between three states. This is distinct from the SU(2) case, where the state of ψ1(r) and ψ−1(r) cannot be coupled directly. In real experiments, the Hamiltonian with SU(3) SOC can be realized using a similar method of Raman dressing as in the SU(2) case, where three laser beams with different polarizations and frequencies, intersecting at an angle of 2π/3, are used for the Raman coupling . The external magnetic field is given by
where the length and time units are chosen as
3 Results and Discussions
In what follows, we will perform a detailed numerical study of the effects of system parameters on the ground-state structure. The richness of the present system lies in the large number of free parameters, which include the spin-dependent and spin-independent contact interactions, SU(3) SOC, rotation, and gradient field. To highlight the effects of the SOC, rotation, and gradient field, we focus on the antiferromagnetic condensate with c2 > 0 and fix the contact coupling parameters c0 = 100 and c2 = 3.
We first consider the system without SOC and fix the rotation frequency as Ω = 0.7. The typical density and phase distributions of such a system are shown in Figure 1 for different gradient magnetic fields, where only component 1 is shown as the other two components show similar behavior. In the absence of both gradient field (B = 0) and SU(3) SOC (κ = 0), the system is located at the center of the external harmonic potential, and no other lumps are formed. In this case, discrete vortices are formed due to the external rotation, but no vortex lattice can be formed due to the small contact interactions, as shown in Figure 1A. In the presence of a gradient magnetic field, the density distribution of the system shows a ring structure, and the size of center hole increases with the strength of the gradient magnetic field, as shown in Figures 1B,C for B = 1 and B = 2, respectively. In these cases, a vortex ring is formed along the ring direction (see Figure 1F for the phase distribution of Figure 1C), and the number of vortex also increases with the gradient magnetic field. Furthermore, if we increase the strength of gradient magnetic to B = 3 (see Figure 1D), most of vortices begin to move into the central hole, and the vortex ring is gradually destroyed. Upon increasing the gradient field to B = 4, all the vortices move to the central hole and form a giant vortex at the center of the trapping potential. In this case, no visible vortex remains, as shown in Figure 1E. Comparing Figures 1A–E, it is found that the radius of the ring increases with the gradient magnetic field.
FIGURE 1. Typical ground-state structures of a rotating BEC as a function of the gradient magnetic field in the absence of SU(3) spin–orbit coupling. (A)–(E) Density distributions of component 1 for different gradient magnetic fields: (A) B = 0, (B) B = 1, (C) B = 2, (D) B = 3, and (E) B = 4. (F) The corresponding phase distribution of (C). Other parameters are fixed as c0 = 100, c2 = 3, Ω = 0.7, and κ = 0, and the scale of each figure is [−12.8,12.8] in units of
Now, we turn our attention to the combined effects of SU(3) SOC and gradient magnetic on the ground-state structure of the system and thus set Ω = 0. Figure 2 shows the typical density and phase distributions of the system for the varying SU(3) SOC and gradient magnetic field. Previous studies on the SU(3) SOC Bose gases have shown that the triangular lattice structure is energy favorable for both homogeneous and confined systems [27, 29]. In the present system, we find a similar triangular lattice structure in the absence of the gradient magnetic field (B = 0), as shown in Figure 2A. In the presence of a gradient magnetic field, it is found that, with the increasing strength of the gradient magnetic field, the triangular lattice structure is gradually destroyed. In addition, the system evolves into three parts in space and eventually into a clover-type structure, as shown in Figures 2B–D for B = 1, 2, 3, respectively. Here, we want to note that a similar clover-type structure has been previously discovered in an SU(3) SO-coupled Bose gas with rotation . However, the clover-type structure found in the present work is induced by the combined effects of both SU(3) SOC and a gradient magnetic field, which give us another way to realize the clover-type structure in a non-rotating system. It is interesting to observe that the distance among such three space parts of the clover-type structure is unchanged if we increase the strength of SU(3) SOC but fix the strength of the gradient magnetic field, as shown in Figures 2D,E for B = 3, κ = 3 and B = 3, κ = 5, respectively. If we look at the phase distribution of the clover-type structure, as shown in Figure 2F for Figure 2D, it is easy to find that the local lump is in the plane wave phase, which is consistent with our previous results on the SU(3) SO-coupled system .
FIGURE 2. Typical ground-state structures of an SU(3) spin–orbit-coupled BEC as a function of SU(3) spin–orbit coupling and gradient magnetic field in the absence of rotation. (A)–(D) Density distributions of component 1 for fixed SU(3) spin–orbit coupling κ = 3 and for different gradient magnetic fields: (A) B = 0, (B) B = 1, (C) B = 2, and (D) B = 3. (E) Density distribution of the system for B = 3 and κ = 5. (F) The corresponding phase distribution of (D). Other parameters are fixed as c0 = 100, c2 = 3, and Ω = 0, and the scale of each figure is [-12.8,12.8] in units of
Finally, we move to the combined effects of SU(3) SOC, rotation, and the gradient magnetic field on the ground-state structure of the system. To highlight the effect of SU(3) SOC, we further fix Ω = 0.7. Figure 3 shows the typical ground-state structures of a rotating BEC as a function of the SU(3) SOC and gradient magnetic field. In the absence of SU(3) SOC, the system shows a ring structure with a vortex ring (see Figure 3A for B = 2 and κ = 0), which is consistent with that reported in Figure 1C. In the presence of SU(3) SOC, it is found that the annular structure is destroyed even for a small SU(3) SOC, as shown in Figure 3B for κ = 0.5 and B = 2, where two asymmetry lumps with visible vortices are formed. Further increasing the strength of SU(3) SOC, this tendency is becoming more and more obvious, and eventually, the system evolves into a single lump, as shown in Figures 3C,D, and F (the corresponding phase distribution of (Figure 3D)) for κ = 1 and κ = 2, respectively. This is different from the phase diagram in , where the density distributions show some symmetry and the center of mass of the system is still around the center of the trap (which also can be seen from their momentum distributions). For the present case, the symmetry of the system is further broken and the center of mass of the system deviates from the center of the trap. Actually, such difference can be understood by the fact that the SU(3) SO-coupled system has three discrete minima in the single-particle energy spectrum, and the number and weight of such three minima are selected for different strengths of SOC. Similar to the former case, it is interesting to observe that the distance between the lump and the center of trapping potential increases with the gradient magnetic field (see Figure 3E for B = 1 and κ = 2). Here, we want to note that we have also calculated other rotation frequencies, and the results show similar behavior.
FIGURE 3. Typical ground-state structures of a rotating BEC as a function of the SU(3) spin–orbit coupling and gradient magnetic field. (A)–(D) Density distributions of component 1 for a fixed gradient magnetic field B = 2 and for different SU(3) spin–orbit coupling: (A) κ = 0, (B) κ = 0.5, (C) κ = 1, and (D) κ = 2. (E) Density distribution of the system for κ = 2 and B = 1. (F) The corresponding phase distribution of (D). Other parameters are fixed as c0 = 100, c2 = 3, and Ω = 0.7, and the scale of each figure is [−12.8,12.8] in units of
To give a clearer understanding of the above results, we employ the variational method and analytically calculate the possible ground state. We begin with the single-particle energy spectrum of the system, which can be obtained by diagonalizing the kinetic energy and SOC terms. It is found that there exist three discrete minima residing on the vertices of an equilateral triangle in the momentum space. As discussed in , a threefold-degenerate plane wave ground state is selected as the ground state for a ferromagnetic condensate with SU(3) SOC, while three discrete minima with unequal (equal) weights are selected for the antiferromagnetic case. Consequently, the variational wave function can be written as Ψ = Ψ1 + Ψ2 + Ψ3, where
where Ψ1, Ψ2, and Ψ3 are the wave functions corresponding to the three minima of the single-particle spectrum, ( −2κ, 0),
It is found that the minimum energies are equal to
In this case, the second equation has the minimum energy
We have investigated the ground-state structure of a harmonically trapped rotating spin-1 BEC with SU(3) SOC subject to a gradient magnetic field. In the absence of SU(3) SOC, the system shows an annular structure, where a vortex ring is formed. In the absence of rotation, it is found that the clover-type structure discovered in the previous work can also be induced by the combined effects of SU(3) SOC and gradient magnetic field, and their distance shows sole dependence on the gradient magnetic field. When the rotation is included, we found that only one lump is formed in the three quadrants, and the distance between the lump and the center of trapping potential increases with the gradient magnetic field. Finally, we have employed the variational method and analytically calculated the possible ground state, which agrees well with our numerical simulations.
Data Availability Statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
All authors listed have made a substantial, direct, and intellectual contribution to the work and approved it for publication.
This work was supported by the National Natural Science Foundation of China under Grant Nos. 12075163, 12175129, 11775253, 12175027, and 11875010; the Dazhou Science and Technology Bureau under Grant No. 2019CYR0001; and the Natural Science Foundation of Chongqing under Grant No. cstc2019jcyj-msxmX0217.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
3. Cheuk LW, Sommer AT, Hadzibabic Z, Yefsah T, Bakr WS, Zwierlein MW. Spin-Injection Spectroscopy of a Spin-Orbit Coupled Fermi Gas. Phys Rev Lett (2012) 109:095302. doi:10.1103/PhysRevLett.109.095302
4. Huang L, Meng Z, Wang P, Peng P, Zhang S-L, Chen L, et al. Experimental Realization of Two-Dimensional Synthetic Spin-Orbit Coupling in Ultracold Fermi Gases. Nat Phys (2016) 12:540–4. doi:10.1038/nphys3672
7. Li J-R, Lee J, Huang W, Burchesky S, Shteynas B, Top FÇ, et al. A Stripe Phase with Supersolid Properties in Spin-Orbit-Coupled Bose-Einstein Condensates. Nature (2017) 543:91–4. doi:10.1038/nature21431
16. Hu H, Ramachandhran B, Pu H, Liu XJ, Drummond PD, Hu H. Spin-orbit Coupled Weakly Interacting Bose-Einstein Condensates in Harmonic Traps. Phys Rev Lett (2012) 108:010402. doi:10.1103/PhysRevLett.108.010402
17. Kato M, Zhang XF, Sasaki D, Saito H. Twisted Spin Vortices in a Spin-1 Bose-Einstein Condensate with Rashba Spin-Orbit Coupling and Dipole-Dipole Interaction. Phys Rev A (2016) 94:043633. doi:10.1103/physreva.94.043633
24. Sakaguchi H, Li B, Malomed BA. Creation of Two-Dimensional Composite Solitons in Spin-Orbit-Coupled Self-Attractive Bose-Einstein Condensates in Free Space. Phys Rev E Stat Nonlin Soft Matter Phys (2014) 89:032920. doi:10.1103/PhysRevE.89.032920
26. Chen H-R, Lin K-Y, Chen P-K, Chiu N-C, Wang J-B, Chen C-A, et al. Spin-Orbital-Angular-Momentum Coupled Bose-Einstein Condensates. Phys Rev Lett (2018) 121:113204. doi:10.1103/physrevlett.121.113204
28. Li CH, Qu C, Niffenegger RJ, Wang SJ, He M, Blasing DB, et al. Spin Current Generation and Relaxation in a Quenched Spin-Orbit-Coupled Bose-Einstein Condensate. Nat Commun (2019) 10:375–02. doi:10.1038/s41467-018-08119-4
29. Peng P, Li G-Q, Wang B-H, Cao Z-Z. Exotic Ground States in a SU(3) Spin-Orbit Coupled Spin-1 Bose-Einstein Condensate under Rotation. Chaos, Solitons & Fractals (2020) 141:110332. doi:10.1016/j.chaos.2020.110332
32. Aidelsburger M, Atala M, Lohse M, Barreiro JT, Barreiro B, Bloch I. Realization of the Hofstadter Hamiltonian with Ultracold Atoms in Optical Lattices. Phys Rev Lett (2013) 111:185301. doi:10.1103/physrevlett.111.185301
33. Kennedy CJ, Siviloglou GA, Miyake H, Burton WC, Ketterle W. Spin-orbit Coupling and Quantum Spin Hall Effect for Neutral Atoms without Spin Flips. Phys Rev Lett (2013) 111:225301. doi:10.1103/physrevlett.111.225301
35. Li J, Yu YM, Zhuang L, Liu WM. Dirac Monopoles with a Polar-Core Vortex Induced by Spin-Orbit Coupling in Spinor Bose-Einstein Condensates. Phys Rev A (2017) 95:043633. doi:10.1103/physreva.95.043633
41. Liu CF, Juzeliūnas G, Liu WM. Spin-orbit Coupling Manipulating Composite Topological Spin Textures in Atomic-Molecular Bose-Einstein Condensates. Phys Rev A (2017) 95:023624. doi:10.1103/physreva.95.023624
42. Liu CF, Liu WM. Spin-orbit-coupling-induced Half-Skyrmion Excitations in Rotating and Rapidly Quenched Spin-1 Bose-Einstein Condensates. Phys Rev A (2012) 86:033602. doi:10.1103/physreva.86.033602
47. Zhang XF, Dong RF, Liu T, Liu WM, Zhang SG. Spin-orbit-coupled Bose-Einstein Condensates Confined in Concentrically Coupled Annular Traps. Phys Rev A (2012) 86:063628. doi:10.1103/physreva.86.063628
Keywords: Bose–Einstein condensate, SU(3) spin–orbit coupling, gradient magnetic field, rotation, GP equation
Citation: Chen G-P, Tu P, Qiao C-B, Zhu J-X, Jia Q and Zhang X-F (2021) SU(3) Spin–Orbit Coupled Rotating Bose–Einstein Condensate Subject to a Gradient Magnetic Field. Front. Phys. 9:768799. doi: 10.3389/fphy.2021.768799
Received: 01 September 2021; Accepted: 29 September 2021;
Published: 02 November 2021.
Edited by:Zhaoxin Liang, Zhejiang Normal University, China
Reviewed by:Qing Sun, Capital Normal University, China
Wei Yi, University of Science and Technology of China, China
Xiang-Fa Zhou, University of Science and Technology of China, China
Copyright © 2021 Chen, Tu, Qiao, Zhu, Jia and Zhang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.