Charge and energy transport in graphene with smooth finite-range disorder

We investigate charge and energy transport in monolayer graphene with smooth finite-range disorder, modeled by soft impurity potentials. Using a continuum Dirac model, we go beyond the Born approximation by computing the exact scattering matrix for individual impurities. This captures the full nonperturbative physics of smooth disorder. From the exact scattering data, we evaluate transport coefficients by solving the Boltzmann equation with energy-resolved phase shifts. We analyze electrical and electronic thermal conductivities versus carrier density and temperature, including deviations from the Wiedemann-Franz law. Our results reveal that finite-range disorder nontrivially modifies charge and heat currents, especially at low energies where perturbative methods fail. These findings provide a more accurate transport characterization for disordered Dirac materials and clarify how smooth disorder governs energy flow in graphene.

💡 Research Summary

This paper presents a comprehensive theoretical study of charge and energy transport in monolayer graphene subject to smooth, finite‑range disorder, modeled as circular “soft‑sphere” impurity potentials of constant strength V₀ and radius R. Starting from the low‑energy Dirac Hamiltonian H₀ = χ v_F σ·p, the authors add a scalar potential V(r)=V₀ Θ(R−r) that is smooth on the lattice scale (R ≫ a) and therefore suppresses intervalley scattering. The impurity ensemble is assumed dilute (n_imp R² ≪ 1), allowing independent scattering events.

The core technical achievement is the exact solution of the single‑impurity scattering problem using a partial‑wave expansion. By separating the Dirac equation into inner (r < R) and outer (r > R) regions, the spinor components are expressed in terms of Bessel functions J_m(k r) and Y_m(k r). Continuity of the wavefunction at r = R yields analytic phase shifts δ_{m j} that depend on the product k_in R, where k_in = s k − V₀/(χħv_F) and k_out = s k. These phase shifts are independent of the valley index χ, reflecting the smoothness of the disorder. The exact transition matrix T(k′,k) is then constructed from the phase shifts, providing a non‑perturbative expression for the scattering amplitude that remains valid for strong potentials and low energies where the Born approximation fails.

Using the exact T‑matrix, the authors evaluate the elastic transition rate W(k′,k)=2πħ n_imp |T|² δ(E_k−E_k′) and, within the relaxation‑time approximation of the Boltzmann equation, define the transport scattering time τ_tr(E) with the usual (1 − cosθ) weighting. The two‑dimensional Dirac density of states ν₁(E)=|E|/(2πħ²v_F²) is employed to perform the energy integrals. Electrical conductivity σ, electronic thermal conductivity κ_e, Lorenz number L = κ_e/(σT), Seebeck coefficient S, and the electronic figure of merit ZT = S²σT/κ_e are obtained from standard semiclassical formulas, with the Fermi‑Dirac distribution providing temperature dependence.

Numerical evaluation reveals several key physical insights. First, the impurity radius R dominates the transport behavior: larger R reduces the effective scattering cross‑section, leading to longer τ_tr and consequently higher σ and κ_e. The potential strength V₀ influences the phase shifts only weakly unless R is very small; in that limit, even modest V₀ can cause strong back‑scattering and suppress conductivity. Near charge neutrality (μ≈0) the conductivity follows a sublinear σ ∝ |μ|^{1/2} dependence, markedly different from the linear behavior predicted by Born‑type theories, reflecting the non‑perturbative nature of the scattering.

Second, the Lorenz number deviates substantially from the universal Sommerfeld value L₀ = π²k_B²/3e². Depending on carrier density and temperature, L can be 2–5 times larger or smaller than L₀, indicating that momentum relaxation and energy relaxation are decoupled for smooth finite‑range scatterers. The Seebeck coefficient, obtained via the Mott‑type relation S ≈ (π²k_B²T/3e)(∂lnσ/∂E)_{E=μ}, inherits the non‑linear σ(E) and therefore shows pronounced variations, including sign changes, as a function of doping.

Third, the electronic figure of merit ZT reaches its maximum for impurity radii of a few nanometers and moderate V₀ (≈0.2 eV), yielding ZT≈0.2–0.4 when only the electronic contribution to thermal conductivity is considered. This represents a notable improvement over models based on point‑like δ‑function scatterers, suggesting that engineering the spatial extent of impurities can be an effective route to enhance thermoelectric performance in graphene.

The authors conclude that smooth finite‑range disorder provides a powerful knob to simultaneously tune charge and heat transport in Dirac materials. The radius R emerges as the crucial design parameter, offering a pathway for impurity‑engineering, substrate‑distance control, or nanostructuring to optimize both conductivity and thermoelectric efficiency. They propose future extensions that incorporate electron‑electron and electron‑phonon interactions, as well as direct comparisons with Kubo‑Greenwood calculations, to fully capture transport in realistic high‑mobility graphene devices.