Arxiv:2104.07559V1 [Cond-Mat.Str-El] 15 Apr 2021

Hydrodynamics of quantum spin liquids

Vir B. Bulchandani,1 Benjamin Hsu,2 Christopher P. Herzog,3 and S. L. Sondhi4 1Princeton Center for Theoretical Science, Princeton University, Princeton, New Jersey 08544, USA 2Amazon Web Services, Inc., Seattle, WA 98109∗ 3Mathematics Department, King’s College London, The Strand, London, WC2R 2LS, UK 4Department of Physics, Princeton University, Princeton, New Jersey 08544, USA (Dated: April 16, 2021) Quantum spin liquids are topological states of matter that arise in frustrated quantum magnets at low temperatures. At low energies, such states exhibit emergent gauge fields and fractionalized quasiparticles and can also possess enhanced global symmetries compared to their parent microscopic Hamiltonians. We study the consequences of this emergent gauge and symmetry structure for the hydrodynamics of quantum spin liquids. Specifically, we analyze two cases, the U(1) spin liquid with a Fermi surface and the SU(4)-symmetric “algebraic” spin liquid. We show that the emergent degrees of freedom in the spin liquid phase lead to a variety of additional hydrodynamic modes compared to the high-temperature paramagnetic phase. We identify a hydrodynamic regime for the internal U(1) gauge field common to both states, characterized by slow diffusion of the internal transverse photon.

I. INTRODUCTION tools based on tensor-network methods13–15 yield quan- titative predictions in remarkable agreement with ex- 16–18 Hydrodynamics provides a universal description of the perimental results . These studies complement new long-wavelength, low-frequency dynamics of physical sys- theoretical methods for deriving hydrodynamics from 19,20 tems as they approach a state of thermal equilibrium. the AdS-CFT correspondence , random quantum cir- 21–23 24,25 This dynamics is typically determined by a handful of cuits and effective field theory . classical partial differential equations, whose structure is In this paper, we explore the linear-response hydrody- dictated by the symmetries of the system at hand. The namics of quantum spin liquids. These unusual topolog- physical insight underlying this enormous reduction in ical states of matter arise in frustrated quantum mag- complexity is the notion of a “conservation law”: only lo- nets at low temperatures and exhibit non-trivial low- cal physical quantities that are conserved in every single energy degrees of freedom26,27. These degrees of free- microscopic collision can give rise to cumulative effects dom include fractionalized quasiparticles with enlarged that are discernible on macroscopic scales. Hydrodynam- flavour symmetries, and internal gauge fields coupled to ics in this sense is equally applicable to both quantum these quasiparticle degrees of freedom. This enhanced and classical systems. Although the emergence of hy- low-energy structure can give rise to an elaborate struc- drodynamics in many-body classical systems has been ture of hydrodynamic slow modes compared to the high- appreciated for centuries, its emergence in many-body temperature paramagnetic phase—even though neither quantum systems is less well understood. Some notable phase breaks any symmetries of the microscopic Hamil- exceptions include the theory of transport in metals and tonian. The extra modes come in part from the emergent liquid helium-3, as based on Landau’s Fermi liquid pic- gauge fields and thus an emergent magnetohydrodynam- ture1,2, and the theory of superfluidity in weakly inter- ics and in part from the new global flavour symmetries. acting Bose gases3,4. We illustrate this behaviour by focusing on two repre- The idea that hydrodynamics might be a useful tool sentative examples of quantum spin liquid states, namely for understanding dynamics in a much broader range of the U(1) spin liquid with a Fermi surface, that is conjec- quantum systems has only come into focus comparatively tured to describe the uniform resonating valence bond recently. Simultaneous advances in various experimen- (RVB) state of the square-lattice Heisenberg antiferro- tal and theoretical directions have provided an intense magnet27, and the SU(4)-symmetric “algebraic” spin liq- stimulus for the study of hydrodynamics in many-body uid, that describes its “d-wave” RVB state28. While quantum systems. One active field of inquiry focuses on arXiv:2104.07559v1 [cond-mat.str-el] 15 Apr 2021 these are much studied spin liquids and hence appro- high-mobility two-dimensional electron gases, for which priate objects of further investigation, our focus on these a long-anticipated regime of viscous electron flow5–7 has two dimensional examples compels us to address an addi- recently been realized in materials such as graphene8 and tional subtlety arising from microscopic compactness of 9 PdCoO2 . Another rapidly growing field is the hydrody- the emergent gauge field. This is the fact that spin liq- namics of one-dimensional quantum systems, for which uid states in two dimensions and at non-zero temperature analytical tools based on integrability10–12 and numerical are unstable at the longest length and time scales, since the possibility of topological excitations of the internal gauge field ultimately leads to spinon confinement2829. Thus the question arises of whether quantum spin liq- ∗ Work done while the author was at Princeton University. uids in two dimensions exhibit a hydrodynamic regime 2 at all. We argue below that such a regime indeed exists, surface is neither Galilean nor Lorentz invariant, and provided that monopole operators in 2+1 dimensions are there are no internal symmetries beyond the emergent sufficiently irrelevant, as has been conjectured for both U(1) symmetry. Meanwhile, the field theory describing the U(1) spin liquid with a Fermi surface30 and Dirac spin the algebraic spin liquid is Lorentz invariant, quantum liquid states with many fermion flavours31. We find that critical, and exhibits multiple quasiparticle flavours with despite their fundamental differences, the internal U(1) an internal SU(4) symmetry. gauge symmetry common to both these states gives rise Despite these differences, both spin liquid states ex- to a diffusive hydrodynamic mode of the internal trans- hibit various hydrodynamic features in common, most verse photon, that is a spin-liquid analogue of magnetic notably a diffusive transverse mode of the internal gauge diffusion in ordinary conductors32,33. photon. In physical terms, this corresponds to an ex- citation of the internal magnetic field b, and therefore couples to the “spin chirality” P = S1 · S2 × S3 about II. GENERALITIES a plaquette, that has been the focus of earlier propos- als for experimental detection of spin liquid states35–38. A. Models We now review the standard hydrodynamics of the high- temperature phase. For concreteness, consider the family of antiferromag- netic, spin-1/2 Heisenberg models on a square, two- dimensional lattice, B. High-temperature hydrodynamics X H = J Si · Sj + ..., (1) In the high-temperature, paramagnetic phase, the hy- hiji drodynamic modes correspond to the na¨ıve conservation laws of the model, namely energy and the three compo- where the ellipsis denotes terms consistent with both lat- nents of spin, tice symmetry and SU(2) symmetry. At high temper- atures, such models will exhibit standard paramagnetic ∂tnE + ∇ · jE = 0, (3) α α behaviour. At low temperatures, they can potentially ∂tnS + ∇ · jS = 0, α = 1, 2, 3. (4) host a variety of spin liquid phases that depend on the microscopic details of the model in question. The field In a given equilibrium state, the expectation values of theory of these states can be studied with the aid of slave- the particle currents vanish, so that the linear response fermion variables34 behaviour of the currents near equilibrium is given by the constitutive relations 1 X Sj = f † σj f , (2) i 2 iα αβ iβ α,β δjE = −DEE∇δnE, (5) α α δjS = −DSS∇δnS, α = 1, 2, 3. (6) which yield an exact rewriting subject to the single- P † occupancy constraint α fiαfiα = 1. Various “spin liq- to leading order in a derivative expansion; these are de- uid” saddle points may then be studied at the level of coupled by symmetry. Thus the paramagnetic phase ex- mean-field theory. Physics in the vicinity of these saddle hibits four diffusive hydrodynamic modes in total, with points is characterized by the emergence of fractionalized dispersion excitations and dynamical gauge fields associated with 2 these excitations. One question that has not been stud- iωE = DEEk , (7) ied systematically is how far these additional excitations 2 iωSα = DSSk , α = 1, 2, 3. (8) and enhanced symmetries lead to novel signatures in the hydrodynamics of quantum spin liquids, compared to the One might ask how far interaction effects, that appear high-temperature paramagnetic regime. in hydrodynamics as nonlinear corrections to the consti- This is the question we address below, focusing on two tutive relations, can modify this na¨ıve picture of normal representative examples of spin liquid states, the U(1) diffusion. For example, the leading nonlinear corrections spin liquid with a Fermi surface, corresponding to the to Eqs. (5) and (6) (note that simple total derivatives, 39 uniform RVB , or “zero flux state”, of the Heisenberg e.g. δnE∇δnE, do not generate hydrodynamic tails ) model27, and the SU(4)-symmetric algebraic spin liquid, have the form corresponding to the d-wave RVB, or “staggered flux 31 α α (2) α state” . For each of these states, there is an effective δjS = −DSS∇δnS − DSESδnE∇δnS + . . . , α = 1, 2, 3. low-energy field theory description consisting of multiple spinon degrees of freedom minimally coupled to a dy- Such nonlinear corrections lead to long-time tails namical, non-compact, U(1) gauge field. However, the hJ (t)J i ∼ 1/td/2+1 in the current autocorrelation func- physics of these two effective theories is very different. tions of diffusing charges39 and tails ∼ 1/td/2 for ballistic The field theory describing the spin liquid with a Fermi ones40. In two dimensions, nonlinear diffusive corrections 3 are nonsingular (∝ ω) as ω → 0+, while nonlinear ballis- try holds approximately, so that the gauge field a may tic corrections give rise to weak (logarithmic) finite-size be modelled as non-compact, while the spinon degrees dependence of transport coefficients40. We will therefore of freedom have sufficient time to reach a state of local disregard all such effects in the remainder of this paper, equilibrium. The question now arises of whether such a −1 −1 as they do not modify the long-time dynamical exponent regime exists at all, i.e. whether the condition τa  τψ measured in a given, finite, two-dimensional system. holds in some regime of temperature and disorder. The We now examine how this simple picture of four, de- temperature dependence of these timescales is not well- coupled, diffusive hydrodynamic modes breaks down in studied. For concreteness, we focus on the clean case, for −4/3 spin liquid states at low temperature. which it has been argued that τψ ∼ τP ∼ T due to scattering off the internal gauge degrees of freedom41,42. In particular, momentum is relaxed in the hydrodynamic III. HYDRODYNAMICS OF THE U(1) SPIN regime. LIQUID WITH A FERMI SURFACE We now put forward a Fermi Golden Rule argument for the temperature dependence of the timescale τa for At low temperatures, the square-lattice Heisenberg fluctuations in the U(1)flux quantum number. Denote by 27 model is conjectured to support a zero-flux state . The H0 the quantum Hamiltonian associated with the lead- effective action for low-lying excitations about this state ing terms in the Lagrangian Eq. (9). As written, the consists of two fermion (“spinon”) fields coupled to an Lagrangian description pertains to a non-compact gauge emergent non-compact U(1) gauge field, which can be field. Compactness of the gauge field can be restored by written as a Lagrangian density26,27 adding “vortex creation” operators  Z   ∗  2 † 2 cl X † ∗ 1 e mx = exp i d x ax · e (10) L = ψ (i∂t + e a0)ψσ − ∇ − ia ψσ 0 0 σ 2m∗ c∗ σ=1,2 cl to the non-compact effective Hamiltonian H0. Here, a 1 x0 + (e2 − b2) + ... (9) denotes a singular classical gauge texture with a 2π-flux 8π vortex centered at x0, so that by virtue of the canoni- 0 0 to leading order in fields. For clarity, we have restored cal commutation relations [ai(x), ej(x )] = iδijδ(x − x ), ∗ the operator m† realizes a vortex in the quantum gauge dimensions for the speed of light c of the internal U(1) x0 ∗ field a. Notice that the operator m† is invariant under field, as well as the effective mass m and internal U(1) x0 charge e∗ of the spinon, but set ~ = 1. We shall also gauge transformations of a, because the only operator use vector notation e, b for the Maxwell fields, where it appearing in its definition is the gauge-invariant quan- is understood that ez = bx = by = 0. Our immediate tum field e. The definition of the classical gauge texture acl is not unambiguous, but the details beyond its 2π goal will be to derive the linear-response hydrodynamics x0 of the effective action given by Eq. (9). flux are not important43. Schematically, this gives rise to the expression Z Z  A. Regime of validity H = H + g d2x cos d2x acl · e + ... (11) 0 0 x0

Let us first address the regime of validity of a hydro- for the true Hamiltonian of the system. Thus we have dynamic description based on the continuum field theory included compactness of the gauge field as a perturbation Eq. (9), which turns out to be quite subtle. Superfi- to the effective Hamiltonian; this is entirely analogous to cially, the action Eq. (1) exhibits two symmetries that how compactness of the target space is encoded in 1 + 1 are not possessed by the underlying lattice model Eq. (1), dimensional Luttinger liquid theory, for example, when namely continuous translation symmetry and U(1) gauge modelling Kosterlitz-Thouless-type transitions44. symmetry. At the level of hydrodynamics, these give rise In the Euclidean time description, the operators m† x0 to conservation of momentum R d2x P and conservation create three-dimensional monopoles, as used by Polyakov R 2 of the total U(1) flux, Φ = d x bz. Both symmetries to argue for confinement in compact 2 + 1 dimensional are approximate; even in the absence of disorder, the Maxwell theory45. In real time, such monopoles are cor- underlying lattice relaxes momentum on a timescale τP rectly viewed as instantons, i.e. tunneling events that associated with Umklapp scattering, while microscopic change the U(1)flux quantum number of the state un- compactness of the emergent gauge field a leads to a vio- der consideration28. The stability of the U(1) spin liquid 28 lation of U(1)flux symmetry , giving rise to fluctuations with a Fermi surface at non-zero temperature T > 0 is in the topological U(1)flux quantum number Φ on some contingent on irrelevance of the vortex creation opera- timescale τ to be determined. For comparison with these tors m† . Following various authors30,46, we shall assume a x0 timescales, it will be useful to introduce a scattering time that the vortex creation operator exhibits scaling with a τψ for the spinon degrees of freedom. highly irrelevant scaling dimension da  3. Neither as- −1 In what follows, we shall focus on the regime τa  sumption is trivial, as the spin liquid with a Fermi sur- −1 ω  τψ . In this regime of frequencies, U(1)flux symme- face is neither scaling invariant nor known rigorously to 4 be stable to confinement in two dimensions28,47; we refer B. Linear-response hydrodynamics to Ref.30 for an argument in favour of both hypotheses. (We note that both assumptions are easier to justify for In the regime of frequencies and wavenumbers given by 28,31 Dirac spin liquids , as discussed in Sec. IV below.) Eqs. (13) and (14), spinon degrees of freedom are in lo- Subject to these assumptions, the timescale for fluc- cal thermal equilibrium with each other and with micro- tuations in Φ can be deduced as follows. By Fermi’s scopic fluctuations of the gauge field, while the gauge field Golden Rule, the timescale for instantons to generate itself can be modelled as non-compact. For non-zero tem- scattering between eigenstates of H0 scales with the peratures T > 0, this means that the effective Lagrangian coupling strength of the vortex creation operator, as Eq. (9) admits a classical magnetohydrodynamic regime −2 τa ∼ g . This represents the shortest timescale on largely analogous to the usual hydrodynamics of charged which the U(1)flux quantum number in a thermal state particles coupled to macroscopic Maxwell fields, except of the effective Hamiltonian H0 can vary due to instan- that momentum is no longer conserved. ton corrections. It remains to deduce the temperature We now derive the linear-response magnetohydrody- dependence of g. To this end, let g0 denote the bare namics of Eq. (9). First note that conservation of energy value of the coupling g at some cut-off scale Λ  T . At and spin yields four hydrodynamic modes, low temperature, the coupling g is renormalized from its d T
a ∂tnE + ∇ · jE = 0, (15) bare value g0 by a factor g = Λ g0. It follows that d α α g ∼ T a , which yields ∂tnS + ∇ · jS = 0, α = 1, 2, 3, (16)

τ ∼ T −2da . (12) just as for the paramagnet. These must, however, be a augmented by the (2 + 1)D Maxwell equations,

Strong irrelevance of the vortex creation operator da  ∇ · b = 0, (17) 3 then implies a clear separation of scales τ −1  τ −1, a ψ 1 giving rise to a hydrodynamic regime ∇ × e = − ∂ b, (18) c∗ t

−1 2da −1 4/3 ∇ · e = 4πn, (19) τa ∼ T  ω  τψ ∼ T . (13) 1 4π ∇ × b = ∂te + j. (20) Associated with these frequency cut-offs is a range c∗ c∗ of allowed wavenumbers. The microscopic motion of Together, these imply conservation of the internal U(1) spinons is expected to be ballistic at a speed fixed by charge (spinon density) the Fermi velocity, leading to a microscopic length-scale l = v τ for spinon scattering. Meanwhile, on fre- ψ F ψ ∂tn + ∇ · j = 0. (21) −1 −1 quency scales τa  ω  τψ , the internal gauge field is both locally equilibrated and approximately non- In the interests of treating gauge and spinon degrees of compact. This means that the magnetohydrodynamic freedom on equal footing, we have combined “field” and description discussed below holds sway, indicating dif- “matter” contributions in the energy density nE. The fusive hydrodynamic transport of the transverse photon advantage of this approach is that effects such as Joule mode with a dynamical exponent z = 2, according to Eq. heating, which consists of a transfer of energy from fields (27). This allows us to self-consistently deduce a length- to matter but does not violate local energy conserva- 1/2 scale la = (Daτa) below which the emergent magnetic tion, are taken into account automatically. Note also field exhibits hydrodynamic behaviour, where the diffu- that because momentum is relaxed on the hydrodynamic −1 −1 sion constant Da is defined below in Eq. (28). Since we timescale, τP ∼ τψ , there are no sound modes (in con- have assumed that microscopic scattering off the gauge trast to the Dirac spin liquid, as discussed in Sec. IV field is the dominant contribution to spinon relaxation, below). −1 −1 4/3 we additionally have Da ∼ σ ∼ τψ ∼ T , implying Let us restrict to the linear response regime, and study a hydrodynamic regime of wavenumbers small perturbations e = e0 + δe and n = n0 + δn about a stationary background charge distribution ∇·e0 = 4πn0. −1/2 da−2/3 −1 4/3 (Daτa) ∼ T  k  (vF τψ) ∼ T (14) We assume general, isotropic, linear-response constitu- tive relations of the form associated with the hydrodynamic regime of frequencies, α α Eq. (13). Finally, we note that above the frequency and δjS = −DSS∇δnS, (22) −1 −1 δj = σδe − D ∇δn − D ∇δn , (23) wavenumber scales τψ and lψ , the spinon degrees of nn nE E 0 freedom do not have time to equilibrate with each other δjE = σ δe − DEn∇δn − DEE∇δnE (24) or with the internal gauge field. This defines a “spinon plasma” regime, whose behaviour is beyond the scope of Note that the local magnetization is decoupled from the this work. The range of possible dynamical regimes is other hydrodynamic modes at linear order, since no other summarized in Fig. 1. gradients transform appropriately under spin rotation. 5

Spinon plasma Spinon hydrodynamics

Confined

FIG. 1. Schematic plot of the hydrodynamic regime in quantum spin liquids. The hydrodynamic frequency range is bounded −1 −1 on either side by the rate of fluctuations in the U(1)flux quantum number τa and the rate of spinon scattering τψ . At timescales smaller than τψ, the spinon degrees of freedom are no longer in equilibrium with each other, and the system is in a “spinon plasma” regime. At timescales larger than τa, compactness of the microscopic gauge field becomes apparent and leads to confinement. A magnetohydrodynamic regime of locally equilibrated spinon degrees of freedom coupled to an approximately non-compact gauge field arises between these two extremes. The origin of the associated wavenumber scales is discussed in the main text. The temperature dependence of these cutoffs is not well-studied in general and depends on the model at hand. For the U(1) spin liquid with a Fermi surface, estimates are provided in Eqs. (13) and (14), while for Dirac spin liquids with a large number of fermion flavour Nf  1, estimates are provided in Eqs. (31) and (32).

Similarly, vector currents cannot couple to gradients of constant the pseudoscalar bz at linear order. c∗2 We find five hydrodynamic modes in total (see Ap- Da = . (28) 4πσ pendix for details). There are three diffusive magnetiza- tion modes, The transverse photon is an excitation of the slow modes {δe⊥, δb, δj⊥, δjE,⊥}. Upon fixing δe⊥, the induced mag- 2 netic field fluctuations δb = (c∗/ω)k × δe by Amp`ere’s iωSα = DSSk , α = 1, 2, 3, (25) ⊥ law, while the transverse current fluctuations are propor- as in the paramagnetic phase. There is again a diffusive tional to δe⊥. The main point is that the transverse pho- energy mode ton can be excited by driving fluctuations of the internal magnetic field b. The latter is related to the microscopic  0  35 σ spin chirality P = S1 · S2 × S3 about a plaquette. This iω = D − D k2 + O(k4) (26) E EE σ nE might provide one physical route towards coupling to the hydrodynamic transverse photon mode36–38. corresponding to fluctuations of δnE. However, in the This mode can be viewed as the spin-liquid analogue of spin liquid phase, this mode is accompanied by fluctua- magnetic diffusion that arises in ordinary conductors32. tions of the longitudinal slow modes {δn, δek, δjk, δjE,k}. In more detail, the hydrodynamic interpretation of the Finally there is an additional hydrodynamic mode “skin effect” for electromagnetic fields in an electrical compared to the paramagnet, which arises due to the conductor is that it describes a crossover from ballistic to emergent U(1) gauge symmetry of the spin liquid phase. diffusive behaviour33. In the present case, the crossover This is the transverse photon mode, scale k∗ such that diffusive behaviour sets in for k  k∗ ∗ ∗ −4/3 is given by (see Appendix) k = 2πσ/c ∼ τψ ∼ T , 2 4 iωa = Dak + O(k ), (27) which lies well beyond the UV cutoff in Eq. (14) at low temperatures. Thus the long-wavelength hydrodynamics which in the presence of a non-zero conductivity σ > 0 of the transverse internal photon is diffusive, rather than is diffusive at long wavelengths, with effective diffusion ballistic as one might na¨ıvely expect. 6

We conclude this section with some speculative re- invariance and the presence of strong spinon-gauge inter- marks on the hydrodynamics of the spinon Fermi surface actions, the spinon scattering rate is expected to exhibit 41,42 51 −1 with conserved momentum. Standard treatments , Planckian scaling τψ ∼ T . We again estimate the based on physically realistic lattice geometries that are timescale for fluctuations of the U(1)flux quantum num- −2da conjectured to yield a spinon Fermi surface state, find ber as τa ∼ T , but for the algebraic spin liquid with a that τP and τψ are the same order of magnitude, so small number of fermion flavours Nf = 4, as in Eq, (29), that momentum is always relaxed in the hydrodynamic it is not known whether instantons are irrelevant31, i.e. regime. However, the effective action Eq. (9) is mean- whether da > 3. For a large number of fermion flavours 28 ingful for Fermi surfaces that do not relax momentum Nf  1, it is expected that da = O(Nf ), so that a efficiently. For example, a suppression of Umklapp scat- hydrodynamic regime tering would imply a momentum relaxation rate expo- 48,49 −1 −1 O(Nf ) −1 nentially large in inverse temperature , so that τP  τa ∼ T  ω  τψ ∼ T (31) −1 τa . In this case, one would expect ballistic sound modes of charge, momentum and energy on generic grounds50, exists. The associated length scales may be determined −1 but the usual analysis is complicated by the possibility of as above. On hydrodynamic timescales ω  τψ ∼ T , momentum exchange between field and matter degrees of the gauge degrees of freedom are again expected to prop- freedom via Poynting flux. In particular, it is not clear agate diffusively with dynamical exponent z = 2, leading 1/2 1/2−da −O(N ) whether a ballistically propagating charge mode is phys- to a length scale la ∼ (Daτa) ∼ T ∼ T f in ical; a definitive resolution of this question seems to lie this case. Meanwhile, the microscopic motion of spinons ∗ beyond linear-response hydrodynamics. is again expected to be ballistic with speed vF = c , leading to a hydrodynamic regime of wavenumbers

−1/2 O(Nf ) ∗ −1 IV. HYDRODYNAMICS OF THE ALGEBRAIC (Daτa) ∼ T  k  (c τψ) ∼ T. (32) SPIN LIQUID Thus we expect a robust magnetohydrodynamic As a more sophisticated illustration of the hydrody- regime in the limit of a large number of fermion flavours, namic approach, we now consider the hydrodynamics of as depicted in Fig. 1. The existence of this regime is less the SU(4)-symmetric “algebraic spin liquid” that arises clear for Nf = 4; nevertheless, the basic Nf = 4 case in the staggered flux phase of the square-lattice Heisen- serves as an instructive template for the hydrodynamics berg model. The low-energy effective action for this state of spin liquid states with higher flavour symmetry and can be expressed as the Lagrangian density31 emergent Lorentz invariance. Finally, we note that the continuum Lagrangian Eq.   ie∗  1 (29) is translationally invariant although the underlying L = Ψ −iγµ ∂ + a Ψ+ (e2 −b2)+... (29) µ c∗ µ 8π lattice model is not. Since there is no Fermi surface, however, the lowest-order processes in the lattice model where Ψ denotes an eight-component spinor, that relax momentum are Umklapp scattering events be- tween Fermi points, which cost an energy E ∼ v |G|, ψ11 U F where |G| = 2π/a is the magnitude of a reciprocal lat- ψ12 Ψ =   , (30) tice vector on the square lattice with lattice spacing a. ψ21 Thus for a clean system, the momentum scattering time ψ22 EU /kB T τP ∼ e is exponentially large at low tempera- 48,49 in which each ψia is a two-component spinor, i denotes a tures , and momentum must be included as a con- valley index and a is the spin index. Following Ref.31, we served mode on the timescales of Eq. (31). denote Pauli matrices acting on the Dirac, valley and spin degrees of freedom by {τ i, µi, σi} respectively, and define Ψ = iΨ†τ 3. The U(1) internal and gauge symmetries B. Linear-response hydrodynamics act similarly to the uniform flux phase considered above. However, the SU(2) spin symmetry now sits within a The hydrodynamics of the algebraic spin liquid can larger, internal SU(4) symmetry flavour group, that is again be modelled by the equations Eq. (15)-(21) in- i i i j generated by the fifteen operators {σ , µ , σ µ }. troduced for the U(1) spin liquid with a Fermi surface, i.e. the Maxwell equations and conservation of energy and spin. However, there are several important differ- A. Regime of validity ences compared to the Fermi surface case. First, since −1 −1 τP  τa due to microscopic suppression of Umklapp As above, the regime of validity of the hydrodynamic scattering, momentum must be treated as a conserved description we propose is bounded on either side by the mode in the putative hydrodynamic regime. This im- scattering rates for instanton and spinon degrees of free- plies an additional hydrodynamic equation −1 −1 dom, with the separation of scales τa  τψ a nec- essary condition for its existence. By emergent Lorentz ∂tPi + ∂jΠij = 0, i, j = x, y, (33) 7 where P = (Px,Py) is the local momentum density and Notice that at low temperatures, the wavenumber scale Π denotes the momentum flux tensor. (As for the energy k∗ = 2πσ/c∗ ∼ T −1 below which the internal transverse density nE in the previous section, we have implicitly photon crosses over from ballistic to diffusive behaviour combined “field” and “matter” contributions in the defi- is once again well beyond the UV cutoff for the hydrody- nitions of Pi and Πij.) Second, the effective Lagrangian namic regime of wavenumbers Eq. (32). Eq. (29) is Lorentz invariant. This means that in the By virtue of the reactive coupling in Eq. (37), the limit that lattice effects are neglected, the equations of energy mode combines with the longitudinal momentum motion associated with Eq. (29) are Lorentz invariant. mode to yield a ballistic sound mode for the algebraic Finally, there are twelve additional flavour modes cor- spin liquid, at least on the timescales t  τP for which responding to the internal SU(4) symmetry. However, momentum appears conserved. This phenomenon has these transform non-trivially under elements of the lat- been noted in other condensed matter systems with emer- tice space group and are therefore decoupled from the gent Lorentz symmetry, such as graphene52,53. The result slow modes of Eqs. (15) - (21) and (33) at the level of (see Appendix) is that the energy and longitudinal mo- linear response, and may be considered separately. mentum modes combine to form a “sound mode” with The modifications due to emergent Lorentz symmetry diffusive broadening, can be derived as follows. First notice that linear dis- 2 sipative corrections include lattice effects and are there- ωs = ±csk + O(k ), (41) fore exempt from the requirement of Lorentz covariance. In fact, the only substantive consequence of emergent q ∂p ∗ where the effective speed of sound cs = c . If Lorentz symmetry is that the energy and momentum op- ∂nE erators must transform as a tensor in the limit that lattice we further assume an ultrarelativistic, two-dimensional effects are neglected. In this limit, the operator equations equation of state p = nE/2, on the grounds that the of motion for energy and momentum associated with Eq. effective action Eq. (29) describes purely massless exci- (29) take the form tations, we recover the universal relation for sound modes in two-dimensional ultrarelativistic plasma52,54,55, ∂ nˆ + c∗2∂ Pˆ = 0, (34) t E j j c∗ ˆ ˆ cs = √ . (42) ∂tPi + ∂jΠij = 0. (35) 2 Hydrodynamic equations at the Euler (ballistic) scale are These modes correspond to joint oscillations of obtained by averaging the operator equations of motion the full set of longitudinal hydrodynamic modes over locally thermal states, and here yield the constitu- {δn , δn, δe , δj , δj , δP }. We emphasize that bal- tive relation E k k E,k k listic behaviour of these modes is, strictly speaking, a ∗2 jE = c P (36) property of the field theory Eq. (29) and not the un- derlying lattice model; in realistic settings, such proxi- for the energy current. Restoring Navier-Stokes (dif- mate Lorentz invariance yields a vanishing Drude weight fusive) corrections and passing to the linear response and an anomalously large thermal conductivity48,49. It is regime, it follows that the linear-response energy current nevertheless consistent to treat these modes as ballistic picks up a reactive coupling to momentum because the crossover timescale to diffusive behaviour, which is the momentum relaxation time τ , exceeds the ∗2 0 P δjE = c δP + σ δe − DEn∇δn − DEE∇δnE, (37) IR cutoff of our hydrodynamic theory in Eq. (31), at least for clean systems in which τ ∼ eEU /kB T  τ ∼ T −1. while the linear-response momentum current is given by P ψ Finally there is a diffusive, transverse momentum mode  ∂p ∂p  due to the shear viscosity, with dispersion relation δΠij = δn + δnE δij ∂n ∂nE ⊥ 2 iωP = νk , (43) − ν(∂iδPj + ∂jδPi − δij∇ · δP) − ζδij∇ · δP. (38) that is decoupled from the other hydrodynamic modes. where p(n, nE) is the equilibrium pressure and ν and ζ define kinematic shear and bulk viscosities respectively. We now turn to the twelve flavour modes that arise The generic form of the other constitutive relations, Eqs. from SU(4) symmetry, excluding the three spin modes (22) and (23), is unchanged. considered above. As noted above, these twelve modes We deduce that of the five hydrodynamic modes transform nontrivially under elements of the lattice space present for the U(1) spinon Fermi surface state, four are group, which prevents them from coupling linearly to any unchanged for the algebraic spin liquid, namely the three of the hydrodynamic modes considered so far. We first components of spin and the transverse photon mode define flavour charge densities

2 i † i iωSα = DSSk , α = 1, 2, 3, (39) nF = hΨ µ Ψi, (44) 2 4 i,α † α i iωa = Dak + O(k ). (40) nF = hΨ σ µ Ψi, (45) 8 through local equilibrium averages of local observables, arising from internal non-Abelian symmetries and diffu- i i,α 25 and associated current densities jF , jF , which satisfy sive hydrodynamic modes was derived recently as a pre- diction of mode-coupling theory, and is consistent with i i ∂tnF + ∇ · jF = 0, i = 1, 2, 3, (46) our treatment of internal symmetries above. To summarize, the SU(4)-symmetric algebraic spin liq- ∂ ni,α + ∇ · ji,α = 0, i, α = 1, 2, 3. (47) t F F uid state has nineteen hydrodynamic modes, of which fifteen arise from the internal SU(4) flavour symmetry These conserved densities correspond to competing or- (this includes the three components of spin), one is the der parameters of the algebraic spin liquid phase, each of transverse photon mode discussed above for the spinon which arises from a generator of the SU(4) flavour sym- Fermi surface state, one is the diffusive transverse mo- 31 1/2 metry . For example, nF acts as an order parameter mentum mode, and two are ballistic, longitudinal energy- 3 for valence-bond-solid (VBS) order, while nF is a local momentum modes arising from emergent Lorentz sym- proxy for the density of N´eelSkyrmions. metry. Now consider the linear-response hydrodynamic be- haviour about a fully disordered equilibrium state, in i i,α which nF,0 = nF,0 = 0 and all such competing orders V. DISCUSSION are on equal footing. The linear response of the system is stratified into symmetry sectors by the action of the lattice space group, yielding general constitutive relations Arguing based on general principles of symmetry and linear response, we have derived the linear-response hy- δj1 = −D ∇δn1 − D ∇δn2 , (48) drodynamics of two representative quantum spin liquid F 11 F 12 F states. We find that the intricate low energy physics of 2 1 2 δjF = −D21∇δnF − D22∇δnF , (49) such states, as made manifest by the gauge and flavour 3 3 symmetries of their effective field theories, leads to a va- δjF = −D33∇δnF , (50) 1,α 0 1,α 0 2,α riety of hydrodynamic modes beyond the quartet of dif- δjF = −D11∇δnF − D12∇δnF , (51) fusive spin and energy modes that characterize ordinary 2,α 0 1,α 0 2,α δjF = −D21∇δnF − D22∇δnF , (52) paramagnets. Possible experimental signatures of nonlinear spinon δj3,α = −D0 ∇δn3,α. (53) F 33 F physics have been the focus of recent theoretical atten- tion56,57. Above, we identified a simple and generic fea- We have already exploited spin rotation symmetry of ture of spin liquids’ linear response, namely the presence the background equilibrium state to reduce the num- of a diffusive transverse mode of the internal gauge pho- ber of independent linear-response coefficients to ten. ton. Physically speaking, our results on the hydrody- In principle, the system Eqs. (48)-(53) is sensitive to namic modes arising from gauge and flavour symmetry VBS-type order of the background equilibrium state. If yield predictions for the finite-temperature dynamics of such order is absent, as we have assumed, then the dif- various autocorrelation functions of local two-, three- and D D  D0 D0  fusion tensors Dˆ = 11 12 , Dˆ 0 = 11 12 in four-spin operators31,35 in the “hydrodynamic” range of D D D0 D0 21 22 21 22 frequencies and wavenumbers depicted in Fig. 1. While the VBS sector cannot pick out a preferred direction the existence of such a hydrodynamic regime for spin in the space of modes, yielding additional constraints liquids is relatively easy to justify in three spatial dimen- D = D = 0, D = D , D0 = D0 = 0 and 21 12 11 22 12 21 sions, for which compactness of the emergent gauge field D0 = D0 , leaving only four independent transport co- 11 22 does not generically destabilize the spin liquid state28, its efficients, {D ,D ,D0 ,D0 }. The resulting hydrody- 11 33 11 33 existence for physically realistic examples of spin liquid namic normal modes are diffusive, and tabulated in Table states in two dimensions is unclear. Pending a more rig- I. orous understanding of the stability of two-dimensional Mode Dispersion Multiplicity spin liquid states, the arguments summarized in Fig. 1 1 2 2 predict a robust hydrodynamic regime at small T > 0 for δnF , δnF ω = −iD11k 2 3 2 the two-dimensional Dirac spin liquid with a large num- δnF ω = −iD33k 1 ber of fermion flavours28,31, and possibly even for the 1,α 2,α 0 2 30,46 δnF , δnF ω = −iD11k 6 U(1) spin liquid with a Fermi surface . 3,α 0 2 δnF ω = −iD33k 3 Our study of the hydrodynamics of quantum spin liq- uids connects to an ongoing theoretical effort to under- TABLE I. Hydrodynamic normal modes of the system Eqs. stand the implications of Lie group symmetry within (48)-(53). hydrodynamics. Various anomalous transport phenom- ena in one-dimensional quantum systems can be traced These multiplicities are in agreement with the lattice to the presence of Lie group symmetry13,15,58–60, and it symmetry classification of the SU(4) flavour generators would be interesting to understand how far ideas useful that was given in previous work31. We also note that a in understanding spin liquid physics, for example, one- one-to-one correspondence between the Noether charges dimensional versions of parton constructions61,62, might 9 shed light on these phenomena. In the context of effec- systems with gauge symmetries. tive field theory approaches to hydrodynamics24, the spin Acknowledgments. We thank Igor Herbut and espe- liquid states studied in this paper motivate an extension cially Andrew Lucas and Joseph Maciejko for several il- of the hydrodynamic formalism recently developed for luminating discussions. SLS acknowledges support from physical systems with global Lie group symmetries25 to the United States Department of Energy via grant No. DE-SC0016244.

[1] E. M. Lifshitz and L. P. Pitaevskii, Statistical physics: teracting 1d bose gases (2020), 2009.06651. theory of the condensed state, vol. 9 (Elsevier, 2013). [18] A. Scheie, N. E. Sherman, M. Dupont, S. E. Nagler, M. B. [2] R. Shankar, Rev. Mod. Phys. 66, 129 (1994), URL Stone, G. E. Granroth, J. E. Moore, and D. A. Tennant, https://link.aps.org/doi/10.1103/RevModPhys.66. Nature Physics (2021), ISSN 1745-2481, URL https:// 129. doi.org/10.1038/s41567-021-01191-6. [3] E. B. Kolomeisky, T. J. Newman, J. P. Straley, and X. Qi, [19] G. Policastro, D. T. Son, and A. O. Starinets, Journal Phys. Rev. Lett. 85, 1146 (2000), URL https://link. of High Energy Physics 2002, 043 (2002), URL https: aps.org/doi/10.1103/PhysRevLett.85.1146. //doi.org/10.1088/1126-6708/2002/09/043. [4] L. Erd˝os,B. Schlein, and H.-T. Yau, Phys. Rev. Lett. [20] P. K. Kovtun, D. T. Son, and A. O. Starinets, Phys. 98, 040404 (2007), URL https://link.aps.org/doi/ Rev. Lett. 94, 111601 (2005), URL https://link.aps. 10.1103/PhysRevLett.98.040404. org/doi/10.1103/PhysRevLett.94.111601. [5] R. Gurzhi, J Exp Theor Phys 17, 521 (1963). [21] C. W. von Keyserlingk, T. Rakovszky, F. Pollmann, and [6] Z. Z. Yu, M. Haerle, J. W. Zwart, J. Bass, W. P. S. L. Sondhi, Phys. Rev. X 8, 021013 (2018), URL https: Pratt, and P. A. Schroeder, Phys. Rev. Lett. 52, //link.aps.org/doi/10.1103/PhysRevX.8.021013. 368 (1984), URL https://link.aps.org/doi/10.1103/ [22] A. Nahum, S. Vijay, and J. Haah, Phys. Rev. X 8, PhysRevLett.52.368. 021014 (2018), URL https://link.aps.org/doi/10. [7] L. Molenkamp and M. de Jong, Solid-State Elec- 1103/PhysRevX.8.021014. tronics 37, 551 (1994), ISSN 0038-1101, URL [23] V. Khemani, A. Vishwanath, and D. A. Huse, Phys. Rev. https://www.sciencedirect.com/science/article/ X 8, 031057 (2018), URL https://link.aps.org/doi/ pii/0038110194902445. 10.1103/PhysRevX.8.031057. [8] D. A. Bandurin, I. Torre, R. K. Kumar, M. Ben Shalom, [24] M. Crossley, P. Glorioso, and H. Liu, Journal of High A. Tomadin, A. Principi, G. H. Auton, E. Khestanova, Energy Physics 2017, 95 (2017), ISSN 1029-8479, URL K. S. Novoselov, I. V. Grigorieva, et al., Science 351, https://doi.org/10.1007/JHEP09(2017)095. 1055 (2016), ISSN 0036-8075, URL https://science. [25] P. Glorioso, L. V. Delacr´etaz,X. Chen, R. M. Nandk- sciencemag.org/content/351/6277/1055. ishore, and A. Lucas, SciPost Phys. 10, 15 (2021), URL [9] P. J. W. Moll, P. Kushwaha, N. Nandi, B. Schmidt, https://scipost.org/10.21468/SciPostPhys.10.1. and A. P. Mackenzie, Science 351, 1061 (2016), 015. ISSN 0036-8075, URL https://science.sciencemag. [26] L. Savary and L. Balents, Reports on Progress in Physics org/content/351/6277/1061. 80, 016502 (2016), URL https://doi.org/10.1088% [10] O. A. Castro-Alvaredo, B. Doyon, and T. Yoshimura, 2F0034-4885%2F80%2F1%2F016502. Phys. Rev. X 6, 041065 (2016), URL https://link.aps. [27] Y. Zhou, K. Kanoda, and T.-K. Ng, Rev. Mod. Phys. org/doi/10.1103/PhysRevX.6.041065. 89, 025003 (2017), URL https://link.aps.org/doi/ [11] B. Bertini, M. Collura, J. De Nardis, and M. Fagotti, 10.1103/RevModPhys.89.025003. Phys. Rev. Lett. 117, 207201 (2016), URL https:// [28] M. Hermele, T. Senthil, M. P. A. Fisher, P. A. link.aps.org/doi/10.1103/PhysRevLett.117.207201. Lee, N. Nagaosa, and X.-G. Wen, Phys. Rev. B 70, [12] V. B. Bulchandani, R. Vasseur, C. Karrasch, and J. E. 214437 (2004), URL https://link.aps.org/doi/10. Moore, Phys. Rev. Lett. 119, 220604 (2017), URL 1103/PhysRevB.70.214437. https://link.aps.org/doi/10.1103/PhysRevLett. [29] Note1, for this reason it would be simpler to consider the 119.220604. hydrodynamic regime for quantum spin liquids in three [13]M. Znidariˇc,Phys.ˇ Rev. Lett. 106, 220601 (2011), URL dimensions, where topological phases can persist to non- https://link.aps.org/doi/10.1103/PhysRevLett. zero temperatures. 106.220601. [30] S.-S. Lee, Phys. Rev. B 78, 085129 (2008), URL https: [14] M. Ljubotina, M. Znidariˇc,andˇ T. Prosen, Phys. Rev. //link.aps.org/doi/10.1103/PhysRevB.78.085129. Lett. 122, 210602 (2019), URL https://link.aps.org/ [31] M. Hermele, T. Senthil, and M. P. A. Fisher, Phys. Rev. doi/10.1103/PhysRevLett.122.210602. B 72, 104404 (2005), URL https://link.aps.org/doi/ [15] M. Dupont and J. E. Moore, Phys. Rev. B 101, 10.1103/PhysRevB.72.104404. 121106 (2020), URL https://link.aps.org/doi/10. [32] J. D. Jackson, Classical electrodynamics (1999). 1103/PhysRevB.101.121106. [33] M. Baggioli, M. Vasin, V. Brazhkin, and K. Tra- [16] M. Schemmer, I. Bouchoule, B. Doyon, and J. Dubail, chenko, Physics Reports 865, 1 (2020), ISSN 0370- Phys. Rev. Lett. 122, 090601 (2019), URL https:// 1573, URL http://www.sciencedirect.com/science/ link.aps.org/doi/10.1103/PhysRevLett.122.090601. article/pii/S0370157320301198. [17] N. Malvania, Y. Zhang, Y. Le, J. Dubail, M. Rigol, and [34] A. A. Abrikosov, Physics Physique Fizika 2, 5 D. S. Weiss, Generalized hydrodynamics in strongly in- (1965), URL https://link.aps.org/doi/10.1103/ 10

PhysicsPhysiqueFizika.2.5. 68, 104401 (2003), URL https://link.aps.org/doi/ [35] X. G. Wen, F. Wilczek, and A. Zee, Phys. Rev. B 10.1103/PhysRevB.68.104401. 39, 11413 (1989), URL https://link.aps.org/doi/10. [50] S. Hartnoll, A. Lucas, and S. Sachdev, Holographic 1103/PhysRevB.39.11413. Quantum Matter, The MIT Press (MIT Press, 2018), [36] B. S. Shastry and B. I. Shraiman, Phys. Rev. Lett. ISBN 9780262348027, URL https://books.google. 65, 1068 (1990), URL https://link.aps.org/doi/10. com/books?id=DW9RDwAAQBAJ. 1103/PhysRevLett.65.1068. [51] K. Damle and S. Sachdev, Phys. Rev. B 56, 8714 (1997), [37] W.-H. Ko and P. A. Lee, Phys. Rev. B 84, URL https://link.aps.org/doi/10.1103/PhysRevB. 125102 (2011), URL https://link.aps.org/doi/10. 56.8714. 1103/PhysRevB.84.125102. [52] T. V. Phan, J. C. W. Song, and L. S. Levitov, Ballis- [38] P. A. Lee and N. Nagaosa, Phys. Rev. B 87, tic heat transfer and energy waves in an electron system 064423 (2013), URL https://link.aps.org/doi/10. (2013), 1306.4972. 1103/PhysRevB.87.064423. [53] A. Lucas, Phys. Rev. B 93, 245153 (2016), URL https: [39] S. Mukerjee, V. Oganesyan, and D. Huse, Phys. Rev. //link.aps.org/doi/10.1103/PhysRevB.93.245153. B 73, 035113 (2006), URL https://link.aps.org/doi/ [54] P. Kovtun, Journal of Physics A: Mathematical and The- 10.1103/PhysRevB.73.035113. oretical 45, 473001 (2012), URL https://doi.org/10. [40] H. Spohn, in Thermal Transport in Low Dimensions 1088/1751-8113/45/47/473001. (Springer, 2016), pp. 107–158. [55] A. Lucas and K. C. Fong, Journal of Physics: Condensed [41] P. A. Lee and N. Nagaosa, Phys. Rev. B 46, Matter 30, 053001 (2018), URL https://doi.org/10. 5621 (1992), URL https://link.aps.org/doi/10. 1088/1361-648x/aaa274. 1103/PhysRevB.46.5621. [56] L. Balents and O. A. Starykh, Phys. Rev. B 101, [42] C. P. Nave and P. A. Lee, Phys. Rev. B 76, 020401 (2020), URL https://link.aps.org/doi/10. 235124 (2007), URL https://link.aps.org/doi/10. 1103/PhysRevB.101.020401. 1103/PhysRevB.76.235124. [57] S. C. Morampudi, F. Wilczek, and C. R. Laumann, Phys. [43] M. Hermele, M. P. A. Fisher, and L. Balents, Phys. Rev. Rev. Lett. 124, 097204 (2020), URL https://link.aps. B 69, 064404 (2004), URL https://link.aps.org/doi/ org/doi/10.1103/PhysRevLett.124.097204. 10.1103/PhysRevB.69.064404. [58] E. Ilievski, J. De Nardis, M. Medenjak, and T. Prosen, [44] T. Giamarchi and O. U. Press, Quantum Physics in Phys. Rev. Lett. 121, 230602 (2018), URL https:// One Dimension, International Series of Monogr (Claren- link.aps.org/doi/10.1103/PhysRevLett.121.230602. don Press, 2004), ISBN 9780198525004, URL https: [59] V. B. Bulchandani, Phys. Rev. B 101, 041411 (2020), //books.google.com/books?id=1MwTDAAAQBAJ. URL https://link.aps.org/doi/10.1103/PhysRevB. [45] A. Polyakov, Nuclear Physics B 120, 429 (1977), 101.041411. ISSN 0550-3213, URL http://www.sciencedirect.com/ [60] E. Ilievski, J. D. Nardis, S. Gopalakrishnan, R. Vasseur, science/article/pii/0550321377900864. and B. Ware, Superuniversality of superdiffusion (2020), [46] K.-S. Kim, Phys. Rev. B 72, 245106 (2005), URL https: 2009.08425. //link.aps.org/doi/10.1103/PhysRevB.72.245106. [61] C. Mudry and E. Fradkin, Phys. Rev. B 49, 5200 (1994), [47] I. F. Herbut, B. H. Seradjeh, S. Sachdev, and G. Murthy, URL https://link.aps.org/doi/10.1103/PhysRevB. Phys. Rev. B 68, 195110 (2003), URL https://link. 49.5200. aps.org/doi/10.1103/PhysRevB.68.195110. [62] C. Mudry and E. Fradkin, Phys. Rev. B 50, [48] A. Rosch and N. Andrei, Phys. Rev. Lett. 85, 11409 (1994), URL https://link.aps.org/doi/10. 1092 (2000), URL https://link.aps.org/doi/10. 1103/PhysRevB.50.11409. 1103/PhysRevLett.85.1092. [49] E. Shimshoni, N. Andrei, and A. Rosch, Phys. Rev. B 11

Appendix A: Linear-response hydrodynamics of quantum spin liquids

1. Hydrodynamic modes of the U(1) spin liquid with a Fermi surface

a. Hydrodynamic equations

We first Fourier transform in space and time. It is helpful to write v = vk + v⊥ for vector components parallel and perpendicular to k. In this notation, the Fourier transformed Maxwell equations read

δbk = 0, (A1) c∗ δb = k × δe , (A2) ⊥ ω ⊥ ik · ek = 4πδn, (A3) 4π e = j , (A4) k iω k  ω2  4πiω k2 − δe = δj . (A5) c∗2 ⊥ c∗2 ⊥

Meanwhile, the constitutive relations for charge and energy current read

δj⊥ = σ⊥δe⊥, (A6)

δjk = σkδek − iDnEδnEk, (A7) 0 δjE,⊥ = σ⊥δe⊥, (A8) 0 δjE,k = σkδek − iDEEδnEk, (A9) with transverse and longitudinal conductivities given by

1 σ = σ, σ = σ + D k2, (A10) ⊥ k 4π nn 1 σ0 = σ0, σ0 = σ0 + D k2. (A11) ⊥ k 4π En

b. Transverse hydrodynamic modes

There is a single transverse hydrodynamic mode corresponding to the transverse photon of the internal gauge field. By the constitutive relations given in Eqs. (A6)-(A9) and Eq. (A5), we deduce that transverse electric field excitations satisfy the equation

 ω2  4πi  k2 − 1 + σ δe = 0, (A12) c∗2 ω ⊥ ⊥

Eq. (A12) is the vanishing of the dielectric function familiar in ordinary conductors, and we deduce that a transverse photon propagates when the condition

4πiσ c∗2k2 1 + − = 0 (A13) ω ω2 is satisfied, which has roots q ω = −2πiσ ± c∗2k2 − (2πσ)2. (A14)

Only the positive branch of this equation is hydrodynamic (i.e. ω → 0 as k → 0), and yields the dispersion relation q ∗2 2 2 ωa = −2πiσ + c k − (2πσ) . (A15) 12

This dispersion relation exhibits the well-known peculiarity of a “k-gap”, or skin effect33, due to its branch point at the critical wavenumber 2πσ k∗ = . (A16) c∗ Thus short-wavelength excitations k  k∗ look like ballistic photon modes with ω ∼ c∗k, while the hydrodynamic regime, k  k∗, exhibits diffusive behaviour with

2 4 ωa ∼ −iDak + O(k ), k → 0, (A17) where the photon diffusion constant

c∗2 D = . (A18) a 4πσ ∗ This diffusive transverse electric field mode induces magnetic fluctuations δb = (c /ω)k × δe⊥ by Amp`ere’slaw Eq. (A2). By the constitutive relations Eqs. (A6) and (A8) it is also accompanied by transverse charge and energy 0 currents, δj⊥ = σ⊥δe⊥ and δjE,⊥ = σ⊥δe⊥. Thus the transverse photon is a simultaneous excitation of the modes {δe⊥, δb, δj⊥, δjE,⊥}.

c. Longitudinal modes

With the aid of Gauss’s law Eq. (A3) and the constitutive relations Eqs. (A6)-(A9), the conservation laws for charge and energy density reduce to the coupled linear system

2 ! ! iω − 4πσk −DnEk δn 0 2 = 0. (A19) −4πσk iω − DEEk δnE

Longitudinal modes of charge and energy exist when the determinant condition

iω − 4πσ −D k2 k nE = 0. (A20) 0 2 −4πσk iω − DEEk is satisfied. The roots of this equation yield dispersion relations

 h  σ0  2 4 i −i 4πσ + Dnn + DnE k + O(k ) ω = σ (A21)  σ0  2 4 −i DEE − σ DnE k + O(k )

Notice that the charge mode decays at a rate 4πσ, and is not hydrodynamic as ω → 0; this is the rapid relaxation of free charge that is familiar in ordinary conductors. The energy mode is diffusive, just as in the paramagnetic phase. From  2  DnE k 4 Eq. (A19), it follows that the energy mode is accompanied by the charge fluctuation δn = − 4πσ + O(k ) δnE; by the constitutive relations (A7) and (A9), these give rise to fluctuations of longitudinal charge and energy current. Thus the diffusive energy mode is a simultaneous excitation of the modes {δnE, δn, δek, δjk, δjE,k}.

d. Summary

We deduce that the U(1) spin liquid with a Fermi surface has five diffusive hydrodynamic modes in total, given to leading order in k by (restoring magnetization modes)

2 iωSα = DSSk , α = 1, 2, 3, (A22)  σ0  iω = D − D k2, (A23) E EE σ nE 2 iωa = Dak , (A24) with linear response transport coefficients defined as above. 13

2. Hydrodynamic modes of the SU(4) algebraic spin liquid

The key difference compared to the hydrodynamics of the U(1) spin liquid with a Fermi surface, beyond the additional hydrodynamic conservation laws arising from flavour symmetry that were discussed in the main text, is the emergence of momentum as a hydrodynamic mode. In Fourier space, the linear-response momentum modes satisfy

2 (iω − νk )δP⊥ = 0, (A25)   2 ∂p ∂p (iω − (ν + ζ)k )δPk = ik δn + δnE . (A26) ∂n ∂nE

a. Transverse modes

The only change in the transverse mode structure compared to the spinon Fermi surface state is an additional transverse mode of momentum, whose dispersion is given by

2 iωP,⊥ = νk (A27)

This mode corresponds to a constant pressure excitation of {δP⊥} and is decoupled from the other hydrodynamic modes.

b. Longitudinal modes

The longitudinal slow modes of the algebraic spin liquid consist of charge, energy and longitudinal momentum, ˆ δPk = k · δPk. Assuming an ultrarelativistic equation of state p = p(nE) = nE/2 at half-filling, the linear-response mode coupling equations are given by

 2    iω − 4πσk −DnEk 0 δn  −4πσ0 iω − D k2 −ic∗2k    = 0, (A28)  k EE  δnE 0 −ik/2 i(ω − ωP,k(k)) δPk

2 where ωP,k = (ν + ζ)k . To understand the effect of the coupling to longitudinal momentum compared to Eq. (A19), it is clearest to first discard dissipative terms (i.e. pass to the ballistic scaling limit), which yields     iω 0 0 δn  ∗2     0 iω −ic k δnE = 0. (A29) 0 −ik/2 iω δPk

The resulting determinant condition reads   ω ω2 − c∗2k2/2 = 0, (A30) from which we deduce that at the Euler scale, fluctuations of energy and longitudinal momentum combine to form ballistic “sound” modes, with speed c∗ cs = √ . (A31) 2 The full determinant condition resulting from Eq. (A28) is rather opaque, but diffusive broadening of the sound modes can be obtained perturbatively in k and reads

i  σ0  ω = ±c k − D − D + ν + ζ k2 + O(k3). (A32) s s 2 EE σ nE