Abstract
Electron hydrodynamics typically emerges in electron fluids with a high electronâelectron collision rate. However, new experiments with thin flakes of WTe2 have revealed that other momentum-conserving scattering processes can replace the role of the electronâelectron interaction, thereby leading to a novel, so-called para-hydrodynamic regime. Here, we develop the kinetic theory for para-hydrodynamic transport. To this end, we consider a ballistic electron gas in a thin three-dimensional sheet where the momentum-relaxing (lmr) and momentum-conserving (lmc) mean free paths are decreased due to boundary scattering from a rough surface. The resulting effective mean free path of the in-plane components of the electronic flow is then expressed in terms of microscopic parameters of the sheet boundaries, predicting that a para-hydrodynamic regime with lmrââ«âlmc emerges generically in ultraclean three-dimensional materials. Using our approach, we recover the transport properties of WTe2 in the para-hydrodynamic regime in good agreement with existing experiments.
Similar content being viewed by others
Introduction
The viscous flow of an interacting electron fluid was predicted a long time ago1; however, experimental evidence for it has remained scarce for many decades2. The advent of ultraclean quantum materials with low carrier density3,4,5 yielded a growing number of cases demonstrating viscous electron flow in the last few years6,7,8,9,10,11. Moreover, it has even become possible to establish ballistic and viscous flow profiles using spatially resolved techniques12,13,14,15,16. Among the coveted properties of hydrodynamic flow is, for example, a negative nonlocal resistance17,18,19,20 as well as other nonlocal transport signatures21,22,23,24,25,26,27. However, the observation of vortical flow (electron whirlpools) still remained elusive28,29,30,31. This situation was upended very recently when a high-fidelity, spatially resolved experiment32 in ultraclean WTe233,34 could unambiguously demonstrate hydrodynamic vortical flow, clearly excluding a ballistic origin of the observed whirlpool pattern.
The observation of this hydrodynamical vortical flow is remarkable for two additional reasons. First, the device is not an effective two-dimensional system but a thin (thickness dâ=â48ânm, width wâ=â550ânm), three-dimensional flake exfoliated and fabricated from WTe2 flakes with a very large bulk mean free path ââââ20âμm. Second, at the measured temperature Tâ=â4.5âK, the electronâelectron interaction leads to a momentum-conserving mean free path âeeâââ10âμmââ«âw. This means that the sample cannot be in the hydrodynamic regime, which is characterized by the condition that âeeââªâwââªââ1.
Based on the unusual properties of WTe2, the authors in ref. 32 suggested that an effectively hydrodynamic flow could instead be induced by almost specular (and thus predominantly momentum-conserving) scattering from the top and bottom surfaces of the three-dimensional WTe2 sheet, a mechanism termed para-hydrodynamics (cf. Fig. 1a). However, no kinetic theory has been put forward which would explain how such a novel type of hydrodynamic flow emerges microscopically.
a The current flows through a thin slab of thickness dââªâwââªââ. Electrons scatter at the microscopically rough top and bottom surfaces with incident angle θ. Most trajectories are reflected almost specularly, leading to angular diffusion in θ, while few trajectories scatter randomly, thus dissipating momentum. bâd Comparison of the non-equilibrium distribution functions f(θ) in different flow regimes (red curve). An approximate cosine that inscribes f(θ) for steep angles θââ â0,âÏ is shown as a blue shaded area, whereas the red shaded area indicates grazing trajectories (θâââ0,âÏ). Different parts of f(θ) have dissimilar scale dependencies in terms of the ratio αâ=âd/2âââªâ1. For ballistic flow, f(θ) decays strongly upon approaching θâ=âÏ/2, while for hydrodynamic flow, it assumes a smooth cosine form everywhere. The intermediate para-hydrodynamic regime carries signatures of both ballistic and hydrodynamic flow, but the smooth part of f(θ) is logarithmically larger than the ballistic part, thus leading to an effectively hydrodynamic current in the limit αâââ0.
In this letter, we consider a kinetic theory for the in-plane flow of an electron fluid in a thin, three-dimensional slab that takes into account weak boundary scattering from the rough top and bottom surfaces. We demonstrate that this setting naturally leads to para-hydrodynamic flow in thin three-dimensional devices as long as the bulk mean free path â of the material is very large compared to the device thickness d. The proposed scattering model is generic, with the only parameters being the amplitude and correlation length of the surface roughness.
The mechanism is depicted schematically in Fig. 1bâd. The emergence of the para-hydrodynamic regime is a result of the conversion of rarely colliding trajectories that impact the top and bottom surfaces at a grazing angle (i.e., ballistic flow) into trajectories that scatter often and at a steep angle from the surface (hydrodynamic flow). This conversion happens due to the slow angular diffusion of the scattered trajectories. The redistribution of statistical weight means that the distribution function and, thus, also the current predominantly resembles hydrodynamic transport. Therefore, while it is not possible to describe the three-dimensional distribution function using a StokesâOhm hydrodynamic approach, the in-plane components of the flow velocity exhibit relaxation properties which resemble viscous flow. Specifically, we find that the effective in-plane momentum-relaxing (âmr) and momentum-conserving (âmc) mean free paths are related as \({\ell }_{mr}/{\ell }_{mc}\propto \log (\ell /d)\, > \,1\). Since the conversion of forward trajectories into steep trajectories is bound to happen whenever angular diffusion is present in the boundary scattering, the only reason why the phenomenon of para-hydrodynamics has so far remained elusive seems to be the logarithmically slow enhancement of this effect with increasing mean free path. In particular, our model posits that neither the precise surface roughness nor the material itself sensitively affects whether a para-hydrodynamic transport regime can emerge in a given material. Instead, the para-hydrodynamic regime merely requires a very large fineness ratio â/d, combined with an appropriate choice of the width w of the slab so that âmrââ«âwââ«ââmc. These conditions were fulfilled in the experiment of ref. 32, which was done with high-quality samples where â/dâ>â500.
Results
Absence of eâe interactions at low T
It has been argued16,32 that WTe2 cannot be in the hydrodynamic regime below 20âK. However, estimates for the effective electronâelectron mean free path âee may vary considerably depending on the employed band structure model and other details of the calculational approach. Therefore, before invoking a surface mechanism, we briefly comment that bulk scattering is unequivocally too weak to matter in these mesoscopic devices. To this end, consider the electronic self-energy for three qualitatively different but realistic candidate band structures of the three-dimensional phase of WTe2. The proposed Fermi surfaces describe (i) a Weyl semimetal phase with very small pockets, (iii) a phase with highly anisotropic pockets, and (iii) a relaxed phase with large pockets (Fig. 2). Employing standard methods to calculate the imaginary part of the self-energy due to the screened Coulomb interaction35 and using a fine momentum grid (cf. Supplementary Information), we obtain the estimates \({\ell }_{ee}^{(i)}=19{{{{{{{\rm{\mu m}}}}}}}}\), \({\ell }_{ee}^{(ii)}=22{{{{{{{\rm{\mu m}}}}}}}}\) and \({\ell }_{ee}^{(iii)}=155{{{{{{{\rm{\mu m}}}}}}}}\) at Tâ=â4.5âK. Therefore, we conclude that details of the calculation, and differing starting assumptions about the qualitative shape of the Fermi pockets do not substantially affect the electronâelectron scattering rate, soundly excluding bulk mechanisms as the source of the hydrodynamic flow. Note that ref. 16 additionally considered phonon-assisted electronâelectron interactions, but they are likewise too weak at low T.
Boundary scattering model
Due to the nature of the para-hydrodynamic transport regime, the derivation cannot rely on a hydrodynamic treatment but has to start from a Boltzmann transport approach. In the following, we consider the nearly ballistic flow in a three-dimensional rectangular geometry where dââªâwââªââ (cf. Fig. 1). In the kinetic approach, the distribution function is denoted by f(r,âk) at real-space position r and for momentum k on a spherical Fermi surface. While the material in consideration has a complicated, non-spherical Fermi surface, the relevant Fermi surface quantity in the Boltzmann equation is the Fermi velocity, which is indeed relatively isotropic. The Boltzmann equation in the steady state is
Here, vF is the Fermi velocity, e is the electron charge, Eâ=â(Ex,â0,â0) is the electric field and \({{{{{{{{\mathcal{I}}}}}}}}}_{0}(f)=|{{{{{{{\bf{v}}}}}}}}|(f-{f}_{0})/\ell\) is the bulk collision integral in relaxation time approximation. \({{{{{{{{\mathcal{I}}}}}}}}}_{b}(f)\) is the collision integral due to boundary scattering from the top and bottom surfaces. Because the thickness d of the flake is much smaller than its width w, we can neglect the spatial dependence along the width of the channel. We note that this approximation becomes exact when the in-plane boundaries at yâ=â±w/2 are completely specular, in which case the distribution function is y-independent. We parametrize \(f({{{{{{{\bf{r}}}}}}}},{{{{{{{\bf{k}}}}}}}})-{f}_{0}=Ah(z,\theta )\cos \theta \cos \phi\), choosingââd/2ââ¤âzââ¤âd/2 along the third dimension,ââÏ/2ââ¤âÏââ¤âÏ/2 as the angle in the plane, andââÏâ<âθââ¤âÏ as the out-of-plane angle, with θâ=âÏâ=â0 pointing alongâ+x (cf. Fig. 1). The dimensionless coefficient Aâ=ââ(âϵâf0)eExâ is chosen such that the solution h becomes normalized (h(z,âθ)â=â1) when boundary scattering is absent, associated with the bulk current density \({j}_{0}=\frac{e{v}_{F}}{32\pi }\int{d}^{3}{{{{{{{\bf{k}}}}}}}}A\).
We now construct h(z,âθ), which solves Eq. (1) for boundary scattering from a rough surface. The microscopic process of boundary scattering has been studied for many decades36,37,38,39,40,41,42,43. It is common to rewrite the boundary collision integral in terms of boundary conditions using a specularity parameter Rθ which may depend on the angle of incidence θ. Here, Rθâ=â0 corresponds to completely diffuse scattering, while Rθâ=â1 is fully specular. Using this parameter, one can express the reflected part of the distribution function h(d/2, ââ£Î¸â£) at the top surface in terms of the incident one as h(d/2, ââ£Î¸â£)â=âRθh(d/2,ââ£Î¸â£), and vice versa at the bottom surface it is h(âd/2, â£Î¸â£)â=âRθh(âd/2, ââ£Î¸â£). However, using a true 2D scattering cross section to describe a boundary scattering event, one has to go back to the full distribution function, expressing the reflected fâ> in terms of the incident fâ< by an integral condition38
where W(k) is the correlation function of the surface scattering potential and the integral runs over the half-sphere (FSâ) of the Fermi surface, which corresponds to trajectories incident to the boundary. We henceforth employ a generic Gaussian-correlated scattering potential38, defined as \(W({{{{{{{\bf{k}}}}}}}})=\pi {a}^{2}{b}^{2}{e}^{-{{{{{{{{\bf{k}}}}}}}}}^{2}{b}^{2}/4}\), with potential depth a and correlation length b.
Historically, Eq. (2) has been solved in two limiting cases, for grazing angles θâââ0 and for steep angles of incidence, θâââÏ/238. For grazing angles, the distribution function changes more rapidly than the scattering cross section W and one obtains a standard form of the boundary in terms of the angle-dependent specularity (1âââqâ£Î¸â£) with parameter \(q=4\sqrt{\pi }{{\Gamma }}(\frac{3}{4}){a}^{2}{k}_{F}^{3/2}/\sqrt{b}\). Conversely, for steep angles, the scattering potential W changes faster than the distribution function. Using a saddle-point approximation therefore yields a FokkerâPlanck equation in the angle of the form \(h(d/2,-|\theta|)={\hat{O}}_{\theta }h(d/2,|\theta|)\) with the operator of the angular diffusion on the Fermi surface being \({\hat{O}}_{\theta }=Q{\sin }^{2}\theta ((2\cot \theta -\tan \theta ){\partial }_{\theta }+{\partial }_{\theta }^{2})\), where Qâ=â8a2/b2. To our knowledge, no attempts have been made to treat the scattering for both limits, θâââ0 and θâââÏ/2, in a unified framework. However, to capture the para-hydrodynamic behavior, we seek a solvable description that holds for all angles of incidence. Indeed, as we demonstrate next, the description for general angle θ is absolutely vital to capture the parametric dependencies of the para-hydrodynamic flow correctly.
In the spirit of Matthiessenâs rule, we propose to add both scattering limits as two different types of scattering processes, which yields for the distribution function at the upper boundary of the slab the boundary condition
Here, we introduced the specularity parameter \({R}_{\theta }=(1-q\sin|\theta|)\), which is the periodic extension of the previously mentioned specularity at small angles. Note that both scattering types can be combined safely because for θâ=â0,âÏ, the momentum-relaxing scattering vanishes (R0â=â1), meaning that the distribution function has no discontinuities anywhere. Equation (3) presents the key innovation for a unified treatment of boundary scattering.
The ordinary differential Eq. (3) exhibits several favorable properties. First, solutions for qâ=â0 and for any Qâ>â0 are non-dissipative, with distribution function h(z,âθ)ââ¡â1, corresponding to a bulk current profile without any stresses. This can be understood as follows. Q parameterizes the relative importance of angular diffusion due to scattering from the boundaries. However, without any momentum loss (i.e., unless qâ>â0), the angular diffusion of momentum will redistribute momenta equally into higher and lower momentum states, thereby conserving momentum exactly. Second, the differential equation becomes stiff both at θâ=â0 and θâ=âÏ/2, both of which constitute singular points. Therefore, the resulting distribution function at either point singularly depends on the initial conditions at the respective other point. Indeed, we find that for Qââ«âq, the solution retains a singular dependence on q for all angles, even though Rθ is dominant over \({\hat{O}}_{\theta }\) only for shallow angles smaller than d/â.
Intermediate regime of para-hydrodynamics
Using a symmetric parametrization in terms of the characteristics of the motion44, the solution of the Boltzmann equation can be written as
c(θ) is symmetric in θ and is determined by the modified boundary scattering condition. Inserting the general solution of Eq. (4) into (3), one obtains for 0â<âθâ<âÏ/2
where αâ=âd/2â and \({h}_{in}^{\pm }=(1-{e}^{\mp \alpha \csc \theta })\). Solving Eq. (5) is rather involved (cf. Supplementary Information). Typical solutions for qâ=â1 and several values of Q are shown in Fig. 3a. For Qâââ0, the distribution function quickly decays with increasing angle, meaning that the current originates almost exclusively from long-lived trajectories in the forward direction (θâ=â0). This is the expected behavior for ballistic flow. In contrast, for Qâ>â1, the profile changes qualitatively, with the trajectories at grazing angles being strongly suppressed, while the distribution function becomes a cosine (corresponding to \(c(\theta )={{{{{{{\rm{const.}}}}}}}}\)) for all other angles. This shape of the distribution function resembles a hydrodynamic distribution function. The flat part is characterized by the asymptotic value c1â=âc(Ï/2), which by Taylor expansion in θâââÏ/2, evaluates to
where \({c}^{\prime{\prime}}(\pi /2)=\frac{1}{2}{\partial }_{\theta }^{2}c(\theta ){|}_{\theta=\pi /2}\). In the ballistic regime (Qâ=â0), we find numerically that câ³(Ï/2)âââα, and thus also c1âââα, which is subleading compared to the forward trajectories at θâââ0, which contribute to the current at order \({{{{{{{\mathcal{O}}}}}}}}(\sqrt{\alpha })\)38. In contrast, for Qâ>â0, we find that \({c}^{\prime{\prime}}(\pi /2)\propto \alpha \log {\alpha }^{-1}\). The weight of trajectories with steep angles is therefore substantially increased compared to the forward trajectories, which are in turn suppressed and contribute to the current only at order \({{{{{{{\mathcal{O}}}}}}}}(\alpha )\). For small values of α, the enhancement of the flat part of the distribution function compared to the forward (ballistic) trajectories is therefore large enough so that a para-hydrodynamic regime emerges (cf. Fig. 1).
a The distribution function c(θ) in the para-hydrodynamic regime, for qâ=â1, three different values of Q and two α. With increasing Q, more weight is accumulated at the steep trajectories away from θâ=â0. For αâ=â0.05, Qâ=â4, the asymptotic value c(Ï/2)â=âc1 is indicated by a gray dashed line. In the two-fluid approximation, the current jh stemming from the area below c1 (gray, hatched) is compared against the current jb induced by the rest of the non-equilibrium distribution. b Comparison of the para-hydrodynamic (jh/j0) and ballistic (jb/j0) current density contributions in the two-fluid approximation for two values of α, where j0 denotes the bulk current density. For all αââªâ1, the contribution of the para-hydrodynamic current increases with Q, quickly overtaking the ballistic contribution.
Since the logarithmic enhancement makes it impossible to construct the size of c1 from local properties around θâ=âÏ/2, we now focus on the limit where αâââ0. To leading order in α, in terms of the variable \(s=\sin \theta\), Eq. (5) can be expanded as
where \({s}_{\perp }=\sqrt{1-{s}^{2}}\). This equation no longer contains any essential singularities for shallow θ and is much easier to handle numerically. Based on the analytical structure of câ³(Ï/2) in terms of higher derivatives and using the numerical solution as a reference, we find by fitting that \({c}^{\prime{\prime}}(\pi /2)\,\approx \,\alpha \log (0.06/\alpha )/(0.4q+0.5Q)\).
Effective mean free paths
The crossover from ballistic to bulk hydrodynamic flow is well studied in two-dimensional electron fluids2,22,42,44,45,46,47,48,49 and is typically described using a dual relaxation time approximation for the scattering integral in terms of the two length scales âmr and âmc. However, it has been shown44,50 that the nature of the transport regime can also be reconstructed from an inspection of the distribution function. Namely, the distribution is smooth and relatively flat in the hydrodynamic regime (âmcââªâwââªââmr), while in the ballistic regime (wââªââmr,ââmc), scattering from the boundaries makes the distribution function strongly angle dependent. This can be more formally restated by considering the angular harmonics of h(θ). Keeping only the first and second terms in such an expansion, the Boltzmann equation simplifies to the StokesâOhm equation44, i.e., the distribution function can be obtained in a hydrodynamic description. On the other hand, if higher angular harmonics are present in the distribution function, this indicates the presence of additional long-lived modes in the flow, as would be expected for a ballistic distribution function.
In the present case, different parts of the distribution function at respectively shallow or steep angles resemble either the ballistic or the hydrodynamic situation. We therefore propose a two-fluid approximation (cf. Figs. 1 and 3), whereby we decompose the distribution function h(θ) into a constant part which constitutes a hydrodynamic current density jhâ=âc1j0, while the remaining strongly angle-dependent parts constitute a ballistic current density jbâ=âjâââjh.
As shown in Fig. 3b, at Qâ=â0, it is jhâ<âjb, but upon increasing Q, there is a crossover into a regime with jhâ>âjb. In the latter regime, one can immediately infer âmr and âmc from the flat part of the distribution function, which creates the dominant contribution jh to the current density. To this end, using the reduced current density, we write for the effective momentum-relaxing mean free path,
The flatness of the distribution function furthermore suggests that the total scattering rate is large and completely dominated by momentum-conserving processes. âmc is therefore expected to approach its maximally possible value, which for a flat distribution is entirely determined geometrically by the normalized travel distance between two successive scatterings from the top and bottom surfaces. In other words, we can estimate that \(\frac{1}{{\ell }_{mc}}\approx \frac{1}{d}\frac{1}{2\pi }d\theta \int\nolimits_{0}^{\pi }\sin|\theta|\), which yields âmcâ=âÏd.
Using these values for âmrââªââ and âmcââªââ, the effective Gurzhi parameter for the para-hydrodynamic flow in thin sheets becomes
In the experiment, the WTe2 samples had the properties dâ=â48ânm, ââ=â20âμm, âmrâ=â530ânm, and Dâ=â155ânm. Using as input the experimental values for d and â, Eq. (10) yields \({D}^{{\prime} }=145{{{{{{{\rm{nm}}}}}}}}\) and âmrâ=â560ânm for the choice (q,âQ)â=â(1,â4), while it is \({D}^{{\prime} }=143{{{{{{{\rm{nm}}}}}}}}\) and âmrâ=â545ânm upon choosing (q,âQ)â=â(0.8,â2). This indicates that the kinetic theory is not strongly sensitive to the precise value of q and Q, and it can explain the experimental findings well for a reasonable range of boundary scattering parameters. Importantly, this versatility implies that para-hydrodynamic flow emerges generically as long as there is a large-scale separation between sheet thickness d and bulk mean free path â and does not require extensive fine-tuning. We remark that q and Q are, in turn, only weakly dependent on the microscopic scattering parameters a and b. For example, using kFâ=â1ânmâ1, we obtain aâ=â0.27ânm, bâ=â0.37ânm for (q,âQ)â=â(1,â4), and aâ=â0.26ânm, bâ=â0.51ânm for (q,âQ)â=â(0.8,â2). These estimates correspond to boundary roughness at the lattice scale and are indeed reasonable for cleaved samples34.
Discussion
We derived a microscopic boundary scattering model which can explain para-hydrodynamic flow in the absence of strong electronâelectron scattering. As the main characteristics of the new regime, we identified singular points in the resulting boundary condition, which lead to a different thickness dependence of the current density that onsets in the presence of small-angle scattering from the boundary. Our findings constitute a new type of ballistic-to-hydrodynamic crossover, where the three-dimensional problem is microscopically ballistic, but the in-plane components of the flow velocity exhibit relaxation properties which are indistinguishable from viscous flows. Our results indicate that the phenomenology which was previously suggested to govern the hydrodynamic-to-ballistic crossover is not universal. Since para-hydrodynamic flow exclusively emerges in the presence of angular diffusion from short-range correlated disorder, conventional approaches which rely solely on a reflectivity coefficient have not been able to capture this mechanism2,51,52. We note that small-angle scattering can also appear from bulk scattering, in which case it typically leads to a ratio âmr/âmcâ=â429,53. For small-angle boundary scattering, we instead found âmrâ=âc1â and âmcâ=âÏd, which yields a scale-dependent ratio \({\ell }_{mr}/{\ell }_{mc}\propto \log (\ell /d)\) that becomes large enough to support hydrodynamic phenomena for very large ratios â/d.
It would be interesting to find additional signatures of the para-hydrodynamic regime for the channel flow. Since the in-plane current density in narrow channels is not suitable for distinguishing between the ballistic and hydrodynamic transport regimes14, this would probably involve the investigation of the Hall viscosity at finite magnetic fields46,54 or optical probes.
Since only very few microscopic parameters other than the bulk mean free path enter into our results, we expect that the para-hydrodynamic flow observed in WTe2 is not unique and that a number of clean, three-dimensional compounds, for example, Weyl and Dirac semimetals4,5 and ultrapure delafossites7,55, can exhibit this new transport regime. The proposed boundary scattering model and the two-fluid approximation used here for extracting the effective mean free paths are generic. Thus, the same methodology can very likely be applied directly to many other three-dimensional materials and also integrated into numerical schemes56.
Methods
The calculation of the electronâelectron scattering rate was done in four steps. The first step in the process is finding the energy bands and wave functions of the Weyl semimetal phase of WTe2 within DFT using VASP and Wannier9057,58. These were evaluated inside the Brillouin zone on a k-mesh of NxâÃâNyâÃâNzâ=â100âÃâ50âÃâ7. The particle-hole bubble Î was calculated with an IR cutoff of 5âmeV, which is approximately a tenth of the Fermi energy. The real and imaginary parts were calculated separately to avoid numerical errors. \({{{{{{{\rm{Re}}}}}}}}\,{{\Pi }}(q,\omega )\) and \({{{{{{{\rm{Im}}}}}}}}\,{{\Pi }}(q,\omega )\) were evaluated on the same k-mesh as the bands and wave functions and on energies between â11 and 11âeV with a resolution of 0.055âeV. On energies in between the grid points, a linear interpolation was used. The UV energy cutoff is chosen such that the low-temperature self-energy corrections are converged, where the main contribution to the relaxation time is by electrons with energy ϵnq in band n residing in an energy window of ϵFâââTâ<âϵnqâ<âϵFâ+âT that are scattered by electrons in band m that likewise satisfy ϵFâââTâ<âϵmkâ<âϵFâ+âT. Thus, the main contribution to the self-energy is generated by electrons whose differences in energies are Ïâ=âϵnqâââϵmkâ<â2T, while the target temperatures are below Tâ<â0.025âeVââªâ11âeV. It is then sufficient to let the band indices run through the first 10 bands above and below the Fermi surface, which corresponds to a maximal energy difference of 4âeV.
The self-energy was calculated with a regularized Bose distribution function, \(b(\epsilon )\to ({e}^{\beta \epsilon }-1)/({({e}^{\beta \epsilon }-1)}^{2}+\lambda )\) (cf. Supplementary Information). A good choice for the regulator is λâ=â10â4, such that the pole at zero frequency is cut off at b(ϵ)â<â104.
We confirmed the convergence of our results by comparing various grid resolutions. First, two k-meshes were chosen at sizes 200âÃâ50âÃâ3 and 100âÃâ100âÃâ3. The difference in the self-energy between both grid resolutions is less than 2%. Additionally, the frequency resolution was also tested by increasing the resolution to 0.02âeV, which yielded a change of less than 1%. Finally, the bosonic occupation number cutoff was increased to λâ=â10â6, which yielded a change of less than 1%. The sensitivity of our results against small changes in the Fermi level was also tested. On two points separated by 0.02âeV, the results changed by about 36%, which is significant but shows that the order of magnitude is robust to small changes in the chemical potential. A calculation with a different number of bands included in the correction of the photon propagator showed no relevant changes.
Data availability
The data that support the findings of this study are available from the corresponding author upon request.
References
Gurzhi, R. N. Reviews of topical problems: hydrodynamic effects in solids at low temperature. Sov. Phys. Usp. 11, 255 (1968).
de Jong, M. J. M. & Molenkamp, L. W. Hydrodynamic electron flow in high-mobility wires. Phys. Rev. B 51, 13389 (1995).
Narozhny, B. N. Electronic hydrodynamics in graphene. Ann. Phys. 411, 167979 (2019).
Yan, B. & Felser, C. Topological materials: Weyl semimetals. Annu. Rev. Condens. Matter Phys. 8, 337 (2017).
Armitage, N. P., Mele, E. J. & Vishwanath, A. Weyl and Dirac semimetals in three-dimensional solids. Rev. Mod. Phys. 90, 015001 (2018).
Krishna Kumar, R. et al. Superballistic flow of viscous electron fluid through graphene constrictions. Nat. Phys. 13, 1182 (2017).
Moll, P. J. W., Kushwaha, P., Nandi, N., Schmidt, B. & Mackenzie, A. P. Evidence for hydrodynamic electron flow in PdCoO2. Science 351, 1061 (2016).
Gooth, J. et al. Thermal and electrical signatures of a hydrodynamic electron fluid in tungsten diphosphide. Nat. Commun. 9, 4093 (2018).
Berdyugin, A. I. et al. Measuring Hall viscosity of grapheneâs electron fluid. Science 364, 162 (2019).
Osterhoudt, G. B. et al. Evidence for dominant phonon-electron scattering in Weyl semimetal WP2. Phys. Rev. X 11, 011017 (2021).
Kumar, C. et al. Imaging hydrodynamic electrons flowing without Landauer-Sharvin resistance. Nature 609, 276 (2022).
Bandurin, D. A. et al. Negative local resistance caused by viscous electron backflow in graphene. Science 351, 1055 (2016).
Marguerite, A. et al. Imaging work and dissipation in the quantum Hall state in graphene. Nature 575, 628 (2019).
Sulpizio, J. A. et al. Visualizing Poiseuille flow of hydrodynamic electrons. Nature 576, 75 (2019).
Ku, M. J. H. et al. Imaging viscous flow of the Dirac fluid in graphene. Nature 583, 537 (2020).
Vool, U. et al. Imaging phonon-mediated hydrodynamic flow in WTe2. Nat. Phys. 17, 1216 (2021).
Levin, A. D., Gusev, G. M., Levinson, E. V., Kvon, Z. D. & Bakarov, A. K. Vorticity-induced negative nonlocal resistance in a viscous two-dimensional electron system. Phys. Rev. B 97, 245308 (2018).
Bandurin, D. A. et al. Fluidity onset in graphene. Nat. Commun. 9, 4533 (2018).
Kim, M. et al. Control of electron-electron interaction in graphene by proximity screenings. Nat. Commun. 11, 2339 (2020).
Gupta, A. et al. Hydrodynamic and ballistic transport over large length scales in GaAs /AlGaAs. Phys. Rev. Lett. 126, 076803 (2021).
Guo, H., Ilseven, E., Falkovich, G. & Levitov, L. S. Higher-than-ballistic conduction of viscous electron flows. Proc. Natl Acad. Sci. USA 114, 3068 (2017).
Alekseev, P. S. Negative magnetoresistance in viscous flow of two-dimensional electrons. Phys. Rev. Lett. 117, 166601 (2016).
Torre, I., Tomadin, A., Geim, A. K. & Polini, M. Nonlocal transport and the hydrodynamic shear viscosity in graphene. Phys. Rev. B 92, 165433 (2015).
Mendoza, M., Herrmann, H. J. & Succi, S. Preturbulent regimes in graphene flow. Phys. Rev. Lett. 106, 156601 (2011).
Stern, A. et al. How electron hydrodynamics can eliminate the Landauer-Sharvin resistance. Phys. Rev. Lett. 129, 157701 (2022).
Di Sante, D. et al. Turbulent hydrodynamics in strongly correlated Kagome metals. Nat. Commun. 11, 3997 (2020).
Raichev, O. E., Gusev, G. M., Levin, A. D. & Bakarov, A. K. Manifestations of classical size effect and electronic viscosity in the magnetoresistance of narrow two-dimensional conductors: theory and experiment. Phys. Rev. B 101, 235314 (2020).
Falkovich, G. & Levitov, L. Linking spatial distributions of potential and current in viscous electronics. Phys. Rev. Lett. 119, 066601 (2017).
Nazaryan, K. G. & Levitov, L. Robustness of vorticity in electron fluids. Preprint at https://arxiv.org/abs/2111.09878 (2021).
Pellegrino, F. M. D., Torre, I., Geim, A. K. & Polini, M. Electron hydrodynamics dilemma: whirlpools or no whirlpools. Phys. Rev. B 94, 155414 (2016).
Alekseev, P. S. et al. Counterflows in viscous electron-hole fluid. Phys. Rev. B 98, 125111 (2018).
Aharon-Steinberg, A. et al. Direct observation of vortices in an electron fluid. Nature 607, 74â80 (2022).
Wang, L. et al. Tuning magnetotransport in a compensated semimetal at the atomic scale. Nat. Commun. 6, 8892 (2015).
Woods, J. M. et al. Suppression of magnetoresistance in thin WTe2 flakes by surface oxidation. ACS Appl. Mater. Interfaces 9, 23175â23180 (2017).
DellâAnna, L. & Metzner, W. Fermi surface fluctuations and single electron excitations near Pomeranchuk instability in two dimensions. Phys. Rev. B 73, 045127 (2006).
Soffer, S. B. Statistical model for the size effect in electrical conduction. J. Appl. Phys. 38, 1710 (1967).
Okulov, V. I. & Ustinov, V. V. Boundary condition for the distribution function of conduction electrons scattered by a metal surface. J. Exp. Theor. Phys. 40, 584 (1975).
Falkovsky, L. A. Transport phenomena at metal surfaces. Adv. Phys. 32, 753 (1983).
Moraga, L. A. Quantum reflectivity of random rough surfaces and the electrical resistivity of thin films and wires. J. Phys. F Metal Phys. 17, 1751 (1987).
Einzel, D., Panzer, P. & Liu, M. Boundary condition for fluid flow: curved or rough surfaces. Phys. Rev. Lett. 64, 2269 (1990).
Akhmerov, A. R. & Beenakker, C. W. J. Boundary conditions for Dirac fermions on a terminated honeycomb lattice. Phys. Rev. B 77, 085423 (2008).
Kiselev, E. I. & Schmalian, J. Boundary conditions of viscous electron flow. Phys. Rev. B 99, 035430 (2019).
Raichev, O. E. Linking boundary conditions for kinetic and hydrodynamic description of fermion gas. Phys. Rev. B 105, L041301 (2022).
Holder, T. et al. Ballistic and hydrodynamic magnetotransport in narrow channels. Phys. Rev. B 100, 245305 (2019).
Shytov, A., Kong, J. F., Falkovich, G. & Levitov, L. Particle collisions and negative nonlocal response of ballistic electrons. Phys. Rev. Lett. 121, 176805 (2018).
Scaffidi, T., Nandi, N., Schmidt, B., Mackenzie, A. P. & Moore, J. E. Hydrodynamic electron flow and Hall viscosity. Phys. Rev. Lett. 118, 226601 (2017).
Ledwith, P., Guo, H., Shytov, A. & Levitov, L. Tomographic dynamics and scale-dependent viscosity in 2D electron systems. Phys. Rev. Lett. 123, 116601 (2019).
Afanasiev, A. N., Alekseev, P. S., Greshnov, A. A. & Semina, M. A. Ballistic-hydrodynamic phase transition in flow of two-dimensional electrons. Phys. Rev. B 104, 195415 (2021).
Li, S., Khodas, M. & Levchenko, A. Conformal maps of viscous electron flow in the Gurzhi crossover. Phys. Rev. B 104, 155305 (2021).
Holder, T., Queiroz, R. & Stern, A. Unified description of the classical Hall viscosity. Phys. Rev. Lett. 123, 106801 (2019).
Alekseev, P. S. & Semina, M. A. Ballistic flow of two-dimensional interacting electrons. Phys. Rev. B 98, 165412 (2018).
Svintsov, D. Hydrodynamic-to-ballistic crossover in Dirac materials. Phys. Rev. B 97, 121405 (2018).
Ledwith, P. J., Guo, H. & Levitov, L. The hierarchy of excitation lifetimes in two-dimensional Fermi gases. Ann. Phys. 411, 167913 (2019).
Delacrétaz, L. V. & Gromov, A. Transport signatures of the Hall viscosity. Phys. Rev. Lett. 119, 226602 (2017).
Mackenzie, A. P. The properties of ultrapure delafossite metals. Rep. Prog. Phys. 80, 032501 (2017).
Varnavides, G., Jermyn, A. S., Anikeeva, P. & Narang, P. Probing carrier interactions using electron hydrodynamics. Preprint at https://arxiv.org/abs/2204.06004 (2022).
W. Kohn, A. D. B. & Parr, R. G. Density functional theory of electronic structure. J. Phys. Chem. 100, 12974â12980 (1996).
Perdew, J. P., Burke, K. & Ernzerhof, M. Generalized gradient approximation made simple. Phys. Rev. Lett. 77, 3865 (1996).
Acknowledgements
We thank J.S. Hofmann and E. Zeldov for fruitful discussions. B.Y. acknowledges the financial support from the European Research Council (ERC Consolidator Grant No. 815869, âNonlinearTopoâ) and Israel Science Foundation (ISF No. 2932/21).
Author information
Authors and Affiliations
Contributions
T.H. conceived the project. T.H. and Y.W. did the analytical and numerical calculations for the boundary model. Y.W. did the DFT and eâe scattering rate calculations. Y.W., A.A.-S., B.Y., and T.H. analyzed the results. T.H. wrote the manuscript with input from all other authors. B.Y. supervised the project.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Peer review
Peer review information
Nature Communications thanks Georgios Varnavides for their contribution to the peer review of this work. Peer reviewer reports are available.
Additional information
Publisherâs note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the articleâs Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the articleâs Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Wolf, Y., Aharon-Steinberg, A., Yan, B. et al. Para-hydrodynamics from weak surface scattering in ultraclean thin flakes. Nat Commun 14, 2334 (2023). https://doi.org/10.1038/s41467-023-37966-z
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41467-023-37966-z
This article is cited by
-
Alternative routes to electron hydrodynamics
Communications Physics (2024)
-
Charge transport and hydrodynamics in materials
Nature Reviews Materials (2023)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.