Wiedemann-Franz law for massless Dirac fermions with implications for graphene

In the 2016 experiment by Crossno et al. [Science 351, 1058 (2016)], electronic contribution to the thermal conductivity of graphene was found to violate the well-known Wiedemann-Franz (WF) law for metals. At liquid nitrogen temperatures, the thermal to electrical conductivity ratio of charge-neutral samples was more than 10 times higher than predicted by the WF law, what was attributed to interactions between particles leading to collective behavior described by hydrodynamics. Here we show, by adapting the handbook derivation of the WF law to the case of massless Dirac fermions, that significantly enhanced thermal conductivity should appear also in few- or even sub-kelvin temperatures, where the role of interactions can be neglected. The comparison with numerical results obtained within the Landauer-B\"uttiker formalism for rectangular and disk-shaped (Corbino) devices in ballistic graphene is also provided.

In the 2016 experiment by Crossno et al. [Science 351, 1058(2016)], electronic contribution to the thermal conductivity of graphene was found to violate the well-known Wiedemann-Franz (WF) law for metals.At liquid nitrogen temperatures, the thermal to electrical conductivity ratio of charge-neutral samples was more than 10 times higher than predicted by the WF law, what was attributed to interactions between particles leading to collective behavior described by hydrodynamics.Here we show, by adapting the handbook derivation of the WF law to the case of massless Dirac fermions, that significantly enhanced thermal conductivity should appear also in few-or even sub-kelvin temperatures, where the role of interactions can be neglected.The comparison with numerical results obtained within the Landauer-Büttiker formalism for rectangular and disk-shaped (Corbino) devices in ballistic graphene is also provided.

I. INTRODUCTION
Soon after the advent of graphene it become clear that this two-dimensional form of carbon shows exceptional thermal conductivity, reaching the room temperature value of ∼ 5000 W/m/K [1], being over 10 times higher than that of copper or silver [2].Although the dominant contribution to the thermal conductivity originates from lattice vibrations (phonons), particularly these corresponding to out-of-plane deformations [3,4] allowing graphene to outperform more rigid carbon nanotubes, the electronic contribution to the thermal conductivity (κ el ) was also found to be surprisingly high [5] in relation to the electrical conductivity (σ) close to the charge-neutrality point [6].One can show theoretically that the electronic contribution dominates the thermal transport at sub-kelvin temperatures [7], but direct comparison with the experiment is currently missing.Starting from a few kelvins, up to the temperatures of about T 80 K, it is possible to control the temperatures of electrons and lattice independently [5], since the electron-phonon coupling is weak, and to obtain the value of κ el directly.Some progress towards extending the technique onto sub-kelvin temperatures has been recently reported [8].
The Wiedemann-Franz (WF) law states that the ratio of κ el to σ is proportional to the absolute temperature [9] where the proportionality coefficient L is the Lorentz number.
For ideal Fermi gas, we have For metals, Eq. (1) with L ≈ L 0 (2) holds true as long as the energy of thermal excitations k B T ε F , with ε F being the Fermi energy.What is more, in typical metals close to the room temperature κ el κ ph , with κ ph being the phononic contribution to the thermal conductivity, and even when approximating the Lorentz number as L ≈ (κ el + κ ph )/σT one restores the value of L 0 (2) with a few-percent accuracy.
In graphene, the situation is far more complex, partly because κ el κ ph (starting from few Kelvins) but mainly be-cause unusual properties of Dirac fermions in this system.Experimental results of Ref. [5] show that the direct determination of κ el leads to L/L 0 = 10 − 20 for T = 50 − 75 K near the charge-neutrality point.Away from the charge-neutrality point, the value of L ≈ L 0 is gradually restored [10].Also, the Lorentz number is temperature-dependent, at a fixed carrier density, indicating the violation of the WF law.
High values of the Lorentz number (L/L 0 > 10) were observed much earlier for semiconductors [11], where the upper limit is determined by the energy gap (∆) to temperature ratio, L max ≈ (∆/2eT ) 2 , but for zero-gap systems strong deviations from the WF law are rather unexpected.Notable exceptions are quasi one-dimensional Luttinger liquids, for which L/L 0 > 10 4 was observed [12], and heavy-fermion metals showing L < L 0 [13].
The peak in the Lorentz number appearing at the charge neutrality point for relatively high temperatures (close to the nitrogen boiling point) can be understood within a hydrodynamic transport theory for graphene [14,15].However, it is worth to stress that for clean samples and much lower temperatures, where the ballistic transport prevails, one may still expect similar peaks with the maxima reaching L max /L 0 ≈ 2 − 3 and the temperature-dependent widths.
In this paper we show how to adapt the handbook derivation of the WF law [9] in order to describe the violation of this law due to peculiar dispersion relation and a bipolar nature of graphene.The quantitative comparison with the Landauer-Büttiker results is also presented, both for toy models of the transmission-energy dependence, for which closed-form formulas for L are derived, and for the exact transmission probabilities following from the mode-matching analysis for the rectangular [16][17][18] and for the disk-shaped [19,20] samples.
The remaining part of the paper is organised as follows.In Sec.II we recall the key points of the WF law derivation for ideal Fermi gas, showing how to adapt them for massless fermions in graphene.In Sec.III, the Landauer-Büttiker formalism is introduced, and the analytical results for simplified models for transmission-energy dependence are presented.The Lorentz numbers for mesoscopic graphene systems, the rectangle and the Corbino disk, are calculated in Sec.IV.The conclusions are given in Sec.V.
The derivation of the WF law for metals [9] starts from the relation between thermal conductivity of a gas with its heat capacity per unit volume (C) derived within kinetic theory of gases [21], which can be written as where d = 1, 2, 3 is the system dimensionality, v is a typical particle velocity, and is the mean-free path (travelled between collisions with boundaries or other particles).For the key points necessary to obtain Eq. ( 3), see Fig. 1.It is worth to notice that the definition of C in Eq. ( 3), used instead of a familiar specific heat (per unit mass), allows to generalize the reasoning onto the massless perticles easily.
Next, the electrical conductivity in Eq. ( 1) is expressed via the Drude formula where n = N/V is the carrier density (to be redefined later for a bipolar system containing electrons and holes), and m is the carrier effective mass.We skip here the detailed derivation of Eq. ( 4), which can be found in Ref. [9]; we only mention that it follows from Ohm's law in the form j = σE, with j the current density and E the electric field, supposing that carriers of the ±e charge and the m mass accelerate freely during the time τ = /v [with the symbols and v same as in Eq. ( 3)].This time, a generalization for massless particles is more cumbersome; we revisit this issue in Sec.II C.
The system volume, referred in definitions of C and n, can be denoted as V = L d , with L being linear dimension of a box of gas.In the SI units, the dimension of C is J/(m d •K), and the unit of thermal conductivity is Similarly, the unit of electrical conductivity is In turn, the unit of length (m) vanishes in the κ/σ ratio occurring in Eq. (1) and the WF law remains valid for arbitrary d (provided that the suppositions given explicitly in Sec.II B are satisfied.)Unfortunately, in the literature on graphene σ is commonly specified in Ω −1 (≡ S), as follows from Eq. ( 6) for d = 2, but the values of κ are reported in W/m/K, as for d = 3 [2].Such an inconsistency can be attributed to the fact that for the thermal conductivity of multilayer graphenes linear scaling with the number of layers remains a reasonable approximation [22], yet the behavior of electrical conductivity is far more complex [23,24] even for bilayers [25].
< l a t e x i t s h a 1 _ b a s e 6 4 = " h L N 8 m e 7 q 6 e r T v 5 Figure 1: Relation between the thermal conductivity (κ), heat capacity per unit volume (C), average particle velocity (v), and the mean-free path ( ).The non-equilibrium heat flow occurs between the interfaces (1) and ( 2), with local temperatures T and T + ∆T , separated by a distance ∆x = vxτ (with vx the mean velocity in x direction and τ the relaxation time) and can be quantified by ∆Q = CS∆x∆T .The corresponding thermal conductivity is κ = ∆Q (Sτ ∆T /∆x) −1 = Cv 2 x τ .Substituting v 2 x = v 2 /d, and = vτ , we obtain Eq. (3) in the main text.

B. The Fermi gas in metals
The calculation of C in Eq. ( 3) employs the free Fermi gas approximation for electrons in a metal.In this approximation, one assumes that leading contributions to thermodynamic properties originate from a thin layer around the Fermi surface.For instance, a contribution to the internal energy can be written as where ε F is the Fermi energy, 2Λ is the relevant energy interval considered (ε F Λ k B T ), D(ε) is the density of states per unit volume (i.e., the number of energy levels lying in the interval of ε, . . ., ε + dε is V D(ε)dε ), and f (ε) is the Fermi-Dirac distribution function In a general case, the chemical potential in Eq. ( 8) is adjusted such that the particle density take a desired value n(µ) ≡ n, defining the temperaturedependent chemical potential µ = µ(T ).Here, the constantdensity of states approximation, D(ε) ≈ D(ε F ) for F − Λ F + Λ imposed in the rightmost expression in Eq. (7), is equivalent to µ ≈ ε F [26].
Definite integral in Eq. ( 7) is equal to where the Riemann zeta function is introduced to be used in forthcoming expressions.Differentiating ∆U el (7) over temperature, one gets approximating expression for the electronic heat capacity In fact, the factor of π 2 /3 in Eq. ( 12) is the same as appearing in the Lorentz number L 0 (2), what is shown in a few remaining steps below.For an isotropic system with parabolic dispersion relation bounded in a box of the volume V = L d with periodic boundary conditions, the wavevector components k = (k i ) take discrete values of k i = 0, ± 2π L , ± 4π L , . . .(with i = x, y, z for d = 3).Calculation of the density of states in d = 1, 2, 3 dimensions is presented in numerous handbooks [27]; here, we use a compact form referring to the particle density on the Fermi level where n(ε F ) = ε F 0 D(ε)dε representing the T → 0 limit of Eq. ( 9).Substituting D(ε F ), given by Eq. ( 14), into Eq.( 12) we obtain Now, taking F = 1 2 m v 2 F with the Fermi velocity and the Fermi wavevector k F = 2m F / 2 , we further set v = v F in Eq. (3), obtaining It is now sufficient to divide Eqs. ( 17) and ( 4) side-by-side to derive the WF law as given by Eqs.
(1) and (2).As mentioned earlier, the result for free Fermi gas is same for arbitrary dimensionality d.More careful analysis also shows that the parabolic dispersion of Eq. ( 13) is not crucial, provided that the Fermi surface is well-defined, with an (approximately) constant D(ε) > 0 in the vicinity of |ε − ε F | k B T , and that the effective mass 0 < m < +∞.In the framework of Landau's Fermi-liquid (FL) theory, the reasoning can be extended onto effective quasiparticles, and the validity of the WF law is often considered as a hallmark of the FL behavior [28,29].
The suppositions listed above are clearly not satisfied in graphene close to the charge-neutrality point.

C. The Dirac gas in graphene
The relation between thermal conductivity and heat capacity given by Eq. (3) holds true for both massive and massless particles.A separate issue concerns the Drude formula (4), directly referring to the effective mass, an adaptation of which for massless Dirac fermions requires some attention.
The Landauer-Büttiker conductivity of ballistic graphene, first calculated analytically employing a basic mode-matching technique [16][17][18] and then confirmed in several experiments [30,31], is given solely by fundamental constants Remarkably, for charge-neutral graphene both the carrier concentration and the effective mass vanish; a finite (and nonzero) value of σ 0 (18) may therefore be in accord with the Drude formula, at least in principle.
In order to understand the above conjecture, we refer to the approximate dispersion relation for charge carriers in graphene, showing up so-called Dirac cones, The value of the Fermi velocity v F ≈ 10 6 m/s is now energyindependent, being determined by the nearest-neighbor hopping integral on a honeycomb lattice (t 0 = 2.7 eV) and the lattice constant (a = 0.246 nm) via Charge carriers in graphene are characterized by an additional (next to spin) quantum number, the so-called valley index.This leads to an additional twofold degeneracy of energy levels, which needs to taken into account when calculating the density of states, Subsequently, the carrier concentration at T = 0 is related to the Fermi energy (and the Fermi wavevector) via In the above we intentionally omitted the F index for symbols denoting the Fermi energy and the Fermi wavevector to emphasize that they can be tuned (together with the concentration) by electrostatic gates, while the Fermi velocity v F (20) is a material constant [32].Despite the unusual dispersion relation, given by Eq. ( 19), the relevant effective mass describing the carrier dynamics in graphene is the familiar cyclotronic mass where A(ε) denotes the area in momentum space (k x , k y ) bounded by the equienergy surface for a given Fermi energy (ε).It is easy to see that for two-dimensional system, with fourfold degeneracy of states, we have ∂A(ε)/∂ε = π 2 D(ε); substituting D(ε) given by Eq. ( 21) leads to the rightmost equality in Eq. (23).Remarkably, the final result is formally identical with the rightmost equality in Eq. ( 16) for free Fermi gas (albeit now the effective mass, but not the Fermi velocity, depends on the Fermi energy).
Assuming the above carrier density n (22), and the effective mass m (23), and comparing the universal conductivity σ 0 (18) with the Drude formula (4), we immediately arrive to the conclusion that mean-free path for charge carriers in graphene is also energy-dependent, taking the asymptotic form Strictly speaking, in the ε → 0 limit is have n → 0, i.e., no free charge carriers, and the transport is governed by evanescent waves [6].The universal value of σ 0 (18) indicates a peculiar version of the tunneling effect appearing in graphene, in which the wavefunction shows a power-law rather then exponential decay with the distance [19], resulting in the enhanced charge (or energy) transport characteristics.Therefore, the mean-free path should be regarded as an effective quantity, allowing one to reproduce the measurable characteristics in the ε → 0 limit.Away from the charge-neutrality point, i.e., for |ε| π v F /L (with the geometric energy quantization ∼ π v F /L), graphene behaves as a typical ballistic conductor, with eff ∼ L. We revisit this issue in Sec.IV, where the analysis starts from actual σ(ε) functions for selected mesoscopic systems, but now the approximation given by Eq. ( 24) is considered as a first.
We further notice that the form of eff (ε) in Eq. ( 24) is formally equivalent to the assumption of linear relaxation time on energy dependence in the Boltzmann equation, proposed by Yoshino and Murata [33].
In the remaining part of this section, we derive explicit forms of the thermal conductivity κ and the Lorentz number L, pointing out the key differences appearing in comparison to the free Fermi gas case (see Sec. II B).
The calculations are particularly simple for charge-neutral graphene (n = ε = 0), which is presented first.Although we still can put v = v F in Eq. ( 3), since the Fermi velocity is energy-independent, the constant-density of states approximation applied in Eq. ( 7) in now invalid.(Also, for T > 0 we cannot put ε F k B T now.)In turn, the expression for heat capacity C needs to be re-derived.
For charge-neutral graphene at T > 0, contributions from thermally excited electron and holes are identical, it is therefore sufficient to calculate the former Again, the integral in the rightmost expression in Eq. ( 25) can be expressed via the Riemann zeta function, and is equal to Differentiating Eq. ( 25) with respect to T , and multiplying by a factor of 2 due to the contribution from holes in the valence band, we obtain the heat capacity It remains now to calculate the effective mean-free path to be substituted to Eq. (3).We use here the asymptotic form of eff (ε) (24), replacing the ε −1 factor by its overage over the grand canonical ensemble, namely Substituting the above, together with the heat capacity C (27) into Eq.( 3), we get and with L 0 being the Fermi-gas result given by Eq. ( 2).A simple reasoning, presented above, indicates that the κ/σ ratio is significantly enhanced in charge-neutral graphene, comparing to the free Fermi gas.However, the WF law is still satisfied, since the Lorentz number given by Eq. ( 30) is temperature-independent.The situation becomes remarkably different for graphene away from the charge-neutrality point, which is studied next.
Without loss of generality, we suppose µ > 0 (the particle hole-symmetry guarantees that measurable quantities are invariant upon µ → − µ).The internal energy U (T ) now consists of contributions from majority carries (electrons), with ε > µ, and minority carriers (holes), with ε < µ, where D(ε) is given by Eq. ( 21).The heat capacity can be written as where we have defined with y = µ/k B T > 0 and Li s (z) being the polylogarithm function [34].
Similarly, the mean-free path can be calculated as where and y = µ/k B T again.Hence, the Lorentz number for µ > 0 is given by with F (y) and G(y) given by Eqs.(33) and (35).The Lorentz number given by Eq. ( 36) is depicted in Fig. 2. It is straightforward to show that in the y → 0 limit one obtains the value given by Eq. ( 30) for µ = 0; also, for y → ∞ we have L → L 0 , restoring a standard form of the WF law for metals.However, for 0 < y < +∞, a fixed value of µ (or n) corresponds to y (and thus L) varying with temperature; namely, the violation of the WL law occurs.

III. LANDAUER BÜTTIKER FORMALISM AND SIMPLIFIED MODELS
A. The formalism essential In the Landauer-Büttiker description transport properties of a mesoscopic system, attached to the leads, are derived from the transmission-energy dependence T (ε), to be found by solving the scattering problem [35][36][37][38].In particular, the Lorentz number can be written as [39] where L n (with n = 0, 1, 2) are given by with g s = g v = 2 denoting spin and valley degeneracies in graphene, and the Fermi-Dirac distribution function f (ε) given by Eq. ( 8).It is easy to show that energy-independent transmission ( T (ε) = const ) leads to L = L 0 (2).

B. Simplified models
Before calculating T (ε) directly for selected systems in Sec.IV, we first discuss basic consequences of some model T (ε) functions for L.
For instance, the linear transmission-energy dependence (i.e., T (ε) ∝ |ε| ), allows one to obtain a relatively short formula for L at arbitrary doping [7], namely with y = µ/k B T .For y = 0, the Lorentz number given by Eq. ( 39) takes the value of being close to that given in Eq. ( 30).The approximation given in Eq. ( 40) was earlier put forward in the context of high-temperature supeconductors also showing the linear transmission-energy dependence [40].Numerical values of L(y) are presented in Fig. 2. Remarkably, L(y) obtained from Eq. (36) [blue line] is typically 20 − 30% higher than obtained Eq. (39) [red line].The deviations are stronger near |µ|/k B T ≈ 4.5, where the latter shows broad minima absent for the former.Above this value, L(y) obtained from Eq. ( 36) approaches L 0 from the top, whereas L(y) obtained from Eq. ( 39) approaches L 0 from the bottom.Also, the right-hand side of Eq. ( 36) converges much faster to L 0 for |µ| k B T than the right-hand side of Eq. ( 39).
< l a t e x i t s h a 1 _ b a s e 6 4 = " L F p t w y u 2 2 p 4 y q c e w 8 X J T s 5 5 U 9 w 0 Figure 3: The Lorentz number for model transmission-energy dependence T (ε) given by Eq. ( 41) with m varied from −0.5 to 2.5 with the steps of 0.5 displayed as a function of the chemical potential.Solid (dashes) lines mark integer (non-integer) m.
In both cases, the Lorentz number enhancement at the charge-neutrality point (µ = 0) is significant, and the violations of the WF law for µ = 0 is apparent.A relatively good agreement between the two formulas is striking: Although both derivations have utilized the linear dispersion of the Dirac cones, being link to D(ε) given by Eq. ( 21) in the first case, or to the T (ε) ∝ |ε| assumption in the second case (see Sec. IV for further explanation), but only the derivation of Eq. (36) incorporates the information about the universal conductivity (σ = σ 0 ).We can therefore argue that the L enhancement occurs in graphene due to the linear dispersion rather then due the transport via evanescent waves (being responsible for σ = σ 0 at µ = 0).
We now elaborate possible effects, on the Lorentz number, of toy-models of transmission-energy dependence where the proportionality coefficient is irrelevant due to the structure of Eq. (37).For some cases, integrals can be calculated analytically, leading e.g. to L = L 0 for m = 0 (the constant transmission case), or to L = L(y) given by Eq. ( 39) for m = 1 (the linear transmission-energy dependence).Numerical results for selected values of m = −0.5 . . .2.5 are displayed in Fig. 3.
The violation of the WF law appears generically for m = 0 away from the charge-neutrality point (i.e., for µ = 0).
For µ = 0, the Lorentz number reaches a global maximum (with L > L 0 ) if m > 0, or a global minimum (with L < L 0 ) if −1 < m < 0. A close-form expression can be derived for L(µ = 0)/L 0 < l a t e x i t s h a 1 _ b a s e 6 4 = " u f d e + J l r 6 f G r y j M l a r 0 < l a t e x i t s h a 1 _ b a s e 6 4 = " t 9 5 < l a t e x i t s h a 1 _ b a s e 6 4 = " t 9 5   42).Insets visualize the T (ε) function given by Eq. ( 41) for m = 0 and m = 1, with contributions from the valence band (p) and the conduction band (n).
both the cases, namely and is visualized in Fig. 4 [41].It is clear that T (ε) models given by Eq. ( 41) may lead to arbitrarily high L max ; in particular, the value of 10 L 0 is exceeded starting from m ≈ 4.1.
Hence, for m > 1, the model grasps the basic features of one-dimensional Luttinger liquids, showing both the power-law transmission energy dependence, with nonuniversal (interaction dependent) exponents, and the significantly enhanced Lorentz numbers [12].
On the other hand, the suppression of L is observed for −1 < m < 0, due to the integrable singularity at ε = 0, constituting an analogy with heavy fermion systems [13].
Both the above-mentioned scenarios were described theoretically for quantum dot systems, which may be tuned from the suppression of L due to Breit-Wigner resonance, to the enhancement of L due to Fano resonance [42][43][44].

C. Gapped systems
For a sake of completeness, we show here how the energy (or transport) gap may enhance the Lorentz number.Instead of T (ε) given by Eq. ( 41), we put < l a t e x i t s h a 1 _ b a s e 6 4 = " < l a t e x i t s h a 1 _ b a s e 6 4 = " The Lorentz number as a function of chemical potential for the limit of an unipolar system, corresponding to T (ε) given by Eq. ( 43) with ∆ kBT and µ ≈ ∆/2 [see also Eq. ( 45)].The exponent m is varied from 0 to 3 with the steps of 1 (solid lines).Dashed line marks the Wiedemann-Franz value (L = L0).
where Θ(x) is the Heaviside step function.
An apparent feature of Eq. ( 44) is that L max shows an unbounded growth with a gap (with the leading term being of the order of ∼ ∆ 2 ), in agreement with the experimental results for semiconductors [11].Similar behaviors can be expected for tunable-gap systems, such as bilayer graphene or silicene, which are beyond the scope of this work.
A different behavior appears near the band boundary, i.e., for µ ≈ ∆/2 (or µ ≈ − ∆/2).Assuming ∆ k B T again, we arrive to the limit of an unipolar system, for which only the contribution from majority carries to integrals L n (38) matters.In effect, the Lorentz number can be approximated as where y = (µ − 1 2 ∆)/k B T and Closed-form expressions for J n (m; y) are not available; a few numerical examples for m = 0 . . . 3 are presented in Fig. 5. Since now L 1 ∝ J 1 = 0 (in contrast to the bipolar case studied before), the Lorentz number is significantly reduced, and relatively close to L 0 , which is approached for y 1. Asymptotic forms of J n (m; y) can be derived for |y| 1, namely where Γ(z) denotes the Euler gamma function, and Substituting the above into Eq.( 45), we obtain or Both limits are closely approached by the numerical data in Fig. 5 for |y| 5.In all the cases considered, the values of L are now much lower than the corresponding L max for a gapless model with the same m (see Fig. 4).
Therefore, it becomes clear from analyzing simplified models of T (ε) that a bipolar nature of the system, next to the monotonically-increasing transmission (the m > 0 case) are essential when one looks for a significant enhancement of the Lorentz number L (compared to L 0 ).
Both these conditions are satisfied for graphene.

IV. EXACTLY SOLVABLE MESOSCOPIC SYSTEMS A. Transmission-energy dependence
The exact transmission-energy dependence T (ε) can be given for two special device geometries in graphene: a rectangular sample attached to heavily-doped graphene leads [16][17][18] and for the Corbino disk [19,20].Although these systems posses peculiar symmetries, allowing one to solve the scattering problem employing analytical mode-matching method (in particular, the mode mixing does not occur), both the solutions were proven to be robust against various symmetrybreaking perturbations [46][47][48].More importantly, several features of the results have been confirmed in the experiments [30,31,50,51] showing that even such idealized systems provide valuable insights into the quantum transport phenomena involving Dirac fermions in graphene.
For a rectangle of width W and length L, the transmission can be written as [17,19] where the transmission probability for n-th normal mode is given by with q n = π(n + 1 2 )/W the quantized transverse wavevector (the constant 1 2 corresponds to infinite-mass confinement; for other boundary conditions, see Ref. [17]), and k = |ε|/( v F ).The two cases in Eq. ( 53) refer to the contributions from propagating waves (k q n , so-called open channels) and evanescent waves (k < q n ).
For the Corbino disk, with its inner (R 1 ) and outer (R 2 ) radii, we have [19] where j the the half-odd integer angular momentum quantum number, with a corresponding transmission probability where k is same as in Eq. ( 53), and (1,2) ν (ρ) the Hankel function of the (first, second) kind.

B. The conductivity
A measurable quantity that provides a direct insight into the T (ε) function is zero-temperature conductivity with the conductance quantum g 0 = 4e 2 /h and a shapedependent factor  For T > 0, Eq. ( 57) needs to be replaced by σ(µ) = e 2 Ω X L 0 , where L 0 is given by Eq. ( 38) with n = 0. Numerical results, for T = 0, are presented in Fig. 6.The data for both systems, displayed versus a dimensionless quantity εL/ v F (with L ≡ R 2 − R 1 for a disk) closely follow each other up to |ε|L/ v F ≈ 3.For larger values of |ε|, the results become shape-dependent and can be approximated, for |ε| v F /L, as  6] match the approximation given by Eq. ( 59) with a few percent accuracy for |ε| 5 v F /L, whereas for a disk with R 2 /R 1 = 2 [dashed red line] a systematic offset of ≈ (1/π)g 0 occurs, signaling an emphasized role of evanescent waves in the Corbino geometry.This observation coincides with a total lack of Fabry-Perrot oscillations in the Corbino case.

C. The Lorentz number
The exact transmission-energy functions T (ε), discussed above, are now substituted to Eq. ( 37) for the Lorentz number.Calculating the relevant integrals numerically, we obtain the results presented in Figs.7 and 8.
Close to the charge-neutrality point, i.e., for |µ| max( v F L −1 , k B T ), both systems show a gradual crossover < l a t e x i t s h a 1 _ b a s e 6 4 = " / D y s < l a t e x i t s h a 1 _ b a s e 6 4 = " e l G 5       (with increasing T ) from the Wiedemann-Franz regime, with a flat L ≈ L 0 , to the linear-transmission regime characterized by L(µ) close to the predicted by Eq. ( 39) [see Fig. 7].For higher µ, some aperiodic oscillations of L(µ) are visible if k B T v F /L, being particularly well pronounced for a rectangular sample.For higher temperatures, the oscilltions are smeared out, leaving only one shallow minimum near |µ|/k B T ≈ 4 − 5, in agreement with Eq. (39).
Maximal values of L for the two systems (reached at µ = 0) are displayed, as functions of temperature, in Fig. 8.It is clear that a crossover between low and high temperature regimes takes place near k B T ∼ v F /L (corresponding to ≈ 6.67 K for L = 1 µm): For lower temperatures (and near µ = 0), thermally-excited carriers appear in the area where T (ε) ≈ const (leading to L ≈ L 0 ), whereas for significantly higher temperatures, the detailed behavior of T (ε) near ε = 0 becomes irrelevant, and the linear-transmission approximation ( T (ε) ∝ |ε| ) applies.Remarkably, the convergence to the value given in Eq. ( 40) is much slower (yet clearly visible) in the Corbino disk case, due to a higher (compared to a rectangular sample) contribution from evanescent waves to the transmission away from the charge-neutrality point.

V. CONCLUSIONS
We have calculated the Lorentz number (L = κ el /σT ) for noninteracting massless Dirac fermions following two different analytic approaches: first, adapting the handbook derivation of the Wiedemann-Franz (WF) law, starting from the relation between thermal conductivity and heat capacity obtained within the kinetic theory of gases, and second, involving the Landauer-Büttiker formalism and postulating simple model of transmission-energy dependence, T (ε) ∝ |ε|.In both approaches, the information about conical dispersion relation is utilized, but the universal value of electrical conductivity, σ ∼ e 2 /h at ε = 0, is referred only in the first approach.Nevertheless, the results are numerically close, indicating the violation of the WF law with maximal Lorentz numbers L max /L 0 ≈ 2.77 and 2.37 (respectively) and L → L 0 = (π 2 /3) k 2 B /e 2 for high doppings (|ε| k B T ).This observation suggests that violation of the WF law, with L max /L 0 ≈ 2 − 3 should appear generically in weakly-doped systems with approximately conical dispersion relation, including multilayers and hybrid structures, even when low-energy details of the band structure alter the conductivity.
Moreover, a generalized model of power law transmissionenergy dependence, T (ε) ∝ |ε| m (with m > −1), is investigated in order to address the question whether the enhancement of L is due to the bipolar band structure or due to the conical dispersion.Since L > L 0 shows up for any m > 0, and the maximal value grows monotonically with m, we conclude that the dispersion relation has a quantitative impact on the effect.On the other hand, analogous discussion of gapped systems, with the chemical potential close to the center of the gap (the bipolar case) or to the bottom of the conduction band (the unipolar case) proves that the bipolar band structure is also important (no enhancement of L is observed in the unipolar case up to m ≈ 2).
Finally, the Lorentz numbers, for different dopings and temperatures, are elaborated numerically from exact solutions available for the rectangular sample and the Corbino (edge-free) disk in graphene, both connected to heavily-doped graphene leads.The results show that L, as a function of the chemical potential µ, gradually evolves (with growing T ) as expected for a model transmission energy dependence, T (ε) ∝ |ε| m , with the exponent varying from m = 0 to m = 1.The upper bound is approached faster for the rectangular sample case, but in both cases L/L 0 > 2 is predicted to appear for T 13 K • µm × L −1 with L the sample length.
Our results complement earlier theoretical study on the topic [33] by including the finite size-effects and the inter-  play between propagating and evanescent waves, leading to the results dependent, albeit weakly, on the sample geometry.

Figure 2 :
Figure 2: The Lorentz number L = κ el /(σT ) for massless Dirac fermions as a function of the chemical potential.Solid lines represent the approximations given by Eq. (36) [blue line] and Eq.(39) [red line].Dashed lines (top to bottom) depict the two corresponding µ = 0 values, and the value of L0 = (π 2 /3) k 2B /e 2 representing the Wiedemann-Franz law restored in the |µ| kBT limit.
t e x i t s h a 1 _ b a s e 6 4 = " Q 8 W 4 z P 7

y 4 r 2 A
e 0 p W T S T B u b S Y b k T q G U / o M b E T c K f o 2 / 4 N + Y t r N p 6 4 H A 4 Z w T 7 j 0 3 T K S w 6 P u / X m 5 r e 2 d 3 L 7 9 f O D g 8 O j 4 p n p 4 1 r E 4 N 4 3 W m p T a t k F o u h e J 1 F C h 5 K z G c x q H k z X D 0 M P e b Y 2 6 s 0

a f 8 O
W 9 e K / e u / e x j O a 8 7 M 8 5 r M D 7 / g N + c I 0 c < / l a t e x i t > " m = 1 < l a t e x i t s h a 1 _ b a s e 6 4 = " 7 4 C 5 R w

y 4 r 2 A
e 0 p W T S T B u b S Y b k T q G U / o M b E T c K f o 2 / 4 N + Y t r N p 6 4 H A 4 Z w T 7 j 0 3 T K S w 6 P u / X m 5 r e 2 d 3 L 7 9 f O D g 8 O j 4 p n p 4 1 r E 4 N 4 3 W m p T a t k F o u h e J 1 F C h 5 K z G c x q H k z X D 0 M P e b Y 2 6 s 0

a f 8 OFigure 4 :
Figure 4: Maximal (solid blue line for m > 0) or minimal (dashed red line for −1 < m < 0) values of the Lorentz number L (reached at µ = 0) obtained from Eq. (42).Insets visualize the T (ε) function given by Eq. (41) for m = 0 and m = 1, with contributions from the valence band (p) and the conduction band (n).
0 5 X t 5 e L / i X N c s m u l K z 7 i 0 K 1 P B 8 h S 4 7 I M S k S m 1 y S K r k j N d I g j D y T V / J J v o w n 4 8 V 4 M 9 5 n 0 Y w x v 3 N I F m B 8 / w K 6 y Z 7 7 < / l a t e x i t > t e x i t s h a 1 _ b a s e 6 4 = " E b o v 7 c q t 9 Z b I n X L A w 0 b z K / u h 7 z c = " > A A A B 9 n i c b V D L S g M x F L 1 T X 7 W + R l 2 4 c J N a B F d l p l R 0 I x T c u H B R w T 6 g U 4 Z M m m l D M w + S j F i G + R U 3 I m 4 U / A x / w b 8 x b c d F W w 8 E T s 4 5 I f c e L + Z M K s v 6 M Q p r 6 x u b W 8 X t 0 s 7 u 3 v 6 B e X j U l l t e x i t s h a 1 _ b a s e 6 4 = " 9 C K C p 0 V 0 j M K p T G P 9 D c e m E 4 1 2 S P o = " > A A A B 4 n i c b

Figure 6 :
Figure 6: Zero-temperature conductivity as a function of the Fermi energy for a rectangular sample with width-to-length ratio W/L = 5 (solid blue line) and the Corbino disk with radii ratio R2/R1 = 2 (dashed red line).Both system are shown schematically.Dashed black line marks the universal conductivity σ0 = (4/π) e 2 /h.
where the number of open channels N open (ε) = kW/π , for rectangle, 2 kR 1 , for disk, (60) with x being the floor function of x, and the average transmission per open channel T open ≈ π/4 < 1 (for the derivation, see Appendix A).Remarkably, numerical values of σ(ε) for a rectangle with W/L = 5 [solid blue line in Fig.

2 < l a t e x i t s h a 1 _ b a s e 6 4 =
9 o S 8 m k m T Y 0 8 y C 5 I 5 T S v R s R N w p + k b / g 3 5 h O Z 9 P W A 4 H D O S f c e 6 6 f K G m Q 0 l + n s L W 9 s 7 t X 3 C 8 d H B 4 d n 5 R P z 1 om T j U X T R 6 r W H d 8 Z o S S k W i i R C U 6 i R Y s 9 J V o + 5 P 7 h d 9 + E d r I O H r G a S L 6 I R t F M p C c o Z W e q F s b l C v U p R n I J v F y U o E c j U H 5 p z e M e R q K C L l i x n Q 9 m m B / x j R K r s S 8 1 E u N S B i f s J G Y Z Q v O y Z W V h i S I t X 0 R k k x d y b H Q m G n o 2 2 T I c G z W v Y X 4 n 9 d N M b j r z 2 S U p C g i v h w U p I p g T B Z t y V B q w V F N L W F c S 7 s h 4 W O m G U d 7 k 5 K t 7 q 0 X 3 S S t m u v d u P T x u l K v 5 k c o w g V c Q h U 8 u I U 6 P E A D m s B h B G / w C V + O 7 7 w 6 7 8 7 H M l p w 8 j / n s A L n + w 9 C h 4 i 9 < / l a t e x i t > 0." v P Q g F / u y C / h U h T 4 y 0 t V N i n M O A 2 0 = " > A A A B 4 H i c b V D L S g M x F L 1 T X 7 W + q i 7 d B I v Q 1 Z A R i y 4 L b l x W t A 9 o S 8 m k m T Y 0 8 y C 5 I 5 T S v R s R N w p + k b / g 3 5 h O Z 9 P W A 4 H D O S f c e 6 6 f K G m Q 0 l + n s L W 9 s 7 t X 3 C 8 d H B 4 d n 5 R P z 1 o m T j U X T R 6 r W H d 8 Z o S S k W i i R C U 6 i R Y s 9 J V o + 5 P 7 h d 9 + E d r I O H r G a S L 6 I R t F M p C c o Z W e q F s b l C v U p R n I J v F y U o E c j U H 5 p z e M e R q K C L l i x n Q 9 m m B / x j R K r s S 8 1 E u N S B i f s J G Y Z Q v O y Z W V h i S I t X 0 R k k x d y b H Q m G n o 2 2 T I c G z W v Y X 4 n 9 d N M b j r z 2 S U p C g i v h w U p I p g T B Z t y V B q w V F N L W F c S 7 s h 4 W O m G U d 7 k 5 K t 7 q 0 X 3 S St a 9 e r u f T x p l K v 5 k c o w g V c Q h U 8 u I U 6 P E A D m s B h B G / w C V + O 7 7 w 6 7 8 7 H M l p w 8 j / n s A L n + w 9 G 9 Y j A < / l a t

Figure 7 :
Figure7: The Lorentz number for a rectangular sample (left) and the Corbino disk (right) displayed as a function of the chemical potential.The temperature, specified in the units of vF /(kBL) ≈ 6.67 K • µm × L −1 , is varied between the lines and same in both panels.Remaining parameters are same as in Fig.6.
l a t e x i t > k B T L/~v F < l a t e x i t s h a 1 _ b a s e 6 4 = " G d r V c s 7 w C R 5 N E r S G u S r u P c L r 8 1 s = " > A A A B / H i c b V D L S s N A F L 3 x W e s r 6 l K E w S J 0V R N R d F l w 4 8 J F B f u A p o T J d N o O n U n C z E Q s I W 7 8 F T c i b h T 8 B n / B v 3 H a B q G t B w Y O 5 5 z h 3 n u C m D O l H e f H W l p e W V 1 b L 2 w U N 7 e 2 d 3 b t v f 2 G i h J J a J 1 E P J K t A C v K W U j r m m l O W 7 G k W A S c N o P h 9 d h v P l C p W B T e 6 1 F M O w L 3 Q 9 Z j B G s j + f Z R 6 h H M 0 W 3 m p 5 4 U S O D H 7 P R P c n y 7 5 F S c C d A i c X N S g h w 1 3 / 7 2 u h F J B A 0 1 4 V i p t u v E u p N i q R n h N C t 6 i a I x J k P c p + l k + Q y d G K m L e p E 0 L 9 R o o s 7 k s F B q J A K T F F g P 1 L w 3 F v / z 2 o n u X X V S F s a J p i G Z D u o l H O k I j Z t A X S Y p 0 X x k C C a S m Q 0 R G W C J i T Z 9 F c 3 p 7 v y h i 6 R x V n E v K s 7 d e a l a z k s o w C E c Q x l c u I Q q 3 E A N 6 k D g G V 7 h A z 6 t J + v Fe r P e p 9 E l K / 9 z A D O w v n 4 B 4 b e T + A = = < / l a t e x i t > L max /L 0 < l a t e x i t s h a 1 _ b a s e 6 4 = " 3 s O Z U 5 X c 7 b u x w D 0 a o / G c U 2 9 3 9 V o = " > A A A B 8 X i c b V D L S s N A F L 2 p r 1 p f s S 7 d D B a h q 5 q I o s u i I C 5 c V O g L 2 h I m 0 0 k 7 d P J g Z l I s I R / i R s S N g j / i L / g 3 T t p s 2 n p g 4 H D O G e 4 9 1 4 0 4 k 8 q y f o 3 C x u b W 9 k 5 x t 7 S 3 f 3 B 4 Z B 6 X 2 z K M B a E t E v J Q d F 0 s K W c B b S m m O O 1 G g m L f 5 b T j T u 4 z v z O l Q r I w a K p Z R A c + H g X M Y w Q r L T l m e e L c J W n z 6 a I / d r F I 0 q n z 4 J g V q 2 b N g d a J n Z M K 5 G g 4 5 k 9 / G J L Y p 4 E i H E v Z s 6 1 I D R I s F C O c p q V + L G m E y Q S P a D L f O E X n W h o i L x T 6 B Q r N 1 a U c 9 q W c + a 5 O + l i N 5 a q X i f 9 5 v V h 5 t 4 O E B V G s a E A W g 7 y Y I x W i r D 4 a M k G J 4 j N N M B F M b 4 j I G A t M l D 5 S S V e 3 V 4 u u k / Z l z b 6 u W c 9 X l X o 1 P 0

Figure 8 :
Figure 8: Maximal Lorentz number (corresponding to µ = 0) for same systems as in Fig. 6 versus temperature.Inset shows the data replotted from main panel with the abscissa scaled logarithmically.Dashed horizontal line marks the prediction given in Eq. (40).
t e x i t s h a 1 _ b a s e 6 4 = "

3 < 4 Figure 9 :
Figure9: Fermi energy dependence of the Fano factor (main panel) and the average transmission per channel (inset), defined in Eq. (A4), for same systems as in Fig.6.