arXiv:1905.03335v1 [cond-mat.mes-hall] 8 May 2019 osuperfluid to locity hsvoaecnevto of conservation violate thus bandwidth the the of of hand, energies order other with excitations the involving On processes scattering system. the of excitations energy by described be should hydrodynamics. flow two-fluid its and motion, macroscopic and a htmvswt t w velocity own its with a moves form that excitations gas elementary the temperature non-zero h unu iudhsafiiemomentum finite a has liquid particles. quantum left-moving the and right- of At numbers the of ence iscmoe ffrin,i hc case which in fermions, of composed con- uids the of meaning quantity physical served The follows. as understood hydrodynamics. two-fluid ftemperature of e,btas nadtoa quantity sys- additional the an of also energy but the and tem, only momentum, not particles, conserve of number perturbations irrelevant these by sarsl utne iud eaea ueflisand, superfluids as that superfluid behave [3–5] to liquids recently shown Luttinger similarly was result It a as [1]. Ref. in introduced vn etrain cl spwr feeg,adthus and rate energy, relaxation of corresponding powers irrel- the as The scale re- equilibrium. perturbations thermal to to evant rise system give the processes of scattering laxation resulting The group which renormalization sense. Hamiltonian, the in liquid perturbations irrelevant from Luttinger are arise the excitations to equi- the corrections not between does Interactions non-interacting system the librate. are approximation this excitations In this of bosons. elementary form so-called simplest the the the In theory of 2]. [1, framework theory liquid the Luttinger in described commonly utne iudter consol o h low- the for only accounts theory liquid Luttinger h esntesse eae sasprudcnbe can superfluid a as behaves system the reason The Introduction. w-uddnmc foedmninlqatmlqisi th in liquids quantum one-dimensional of dynamics Two-fluid J m vni h bec feeetr excitations elementary of absence the in even 0 6= u and 0 = p N bu he e rnpr offiins eoti expression obtain We coefficients. rate. the backscattering transport However, liq new quantum hydrodynamics. three such single-fluid about of theory to two-fluid reduces the theory one-dimensiona develop we of Here realizations sup ance. with experimental ap liquids most present at Luttinger whereas available endows description This hydrodynamic two-fluid particles. of tering F 4 J/mN e utne iudehbt w ye of types two exhibits liquid Luttinger He, r h asadnme fprils At particles. of number and mass the are utne iudter foedmninlqatmsystem quantum one-dimensional of theory liquid Luttinger n-iesoa unu iud are liquids quantum One-dimensional — T 1 J aeil cec iiin ron ainlLaboratory, National Argonne Division, Science Materials motnl,tecliin generated collisions the Importantly, . 2 sms rnprn o unu liq- quantum for transparent most is eateto hsc,Uiest fWsigo,Seattle, Washington, of University Physics, of Department Here . 4 D e[,7,cnb ecie by described be can 7], [6, He a akcte emosand fermions backscatter can p J F h otecetpro- efficient most The . steFrimomentum, Fermi the is τ ex − 1 u J .A Matveev A. K. clsa power a as scales ex hc a first was which , hs similar Thus, . J p stediffer- the is F J Dtd a ,2019) 8, May (Dated: n ve- and invariance 1 n .V Andreev V. A. and ad[] n hshv aeta sepnnilysmall the exponentially of is temperatures, that bottom low rate the a at have near thus holes and involve [8], band type this of cesses iud slmtdt h rqec range frequency the to limited quantum is one-dimensional liquids of behavior superfluid the that ierapoiaintedsiaieprsof parts dissipative the approximation linear httepril ubrcurrent number particle the that is regime this goal. in main liquids our quantum one-dimensional of the- Galilean ory hydrodynamic of two-fluid violation of to Development an leads invariance. plays usually lattice and solid role in crystalline important realized underlying liquids the quantum systems Galilean-invariant. for state is hand, other system the the On that approach assumption this the of limitation is important An dimension. one oetm ontaffect conserve not which do particles, the momentum, between collisions the fore, density mentum n h unu iudbhvsa omlfluid. normal a as behaves liquid quantum the and quencies essla oeulbaino h velocities, the of equilibration to lead cesses aduster fsprudt fliquid of the superfluidity adapting of by theory obtained was Landau’s which [3–5], theory namic inpootoa to proportional tion current energy the atceadeeg urnsapa nrsos ogra- to response potential in chemical appear of dients currents energy and particle to similar and, sity n eprtr rdet nteasneo Galilean of absence the In invariance conductivity gradient. thermal the temperature of and terms in expressed typically ttelws frequencies lowest the At r eae otegradients the to related are st xrs hs rnpr offiinsi em fthe of terms in coefficients time transport enables relaxation these theory two-fluid express Our to us [12]. conductors coef- in thermoelectric ficients to similar coefficients, transport four nipratpoet fGlla-nain ytm is systems Galilean-invariant of property important An yaiso n-iesoa unu iud tfre- at liquids quantum one-dimensional of Dynamics le nyt aienivratsystems, Galilean-invariant to only plies unu iud akGlla invari- Galilean lack liquids quantum l rudpoete.Tecorresponding The properties. erfluid goe xoetal ekbackscat- weak exponentially ignores s o hs offiinsi em fthe of terms in coefficients these for s ω bec fGlla nainebrings invariance Galilean of absence j is ntelwfeunylmtthe limit low-frequency the In uids. ≪ n ron,Ilni 03,USA 60439, Illinois Argonne, sntuiul endb oetmden- momentum by defined uniquely not is τ ahntn915 USA 98195, Washington ex − τ p 2 1 . iie ytepril mass particle the by divided sdsrbdb h w-udhydrody- two-fluid the by described is j j p ε ε a isptv opnn.The component. dissipative a has , hc nadto oacontribu- a to addition in which , otisadsiaiecontribution, dissipative a contains τ − 1 bec fGalilean of absence e ω ∝ j n ∂ ≪ µ hsi ncnrs with contrast in is This . x e µ − n temperature and τ D/T and − j 1 n akcteigpro- backscattering , sgvnb h mo- the by given is ∂ 91] eas of Because [9–11]. x T τ − yamti of matrix a by 1 4 u ≪ e[,7 to 7] [6, He ex m j − ω n There- . T u ≪ and 0 na In . → τ ex − j 0, 1 ε . 2

Luttinger liquid.—For simplicity, we consider a quan- In most applications the Hamiltonian of the Luttinger tum liquid of spinless particles. In the Luttinger liquid liquid is approximated by the first line of Eq. (5). In approximation it is described in terms of the bosonic this case the elementary excitations are noninteracting fields ϕ(x) and θ(x), satisfying the commutation relation bosons. The energy of the boson with momentum q is ′ ′ ǫq = v|q| [1, 2]. We will need to account for a correction [ϕ(x), ∂x′ θ(x )] = iπδ(x − x ). (1) to this expression arising at small but finite χ. To this The fields are subject to the boundary conditions end we substitute Eq. (3) into the cubic term in Eq. (5) proportional to αθ. To linear order in χ this generates a ϕ(L)= ϕ(0), θ(L)= θ(0) − πJ, (2) correction to the quadratic Hamiltonian in the first line of Eq. (5) proportional to the momentum of the elementary where L is the size of the system. The second condi- excitations. As a result, the energy spectrum takes the tion formally defines the conserved quantity J, which for form systems of fermions is the difference of the numbers of χ ǫ = v|q| + q. (7) right- and left-moving particles. It will be convenient to q m∗ also define a field ϑ(x), which satisfies periodic boundary conditions ϑ(L)= ϑ(0) and is related to θ(x) by Here we have introduced the effective mass ∗ π~ χ J m = . (8) ∂ θ = − + ∂ ϑ, χ = π~ . (3) ∂n(vK) x ~ x L In the Galilean-invariant case the parameters v and K are Here ~ is the Planck’s constant, while χ has the dimen- not independent. Specifically, vK = v where the Fermi sion of momentum and is proportional to the density of F velocity v = π~n/m. In this case the effective mass (8) the conserved quantity J. F coincides with the mass m of the particles comprising the The momentum of the Luttinger liquid can be ex- Luttinger liquid. Equation (7) can then be interpreted pressed in terms of the bosonic fields as [1, 13] as the Galilean transformation of the original spectrum

∂xϕ ~ ǫq = v|q| to the frame moving with velocity χ/m. P = −~ dx n + ∂xθ = Nχ − dx ∂xϕ ∂xϑ. Z  π  π Z Thermal equilibrium.—We now consider an equilib- (4) rium state of the Luttinger liquid. It is described by Here we denote the total number of particles by N, while the number of particles N, the value of J, and the boson n = N/L is the average density. The first form of Eq. (4) occupation numbers Nq, which take the usual Bose form is clear from the physical meaning of the bosonic fields: −1 ǫq − uq ∂xϕ(x)/π is the density of particles at point x, measured N = exp − 1 . (9) q   T   from n, whereas −~ ∂xθ(x) is the momentum of the liquid per particle. In the second form of Eq. (4) for a liquid The parameter u appears as a consequence of the con- of spinless fermions Nχ can be expressed in terms of the servation of momentum and can be thought of as the ~ Fermi momentum pF = π n as pF J and represents the velocity of the gas of excitations. It is worth noting that momentum of the ground state with unequal numbers the occupation numbers (9) implicitly depend on N and of the right- and left-moving particles. The remaining J, because the velocity of the excitations v in Eq. (7) is ~ integral of −( /π)∂xϕ ∂xϑ accounts for the momentum a function of density, while χ ∝ J. of the elementary excitations. Using the distribution function (9) it is straightforward The Hamiltonian of the Luttinger liquid can be pre- to obtain the thermodynamic properties of the system to sented in the form [1, 2] leading order in T/D. Below we will need expressions ~ for the energy and momentum densities of the Luttinger v 2 1 2 H = dx K(∂xθ) + (∂xϕ) liquid. We will limit our consideration to low-velocity Z 2π  K  flows and neglect contributions beyond linear order in 2 3 u and χ. The energy density ε is obtained by adding +αθ∂xϕ(∂xθ) + αϕ(∂xϕ) + ... , (5)  contributions ǫqNq of all excitations, resulting in where v is the velocity of the elementary excitations and πT 2 K is the Luttinger liquid constant. The coefficients α ε = . (10) θ 6~v and αϕ can be expressed in terms of the dependences of v and K on density as [11] We obtain the momentum density p using the second form of Eq. (4), in which the second term is given by the ~ ~ v sum of the contributions qNq of all bosonic states. This αθ = 2 ∂n (vK) , αϕ = 2 ∂n . (6) 2π 6π K  yields We omitted from the Hamiltonian (5) terms with scaling χ πT 2 p = nχ + ρ u − , ρ = . (11) dimensions exceeding three. ex ∗ ex ~ 3  m  3 v 3

We earlier interpreted u as the velocity of the gas of The evaluation of the momentum current jp starts with bosonic excitations of the Luttinger liquid. Given that the equation of motion for the operator of momentum the corresponding contribution to the momentum den- density −~(n + ∂xϕ/π)∂xθ, cf. Eq. (4). The result is sity is ρexu, one can interpret ρex as the mass density of 2 ∂n(nv) the gas of excitations. jp = m0v δn + ε . (20) Hydrodynamics of Luttinger liquids.—Let us now con- v sider the dynamics of the Luttinger liquid at time scales Finally, the current jχ is obtained by averaging the equa- much longer than τex. In this regime the liquid is in local tion of motion for the operator ∂xθ and using Eq. (3). thermal equilibrium at every point. As a result, its state This yields is fully described by the densities of the four conserved 2 quantities: number of particles, energy, momentum, and m0v ∂nv jχ = δn + ε . (21) J. Time evolution of these densities is described by four n v continuity equations expressing the conservation laws: It is easy to show that the above expressions for jp and j coincide with the pressure Π and chemical potential µ ∂tn + ∂xjn =0, (12) χ of the Luttinger liquid, respectively. ∂ ε + ∂ j =0, (13) t x ε Compared to conventional hydrodynamics, the above ∂tp + ∂xjp =0, (14) theory includes an additional equation (15) originating ∂tχ + ∂xjχ =0, (15) from the conservation of J. In the Galilean-invariant case the resulting theory was interpreted [3–5] as two-fluid hy- where we used χ defined in (3) instead of the density J/L. drodynamics, fully analogous to that of superfluid 4He In order to complete the hydrodynamic description of the in three dimensions [6, 7]. The normal component of the system we need to evaluate the currents jn, jε, jp, and fluid in this analogy is the gas of excitations characterized jχ. We will restrict ourselves to terms up to linear order by velocity u, while the superfluid component describes in small parameters χ(x), u(x), and δn(x)= n(x)−N/L. the Fermi surface and moves with velocity χ/m. In the To evaluate the particle current jn we write the Heisen- absence of Galilean invariance the full set of two-fluid berg equation of motion for the operator of particle den- hydrodynamic equations of the Luttinger liquid is given sity n(x) = N/L + ∂xϕ(x)/π. This yields the evolution by the four continuity equations (12)–(15) along with the equation (12) in the operator form with constitutive relations (18)–(21). The latter were derived i vK 2α to linear order in the deviations from the static equilib- j (x)= [ϕ(x),H]= − ∂ θ − θ ∂ ϕ ∂ θ. (16) n π~ π x ~ x x rium state, in which δn, χ and u vanish. To this approx- imation the system is described by two effective masses, Next we substitute ∂xθ in the form (3) and evaluate the ∗ m and m0. The Galilean invariant limit is recovered at expectation value of the operator (16) in the equilibrium ∗ m = m0 = m. state of the Luttinger liquid. Introducing a new effective Hydrodynamics of one-dimensional quantum liquids.— mass As discussed above, Luttinger liquid approximation ne- π~n glects the exponentially weak backscattering processes, m = (17) 0 vK which relax J to its equilibrium value. The latter is determined by the velocity u and can be found as fol- ~ and noting that the operator (− /π)∂xϕ∂xϑ is the mo- lows. To leading order at T → 0 the value of J is de- mentum density of the bosonic excitations (see above), −(E−uP )/T termined by the Gibbs distribution wJ ∝ e , we obtain where E = π~vKJ 2/2L is the ground state energy of the nχ ρex χ Luttinger liquid, obtained by substituting Eq. (3) into the jn = + ∗ u − ∗ . (18) m0 m  m  first line of Eq. (5), and the momentum P = π~NJ/L is obtained from Eq. (11). The maximum of this dis- For Galilean-invariant systems vK = vF , and thus m0 = tribution gives the equilibrium value J = uN/vK or, m. In this case jn must coincide with the ratio of mo- mentum density p and the particle mass m. This is easily equivalently, χ = m0u. ∗ At time scales much longer than τex the deviation of verified by substituting m0 = m = m into Eqs. (11) and (18). the system from equilibrium is described by χ−m0u 6= 0. The above procedure can be extended to the evalua- For small deviations, relaxation to equilibrium proceeds with the exponentially small rate τ −1 [9–11], following tion of the remaining three currents. To obtain jε one can write the equation of motion for the operator of Hamilto- the usual relaxation law nian density, keeping both the quadratic and cubic terms d 1 (χ − m u)= − (χ − m u). (22) in Eq. (5). After thermal averaging one obtains dt 0 τ 0

n∂nv χ The relaxation processes obey the remaining three con- jε =2εu + ε . (19) v m0 servation laws. In particular, the momentum density is 4 conserved, dp/dt = 0. Using the expression (11), we ob- expressing them instead in terms of the momentum den- tain a linear relation between the dχ/dt and du/dt. This sity (11) and δχ = χ − m0u. Due to conservation of enables us to find dχ/dt. We then find the time evolution momentum, the parts of jn and jε proportional to p are equation for χ by substituting dχ/dt into the right-hand not affected by the relaxation processes and represent side of Eq. (15), the equilibrium contributions to the currents. The com- 1 ρ ponents of the currents proportional to δχ correspond to ∂ χ + ∂ j = − ex (χ − m u). (23) t x χ τ m n 0 dissipative contributions. In the low-frequency limit, the 0 latter can be expressed in terms of the gradients of tem- Here we accounted only for the leading at small T con- perature and chemical potential with the aid of Eq. (24). tribution in the right-hand side. This yields Equation (23) completes our generalizaiton of two-fluid hydrodynamics of one-dimensional quantum liquids to p ∂xT jn = − γ11∂xµ − γ12 , (25) the non-Galilean-invariant case. The dynamics of the m0 T − 1 2 system at frequencies ω ≪ τex is fully described by p ∂ (n v) ∂ T j = n ε − γ ∂ µ − γ x , (26) Eqs. (12)–(14) and (23) along with the constitutive rela- ε m n2v 21 x 22 T −1 −1 0 tions (18)–(21). At τ ≪ ω ≪ τex the right-hand side of Eq. (23) can be neglected and the system is described where by Eqs. (12)–(15). In this regime the system supports 2 1 1 2 two sound modes, similar to the Galilean-invariant case γ11 = τρex ∗ − , γ22 =2τεv , − m m  [3–5]. On the other hand, at low frequencies ω ≪ τ 1 0 qualitatively new physics arises. 1 1 γ12 = γ21 =2τε ∗ − . Dissipation in non-Galilean-invariant fluids.—In the m m0  low-frequency limit the backscattering processes are very Hydrodynamics of one-dimensional quantum liquids at effective at bringing the system to equilibrium charac- − ω ≪ τ 1 is fully described by the set of three continuity terized by only three conserved quantities: number of equations (12)–(14) along with the expressions (20), (25), particles, energy, and momentum. The quantity J is no and (26) for the currents. longer conserved; in a uniform fluid it takes the equilib- ∗ In the Galilean-invariant case, m = m = m, the dis- rium value uN/vK. In a non-uniform fluid, the system 0 sipative contribution to the particle current vanishes, and approaches a local equilibrium at every point in space, Eq. (25) recovers the expected relation j = p/m. The with χ(x) = m u(x). However, full local equilibrium n 0 energy current contains both the equilibrium and dissi- cannot be achieved in a non-uniform system, i.e., the de- pative contributions, with the latter defining the thermal viation from equilibrium δχ = χ − m u acquires a finite 0 conductivity of the liquid κ = γ /T , cf. Ref. [14]. In value proportional to the gradients of physical parame- 22 the absence of Galilean invariance the dissipation is de- ters of the system, such as the temperature and chemical scribed by the matrix of four coefficients γ . The latter potential. As a result, the system behaves as a conven- ij relates dissipative components of j and j to ∂ µ and tional fluid described the hydrodynamic equations (12)– n ε x ∂ T/T , and is analogous to the martix of thermoelectric (14), with the constitutive relations containing dissipa- x coefficients in conductors [12]. The coefficient γ is anal- tive contributions proportional to ∂ T and ∂ µ. 11 x x ogous to electrical conductivity. The off-diagonal martix We now apply our two-fluid description to evaluate dis- elements satisfy the Onsager relation, γ = γ . sipative components of the currents. We start by using 12 21 Summary.—We have developed a hydrodynamic the- Eq. (11) to bring δχ to the form δχ = (m0/ρex){[n + −1 ∗−1 ory of non-Galilean-invariant one-dimensional quantum ρ (m − m )]χ − p}. We then use Eqs. (14) and − ex 0 liquids, which applies at frequencies ω ≪ τ 1. The fluid (23) to express ∂ δχ as a linear combination of δχ, ex t is described by the three usual continuity equations (12)– ∂ j = ∂ Π, and ∂ j = ∂ µ. At low frequencies the x p x x χ x (14), expressing conservation of particle number, energy, time derivative ∂ δχ ∼ ωδχ ≪ δχ/τ and can be ne- t and momentum, and the additional evolution equation glected. This enables us to express δχ as (23). At small deviations from equilibrium the dynamics m s m ∗ δχ = τ 0 ∂ T − τ 1 − 0 ∂ µ. (24) of the system is characterized by two effective masses, m ρ x m∗ x ex   and m0, both of which coincide with the particle mass in Here we kept only the leading at T → 0 terms and used the Galilean-invariant case. At the lowest frequencies, −1 the thermodynamic relation ∂xΠ = n∂xµ + s∂xT to ex- ω ≪ τ , our theory reduces to ordinary hydrodynam- press the gradient of pressure in terms of the gradients ics. However, the absence of Galilean invariance results of temperature and chemical potential, with s = πT/3~v in the emergence of new transport coefficients, analogous being the entropy density of the Luttinger liquid. to conductivity and thermoelectric coefficients. The particle and energy currents given by Eqs. (18) The authors are grateful to I. L. Aleiner and W. De- and (19) are linear in u and χ. One can get further in- Gottardi for insightful discussions. Work at Argonne Na- sight into the physics associated with these currents by tional Laboratory was supported by the U.S. Department 5 of Energy, Office of Science, Materials Sciences and Engi- Phys. Rev. B 98, 155441 (2018). neering Division. Work at the University of Washington [6] L. D. Landau, The Theory of Superfluidity of Helium II, was supported by the U.S. Department of Energy Of- J. Phys. USSR 5, 71 (1941). fice of Science, Basic Energy Sciences under Award No. [7] I. M. Khalatnikov, An introduction to the theory of su- perfluidity (Perseus, New York, 2000). DE-FG02-07ER46452. [8] A. M. Lunde, K. Flensberg, and L. I. Glaz- man, Three-particle collisions in quantum wires: Corrections to thermopower and conductance, Phys. Rev. B 75, 245418 (2007). [9] T. Micklitz, J. Rech, and K. A. Matveev, Trans- [1] F. D. M. Haldane, ’Luttinger liquid theory’ of port properties of partially equilibrated quantum wires, one-dimensional quantum fluids. I. Properties Phys. Rev. B 81, 115313 (2010). of the Luttinger model and their extension to [10] K. A. Matveev and A. V. Andreev, Equi- the general 1d interacting spinless Fermi gas, libration of a spinless Luttinger liquid, J. Phys. C: Solid State Phys. 14, 2585 (1981). Phys. Rev. B 85, 041102(R) (2012). [2] T. Giamarchi, Quantum physics in one dimension [11] K. A. Matveev and A. V. Andreev, Scattering of hole (Clarendon, Oxford, 2004). excitations in a one-dimensional spinless quantum liquid, [3] K. A. Matveev and A. V. Andreev, Second Phys. Rev. B 86, 045136 (2012). Sound in Systems of One-Dimensional Fermions, [12] A. A. Abrikosov, Fundamentals of the Theory of Metals Phys. Rev. Lett. 119, 266801 (2017). (North-Holland, Amsterdam, 1988). [4] K. A. Matveev and A. V. Andreev, Hybrid Sound [13] Henceforth we assume normal ordering of bosonic fields. Modes in One-Dimensional Quantum Liquids, [14] W. DeGottardi and K. A. Matveev, Electrical and Phys. Rev. Lett. 121, 026803 (2018). Thermal Transport in Inhomogeneous Luttinger Liquids, [5] K. A. Matveev and A. V. Andreev, Propagation and at- Phys. Rev. Lett. 114, 236405 (2015). tenuation of sound in one-dimensional quantum liquids,