Hydrodynamic Electron Flow and Hall Viscosity

Hydrodynamic Electron Flow and Hall Viscosity

Thomas Scaffidi Department of Physics, University of California, Berkeley, CA 94720, USA    Nabhanila Nandi Max Planck Institute for Chemical Physics of Solids, Nöthnitzer Strasse 40, 01187 Dresden, Germany    Burkhard Schmidt Max Planck Institute for Chemical Physics of Solids, Nöthnitzer Strasse 40, 01187 Dresden, Germany    Andrew P. Mackenzie Max Planck Institute for Chemical Physics of Solids, Nöthnitzer Strasse 40, 01187 Dresden, Germany Scottish Universities Physics Alliance, School of Physics and Astronomy, University of St. Andrews, St. Andrews, Fife KY16 9SS, UK    Joel E. Moore Department of Physics, University of California, Berkeley, CA 94720, USA Materials Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, CA 94720
September 12, 2019
Abstract

In metallic samples of small enough size and sufficiently strong momentum-conserving scattering, the viscosity of the electron gas can become the dominant process governing transport. In this regime, momentum is a long-lived quantity whose evolution is described by an emergent hydrodynamical theory. Furthermore, breaking time-reversal symmetry leads to the appearance of an odd component to the viscosity called the Hall viscosity, which has attracted considerable attention recently due to its quantized nature in gapped systems but still eludes experimental confirmation. Based on microscopic calculations, we discuss how to measure the effects of both the even and odd components of the viscosity using hydrodynamic electronic transport in mesoscopic samples under applied magnetic fields.

pacs:

The semiclassical theory of electronic conduction, based on relaxation of total momentum by impurities, phonons and umklapp scattering, occupies a central place in condensed matter physics. It is therefore of particular interest to study the cases for which it fails. One case that has attracted much interest is the possibility of a hydrodynamic regime, where transport is dominated by viscous effects Gurzhi (1963, 1968); Black (1980); Yu et al. (1984); Gurzhi et al. (1995); Molenkamp and De Jong (1994); Molenkamp and de Jong (1994); de Jong and Molenkamp (1995); Yu et al. (1984); Dyakonov and Shur (1993); Chow et al. (1996); Spivak and Kivelson (2006); Andreev et al. (2011); Torre et al. (2015); Alekseev (2016); Tomadin et al. (2014); Levitov and Falkovich (2016); Torre et al. (2015); Sherafati et al. (2016); Guo et al. (2016a, b); Lucas (2016); Levchenko et al. (2016); Ganeshan and Abanov (2017). One needs a large separation of scales between momentum-relaxing and momentum-conserving scattering in order to see these effects. This was recently achieved in grapheneBandurin et al. (2016); Crossno et al. (2016) and PdCoOMoll et al. (2016); Mackenzie (2017).

Interest in such a hydrodynamic regime also emanated from a conjectured bound on diffusion constants for the hydrodynamics of strongly interacting quantum systems Kovtun et al. (2005); Hartnoll (2015). Even though the physics described in this work is semiclassical and probably still quite far from these quantum-mechanical bounds, the observations that we hope to stimulate would constitute an important first step towards the understanding of emergent hydrodynamical regimes in electronic systems.

A further motivation for the work is that reaching a viscous regime for a charged fluid enables one to break time-reversal symmetry by adding a magnetic field and hence to study a non-dissipative component to the viscosity tensor called the Hall viscosity. The recent interest in this Hall viscosity emanates from the fact that it is topologically quantized in gapped systems Read (2009). In order to study this effect experimentally, in analogy with the Hall conductivity, the first step would obviously be to measure the classical Hall viscosity. We show in this letter how this measurement could be done by describing specific size effects from Hall viscosity in transport in restricted 2D channels under transverse magnetic fields.

This paper is organized as follows. We start by assuming a perfect hydrodynamic regime and calculate and . We show that the component of is proportional to the Hall viscosity, thereby providing a way of measuring it. In order to have realistic predictions to compare with experiments, one should also take into account other, non-viscous effects that can lead to a size-dependent resistivity. We thus perform a kinetic Boltzmann calculation in which the effect of diffuse boundaries, gradient along the section of the wire, momentum-conserving scattering, and magnetic field are taken into account. We show that the size effects in resistivities with and without momentum-conserving scattering are markedly different, thereby making it possible to distinguish hydrodynamic and non-hydrodynamic size effects. Finally, we comment on how these measurements could also enable one to measure the quantum Hall viscosity in the quantum Hall regime and establish a relation similar to Hoyos-Son Hoyos and Son (2012).

I Fluid equation

Even though a finite amount of momentum relaxation is always present, it is instructive to first look at the limit where it is zero. In this limit, the momentum density of the electron gas is conserved and one can write a hydrodynamic equation to model its dynamics Alekseev (2016) 111Following Ref. Alekseev (2016), we consider the case of an isochoric flow, . Even though our solution exhibits non-zero variations in the local electron density compared to the average density, the pressure contribution coming from those variations is negligible.:

(1)

where is the electron mass and and are the regular and Hall components of the kinematic viscosity tensor. As mentioned previously, the two diffusion constants in this equation are interesting from a fundamental point of view: (1) from a holographic argument, a bound on was conjectured (coming from the bound on the dissipative time scale ) and (2) is non-dissipative (and therefore not subject to this bound) but was shown to be quantized in gapped systems Avron et al. (1995); Read and Rezayi (2011); Hughes et al. (2013); Hoyos (2014).

We consider the case of a two-dimensional channel of width along the direction, with a uniform applied electric field along , uniform magnetic field along , and zero current along . In the stationary regime, one finds

(2)

where . Using the conventional no-slip boundary condition (we will treat the more realistic case of diffuse boundaries within a kinetic formalism in the next section), one finds

(3)

where . In the wide sample limit (), one has , since the only source of resistance comes from the boundaries, and , which is indeed its bulk value. We propose to measure and by measuring the size dependence of and in restricted channels of varying size and under varying magnetic fields.

In the next section, we will compare these results with the results of a kinetic Boltzman formalism. In order to do so, it is convenient to inject a microscopically derived magnetic field dependence of the viscosities in the above hydrodynamic solutions. The following dependence was found by Alekseev Alekseev (2016):

(4)

where is the cyclotron radius, is the momentum-conserving scattering length and where Alekseev (2016). This leads to

(5)

The size effect on is maximal at zero field due to the Lorentzian factor. In this limit, one obtains

(6)

In order to measure this effect, should be as small as possible for it to be sizable, but still somewhat larger than in order to remain in the hydrodynamic regime. For example, for , one expects a relative change of the Hall slope at zero field of the order of , which should be measurable.

Ii Kinetic theory

As mentioned previously, in order to make a quantitative comparison with experiments, it is crucial to go beyond a purely hydrodynamical theory and take into account several other effects like the non-zero momentum relaxation and the diffuse scattering at the boundaries. In order to do this, we perform kinetic Boltzmann calculations Ditlefsen and Lothe (1966); de Jong and Molenkamp (1995); Moll et al. (2016):

(7)

where is the external magnetic field, is the semiclassical occupation number for a wavepacket at position and velocity and where

(8)

with the local charge density, the local current, the momentum average and . For the sake of simplicity, we consider the case of a circular Fermi surface with with the radial unit vector. Scattering lengths are then simply defined as . The term proportional to is the most simple momentum-conserving scattering term that can be written assuming that the electrons relax to a Fermi-Dirac distribution shifted by the drift velocity Molenkamp and de Jong (1994); Moll et al. (2016). The boundary conditions are given by

(9)

which imposes, respectively, no current in the direction and diffuse scattering at the boundaries. Equation 7 is supplemented by Gauss’s law with a charge density given by . The resulting integrodifferential equation is advantageously solved numerically by using the method of characteristics Courant and Hilbert (2008).

Three limiting regimes can be identified de Jong and Molenkamp (1995): Ohmic, hydrodynamic and ballistic. In the Ohmic case, one has , and transport is therefore dominated by momentum-relaxing scattering, leading to bulk values for transport coefficients: and . In the ballistic case, one has , and transport is dominated by scattering with the boundaries. This leads to strong size effects, but with a different size dependence from the one arising in the hydrodynamic regime. Finally, in the hydrodynamic case, one has , and transport is dominated by the diffusion of momentum at the boundaries through viscosity. Note that this regime will only appear if .

Using the above kinetic theory, we look at transport in and between these different regimes and identify clear signatures of a hydrodynamic regime. In figures 1 and 2, we give respectively the magnetoresistance and the Hall resistivity in a typical ballistic case (, ) and a typical hydrodynamic case (, ).

In the ballistic regime (Fig 1.a), the magnetoresistance shows a maximum around and a rapid change of slope at , as previously reported Ditlefsen and Lothe (1966); Beenakker and van Houten (1991). In this ballistic regime, transport is dominated by electrons with a velocity close to the longitudinal () direction for which scattering on boundaries are very rare. These trajectories are bent by the field, leading to an increase in boundary scattering and therefore an increase of at low fields. At higher fields, when the cyclotron radius becomes of the order of , a larger and larger fraction of electrons present near the middle of the wire stop seeing the boundaries at all, and at large fields the bulk resistivity is therefore recovered.

In contrast, in the hydrodynamic regime (Fig 1.b), the longitudinal magnetoresistance decays monotonically and there is no sharp slope change at . The magnetoresistance follows a Lorentzian-like curve, in agreement with the hydrodynamic calculation of Eq. 5. This absence of sharp behavior at shows that, in this case, momentum-conserving scattering is dominant and no electron can ever go over a full cyclotron orbit without being scattered. As explicit in Eq. 5, the relevant dimensionless parameter is therefore now , as opposed to in the ballistic case. The absence of a maximum of the magnetoresistance and of a sharp slope change at can be used as clear signatures of the hydrodynamic regime. Note that the different scaling of with in zero applied field can also be used to differentiate the two regimes, as reported in Moll et al. (2016).

The Hall resistivity also shows strikingly different behaviors in the ballistic and hydrodynamic regimes. As seen in Fig. 2.b, in the ballistic case, exhibits a minimum at and a sharp slope change at , in analogy with and in agreement with the perturbative calculation given in Ref. Sondheimer (1950). While, in the ballistic case, is positive at low fields and negative at large fields, it is always negative in the hydrodynamic case (see Fig 2.b). The sign of can therefore be used a clear signature of hydrodynamic effects: one would expect a smaller (resp. larger) Hall slope at small fields than at large fields in the hydrodynamic (resp. ballistic) regime.

It was checked that, in both the cases of and , if becomes larger than , Ohmic behavior is recovered. It was also checked that, in the case of , if becomes smaller than , ballistic behavior is recovered.

a)
b)
Figure 1: Magnetoresistance for the ballistic () (a) and hydrodynamic case () (b).
a)
b)
Figure 2: Hall resistivity for the ballistic () (a) and hydrodynamic case () (b). The inset of panel (a) shows the saturation of the Hall resistivity at low fields in the ballistic case.

Iii Quantum limit

While the calculations up to this point have been classical, it is instructive to look at the large magnetic field limit, . In this limit, one finds from Eq. 4 that

(10)

which, crucially, does not depend on . Now, using , with and , one finds

(11)

which gives the quantized value of the dynamical viscosity for a quantum Hall system at filling . Inverting the resistivity tensor obtained in Eq. 3 and taking the large field limit leads to the conductivity:

(12)

The large magnetic field extrapolation of the classical calculation found in this paper is therefore consistent with the form of the finite- correction found by Hoyos and Son in the quantum case Hoyos and Son (2012) 222Our classical calculation deals with compressible fluids and is therefore not a priori expected to apply to incompressible quantum Hall fluids. Yet, out compressible calculation could still apply to a quantum Hall system with , which is effectively a gapless system because of finite-size effects. Interestingly, this is the regime in which our calculations predict a sizable viscous contribution to .. In the quantum Hall regime, the relative deviation of from its bulk value is of the order of . For realistic system sizes, the factor would probably be too small to lead to a measurable deviation. One could then turn to optical probes where could be chosen in the X-ray range.

Iv Conclusion

In conclusion, we have identified clear, qualitative signatures of hydrodynamic behavior in transport measurements of mesoscopic metallic samples under magnetic fields. This type of measurement is possible with readily available experimental techniques, and would make it possible to measure the classical Hall viscosity of the electron gas, which to the best of our knowledge has never been measured in a solid-state system. This would both further the evidence of a hydrodynamic regime in electronic transport and constitute an important step towards an experimental understanding of the quantum Hall viscosity.

V Acknowledgements

Helpful conversations with Tankut Can and Sriram Ganeshan are acknowledged. The authors acknowledge support from the Emergent Phenomena in Quantum Systems initiative of the Gordon and Betty Moore Foundation (T. S.) and NSF DMR-1507141 and a Simons Investigatorship (J.E.M.). We also acknowledge the support of the Max Planck Society and the UK Engineering and Physical Sciences Research Council under grant EP/I032487/1.

References

Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
""
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
   
Add comment
Cancel
Loading ...
338468
This is a comment super asjknd jkasnjk adsnkj
Upvote
Downvote
""
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters
Submit
Cancel

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test
Test description