Abstract
A current challenge in condensed matter physics is the realization of strongly correlated, viscous electron fluids. These fluids can be described by holography, that is, by mapping them onto a weakly curved gravitational theory via gauge/gravity duality. The canonical system considered for realizations has been graphene. In this work, we show that Kagome systems with electron fillings adjusted to the Dirac nodes provide a much more compelling platform for realizations of viscous electron fluids, including non-linear effects such as turbulence. In particular, we find that in Scandium Herbertsmithite, the fine-structure constant, which measures the effective Coulomb interaction, is enhanced by a factor of about 3.2 as compared to graphene. We employ holography to estimate the ratio of the shear viscosity over the entropy density in Sc-Herbertsmithite, and find it about three times smaller than in graphene. These findings put the turbulent flow regime described by holography within the reach of experiments.
Similar content being viewed by others
Introduction
Electrons in solids typically interact not only with impurities and phonons, but also with each other via the Coulomb interaction. If the momentum relaxing effects of impurities and phonons are weak, the Coulomb interaction can become dominant and lead to local thermalization and the formation of an electronic fluid1,2. Thus, the length and time scales over which thermalization occurs are controlled by the strength of the Coulomb interaction. The regime of electron hydrodynamics has been realized in several systems3,4,5,6,7, giving rise to new transport properties8,9,10,11,12 distinctly different from the ballistic regime. Therefore, characterizing the transport properties of viscous electronic fluids is tantamount to determining the strength of the Coulomb coupling α.
The Coulomb interaction also mediates energy and momentum transfer within the fluid ensuring a local thermal equilibrium, and thereby controls transport coefficients such as the shear viscosity η. Since direct access to the Coulomb coupling α proves difficult, we focus in this work on the easily accessible shear viscosity, or more precisely on the ratio between η and the entropy density of the fluid, s. As explained in the âMethodsâ section, the shear viscosity η is straightforwardly obtained from a Kubo formula (see Eq. (16) in âMethodsâ section), and the entropy density s from thermodynamics (see Eq. (8) in âMethodsâ section). For a Dirac fluid, the ratio η/s is equal to the temperature multiplied by the kinematic viscosity ν, which is experimentally accessible11,12, η/sâ=âTν. The physical relevance of η/s is that it yields the shear viscosity per effective degree of freedom participating in the momentum diffusion in the fluid. The ratio η/s depends sensitively on α, as shown in Fig. 1. In the weakly Coulomb interacting regime, i.e. for αââªâ1, first order perturbation theory such as Boltzmannâs kinetic theory predicts a fast fall off, η/sâ~âαâ2, as shown by the black line in Fig. 1. For intermediate Coulomb couplings, perturbative approaches lose their validity. Holographic gauge/gravity duality13 provides a nonperturbative approach to predict the coupling dependence of η/s. In the limit of infinitely strong coupling and for systems with a large number of degrees of freedom, it predicts the universal value14,15,16 (see the red dashed line in Fig. 1)
The essential feature of this result is that it is significantly smaller than any value of η/s obtained within weak-coupling perturbation theory. Beyond the infinite coupling limit, gauge/gravity duality allows to include finite-coupling corrections to the infinite coupling result in an expansion in the inverse coupling. As we explain in the âMethodsâ section, the leading-order finite-coupling correction relevant to the materials considered is
where α is the fine-structure constant and the constant \({\mathcal{C}}\) parametrizes the class of gauge/gravity duals considered.
Black line: prediction in the weak coupling regime, \(\left(\eta /s\right)\propto {\alpha }^{-2}\), reliable for small values of α. The red dashed line corresponds to the universal holographic value, \(\eta /s = \, \hbar \,/4\pi {k}_{\mathrm{B}}\). The blue shaded region, for which η/s is given by Eq. (2), represents the holographic prediction. This shaded region parametrizes the class of extrapolating models beyond the α â â limit (see âMethodsâ section). Notice that Sc-Herbertsmithite (Sc-Hb) shows a much smaller variance than graphene (given by the vertical blue bars at α = 2.9 for Sc-Hb and α = 0.9 as a representative value for graphene, respectively), providing further support for the applicability of holographic methods to correlated Kagome metals.
So far, the material of choice to investigate hydrodynamics of electronic systems has been graphene17. Here we show that in certain Kagome materials, hydrodynamic behavior will be significantly enhanced. Specifically, we focus on Kagome materials at filling levels such that the chemical potential is located at the Dirac point. These materials are particularly suited for hydrodynamic studies since not only the Coulomb interaction is enhanced as compared to graphene (see below), but also because the Kagome lattice structure suppresses the formation of ordered phases. The reason is that, in contrast to graphene, the Dirac cones on the Kagome lattice are located far away from half filling. As explained in the âMethodsâ section, combined with the small low energy density of states, this suggests a strong resilience of the metallic Dirac state against ordering instabilities18,19, which implies that it sustains stronger Coulomb coupling than a Dirac metal on the honeycomb lattice20. For the explicit candidate Kagome metal Scandium-substituted Herbertsmithite ScCu3(OH)6Cl2 (Sc-Herbertsmithite hereafter, see Fig. 2)21, we further calculate the phonon spectrum and find that the optical phonons decouple from the electronic degrees of freedom for temperatures below ~80âK, ensuring that a Kagome Dirac metal with electronic interactions is the appropriate microscopic description.
a Top view of the crystal structure, where the CuO4 plaquettes form a Kagome lattice, as highlighted by the orange lines. b Ab-initio band structure of the low-energy manifold of Sc-Herbertsmithite along the high-symmetry directions of the conventional hexagonal Brillouin zone21. The dashed black lines refer to a linear fit around the Dirac point of Sc-Herbertsmithite, whereas the solid gray lines denote the counterpart for graphene. c Phonon dispersion of Sc-Herbertsmithite (gray lines) and relative distribution of the electron-phonon coupling strengths λνq (blue circles, ν is the branch index, q the phonon momentum). The horizontal dashed line marks the temperature Tph above which optical phonon modes with sizeable λνq are thermally activated.
Results
The fine-structure constant
The suppression of gapped ordered phases and the decoupling of phonons in the proposed Kagome materials allows the electrons to form a fluid by Coulomb interactions up to very low temperatures. The electronic Dirac fluid of particle-hole excitations around the Dirac point then has emergent relativistic symmetry, is parity and time reversal invariant, and particle-hole symmetric. With these symmetries in place, a relativistic fluid is described by relativistic hydrodynamic equations of motion, which depends on the following key parameters: The Fermi velocity vF of the relativistic dispersion relation playing the role of the speed of light, the relative dielectric constant ϵr in the medium, and the shear viscosity η of the fluid. As explained in the âMethodsâ section, the two other transport coefficients bulk viscosity ξ and interaction induced conductivity ÏQ are not relevant for our arguments. The quantities vF and ϵr set the value of the fine-structure constant (effective Coulomb coupling) α via
where ϵ0 is the dielectric constant in vacuum. In turn, η depends on α. We calculate vF and ϵr, and hence α, within the framework of the constrained Random Phase Approximation (cRPA) (see the âMethodsâ section), and compare the results for Sc-Herbertsmithite with (hBN encapsulated) graphene. While our findings are likely to be applicable to a broad class of Kagome metals, Sc-Herbertsmithite suitably underlines a prime motif of how to accomplish Dirac fillings in Kagome systems. In pristine Herbertsmithite, Zn2+ acts as a charge donor to the Cu \({d}_{{{\rm{x}}}^{2}-{{\rm{y}}}^{2}}\) orbitals which dominate the low energy theory, yielding a half filled Kagome lattice setting. Synthesis of the otherwise identical compound with Zn2+ replaced by Sc3+ provides the precise stoichiometry for the electrochemical potential to coincide with the Dirac points. While our ab-initio studies can only confirm the chemical stability of the ultimate Sc-Herbertsmithite crystal, the similar atomic radius of Sc in comparison to Zn at least opens the possibility that a related chemical synthesis may be possible. Note that in ref. 21 the similar motif was pursued with Ga-substituted Herbertsmithite. For chemical reaons, we believe Sc to be more suited, because Ga quickly turns fluid and may more easily leak out within a chemical synthesis procedure at higher temperature.
Our first result (c.f. Table 1) is that the fine-structure constant α in Sc-Herbertsmithite is more than three times larger than in graphene, implying a strong enhancement of applicability of viscous hydrodynamics in this material. With Sc-Herbertsmithite being our candidate material for the realization of holographic hydrodynamics, we refer to graphene as a benchmark for our theoretical analysis, as well as a reference point at weaker coupling, which helps us to identify expected trends in observables in more strongly coupled Dirac materials.
Corrections to the infinite coupling result
Our second result is an estimate of the leading finite coupling correction in Eq. (2) within a broad class of holographic models, as explained below and in the âMethodsâ section. This estimate leads to the blue band in Fig. 1, allowing for predictions for the possible range of η/s for both graphene and Sc-Herbertsmithite. Our predictions are indicated by the blue bars and show not only a considerably smaller value η/s for Sc-Herbertsmithite as compared to graphene, but also a smaller variance. In particular, our holographic estimate shows that η/s for the Dirac fluid in Sc-Herbertsmithite lies significantly closer to the infinite coupling value of Eq. (1) predicted by gauge/gravity duality, than in graphene. This makes the Dirac fluid in Sc-Herbertsmithite an interesting candidate for experimental realizations of holographic hydrodynamics.
Robustness of the hydrodynamic regime
We can understand the difference in the scales for Sc-Herbertsmithite as compared to graphene along the following lines. In order for an electron system to display hydrodynamic behavior, the electron-electron mean free path must obey1
where âph is the electron-phonon mean free path, âimp is the electron-impurity mean free path, and W the width of the sample. From Eq. (4), we infer that the emergence of hydrodynamic flow in solids is closely related to the relative value of the characteristic scales of the system3,4. This implies that the formation of viscous flow in solids is restricted to specific, and sometimes narrow, regimes. For Dirac fluids satisfying μââªâkBT, the emergence of hydrodynamic behavior is greatly enhanced2, since in this large temperature regime the available phase space for electron-electron collisions is vastly enhanced. In turn, Landau quasiparticles are short lived, which yields a decrease of âee. In addition, since âeeâââ1/α2, increasing the Coulomb interaction strength leads to shorter thermalization length and time scales. This leads to an enhancement of hydrodynamic behavior and a higher tendency to exhibit viscous flow2. Moreover, the emergence of additional non-hydrodynamic modes sets a critical length scale âc, below which the standard hydrodynamic approach ceases to be valid. For the particular case of Sc-Herbertsmithite, we find (see âMethodsâ section)
This shows that the hydrodynamic regime is more robust in Sc-Herbertsmithite than in graphene.
Turbulent hydrodynamics
Note that the shear viscosity controls the turbulent behavior of the fluid22,23, an aspect of electron hydrodynamics that still remains largely unexplored. Within the hydrodynamic regime, the relative size of dissipative forces to nonlinear inertial forces in the fluid distinguishes the laminar flow regime, in which viscous forces dominate, from the turbulent flow regime, where non-linear dynamics in the Navier-Stokes equations dominate. More specifically, the ratio of inertial to dissipative/viscous forces is characterized by the Reynolds number of the flow. For a 2â+â1 dimensional Dirac fluid at charge neutrality flowing within a channel of width W, the Reynolds number is given by
where s the entropy density of the fluid and utyp is the typical velocity of the fluid, which strongly increases as η/s decreases24. The fluid exhibits turbulent behavior when \(Re={\mathcal{O}}(1000)\)25.
The onset of turbulence is related to the fluid flow, which is sensitive not only to the sample geometry, but also to the intrinsic transport properties such as its shear viscosityâfrom Eq. (6), Reâââ(η/s)â1utyp(η/s). Since the shear viscosity is much smaller for strongly coupled fluids (see Fig. 1), the experimental realization of turbulence will be greatly enhanced for materials exhibiting large Coulomb interactions. This conclusion holds regardless of the particular geometry employed. Furthermore, even laminar flows such as the Poiseuille channel flow will, compared to graphene or other even more weakly coupled materials, exhibit enhanced hydrodynamic transport properties such as larger typical fluid velocities and smaller differential resistances24.
The third and final result of this work is an estimate of the Reynolds number for the flow of our proposed holographic Dirac fluid in Sc-Herbertsmithite through a long and straight channel flow as described by Eq. (6). In graphene, the Reynolds number has been found to be sufficiently large for preturbulent behavior such as vortex production, but not for fully developed turbulence26. For Sc-Herbertsmithite, we find an enhancement of the Reynolds number compared to graphene by a factor of order 100, mostly due to the Fermi velocity νF that is six times smaller (see Fig. 2b and Table 1, as well as âMethodsâ for the precise estimate of the enhancement factor), leading to Reynolds numbers sufficiently large for fully turbulent channel flows.
Discussion
We hence predict that Sc-Herbertsmithite, and correlated Kagome metals at Dirac filling in general, bring the turbulent flow regime in two-dimensional electron hydrodynamics within experimental reach. As alternative routes, turbulent flow is predicted for three-dimensional hydrodynamic Weyl semimetals, from anisotropic two-dimensional Dirac fluids setting in at a possibly small Reynolds number27,28, or particular superconductors in a fluctuation regime around Tc29. Interestingly, yet another promising pursuit of a turbulent electronic regime is to consider electronic fluids with particularly low entropy. This might be the case for the strange metal regime in cuprate superconductors30.
We note that our theoretical results apply more generally than to the specific material candidate Sc-Herbertsmithite, since they provide a new viewpoint on viscous electron fluids in strongly correlated materials. Further examples they may apply to include copper-based metal-organic frameworks (MOFs) in which electrons live on a Kagome lattice. These were recently synthesized successfully. They exhibit strong electronic correlations31 and even superconductivity32. The fact that the number of available strongly correlated Kagome metals is growing significantly is encouraging in view of experimentally realizing our results, and outweighs the difficulties so far encountered in achieving conducting Herbertsmithite materials33,34,35.
Methods
Properties of Herbertsmithite
Kagome materials, such as Herbertsmithite combine the features of Dirac fermions and strong correlations. ScCu3(OH)6Cl2, which we refer to as Sc-Herbertsmithite, consists of CuO4 plaquettes forming a Kagome lattice (see Fig. 2a). The crystal field is such that the low-energy physics is correctly captured by a single \({d}_{{{\rm{x}}}^{2}-{{\rm{y}}}^{2}}\) orbital on each Cu site, and the resulting low-energy band structure is qualitatively similar to that of a one-orbital model at n = 4/3 electron filling, where, as shown in Fig. 2b, the Fermi level is pinned at the Dirac points.
As compared to graphene, where the underlying Dirac spectrum originates from weakly correlated pz orbitals, Sc-Herbertsmithite is expected to show a larger degree of electronic correlations. This is confirmed by our constrained random phase approximation (cRPA) estimates of the coupling strength α. For a linear band dispersion, the strength of the Coulomb interaction is characterized by the effective fine-structure constant in Eq. (3). According to the values summarized in Table 1, αScâHbâ=â 2.9â±â0.3 as compared to αGrâ=â0.5âââ1.0 for (hBN-encapsulated) graphene. We use αGrâ=â0.9 as a representative example in the main text. The latter value was evaluated using an ab-initio estimate of the dielectric constant36,37.
First principles calculations
For our numerical study of Sc-Herbertsmithite, we employed state-of-the-art first-principle calculations based on the density functional theory as implemented in the Vienna ab initio simulation package (VASP)38 following the projector-augmented-plane-wave (PAW) method39,40. We use the generalized gradient approximation as parametrized by the PBE-GGA functional for the exchange-correlation potential41 by expanding the Kohn-Sham wave functions into plane waves up to an energy cutoff of 400 eV and sampling the Brillouin zone on an 6âÃâ6âÃâ6 regular mesh. We obtain the phonon dispersion Ïνq within the context of the density functional perturbation theory (DFPT)42 as implemented in the Quantum ESPRESSO suite43 with a 2âÃâ2âÃâ2 supercell. The electron-phonon coupling strengths λνqâ=â Î â³Î½q/ÏN(ϵF)Ïνq are computed from the imaginary part Î â³Î½q of the phonon self-energy within the Migdal approximation
where \({g}_{{\rm{mn}},\nu }^{{\bf{k}},{\bf{q}}}\) are the electron-phonon matrix elements and the electron momentum integration is performed on an extreme dense mesh of Nkâ=â503 points via Wannier interpolation44.
We subsequently determine the Fermi velocity by a fit to the band structure along the Î-X-J-Î path. The low-energy theory of Herbertsmithite consists of the Cu \({d}_{{{\rm{x}}}^{2}-{{\rm{y}}}^{2}}\) Dirac electrons and their hydrodynamic properties are characterized by the Coulomb coupling. This effective low-energy theory needs to take into account the effects of all other electronic states, in particular the dielectric screening. In other words, the hydrodynamic theory starts from partially dressed Cu \({d}_{{{\rm{x}}}^{2}-{{\rm{y}}}^{2}}\) electrons, with an effective dielectric constant ϵr that captures the screening by all other electronic states. Since these other states are far away from the Fermi level, they are less affected by correlation effects and the ab-initio calculation of their dielectric constant is feasible. To determine ϵr, we start from the density functional theory band structure and fix (constrain) the occupation of all electronic states in the three Kagome bands to be nâ=â4/3 (including a factor of 2 from spin degeneracy), thereby excluding intraband screening processes45. To first approximation, the resulting dielectric tensor is diagonal and given by a single constant ϵrâ=â5 (see Table 1), with relative deviations of up to 10% from this constant value. We use this 10% as an estimate of the methodological uncertainty. The cRPA dielectric constant and the Fermi velocity from density functional theory define the partially dressed electronic states entering the hydrodynamic theory, which captures the low-energy correlations in the material and determines the resulting viscosity.
Phonons
A number of caveats need to be addressed. Interactions of the electron fluid with lattice vibrations are detrimental to the hydrodynamic behavior. As shown in Fig. 2c, however, the optical phonon modes in Sc-Herbertsmithite that sensibly couple to the electronic states, and indeed contribute to the electron-phonon coupling, are populated for temperatures above Tphâ~â80âK. This analysis indicates that the hydrodynamic regime can be expected to extend over a range of temperatures within current experimental reach.
Instabilities
The Kagome lattice itself plays a crucial role. In contrast to grapheneâs Dirac spectrum, which occurs precisely at half-filling, and as such it is prone to magnetic instabilities at comparably small couplings20, the Dirac fluid in Kagome Sc-Herbertsmithite is reached at the stoichiometric filling nâ=â4/3. Therefore, tendencies towards a correlation driven magnetic groundstate are strongly suppressed, and the linear dispersion of the Dirac fluid is expected to be highly robust against the opening of a band gap. While the geometric frustration inherent in the Kagome lattice precludes a rigorous treatment, this conjecture is supported both from weak coupling18,46 as well as strong coupling studies19.
Hydrodynamics
Electronic fluids are in local thermal equilibrium on length scales larger than the interaction mean free path âee, which needs to be the smallest length scale for hydrodynamics to apply (see Eq. (4)). Locally, the laws of thermodynamics such as the Gibbs-Duhem relation
holds for the densities of energy ϵ(T, μ), entropy s(T, μ) and charge Ï(T, μ). Given the equation of state P(ϵ, Ï), Eq. (8) can be used to determine the entropy density s(μ, T). A two-dimensional Dirac fluid, i.e. a fluid of electrons with relativistic dispersion and the chemical potential pinned to the Dirac point, μâ=â0, is then described by relativistic hydrodynamics22. The hydrodynamic equations are the conservation equations of the energy-momentum tensor Tμν and the electric current Jμ,
where we have also allowed for the coupling of the fluid to an external electromagnetic field, Fνμ.
The hydrodynamic derivative expansion expresses Tμν and Jμ order by order in derivatives of the local temperature T(xν), the chemical potential μ(xν), and the relativistic 3-velocity uμ(xν) (uμuμâ=ââc2). For a parity and time reversal invariant, relativistic two-dimensional fluid, one obtains to first order in the derivatives
where Îμνâ=âuμuν/c2 + ημν is the projection matrix in the directions transverse to uμ and A(μν)â=â(Aμνâ+âAνμ)/2 denotes symmetrization of any two-tensor.
The shear viscosity η is the main observable in this context. It can be calculated from kinetic theory in the perturbative regime αââªâ147, and will be derived from a holographic model below in the nonperturbative regime αââ«â1. The bulk viscosity ξ is negligible due to the approximate scale invariance of the linear dispersion relation in a Dirac fluid, ξâââ0. For incompressible flows (âμuμâ=â0) such as the Poiseuille flow, a finite bulk viscosity has no effect either. The interaction-induced intrinsic conductivity ÏQ48 affects electric transport and controls the rate of Joule heating, but does not affect the flow of the fluid.
The three transport coefficients (η, ξ, ÏQ) have to be calculated from a microscopic theory or an effective model, such as the kinetic theory at weak coupling or holography at strong coupling. In linear response theory around a global equilibrium (i.e., constant temperature T and chemical potential μ, and uμâ=â(1, 0, 0)T), the shear viscosity η can be calculated with the following Kubo formula at zero frequency Ï and momentum k,
In the strong coupling regime, the AdS/CFT correspondence provides a framework to calculate η via a semi-classical analysis of the dual gravitational action of the system49,50. The entropy density of the fluid is then given by the entropy density of the black hole horizon in the bulk of spacetime51.
Holographic model
At infinite coupling and for a large number of degrees of freedom, the holographic dual is given by the Einstein-Hilbert action (in 2+1d dimensional field theory the coupling is dimensionfull and hence runs. This is usually reflected in holography by including a non-trivial dilaton scalar field, Ï, into the gravitational action. We have neglected the dilaton contribution in our analysis, because it leads to sub-leading derivative corrections to η/s)52.
yielding Eq. (1) for the ratio η/s, independently of the coupling constant and number of degrees of freedom. Our corrections to this limit are computed by adding higher powers of the curvature R to Eq. (15). In 3â+â1d gravitational duals, the next-to-leading order R2 terms contribute a topological Gauss-Bonnet term to the gravitational action53 and therefore do not alter the value of η/s. Furthermore, type-II supergravity, the ten-dimensional parent theory of our four-dimensional gravity dual, does not contain R3 corrections, as was shown in ref. 54. The next higher derivative corrections contain four powers of R14,55,56 and are induced in type-II supergravity by quartic terms involving the Weyl tensor. These terms yield
In top-down constructions of holography originating from string theory, the holographic dual theory is typically a non-Abelian gauge theory with gauge group rank N and ât Hooft coupling λâ=âαN, related to the fine-structure constant α. N parametrizes the number of degrees of freedom in the dual gauge theory. Dimensional analysis implies that the correction at order R4 to η/s scales as λâ3/2 in any spacetime dimension, up to a multiplicative constant that depends on the details of the holographic model49. Eq. (16) in particular universally describes the coupling dependence of the ratio η/s for all (3â+â1)d-holographic duals with relativistic symmetry at leading order in the inverse coupling expansion and in the large N limit near charge neutrality. However, the prefactors of the allowed subleading curvature corrections are model dependent. We parametrize the model dependence of the R4 correction through the prefactor \({\mathcal{C}}^{\prime}\) in Eq. (16). For the original example of a theory with holographic dual, \({\mathcal{N}}=4\) SYM theory13, \({\mathcal{C}}^{\prime} =135\zeta (3)/8\)55. The unknowns in Eq. (16) are N and \({\mathcal{C}}^{\prime}\), with \({\mathcal{C}}^{\prime}\) depending on the particular holographic dual considered. We absorb Nâ3/2 into \({\mathcal{C}}^{\prime}\) via \({\mathcal{C}}\equiv {\mathcal{C}}^{\prime} {N}^{-3/2}\), which brings Eq. (16) to the form of Eq. (2). Thus, the details of the holographic model suitable for describing Dirac fluids are encoded in a single coefficient, namely \({\mathcal{C}}\). Using the parametrization in terms of \({\mathcal{C}}\), we generate the blue band in Fig. 1 as follows: We vary \({\mathcal{C}}\) from \({\mathcal{C}}=0.0005\) to \({\mathcal{C}}=5\) to generate the blue band in Fig. 1. For \({\mathcal{N}}=4\) SYM, this range of variation of \({\mathcal{C}}\) corresponds to formally varying N from Nâââ103 to Nâââ2. As shown in Fig. 1, even varying \({\mathcal{C}}\) over four orders of magnitude changes the value of η/s of Sc-Hb by at most a factor of two. In principle, there may be corrections of even higher order in 1/λ to Eq. (16), corresponding to terms involving even higher orders in the curvature. However, we do not expect these to alter the predicted estimate of the value of η/s for values of the coupling not accessible to weak-coupling perturbation theory. Furthermore, the absence of R2 and R3 corrections to the gravitational action imply that Eq. (16) does not receive 1/N corrections53. Hence, corrections to η/s independent of α can appear at order 1/N2 or higher, however, we expect these to be subleading as compared to the \({\mathcal{C}}/{\alpha }^{3/2}\) corrections at αScâHbâ=â2.9, for the range of \({\mathcal{C}}\) chosen. Moreover, our main argument that the applicability of hydrodynamics and the appearance of turbulence are more likely at larger values of α is independent of the value of N. The black, fading, line in Fig. 1 shows the perturbative prediction at weak coupling47,48, which cannot be extrapolated to stronger coupling without violating the condition ηââ¥â0 that follows from local application of the second law of thermodynamics
Turbulence
Using the extrapolation presented in Fig. 1, we infer that materials which display a stronger Coulomb interaction, such as Sc-Herbertsmithite, yield a much larger Reynolds number than that reported for graphene. In particular, using Eq. (6), Fig. 1 and the results for the Fermi velocity vF from Table 1, we find that Sc-Herbertsmithite will exhibit a 63 to 156 times larger Reynolds number as compared to graphene, depending on whether we use the value of η/s at the bottom or the top of the blue band, respectively. This enhancement is large enough to expect fully developed turbulence in a constriction setup26.
Coupling dependence of the regime of validity of hydrodynamics
In this section, we derive Eq. (5) for the relative ratio between the length scales where hydrodynamics is expected to break down for Sc-Herbertsmithite and graphene. One can characterize the breakdown of hydrodynamics by the existence of a diffusive pole, Ïâ=ââiDk2, in the energy-momentum tensor self-correlation in Eq. (14). Through the calculation of the poles of the mentioned self-correlation function via holography, the diffusive pole was shown57 to disappear at a critical wavelength âc as the coupling strength of the system is decreased. The pole moves down the imaginary frequency axis and collides with an upward-moving pole at ââ= âc. After the collision, the two poles split into a conjugate pair and the purely imaginary diffusion pole disappears. An approximate analytic result for lc is given by57
where λ is the ât Hooft coupling defined below Eq. (16). Within the model employed57, the analytical approximation in Eq. (17) lies closer to the numerical value of âc for smaller values of λ. To our knowledge, âc was derived only through holographic methods. At weak coupling, it is only known that hydrodynamics has a finite radius of convergence. For the particular case of graphene and Sc-Herbertsmithite, with αGrâ=â0.9 and αScâHbâ=â2.9, we find the ratio given in Eq. (5). Note that we calculated only the ratio of critical wavelengths and not their absolute values. We did so because the two materials have the same low-lying energy spectra and thus Eq. (5) is independent of the unknowns \(N,{\mathcal{C}}\) of Eq. (16). Eq. (5) shows that if both graphene and Sc-Herbertsmithite are described by holography, then the hydrodynamic approximation is more robust for Sc-Herbertsmithite. This argument presumes that graphene can be described holographically. Because of grapheneâs intermediate coupling strength, however, αGrâ=â0.9, it is more likely to lie in the intermediate regime, where neither simple truncated holographic models nor perturbative methods are applicable. In this case Eq. (5) entails that Sc-based Herbertsmithite is a better candidate for the realization of holographic hydrodynamics.
Data availability
All data are available from the corresponding author upon request.
Code availability
All codes used in this work, with the exception of the licensed VASP package, are either publicly available or available from the authors upon request.
References
Andreev, A. V., Kivelson, S. A. & Spivak, B. Hydrodynamic description of transport in strongly correlated electron systems. Phys. Rev. Lett.106, 256804 (2011).
Polini, M. & Geim, A. K. Viscous electron fluids. Physics Today73, 28 (2020).
Molenkamp, L. W. & de Jong, M. J. M. Electron-electron-scattering-induced size effects in a two-dimensional wire. Phys. Rev. B49, 5038 (1994).
Molenkamp, L. W. & de Jong, M. J. M. Observation of Knudsen and Gurzhi transport regimes in a two-dimensional wire. Solid State Electron.37, 551 (1994).
Moll, P. J. W., Kushwaha, P., Nandi, N., Schmidt, B. & Mackenzie, A. P. Evidence for hydrodynamic electron flow in PdCoO2. Science351, 1061 (2016).
Gooth, J. et al. Thermal and electrical signatures of a hydrodynamic electron fluid in tungsten diphosphide. Nat. Commun.9, 4093 (2018).
Bandurin, D. A. et al. Fluidity onset in graphene. Nat. Commun.9, 4533 (2018).
Torre, I., Tomadin, A., Geim, A. K. & Polini, M. Nonlocal transport and the hydrodynamic shear viscosity in graphene. Phys. Rev. B92, 165433 (2015).
Levitov, L. & Falkovich, G. Electron viscosity, current vortices and negative nonlocal resistance in graphene. Nat. Phys.12, 672 (2016).
Bandurin, D. et al. Negative local resistance caused by viscous electron backflow in graphene. Science351, 1055 (2016).
Pellegrino, F. M. D., Torre, I., Geim, A. K. & Polini, M. Electron hydrodynamics dilemma: whirlpools or no whirlpools. Phys. Rev. B94, 155414 (2016).
Kumar, R. K. et al. Superballistic flow of viscous electron fluid through graphene constrictions. Nat. Phys.13, 1182 (2017).
Maldacena, J. M. The large N limit of superconformal field theories and supergravity. Adv. Theor. Math. Phys.2, 231 (1998).
Policastro, G., Son, D. T. & Starinets, A. O. The shear viscosity of strongly coupled N = 4 supersymmetric Yang-Mills plasma. Phys. Rev. Lett.87, 081601 (2001).
Buchel, A. & Liu, J. T. Universality of the shear viscosity from supergravity duals. Phys. Rev. Lett.93, 090602 (2004).
Kovtun, P. K., Son, D. T. & Starinets, A. O. Viscosity in strongly interacting quantum field theories from black hole physics. Phys. Rev. Lett.94, 111601 (2005).
Novoselov, K. S. et al. Two-dimensional gas of massless Dirac fermions in graphene. Nature438, 197 (2005).
Kiesel, M. L., Platt, C. & Thomale, R. Unconventional fermi surface instabilities in the Kagome Hubbard model. Phys. Rev. Lett.110, 126405 (2013).
Indergand, M., Läuchli, A., Capponi, S. & Sigrist, M. Modeling bond-order wave instabilities in doped frustrated antiferromagnets: valence bond solids at fractional filling. Phys. Rev. B74, 064429 (2006).
Assaad, F. F. & Herbut, I. F. Pinning the order: the nature of quantum criticality in the Hubbard model on honeycomb lattice. Phys. Rev. X3, 031010 (2013).
Mazin, I. I. et al. Theoretical prediction of a strongly correlated Dirac metal. Nat. Commun.5, 4261 (2014).
Landau, L. & Lifshitz, E. Fluid Mechanics. Course of Theoretical Physics (Butterworth-Heinemann, 1987).
Adams, A., Chesler, P. M. & Liu, H. Holographic turbulence. Phys. Rev. Lett.112, 151602 (2014).
Erdmenger, J., Matthaiakakis, I., Meyer, R. & RodrÃguez Fernández, D. Strongly coupled electron fluids in the Poiseuille regime. Phys. Rev. B98, 195143 (2018).
Tuckerman, L. S., Chantry, M. & Barkley, D. Patterns in wall-bounded shear flows. Annu. Rev. Fluid Mech.52, 1 (2020).
Gabbana, A., Polini, M., Succi, S., Tripiccione, R. & Pellegrino, F. M. D. Prospects for the detection of electronic preturbulence in graphene. Phys. Rev. Lett.121, 236602 (2018).
Galitski, V., Kargarian, M. & Syzranov, S. Dynamo effect and turbulence in hydrodynamic weyl metals. Phys. Rev. Lett.121, 176603 (2018).
Link, J. M., Narozhny, B. N., Kiselev, E. I. & Schmalian, J. Out-of-bounds hydrodynamics in anisotropic dirac fluids. Phys. Rev. Lett.120, 196801 (2018).
Liao, Y. & Galitski, V. Critical viscosity of a fluctuating superconductor. Phys. Rev. B100, 060501 (2019).
Zaanen, J. Planckian dissipation, minimal viscosity and the transport in cuprate strange metals. SciPost Phys.6, 61 (2019).
Fuchs, M. et al. Kagome metal-organic frameworks as a platform for strongly correlated electrons. J. Phys. Mater.3, 025001 (2020).
Huang, X. et al. Superconductivity in a copper(ii)-based coordination polymer with perfect kagome structure. Angew. Chem., Int. Ed.57, 146 (2018).
Kelly, Z. A., Gallagher, M. J. & McQueen, T. M. Electron doping a kagome spin liquid. Phys. Rev. X6, 041007 (2016).
Liu, Q. et al. Electron doping of proposed kagome quantum spin liquid produces localized states in the band gap. Phys. Rev. Lett.121, 186402 (2018).
Puphal, P. et al. Tuning of a kagome magnet: Insulating ground state in ga-substituted Cu4(OH)6Cl2. Phys. Stat. Sol. (b)256, 1800663 (2019).
Rösner, M., ÅaÅÃoÄlu, E., Friedrich, C., Blügel, S. & Wehling, T. O. Wannier function approach to realistic Coulomb interactions in layered materials and heterostructures. Phys. Rev. B92, 085102 (2015).
Wehling, T. O. et al. Strength of effective Coulomb interactions in graphene and graphite. Phys. Rev. Lett.106, 236805 (2011).
Kresse, G. & Furthmüller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. Rev. B54, 11169 (1996).
Kresse, G. & Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. Rev. B59, 1758 (1999).
Blöchl, P. E. Projector augmented-wave method. Phys. Rev. B50, 17953 (1994).
Perdew, J. P., Burke, K. & Ernzerhof, M. Generalized gradient approximation made simple. Phys. Rev. Lett.77, 3865 (1996).
Baroni, S., de Gironcoli, S., Dal Corso, A. & Giannozzi, P. Phonons and related crystal properties from density-functional perturbation theory. Rev. Mod. Phys.73, 515 (2001).
Giannozzi, P. et al. QUANTUM ESPRESSO: a modular and open-source software project for quantum simulations of materials. J. Phys. Cond. Matter21, 395502 (2009).
Poncé, S., Margine, E., Verdi, C. & Giustino, F. Epw: Electronâphonon coupling, transport and superconducting properties using maximally localized wannier functions. Comp. Phys Commun.209, 116 (2016).
Aryasetiawan, F. et al. Frequency-dependent local interactions and low-energy effective models from electronic structure calculations. Phys. Rev. B70, 195104 (2004).
Kiesel, M. L. & Thomale, R. Sublattice interference in the kagome Hubbard model. Phys. Rev. B86, 121105 (2012).
Müller, M., Schmalian, J. & Fritz, L. Graphene - a nearly perfect fluid. Phys. Rev. Lett.103, 025301 (2009).
Fritz, L., Schmalian, J., Müller, M. & Sachdev, S. Quantum critical transport in clean graphene. Phys. Rev. B78, 085416 (2008).
Cremonini, S. The shear viscosity to entropy ratio: a status report. Mod. Phys. Lett. B25, 1867 (2011).
Buchel, A., Myers, R. C. & Sinha, A. Beyond η/sâ=â1/4Ï. J. High Energy Phys.03, 084 (2009).
Gubser, S. S., Klebanov, I. R. & Tseytlin, A. A. Coupling constant dependence in the thermodynamics of Nâ=â4 supersymmetric Yang-Mills theory. Nuc. Phys. B534, 202 (1998).
Itzhaki, N., Maldacena, J. M., Sonnenschein, J. & Yankielowicz, S. Supergravity and the large N limit of theories with sixteen supercharges. Phys. Rev.D58, 046004 (1998).
Brigante, M., Liu, H., Myers, R. C., Shenker, S. & Yaida, S. Viscosity bound violation in higher derivative gravity. Phys. Rev.D77, 126006 (2008).
Metsaev, R. R. & Tseytlin, A. A. Curvature cubed terms in string theory effective actions. Phys. Lett. B185, 52 (1987).
Buchel, A., Liu, J. T. & Starinets, A. O. Coupling constant dependence of the shear viscosity in N = 4 supersymmetric Yang-Mills theory. Nucl. Phys. B707, 56 (2005).
Benincasa, P. & Buchel, A. Transport properties of N = 4 supersymmetric Yang-Mills theory at finite coupling. J. High Energy Phys.01, 103 (2006).
Grozdanov, S., Kaplis, N. & Starinets, A. O. From strong to weak coupling in holographic models of thermalization. J. High Energy Phys.2016, 151 (2016).
Acknowledgements
We thank Andy Lucas and Fabio Caruso for useful discussions. The work in Würzburg is funded by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) through Project-ID 258499086âSFB 1170 and through the Würzburg-Dresden Cluster of Excellence on Complexity and Topology in Quantum Matter âct.qmat Project-ID 390858490âEXC 2147. We acknowledge the Gauss Centre for Supercomputing e.V. for funding this project by providing computing time on the GCS Supercomputer SuperMUC at Leibniz Supercomputing Centre. Open access funding provided by Projekt DEAL.
Author information
Authors and Affiliations
Contributions
J.E., M.G., R.M., and R.T. initialized the project. J.E. and R.T. coordinated and supervised the investigation. I.M., R.M., and D.R.F. developed and performed the AdS-CFT analysis, while D.D.S., E.v.L., and T.W. performed the ab-initio and constrained RPA calculations. I.M., D.D.S., J.E., M.G., R.M., D.R.F., and R.T. wrote the paper. The paper reflects contributions from all authors.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Peer review informationNature Communications thanks the anonymous reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.
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
Di Sante, D., Erdmenger, J., Greiter, M. et al. Turbulent hydrodynamics in strongly correlated Kagome metals. Nat Commun 11, 3997 (2020). https://doi.org/10.1038/s41467-020-17663-x
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41467-020-17663-x
This article is cited by
-
Charge transport and hydrodynamics in materials
Nature Reviews Materials (2023)
-
Para-hydrodynamics from weak surface scattering in ultraclean thin flakes
Nature Communications (2023)
-
Flat band separation and robust spin Berry curvature in bilayer kagome metals
Nature Physics (2023)
-
Dynamics and resilience of the unconventional charge density wave in ScV6Sn6 bilayer kagome metal
Communications Materials (2023)
-
Topological kagome magnets and superconductors
Nature (2022)
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.