A 2007 finite-size Kubo-formula study of two-dimensional massless Dirac fermions (the low-energy model of graphene) with screened Coulomb disorder found a finite minimum conductivity σmin⁡∼e2/h\sigma^{\min}\sim e^2/h at the Dirac point, and identified as open questions the ultimate low-temperature, infinite-size transport state: whether symplectic-class weak antilocalization drives the conductivity to diverge, whether a finite universal value holds, or whether localization sets in once intervalley scattering becomes relevant, together with the quantitative role of electron-electron screening beyond Thomas-Fermi theory. Later work resolved much of the practical phenomenology (nonuniversal, impurity-density-dependent σmin⁡\sigma^{\min}, percolating electron-hole puddles, quantum-mechanical screening, weak-localization crossovers), but the strict T→0T\to 0, L→∞L\to\infty classification of Dirac-point transport and a nonperturbative value of the ground-state conductivity remain unsettled, with active numerical work as recent as 2026.

Background

Graphene's low-energy carriers are described by a two-dimensional massless Dirac fermion (MDF) model: the honeycomb band structure has conduction and valence bands meeting linearly at two inequivalent valleys (KK and K′K'), so that the Hamiltonian is H=ℏv σ⋅kH=\hbar v\,\boldsymbol{\sigma}\cdot\mathbf{k} with vv the Fermi velocity and σ\boldsymbol{\sigma} Pauli matrices acting on the sublattice degree of freedom. At the charge-neutrality (Dirac) point the clean density of states vanishes, yet experiments find that graphene's conductivity never falls below a minimum value σmin⁡\sigma^{\min} of order 4e2/h4e^2/h per layer, larger than the value (1/π)(4e2/h)(1/\pi)(4e^2/h) predicted by the self-consistent Born approximation (SCBA, a self-consistent lowest-order perturbative treatment of disorder) for short-range scatterers. Nomura and MacDonald [Phys. Rev. Lett. 98, 076602 (2007)] evaluated the Kubo formula (the quantum-mechanical linear-response expression for conductivity) numerically for finite systems with either short-range or screened Coulomb disorder. They found that Coulomb scatterers, which do not scatter between the two valleys because their potential is smooth on the lattice scale, give σmin⁡∼e2/h\sigma^{\min}\sim e^2/h per channel and a conductivity that increases only very slowly with system size LL at the Dirac point.

Three notions are needed to state what they left open. First, smooth disorder that suppresses intervalley (valley-changing) scattering places the system in the symplectic symmetry class of Anderson localization theory; in two dimensions this class shows weak antilocalization — a quantum interference effect, driven by the Berry phase accumulated along the honeycomb-lattice momentum-space paths, that makes the conductivity grow logarithmically with system size, so that scaling arguments suggest an ultimately divergent ('perfect metal') conductivity — whereas with intervalley scattering the system crosses to the orthogonal class, in which all states localize at sufficiently large sizes. The long-wavelength field theory of this quantum-interference physics is the nonlinear sigma model (the standard effective field theory of diffusive, noninteracting disordered electrons); in the symplectic case it acquires a topological term with angle θ=π\theta=\pi, a Berry-phase consequence of the honeycomb-lattice band structure, whose presence underlies the competing prediction of a quantum critical point with a finite universal conductivity. Second, the disorder potential of charged impurities is screened by the graphene electrons; Nomura and MacDonald used the Thomas-Fermi approximation (a local, semiclassical screening treatment), which strictly fails at the Dirac point where the clean density of states vanishes, so the value of σmin⁡\sigma^{\min} may depend on a more faithful, self-consistent treatment of electron-electron interactions. Third, in a spatially fluctuating long-range potential the local carrier density does not vanish on average: graphene near charge neutrality breaks up into electron-rich and hole-rich regions ('electron-hole puddles') through which current percolates, which is the now-standard competing picture of finite Dirac-point conductivity. The open issue is which description wins in the true zero-temperature, infinite-system limit, and what the quantitative value and universality of the minimum conductivity are.

Problem Statement

For the two-dimensional massless Dirac fermion description of graphene at the charge-neutrality (Dirac) point with disorder, determine the zero-temperature, infinite-system-size transport state and the quantitative value of the minimum conductivity σmin⁡\sigma^{\min}. Specifically: (1) establish the asymptotic scaling of the Dirac-point conductivity σ(E=0,L)\sigma(E=0, L) as L→∞L\to\infty for smooth (intravalley, screened Coulomb) disorder — does it diverge as expected for the symplectic symmetry class with weak antilocalization, flow to a finite universal critical value of order e2/he^2/h as predicted by the topological (θ=π\theta=\pi) nonlinear sigma model, or is the finite-size puddle-percolation picture the correct thermodynamic-limit description; (2) determine at what length scales and temperatures intervalley scattering becomes relevant for realistic disorder, whether it drives localization of the Dirac-point states in the true ground state, and what the conductivity is in that intervalley-dominated regime; and (3) fix, within a single self-consistent treatment of disorder and electron-electron interactions valid at the Dirac point (beyond Thomas-Fermi screening), the quantitative value and degree of universality of σmin⁡\sigma^{\min} as a function of the impurity density and disorder range, and reconcile it with the observed nonuniversal sample-dependent values. A complete answer must address all three parts; each part may be answered independently, and establishing any one of the three alternatives in (1) for the L→∞L\to\infty limit resolves part (1). The answer must address the thermodynamic-limit fate and not only accessible finite sizes.

References

  1. K. Nomura and A. H. MacDonald, "Quantum Transport of Massless Dirac Fermions", Phys. Rev. Lett. 98, 076602 (2007). DOI: 10.1103/PhysRevLett.98.076602; https://doi.org/10.1103/PhysRevLett.98.076602 (arXiv: cond-mat/0606589)
  2. P. M. Ostrovsky, I. V. Gornyi, and A. D. Mirlin, "Electron transport in disordered graphene", Phys. Rev. B 74, 235443 (2006). DOI: 10.1103/PhysRevB.74.235443; https://doi.org/10.1103/PhysRevB.74.235443 (arXiv: cond-mat/0609617)
  3. E. McCann, K. Kechedzhi, Vladimir I. Fal'ko, H. Suzuura, T. Ando, and B. L. Altshuler, "Weak-Localization Magnetoresistance and Valley Symmetry in Graphene", Phys. Rev. Lett. 97, 146805 (2006). DOI: 10.1103/PhysRevLett.97.146805; https://doi.org/10.1103/PhysRevLett.97.146805
  4. S. V. Morozov, K. S. Novoselov, M. I. Katsnelson, F. Schedin, L. A. Ponomarenko, D. Jiang, and A. K. Geim, "Strong Suppression of Weak Localization in Graphene", Phys. Rev. Lett. 97, 016801 (2006). DOI: 10.1103/PhysRevLett.97.016801; https://doi.org/10.1103/PhysRevLett.97.016801
  5. P. M. Ostrovsky, I. V. Gornyi, and A. D. Mirlin, "Quantum criticality and minimal conductivity in graphene with long-range disorder", Phys. Rev. Lett. 98, 256801 (2007). DOI: 10.1103/PhysRevLett.98.256801; https://doi.org/10.1103/PhysRevLett.98.256801 (arXiv: cond-mat/0702115)
  6. S. Adam, E. H. Hwang, V. M. Galitski, and S. Das Sarma, "A self-consistent theory for graphene transport", Proc. Natl. Acad. Sci. USA 104, 18392 (2007). DOI: 10.1073/pnas.0704772104; https://doi.org/10.1073/pnas.0704772104
  7. J. Martin, N. Akerman, G. Ulbricht, T. Lohmann, J. H. Smet, K. von Klitzing, and A. Yacoby, "Observation of electron-hole puddles in graphene using a scanning single-electron transistor", Nat. Phys. 4, 144 (2007). DOI: 10.1038/nphys781; https://doi.org/10.1038/nphys781
  8. E. Rossi and S. Das Sarma, "Ground State of Graphene in the Presence of Random Charged Impurities", Phys. Rev. Lett. 101, 166803 (2008). DOI: 10.1103/PhysRevLett.101.166803; https://doi.org/10.1103/PhysRevLett.101.166803
  9. I. S. Terekhov, A. I. Milstein, V. N. Kotov, and O. P. Sushkov, "Screening of Coulomb Impurities in Graphene", Phys. Rev. Lett. 100, 076803 (2008). DOI: 10.1103/PhysRevLett.100.076803; https://doi.org/10.1103/PhysRevLett.100.076803
  10. F. V. Tikhonenko, D. W. Horsell, R. V. Gorbachev, and A. K. Savchenko, "Weak Localization in Graphene Flakes", Phys. Rev. Lett. 100, 056802 (2008). DOI: 10.1103/PhysRevLett.100.056802; https://doi.org/10.1103/PhysRevLett.100.056802
  11. S. Das Sarma, Shaffique Adam, E. H. Hwang, and Enrico Rossi, "Electronic transport in two-dimensional graphene", Rev. Mod. Phys. 83, 407 (2011). DOI: 10.1103/RevModPhys.83.407; https://doi.org/10.1103/RevModPhys.83.407
  12. V. N. Kotov, B. Uchoa, V. M. Pereira, F. Guinea, and A. H. Castro Neto, "Electron-Electron Interactions in Graphene: Current Status and Perspectives", Rev. Mod. Phys. 84, 1067 (2012). DOI: 10.1103/RevModPhys.84.1067; https://doi.org/10.1103/RevModPhys.84.1067
  13. Weiwei Chen, Yedi Shen, Bo Fu, Qinwei Shi, and W. Zhu, "On the sample-dependent minimal conductivity in weakly disordered graphene", arXiv:2108.13787 (2021); https://arxiv.org/abs/2108.13787
  14. Yositake Takane, "Finite-Size Scaling Analysis of the Conductivity of Dirac Electrons on a Surface of Disordered Topological Insulators", J. Phys. Soc. Jpn. 85, 094715 (2016). DOI: 10.7566/JPSJ.85.094715; https://doi.org/10.7566/JPSJ.85.094715 (arXiv: 1608.07873)
  15. Swadeepan Nanda and Pavan Hosur, "Transport scaling and critical tilt effects in disordered two-dimensional Dirac fermions", Phys. Rev. B 113, 214203 (2026). DOI: 10.1103/6qrf-2ffm; https://doi.org/10.1103/6qrf-2ffm (arXiv: 2512.09133)