# Tricriticalities and Quantum Phases in Spin-Orbit-Coupled Spin- Bose Gases

###### Abstract

We study the zero-temperature phase diagram of a spin-orbit-coupled Bose-Einstein condensate of spin , with equally weighted Rashba and Dresselhaus couplings. Depending on the antiferromagnetic or ferromagnetic nature of the interactions, we find three kinds of striped phases with qualitatively different behaviors in the modulations of the density profiles. Phase transitions to the zero-momentum and the plane-wave phases can be induced in experiments by independently varying the Raman coupling strength and the quadratic Zeeman field. The properties of these transitions are investigated in detail, and the emergence of tricritical points, which are the direct consequence of the spin-dependent interactions, is explicitly discussed.

###### pacs:

67.85.Bc, 03.75.Mn, 05.30.RtIntroduction — Ultracold atoms with spin-orbit (SO) coupling represent a very active field of research within the physics of quantum gases. They are characterized by highly nontrivial phase structures, arising when their peculiar single-particle properties are combined with the effects of the interparticle interactions. In the case of Bose gases, this results, e.g., in the appearance of spin-polarized phases as well as of spatially modulated (striped) configurations with supersolid-like features Boninsegni2012 () (see the recent reviews Dalibard2011review (); Galitski2013review (); Zhou2013review (); Goldman2014review (); Zhai2015review (); Li2015review (); Zhang2016review () and references therein).

Since the first achievement of a SO-coupled Bose-Einstein condensate (BEC) by the NIST group Lin2011 (), most studies have focused on spin- systems. Recent experiments in a gas of Rb atoms Campbell2016 (); Luo2016 () have succeeded in implementing spin-orbit coupling on spin- BECs, thus paving the way towards the exploration of the interplay between the spinor character and the SO coupling in determining the phase diagram of the system. Some properties of SO-coupled BECs with higher spin have been investigated in Wang2010 (); Xu2012 (); Ruokokoski2012 (); Lan2014 (); Song2014 (); Gautam2014 (); Gautam2015a (); Natu2015 (); Gautam2015b (); Pixley2016 (); Chen2016a (); Chen2016b (); Hurst2016 ().

The purpose of this Letter is to analyze the ground state of a SO-coupled BEC of spin in uniform matter and to point out the emergence of quantum tricriticalities, which are approachable in a highly controllable way by varying the Raman coupling and the quadratic Zeeman field. Tricriticality is an ubiquitous phenomenon in multicomponent systems Widom1996 (), like superconductors Kleinert1989 (); Green2005 (); Giovannetti2011 (), antiferromagnets Kato2015 () and QCD configurations Stephanov1998 (). Experimental evidence of quantum tricriticality has been recently reported in ferromagnetic superconductors Tokunaga2015 (). Controllability is crucial to reveal experimentally the predicted unconventional critical behavior of systems close to tricriticality, including the temperature dependence of the observables and their critical exponents Mercaldo2011 (); Kato2014 (); Kato2015 (), which can lead to the identification of novel universality classes Kato2014 ().

Single-particle physics — Let us consider a spin- BEC of particles in a volume , characterized by equally weighted Rashba Bychkov1984 () and Dresselhaus Dresselhaus1955 () SO couplings. The single-particle Hamiltonian reads (we set )

(1) |

where and is the spin- operator. Hamiltonian (1) describes an atom coupled to two Raman laser beams (see, e.g., the experiment Campbell2016 ()), within the rotating-wave approximation and in a frame connected to the laboratory frame by a unitary transformation , depending on space, time and Note1 (), and consisting of a proper generalization of the transformation employed in spin- gases Martone2012 (). This transformation, with an additional renormalization of the quadratic Zeeman term, yields the time-independent and translationally invariant Hamiltonian (1) Campbell2016 (); Campbell2015thesis (). The Raman coupling, quantified by , induces transitions between the three spin components of the spin- BEC, providing at the same time a momentum transfer along . The spinor nature of the system enables one to introduce effective linear and quadratic Zeeman terms, with strength and , respectively, which can be tuned independently from each other by varying the frequencies of the Raman lasers. Henceforth, we will focus on the case and study the interplay of Raman and quadratic Zeeman terms with interactions.

The single-particle spectrum of (1) has been analyzed by Lan and Öhberg in Lan2014 () and is characterized by three branches. At small and for large positive values of , the lowest branch has a single minimum at . At intermediate values of two local minima appear at , with , which become global degenerate minima when is further decreased. The minimum at eventually disappears for large negative . By increasing , the two curves delimiting the three-minima regime in the - plane and the line where the three minima are degenerate eventually merge in a critical point , beyond which a direct transition from the single-minimum to the double-minimum structure occurs.

Interacting gas — Two-body interparticle interactions introduce important novel features in the phase diagram. To investigate the properties of the interacting gas, we resort to the Gross-Pitaevskii (GP) mean-field approach. The ground-state properties of the system will be analyzed by minimizing the energy functional

(2) |

where is the three-component condensate wave function, with each component corresponding to the eigenvalues of . The density of the gas obeys the normalization constraint . The coupling constants and are related to the -wave scattering lengths in the spin- () and the spin- () channels Ho1998 (); Ohmi1998 (); StamperKurn2013 (). The form of the two-body interaction in (2) is characteristic of spin- BECs and reflects their rotational invariance. In the following, we shall focus on the ratios typical of Na () and Li () StamperKurn2013 ().

The stationary condensate wave functions can be found by solving the GP equation , with the chemical potential. Since we are interested in the ground state, we will consider solutions represented by a superposition of plane waves of the form

(3) |

where is the average density, the condensation momentum, are real three-component normalized spinors, and are complex coefficients obeying the normalization condition . The reality of is related to the fact that the expectation value in a stationary state of (2) vanishes at Note2 (). To motivate the form (3), first notice that it reproduces the ground-state wave function for a noninteracting gas, obtained by solving the linear Schrödinger equation for the single-particle Hamiltonian (1). In that case, only the terms with can be nonvanishing (see above); since the equation is linear, in the regimes where degenerate minima are present, the relative values of the corresponding coefficients are arbitrary.

As we shall see, interactions lift this ground-state degeneracy, just as in spin- SO-coupled systems Ho2011 (); Li2012PRL (). We have studied the ground state of the system investigating stationary solutions of the GP equation in the form (3) as functions of the relevant parameters , , , and . We find quantum phases characterized by condensation both in a single momentum state (plane-wave phases) and in a superposition of states with momenta (striped phases), whose coefficients satisfy . The contribution of the terms with in Eq. (3), appearing as a consequence of the nonlinearity of the GP theory Li2013 (), is negligible for low-density systems, and becomes more and more relevant as increases. Our results for the striped phases agree with Lan2014 ().

The global spin configuration of a spin- system can generally be described by means of two tensors, and , related to two competing kinds of spin order: ferromagnetic and nematic order, respectively Mueller2004 (); StamperKurn2013 (). The expectation value of the spin operator is the magnetization vector, while the expectation value of the spin-quadrupole tensor , with , is the nematicity tensor. The eigenvectors belonging to positive eigenvalues of identify the nematic directors. One can also define magnetization and nematicity densities and . The sign of the spin-dependent coupling plays a crucial role in establishing which kind of order characterizes the ground state of the system Ho1998 (); Ohmi1998 (); StamperKurn2013 ().

Antiferromagnetic interactions — Antiferromagnetic (or polar) spin-dependent interactions () tend to minimize the spin term of Eq. (2). Without SO coupling, if , the ground state of the system is represented by the longitudinal polar phase with wave function , while for one has the transverse polar phase , with an arbitrary phase factor Ho1998 (); Ohmi1998 (); StamperKurn2013 (). The two phases are separated by a first-order transition at . They both have and a single nematic director, which is parallel to the axis for the longitudinal polar phase, and lies in the - plane for the transverse polar state, with its direction fixed by .

The phase diagram in the presence of SO coupling is shown in Fig. 1, as a function of the Raman coupling and the quadratic Zeeman field . The upper part of the diagram is dominated by a zero-momentum (ZM) phase, with . Typical features of this phase are uniform density and vanishing longitudinal magnetization . At small and positive , the ground state approaches the above longitudinal polar state.

The lower part of the diagram of Fig. 1 corresponds to the region where the lowest branch of the single-particle spectrum has two degenerate minima. The competition between density and spin-density interactions in (2) breaks this degeneracy. At low , a polar striped (PS1) phase exists, characterized by spatial modulations and vanishing longitudinal magnetization . In this phase, only harmonic terms with odd are present in the wave function (3), and the two states with momenta are predominantly populated. Stripe modulations with wavelength appear in the total density , in the magnetization densities , , and in the nematicity densities , , , and . Notice that, in the limit, the polar striped phase approaches the transverse polar phase discussed above. The PS1 phase keeps a dominantly transverse polar character also at finite , since the nematicity densities , , and are the observables that oscillate with the largest contrast (see Supplemental Material Supplemental ()). However, the oscillations in and can augment significantly at large . All these quantities are not invariant under , and exhibit a different space and time dependence in the laboratory frame Martone2012 (). The total density and the nematicity density are instead unaffected by , and their behavior can be safely used to reveal the PS1 phase in experiments: the contrast of their modulations grows with and reduces at increasing , vanishing asymptotically at large negative values of the quadratic Zeeman field.

The energetic cost of density modulations increases with , until a transition occurs towards a plane-wave (PW) phase, where all particles condense in a single minimum with momentum . In this phase the density is uniform and the longitudinal magnetization is fixed by the momentum, as . The degeneracy with the state having opposite momentum (and magnetization along ) persists.

In a very narrow region between the ZM phase (above) and the PS1 and PW phases (below), a second polar striped (PS2) phase appears, featuring occupation of both odd and even- momentum states in (3). In the PS2 phase, having , the densities , , and oscillate with period , while , , , , , , and have oscillation wavelength as in the PS1 phase (see Supplemental Material Supplemental ()).

The PS1-PS2 and PS2-ZM phase transitions have a second-order nature, involving a smooth variation of the contrast of the stripes and of the momentum distribution. The PS1-PW and PS2-PW transitions are instead of first order, and both the contrast of the stripes and the magnetization can be used as experimentally testable order parameters.

The PS1-PS2, PS1-PW, and PS2-PW transition lines meet at the tricritical point , where the three phases coexist. The appearance of is an effect of the antiferromagnetic spin-dependent interactions and represents a key feature predicted by our work. Interactions also turn , which in the noninteracting model connects three different single-particle regimes Lan2014 (), into a second tricriticality, where the PS2-ZM and PS2-PW lines intersect. After , a direct PW-ZM transition occurs, whose second-order nature is witnessed by a smooth behavior of the magnetization . Up to corrections due to the spin-dependent interaction, the PW-ZM transition follows the line separating the nondegenerate () and twofold-degenerate () regimes of the single-particle ground state Lan2014 ().

For a fixed , as the point shifts toward the origin , while , whose position generally depends on interactions, approaches its noninteracting counterpart. In the same limit, the region of the PS1 phase collapses into the negative half-axis, and the PW phase, which is uniform and favored by density-density interactions, extends in the whole lower region of the phase diagram. The PS2 phase instead shrinks into the same line which, in the noninteracting limit, separates the nondegenerate and twofold-degenerate single-particle regimes, and connects the origin with Lan2014 ().

For a given value of the critical Raman coupling characterizing the PS1-PW transition approaches, as , a density-independent value, fixed by the ratio , as in the spin- case Ho2011 (). Remarkably, also the tricritical point approaches a finite value as . The behavior of the PS2 phase in the low-density limit is the same as for . A different scenario takes place for large densities: as increases, approaches while at the same time the PW phase shrinks and eventually disappears from the phase diagram, again as in the spin- case Li2012PRL (). However, the merging of the two tricritical points requires exceedingly high densities.

Ferromagnetic interactions — In the case of ferromagnetic spin-dependent interactions (), three phases are possible without SO coupling Ho1998 (); Ohmi1998 (); StamperKurn2013 (). For , the ground state is again the longitudinal polar state discussed above, while at the two degenerate longitudinal ferromagnetic states and are energetically favored. In the intermediate regime , the ground-state wave function reads , with . This state is characterized by ferromagnetic order in the transverse - plane, with the magnetization direction fixed by the arbitrary phase factor , and continuously approaches the polar state at , in a second-order phase transition. The transverse and longitudinal ferromagnetic states are instead separated by a first-order transition.

The phase diagram in presence of SO coupling is shown in Fig. 2. At low , the ZM and PW phases are still present in the upper and the lower part of the diagram, respectively. Their properties are essentially the same as in the antiferromagnetic case; however, for the ground state is in the PW phase also at small Raman couplings, and approaches the above longitudinal ferromagnetic phase as .

At low and another class of ground states is found, stemming from the transverse ferromagnetic state at . In this ferromagnetic striped (FS) phase, the harmonic terms with in (3) are predominantly populated, with . The FS phase results from the competition between the spin-dependent interaction term, proportional to , and the quadratic Zeeman field with strength . It has and, due to its connection with the transverse ferromagnetic phase at , the contrast of the density modulations is larger for the magnetization densities and than for the total density and the nematicity densities , , , and . All these quantities oscillate with wavelength , twice larger than in the PS1 phase for Lan2014 (), with relevant anharmonicities for the nematicity densities (see Supplemental Material Supplemental ()). As for the PS1 phase, in the FS phase the maximum contrast in at fixed is achieved by decreasing , i.e., going away from the upper transition line. However, the behavior of the contrast versus is not monotonic. Concerning the nematicity density , its contrast is generally neither monotonic in nor in .

The PW-FS transition has a first-order nature, which can be revealed by looking at the jumps in and in the contrast of the density modulations of the FS phase. The behavior of the contrast also witnesses the order of the FS-ZM transition. At small values of , the transition is expected to be second order as for . However, for our values of the parameters we find that, at larger , the contrast suddenly jumps from zero to a finite value, indicating the presence of a first-order phase transition.

The PW-FS and FS-ZM transition lines meet at the tricritical point , which exists as a consequence of the ferromagnetic spin-dependent interactions, and whose distance from the origin increases with . Passing one has a direct PW-ZM transition, which is of first order at small Raman couplings. However, the jump in the order parameter decreases when moving towards higher , while the region of the - plane where both the PW and ZM phases are present as local minima of the energy shrinks (see Fig. 2). Eventually, at the critical point , the PW-ZM transition becomes second order. Notice the difference in the nature of with respect to both the antiferromagnetic case and the noninteracting model Lan2014 (). The change in the order of the PW-ZM transition has already been confirmed in experiments with ferromagnetic Rb BECs Campbell2016 ().

For a fixed ratio , as the average density increases, the region of the FS phase in the phase diagram grows indefinitely along , and the position of shifts to the right. In the large- regime can meet and cross , which moves leftward, and the whole PW-ZM transition becomes second order. In the opposite limit, the FS phase shrinks into a finite arc of the curve connecting the origin with , similar to the PS2 phase in the antiferromagnetic case. At fixed , and for , the region of the FS phase collapses, together with , into the origin of the - plane.

Conclusions — Let us now comment on the experimental relevance of our results. The phase diagram for antiferromagnetic interactions could be explored in experiments with Na spin- BECs. The transition from the ZM to the striped phases can be studied by observing the change in the momentum distribution of the condensate, while the PW-ZM, the PS1-PW, and the PS2-PW transitions can be revealed measuring the longitudinal magnetization . The existence of the tricriticalities and can be deduced from the simultaneous study of these phase transitions. The possibility of tuning and arbitrarily close to such points opens realistic perspectives for the first experimental detection of tricriticality in ultracold atomic gases.

Similar considerations as above hold for ferromagnetic interactions. However, in current experiments with Rb BECs the investigation of the FS phase is made difficult by the small ratio StamperKurn2013 (), which compresses the stripes in a narrow region of the - plane. Possible strategies to improve the situation could involve atomic species with larger or tunable values of , or the implementation of properly chosen spin-dependent trapping conditions Martone2014 (); this would also make the merging of or with of more realistic achievement.

Further developments concern the study of the dynamic properties of spin- configurations, including the behavior of the roton minima in the excitation spectrum Martone2012 (); Zheng2013 (); Khamehchi2014 (); Ji2015 () and the emergence of a double gapless band structure in the striped phases Li2013 ().

###### Acknowledgements.

Insightful discussions with Giacomo Lamporesi, Yun Li, Chunlei Qu, and Zeng-Qiang Yu are acknowledged. This work was partially supported by PRIN Grant No. 2010LLKJBX “Collective quantum phenomena: from strongly correlated systems to quantum simulators,” by INFN through the project “QUANTUM,” by the Italian National Group of Mathematical Physics (GNFM-INdAM), by the QUIC grant of the Horizon2020 FET program, and by Provincia Autonoma di Trento. Note added — During the final preparation of the present work we became aware of two papers Sun2016 (); Yu2016 () discussing spin- SO-coupled BECs, whose predictions, when comparable, are in agreement with our findings.## References

- (1) M. Boninsegni and N. V. Prokof’ev, Rev. Mod. Phys. 84, 759 (2012).
- (2) J. Dalibard, F. Gerbier, G. Juzeliūnas, and P. Öhberg, Rev. Mod. Phys. 83, 1523 (2011).
- (3) V. Galitski and I. B. Spielman, Nature (London) 494, 49 (2013).
- (4) X. Zhou, Y. Li, Z. Cai, and C. Wu, J. Phys. B 46, 134001 (2013).
- (5) N. Goldman, G. Juzeliūnas, P. Öhberg, and I. B. Spielman, Rep. Prog. Phys. 77, 126401 (2014).
- (6) H. Zhai, Rep. Progr. Phys. 78, 026001 (2015).
- (7) Y. Li, G. I. Martone, and S. Stringari, in Annual Review of Cold Atoms and Molecules (World Scientific, Singapore, 2015), Vol. 3, chap. 5, pp. 201–250.
- (8) Y. Zhang, M. E. Mossman, T. Busch, P. Engels, and C. Zhang, Front. Phys. 11, 118103 (2016).
- (9) Y.-J. Lin, K. Jimenez-Garcia, and I. B. Spielman, Nature (London) 471, 83 (2011).
- (10) D. L. Campbell, R. M. Price, A. Putra, A. Valdés-Curiel, D. Trypogeorgos, and I. B. Spielman, Nat. Commun. 7, 10897 (2016).
- (11) X. Luo, L. Wu, J. Chen, Q. Guan, K. Gao, Z.-F. Xu, L. You, and R. Wang, Sci. Rep. 6, 18983 (2016).
- (12) C. Wang, C. Gao, C.-M. Jian, and H. Zhai, Phys. Rev. Lett. 105, 160403 (2010).
- (13) Z. F. Xu, Y. Kawaguchi, L. You, and M. Ueda, Phys. Rev. A 86, 033628 (2012).
- (14) E. Ruokokoski, J. A. M. Huhtamäki, and M. Möttönen, Phys. Rev. A 86, 051607 (2012).
- (15) Z. Lan and P. Öhberg, Phys. Rev. A 89, 023630 (2014).
- (16) S.-W. Song, Y.-C. Zhang, H. Zhao, X. Wang, and W.-M. Liu, Phys. Rev. A 89, 063613 (2014).
- (17) S. Gautam and S. K. Adhikari, Phys. Rev. A 90, 043619 (2014).
- (18) S. Gautam and S. K. Adhikari, Phys. Rev. A 91, 013624 (2015).
- (19) S. S. Natu, X. Li, and W. S. Cole, Phys. Rev. A 91, 023608 (2015).
- (20) S. Gautam and S. K. Adhikari, Phys. Rev. A 91, 063617 (2015).
- (21) J. H. Pixley, S. S. Natu, I. B. Spielman, and S. Das Sarma, Phys. Rev. B 93, 081101(R) (2016).
- (22) Li Chen, H. Pu, and Y. Zhang, Phys. Rev. A 93, 013629 (2016).
- (23) G. Chen, Li Chen, and Y. Zhang, New J. Phys. 18, 063010 (2016).
- (24) H. M. Hurst, J. H. Wilson, J. H. Pixley, I. B. Spielman, and S. S. Natu, arXiv:1608.01346.
- (25) B. Widom, J. Phys. Chem. 100, 13190 (1996).
- (26) H. Kleinert, Gauge Fields in Condensed Matter (World Scientific, Singapore, 1989), Vol. 1.
- (27) A. G. Green, S. A. Grigera, R. A. Borzi, A. P. Mackenzie, R. S. Perry, and B. D. Simons, Phys. Rev. Lett. 95, 086402 (2005).
- (28) G. Giovannetti, C. Ortix, M. Marsman, M. Capone, J. van den Brink, and J. Lorenzana, Nat. Commun. 2, 398 (2011).
- (29) Y. Kato and T. Misawa, Phys. Rev. B 92, 174419 (2015).
- (30) M. Stephanov, K. Rajagopal, and E. Shuryak, Phys. Rev. Lett. 81, 4816 (1998).
- (31) Y. Tokunaga, D. Aoki, H. Mayaffre, S. Krämer, M.-H. Julien, C. Berthier, M. Horvatić, H. Sakai, S. Kambe, and S. Araki, Phys. Rev. Lett. 114, 216401 (2015).
- (32) M. T. Mercaldo, I. Rabuffo, A. Naddeo, A. Caramico D’Auria, and L. De Cesare, Eur. Phys. J. B 84, 371 (2011).
- (33) Y. Kato, D. Yamamoto, and I. Danshita, Phys. Rev. Lett. 112, 055301 (2014).
- (34) Y. A. Bychkov and E. I. Rashba, J. Phys. C 17, 6039 (1984).
- (35) G. Dresselhaus, Phys. Rev. 100, 580 (1955).
- (36) For the laser geometry employed in Campbell2016 (), characterized by two different Raman transitions and , the unitary transformation reads .
- (37) G. I. Martone, Y. Li, L. P. Pitaevskii, and S. Stringari, Phys. Rev. A 86, 063621 (2012).
- (38) D. L. Campbell, Ph. D. thesis, University of Maryland, 2015.
- (39) T.-L. Ho, Phys. Rev. Lett. 81, 742 (1998).
- (40) T. Ohmi and K. Machida, J. Phys. Soc. Jpn. 67, 1822 (1998).
- (41) D. M. Stamper-Kurn and M. Ueda, Rev. Mod. Phys. 85, 1191 (2013).
- (42) In the single-particle picture this follows from the identity , holding for any eigenstate of the single-particle Hamiltonian (1). This result can be immediately extended to many-body configurations with rotationally invariant interparticle interactions.
- (43) T.-L. Ho and S. Zhang, Phys. Rev. Lett. 107, 150403 (2011).
- (44) Y. Li, L. P. Pitaevskii, and S. Stringari, Phys. Rev. Lett. 108, 225301 (2012).
- (45) Y. Li, G. I. Martone, L. P. Pitaevskii, and S. Stringari, Phys. Rev. Lett. 110, 235302 (2013).
- (46) E. J. Mueller, Phys. Rev. A 69, 033606 (2004).
- (47) See Supplemental Material for additional information about the behavior of the various kinds of densities in the striped phases.
- (48) G. I. Martone, Y. Li, and S. Stringari, Phys. Rev. A 90, 041604 (2014).
- (49) W. Zheng, Z.-Q. Yu, X. Cui, and H. Zhai, J. Phys. B 46, 134007 (2013).
- (50) M. A. Khamehchi, Y. Zhang, C. Hamner, T. Busch, and P. Engels, Phys. Rev. A 90, 063624 (2014).
- (51) S.-C. Ji, L. Zhang, X.-T. Xu, Z. Wu, Y. Deng, S. Chen, and J.-W. Pan, Phys. Rev. Lett. 114, 105301 (2015).
- (52) K. Sun, C. Qu, Y. Xu, Y. Zhang, and C. Zhang, Phys. Rev. A 93, 023615 (2016).
- (53) Z.-Q. Yu, Phys. Rev. A 93, 033648 (2016).