Arxiv:1909.08936V5

arXiv:1909.08936v5 [cond-mat.str-el] 23 Jun 2020 etLrnzsmer otelna re u h linear the but at order (1) linear the expansion to symmetry Lorentz gent h ee fteHmloin h hf fteWy node Weyl the of shift the Hamiltonian, p the of level the in ftecia nml nWy eiea n Weyl and semimetal Weyl e.g. in see anomaly applica- superfluids/superconductors, chiral the the For of tions quasiparticles. massless low-energy the in can potential fermions vector the coupling. electromagnetic minimal charged, the via are to fermions couple the addition If ever. h o-rva oriaedpnec trin related tetrads (torsion) the dependence to coordinate non-trivial the rudsaeietta ed ihdmnin funity of indices dimensions temporal with for fields tetrad spacetime ground urdb h aoia iesosof dimensions canonical the by quired utedet rsneo iesoflutavoe (UV) ultra-violet and Λ dimensionful a contentious scale of remained presence to has due anomaly subtle the to tribution anomaly (NY) Nieh-Yan lsia fields classical P tr nteoperator the in eters iersetu hnmn a emerge may phenomena with spectrum quasiparticles linear low-energy of description relativistic ls oteWy oes at node(s) Weyl the to close a,i he pta iesosa eei (two-fold) generic a momentum at at crossing dimensions band fermion spatial degenerate three in lar, aitna so h elform Weyl the of is Hamiltonian hr the where inso h emta al matrices Pauli Hermitean the of cients W hs akrudfilsipytecia nml for anomaly chiral the imply fields background These nnnrltvsi oooia atr ffcieyquasi- effectively matter, topological non-relativistic In { p csa neegn ail ag edwt emer- with field gauge (axial) emergent an as acts W } 2 sgn(det 1 o eprtr aoaoy eateto ple Physics Applied of Department Laboratory, Temperature Low ihcnncldmnin fmmnu,a re- as momentum, of dimensions canonical with e h yrdnmc fachiral a of hydrodynamics the ASnumbers: syste PACS Weyl ferm in Weyl transport relativistic thermal to to comparison anomalies low gravitational in emergent results with our parameters chemical discuss chiral fundamental the dimensionless and of corrections Fermi-liquid other e.g. from how show and discuss uefli re aaee.A itntfo h dimensionf the from distinct As parameter. order superfluid emty oooyadnme fcia unu ed nthe in fields quantum fund corrections, chiral low-temperature be with of may superfluid number Weyl parameter and dimensionless topology such geometry, that suggests This h ovninltrinlNe-a nml,teparamete the anomaly, Nieh-Yan torsional conventional the nWy ueflis where superfluids, Weyl in a i H edsustepsiiiyo h osoa ihYnanomal Nieh-Yan torsional the of possibility the discuss We e e e a µ .INTRODUCTION I. a µ = a i W ( p ( vanishes ) x x W n()cnla otegravitational the to lead can (1) in ) = = ) ∂ e p svlda uhlwrsae,how- scales, lower much at valid is σ i a 0 σ a H 2 ntemlNe-a nml nWy superfluids Weyl in anomaly Nieh-Yan thermal On a e aduIsiuefrTertclPyis 442Chernogo 142432 Physics, Theoretical for Institute Landau { n eoiyfrspatial for velocity and 8–17 e a i ( e a µ p ( a 0 p i∂ 5 ) e , eetees h Ycon- NY the Nevertheless, . o lwyvrigparam- varying slowly For . | − p µ a i W } p p ≡ W 1,3–5 r rmtdt back- to promoted are W r h iercoeffi- linear the are ) i∂ i h e chirality net The . + 1,5–7 t T − e · · · T µ a p stmeaueand temperature is σ 1,2 H 2 .Nissinen J. e,however, See, . nparticular, In . wv uefli,sc as such superfluid, -wave a em fsmlrfr n ftesm re ngains comi gradients, in order same the of and form similar of terms W , nparticu- In . h semi- the , a 1 = Dtd ue2,2020) 24, June (Dated: p e W a i , At . the , 2 1 (1) , 3, n ..Volovik G.E. and T T uecnutr swl,we h lcrmgei po- electromagnetic added the chiral is the when tential for well, apply as results The superconductors of transport superfluids. momentum Weyl tempera- hydrodynamic chiral the finite in different terms and ture anomaly NY gravitational below. discussion the and 15 Ref. ntedtiso h irsoi hsc,btaefully ef- are of but number species the physics, fermionic and fective microscopic topology depend geometry, the by not of determined do details i.e. the dimensionless, on are terms responding o l uhlweeg eprtr orcin propor- corrections to temperature tional low-energy anomaly. Nieh-Yan such torsional all momen- the For from the terms, originating to similar term other corrections a various among these include, that currents tum show We gauge. Yan i elfrin ihpstv n eaiebace at relativis- branches negative to momentum. and results zero positive these with fermions compare NY Weyl Weyl we tic general temperature in Finally finite universal similarly the systems. be that can term conjecture anomaly to us which anomaly superfluid, Nieh-Yan leads Weyl chiral quasirelativistic low-energy the the from in iden- contribution and superfluid, the chiral gradi- the tify order of free-energy lowest in the terms ent to corrections temperature finite where T h erdadspin-connection and tetrad the a ≪ a eew ics h eprtr orcin othe to corrections temperature the discuss we Here steeetv reegn oso rmthe from torsion emergent or effective the is o pctmswt oso adcraue,Ne and Nieh curvature), (and torsion with spacetimes For = 8–10 ∆ 3 at nvriy I006ALO Finland AALTO, FI-00076 University, Aalto , eA racia uecnutr ealso We superconductor. chiral a or He-A, 0 de oeta,cnas eepesdi terms in expressed be also can potential, osadtecneto ftetorsional the of connection the and ions ms. eso htsc emde xs in exist does term a such that show we , e r a a nrdcdte4dmninlinvariant 4-dimensional the introduced ftetype the of y eeg eaiitcfils aty we Lastly, fields. relativistic -energy γ T lutailt(V aaee Λ parameter (UV) ultraviolet ul + = 2 sdmninesi aoia units. canonical in dimensionless is mna,bigdtrie ythe by determined being amental, I OSOA ANOMALY TORSIONAL II. ω ,2 1, ntefe nry h rfcosi h cor- the in prefactors the energy, free the in ytm ycmaigti oa to this comparing By system. e a N b µ a ∧ dx = e v b µ s T and ok,Russia. lovka, stelcltta -omfil and field 1-form tetrad local the is → a ∂ T ∧ µ R 1,18 ( v ej ab s a 5 µ − pcfial,w opr the compare we Specifically, . − = ) = e e dω a A γT ω ∧ /m b a a b 2 e + ( b = T ω ∧ nsm convenient some in a a ω c R T ∧ µb a ∧ 2 ab ω ng dx a in c ) b µ ntrsof terms in , susual, As . (2) 2

a b the spacetime metric follows as gµν = eµeν ηab, where ηab CS action was recently generalized to 3+1d (and higher is the local orthonormal (Lorentz) metric. This invariant odd space dimensions) crystalline topological insulators, can be written, using the associated Bianchi identities, using so-called dimensionful elasticity tetrads E with di- as mension [E] = [1/L] = [M]. Topological polarization terms can also be written down22 by dimensional descent. N = dQ, Q = ea . (3) ∧ Ta The ensuing higher dimensional Chern-Simons and po- N is a locally exact 4-from independent from tr(R R) larization terms are expressed as the mixed responses ∧ E A dA and E E dA with quantized dimensionless and the dual of the scalar curvature √g in the presence ∧ ∧ 22–24 ∧ ∧ of non-zero torsion. It can be associatedR with a difference coefficients . of two topological terms, albeit in terms of an embedding Another such example is the temperature correction to 2 to five dimensions12–14. In terms of four-dimensional chi- curvature effects, with δSeff = T in the low-energy action18. This represents the analogR of the gravitational ral fermions on such a spacetime, it has been suggested R that this invariant contributes to the axial anomaly, i.e. coupling (inverse Newton constant) in the low-energy ac- the anomalous production of the chiral current: tion, where the curvature scalar is the effective scalar spacetime curvature. Since [T ]2R[ ] = [M]4, the coef- 2 R µ Λ ficient of this term is dimensionless, and can be given ∂ (ej )= N(r,t) , (4) µ 5 4π2 in terms of universal constants: it is fully determined by the number of the fermionic and bosonic species in where e = det ea and ejµ is the axial current (pseudoten- µ 5 the effective theory on flat background, and thus works sor) density and the non-universal parameter Λ has di- both in relativistic and non-relativistic systems18. The mension of relativistic momentum (mass) [Λ] = [1/L]= same universal behavior takes place with the terms de- [M] and is determined by some ultraviolet (UV) energy scribing the chiral magnetic and chiral vortical effects scale. in Weyl superfluid 3He-A, where the coefficients are Given the anomaly term (4), there has been sev- dimensionless1,25,26. Similarly, it has been observed that eral attempts to consider the Nieh-Yan anomaly in con- the coefficient of the tr(R R) graviational anomaly in densed matter systems with Weyl fermions like (1), see chiral Weyl systems affects∧ the thermal transport coef- e.g.15–17,19. However, in non-relativistic systems the rela- ficients in flat space27–33. These coefficients are funda- tivistic high-energy cut-off Λ is not a well defined parame- mental, being determined by the underlying degrees of ter and, moreover, can be anisotropic. The complete UV freedom in addition to symmetry, topology and geome- theory is, of course, non-Lorentz invariant and the lin- try. From this perspective especially, since the NY form ear, quasirelativistic Weyl regime is valid at much lower is second order in derivatives and can be computed from scales. Moreover, the anomalous hydrodynamics of su- linear response, our findings are very interesting and war- perfluid 3He at zero temperature suggests that the chiral rant further research. anomaly is either completely exhausted by the emergent Our goal in the present paper is to separate different T 2 axial gauge field corresponding to the shift of the node contributions in the hydrodynamic free-energy of Weyl or, conversely, the gravitational NY anomaly term aris- superfluids in order to identify the terms responsible for ing from the tetrad and spin-connection for local Lorentz different relativistic phenomena, including a term from invariance along the uniaxial symmetry direction. It was the thermal Nieh-Yan anomaly, as well as gravitational recently shown in Ref. 15 that the low-energy theory terms of the form T 2( + µ2), where µ is the chiral satisfies the symmetries and conservations laws related chemical potential µ TR. to an emergent quasirelativistic spacetime with torsion R ≪ and Λ is determined from the UV-scale where the linear Weyl approximation breaks down, as dictated by the underlying p-wave BCS Fermi superfluid. Here we will con- III. TEMPERATURE CORRECTION TO THE sider the temperature corrections to the anomaly and the NIEH-YAN TERM hydrodynamic free-energy in the chiral Weyl superfluid, with the expectation that Λ = T , kB = 1, with some di- The relativistic zero-temperature anomaly term in the 2 a a b mensionless prefactor γ, in units where the tetrads have axial current production Λ ( a + e e Rab) canonical dimensions of velocity, or the relevant “speed is still not confirmed in general,T ∧ see T however∧ Ref.∧ 15. of light” is set unity. On one hand the UV cut-off parameter Λ is not well- In terms of effective low-energy actions, the fully rel- defined in relativistic field theory with fundamental chi- ativistic analogs work unambiguously only for terms ral fermions. On the other hand, such a cut-off is in the effective action with dimensionless coefficients. not in general available in non-relativistic matter with Perhaps the most well-known example being the 2+1- quasi-relativistic low-energy chiral fermions and can be dimensional topological Chern-Simons (CS) terms de- anisotropic15 or even zero. However, the term of the form 2 a a b scribing the Quantum Hall effect. Gravitational Chern- γT ( a +e e Rab) has the proper dimensionality Simons terms can be similarly quantized in terms of chi- [M]4,T and∧T its prefactor∧ ∧ γ could be a universal constant ral central charge which has relation to thermal trans- in canonical units, being expressed via some invariant port and the boundary conformal field theory20,21. The related to the degrees of freedom. 3

For concreteness, we focus on the ﬁnite temperature where n5 = n+ n− is the chiral density. The corre- Nieh-Yan anomaly in chiral p-wave Weyl superﬂuid (such sponding quasirelativistic− momentum density of the Weyl as 3He-A) with fermions, valid in the vicinity of the node, is15

µ 2 ˆ 0 ∂ (ej )= γT N(r,t) , (5) PNY-node = pF lej5 . (11) µ 5 − and check whether the dimensionless parameter γ can be Thus Eq. (8) gives the temperature correction to this universal. We now use the result obtained by Khaidukov anomalous momentum production, leading to a mass cur- and Zubkov34 and Imaki and Yamamoto35 for the finite rent due to thermal Nieh-Yan anomaly. Next we discuss temperature contribution to the chiral current. For a this in detail for the quasiparticles and the superfluid vacuum in the non-relativistic chiral Weyl p-wave super- single (complex) chiral fermion, one has for the chiral 15,36 current fluid, using a Landau level model for the currents . For relativistic Weyl fermions, see37. For considerations 2 k T kij 0 of similar temperature effects in chiral Weyl superfluid, j = ǫ . (6) 49 5 − 24 Tij see Ref. We assume that this current can be covariantly generalized to the 4-current: A. Nieh-Yan term from the hydrodynamics of chiral Weyl superfluid 2 µ T µναβ a ej5 = ǫ eνa αβ . (7) − 24 T It is known that the hydrodynamics of chiral gapless 3 Then one obtains the divergence He-A experiences momentum anomalies related to the spectral flow through the Weyl points. Let us express 2 Eqs. (10), (11) in terms of the hydrodynamic variables µ T µναβ a ∂µ(ej )= ǫ aµν . (8) and quasiparticles of the chiral superfluid. The Weyl 5 − 48 T Tαβ fermions arise from the BdG Hamiltonian close to the In the presence of curvature R(ω), this becomes the tem- nodes, perature correction to the full Nieh-Yan term in Eq. (5), 1 i ǫ( i∂) 2i ∂i, ∆ where now the non-universal cut-off Λ is substituted by HBdG( i∂)= 1 − ∗i { } (12) − ∂i, ∆ ǫ(i∂) the well defined temperature T , and the dimensionless 2i { } − 2 2 parameter γ =1/12: p −pF Here ǫ(p)= 2m∗ is the normal state dispersion minus T 2 the Fermi level µF ; ∆i = c⊥(mˆ + inˆ) is the order param- ∂ (ejµ)= N(r,t) . (9) ˆ µ 5 − 12 eter in the p + ip chiral superfluid; l = mˆ nˆ is the unit vector in the direction of the orbital angular× momentum Note it is possible that the local relativistic (Tolman) of Cooper pairs; anisotropic node velocities are ck = vF temperature T = T / e0 enters the local anomaly, while and c = ∆ /p , constituting effective speeds of light in 0 | t | ⊥ 0 F the constant T0 is the global equilibrium temperature of Weyl equation along ˆl and in transverse directions, re- 18 the condensed matter system . In (1) we have simply spectively. In the weak coupling BCS theory c⊥ ck, e0 = 1. 3 −3 ≪ t − in He-A their ratio is of order 10 . In a chiral superconductor we in addition perform the minimal coupling ǫ( i∂i) ǫ(iDi)+ eA0, where Di = ∂i eAi is the IV. FROM RELATIVISTIC PHYSICS TO − → − gauge covariant derivative and Aµ is the electromagnetic CHIRAL WEYL SUPERFLUID eA potential. Equivalently this amounts to vs vs m in the free-energy. In what follows we ignore→ the super-− In the presence of a finite Weyl node pW = 0, the chi- fluid velocity (giving rise to spin-connection in addition 6 ral anomaly for the chiral current leads to the anomalous to tetrads) in the anomaly considerations. Then the only production of the linear momentum1,15. Even though hydrodynamic variable appearing in the torsion is the chiral current is not well-defined at high-energy, the spec- unit vector of the orbital momentum ˆl15. Our results tral flow of chiral quasiparticles is accompanied by the can afterwards then generalized to include the superfluid spectral flow of the linear momentum pW of the Weyl velocity appearing in the free-energy. point. In 3He-A there are two (spin-degenerate) Weyl The Weyl nodes are at pW ± = pFˆl. Expanding the points with opposite chirality and opposite momenta, Hamiltonian as H σaei (ˆp ±p ) , the vierbein ei ˆ ˆ BdG a W i a pW ± = pF l, where l is the unit vector of the orbital in the vicinity of the Weyl≃ point− take the form: momentum± of the superfluid. In particular, the anomalous production of the linear momentum density n± at 1 0 nodes of the two opposite Weyl points sum-up leading to µ t i 0 c⊥mˆ ea = ea, ea =   , a =0, 1, 2, 3. (13) { } 0 c⊥nˆ ˆ P˙ anom = pF l (∂tn5), (10)  0 c ˆl  −  k    4

i a i For the inverse vierbein eaej = δj we have whereas the finite temperature contribution to the quasiparticle momentum is 10 0 0 2 a a a qp pF T eµ = et , ei = 1 1 1 ˆ . (14) j (T )= ˆl(ˆl ˆl). (18) { } 0 c mˆ c nˆ c l anomk 2 ⊥ ⊥ k ! − 6c⊥ ·∇× In order to compute the thermal contribution from the and arises due to thermal normal component close to the quasiparticles in the presence of torsion, we assume that nodes. These contributions to the anomalous vacuum the tetrad gives rise to a constant torsion via mˆ = xˆ, current from the perspective of the superfluid are ana- 15,36 nˆ = yˆ TBxˆz, where TB 1 is a perturbation . In lyzed in the next subsection. this case,− the 3d spectrum organizes≪ into one-dimensional The relativistic anomaly results in Sec. III and the states on 2d Landau levels (LLs)1,16,36–38, where the rel- Landau level argument for the chiral superfluid suggest evant spatial torsion in addition the following temperature correction to the anomalous momentum in 3He-A, in the vicinity of the 15 1 ˆl node, from the thermal Nieh-Yan anomaly : ǫijk 3 ei = c ˆl ( ) T , (15) jk 3 k B 2 2 T · ∇× ck ≡ p T P (T )= p ˆl n (T )= F ˆl(ˆl ( ˆl)). (19) NY-node − F 5 12c2 · ∇× ˆ ˆ z ⊥ where TB = l l = xy, is playing the role of effective magnetic·∇× field. OnlyT the gapless lowest LLs are With the Landau level approximation, the quasiparticle density n (T ) between the two Weyl nodes is computed relevant, while the gapped levels, with gap c⊥√pF TB 5 are particle-hole symmetric and cancel out36,39∼ . We ap- from a similar integral as Eq. (16), 2 2 2 2 p −pF pz −pF ∞ proximate ǫ(p) = 2m 2m = ǫ(pz), strictly valid when p mc which≃ coincides with the linear Weyl n5(T )= 2 NLL(pz)dpz nF (ǫpz µF ) Θ(µ ǫpz ) ⊥ ≪ ⊥ − 0 − − − regime ǫ = v (p p ). However, we can fix the Z ∞ pz F z F T anisotropic dispersion± to− ǫ(p ) for all momenta to obtain B z = 2 dx xnF (x) a convenient UV-complete model, for which the following −2π 0 2 Z analysis for the total vacuum current is valid, as long as T ˆ ˆ = 2 (l l). (20) the correct cutoff for the linear Weyl regime is maintained −12c⊥ ·∇× for 3He-A, see Appenix A 2. The dispersion of the low- Similar to the total current, this is the contribution est LL becomes En=0 = sgn(pzTB)ǫ(pz), while density − from the linear spectrum close to the nodes with ∆pk of states per lowest LL per three-momentum pz becomes c⊥ T ∼ |pzTB | ( c ) c due to thermal fluctuations, which renormalizes nLL(pz) = 4π2 . The lowest LL is particle-like (hole- k ⊥ ˆ the velocity coefficient and requires the cutoff of the Weyl like) for pW = pF l. For more on this non-relativistic 3 model with torsion,∓ see15,36,38,39 regime relevant for the full dispersion, such as He-A. i The thermal fluctuations for the superfluid are discussed The total (chiral) quasiparticle current along ˆl = e3 vF ≈ in Appendix A 1. ˆz becomes What have been calculated above are actually the tem- ∞ perature dependence of the anomalous torsional conduc- qp ˆ P l = 2 NLL(pz)dpz pznF (ǫpz µF ) tivities of the superfluid, i.e. then chiral momentum- · − 0 − Z density and number densities in response to TB, see p3 p T 2 30 = F + F ˆl ˆl (16) e.g. . Note in particular that the leading torsional − 6π2 6c2 ·∇× contribution from the NY anomaly in Eqs. (16), (20) ⊥ jvac jqp from thermal fluctuations is suppressed by the the = anomk + anomk(T ). −2 factor c⊥ , in contrast to the antisymmetric torsion 1 ijk 3 3 TB βx −1 2 where n (x) = (e + 1) is the quasiparticle distribu- S = 2 ǫ ei jk = c from (14), as was found in F T k tion function and a factor of two comes from the spin- Ref.15 without the Landau level approximation to the degeneracy. In order to compute the integral, we have anomaly. It is also intriguing to see that the non- used T ∆0 µF as well as the linear expansion relativistic thermal integrals in Eqs. (16),(20) with ex- ≪ ≪ ǫ(∆pz) vF (∆pzc⊥/vF ) T with anisotropic momen- plicit non-relativistic vacuum regulatization due to the ≈ ≈ T T tum scaling, where close to the nodes ∆pz pF = ∆0 c⊥ filled quasiparticle states coincide with similar expression 2 ∼ for both T,ǫ mc⊥, the cutoff of the linear Weyl regime for relativistic Weyl fermions with positive and negative in 3He-A. See≪ Appendix A 2 and15,38 for more details. branches26,28–30,34,35, see Appendix A 3. vac From the above we conclude that janom is the anomalous superfluid current from filled quasiparticle states, see Eq. (21) below, B. Anomalous vacuum current

p3 C jvac = F ˆl(ˆl ˆl)= 0ˆl(ˆl ˆl), (17) The hydrodynamic anomaly in momentum conserva- anomk −6π2 ·∇× − 2 ·∇× tion arises between the quasiparticles and vacuum mass 5 current P = j of the superfluid, which at T = 0 has the which is the Fermi-liquid corrected result Eq. (25) and following general form1,40,41: corresponds to the result (18) from the reduction of superfluid density when the Fermi-liquid corrections are ig- 1 C0 j = ρvs + (ρˆl) ˆl(ˆl ( ˆl)) , (21) nored. 4m∇× − 2 · ∇× We summarize these findings as follows. While Eq. where the last term is anomalous with the parameter (23) with temperature corrections Eq. (24) looks similar C0 from the combined orbital-gauge symmetry of the to Eq. (19), Panom(T ) represents the consistent anomaly superfluid42, that fully determines the axial anomaly vacuum momentum density from all filled states which in the system due to the Weyl quasiparticles of the depends on the non-fundamental parameters m∗ and m superfluid1 via the Fermi-liquid corrections, while we expect that the T 2-contribution to the (covariant) NY anomaly arise C (T =0)= p3 /3π2 = ρ, (22) 0 F from contributions close to the linear Weyl nodes and 2 2 should contain fundamental and universal prefactors. where we ignore corrections of the order of (∆0/EF ) = 2 (c⊥/ck) 1 at zero temperature to the density from Concerning the vacuum momentum current of the su- ≪ 2 2 perfluid, one reason is that there are several different T 2 pairing (ρ C0 ρ(∆0/EF )). Notably, this term ex- ists only on− the∼ weak coupling side of the topological contributions to the current and free energy in the chiral BEC-BCS Lifshitz transition, where the pair of the Weyl superfluid, and they correspond to phenomemena with different origins and scales. Although of similar form in points with pW ± = pFˆl appears in the quasiparticle spectrum1,42. ± terms of the low-energy Goldstone variables of the su- Here we are interested in the temperature correction to perfluid, they can expressed in relativistic form with fun- the anomaly, which may come from the Nieh-Yan term. damental parameters only when carefully keeping track The extension of the anomalous current to nonzero tem- of each individual contribution to avoid double counting. perature gives In particular, in the next section we shall see that the non-fundamental parameters m∗ and m do not enter the 1 P (T )= C (T )ˆl(ˆl ( ˆl) , (23) free energy or the final results, when experessed in terms anom −2 0 · ∇× of the correct relativistic variables valid in the quasirelativistic low-energy theory. where according to Cross40,41, the anomalous parameter C0(T ) has the following temperature dependence at low T Tc: ≪ V. RELATIVISTIC CORRECTIONS TO p m∗ C (T )= C (0) T 2 F 1+ . (24) FREE-ENERGY 0 0 − 6c2 m ⊥ 2 Here m∗ is the effective mass of quasiparticles in the nor- Let us consider the T -corrections in the free energy mal Fermi liquid, which differs from the bare mass m which are second order in derivatives containing combi- 3 ˆ ˆ ˆ of the He atom due to the Fermi liquid corrections. nations of (l vs) and (l ( l)), neglecting all higher 4 · 3 · ∇× In Eqs. (17) and (18) we have neglected the Fermi- order O(T , ∂ ) terms. These terms can be distributed liquid corrections due to interactions: An additional fac- into three groups, which have different dependence on m∗ 1 m∗ (0) and m: tor 2 ( m 1)(ρn /ρ) arises from the reduced quasiparticle momentum− flow, due to Galilean invariance. The current becomes41 F = F1 + F2 + F3 , (28) ∗ p T 2 1 m (0) F [ˆl, v ]= F (ˆl v )(ˆl ( ˆl)), (29) janom = janom (25) 1 s 2 s 1 s (0) m 12 c⊥ · · ∇× 1+ 3 F1 (ρnk /ρ) ∗ 2 2 pF m T ∗ F [ˆl, v ]= 4m(ˆl v ) (ˆl ( ˆl)) , where the Landau parameter 1 F s = m 1 and ρ(0) (T ) 2 s − 96m2 c2 · s − · ∇× 3 1 m − nk,⊥ ⊥ are the bare thermal quasiparticle densities without (30) Fermi-liquid corrections along and perpendicular to ˆl. p T 2 ˆl F ˆl ˆl 2 They are given by, see the Appendix A 1, F3[ ]= ∗ 2 ( ( )) , (31) −288m c⊥ · ∇× 2 2 4 4 (0) π T (0) 7π T ρ = ρ, ρ = ρ, (26) The relativistic form of each of these free-energy contri- nk ∆2 n⊥ 15∆4 0 0 butions arises separately as follows. where by Galilean invariance, ρδij = ρsij + ρnij at all The term F3 in Eq. (31) describes the universal tem- temperatures. From Eqs. (25), (26) we gather perature correction to the Newton gravitational coupling, which depends on the number of fermionic species18: C (T ) p3 p T 2 p T 2 m∗ 0 = F F F 1 + O(T 4), 2 6π2 − 6c2 − 12c2 m − 2 2 ⊥ ⊥ vF T ˆ ˆ 2 T (27) F3[ ,T ]= 2 (l ( l)) = √ g . (32) R −288 c⊥ · ∇× 144 − R 6

Note that being exressed in terms of the scalar curva- presence of an explicit vacuum of filled quasiparticles, al- 2 ture and metric determinant √ g = 1/vF c⊥, this term though the end results are similar to relativistic fermions becomes universal: it does not− contain the microscopic when interpreted carefully in terms of the anisotropy and ∗ parameters of the system: pF , m and m . The prefactor cutoff of the quasirelativistic Weyl regime. is fully determined by the number of fermionic species. What we calculated, via the anisotropic Landau level The term F2 in Eq. (30) is expressed in terms of the model with non-relativistic symmetries in Sec. IV, are ˆ ˆ 1 ˆ ˆ actually the temperature dependence of the anomalous combination l v = l vs 4m l ( l), proportional to the · · − · ∇× a ground state current which does not receive corrections torsional conductivities σT of the quasiparticle system, due to Galilean invariance. Here the velocity v = j/ρ i.e. the chiral momentum-density and number densities is the velocity of the “total quantum vacuum”, where (at the nodes at finite momenta pF ) in response to the 1 ijk 3 3± j j(T = 0) is the total vacuum current in Eq. (21) at spacelike torsion TB = ǫ ei jk. Namely, for example, ≡ 2 T T = 0. We conclude that the F2 contribution gives the temperature T and chemical potential µ correction to the a a 0i a σT b (T ) 0ijk b a3 P = e T = e ǫ T = δ σTB TB, (36) free energy of the gas of chiral fermionic particles in the h i h i i i 2 jk 1 limit µ T (see Eqs. (9.12) and (10.42) in Ref. ): i 0i µν | | ≪ P = T , where T is the stress-tensor, which are of second order in derivates and can be calculated in p m∗ T 2 2 F ˆl v ˆl ˆl linear response, e.g. Kubo formulas, to the back- F2[µ,T ]= 2 2 4m s ( ) = − 96m c⊥ · − · ∇× 27,30,44,45 ground tetrads . A similar momentum anomaly T 2 in anisotropic system with non-relativistic symmetries = √ gµ2 , (33) − 6 − was also considered in46. For the zero temperature case, see15,38 as well as37,47 for a relativistic model related to where the chiral chemical potential of the superfluid is Weyl semi-metals. We, however, stress that the universal determined by the Doppler shift gravitational NY anomaly and thermal physics we discussed arise in flat space from the geometric background µ = µ = p v ˆl = p ˆl v (p /4m)ˆl ( ˆl) . (34) R − L F · F · s − F · ∇× fields in the low-energy quasiparticle Hamiltonian. These tetrads arise universally in all Weyl systems (1) and cou- The final version of Eq. (33) does not contain microscopic ple to the momentum16,48 as in gravity. The connection parameters. The Eq. (33) also gives rise to the mass of of our results and the relation of the gravitational NY “photon” in 3He-A43. In principle, the chiral chemical anomaly, with the coefficient γ, to thermal transport in potential may also serve as the Λ parameter in the Nieh- Weyl system should be further elucidated27–29,31,33,49,50. Yan term, see however Appendix A 3. Detailed consideration of the temperature depen- Finally F in Eq.(29) is the term in free energy, which 1 dent anomaly terms in the hydrodynamics of the non- gives rise to the vacuum current without any factors ∗ relativistic p-wave chiral superfluid with quasirelativistic of m with contribution from the Nieh-Yan thermal m Weyl fermions demonstrates that in the hydrodynamics anomaly in Eqs. (19), (24), (27): of this liquid there are several T 2 terms, which can be δF p T 2 assigned to different emergent relativistic phenomena, P (T )= 1 = F ˆl(ˆl ( ˆl)) (35) both anomalous and non-anomalous. In particular, we anom δv 12 c2 s ⊥ · ∇× identified and discussed the term in the vacuum momen- pF n5(T )ˆl, tum corresponding to the (consistent) thermal Nieh-Yan ≃− anomaly. As expected this term originates from ther- where the last equality was derived in Eq. (20). We stress mal fluctuations close to the linear nodes with emergent however that Panom(T ) represents the total anomalous quasirelativistic torsion with anisotropy. Note that in superfluid current, whereas the right-hand side is equal terms of the superfluid and Weyl fermions, the T = 0 to the quasiparticle momentum density due to thermal anomalous vacuum contribution to the current can be as- fluctuations close to the node. signed only to a non-local action51. We showed how the various T 2 low-temperature corrections can be written aslow-energy relativistic terms with dimensionless pref- VI. CONCLUSIONS actors, which do not seem to depend on microscopic physics, but are fully determined by geometry, topol- We discussed the possibility of thermal Nieh-Yan ogy and the number of fermionic and bosonic quantum anomaly where the role of the non-universal dimension- fields. Detailed comparison of the finite temperature su- ful UV cut-off Λ is played by the temperature IR scale, perfluid hydrodynamics with Fermi-liquid corrections to and the dimensionless prefactor γ in the anomaly is uni- the anomalous quasiparticle axial current production in versal. We identified a contribution from this anomaly in the presence of arbitrary textures and superfluid velocity the known low-temperature corrections of non-relativistic remains to be identified36,52 . See however Ref. 15 for chiral p-wave Weyl superfluid (or superconductor). In the zero temperature case. this system, the anomaly results from thermal effects Note added: After the initial submission of this of the linear Weyl spectrum at finite momentum in the manuscript as a preprint, arXiv:1909.08936v1, with the 7

2 predicted T -contribution to the NY anomaly, Eqs. (9), anisotropic contributions from, with ωn = πT (2n + 1), (35), the recent preprints37,47,53 discussing related torsional anomaly phenomena at finite temperatures ap- ∞ ∂n ∞ ∞ ω2 E2 F n − kF peared. General aspects of the temperature anomaly dǫ = dǫ T 2 2 −∂Ek (ω + E F ) −∞ kF −∞ n=−∞ n k in Weyl materials were further discussed in the short Z Z X paper54. In particular, the result Eq. (9) has been con- ∞ ∆ 2 = πT − | kF | (A2) firmed in Ref. 37 by a direct calculation of the spectral 2 2 3/2 n=−∞ (ωn + ∆kF ) flow of relativistic Landau levels in the presence of a con- X | | stant torsional magnetic field T 3 16,36,39 at finite temper- µν πT ature. Here similar computations for the non-relativistic The dimensionless summation variable is xn = ∆0 (2n + Weyl superfluid in Eqs. (16) and (20) give corresponding 1), giving for T ∆0 ≪ results. While the current manuscript was being final- 55 (0) ized, also the paper appeared discussing the anomaly ∞ 1 2 2 ρnk πT u (1 u ) for relativistic fermion at finite temperature and chemical = du − − 3ρ ∆ (x2 + (1 u2))3/2 potential. 0 n=−∞ −1 n X Z − Acknowledgements: GEV thanks Mike Zubkov for dis- πT ∞ 1 = 3 x (3x2 + 1)arctan cussions. JN thanks Z.-M. Huang for correspondence 2∆ | n|− n x 0 n=−∞ n and T. Ojanen for discussions. This work has been sup- X | | ported by the European Research Council (ERC) under πT ∞ π = 4 x (3x2 1) + O(x4 ) (A3) the European Union’s Horizon 2020 research and inno- ∆ | n|− n − 2 n 0 n=−∞ vation programme (Grant Agreement No. 694248). X h i and

(0) ∞ Appendix A: Vacuum regularization and thermal ρ πT 1 du (1 u2)2 n⊥ = − − integrals 3ρ ∆ 4 (x2 + (1 u2))3/2 0 n=−∞ −1 n X Z − πT ∞ 2 x In this appendix we shortly review the low- = 3 x | n| 4∆ | n|− 1+ x2 temperature corrections to the chiral superfluid and the 0 n=−∞ n X quasirelativistic thermal integrals utilized in the main 1 text. Although the temperature and Fermi-liquid correc- +(1 3x2 ) arctan − n x tions in Sec. A 1 chiral superfluid are well-known1,40,41, n ∞ | | it is of interest to compare them to the chiral quasirela- πT 16 π = x 3 + (3x2 1) + O(x5 ) tivistic Landau level model with torsion and explicit UV 4∆ − 3 | n| n − 2 n 0 n=−∞ completion15,36,38 Sec. A 2, as well to relativistic Weyl X fermions with positive and negative branches26,28,30 Sec. We use the regularizations, ζ (s)= ∞ (n + a)s, A 3. The last system has particles and antiparticles at a n=0 finite chemical potential counted from the node and fi- ∞ ∞ P nite temperature with vacuum contributions subtracted, =1+2 =1+2ζ(0) = 0, in contrast to non-relativistic systems at finite chemical n=−∞ n=1 potential. ∞ X X 4πT πT x = ζ ( 1) = , | n| ∆ 1/2 − 6∆ n=−∞ 0 0 X∞ 1. Temperature corrections to normal and 8π2T 2 x2 = ζ ( 2)=0, (A4) superfluid density in chiral superfluid n ∆2 1/2 − n=−∞ 0 ∞X 16π3T 3 7π3T 3 We start with the corrections to the superfluid and x 3 = ζ ( 3) = | n| ∆3 1/2 − − 60∆3 normal density. The anisotropic chiral p-wave superfluid n=−∞ 0 0 density is X

These give the result in Eq. (26). Similarly ρsij (T ) = ∞ (0) ˆ ˆ ∂nF (Ek) ρδij ρnij (T ) by Galilean invariance, valid at all temper- ρ =3ρ dΩ kikj dǫ (A1) − nij − ∂E atures with or without Fermi liquid corrections. Z Z−∞ k k=kF From the superﬂuid and normal densities, the conclu- βx −1 sion is that only anisotropic momenta of the order of where nF (x) = (e + 1) with β = 1/T and Ek = T T c⊥ T ∆p⊥ pF = and ∆pk = ( ) , ck = vF con- 2 2 ∼ ∆0 c⊥ ck c⊥ ǫk + ∆k is the quasiparticle energy with normal | | 2 tribute in the vicinity of the gap nodes to the thermal k qstate dispersion ǫk = µ . We compute the ﬂuctuations ǫ T . 2m − F ∼ 8

2. Quasirelativistic chiral fermions with torsion whereas the quasiparticle current due to the thermal ﬂuc- tuations is

The quasirelativistic model for the chiral superfluid 2 ˆ qp pF T ˆ ˆ ˆ with anisotropic dispersion is, with l ˆz, janomk(T )= 2 l(l l). (A10) ≃ − 6c⊥ ·∇× 2 2 p pz ǫp = µF ǫpz µF (A5) Similarly, the thermal density of states is the vicinity 2m − → ≡ 2m − of the nodes for T ∆0 is which is valid for real 3He-A for E mc2 , c ∆ /p . ≪ ≪ ⊥ ⊥ ≡ 0 F ∞ This limit coincides with linear Weyl regime, where ǫp = n5(T )= 2 NLL(pz)dpz nF (ǫpz µF ) Θ(µ ǫpz ) c (p p ), c v , but the model (A5) allows calcu- − − − − ± k z − F k ≡ F Z0 lations of Landau levels and the filled quasiparticle states T ∞ = B dx xn (x) up till pz = 0 that contribute to the superfluid vacuum. −2π2 F In case of the anomalous current, only the momentum Z0 T 2 along ˆl ˆz contributes justifying the approximation as ˆ ˆ = 2 (l l). (A11) ≃ 2 c⊥ −12c ·∇× long as the cutoff T ∆0 and E mc = ∆0 are ⊥ ≪ ≪ ⊥ ck maintained. In particular, this requires that where the T = 0 vacuum contribution is substracted, since only contribution from the vicinity of the nodes to c⊥ T ǫp(pF + ∆p) ǫpz (pF + ∆p ) c T the chiral density is well-defined. The relevant thermal → k ∼ k c c ∼ k ⊥ integral is similar to Eq. (A2) (A6) ∞ 1 ∞ ∂n (x) by the equipartition principle. The momenta (and en- dx xn (x)= x2 F (A12) F 4 − ∂x ergies) should be anisotropically scaled in the simplified Z0 Z−∞ c⊥ 3 2 2 model, or equivalently ǫpz ( c )ǫpz for He-A with tor- 2 2 π T → k =4π T ζ1/2( 1) = (A13) sion. The density of state per pz per Landau level is − 6 p T N (p )= | z B| (A7) LL z 4π2 3. Relativistic chiral fermions ˆ ˆ where TB = l l. For more on this anisotropic chiral Now we compare these results with relativistic fermion with·∇× torsion, see Ref. 15 and the forthcoming 38 Weyl fermions with both positive and negative energy paper . branches, i.e. particles and antiparticles with the node In case of the Landau level calculation, the quasirel- at pz = 0. At chemical potential µ at the node, the ativistic integrals are given by the lowest Landau level corresponding momentum along ˆl is (LLL) contribution, which is particle- and hole-like for p < 0 and p > 0, respectively. In the chiral superfluid, ∞ z z ˆ the quasiparticles are spin-degerate Majorana-Weyl ex- PW l = NLL(pz)dpzpz nF (ǫpz µ) Θ( pz) · −∞ − − − citations. It follows that the quasiparticle momentum Z ∞ ˆ TB along is the l = zˆ direction is = dp p2 n (ǫ µ) n (ǫ + µ) −4π2 z z F pz − − F pz ∞ Z0 ˆ µ ∞ Pqp l = 2 NLL(pz)dpzpznF (ǫpz µF ) TB 2 x · − 0 − = 2 dpzpz +4µ dx xnF ( ) Z −4π 0 0 T 2T ∞ Z Z = B dp p2n (ǫ µ ) µ3 µT 2 2 z z F pz F ˆl ˆl − 4π 0 − = 2 + , (A14) Z pF ∞ − 12π 12 ·∇× TB 2 4pF x = 2 dpzpz + 2 dx xnF ( ) −2π c T per single node, in units where ǫpz = pz. The thermal Z0 ⊥ Z0 p3 p T 2 chiral density, corresponding to the left-handed node, = F + F ˆl ˆl − 6π2 6c2 ·∇× ∞ ⊥ vac qp n (T )= N (p )dp n (ǫ z µ) Θ( p ) = j + j (T ). (A8) W − − LL z z F p − − − z anomk anomk −∞ Z ∞ TB where we have used T ∆ µF and the linearization = dp p n (ǫ µ)+ n (ǫ + µ) ≪ ≪ −4π2 z z F pz − F pz of the spectrum close to µF with the anisotropic scal- Z0 ing (A6) which applies for (weak-coupling) 3He-A. The µ2 T 2 = ( )ˆl ˆl. (A15) vacuum contribution is −8π2 − 24 ·∇× ∞ TB jvac 2 In order to compute the integral, we have used nF (x)= anomk = 2 dpzpzΘ(pF pz) (A9) −2π 0 − 1 n ( x) as well as the linear expansion close to the Z − F − 9 nodes for T ∆0. Intriguinly, the structure of this cles and that the contribution of filled vacuum states is momentum-density≪ and density contribution are exactly substracted. Lastly, we note the dissimilarity of these rel- of the same form26,28,30,34,35. The difference is, of course, ativistic finite chiral chemical potential and temperature that µ is counted from the node and has different sign for terms to Eqs. (32), (33), and (34). positive energy solutions corresponding to the antiparti-

1 G.E. Volovik, The Universe in a Helium Droplet, Claren- 21 N. Read, D. Green, Paired states of fermions in two dimen- don Press, Oxford (2003). sions with breaking of parity and time reversal symmetries 2 P. Horava, Stability of Fermi Surfaces and K Theory, Phys. and the fractional quantum Hall effect, Phys. Rev. B 61, Rev. Lett. 95, 016405 (2005) 10267 (2000). 3 C. Herring, Accidental degeneracy in the energy bands of 22 Xue-Yang Song, Yin-Chen He, Ashvin Vishwanath, Chong crystals, Phys. Rev. 52, 365–373 (1937). Wang, Electric polarization as a nonquantized topological 4 A.A. Abrikosov and S.D. Beneslavskii, Possible existence response and boundary Luttinger theorem, arXiv preprint of substances intermediate between metals and dielectrics, arXiv:1909.08637 JETP 32, 699–798 (1971). 23 J. Nissinen and G.E. Volovik, Tetrads in solids: from 5 H.B. Nielsen and M. Ninomiya, The Adler-Bell-Jackiw elasticity theory to topological quantum Hall systems and anomaly and Weyl fermions In a crystal, Phys. Lett. B130, Weyl fermions, ZhETF 154, 1051 (2018) [JETP 127, 948 389 (1983). (2018)]. 6 G.E. Volovik, Chiral anomaly and the law of conservation 24 J. Nissinen and G.E. Volovik, Elasticity tetrads, mixed of momentum in 3He-A, Pis’ma ZhETF 43, 428 (1986); axial-gravitational anomalies, and (3+1)-d quantum Hall JETP Lett. 43, 551 (1986). effect, Physical Review Research 1, 023007 (2019), 7 A.A. Zyuzin and A.A. Burkov, Topological response in arXiv:1812.03175. Weyl semimetals and the chiral anomaly, Phys. Rev. B 25 G.E. Volovik and A. Vilenkin, Macroscopic parity violat- 86, 115133 (2012) ing effects and 3He-A, Phys. Rev. D 62, 025014 (2000); 8 H.T. Nieh and M.L. Yan, An identity In Riemann-Cartan hep-ph/9905460. geometry, J. Math. Phys. 23, 373 (1982). 26 Gök¸ce Ba¸sar, Dmitri E. Kharzeev, and Ho-Ung Yee, Trian- 9 H.T. Nieh and M.L. Yan, Quantized Dirac field in curved gle anomaly in Weyl semimetals Phys. Rev. B 89, 035142 Riemann-Cartan background. I. Symmetry properties, (2014). Green’s function, Ann. Phys. 138, 237 (1982). 27 K. Landsteiner, E. Megas, and F. Pena-Benitez, Gravita- 10 H.T. Nieh, A torsional topological invariant, Int. J. Mod. tional anomaly and transport phenomena, Phys. Rev. Lett. Phys. A 22, 5237 (2007). 107, 021601 (2011). 11 S. Yajima, Evaluation of the heat kernel in Riemann - Car- 28 R Loganayagam, P Surówka, Anomaly/transport in an tan space, Class. Quantum Grav. 13, 2423 (1996). ideal Weyl gas Journal of High Energy Physics 2012 (4), 12 O. Chandia and J. Zanelli, Topological invariants, instan- 97. tons, and the chiral anomaly on spaces with torsion, Phys. 29 Kristan Jensen, R. Loganayagam and Amos Yarom, Ther- Rev. D 55, 7580 (1997) modynamics, gravitational anomalies and cones, J. High 13 O. Chandia and J. Zanelli, Torsional topological invariants Energ. Phys. 2013, 88 (2013) (and their relevance for real life), arXiv:hep-th/9708139. 30 K. Landsteiner, Anomalous transport of Weyl fermions in 14 O. Chandia and J. Zanelli, Supersymmetric particle in a Weyl semimetals, Phys. Rev. B 89, 075124 (2014). spacetime with torsion and the index theorem, Phys. Rev. 31 A. Lucas, R.A. Davison, and S. Sachdev, Hydrodynamic D 58, 045014 (1998). theory of thermoelectric transport and negative magne- 15 J. Nissinen, Emergent spacetime and gravitational Nieh- toresistance in Weyl semimetals, PNAS 113, 9463 (2016). Yan anomaly in chiral p + ip Weyl superfluids and su- 32 M. Stone, Y. Kim, Mixed anomalies: Chiral vortical effect perconductors, Phys. Rev. Lett. 124, 117002 (2020). and the Sommerfeld expansion, Phys. Rev. D 98, 025012 arXiv:1909.05846. (2018). 16 O. Parrikar, T.L. Hughes, R.G. Leigh, Torsion, parity-odd 33 J. Gooth, A.C. Niemann, T. Meng, A. G. Grushin, K. response and anomalies in topological states, Phys. Rev. Landsteiner, B. Gotsmann, F. Menges, M. Schmidt, C. D 90, 105004 (2014). Shekhar, V. Süss, R. Hühne, B. Rellinghaus, C. Felser, 17 Y. Ferreiros, Y. Kedem, E.J. Bergholtz, and J.H. Bardar- Binghai Yan, K. Nielsch, Experimental signatures of the son, Mixed axial-torsional anomaly in Weyl semimetals, mixed axial-gravitational anomaly in the Weyl semimetal Phys. Rev. Lett. 122, 056601 (2019) NbP, Nature 547, 324 (2017) 18 G.E. Volovik and A.I. Zelnikov, Universal temperature 34 Z.V. Khaidukov and M.A. Zubkov, Chiral torsional effect, corrections to the free energy for the gravitational field, Pisma v ZhETF, 108, 702–703 (2018) [JETP Lett. 108, Pisma ZhETF 78, 1271–1276 (2003), JETP Lett.78, 751– 670 (2018)], arXiv:1812.00970. 756 (2003); gr-qc/0309066. 35 S. Imaki and A. Yamamoto, Lattice field theory with tor- 19 Liang Sun and Shaolong Wan, Chiral viscoelastic response sion, Phys. Rev. D 100, 054509 (2019). in Weyl semimetals, EPL 108, 37007 (2014) 36 G.E. Volovik, Normal Fermi liquid in a superfluid 3He-A 20 G. E. Volovik, The gravitational-topological Chern-Simons at T = 0 and the anomalous current, JETP Lett. 42, 363 term in a film of superfluid 3He-A, Pis’ma Zh. Eksp. Teor. (1985) Fiz. 51, 111 (1990) [JETP Lett. 51, 125 (1990)]. 37 Ze-Min Huang, Bo Han and M. Stone, The Nieh-Yan 10

anomaly: torsional Landau levels, central charge and Lifschitz theories, JHEP 2020, 190 (2020). anomalous thermal Hall effect, Phys. Rev. B 101, 125201 47 Ze-Min Huang, Bo Han and M. Stone, Hamiltonian ap- (2020). arXiv:1911.00174 (2019). proach to the torsional anomalies and its dimensional lad- 38 S. Laurila and J. Nissinen, Torsional Landau levels and der, Phys. Rev. B 101, 165201 (2020), arXiv:1912.06051. geometric anomalies in condensed matter Weyl systems, 48 H. Shapourian, T.L. Hughes, and S. Ryu, Viscoelastic re- to appear, (2020). sponse of topological tight-binding models in two and three 39 A.V. Balatskii, G.E. Volovik and V.A. Konyshev, On the dimensions, Phys. Rev. B 92, 165131 (2015). chiral anomaly in superfluid 3He-A, Zh.Eksp.Teor.Fiz. 90, 49 T. Kobayashi, T. Matsushita, T. Mizushima, A. Tsuruta, 2038 (1986), [Sov.Phys. JETP 63, 1194 (1986)]. and S. Fujimoto, Negative thermal magnetoresistivity as 40 D. Vollhardt and P. Wlfle, The Superfluid Phases of He- a signature of a chiral anomaly in Weyl superconductors, lium 3 (Taylor & Francis, London, 1990). Phys. Rev. Lett. 121, 207002 (2018). 41 M.C. Cross, A generalized Ginzburg-Landau approach to 50 J. M. Luttinger, Theory of thermal transport coefficients, the superfluidity of helium 3, J. Low Temp. Phys. 21, 525– Phys. Rev. 135, A1505 (1964). 534 (1975). 51 G.E. Volovik, Wess-Zumino action for the orbital dynamics 42 G. E. Volovik and V. P. Mineev, Orbital angular momen- of 3He-A, Pis’ma ZhETF 44, 144-148 (1986) [JETP Lett. tum and orbital dynamics: 3He-A and the Bose liquid, Zh. 44 185 (1986)] Eksp. Teor. Fiz. 81, 989 (1981) [Sov. Phys. JETP 54, 524 52 R. Combescot and T. Dombre, Twisting in superfluid 3He- (1981)], A and consequences for hydrodynamics at T = 0, Phys. 43 G.E. Volovik, Axial anomaly in 3He-A: Simulation of Rev. B 33, 79 (1986). baryogenesis and generation of primordial magnetic field in 53 Long Liang and Teemu Ojanen, Topological magnetotor- Manchester and Helsinki,” Proceedings of Symposium on sional effect in Weyl semimetals, Phys. Rev. Research 2, QUANTUM PHENOMENA at LOW TEMPERATURES, 022016(R) (2020), arXiv:1912.07732. Lammi (Finland), 7-11 January, 1998, Physica B 255, 86- 54 J. Nissinen and G.E. Volovik, On thermal Nieh-Yan 107 (1998); cond-mat/9802091. anomaly in topological Weyl materials, Pis’ma ZhETF 44 Barry Bradlyn and N. Read, Low-energy effective theory 110, 797–798 (2019); JETP Lett. 110, 1–4 (2019), in the bulk for transport in a topological phase, Phys. Rev. arXiv:1911.03382. B 91, 125303 (2015). 55 S. Imaki and Zebin Qiu, Chiral torsional effect with fi- 45 A. Gromov and A.G. Abanov, Thermal Hall effect and ge- nite temperature, density and curvature, arXiv:2004.11899 ometry with torsion, Phys. Rev. Lett. 114, 016802 (2015) (2020). 46 C. Copetti, Torsion and anomalies in the warped limit of