Phonon-limited transport coefficients in extrinsic graphene

Abstract: The effect of electron-phonon scattering processes over the thermoelectric properties of extrinsic graphene was studied. Electrical and thermal resistivity, as well as the thermopower, were calculated within the Bloch theory approximations. Analytical expressions for the different transport coefficients were obtained from a variational solution of the Boltzmann equation. The phonon-limited electrical resistivity ρ{e-ph} shows a linear dependence at high temperatures, and follows ρ{e-ph} \sim T^{4} at low temperatures, in agreement with experiments and theory previously reported in the literature. The phonon-limited thermal resistivity at low temperatures exhibits a \sim T dependence, and achieves a nearly constant value at high temperatures. The predicted Seebeck coefficient at verylow temperatures is Q(T) \sim -π2 k_B T /(3 e E_F), which shows a n^{-1/2} dependence with the density of carriers, in agreement with experimental evidence. Our results suggest that thermoelectric properties can be controlled by adjusting the Bloch-Gruneisen temperature through its dependence on the extrinsic carrier density in graphene.

• The paper derives explicit analytic formulas for electrical resistivity, thermal resistivity, and thermopower in doped graphene using a variational Boltzmann approach.
• It demonstrates a transition in electrical resistivity from a T^4 dependence at low temperatures to a linear behavior at higher temperatures, which aligns well with experimental data.
• The study highlights how carrier density modulates transport coefficients, establishing design principles for graphene-based thermal and thermoelectric devices.

Phonon-Limited Transport Coefficients in Extrinsic Graphene

Introduction

This work presents a comprehensive theoretical analysis of phonon-limited transport coefficients—electrical and thermal resistivity, as well as thermopower—in extrinsic graphene under high carrier densities. The methodology is rooted in a variational solution to the Boltzmann equation within the Bloch theory framework, employing the deformation potential approximation for electron–longitudinal acoustic phonon interactions. By deriving explicit analytical expressions for each transport coefficient and benchmarking them with experimental data, the paper elucidates both the parameter dependencies and the underlying physical mechanisms that dictate charge and heat transport in graphene in the phonon-limited regime.

Methodology

The central model treats graphene as a two-dimensional Dirac electron system, with the Fermi level situated above the Dirac point due to carrier doping. The electron-phonon interaction is described by the deformation potential, justified by the near-spherical (circular) Fermi surface at high carrier concentration. The Boltzmann transport equation is approached variationally, a method that transcends the quasi-elasticity constraints of relaxation time approximations and converges to exact solutions within the variational ansatz.

Charge and heat currents are formally coupled via the Onsager reciprocal relations, leading to interrelations between the conductivities/resistivities and the Seebeck coefficient in terms of the Onsager matrix elements. The variational approach is applied using a minimal set of trial functions, allowing derivation of closed-form temperature- and density-dependent formulas for electrical and thermal resistivity and thermopower, all incorporating the relevant graphene band structure physics and electron-phonon scattering matrix elements.

Electrical Resistivity

Phonon-limited electrical resistivity $\rho_{e-ph}(T)$ exhibits a crossover from a low-temperature $T^4$ power law to a linear-in-$T$ dependence at elevated temperatures. The controlling energy scale is the Bloch-Grüneisen temperature $\Theta_{BG}$, which scales with carrier density as $\Theta_{BG} \propto n^{1/2}$. This behavior is captured through explicit analytic expressions involving dimensionless integrals $\mathcal{J}_p(z)$ over the phonon spectrum. The results align quantitatively with high-quality experimental data, accurately predicting both the magnitude and temperature dependence of $\rho_{e-ph}$ in the regime where electron–phonon scattering dominates (1201.1057).

The zero-temperature offset in resistivity, attributable to impurity scattering, is shown to depend on carrier concentration, displaying a sum of a constant and a $n^{-1}$ term, consistent with expectations from both short-range and long-range (Coulomb) disorder models.

Thermal Resistivity

The phonon-limited thermal resistivity $[\kappa_{e-ph}]^{-1}$ in extrinsic graphene displays a linear temperature dependence at low $T$. This is in marked contrast to the quadratic $T^2$ scaling characteristic of three-dimensional metals, and arises from the reduced phase space for electron-phonon scattering in two dimensions and the specifics of the Dirac spectrum. At high temperatures, $[\kappa_{e-ph}]^{-1}$ saturates to an approximately constant value.

When including impurity (Coulomb) scattering, the total thermal conductivity at room temperature is estimated to be $$\sim 400~\mathrm{W}\,\mathrm{m}^{-1}\mathrm{K}^{-1}$$, which is significantly lower than the phononic contribution ($$\sim 4300~\mathrm{W}\,\mathrm{m}^{-1}\mathrm{K}^{-1}$$). This hierarchy suggests dominant phonon-mediated thermal transport in graphene except at very high carrier densities or low temperatures.

Thermopower (Seebeck Coefficient)

The analytical formula for the Seebeck coefficient $Q(T)$ reveals a leading-order linear-in-$T$ dependence at low temperature:

$Q(T) \sim -\frac{\pi^2}{3e} \frac{k_B^2 T}{E_F}$

This yields a carrier density dependence $Q \propto n^{-1/2}$, consistent with the experimental findings, as $E_F \propto n^{1/2}$. Phonon-drag contribution to thermopower is subdominant at low temperatures compared to the diffusion part, which the theory shows agrees with observed trends.

The model predictions require only the deformation potential $D$ and the sound velocity $v_s$ as free parameters. The extracted values are in close agreement with independent theoretical and experimental characterizations, reinforcing the consistency of the approach.

Implications

The results demonstrate that electrical resistivity, thermal resistivity, and thermopower in extrinsic graphene can be manipulated by carrier density modulation, primarily via electrostatic gating. This tunability is mediated by the density dependence of the Bloch-Grüneisen temperature. The high fidelity of the theoretical framework in matching experimental data confirms its suitability for interpreting device physics in doped graphene, especially when engineering device thermal or thermoelectric properties is required.

The distinctive transport scalings (e.g., linear-$T$ thermal resistivity, $T^4$ to $T$ crossover in electrical resistivity) are direct consequences of graphene's linear band structure and two-dimensional phonon dynamics. The work also underscores the utility of the variational Boltzmann approach in systems where the relaxation time approximation is only marginally justified, particularly for inelastic scatterers like acoustic phonons.

Future Directions

The established theory forms a baseline for investigating transport properties in other two-dimensional Dirac materials, including those with modified phonon spectra or under various strain and substrate conditions. Extending these methods to multiband systems or incorporating strong electron-electron and phonon-phonon interactions presents avenues for deeper theoretical insights.

In practical terms, the demonstrated carrier density control over thermoelectric properties suggests design principles for graphene-based sensors, thermal switches, and low-dimensional thermoelectric converters. Further advances could stem from integrating this theory with ab initio calculations for explicit device geometries and disorder configurations.

Conclusion

This paper provides an exacting and validated theoretical treatment of phonon-limited transport coefficients in extrinsic graphene at high carrier densities, yielding analytic temperature- and density-dependent formulas. The results not only exhibit quantitative agreement with a range of experimental data but also clarify the mechanisms governing thermal and charge transport under phonon-limited conditions. The demonstrated density-tunability of these coefficients establishes a robust platform for both further experimental explorations and practical device engineering in graphene and related Dirac materials.