Emergent gauge field and the Lifshitz transition of spin-orbit coupled bosons in one dimension.

In the presence of strong spin-independent interactions and spin-orbit coupling, we show that the spinor Bose liquid confined to one spatial dimension undergoes an interaction- or density-tuned quantum phase transition similar to one theoretically proposed for itinerant magnetic solid-state systems. The order parameter describes broken Z2 inversion symmetry, with the ordered phase accompanied by non-vanishing momentum which is generated by fluctuations of an emergent dynamical gauge field at the phase transition. This quantum phase transition has dynamical critical exponent z ≃ 2, typical of a Lifshitz transition, but is described by a nontrivial interacting fixed point. From direct numerical simulation of the microscopic model, we extract previously unknown critical exponents for this fixed point. Our model describes a realistic situation of 1D ultracold atoms with Raman-induced spin-orbit coupling, establishing this system as a platform for studying exotic critical behavior of the Hertz-Millis type.

Introduction

Perhaps the first example of a quantum phase transition (QPT) was Stoner’s identification of a zero-temperature critical point distinguishing between unpolarized and spin-imbalanced Fermi liquids, and magnetic transitions in Fermi liquids have remained a rich subject since. These transitions for gapless, itinerant magnets belong to a class that is qualitatively distinct from transitions between gapped phases of matter and still remain mysterious despite the seminal works by Hertz and Millis[1-3]. While quantum criticality in Fermi liquids is a pervasive phenomenon in strongly correlated phases of matter, simple realizations of the ferromagnetic transition appear to be rare, especially in low-dimensional fermion systems where there is hope of a detailed theoretical understanding[4,5]. Nonetheless, a different paradigm for itinerant ferromagnetism appears if our initial degrees of freedom are bosons. In fact, since the ground state of an interacting (weak or strong) spinor Bose gas is a spin-polarized superfluid (SF)[6], ultracold bosons already provide a more natural realization of an itinerant ferromagnetic liquid compared to electrons in solid state systems requiring a Stoner instability. This work explores spin dynamics arising from the interplay of spin-charge separation, a concomitant emergent gauge field, and spin-momentum locking near the ferromagnetic QPT of a spin-orbit coupled interacting 1D Bose liquid[7-12].

The strongly interacting Bose liquid without spin-orbit coupling can be well understood by separating the excitations into a quadratically dispersing (but gapless) spin degree of freedom and a massless acoustic mode. Spin-charge separation in itinerant fermion ferromagnets has been profitably formulated in terms of emergent gauge fields, for example in the context of solid-state spintronics[13]. Likewise, in the Bose liquid, fluctuations of the spin degree of freedom behave as an emergent dynamical gauge field for the SF sound mode, the former coupled to the latter by an emergent electric field, as we show in this work. However, without any spin-dependent perturbations, the ferromagnetic ground state is also a fully-polarized spin eigenstate; therefore, in the absence of spin fluctuations the emergent field vanishes.

Spin fluctuations can be induced in the otherwise static spin-polarized gas by the addition of a helical Zeeman field. Qualitatively, a sufficiently strong field polarizes the local magnetization to be entirely parallel to it. However, the spin stiffness of the ferromagnetic liquid leads to an energy cost associated with the spatial variation of the magnetization, and this competes with the energy gain from precisely following the spatially rotating field. The result is that at intermediate values of the Zeeman field the system reduces its energy by developing an axial component of magnetization either parallel or antiparallel to the axis of the helical Zeeman field. Because the axial component of the magnetization is spatially uniform it does not contribute to the energy cost associated with the spin stiffness. By transforming to the rotating frame of the helical Zeeman field (yielding a frame with uniform spin-orbit coupling and a uniform transverse Zeeman field), one can see that this transition breaks the same symmetries as observed in the mean-field treatment of the spin-orbit coupled Bose gas[8,9,14,15]. However, in the mean-field case this transition is tuned by the intensity of the Zeeman field and has been understood almost exclusively in terms of single-particle physics: for a weak Zeeman field, there are two degenerate single-particle minima related by a Z2 inversion symmetry, and weak repulsive interactions favor condensation in one of these states, breaking the symmetry. At a sufficiently strong Zeeman field, the single-particle band structure changes such that there is a unique lowest-energy state in which to condense. In contrast, in the strongly interacting limit of interest to us it is more natural to understand the transition in terms of the competition between Zeeman energy and spin stiffness, the generalized rigidity associated with the interacting ferromagnetic Bose liquid.

We analyze the two-component strongly interacting Bose liquid subjected to a helical Zeeman field described by the Hamiltonian density (in units with ℏ = 1)where b = (b↑, b↓) represents the two bosonic fields describing the physical microscopic degrees of freedom, is the vector of Pauli matrices, and is the Raman-induced helical magentic field[7,8]. From the wavevector α we can define a natural unit of energy 4E = α2/2m, the recoil energy for the Raman laser[8]. We are interested here exclusively in the effect of spin-isotropic interactions g = g.

In the following, we first start from the general low-energy effective action for a spinor Bose liquid[16] and derive a low-energy effective action corresponding to Eq. (1) that describes the Bose liquid in terms of phase (i.e., number density) and magnetization degrees of freedom. We show that magnetization fluctuations are responsible for an emergent dynamical gauge field that influences the momentum of the superfluid through spin-orbit coupling. In the presence of a helical Zeeman field we find that there is a curve of quantum critical points of the system that can be reached by tuning interaction or density; in other words, sufficiently strong interactions can disorder the Z2 symmetry broken state at arbitrarily small Ω0, even though they cannot disorder the isotropic ferromagnet in the absence of the helical field. Finally we show that spin fluctuations become locked to density fluctuations, and as a result the latter are described by an effective field theory similar to the one-dimensional Lifshitz magnet proposed originally for ferromagnetic Fermi liquids[4]. The resulting critical fluctuations at the transition are qualitatively reminiscent of the non-Luttinger liquid predicted previously at the single-particle “flat band” transition in non- or weakly-interacting spin-orbit coupled bosons[15]. However, we uncover a previously unappreciated relevent interaction that we find drives the system to the much less well-understood interacting fixed point suggested in ref.[4]. We also find, fortuitously, that large interaction strength decreases the length scale and increases the temperature where these critical fluctuations can be observed in experiment. Following these analytic results, we validate our field theory analysis with detailed density matrix renormalization group (DMRG) simulations of the microscopic model Eq. (1). We confirm the presence of the interaction-tuned critical point, and we find new critical exponents for the order parameter and correlation length which differ from all previously studied transitions in this model, and from the zeroth- and the first-order ε-expansion predictions derived previously by Yang[4] and by Senthil and Sachdev[17]. We also obtain the dynamical exponent and find it consistent with the interacting Lifshitz transition value of z almost, but not exactly, 2.

Results

Effective Lagrangian approach for spinor bosons

Our goal is to analytically understand the properties of the strongly interacting Bose gas in a long-wavelength (compared to mean inter-particle spacing) helical magnetic field. Here, we restrict to this limit (analogous to Landau-Ginzburg theory) because the system is not integrable, so only universal properties such as the long wave-length limit are amenable to analytic treatment. As an added bonus, these universal properties should be relevant to other systems such as ferromagnetic Luttinger liquids.

To derive the low-energy long-wavelength properties we follow the effective field theory approach, which allows one to determine correlation functions from an effective Lagrangian that is determined solely by the symmetry of the system. The traditional classification of effective field theories based on symmetry was restricted to relativistic systems, where the symmetry group includes the Lorentz group[18,19]. Recent work by Watanabe and Murayama[16] has extended the classification of effective Lagrangians based on symmetry breaking to non-relativistic systems, including systems which are in the symmetry class of the spinor Bose liquid. In Methods, we translate their results for spinor bosons, first without spin-orbit coupling, to an effective Lagrangian (up to second order in derivatives) described by a phase and a spin degree of freedom which contains only four parameters,where ρ0, v, and κ are the average charge density, charge velocity, and spin-stiffness, respectively. Although ν has dimensions of velocity, it should not be interpreted as the spin velocity, as the spin excitations in the ferromagnetic Bose liquid disperse quadratically[20]. The charge density and spin-stiffness together determine the interaction-renormalized magnon mass m* = ρ0/κ. The degrees of freedom are a normalized vector field that represents the space-time texture of the spin and a phase ϕ whose space-time gradient (∂ϕ, ∂ϕ) is proportional to the fluctuations around the average density (ρ0) and to the momentum density respectively. The fields (a, a) are emergent gauge potentials satisfying the equation , where[21]

Since a gauge transformation may be off-set by a redefinition of the phase ϕ → ϕ + Λ, any choice of the vector potential (a, a) that has a gauge-invariant flux given by is sufficient.

The rather formal rigorous results (i.e. Eqs (2) and (3)) derived in Watanabe et al.[16] can also be understood intuitively on symmetry grounds. For example, the appearance of and ϕ as the appropriate degrees of freedom is expected from the symmetry of the ferromagnetic SF ground state[6] of spinor bosons, which forms a quasi Bose condensate (hence ϕ) and reduces the rotation symmetry from SU(2) (hence ) to O(2). If we ignore the vector potential terms a, then Eq. (2) is the usual Lagrangian for the phase degree of freedom for a Bose condensate while Eq. (3) is the gradient part of the field theory description (i.e. non-linear sigma model[22]) of a ferromagnet. The gauge potential (a, a), which is essential to obtain the correct dispersion for the ferromagnetic spin waves, arises from a complication that a local rotation of the condensate about the magnetization direction by δϕ(x, t) advances the phase ϕ(x, t) → ϕ(x, t) + δϕ(x, t). This is similar to how the application of a potential leads to a winding of the phase according to the Josephson relation. Thus, the winding of the phase in time is a combination of the applied external potential as well as the rotation of the condensate about the magnetization direction. To predict the effect of an external potential on the phase, we must therefore keep track of how much the condensate rotates along the condensate direction.

The simplest solution to this problem would be to avoid rotating the condensate around the magnetization direction . This could be accomplished by defining a vector orthogonal to (i.e. ) and ensuring that remains parallel as position is varied. Since we can think of as lying on a surface (parametrized by (x, t)) which is normal (locally) to the vector field . The problem of choosing to be locally parallel at nearby points is exactly that of parallel transport of the vector . This turns out not to be possible because of the holonomy associated with the curvature of the surface defined by the magnetization . Attempting to parallel transport in a small rectangle of size (δx, δt) leads to a rotation of the vector by an angle proportional to the Gaussian curvature . This suggests that attempting to avoid rotating the condensate when position is changed first along t and then along x leads to a net rotation about the magnetization direction compared to when the position is changed in the other order. This ambiguity in the net rotation leads to an ambiguity in the Berry phase that must be accounted for by a gauge potential whose curvature is given by Eq. (4).

Phase diagram in a helical Zeeman field

Given our understanding of the effective Lagrangians Eqs (2) and (3) of the fully symmetric spinor Bose gas, we now investigate the effects of the helical Zeeman field. This enters into the Lagrangian as the Zeeman potential induced by coupling the external field to the spin density, which is here the product of the normalized spin field and the number density , resulting finally in . This helical Zeeman field can be unwound using a position-dependent rotation of the spin vector around . This also therefore modifies the spin gradient as and thus also the gauge potential as a → a − αs3. Finally the effective Lagrangian with the helical field, in the rotating frame, becomes

The spin part of the Lagrangian () is now spatially uniform, and includes an isotropic ferromagnetic exchange, a Dzyaloshinskii-Moriya term, and an easy-axis anisotropy along , along with a uniform Zeeman field along . Similar 1D spin models have been studied previously in this context, though typically in the Mott-insulating limit on a lattice, where there is no back-action on the spin from charge fluctuations[23-28]. The charge part of the Lagrangian () now depends explicitly on spin from a dynamical vector potential αs3.

We now study the ground state of the system within the saddle point approximation in the limit of small α or large Ω, where the spin aligns along up to long-wavelength fluctuations. Neglecting these fluctuations, a zeroth-order saddle point approximation to the effective spin Lagrangian is

This is essentially the (real time) Landau-Ginzburg action for an Ising transition, favoring a nonzero s3 when the applied helical field satisfies Ω0 < α2/2m*, whether tuned by field strength, pitch angle, density, or effective spin stiffness. The magnon mass m* (or correspondingly, the spin stiffness) can be related to the dimensionless interaction strength γ as m* = m/f(γ) for a known function f[29], such that the ferromagnetic phase in Fig. 1(a) is given by the conditionwhich recovers the mean-field result in the limit m* = m.

Figure 1

(a) Phase diagram as derived from our effective Lagrangian approach. The critical Zeeman coupling Ω0 separating the Z2 ordered phase from the disordered phase depends strongly on the dimensionless interaction strength. (b) Evolution of a localized spin fluctuation. The time evolution was computed using the semiclassical equations of motion Eq. (13) on the critical curve. The initially localized fluctuation spreads as (red dashed lines), consistent with an expected dynamical critical exponent z ≃ 2.

The collective excitations of these phases are determined by the derivatives in Eqs (5) and (6). To understand these excitations we fix a gauge for the vector potential using the Wess-Zumino method[22] of extending the spin-texture field into an extra fictitious dimension parametrized by λ ∈ [0, 1], so that the field is the microscopic spin field and . Given this extension, the vector potential is

Using the fact that for a normalized field , , we can recover . We use the fact that the vector potential for small fluctuations simplifies from the full Wess-Zumino form (Eq. 9) to . The Zeeman field gaps out the magnon modes, however the mass of the s3 term vanishes near the transition and one can integrate out the massive s2 degree of freedom (with mass ). The dynamical term ρ0s2∂s3 leads to the massive phase with s2 ≃ ∂s3. After these manipulations, we obtain a simplified effective Lagrangian (exact as usual at sufficiently low energies) near the phase transition,where the spin velocity . The term from Eq. (5) that was linear in ∂ϕ has been eliminated for convenience by a constant shift ϕ(t) → ϕ(t) − At. We can now interpret the effective model as an Ising field theory gauge-coupled to a scalar boson. The remaining discrete symmetry of our model, however, is crucially not the usual Ising s3 → −s3, instead it is (s3, ∂) → (−s3, −∂) because of the coupling term.

Quantum critical point

While the phase diagram clearly suggests a Z2-breaking transition where the magnetization s3 orders, the itinerant nature of the magnet is expected to modify the critical properties of the transition. To understand the critical properties, we shift s3 as s3 → s3 − ∂ϕ/α (reflecting the symmetry above) and notice that this shifted field is gapped in the ordered phase, and so the charge and spin are locked as s3 = α−1∂ϕ. Substituting in this locking (and scaling length by α−1 and time by (vα)−1), the effective Lagrangian for ϕ alone is

The third term scales to zero in the long-wavelength limit compared to the first term, so in the long-wavelength limit near the critical point we neglect it in writing down the semiclassical equations of motion

We solve this numerically at the critical point and plot the resulting magnetization dynamics in Fig. 1(b). The simulation domain is a line of length 80/α with periodic boundary conditions, and our initial condition is a gaussian magnetization excess of height 0.02 and width α−1. We observe that the initial defect spreads with an envelope x ∝ t1/ with z ≃ 2.

Finally, we estimate the length and temperature scales where we expect this treatment to be valid. The length-scale cutoff for this phase-only action near the critical point comes from the Ising model, so both x, t must be larger than the inverse-gap in the spin sector, i.e. . Since f decreases with increasing interaction, the length and inverse-temperature scales required to observe critical behavior correspondingly decrease, a boon for experimental realization.

As suggested by the estimate z ≃ 2 for the dynamical critical exponent, this low energy effective model we obtained is actually a “Lifshitz” field theory and, including the relevent (∂ϕ)4 interaction, is identical to the effective action derived by Yang for itinerant Fermi liquid ferromagnets[4]. The upper critical dimension of this interaction is d = 2, and the low-order epsilon expansion (ε = 2 − d) predictions for critical exponents[4,17] are reproduced in Table 1, along with a summary of our numerical estimates presented in the following. It was also established in previous analysis that this novel critical point has a temperature dependence for the correlation length in the quantum critical regime of ξ ∝ T−2 and spin susceptibility χ ∝ T−1. Since the spin and density degrees of freedom are locked, the density susceptibility and correlation length also obey “non-Luttinger liquid” scaling at nonzero temperature[15].

Table 1

Scaling exponents for Eq. (12) near the critical point from epsilon expansion[4,17] compared to our numerical results.

ν	β	η	z
12(1+ε6) ~13	12(1−ε3) ~16	4ε^22250–0.2	2−η21.9–2

Numerical results

We next verify and extend the conclusions of the last section by direct numerical simulation of the ground state properties of the Hamiltonian in Eq. (1) using DMRG[30]. Specifically, we show that the continuous quantum phase transition shown in the phase diagram Fig. 1(a) indeed can be accessed by varying the interaction strength at a density where the mean interparticle distance is longer than the pitch of the helical magnetic field. We extract several critical exponents and substantiate our expectation of a Lifshitz-like critical point. As a practical limit in simulating Eq. (1), we discretize the Hamiltonian on a lattice and restrict our Hilbert space to states with at most two bosons per site. This should not affect our results significantly since we consider low density (1/5 boson per site) and relatively strong interactions, described by an isotropic Hubbard interaction U. We emphasize that the lattice discretization is for numerical convenience only; we are not introducing a physical optical lattice potential. For most of the calculations we kept up to 800 states to keep the truncation error per step around 10−12. However, when the interaction is very close to its critical value, we need to include more states (up to 2000) to achieve similar truncation error.

Phase diagram

The predicted Z2 phase transition is identified from the magnetization expectation value 〈s3〉 (In any finite system, the exact ground state is a “cat state” superposition of symmetry broken states with identically. However, the DMRG truncation procedure favors low entanglement, and since the energy difference compared to the exact ground state is exponentially small in system size, the DMRG converges on one of the symmetry broken quasi-ground states with nonzero ) and correlation function 〈s3(x)s3(x′)〉 on open chains while tuning the interaction strength U. Figure 2(a) shows the ground state spin-density expectation value for 300, 450, and 600-site chain systems with boson density of 1/5. The system undergoes a phase transition from a ferromagnetic (〈s3〉≠0) to a paramagnetic (〈s3〉 = 0) phase with increasing interaction strength; this is the same phase transition as a horizontal cut of Fig. 1(a). The weak dependence of the magnetization on the system size suggests that the interaction-tuned quantum phase transition is indeed continuous.

Figure 2

(a) Expectation value of magnetization as a function of interaction strength on open chains of different length. The boson density is fixed to 1/5. An interaction driven second-order phase transition is apparent. Data in the shaded region is less-converged due to a low energy charge modulated state. (Inset) Spin-momentum locking (Eq. (16)) is shown for N = 450. The appropriately scaled momentum tracks the magnetization exactly. (b) For N = 450, the critical interaction strength U is obtained by finding the exponent that gives the best linear fit for 〈s3〉1/(U), with U the extrapolated intercept. We find U ≃ 2.25, indicated as a red dot. (c) Using the obtained U, we confirm the magnetization critical exponent to be β ≃ 1/6 by a linear fit on a log-log plot.

In the numerical calculations, we observe a slightly density-modulated excited state in the paramagnetic region. Such states are an unavoidable artifact of the lattice discretization we employed, and cause convergence difficulty in the calculation. We checked that running an extensive number of sweeps decreases the amplitude of the modulation and the system eventually does converge to the expected paramagnetic state with uniform density. However, in the vicinity of the critical point, where the convergence is slowest, the wave functions we obtained were not converged enough to completely eliminate the density modulation; this less-converged region gives an overestimated value of 〈s3〉 and is shaded gray in Fig. 2(a,b).

Spin-momentum locking as signature of gauge coupling

The role of the magnetization generated by dynamical gauge fluctuations becomes clear from the spin-momentum locking property of the ground state. To see the spin-momentum locking notice that the canonical momentum operator derived from the microscopic Hamiltonian is given bywhile the number current operator derived from that same Hamiltonian is given by

A general theorem[31] rules out nonzero current density in the ground state (so 〈j〉 = 0). Given this constraint, the momentum and spin-density must be related by

Thus, the gauge coupling results in all Ising-symmetry-broken ground states having finite canonical momentum. The change in momentum as one crosses the phase transition from paramagnetic to the ferromagnetic must be attributed to the effective electric field Eq. (4) generated from gauge field fluctuations.

To observe the spin-momentum locking numerically, we also calculated the momentum expectation value 〈Π〉. The result, normalized by the factor α, is plotted in the inset of Fig. 2(a). The perfect agreement of the magnetization and normalized momentum explicitly shows the spin-momentum locking, Eq. (16), of the system.

Magnetization exponent

Next we demonstrate that the transition is not characterized by mean-field (i.e, ε = 0) critical exponents. For this purpose, an in-depth finite-size scaling is not necessary; instead, we estimate the critical point and exponents for a sufficiently large system.

From Fig. 2(a), one can see that finite size or boundary effects are already quite small for the N = 450 chain. Therefore we pinpoint the critical point U by finding the best linear fit (i.e., minimum residual) of 〈s3〉1/(U) over all values of 1/β (representative lines are shown in Fig. 2(b)), and then U is the extrapolated intercept. We then confirm the magnetization critical exponent β from a log-log plot using that value of U in Fig. 2(c). The critical point and exponent we identified are U ≃ 2.25 and β ≃ 1/6. Note that β ≃ 1/6 is both different from both the Ising critical point (β = 1/8) and the first-order ε-expansion result for the Lifshitz critical point in ref.[4] (β = 1/3). However, we do not take β ≠ 1/3 as an indication that this is not the interacting Lifshitz fixed point, but rather we suspect that the first-order ε-expansion is unreliable for such a large ε.

Dynamical critical exponent

Now we compute the dynamical critical exponent numerically to verify whether the transition remains essentially Lifshitz (i.e. z ≃ 2) as estimated from the semiclassical limit.

For this purpose, we compute the equal-time connected correlation function of s3. Near the critical point, i.e. when ξ~N, we expect the connected correlation function to decay as

The log-log plot in Fig. 3(a) shows that our estimated value for d + z − 2 + η ≃ 1.1, and we obtain the dynamical exponent to be z ≃ 2.1 − η.

Figure 3

(a) Log-log plot of the correlation function 〈s3(x0)s3(x)〉 as a function of x. N = 300 sites and U = 2.2, which is close to the critical interaction. As per Eq. (17), the linear slope of the log-log plot yields d + z − 2 + η ≃ 1.1. (b) Data collapse of the s3 correlation function. For different values of interaction (from U = 2.4 to U = 2.8) we rescale x − x0 by the correlation length ξ(U), to achieve data collapse. The dashed line is with a small shift, as a guide to the eye. (c) Correlation length varying U. The solid line is a fit yielding the exponent ν ≃ 1/3. The shaded region represents data collected from the less-converged wavefunctions, as in Fig. 2.

A secondary check for the value of d + z − 2 + η makes use of the scaling relation

The connected correlation function in Eq. (17) has a scaling form of ~F(x/ξ). We collapse the correlation function data at different values of U by scaling the distance (x), and from the scaling values we used we are able to extract the correlation lengths at different interaction strengths. The data collapse is shown in Fig. 3(b), and in Fig. 3(c) we plot the correlation lengths as a function of U. The correlation length decays as ξ ~ |U − U|−, and from the log-log plot in Fig. 3(c) we read off the correlation length exponent ν as ν ≃ 1/3. Plugging ν into Eq. (18), together with β ≃ 1/6 from the previous section, again gives d + z − 2 + η ≃ 1 which is close to the result we obtained from the s3 correlation function directly.

To first order in ε-expansion, η was predicted to be zero, but at second order the correction z = 2 − η/2 was derived in ref.[17] for this critical theory. From our estimate z ≃ 2 − η (using our estimates of β and ν) or z ≃ 2.1 − η (using the critical correlation function directly), we therefore infer a value η between 0 and 0.2, and corresponding z between 2 to 1.9. The small but non-zero anomalous dimension, along with the other estimated critical exponents that differ from their predicted values above the upper critical dimension, indicate that this is a distinct interacting fixed point of the one-dimensional model.

Discussion

In this work we analyzed universal properties of a strongly interacting spinor Bose gas in a helical magnetic field. We found that the rigorous description of the low-energy dynamics of the spinor Bose liquid[16] in terms of a scalar SF that interacts with the fluctuating magnetization as an emergent dynamical gauge field continues to apply with the addition of the helical magnetic field. The helical magnetic field is then responsible for an interaction-tuned quantum critical point, where sufficiently strong interactions can disorder an ordered Ising-like phase with a broken Z2 symmetry. The effective field theory expectation of a continuous quantum phase transition was then verified with detailed DMRG simulations. Although the results are quite different, we note that as a minimal model our study also complements recent work on the 1D Ising field with an interaction (as opposed to gauge) coupling to an acoustic mode[32].

Our effective field theory analysis yields a long wave-length effective Lagrangian valid near the quantum critical point identical to one proposed for one-dimensional ferromagnetic Fermi liquids[4] with a Lifshitz-like dynamics (i.e. z ≃ 2) and similar to that that proposed for “flat band” condensates of weakly interacting spin-orbit coupled bosons[15]. In the latter case, their expectation of a collective mode with Lifshitz-like dynamics (i.e. ω ∝ k2) arises straightforwardly from the underlying k4 spectrum for noninteracting bosons exactly at the flat band point. In our case, strong spin-independent interactions drive the transition even for weak Ω0, far from the flat band point. The ε-expansion results in ref.[4] suggest that interaction-induced fluctuations modify the quantum critical properties from classical, mean-field estimates, which we verify with exact numerics.

In addition to Fermi liquid ferromagnetism and the spin-orbit coupled Bose liquid studied here, Eq. (1) can also be the starting point to describe spinless bosons in flux ladders[33,34], where the “leg” of the ladder plays the role of pseudospin and the flux is a leg-dependent hopping phase (i.e., a pseudospin-orbit coupling). This platform is subtly different because of the necessary presence of an underlying lattice, but also supports an incredibly rich landscape of quantum phases and QPTs. In future work it would be interesting to determine if the interacting fixed point we uncover here is naturally realized there as well. More generally, Lifshitz critical points have also garnered substantial recent interest in higher spatial dimensions and in frustrated spin chains[35-39], and similar physics could even be relevant to superconductor-helimagnet heterostructures[40], with the superconductor providing the phase field and the direction of chirality of the helimagnet providing the Ising-like field.

The continuous phase transition in a strongly interacting gapless itinerant magnet demonstrated in this work shows that the spinor Bose liquid can be used as an experimentally realistic platform to study such quantum critical points. The combination of analytic and numerical results presented here show that the critical dynamics of the interaction-tuned Z2 ferromagnetic transition differ qualitatively from mean field. Our results for a simpler relative of itinerant quantum critical points in Fermi liquids may yield more general insight into those problems. The strongly interacting 1D spinor boson system in a helical magnetic field that is proposed in this work is already accessible in experiments on ultracold Rb and presents the ideal venue for the study of this class of criticality in the near future. Critical properties should be easily accessible from the temperature dependence of spatially resolved correlations in the system. Furthemore, given the slow timescale of the dynamics, this system (similar to the superfluid-Mott transition[41]) might provide an ideal platform to observe the surprising dynamics of this critical point.

Methods

Derivation of the effective Lagrangian

Following refs[18,19] the dynamics of the local order parameter of a system at sufficiently long wave-lengths and low frequencies is governed by a unitary operator U = exp(iπT), with the local order parameter encoded in π(x, t). This has been recently extended to non-relativistic systems by Watanabe and Murayama[16] for general classes of symmetry groups. Here we are interested in the low-energy, long-wavelength Lagrangian of interacting spinor bosons, possessing SU(2) rotation invariance of the spinor as well as Galilean invariance. As mentioned in the main text, the ground state spontaneously breaks Galilean invariance and reduces the rotation symmetry from SU(2) to O(2). From the effective field theory perspective, the low-energy dynamics is completely determined by an effective Lagrangian constrained by these symmetry (and symmetry-breaking) considerations.

We obtain the explicit form of the effective Lagrangian from the more general group-theoretic form given in ref.[16]. We start by recalling their results for the symmetry class of the spinor Bose gas, where the effective Lagrangian in terms of Galilean covariant derivatives contains only four parameters, up to second order in derivatives, and is written as

For convenience, we have combined π1,2 into π⊥. The Galilean covariant derivatives (containing the atomic mass m) are given asfor π⊥, while for π3 we have

Unpacking this notation, these covariant derivatives are constructed from components of a Maurer-Cartan form , so . Next, the covariant derivatives contain time and space derivatives of the π(x, t). Usingsuggests the shorthand notation for μ = x, t. (In terms of the notation of[16], , ). Substituting this into and retaining only terms up to second order in derivatives yields

Now we translate , via the spacetime order parameter derivatives encoded in , into experimentally relevant quantities: the magnetization and phase. We find it useful, following ref.[16], to decompose the transformation U as into where U0 is now a pure SU(2) rotation. (U need not include the Galilean generators at this point, which were used in ref.[16] to produce the Galilean covariant derivatives). In turn, we also find it useful to express U0 in terms of Euler angles (α, β, γ) as . However, π3 appears explicitly as the phase of U and must therefore be viewed as a function π3 ≡ π3(α, β, γ). The ρ can now be written as ρ = Ω + ∂π3σ0, where

The magnetization of the spinor gas in the absence of the field points along z so that s ∝ δ. Any quantity which transforms like a vector and takes the value (s1, s2, s3) = (0, 0, 1) for the ferromagnetic Bose gas must then be proportional to the magnetization. Following this argument, the magnetization can be taken as This is consistent with an intuitive picture where the state of the uniform ferromagnet is taken to be |Ψ0〉, and we are interested in the magnetization of a rotated state U0|Ψ0〉,where the trace including σ3 reflects that only the local z-component of has non-vanishing expectation in the fully-polarized state |Ψ0〉, as Tr[σσ3] = δ. In terms of the Euler angles, we can first define the orthogonal transformation , and then s = R, or (s1, s2, s3) = (sinα sinβ, −cosα sinβ, cosβ), with no dependence on γ.

We can also calculate derivatives of the magnetization s directly,

Now, using the identity and that ρ⊥ = Ω⊥ we get , and since R is an orthogonal matrix that preserves inner products,

The remaining terms of come from the scalar field ρ3. From the above, this is explicitlywith the vector potential a introduced to allow for the possibility of a curvature, i.e. . To determine a, we explicitly compute

Therefore the vector potential a, which satisfies can be chosen to depend only on and is therefore independent of γ. That is, we can choose ϕ(α, β, γ) as our third variable, with the only subtlety being that a must be given by the Wess-Zumino expression (i.e. Eq. 9).

Substituting the various ρ, we obtain the final form of the low-energy long-wavelength effective Lagrangian

The four previously-unassigned parameters are now given their physical significance: is the real time Lagrangian of an acoustic mode with sound velocity v in a liquid with average density ρ0. is the real time NL σM Lagrangian for a spin-1/2 ferromagnet with spin-stiffness κ. ν has velocity dimensions reflecting the spacetime anisotropy, although the spin excitation spectrum is quadratic.