Elementary excitations in a spin-orbit-coupled spin-1 Bose-Einstein condensate

While a spin-orbit-coupled spin-1 Bose-Einstein condensate has been experimentally observed, its elementary excitations remain unclear in the stripe phase. Here, we systematically study the elementary excitations in three distinct phases of a spin-orbit-coupled spin-1 Bose-Einstein condensate. We find that the excitation spectrum as well as the corresponding static response function and structure factor depend strongly on spin-orbit coupling parameters such as the quadratic Zeeman field and the Rabi frequency. In the stripe phase, besides two gapless Goldstone modes, we show the existence of roton excitations. Finally, we demonstrate that quantum phase transitions between these different phases including the zero-momentum, plane wave and stripe phases are characterized by the sound velocities and the quantum depletion.


Introduction
Experimental realizations of spin-orbit-coupling in ultracold neutral atomic gases [1][2][3][4][5][6][7] open an avenue for studying exotic quantum phases and nonlinear dynamics [8][9][10][11][12][13] or simulating topological phases [14][15][16][17]. In neutral atoms, the so-called spins are pseudospins which may come from hyperfine states. The minimum system to support spin-orbit coupling is a spin-1/2 one which involves two hyperfine states. These hyperfine states are coupled by two oppositely propagating Raman lasers via a two-photon transition so that the momentum exchange between the Raman lasers and two hyperfine states generates atomic Raman-induced spin-orbit coupling. In such a system, intriguing ground-state phases for Bose-Einstein condensates (BECs) including the stripe, planewave and zero-momentum phases have been revealed. Quantum phase transitions between these three phases depend on spin-orbit coupling parameters and interaction strengths [18][19][20]. The stripe phase, which spontaneously breaks both the gauge symmetry and continuous translational symmetry, has been identified as possessing supersolid properties [21][22][23][24][25]. The plane-wave phase spontaneously breaks the timereversal symmetry and features nonzero magnetizations, and the zero-momentum phase is analogous to a conventional BEC without spin-orbit coupling. Such exotic BECs support intriguing elementary excitations: roton excitations in the plane-wave phase and two gapless Goldstone excitations in the stripe phase [18,20,26]. The roton excitations have been experimentally observed in the plane-wave phase [27,28].
In stark contrast to a solid-state system, spin-orbit coupling can be synthesized in a spin-1 cold atom system consisting of three hyperfine states. In fact, spin-orbit coupled spin-1 BECs have been realized in a recent experiment using three Raman lasers [29]. The BEC exhibits a richer ground-state phase diagram due to spinor interactions and tunable quadratic Zeeman shifts [29][30][31][32][33][34][35][36]. It has also been shown that there are multiple rotons in elementary excitations of the spin-1 plane-wave phase [32,33]. However, elementary excitations in the stripe phase of a spin-1 BEC have not yet been investigated (possibly due to the fact that it is much harder to study these excitations), and their fundamental properties such as Goldstone modes are still unknown.
In this work, we systematically study elementary excitations in a spin-orbit-coupled spin-1 BEC. Exotic ground states and phase transitions in this system make elementary excitations more interesting. In the stripe ground state, the existence of density modulations gives rise to a Bloch band-gap structure for elementary excitations. Such an excitation Bloch spectrum has two gapless Goldstone modes associated with the spontaneous breaking of two different continuous symmetries. We also find roton excitations at finite quasimomenta in the stripe phase, which is a unique feature of the stripe phase of a spin-1 BEC; such roton excitations have not been found in the stripe phase of a spin-1/2 BEC. In addition, we use the static response function and structure factor to characterize the elementary excitations. Both the static response function and structure factor diverge at a particular quasimomentum in the stripe phase, which is a salient feature of a stripe ground state. In the plane-wave ground state, apart from gapless excitations, there appear multiple rotons. The static response function presents sharp peaks around the location of rotons when a roton excitation is sufficiently softened. From the linear dispersion of the lowest excitation in all these phases, we can determine the sound velocities. We find that the sound velocities and the quantum depletion can be used as signals to reflect quantum phase transitions between different ground states.
This paper is organized as follows. In Sec. 2, we describe our theoretical frame for analyzing elementary excitations, response functions and quantum depletion. In Sec. 3, we show numerical results of elementary excitation for both spin-1 antiferromagnetic and ferromagnetic BECs. The conclusion is finally presented in Sec. 4.

Theoretical model
We consider a three-dimensional spin-1 spinor BEC with a Raman-induced spin-orbit coupling. Such spin-orbit coupling can be realized by applying three Raman lasers along the x direction to atoms with hyperfine states m F = 0, ±1. These lasers induce a two-photon Raman transition between two hyperfine states |m F = 0 and |m F = +1 (|m F = −1 ). The Hamiltonian of the system iŝ Here,Ψ = (Ψ 1 ,Ψ 2 ,Ψ 3 ) T represent the bosonic annihilation operators of the atoms with different pseudo-spins, andF = (F x ,F y ,F z ) are the spin-1 Pauli matrices.Ĥ SOC is the single-particle spin-orbit-coupled Hamiltonian, where ∇ ⊥ = ∂ 2 y + ∂ 2 z , Ω is the Rabi frequency of the Raman lasers which depends on the laser intensity, and is the quadratic Zeeman shift which can be tuned by controlling the detuning of the lasers [29]. In Eq. (1) and the following calculations, we set all quantities dimensionless. The units of the momentum, length, density, and energy are k Ram , 1/k Ram , k 3 Ram , and 2 k 2 Ram /m, respectively, where m is the atom mass and k Ram = 2π/λ Ram is the wave vector of the Raman lasers with λ Ram being the wavelength. g 0 and g 2 are spin-independent and spin-dependent interaction strengths, respectively, which depend on the atom number and the s-wave scattering lengths in the total spin channels [37]. In this work, the spin-independent interaction is always assumed to be repulsive (g 0 > 0), and the spin-dependent interaction could be either antiferromagnetic (g 2 > 0) or ferromagnetic (g 2 < 0).
In the mean-field approximation, the perturbations can be omitted and we have Ψ σ = ψ σ . The mean-field energy functional can be written as By minimizing the energy functional, we obtain the ground-state wave functions. Once we know the ground state, we can study elementary excitations of the system. The bosonic field operators can be written as the sum of the ground-state wave functions and fluctuation operators,Ψ with Here, q is the quasimomentum of the excitation, ω l is the excitation energy with l being the label of the lth band of the excitation spectrum, and µ is the chemical potential. u (l) q,σ (r) and v (l) q,σ (r) (σ = 1, 2, 3) are the perturbation amplitudes, which satisfy the normalization condition, a q,l andâ † q,l are the annihilation and creation operators of quasiparticles of the lth excitation band, respectively. Then we write the Heisenberg motion equation for the operatorΨ σ , By substituting Eqs. (1) and (4) into the above equation, we obtain the equations for ψ and δΨ. The equations for the ground-state wave functions are withĤ q,σ satisfy the Bogoliubov-de Gennes (BdG) equations [38], with X = X 1 + (g 0 + g 2 )X 2 + 2g 2 X 3 , The matrices X j and Y j (j = 1, 2, 3) are and with Atom density fluctuations can be induced when perturbations are applied on the system, which can be characterized by the density response function. In the linear response theory, the static density response function can be written as [20,39,40] where N is the total atom number. Here, the summation should include both the positive and negative excitation bands obtained by solving Eq. (9). Experimentally, excitations in a BEC can be probed by Bragg spectroscopy [27,41,42]. The Bragg beams transfer momentum and energy to the system. The excitation probability of Bragg scattering with momentum q can be described by the static structure factor [20], which can be treated as the sum of the excitation probabilities to different excitation bands. Static structure factor of a spin-orbit-coupled spin-1/2 BEC has been studied theoretically in Refs. [18,20], and experimental measurement in the plane-wave phase has also been reported [27].
The above theoretical frame belongs to the mean-field theory, which is valid when quantum fluctuations are negligible. The quantum fluctuations induce a fraction of the condensation depleted even at zero temperature. The quantum depletion in a spinorbit-coupled BEC can be calculated if the excitation spectrum and amplitudes are known [26,43,44]. The density of depleted atoms can be calculated as [26] where V is the volume of the system, and the perturbation amplitude v (l) q,σ depends on the quasimomentum q. If the number of depleted atoms is much smaller than the total atom number, then the mean-field approximation is valid.
There are plane-wave, zero-momentum and stripe phases for the gound state of a spin-orbit-coupled spin-1 BEC [32,33], which have different excitation properties. In the following, we will describe the theoretical treatments for different phases.

Plane-wave and zero-momentum phases
In the plane-wave or zero-momentum phase, the system has a uniform density profile. The ground-state wave functions of these two phases can be assumed as where n = N/V is the atom density, and k is the quasimomentum. If k = 0, then the system is in the plane-wave phase, and if k = 0, then the system is in the zero-momentum phase. The wave functions satisfy the normalization condition, ϕ 2 1 + ϕ 2 2 + ϕ 2 3 = 1. With the above wave functions, the energy functional in Eq. (3) becomes The parameters k, ϕ 1,2,3 and θ 1,3 can be determined by minimizing E. The excitation spectrum from Eq. (9) can be calculated by using the ansatz of plane-wave excitations [32] where U (l) q,σ and V (l) q,σ are plane-wave amplitudes. The excitation spectrum of the planewave phase exhibits phonon modes in the long-wavelength limit and roton modes at a finite momentum, which will be discussed in Sec. 3. We will also show that the excitation spectrum can be tuned by changing spin-orbit coupling parameters.

Stripe phase
In the plane-wave phase, roton instability emerges when the roton energy becomes zeros or even negative by tuning relevant parameters. In this situation, the plane wave is no longer the ground state of the system. The ground state turns out to be the stripe phase with a spatially modulated density profile [20].
Wave functions of the stripe phase can be constructed as the superposition of plane waves with different wave numbers, with the normalization condition L is the cutoff of the plane-wave modes, and K determines the period of the stripes. ϕ (j) σ and K are determined by substituting Eq. (20) into Eq. (3) followed by minimizing the resultant energy functional. This procedure shows that only coefficients of plane waves with j being an odd number (i.e., j = ±1, ±3, · · · .) are nonzero. We find that the period of the density is π/K, which is similar to that of the spin-1/2 case [20,26].
With the ground states, we can calculate the elementary excitation spectrum. Since the wave functions of the stripe phase is periodic, the BdG equations is also periodic and has the same period. Therefore, the perturbation amplitudes u (l) q,σ (r) and v (l) q,σ (r) are Bloch waves, which can be assumed to be Here, q is the quasimomentum and the superscript l indicates the lth band. By substituting Eq. (21) into the BdG equations Eq. (9), we obtain the excitation spectra and the corresponding excited states, from which the response function and structure factor can be computed.

Results and discussions
Now we are at the stage to study elementary excitations in the spin-orbit-coupled spin-1 BEC. In the following, we will study both the antiferromagnetic BEC (g 2 > 0) and the ferromagnetic BEC (g 2 < 0). (a1)-(a3) Excitation spectrum, (b1)-(b3) static response function and (c1)-(c3) static structure factor of a spin-orbit-coupled spin-1 BEC with an antiferromagnetic interaction. The first, second and third columns correspond to the zero-momentum phase, plane-wave phase and stripe phase, respectively. In (a1)-(a3), different lines represent different bands. In (a1) and (a2), there are three bands, and in (a3), only four lowest bands are plotted. In (a2), the vertical magenta-dashed line marks out the location of the roton at q x = −3. In (b1) and (b2), the red-dot-dashed, blue-dotted, and green-dashed lines represent response functions of the three bands, and the black-solid line denotes the total response function obtained by adding together the contributions of all positive and negative bands. In (b3), the red-dot-dashed and blue-dotted lines refer to the static response functions contributed by the lowest two gappless bands, and the black-solid line refers to the total static response function obtained by adding together the contributions of all bands in (a3). In (c1) and (c2), the red-dot-dashed, blue-dotted, and green-dashed lines describe the structure factor of different bands; their sum gives rise to the total structure factor denoted by the blacksolid line. In (c3), the red-dot-dashed and blue-dotted lines describe the structure factor of the lowest two bands, and the black-solid line describes the total one. π/K [K is defined in Eq. (20)] denotes the period of the stripe phase. The quadratic Zeeman shift is set as = −0.2, −0.7, and −1.5 in the first, second, and third columns, respectively. Here, Ω = 1.8, g 0 n = 1 and g 2 n = 0.1.

Antiferromagnetic BEC
The spin-orbit-coupled antiferromagnetic gas can undergo phase transitions between the the stripe phase, the plane-wave phase and the zero-momentum phase, which depend on the quadratic Zeeman shift and the Rabi frequency. It is also found that the stripe phase occupying the |m F = ±1 states can exist in a large parameter region [29,32], which makes it experimentally accessible to study its excitation properties. Based on previous works on the ground states of such a system [32,33,45], we present our results on elementary excitations in the following.
In Fig. 1, we show the excitation spectra of the three distinct phases (corresponding to different quadratic Zeeman shifts ) and their corresponding static response functions and structure factors. Here, the Rabi frequency and interaction strength are fixed at Ω = 1.8, g 0 n = 1 and g 2 n = 0.1, respectively. The three different phases can emerge by tuning [32]. Since the excitations along the y and z directions are irrelevant to spin-orbit coupling and thus have similar results as the conventional BEC, we here only consider the excitations along the x direction with spin-orbit coupling and set q y = q z = 0. When = −0.2, the ground state is the zero-momentum phase with k = 0. Our calculation shows the existence of three branches of excitation spectra, which is symmetric respect to the quasimomentum q x , as shown in Fig. 1(a1). The static response function χ(q x ) has a single local maximum at q x = 0, which corresponds to the ground state [see the black line Fig. 1(b1)]. Here, the red-dot-dashed, bluedotted, and green-dashed lines represent the contributions from the first, second, and third excitation bands in Fig. 1(a1), respectively. We note that the negative-band contributions to χ(q x ) are not shown here since we find that the response function of lth band possess the symmetry χ l (q x ) = χ −l (−q x ), where χ ±l (q x ) correspond to the energy ±ω l (q x ). Similarly, the static structure factor S(q x ) is plotted in Fig. 1(c1) with contributions from different excitation bands. The black-solid line refers to the total structure factor obtained by adding together the contributions of the three positive bands. We see that the three lines are symmetric with respect to q x arising from the symmetric excitation spectra.
When we change to −0.7, the system is in the plane wave phase. Figure 1(a2) shows the existence of one linear excitation and a roton excitation, which is consistent with previous results [18,32,33]. The roton excitation is a part of the lowest branch of the excitation spectrum, featuring a local minimum around a finite quasimomentum q rot ≈ −3, the location of which is marked out by the vertical magenta-dashed line. The minimum of the roton energy is so low that two peaks around q x = ±q rot emerge in the response function [see Fig. 1(b2)]. We note that the peak on the right side originates from unphysical negative-band contributions, which are not shown here. Different from the structure factor in the zero-momentum phase, the contributions of the three excitation bands are asymmetric; however, the total S(q x ) (the black-solid line) is always symmetric with respect to q x = 0. We also see that the contributions of the first and second bands have a local maximum, which corresponds to the local maximum of the corresponding excitation band. We also note that the structure of the static structure factor for the plane-wave phase in the spin-1 case is vary similar to that of the spin-1/2 case, which has been measured in a recent experiment [27]. In Fig. 1(a2), the low roton energy leads to two peaks in the curve of χ(q x ), which are higher than the peak at q x = 0. If we further decrease , we expect the occurrence of the roton instability induced by roton softening and the divergence of the response function. When this happens, the ground state becomes the stripe phase with spatially modulated density profile instead of the plane-wave phase. Because of the density pattern, the BdG equation in Eq. (9) is periodic with the same period π/K, leading to the excitation spectrum featuring Bloch band-gap structures as shown in Fig. 1(a3). In stark contrast to the zero-momentum and plane wave phases with one gapless mode, in the stripe phase, there are two gapless Goldstone modes arising from simultaneously spontaneous breaking of the gauge symmetry and the continuous translational symmetry [20]. The response function and the structure factor are shown in Figs. 1(b3) and 1(c3), respectively. At q x = 2K, we find σ,j |U (l,j) q,σ | 2 − |V (l,j) q,σ | 2 = 0 for l = 1, 2, so that the normalization condition for u (l) q,σ and v (l) q,σ in Eq. (6) is invalid. This leads to the divergence of χ(q x ) and S(q x ) at q x = 2K. A similar divergence has been reported in a spin-1/2 BEC [20]. The divergence of the static structure factor at the boundary of the Brillouin zone is a character of the spatially modulated density profile.
More interesting phenomena can be observed when we choose a small Rabi frequency so that the dispersion of the single-particle spin-orbit-coupled Hamiltonian H SOC has a three-well shape in the lowest branch. Figure 2 shows the excitation spectrum of the zero-momentum phase and the stripe phase for Ω = 0.2 and different . When = 0.2, the ground state of the system is the zero-momentum phase. In the excitation spectrum, there emerge two symmetric rotons, which is very different from Fig. 1(a1) where rotons do not exist. We note that by increasing Ω, the two rotons stiffen with enlarged roton gaps, leading to the disappearance of the rotons for a sufficiently large Ω. In this case, the ground state occupies the middle well of the three symmetric wells in the single-particle dispersion [32]. Rotons are excited at the positions corresponding to the other two wells. If we tune to a negative value, the ground state becomes the stripe phase [32]. There, compared with Fig. 1(a3), two gapless modes still exist in the excitation spectrum; interestingly, we also find roton excitations in the two lowest bands [see Fig. 2(b)]. Due to the periodicity of Bloch bands, these rotons are distributed periodically in the repeated Brillouin zone representation, reminiscent of the roton excitation of the spin-orbit coupled BEC in optical lattices [46,47]. In the latter case, rotons originates from the ground state occupying only one Bloch state with a nonzero quasimomentum [46,47]. Here, the emergence of roton structures in periodic excitation bands can be understood with the assistance of the single-particle dispersion from the Eq. (2). In this parameter regime, the single-particle dispersion have two symmetric wells located at k ± and a well at k = 0 with a higher energy. The ground state occupies the two symmetric wells, leading to the stripe phase [32]. If atoms can be excited to the middle well, roton structures emerge in the excitation spectrum.  Based on the fact that under the long wavelength limit, the lowest elementary excitations are always gapless with linear dispersion and the excitation spectrum depends strongly on spin-orbit coupling parameters, we expect that the phase transitions between these different phases can be characterized by the sound velocities. Figure 3(a) plots the sound velocities with respect to the quadratic Zeeman shift . When > −0.5, the ground state is the zero-momentum phase, and the sound velocity declines as we decrease . It exhibits discontinuity across the phase transition point around ≈ −0.5: In the plane wave phase, the sound velocity approaches zero whereas in the zeromomentum phase, it approaches a finite value. The discontinuity of the sound velocity indicates that the transition between the zero-momentum and plane-wave phases driven by is a first order phase transition, consistent with the previous results [32]. The sound velocity in the plane-wave phase rises as further decreases. At ≈ −0.8, the system enters into the stripe phase with two gapless Goldstone modes in the excitation spectrum, resulting in two different sound velocities [see Fig. 3(a)]. The discontinuity of the sound velocity at the transition point indicates that the transition between the plane wave phase and the strip phase is of the first order. In addition, we display the sound velocity with respect to the Rabi frequency in Fig. 3(b). When we increase Ω from a very small value, we see that the system changes from the stripe phase to the plane-wave phase and finally to the zero-momentum phase. Similar to the case shown in Fig. 3(a), the sound velocity experiences a discontinuity across the stripe and plane wave phases, implying that the phase transition is of the first order. Yet, different from the preceding case, the sound velocities approach zero as Ω approaches the transition point from either the plane wave phase or the zero-momentum phase, indicating that the transition between the plane wave and zero-momentum phases is of the second order, which is in agreement with previous results in Ref. [32].  We also plot the quantum depletion n ex as a function of in Fig. 4 [where the same parameters are chosen as in Fig. 3(a)]. The quantum depletion in the stripe and plane-wave phases rises with increasing, while it declines with increasing in the zero-momentum phase, which are similar to those in the spin-1/2 case [26]. The first order phase transitions between these three distinct phases have also been revealed by the discontinuous change of the quantum depletion. We note that in our parameter regimes, n ex /n is less than 5%, which verifies the validity of the mean-field treatments of our system. We also set g 0 n = 1 and g 2 n = −0.1.

Ferromagnetic interaction
Now, we turn to study the ferromagnetic case with g 2 < 0. In this scenario, the stripe phase can only exist within a small parameter region, which is very hard to be observed in experiments [29,32]. We also find that it is very difficult to obtain the accurate ground-state wave function for the stripe phase in the parameter regime, rendering the calculated excitation spectrum unreliable. We thus only consider the elementary excitations of the plane-wave and zero-momentum ground states. In Fig. 5, we plot the elementary excitation spectra of the ferromagnetic gas and its corresponding response functions and structure factors. For = −0.01, a doubleroton structure emerges in the excitation spectrum [see Fig. 5(a1)], as a result of the triple-well shape in the dispersion ofĤ SOC . The roton at q x ≈ −4 has a higher energy and the roton at q x ≈ −2 has a lower energy. The energy of the latter roton is so low that two peaks at q x ≈ ±2 emerges in the response function, as shown in Fig. 5(b1). In addition, the two maxons in the excitation spectrum induce a maximum in the curves of the first-band contribution to the structure factor [see Fig. 5(c1)]. By decreasing , we find that the roton at q x ≈ −4 softens and the roton at q x ≈ −2 stiffens [see Fig. 5(a2) and (a3)]. Consequently, the two peaks at q x ≈ ±2 in χ(q x ) disappear. If we choose a large negative and a large Ω, the dispersion of the single-particle spinorbit-coupled HamiltonianĤ SOC has a double-well shape, leading to the elementary excitation spectrum with a single roton, as shown in Fig. 5(a4). In this case, χ(q x ) and S(q x ) in Fig. 5(b4) and (c4) are similar to those of the antiferromagnetic gas, as shown in Fig. 1(b2) and (c2).  Once we know the excitation spectrum of the ferromagnetic gas, we can calculate the sound velocity. The sound velocity as a function of the quadratic Zeeman shift is demonstrated in Fig. 6. Similar to the antiferromagnetic case, the ferromagnetic gas can also undergo first order or second order phase transitions [32]. When Ω = 1.5, the ground state experiences a phase transition from the zero-momentum phase to the plane-wave phase around ≈ −0.4. The sound velocity is discontinuous at this critical value, indicating the occurrence of a first order phase transition [see Fig. 6(a)]. However, the phase transition is of the second order when we choose a large Ω, as illustrated in Fig. 6(b). These results of phase transitions are in accordance with the previous work where the phase transitions are directly determined from the ground-state energy [32]. We also plot the quantum depletion as a function of the quadratic Zeeman shift in Fig. 7 [where we choose the same parameters as in Fig. 6(a)]. Compared with the g 2 > 0 case, the quantum depletion in the case of g 2 < 0 is smaller, possibly due to the fact that the negative spin-dependent interaction lowers the system's energy. The first order phase transition between the plane-wave phase and the zero-momentum phase is also revealed by the sudden change in n ex across the phase transition, which is in consistence with Fig. 6(a).

Conclusion
In summary, we have systematically studied the elementary excitations in a spin-orbitcoupled spin-1 BEC for both antiferromagnetic and ferromagnetic interactions. We use the static response function and structure factor to characterize the excitation properties of the system. We further show that the sound velocity and quantum depletion can be used to identify the phase transitions between three different phases. Specifically, for ferromagnetic interactions, we reveal rich roton structures in the excitation spectrum of the plane-wave and zero-momentum phases. For antiferromagnetic interactions, we find that the stripe phase can not only support two gapless modes in the long wave regime, but also can support roton excitations at finite momenta. The excitation spectrum and structure factors can be measured in experiments by using Bragg spectroscopy in a 87 Rb BEC for the ferromagnetic case and in a 23 Na BEC for the antiferromagnetic case. Our results provide new insights in the exploration of the stripe phase in a spinor BEC and thus may be helpful for studying dynamics in such a system.