Skip to main content
  • Submit
  • About
    • Editorial Board
    • PNAS Staff
    • FAQ
    • Accessibility Statement
    • Rights and Permissions
    • Site Map
  • Contact
  • Journal Club
  • Subscribe
    • Subscription Rates
    • Subscriptions FAQ
    • Open Access
    • Recommend PNAS to Your Librarian
  • Log in
  • My Cart

Main menu

  • Home
  • Articles
    • Current
    • Special Feature Articles - Most Recent
    • Special Features
    • Colloquia
    • Collected Articles
    • PNAS Classics
    • List of Issues
  • Front Matter
  • News
    • For the Press
    • This Week In PNAS
    • PNAS in the News
  • Podcasts
  • Authors
    • Information for Authors
    • Editorial and Journal Policies
    • Submission Procedures
    • Fees and Licenses
  • Submit
  • About
    • Editorial Board
    • PNAS Staff
    • FAQ
    • Accessibility Statement
    • Rights and Permissions
    • Site Map
  • Contact
  • Journal Club
  • Subscribe
    • Subscription Rates
    • Subscriptions FAQ
    • Open Access
    • Recommend PNAS to Your Librarian

User menu

  • Log in
  • My Cart

Search

  • Advanced search
Home
Home

Advanced Search

  • Home
  • Articles
    • Current
    • Special Feature Articles - Most Recent
    • Special Features
    • Colloquia
    • Collected Articles
    • PNAS Classics
    • List of Issues
  • Front Matter
  • News
    • For the Press
    • This Week In PNAS
    • PNAS in the News
  • Podcasts
  • Authors
    • Information for Authors
    • Editorial and Journal Policies
    • Submission Procedures
    • Fees and Licenses

New Research In

Physical Sciences

Featured Portals

  • Physics
  • Chemistry
  • Sustainability Science

Articles by Topic

  • Applied Mathematics
  • Applied Physical Sciences
  • Astronomy
  • Computer Sciences
  • Earth, Atmospheric, and Planetary Sciences
  • Engineering
  • Environmental Sciences
  • Mathematics
  • Statistics

Social Sciences

Featured Portals

  • Anthropology
  • Sustainability Science

Articles by Topic

  • Economic Sciences
  • Environmental Sciences
  • Political Sciences
  • Psychological and Cognitive Sciences
  • Social Sciences

Biological Sciences

Featured Portals

  • Sustainability Science

Articles by Topic

  • Agricultural Sciences
  • Anthropology
  • Applied Biological Sciences
  • Biochemistry
  • Biophysics and Computational Biology
  • Cell Biology
  • Developmental Biology
  • Ecology
  • Environmental Sciences
  • Evolution
  • Genetics
  • Immunology and Inflammation
  • Medical Sciences
  • Microbiology
  • Neuroscience
  • Pharmacology
  • Physiology
  • Plant Biology
  • Population Biology
  • Psychological and Cognitive Sciences
  • Sustainability Science
  • Systems Biology
Research Article

Hydrodynamic theory of thermoelectric transport and negative magnetoresistance in Weyl semimetals

Andrew Lucas, Richard A. Davison, and View ORCID ProfileSubir Sachdev
PNAS August 23, 2016 113 (34) 9463-9468; first published August 10, 2016; https://doi.org/10.1073/pnas.1608881113
Andrew Lucas
aDepartment of Physics, Harvard University, Cambridge, MA 02138;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • For correspondence: sachdev@g.harvard.edu lucas@fas.harvard.edu rdavison@physics.harvard.edu
Richard A. Davison
aDepartment of Physics, Harvard University, Cambridge, MA 02138;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • For correspondence: sachdev@g.harvard.edu lucas@fas.harvard.edu rdavison@physics.harvard.edu
Subir Sachdev
aDepartment of Physics, Harvard University, Cambridge, MA 02138;
bPerimeter Institute for Theoretical Physics, Waterloo, Ontario N2L 2Y5, Canada
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Subir Sachdev
  • For correspondence: sachdev@g.harvard.edu lucas@fas.harvard.edu rdavison@physics.harvard.edu
  1. Contributed by Subir Sachdev, June 30, 2016 (sent for review June 3, 2016; reviewed by Anton Burkov and Kristan Jensen)

  • Article
  • Figures & SI
  • Info & Metrics
  • PDF
Loading

Significance

Weyl semimetals are exotic materials with negative electrical magnetoresistance: when an electric and magnetic field are applied in parallel, the induced electrical current increases upon increasing magnetic field strength. This is due to an emergent axial quantum anomaly in Weyl semimetals. We present a universal description of thermoelectric transport in weakly disordered Weyl semimetals where electron–electron interactions are faster than electron–impurity scattering. We predict negative thermal magnetoresistance: upon applying a parallel temperature gradient and magnetic field, the induced heat current increases with increasing magnetic field strength. This is caused by a distinct emergent quantum anomaly—the axial–gravitational anomaly. Measuring this effect may be the most practical route to experimentally observing this anomaly in any branch of physics.

Abstract

We present a theory of thermoelectric transport in weakly disordered Weyl semimetals where the electron–electron scattering time is faster than the electron–impurity scattering time. Our hydrodynamic theory consists of relativistic fluids at each Weyl node, coupled together by perturbatively small intervalley scattering, and long-range Coulomb interactions. The conductivity matrix of our theory is Onsager reciprocal and positive semidefinite. In addition to the usual axial anomaly, we account for the effects of a distinct, axial–gravitational anomaly expected to be present in Weyl semimetals. Negative thermal magnetoresistance is a sharp, experimentally accessible signature of this axial–gravitational anomaly, even beyond the hydrodynamic limit.

  • Weyl metals
  • hydrodynamics
  • thermoelectric effects
  • anomalies

The recent theoretical predictions (1⇓–3) and experimental discoveries (4⇓–6) of Weyl semimetals open up an exciting new solid state playground for exploring the physics of anomalous quantum field theories. These anomalies can lead to striking signatures in simple transport measurements. Upon applying a magnetic field B=Bz^ and measuring the electrical conductivity σzz parallel to B, one predicts σzz has a contribution that grows as B2 (7⇓–9). This longitudinal negative magnetoresistance is a direct signature of the anomaly associated with the Weyl points in momentum space. Similar results have also been predicted for thermal and thermoelectric transport (10, 11). Negative magnetoresistance in σ, with the predicted B2 dependence, has been observed experimentally in many different materials (12⇓⇓⇓⇓⇓–18).

So far, the theories of this negative magnetoresistance assume two facts about the dynamics of the quasiparticles of the Weyl semimetal. First, it is assumed that the quasiparticles are long lived, and that a kinetic description of their dynamics is valid. Second, it is assumed that the dominant scattering mechanism is between quasiparticles and impurities or phonons. In most simple crystals—including Weyl semimetals—it is likely that this description is reasonable.

However, there are exotic metals in which the quasiparticle–quasiparticle scattering time is much smaller than the quasiparticle–impurity/phonon scattering time. In such a finite temperature metal, the complicated quantum dynamics of quasiparticles reduces to classical hydrodynamics of long-lived quantities—charge, energy, and momentum—on long time and length scales. Most theoretical (19⇓⇓⇓⇓⇓–25) and experimental (26⇓–28) work on such electron fluids studies the dynamics of (weakly interacting) Fermi liquids in ultrapure crystals. As expected, the physics of a hydrodynamic electron fluid is qualitatively different from the kinetic regime where quasiparticle–impurity/phonon scattering dominates, and there are qualitatively distinct signatures to look for in experiments.

Experimental evidence for a strongly interacting quasirelativistic plasma of electrons and holes has recently emerged in graphene (29, 30). The relativistic hydrodynamic theories necessary to understand this plasma are different from ordinary Fermi liquid theory (31), and lead to qualitatively different transport phenomena (32, 33). The hydrodynamics necessary to describe an electron fluid in a Weyl material, when the Fermi energy is close to a Weyl node, is similar to the hydrodynamics of the graphene plasma, though with additional effects related to anomalies (34, 35). Such a quasirelativistic regime is where negative magnetoresistance is most pronounced (9), and also where interaction effects can be strongest, due to the lack of a large Fermi surface to provide effective screening.

In this paper, we develop a minimal hydrodynamic model for direct current (dc) thermoelectric transport in a disordered, interacting Weyl semimetal, where the Fermi energy is close to the Weyl nodes. The first hydrodynamic approach to transport in a Weyl semimetal may be found in ref. 36 (see also refs. 37, 38). In contrast to these, our approach ensures that the conductivity matrix is positive semidefinite and Onsager reciprocal. We apply an infinitesimal electric field Ei and temperature gradient ∂iT to a Weyl semimetal, and compute the total charge current Ji and heat current Qi using hydrodynamics. We then read off the thermoelectric conductivity matrix defined by(JiQi)=(σijαijTα¯ijκ¯ij)(Ej−∂jT).[1]In the limit where disorder, magnetic field, and intervalley scattering are perturbatively weak, we show that all conductivities may be written as the sum of a Drude conductivity for each valley fluid, and a correction due to intervalley scattering: e.g., σij=σijDrude + σijanom. We present a general formula for the coefficient of B2 in σzzanom: the quantitative dependence of this coefficient on temperature and electron density can be different from quasiparticle-based methods.

Although the qualitative form of our results (e.g., σijanom∼BiBj) is very similar to that found using kinetic theory approaches (8⇓⇓–11), we strongly emphasize that the physical interpretations are often quite different. For example, the emergence of Drude conductivities in our model is not due to the existence of long-lived quasiparticles, but due to the fact that momentum relaxation is a perturbatively slow process (31, 39). Furthermore, distinct anomalies are responsible for the negative magnetoresistance in electrical vs. thermal transport. This remains true even beyond our strict hydrodynamic limit.

In this paper, we work in units where ℏ = kB = e = 1. We will also generally set the Fermi velocity vF = 1. In our relativistic formalism, the effective speed of light is set by vF.

Weyl Hydrodynamics

We begin by developing our hydrodynamic treatment of the electron fluid, assuming the chemical potential lies close to the charge neutrality point for every node. For simplicity, we assume that the Weyl nodes are locally isotropic to reduce the number of effective parameters. It is likely straightforward, although tedious, to generalize and study anisotropic systems.

We will firstly review the hydrodynamic theory of a chiral fluid with an anomalous axial U(1) symmetry, derived in refs. 34, 35. Neglecting intervalley scattering, this theory describes the dynamics near one Weyl node. The equations of relativistic chiral hydrodynamics are the conservation laws for charge, energy, and momentum, modified by the external electromagnetic fields, which we denote with Fμν. On a curved space with Riemann tensor Rαβδγ, they read∇μTμν=FνμJμ−G16π2∇μ[ερσαβFρσRαβνμ] ,[2a]∇μJμ=−C8εμνρσFμνFρσ−G32π2εμνρσRβμναRαρσβ,[2b]where C is a coefficient related to the standard axial anomaly and G is a coefficient related to an axial–gravitational anomaly (40). For a Weyl fermionC=k4π2, G=k24,[3]with k∈ℤ the Berry flux associated with the Weyl node (41). Jμ and the energy–momentum tensor Tμν are related to the hydrodynamic variables of chemical potential μ, temperature T, and velocity uμ in a tightly constrained way (34, 35), which we review in the SI Appendix. We will take the background electromagnetic field to beF=Bdx∧dy+∂iμ0dxi∧dt,[4]with B as a constant. Constant B is required by Maxwell’s equations for the external electromagnetic field in equilibrium, at leading order.

A single chiral fluid cannot exist in a Weyl material. Instead, enough Weyl nodes must exist so that the “net” C for the material vanishes. This follows mathematically from the fact that the Brillouin zone of a crystal is necessarily a compact manifold and so the sum of the Berry fluxes associated with each node must vanish—this is the content of the Nielsen–Ninomiya theorem (7). Hence, we must consider the response of multiple chiral fluids when developing our theory of transport.

One might hope that so long as each chiral fluid has a well-behaved response, then the net conductivities are simply additive. This is not so: the transport problem is ill-posed for a single chiral fluid, once we apply a background magnetic field. To see this, suppose that we apply an electric field such that E⋅B≠0. Then, the total charge in the sample obeysdQtotdt=∫d3x ∂μJμ=CE⋅BV3,[5]with V3 the spatial volume of the metal. Even at the linear response level, we see that there is a necessary O(E) time dependence to any solution to the hydrodynamic equations (with spatial directions periodically identified). If there is no static solution to the equations of motion, then any dc conductivity is an ill-posed quantity to compute. There is also energy production in a uniform temperature gradient, proportional to G∇T⋅B, even when C=0 (SI Appendix).

The physically relevant solution to this issue is that multiple Weyl nodes exist in a real material, and this means that we must consider the coupled response of multiple chiral fluids. Rare intervalley processes mediated by phonons and/or impurities couple these chiral fluids together (8) and make the transport problem far richer for Weyl fluids than for simpler quantum critical fluids, including the Dirac fluid (32).

We label each valley fluid quantity with the labels ab…. For example, uaμ is the velocity of valley fluid a. To avoid being completely overwhelmed with free parameters, we only include coefficients at zeroth order in derivatives coupling distinct fluids together. In fact, this will be sufficient to capture the negative magnetoresistance, as we explain in the next section. Accounting for this coupling modifies the conservation equations to∇μJaμ=−Ca8εμνρσFμνFρσ−Ga32π2εμνρσRβμναRαρσβ−∑b[ℛabνb+Sabβb] [6a]∇μTaμν=FνμJμa−Ga16π2∇μ[ερσαβFρσRαβνμ]+uaν∑b[Uabνb+Vabβb], [6b]where we have defined βa≡1/Ta and νa≡βaμa. The transport problem is well-posed if∑aCa=∑aGa=0.[7]

The new coefficients ℛ,S,U, and V characterize the rate of the intervalley transfer of charge, energy, and momentum due to relative imbalances in chemical potential or temperature. In writing [6], we have chosen the intervalley scattering of energy and momentum to be relativistic. This makes the analysis easier as it preserves Lorentz covariance, but will not play an important role in our results. In particular, the intervalley momentum transfer processes are subleading effects in our theory of transport.

The gradient expansion may be different for each fluid, but we will assume that Jaμ and Taμν depend only on fluid a. We require that∑a or bℛab=∑a or bSab=∑a or bUab=∑a or bVab=0.[8]This ensures that globally charge and energy are conserved, as well as that uniform shifts in the background chemical potential and/or temperature, for all fluids simultaneously, are exact zero modes of the equations of motion.

For simplicity in [6], we have implicitly assumed that the Weyl nodes are all at the same chemical potential in equilibrium. This is generally not true for realistic Weyl materials. As nontrivial issues in hydrodynamics already arise without making this generalization, we will stick to the case where all Weyl nodes are at the same chemical potential in equilibrium in this paper.

For the remainder of this paper, we will be interested in transport in flat spacetimes where Rμναβ=0. Except where otherwise stated, we will assume Minkowski space from now on. Hence, for most purposes, we write partial derivatives ∂μ rather than covariant derivatives ∇μ. However, we will continue to use the covariant derivative ∇μ at intermediate steps of the calculations where it is necessary.

Thermodynamic Constraints

We will now derive the constraints on our hydrodynamic parameters which are imposed by demanding that the second law of thermodynamics is obeyed locally. Without intervalley coupling processes, and at the ideal fluid level (derivative corrections, including Fμν, are neglected), the second law of thermodynamics implies that the total entropy current sμ (where sa is the entropy density of fluid a) obeys (42, 43)∂μsμ=∂μ(∑asauaμ)=0.[9]In the more generic, nonideal, case the right-hand side of [9] must be nonnegative. In our theory of coupled chiral fluids, the right-hand side of [9] does not vanish already at the ideal fluid level:∂μsμ=∑ab(βa[Uabνb+Vabβb]+νa[ℛabνb+Sabβb])≥0.[10]There is no possible change we can make to the entropy current that is local and that can subtract off the right-hand side of [10]. Hence, we demand that the matrixA≡(ℛ −S−U V)[11]is positive semidefinite.

Using standard arguments for Onsager reciprocity in statistical mechanics (44), one can show that A=AT. In the SI Appendix, we will show using the memory matrix formalism (39, 45) that whenever the quantum mechanical operators na and ϵa are naturally defined:AIJ=Tlimω→0Im(Gx˙Ix˙JR(ω))ω,[12]where xI denotes (na,ϵa) and dots denote time derivatives. We also prove [8], and the symmetry and positive semidefiniteness of A through the memory matrix formalism, at the quantum mechanical level.

Equilibrium Fluid Flow

We now find an equilibrium solution to [6]. Beginning with the simple case of B=0, it is straightforward to see following ref. 46 that an equilibrium solution isμa=μ0(x),[13a]Ta=T0=constant,[13b]uaμ=(1,0).[13c]Indeed as pointed out in refs. 32, 46, this exactly satisfies [6] neglecting the intervalley and anomalous terms. Using [8] it is straightforward to see that the intervalley terms also vanish on this solution. If B=0 then CεμνρσFμνFρσ=0, and hence this is an exact solution to the hydrodynamic equations. We define the parameter ξ as the typical correlation length of μ(x): roughly speaking, ξ∼|μ0|/|∂xμ0|.

Following ref. 43, we can perturbatively construct a solution to the equations of motion when B≠0, assuming thatB≪T2 and 1≪ξT.[14]Both of these assumptions are necessary for our hydrodynamic formalism to be physically sensible. Using these assumptions, it is consistent at leading order to only change via≠0, but to keep μa and Ta the same:vza=Caμ02B2(ϵa+Pa)+GaT02Bϵa+Pa≡v(μ(x),T,Bi).[15]

It may seem surprising that in a single chiral fluid, there would be a nonvanishing charge current. This is a well-known phenomenon called the chiral magnetic effect (for a recent review, see ref. 47). In our model, the net current flow is the sum of the valley contributions:Jz=∑aJaz=∑aCaμ0B=0,[16]and so indeed, this complies with the expectation that the net current in a solid-state system will vanish in equilibrium, as discussed (in more generality) in refs. 41, 48.

Thermoelectric Conductivity

We now linearize the hydrodynamic equations around this equilibrium solution, applying infinitesimally small external electric fields E˜i, and temperature gradients ζ˜i=−∂i⁡log⁡T to the fluid. Although we have placed an equal sign in this equation, we stress that we will apply ζ˜i in such a way we may apply a constant temperature gradient on a compact space (with periodic boundary conditions). Applying a constant E˜i is simple, and corresponds to turning on an external electric field in Fμν. Applying a constant ζ˜i is more subtle, and can be done by changing the spacetime metric to (49):ds2=ημνdxμdxν−2e−iωt−iωζ˜idxidt,[17]ω is a regulator, which we take to 0 at the end of the calculation. This spacetime is flat (Rαβγδ=0). To account for both E˜i and ζ˜i, the external gauge field is modified to A+A˜, whereA˜i=−(E˜i−μ0(x)ζ˜i)e−iωt−iω.[18]

The hydrodynamic Eq. 6 must then be solved in this modified background. In linear response, the hydrodynamic variables becomeμa=μ0(x)+μ˜a(x),[19a]Ta=T0+T˜a(x),[19b]uaμ≈(1+viaζ˜ie−iωtiω,via+v˜ia).[19c]Note that tilded variables represent objects that are first order in linear response. The correction to uat is necessary to ensure that uμuμ=−1 is maintained. In general we cannot solve the linearized equations analytically, except in the limit of perturbatively weak disorder and magnetic field strength.

We assume that the inhomogeneity in the chemical potential is small:μ0=μ¯0+uμ^0(x),[20]with u≪μ¯0 and T; u is our perturbative parameter, and we assume that μ^0 is a zero-mean random function with unit variance. We assume the scalingsB∼u2  and   ℛ,S,U,V∼u6.[21]The hydrodynamic equations can be solved perturbatively, and the charge and heat currents may be spatially averaged on this perturbative solution. The computation is presented in the SI Appendix, and we present highlights here. At leading order, the linearized hydrodynamic equations reduce to∂i[naw˜ia+σQa(E˜i−∂iμ˜a−μ0T(Tζ˜i−∂iT˜a))]=CaE˜iBi−∑b[ℛabν˜b+Sabβ˜b][22a]∂i[Tsaw˜ia−μ0σQa(E˜i−∂iμ˜a−μ0T(Tζ˜i−∂iT˜a))]=2GaT02ζ˜iBi+∑b[(ℛabμ0+Uab)ν˜b+(Sabμ0+Vab)β˜b][22b]na(∂iμ˜a−E˜i)+sa(∂iT˜a−Tζ˜i)=εijkw˜janaBk.[22c]We have definedv˜ia=w˜ia+∂vi∂μμ˜a+∂vi∂TT˜a.[23]Here, w˜ia represents the fluid velocity after subtracting the contribution coming from [15] in local thermal equilibrium.

Eq. 22 depends on Ga, despite the fact that our spacetime is flat. This follows from the subtle fact that thermodynamic consistency of the anomalous quantum field theory on curved spacetimes requires that the axial–gravitational anomaly alters the thermodynamics of fluids on flat spacetimes (40).

The total charge current is J˜i=∑aJ˜ai, and the total heat current is Q˜i=∑aT˜ati−μ0J˜i. At leading order in perturbation theory, we find that the charge current in each valley fluid may be written asJ˜ai=naV˜ia+CaM˜aBi,[24]and the heat current per valley, T˜ati−μ0J˜ai, may be written asQ˜ai=TsaV˜ia+2GaT0T˜aBi.[25]In the above expressions V˜ia is a homogeneous O(u−2) contribution to w˜ia, and M˜a and T˜a are O(u−4) homogeneous contributions to μ˜a and T˜a, respectively.

The thermoelectric conductivity matrix is:σxx=σyy=∑ana2ΓaΓa2+B2na2,[26a]σxy=∑aBna3Γa2+B2na2,[26b]σzz=∑ana2Γa+sB2,[26c]κ¯xx=κ¯yy=∑aTsa2ΓaΓa2+B2na2,[26d]κ¯xy=∑aBTnasa2Γa2+B2na2,[26e]κ¯zz=∑aTsa2Γa+hB2,[26f]αxx=αyy=∑anasaΓaΓa2+B2na2,[26g]αxy=−αyx=∑aBna2saΓa2+B2na2,[26h]αzz=∑anasaΓa+aB2,[26i]where we have defined the four parameterss≡T(Ca Caμ)(ℛab −Sab−Uab Vab)−1(CbCbμ),[27a]h≡4T4(0 Ga)(ℛab −Sab−Uab Vab)−1(0Gb),[27b]a≡2T2(0 Ga)(ℛab −Sab−Uab Vab)−1(CbCbμ),[27c]Γa≡T02(sa(∂na/∂μ)−na(∂sa/∂μ))23σQa(ϵa+Pa)2u2.[27d]In these expressions, sums over valley indices are implicit. Coefficients odd under z→−z (such as σxz) vanish. Note that all of the contributions to the conductivities listed above are of the same order O(u−2) in our perturbative expansion, explaining the particular scaling limit (21) in u that was taken.

We have not listed the full set of transport coefficients. The unlisted transport coefficients are related to those in ref. 26 by Onsager reciprocity:σij(B)=σji(−B),[28a]κ¯ij(B)=κ¯ji(−B),[28b]αij(B)=α¯ji(−B).[28c]The symmetry of A is crucial in order for the final conductivity matrix to obey [28].

Evidently, the conductivities perpendicular to the magnetic field are Drude-like. This follows from principles, which are by now very well understood (31, 39). In these weakly disordered fluids, the transport coefficients are only limited by the rate at which momentum relaxes due to the disordered chemical potential Γa/(ϵa+Pa), and/or the rate at which the magnetic field relaxes momentum (by “rotating” it in the xy plane), Bna/(ϵa+Pa). This latter energy scale is the hydrodynamic cyclotron frequency (31). In our hydrodynamic theory, we can see this momentum “bottleneck” through the fact that the components of the charge and heat currents in [24] and [25], perpendicular to Bi, are proportional to the same fluid velocity V˜ia∼u−2 at leading order. The transport of the fluid is dominated by the slow rate at which this large velocity can relax. Because the heat current and charge current are proportional to this velocity field, the contribution of each valley fluid to σ, α, and κ¯ are all proportional to one another in the xy plane.

The remaining nonvanishing transport coefficients are σzz, αzz, and κ¯zz. From [26], we see that these conductivities are a sum of a Drude-like contribution (because this is transport parallel to the magnetic field, there is no magnetic momentum relaxation) from each valley, as before, and a new “anomalous” contribution that couples the valley fluids together. This anomalous contribution has a qualitatively similar origin as that discovered in refs. 8, 36. It can crudely be understood as follows: the chemical potential and temperature imbalances M˜ and T˜ are proportional to B and inversely proportional to A, as the homogeneous contributions to the right-hand side of [22] cancel. Such thermodynamic imbalances lead to corrections to valley fluid charge and heat currents, analogous to the chiral magnetic effect—these are the linear in B terms in [24] and [25]. Combining these scalings together immediately gives us the qualitative form of the anomalous contributions to the conductivity matrix.

The positive semidefiniteness of the thermoelectric conductivity matrix is guaranteed. Thinking of the conductivity matrices as a sum of the anomalous contribution and Drude contributions for each valley, it suffices to show each piece is positive-definite individually. The Drude pieces are manifestly positive definite, as is well known (it is an elementary exercise in linear algebra to confirm). To show the anomalous pieces are positive semidefinite, it suffices to show that sh≥Ta2. This follows from [27], and the Cauchy–Schwarz inequality (v1TAv1)(v2TAv2)≥(v1TAv2)2 for any vectors v1,2, and a symmetric, positive-semidefinite matrix A. These arguments also guarantee s,h>0.

Our expression for the conductivity may seem ill posed—it explicitly depends on the matrix inverse A−1, but A is not invertible. In fact, the kernel of A has two linearly independent vectors: (1a,1a) and (1a,−1a), with 1a denoting a vector with valley indices with each entry equal to 1. However, in the final formula for the thermoelectric conductivities, which sums over all valley fluids, we see that A−1 is contracted with vectors that are orthogonal to the kernel of A due to [7]. The expression for the conductivities is therefore finite and unique.

We present a simple example of our theory for a fluid with two identical Weyl nodes of opposite chirality in the SI Appendix, along with a demonstration that the equations of motion are unchanged when we account for long-range Coulomb interactions, or impose electric fields and temperature gradients through boundary conditions in a finite domain. Hence, the transport coefficients we have computed above are in fact those which will be measured in experiment.

In this paper, we used inhomogeneity in the chemical potential to relax momentum when B=0. By following the hydrodynamic derivation in ref. 46, other mechanisms for disorder likely lead to the same thermoelectric conductivities as reported in ref. 26, but with a different formula for Γa.

Violation of the Wiedemann–Franz Law

The thermal conductivity κij usually measured in experiments is defined with the boundary conditions J˜i=0 (as opposed to κ¯ij, which is defined with E˜i=0). This thermal conductivity is related to the elements of the transport matrix [1] byκij=κ¯ij−Tα¯ikσkl−1αlj.[29]In an ordinary metal, the Wiedemann–Franz (WF) law states that (50)ℒij≡κijTσij=π23.[30]The numerical constant of π2/3 comes from the assumption that the quasiparticles are fermions, and that the dominant interactions are between quasiparticles and phonons or impurities, but otherwise is robust to microscopic details.

In general, our model will violate the WF law. Details of this computation are provided in the SI Appendix. In general, the WF law is violated by an O(1) constant, which depends on the magnetic field B. However, in the special case where we have valley fluids of opposite chirality but otherwise identical equations of state, we find that the transverse Lorenz ratios ℒxx, ℒxy, ℒyy are all parametrically smaller than 1 (in fact, they vanish at leading order in perturbation theory). In contrast, we find that ℒzz∼B2 at small B, and saturates to a finite number as B becomes larger (but still ≪T2). This dramatic angular dependence of the WF law would be a sharp experimental test of our formalism in a strongly correlated Weyl material.

If the intervalley scattering rate is almost vanishing, the anomalous conductivities of a weakly interacting Weyl gas are still computable with our formalism. Weak intravalley scattering processes bring the “Fermi liquid” at each Weyl node to thermal equilibrium, and A may be computed via semiclassical kinetic theory. Assuming that intervalley scattering occurs elastically off of pointlike impurities, we compute s, a, and h in the SI Appendix. We find that ℒzzanom<π2/3, asymptotically approaching the WF law when μ≫T. This is in contrast to the nonanomalous conductivities of a semiconductor, where under similar assumptions ℒzz>π2/3 (51). The Mott relation between αzzanom and σzzanom differs by an overall sign from the standard relation. These discrepancies occur because we have assumed that elastic intervalley scattering is weaker than intravalley thermalization. In contrast, [11] makes the opposite assumption when μ≫T, and so recovers all ordinary metallic phenomenology. Even in this limit, negative thermal magnetoresistance is a consequence of nonvanishing Ga.

Outlook

In this paper, we have systematically developed a hydrodynamic theory of thermoelectric transport in a Weyl semimetal where quasiparticle–quasiparticle scattering is faster than quasiparticle–impurity and/or quasiparticle–phonon scattering. We have demonstrated the presence of longitudinal negative magnetoresistance in all thermoelectric conductivities. New phenomenological parameters introduced in our classical model may be directly computed using the memory matrix formalism given a microscopic quantum mechanical model of a Weyl semimetal. Our formalism is directly applicable to microscopic models of interacting Weyl semimetals where all relevant nodes are at the same Fermi energy. Our model should be generalized to the case where different nodes are at different Fermi energies, though our main results about the nature of negative magnetoresistance likely do not change qualitatively.

Previously, exotic proposals had been put forth to measure the axial–gravitational anomaly in an experiment. Measurements involving rotating cylinders of a Weyl semimetal have been proposed in refs. 41, 52, and it is possible that the rotational speed of neutron stars is related to this anomaly (53). A nonvanishing negative magnetoresistance in either αzz or κ¯zz is a direct experimental signature of the axial–gravitational anomaly. It is exciting that a relatively mundane transport experiment on a Weyl semimetal is capable of detecting this novel anomaly.

Acknowledgments

R.A.D. is supported by the Gordon and Betty Moore Foundation EPiQS Initiative through Grant GBMF4306. A.L. and S.S. are supported by the National Science Foundation under Grant DMR-1360789 and Multi-University Research Initiative Grant W911NF-14-1-0003 from Army Research Office. Research at Perimeter Institute is supported by the Government of Canada through Industry Canada and by the Province of Ontario through the Ministry of Research and Innovation. S.S. also acknowledges support from Cenovus Energy at Perimeter Institute.

Footnotes

  • ↵1To whom correspondence may be addressed. Email: sachdev{at}g.harvard.edu, lucas{at}fas.harvard.edu, or rdavison{at}physics.harvard.edu.
  • Author contributions: A.L., R.A.D., and S.S. designed research; A.L., R.A.D., and S.S. performed research; and A.L. and R.A.D. wrote the paper.

  • Reviewers: A.B., University of Waterloo; and K.J., San Francisco State University.

  • The authors declare no conflict of interest.

  • This article contains supporting information online at www.pnas.org/lookup/suppl/doi:10.1073/pnas.1608881113/-/DCSupplemental.

View Abstract

References

  1. ↵
    1. Wan X,
    2. Turner AM,
    3. Vishwanath A,
    4. Savrasov SY
    (2011) Topological semimetal and Fermi-arc surface states in the electronic structure of pyrochlore iridates. Phys Rev B 83:205101.
    .
    OpenUrlCrossRef
  2. ↵
    1. Fang C,
    2. Gilbert MJ,
    3. Dai X,
    4. Bernevig BA
    (2012) Multi-Weyl topological semimetals stabilized by point group symmetry. Phys Rev Lett 108(26):266802.
    .
    OpenUrlCrossRefPubMed
  3. ↵
    1. Weng H,
    2. Fang C,
    3. Fang Z,
    4. Bernevig A,
    5. Dai X
    (2015) Weyl semimetal phase in non-centrosymmetric transition metal monophosphides. Phys Rev X 5:011029.
    .
    OpenUrl
  4. ↵
    1. Lu L, et al.
    (2015) Topological matter. Experimental observation of Weyl points. Science 349(6248):622–624.
    .
    OpenUrlAbstract/FREE Full Text
  5. ↵
    1. Xu SY, et al.
    (2015) Topological matter. Discovery of a Weyl fermion semimetal and topological Fermi arcs. Science 349(6248):613–617.
    .
    OpenUrlAbstract/FREE Full Text
  6. ↵
    1. Lv BQ, et al.
    (2015) Experimental discovery of Weyl semimetal TaAs. Phys Rev X 5:031013.
    .
    OpenUrlCrossRef
  7. ↵
    1. Nielsen HN,
    2. Ninomiya M
    (1983) The Adler-Bell-Jackiw anomaly and Weyl fermions in a crystal. Phys Lett B 130:389–396.
    .
    OpenUrlCrossRef
  8. ↵
    1. Son DT,
    2. Spivak BZ
    (2013) Chiral anomaly and classical negative magnetoresistance of Weyl metals. Phys Rev B 88:104412.
    .
    OpenUrlCrossRef
  9. ↵
    1. Burkov AA
    (2015) Negative longitudinal magnetoresistance in Dirac and Weyl metals. Phys Rev B 91:245157.
    .
    OpenUrlCrossRef
  10. ↵
    1. Lundgren R,
    2. Laurell P,
    3. Fiete GA
    (2014) Thermoelectric properties of Weyl and Dirac semimetals. Phys Rev B 90:165115.
    .
    OpenUrlCrossRef
  11. ↵
    1. Spivak BZ,
    2. Andreev AV
    (2016) Magneto-transport phenomena related to the chiral anomaly in Weyl semimetals. Phys Rev B 93:085107.
    .
    OpenUrlCrossRef
  12. ↵
    1. Kim HJ, et al.
    (2013) Dirac versus Weyl fermions in topological insulators: Adler-Bell-Jackiw anomaly in transport phenomena. Phys Rev Lett 111(24):246603.
    .
    OpenUrlCrossRefPubMed
  13. ↵
    1. Xiong J, et al.
    (2015) Evidence for the chiral anomaly in the Dirac semimetal Na3Bi. Science 350(6259):413–416.
    .
    OpenUrlAbstract/FREE Full Text
  14. ↵
    1. Huang X, et al.
    (2015) Observation of the chiral anomaly induced negative magnetoresistance in 3D Weyl semimetal TaAs. Phys Rev X 5:031023.
    .
    OpenUrl
  15. ↵
    1. Li CZ, et al.
    (2015) Giant negative magnetoresistance induced by the chiral anomaly in individual Cd3As2 nanowires. Nat Commun 6:10137.
    .
    OpenUrlCrossRefPubMed
  16. ↵
    1. Li H, et al.
    (2016) Negative magnetoresistance in Dirac semimetalCd3As2. Nat Commun 7:10301.
    .
    OpenUrlCrossRefPubMed
  17. ↵
    1. Zhang CL, et al.
    (2016) Signatures of the Adler-Bell-Jackiw chiral anomaly in a Weyl fermion semimetal. Nat Commun 7:10735.
    .
    OpenUrlCrossRefPubMed
  18. ↵
    1. Hirschberger M, et al.
    (2016) The chiral anomaly and thermopower of Weyl fermions in the half-Heusler GdPtBi. arXiv:1602.07219.
    .
  19. ↵
    1. Gurzhi RN
    (1963) Minimum of resistance in impurity-free conductors. J Exp Theor Phys 17:521.
    .
    OpenUrl
  20. ↵
    1. Spivak B,
    2. Kivelson SA
    (2006) Transport in two-dimensional electronic micro-emulsions. Ann Phys 321:2079.
    .
    OpenUrl
  21. ↵
    1. Andreev AV,
    2. Kivelson SA,
    3. Spivak B
    (2011) Hydrodynamic description of transport in strongly correlated electron systems. Phys Rev Lett 106(25):256804.
    .
    OpenUrlCrossRefPubMed
  22. ↵
    1. Tomadin A,
    2. Vignale G,
    3. Polini M
    (2014) Corbino disk viscometer for 2D quantum electron liquids. Phys Rev Lett 113(23):235901.
    .
    OpenUrlCrossRefPubMed
  23. ↵
    1. Principi A,
    2. Vignale G
    (2015) Violation of the Wiedemann-Franz law in hydrodynamic electron liquids. Phys Rev Lett 115(5):056603.
    .
    OpenUrlCrossRefPubMed
  24. ↵
    1. Torre I,
    2. Tomadin A,
    3. Geim AK,
    4. Polini M
    (2015) Nonlocal transport and the hydrodynamic shear viscosity in graphene. Phys Rev B 92:165433.
    .
    OpenUrlCrossRef
  25. ↵
    1. Levitov L,
    2. Falkovich G
    (2016) Electron viscosity, current vortices and negative nonlocal resistance in graphene. Nat Phys 12:672–676.
    .
    OpenUrlCrossRef
  26. ↵
    1. Molenkamp LW,
    2. Molenkamp LW,
    3. de Jong MJ
    (1995) Hydrodynamic electron flow in high-mobility wires. Phys Rev B Condens Matter 51(19):13389–13402.
    .
    OpenUrlCrossRefPubMed
  27. ↵
    1. Bandurin DA, et al.
    (2016) Negative local resistance caused by viscous electron backflow in graphene. Science 351(6277):1055–1058.
    .
    OpenUrlAbstract/FREE Full Text
  28. ↵
    1. Moll PJW,
    2. Kushwaha P,
    3. Nandi N,
    4. Schmidt B,
    5. Mackenzie AP
    (2016) Evidence for hydrodynamic electron flow in PdCoO₂. Science 351(6277):1061–1064.
    .
    OpenUrlAbstract/FREE Full Text
  29. ↵
    1. Crossno J, et al.
    (2016) Observation of the Dirac fluid and the breakdown of the Wiedemann-Franz law in graphene. Science 351(6277):1058–1061.
    .
    OpenUrlAbstract/FREE Full Text
  30. ↵
    1. Ghahari F, et al.
    (2016) Enhanced thermoelectric power in graphene: Violation of the Mott relation by inelastic scattering. Phys Rev Lett 116(13):136802.
    .
    OpenUrlCrossRefPubMed
  31. ↵
    1. Hartnoll SA,
    2. Kovtun PK,
    3. Müller M,
    4. Sachdev S
    (2007) Theory of the Nernst effect near quantum phase transitions in condensed matter, and in dyonic black holes. Phys Rev B 76:144502.
    .
    OpenUrlCrossRef
  32. ↵
    1. Lucas A,
    2. Crossno J,
    3. Fong KC,
    4. Kim P,
    5. Sachdev S
    (2016) Transport in inhomogeneous quantum critical fluids and in the Dirac fluid in graphene. Phys Rev B 93:075426.
    .
    OpenUrlCrossRef
  33. ↵
    1. Lucas A
    (2016) Sound waves and resonances in electron-hole plasma. Phys Rev B 93:245153.
    .
    OpenUrlCrossRef
  34. ↵
    1. Son DT,
    2. Surówka P
    (2009) Hydrodynamics with triangle anomalies. Phys Rev Lett 103(19):191601.
    .
    OpenUrlCrossRefPubMed
  35. ↵
    1. Neiman Y,
    2. Oz Y
    (2011) Relativistic hydrodynamics with general anomalous charges. J High Energy Phys 2011(3):023.
    .
    OpenUrl
  36. ↵
    1. Landsteiner K,
    2. Liu Y,
    3. Sun YW
    (2015) Negative magnetoresistivity in chiral fluids and holography. J High Energy Phys 2015(3):127.
    .
    OpenUrl
  37. ↵
    1. Sadofyev AV,
    2. Isachenkov MV
    (2011) The Chiral magnetic effect in hydrodynamical approach. Phys Lett B 697:404–406.
    .
    OpenUrlCrossRef
  38. ↵
    1. Roychowdhury D
    (2015) Magnetoconductivity in chiral Lifshitz hydrodynamics. J High Energy Phys 2015(9):145.
    .
    OpenUrl
  39. ↵
    1. Lucas A,
    2. Sachdev S
    (2015) Memory matrix theory of magnetotransport in strange metals. Phys Rev B 91:195122.
    .
    OpenUrlCrossRef
  40. ↵
    1. Jensen K,
    2. Loganayagam R,
    3. Yarom A
    (2013) Thermodynamics, gravitational anomalies and cones. J High Energy Phys 2013(2):088.
    .
    OpenUrl
  41. ↵
    1. Landsteiner K
    (2014) Anomaly related transport of Weyl fermions for Weyl semi-metals. Phys Rev B 89:075124.
    .
    OpenUrlCrossRef
  42. ↵
    1. Jensen K, et al.
    (2012) Towards hydrodynamics without an entropy current. Phys Rev Lett 109(10):101601.
    .
    OpenUrlCrossRefPubMed
  43. ↵
    1. Banerjee N, et al.
    (2012) Constraints on fluid dynamics from equilibrium partition functions. J High Energy Phys 2012(9):046.
    .
    OpenUrl
  44. ↵
    1. Landau LD,
    2. Lifshitz EM
    (1980) Statistical Physics Part 1 (Butterworth Heinemann, Oxford).
    .
  45. ↵
    1. Forster D
    (1975) Hydrodynamic Fluctuations, Broken Symmetry and Correlation Functions (Perseus Books, Reading, MA).
    .
  46. ↵
    1. Lucas A
    (2015) Hydrodynamic transport in strongly coupled disordered quantum field theories. New J Phys 17:113007.
    .
    OpenUrlCrossRef
  47. ↵
    1. Kharzeev DE
    (2014) The chiral magnetic effect and anomaly-induced transport. Prog Part Nucl Phys 75:133–151.
    .
    OpenUrlCrossRef
  48. ↵
    1. Vazifeh MM,
    2. Franz M
    (2013) Electromagnetic response of Weyl semimetals. Phys Rev Lett 111(2):027201.
    .
    OpenUrlCrossRefPubMed
  49. ↵
    1. Hartnoll SA
    (2009) Lectures on holographic methods for condensed matter physics. Class Quantum Gravity 26:224002.
    .
    OpenUrlCrossRef
  50. ↵
    1. Ashcroft NW,
    2. Mermin ND
    (1976) Solid-State Physics (Brooks Cole, Belmont, CA).
    .
  51. ↵
    1. Nolas GS
    (2004) Thermal conductivity of semiconductors. Thermal Conductivity: Theory, Property and Applications, ed Tritt TM (Kluwer Academic, New York), pp 105–121.
    .
  52. ↵
    1. Chernodub MN,
    2. Cortijo A,
    3. Grushin AG,
    4. Landsteiner K,
    5. Vozmediano MAH
    (2014) Condensed matter realization of the axial magnetic effect. Phys Rev B 89:081407.
    .
    OpenUrlCrossRef
  53. ↵
    1. Kaminski M,
    2. Uhlemann CF,
    3. Bleicher M,
    4. Schaffner-Bielich J
    (2014) Anomalous hydrodynamics kicks neutron stars. arXiv:1410.3833.
    .
PreviousNext
Back to top
Article Alerts
Email Article

Thank you for your interest in spreading the word on PNAS.

NOTE: We only request your email address so that the person you are recommending the page to knows that you wanted them to see it, and that it is not junk mail. We do not capture any email address.

Enter multiple addresses on separate lines or separate them with commas.
Hydrodynamic theory of thermoelectric transport and negative magnetoresistance in Weyl semimetals
(Your Name) has sent you a message from PNAS
(Your Name) thought you would like to see the PNAS web site.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Citation Tools
Hydrodynamic theory of Weyl semimetals
Andrew Lucas, Richard A. Davison, Subir Sachdev
Proceedings of the National Academy of Sciences Aug 2016, 113 (34) 9463-9468; DOI: 10.1073/pnas.1608881113

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Request Permissions
Share
Hydrodynamic theory of Weyl semimetals
Andrew Lucas, Richard A. Davison, Subir Sachdev
Proceedings of the National Academy of Sciences Aug 2016, 113 (34) 9463-9468; DOI: 10.1073/pnas.1608881113
Digg logo Reddit logo Twitter logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Mendeley logo Mendeley
Proceedings of the National Academy of Sciences: 113 (34)
Table of Contents

Submit

Sign up for Article Alerts

Article Classifications

  • Physical Sciences
  • Physics

Jump to section

  • Article
    • Abstract
    • Weyl Hydrodynamics
    • Thermodynamic Constraints
    • Equilibrium Fluid Flow
    • Thermoelectric Conductivity
    • Violation of the Wiedemann–Franz Law
    • Outlook
    • Acknowledgments
    • Footnotes
    • References
  • Figures & SI
  • Info & Metrics
  • PDF

You May Also be Interested in

Abstract depiction of a guitar and musical note
Science & Culture: At the nexus of music and medicine, some see disease treatments
Although the evidence is still limited, a growing body of research suggests music may have beneficial effects for diseases such as Parkinson’s.
Image credit: Shutterstock/agsandrew.
Large piece of gold
News Feature: Tracing gold's cosmic origins
Astronomers thought they’d finally figured out where gold and other heavy elements in the universe came from. In light of recent results, they’re not so sure.
Image credit: Science Source/Tom McHugh.
Dancers in red dresses
Journal Club: Friends appear to share patterns of brain activity
Researchers are still trying to understand what causes this strong correlation between neural and social networks.
Image credit: Shutterstock/Yeongsik Im.
Yellow emoticons
Learning the language of facial expressions
Aleix Martinez explains why facial expressions often are not accurate indicators of emotion.
Listen
Past PodcastsSubscribe
Goats standing in a pin
Transplantation of sperm-producing stem cells
CRISPR-Cas9 gene editing can improve the effectiveness of spermatogonial stem cell transplantation in mice and livestock, a study finds.
Image credit: Jon M. Oatley.

Similar Articles

Site Logo
Powered by HighWire
  • Submit Manuscript
  • Twitter
  • Facebook
  • RSS Feeds
  • Email Alerts

Articles

  • Current Issue
  • Special Feature Articles – Most Recent
  • List of Issues

PNAS Portals

  • Anthropology
  • Chemistry
  • Classics
  • Front Matter
  • Physics
  • Sustainability Science
  • Teaching Resources

Information

  • Authors
  • Editorial Board
  • Reviewers
  • Librarians
  • Press
  • Site Map
  • PNAS Updates

Feedback    Privacy/Legal

Copyright © 2021 National Academy of Sciences. Online ISSN 1091-6490