## Abstract

Interactions between particles in quantum many-body systems can lead to collective behavior described by hydrodynamics. One such system is the electron-hole plasma in graphene near the charge neutrality point, which can form a strongly coupled Dirac fluid. This charge neutral plasma of quasi-relativistic fermions is expected to exhibit a substantial enhancement of the thermal conductivity, thanks to decoupling of charge and heat currents within hydrodynamics. Employing high sensitivity Johnson noise thermometry, we report an order of magnitude increase in the thermal conductivity and the breakdown of the Wiedemann-Franz law in the thermally populated charge neutral plasma in graphene. This result is a signature of the Dirac fluid, and constitutes direct evidence of collective motion in a quantum electronic fluid.

Understanding the dynamics of many interacting particles is a formidable task in physics. For electronic transport in matter, strong interactions can lead to a breakdown of the Fermi liquid (FL) paradigm of coherent quasiparticles scattering off of impurities. In such situations, provided certain conditions are met, the complex microscopic dynamics can be coarse-grained to a hydrodynamic description of momentum, energy, and charge transport on long length and time scales (*1*). Hydrodynamics has been successfully applied to a diverse array of interacting quantum systems, from high mobility electrons in conductors (*2*), to cold atoms (*3*) and quark-gluon plasmas (*4*). Hydrodynamic effects are expected to greatly modify transport coefficients as compared to their FL counterparts, as has been argued for strongly interacting massless Dirac fermions in graphene at the charge-neutrality point (CNP) (*5*–*8*).

Many-body physics in graphene is interesting because of electron-hole symmetry and a linear dispersion relation at the CNP (*9*, *10*). Together with the vanishing Fermi surface, the ultrarelativistic spectrum leads to ineffective screening (*11*) and the formation of a strongly-interacting quasi-relativistic electron-hole plasma, known as a Dirac fluid (*12*). The Dirac fluid shares many features with quantum critical systems (*13*): most importantly, the electron-electron scattering time is fast (*14*–*17*), and well suited to a hydrodynamic description. Because of the quasi-relativistic nature of the Dirac fluid, this hydrodynamic limit is described by equations (*18*) quite different from non-relativistic counterparts. A number of exotic properties have been predicted including nearly perfect (inviscid) flow (*19*) and a diverging thermal conductivity resulting in the breakdown of the Wiedemann-Franz law at finite temperature (*5*, *6*).

Away from the CNP, graphene has a sharp Fermi surface and the standard FL phenomenology holds. By tuning the chemical potential, we may measure thermal and electrical conductivity in both the Dirac fluid (DF) and the FL in the same sample. In a FL, the relaxation of heat and charge currents is closely related as they are carried by the same quasiparticles. The Wiedemann-Franz (WF) law (*20*) states that the electronic contribution to a metal’s thermal conductivity κ_{e} is proportional to its electrical conductivity σ and temperature *T*, such that the Lorenz ratio satisfies(1)where *e* is the electron charge, *k*_{B} is the Boltzmann constant, and is the Sommerfeld value derived from FL theory. depends only on fundamental constants, and not on specific details of the system such as carrier density or effective mass. As a robust prediction of FL theory, the WF law has been verified in numerous metals (*20*). At high temperature, the WF law can be violated due to inelastic electron-phonon scattering or bipolar diffusion in semiconductors, even when electron-electron interactions are negligible (*21*). In recent years, several non-trivial violations of the WF law have been reported in strongly interacting systems such as Luttinger liquids (*22*), metallic ferromagnets (*23*), heavy fermion metals (*24*), and underdoped cuprates (*25*), all related to the emergence of non-Fermi liquid behavior.

The WF law is expected to be violated at the CNP in a DF thanks to the strong Coulomb interactions between thermally excited charge carriers. An electric field drives electrons and holes in opposite directions; collisions between them introduce a frictional dissipation, resulting in a finite conductivity even in the absence of disorder (*26*). In contrast, a temperature gradient causes electrons and holes to move in the same direction inducing an energy current, which grows unimpeded by inter-particle collisions as the total momentum is conserved. The thermal conductivity is therefore limited by the rate at which momentum is relaxed byresidual impurities.

Realization of the Dirac fluid in graphene requires that the thermal energy be larger than the local chemical potential μ(**r**), defined at position **r**: . Impurities cause spatial variations in the local chemical potential, and even when the sample is globally neutral, it is locally doped to form electron-hole puddles with finite μ(**r**) (*27*–*30*). At high temperature, formation of the DF is complicated by phonon scattering, which can relax momentum by creating additional inelastic scattering channels. This high temperature limit occurs when the electron-phonon scattering rate becomes comparable to the electron-electron scattering rate. These two temperatures set the experimental window in which the DF and the breakdown of the WF law can be observed.

To minimize disorder, the monolayer graphene samples used in this report are encapsulated in hexagonal boron nitride (hBN) (*31*). All devices used in this study are two-terminal to keep a well-defined temperature profile (*32*) with contacts fabricated using the one-dimensional edge technique (*33*) in order to minimize contact resistance. We employ a back gate voltage applied to the silicon substrate to tune the charge carrier density *n* = *n*_{e} – *n*_{h}, where *n*_{e} and *n*_{h} are the electron and hole density, respectively (*21*). All measurements are performed in a cryostat controlling the temperature *T*_{bath}. Figure 1A shows the resistance *R* versus *V*_{g} measured at various fixed temperatures for a representative device [see (*21*) for all samples]. From this, we estimate the electrical conductivity σ (Fig. 1B) using the known sample dimensions. At the CNP, the residual charge carrier density *n*_{min} can be estimated by extrapolating a linear fit of log(σ) as a function of log(*n*) out to the minimum conductivity (*34*). At the lowest temperatures we find *n*_{min} saturates to ~8 × 10^{9} cm^{–2}. We note that the extraction of *n*_{min} by this method overestimates the charge puddle energy, consistent with previous reports (*31*). Above the disorder energy scale *T*_{dis} ~ 40 K, *n*_{min} increases as *T*_{bath} is raised, suggesting thermal excitations begin to dominate and the sample enters the non-degenerate regime near the CNP.

The electronic thermal conductivity is measured using high sensitivity Johnson noise thermometry (JNT) (*32*, *35*). We apply a small bias current through the sample that injects a joule heating power *P* directly into the electronic system, inducing a small temperature difference between the graphene electrons and the bath. The electron temperature *T*_{e} is monitored independent of the lattice temperature through the Johnson noise power emitted at 100 MHz with a 20 MHz bandwidth defined by an LC matching network. We designed our JNT to be operated over a wide temperature range 3–300 K (*35*). With a precision of ~10 mK, we measure small deviations of *T*_{e} from *T*_{bath}, i.e., . In this limit, the temperature of the graphene lattice is well thermalized to the bath (*32*) and our JNT setup allows us to sensitively measure the electronic cooling pathways in graphene. At low enough temperatures, electron and lattice interactions are weak (*35*, *36*), and most of the Joule heat generated in graphene escapes via direct diffusion to the contacts (*21*). As temperature increases, electron-phonon scattering becomes appreciable and thermal transport becomes limited by the electron-phonon coupling strength (*36*–*38*). The onset temperature of appreciable electron-phonon scattering, *T*_{el–ph}, depends on the sample disorder and device geometry: *T*_{el–ph} ~ 80 K (*35*, *36*, *39*, *40*) for our samples. Below this temperature, the electronic contribution of the thermal conductivity can be obtained from *P* and Δ*T* using the device dimensions (*21*).

Figure 1C plots κ_{e}(*V*_{g}) alongside the simultaneously measured σ(*V*_{g}) at various fixed bath temperatures. Here, for a direct quantitative comparison based on the WF law, we plot the scaled electrical conductivity as in the same units as κ_{e}; in a FL these two values will coincide as given by Eq. 1. At low temperatures, *T* < *T*_{dis} ~ 40 K, where the puddle induced density fluctuations dominate, we find , monotonically increasing as a function of carrier density with a minimum at the neutrality point, confirming the WF law in the disordered regime. As *T* increases (*T* > *T*_{dis}), however, we begin to see the violation of the WF law. This violation appears only close to the CNP, with the measured thermal conductivity maximized at *n* = 0 (Fig. 1C). The deviation is the largest at 75 K, where κ_{e} is over an order of magnitude larger than the value expected for a FL. The non-monotonicity of κ_{e}(*T*) is consistent with acoustic phonons relaxing momentum more efficiently than impurities as *T* increases (*41*); for in our samples, activation of optical phonons introduces an additional electron-phonon cooling pathway (*35*), and the measured thermal conductivity is larger than κ_{e}. This non-FL behavior quickly disappears as increases; κ_{e} returns to the FL value and restores the WF law. In fact, away from the CNP, the WF law holds for a wide temperature range, consistent with previous reports (*35*, *36*, *39*) (Fig. 1E). For this FL regime, we verify the WF law up to *T* ~ 80 K.

Our observation of the breakdown of the WF law in graphene is consistent with the emergence of the DF. Figure 2 shows the full density and temperature dependence of the experimentally measured Lorenz ratio highlighting the presence of the DF. The blue colored region denotes , suggesting the carriers in graphene exhibit FL behavior. The WF law is violated in the DF (yellow-red) with a peak Lorenz ratio 22 times larger than . The green dotted line shows the corresponding *n*_{min}(*T*) for this sample; the DF is found within this regime, indicating the coexistence of thermally populated electrons and holes. Disorder and electron-acoustic phonon scattering serve as lower and upper limits respectively on the temperature range over which the DF can be observed.

We investigate the effect of impurities on hydrodynamic transport by comparing the results obtained from samples with varying disorder. Figure 3A shows *n*_{min} as a function of temperature for three samples used in this study. *n*_{min}(*T *= 0) is estimated as 5, 8, and 10 × 10^{9} cm^{–2} in samples S1, S2, and S3, respectively. All devices show qualitatively similar Dirac fluid behavior; the largest value of measured in the Dirac fluid regime is 22, 12, and 3 in samples S1, S2, and S3, respectively (Fig. 3B). Cleaner samples not only have a more pronounced peak but also a narrower density dependence (Fig. 3C), as predicted (*5*, *6*).

More quantitative analysis of in our experiment can be done by employing a quasi-relativistic hydrodynamic theory of the DF incorporating the effects of weak impurity scattering (*5*, *6*, *18*)(2)where(3)Here *v*_{F} is the Fermi velocity in graphene, σ_{min} is the electrical conductivity at the CNP, is the fluid enthalpy density, and *l*_{m} is the momentum relaxation length from impurities. Two parameters in Eq. 2 are undetermined for any given sample: *l*_{m} and . For simplicity, we assume we are well within the DF limit where *l*_{m} and are approximately independent of *n*. We fit Eq. 2 to the experimentally measured for all temperatures and densities in the Dirac fluid regime to obtain *l*_{m} and for each sample (see, e.g., Fig. 3C). *l*_{m} is estimated to be 1.5, 0.6, and 0.034 μm for samples S1, S2, and S3, respectively. For the system to be well described by hydrodynamics, *l*_{m} should be long compared to the electron-electron scattering length of ~0.1 μm expected for the Dirac fluid at 60 K (*19*). This is consistent with the pronounced signatures of hydrodynamics in S1 and S2, but not in S3, where only a glimpse of the DF appears in this more disordered sample. We also observe in S1 that dips substantially below : its minimum is . occurs in Eq. 2 for . The inset to Fig. 3C shows the fitted enthalpy density as a function of temperature compared to that expected in clean graphene (dashed line) (*19*), excluding renormalization of the Fermi velocity. In the cleanest sample varies from 1.1-2.3 eV/μm^{2} in the hydrodynamic regime. This enthalpy density corresponds to ~20 meV or ~4*k*_{B}T per charge carrier—about a factor of 2 larger than the model calculation without disorder (*19*). The sharp temperature and impurity dependence observed in is a prediction of these hydrodynamic models. These effects and the magnitude of are inconsistent with alternative models for WF violations, including bipolar diffusion in graphene (*21*, *42*). Furthermore, recent experiments report monotonic behavior in thermopower as a function of *T* (*43*), implying phonon drag is not responsible for the peak in κ_{e} that we observe as a function of *T*.

To fully incorporate the effects of disorder, a hydrodynamic theory treating inhomogeneity non-perturbatively is necessary (*41*, *44*). The enthalpy densities reported here are larger than the theoretical estimation obtained for disorder free graphene, consistent with the picture that chemical potential fluctuations prevent the sample from reaching the Dirac point. While we find thermal conductivity well described by (*5*, *6*), electrical conductivity increases slower than expected away from the CNP, a result consistent with hydrodynamic transport in a viscous fluid with charge puddles (*41*).

In a hydrodynamic system, the ratio of shear viscosity η to entropy density *s* is an indicator of the strength of the interactions between constituent particles. It is suggested that the DF can behave as a nearly perfect fluid (*19*): η/*s* approaches a conjecture by Kovtun-Son-Starinets: for a strongly interacting system (*45*). A non-perturbative hydrodynamic framework can be employed to estimate η (*41*). A direct measurement of η is of great interest.

Beyond a diverging thermal conductivity and an ultra-low viscosity, other peculiar phenomena are expected to arise in this plasma. The massless nature of the Dirac fermions is expected to result in a large kinematic viscosity, despite a small shear viscosity η. Observable hydrodynamic effects have also been predicted to extend into the FL regime (*46*). The study of magnetotransport in the DF will lead to further tests of hydrodynamics (*5*, *18*).

## Supplementary Materials

## References and Notes

- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵Materials and methods are available as supplementary materials on
*Science*Online. - ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
**Acknowledgments:**We thank Matthew Foster, Dmitri Efetov and Gil-Ho Lee for helpful discussions. The major experimental work at Harvard is supported by DOE (DE-SC0012260) and at Raytheon BBN Technologies is supported by IRAD. J.C. thanks the support of the FAME Center, sponsored by SRC MARCO and DARPA. K.W. is supported by ARO MURI (W911NF-14-1-0247). J.K.S. is supported by ARO (W911NF-14-1-0638) and AStar. P.K. acknowledges partial support from the Gordon and Betty Moore Foundation’s EPiQS Initiative through Grant GBMF4543 and Nano Material Technology Development Program through the National Research Foundation of Korea (2012M3A7B4049966). A.L. and S.S. are supported by the NSF under Grant DMR-1360789, the Templeton foundation, and MURI grant W911NF-14-1-0003 from ARO. Research at Perimeter Institute is supported by the Government of Canada through Industry Canada and by the Province of Ontario through the Ministry of Research and Innovation. K.W. and T.T. acknowledge support from the Elemental Strategy Initiative conducted by the MEXT, Japan. T.T. acknowledges support from a Grant-in-Aid for Scientific Research on Grant 262480621 and on Innovative Areas “Nano Informatics” (Grant 25106006) from JSPS. T.A.O. and K.C.F. acknowledge Raytheon BBN Technologies’ support for this work. This work was performed in part at the Center for Nanoscale Systems (CNS), a member of the National Nanotechnology Infrastructure Network (NNIN), which is supported by the National Science Foundation under NSF award no. ECS-0335765. CNS is part of Harvard University.