Hopping-Mediated Charge Transport in Graphene Beyond the Ballistic Regime

Douglas S. Galvao; J. P. Dadario Pereira; Luiz A. Ribeiro Junior; Raphael Tromer

arxiv: 2604.16152 · v1 · submitted 2026-04-17 · ❄️ cond-mat.mes-hall · cond-mat.mtrl-sci

Hopping-Mediated Charge Transport in Graphene Beyond the Ballistic Regime

J. P. Dadario Pereira , Raphael Tromer , Luiz A. Ribeiro Junior , Douglas S. Galvao This is my paper

Pith reviewed 2026-05-10 07:32 UTC · model grok-4.3

classification ❄️ cond-mat.mes-hall cond-mat.mtrl-sci

keywords graphenecharge transportkinetic Monte Carlohopping transportballistic regimevacanciesmobilitydisorder

0 comments

The pith

A kinetic Monte Carlo framework describes charge transport in graphene by hopping on atomic sites, revealing nearly ohmic behavior in pristine samples.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper develops a computational method to simulate how electrons move through graphene when not in the perfect ballistic regime. It uses random hops between carbon atoms to account for temperature, defects like vacancies, strain, and magnetic fields. In clean graphene, the simulations produce currents around 7 to 8 microamperes at 0.1 volts with transmittance close to 1, indicating almost ohmic conduction. The same approach calculates diffusion constants and mobilities from how far carriers travel in given times. This matters for understanding real graphene devices where perfect coherence does not hold and imperfections are present.

Core claim

The authors establish that charge transport beyond the ideal ballistic and coherent limits in graphene can be modeled through kinetic Monte Carlo hopping on a fixed atomic lattice. This method extracts current and transmittance directly from simulated carrier paths and incorporates effects from bias, temperature, magnetic fields, strain, and vacancies. Application to pristine graphene yields an almost ohmic response characterized by currents of approximately 7-8 μA, transmittance values near 0.98-1.00, and conductance in the range (5.8-7.8) × 10^{-5} S at 0.10 V bias, varying slightly with direction.

What carries the argument

kinetic Monte Carlo hopping on a predefined atomic lattice, which generates stochastic carrier trajectories to compute transport properties without phenomenological coefficients

Load-bearing premise

That persistent quantum coherence and interference effects remain unimportant, allowing classical kinetic Monte Carlo hopping to capture the essential transport physics.

What would settle it

Direct experimental measurement showing that conductance in pristine graphene at 0.1 V deviates substantially from the predicted 5.8-7.8 × 10^{-5} S, or exhibits strong non-ohmic characteristics at low temperatures.

Figures

Figures reproduced from arXiv: 2604.16152 by Douglas S. Galvao, J. P. Dadario Pereira, Luiz A. Ribeiro Junior, Raphael Tromer.

**Figure 2.** Figure 2: Bias-dependent transport response of graphene under randomly distributed va [PITH_FULL_IMAGE:figures/full_fig_p016_2.png] view at source ↗

**Figure 3.** Figure 3: Temperature-dependent transport response of graphene with randomly distributed [PITH_FULL_IMAGE:figures/full_fig_p018_3.png] view at source ↗

**Figure 4.** Figure 4: Magnetic-field-dependent current–voltage characteristics of graphene with ran [PITH_FULL_IMAGE:figures/full_fig_p021_4.png] view at source ↗

**Figure 5.** Figure 5: Current–voltage characteristics of graphene under combined structural disorder [PITH_FULL_IMAGE:figures/full_fig_p025_5.png] view at source ↗

**Figure 6.** Figure 6: Bias dependence of the effective diffusion coefficient [PITH_FULL_IMAGE:figures/full_fig_p029_6.png] view at source ↗

**Figure 7.** Figure 7: Bias-dependent effective transmitance T(V ) of graphene under a perpendicular magnetic field for different vacancy concentrations. Results are shown for vacancy concentrations of 0%, 5%, and 10%, with transport resolved along the X and Y directions, for magnetic field strengths B = 0, 2, 5, 7, and 10 T. Vacancies are randomly distributed within the lattice, while temperature, electrostatic bias profile, a… view at source ↗

**Figure 8.** Figure 8: Bias-dependent effective transmitance T(V ) of strained graphene with randomly distributed vacancies. Results are shown for vacancy concentrations of 0%, 5%, and 10% (rows), strain magnitudes of 2%, 7%, and 10% (columns), and three strain configurations: uniaxial strain along X, uniaxial strain along Y , and biaxial strain (XY ). Transport is resolved along the X and Y directions, as indicated by open and … view at source ↗

**Figure 9.** Figure 9: Temperature-dependent transport response of phagraphene. The top panels show [PITH_FULL_IMAGE:figures/full_fig_p035_9.png] view at source ↗

**Figure 10.** Figure 10: Magnetic-field-dependent current–voltage characteristics of phagraphene. The [PITH_FULL_IMAGE:figures/full_fig_p036_10.png] view at source ↗

**Figure 11.** Figure 11: Bias-dependent effective transmittance [PITH_FULL_IMAGE:figures/full_fig_p037_11.png] view at source ↗

**Figure 12.** Figure 12: Strain-dependent transport response of phagraphene. The top panels show the [PITH_FULL_IMAGE:figures/full_fig_p038_12.png] view at source ↗

read the original abstract

We present a trajectory-resolved framework for charge transport in graphene and related two-dimensional carbon systems beyond the ideal ballistic and fully coherent limits. Transport is described by kinetic Monte Carlo hopping on a predefined atomic lattice, allowing the combined treatment of disorder, thermal activation, and external fields. Current and effective transmittance are extracted directly from stochastic carrier trajectories, without phenomenological transport coefficients. We apply the method to graphene under bias voltage (0-0.10 V), temperature (300-900 K), magnetic field (0-10 T), in-plane strain (2-10%, uniaxial and biaxial), and vacancy concentration (0-10%). Pristine graphene shows an almost ohmic response, with currents of about 7-8 uA, effective transmittance near 0.98-1.00, and conductance of about (5.8-7.8) x 10^-5 S at 0.10 V, depending on direction. Vacancies strongly suppress transport, reducing transmittance to about 0.45-0.75 at 10% vacancy. Higher temperature accelerates hopping and partly restores transport, but cannot overcome severe connectivity loss. Magnetic fields further reduce transport, especially in disordered networks. The framework provides a unified computational scheme for realistic two-dimensional carbon materials and also yields diffusion coefficients and effective mobilities from carrier displacements and transit times.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

This KMC hopping framework gives a workable tool for disordered graphene under multiple fields but cannot credibly reproduce near-ballistic transmittance in the clean limit.

read the letter

The paper's main offering is a trajectory-resolved kinetic Monte Carlo scheme that lets carriers hop on a fixed atomic lattice while incorporating vacancies, temperature, bias, magnetic fields, and strain in one run. Current and transmittance come straight from the paths, along with diffusion coefficients and mobilities derived from displacements and transit times. The authors run sweeps across 0-0.10 V, 300-900 K, 0-10 T, 2-10% strain, and 0-10% vacancies, reporting nearly ohmic behavior in pristine graphene with transmittance 0.98-1.00 and conductance around 6-8 x 10^-5 S at 0.1 V, plus strong suppression by vacancies that temperature only partly offsets. That unified setup is the practical advance for anyone modeling realistic 2D carbon devices where defects and external influences coexist. It avoids phenomenological coefficients and works directly from stochastic trajectories, which is cleaner than some earlier approaches. The clean-limit numbers are the clear weak point. Classical hopping assumes loss of phase coherence, yet the reported transmittance near unity in pristine sheets matches what one expects only when Klein tunneling and interference survive. The abstract gives no derivation of the hopping rates, no explicit check against ballistic NEGF or Landauer results, and no error bars, so it is impossible to tell whether the crossover is physical or just an artifact of rate choice. The stress-test concern lands: without coherence built in, the model cannot self-consistently justify the pristine results it claims. Readers who need quick parameter studies of disordered graphene will find the code-like structure useful even if they later add quantum corrections. The work is coherent on its own terms and shows honest engagement with the parameter space, so it merits referee time. I would send it to review with instructions to add rate validation and ballistic benchmarks.

Referee Report

2 major / 3 minor

Summary. The manuscript introduces a trajectory-resolved kinetic Monte Carlo (KMC) framework for charge transport in graphene and related 2D carbon materials that operates beyond ideal ballistic and fully coherent limits. Transport is modeled via incoherent hopping on a predefined atomic lattice, with current and effective transmittance extracted directly from stochastic carrier trajectories under bias (0-0.10 V), temperature (300-900 K), magnetic field (0-10 T), strain (2-10%), and vacancy disorder (0-10%). For pristine graphene the method yields nearly ohmic response, currents of 7-8 μA, transmittance 0.98-1.00, and conductance (5.8-7.8)×10^{-5} S at 0.10 V; vacancies suppress transmittance to 0.45-0.75 at 10% concentration. The framework also reports diffusion coefficients and effective mobilities derived from carrier displacements and transit times.

Significance. If validated, the approach supplies a single computational scheme capable of treating realistic disorder, thermal activation, and external fields in 2D carbon systems while directly yielding transport observables without fitted phenomenological coefficients. The reported near-unity transmittance in the pristine case and the extraction of mobility/diffusion quantities from the same trajectories would be useful for device-scale modeling if the underlying model can be shown to be consistent across regimes.

major comments (2)

[Abstract] Abstract: the reported effective transmittance of 0.98-1.00 (and conductance (5.8-7.8)×10^{-5} S) for pristine graphene at 0.10 V cannot be reconciled with a purely classical, incoherent KMC hopping model on a fixed lattice. Ballistic graphene transport is dominated by phase-coherent effects (Klein tunneling, Fabry-Pérot resonances, weak localization) that are absent from stochastic hopping trajectories; the manuscript therefore provides no internal justification for recovering the ballistic limit from the same hopping-rate and connectivity rules used for the disordered cases.
[Abstract] Abstract and method description: the hopping rates themselves are not derived or validated against known ballistic or Landauer-Büttiker results. Without an explicit demonstration that the chosen rates and lattice connectivity reproduce the expected ballistic conductance (or at least the correct scaling with system size) when disorder is removed, the claim of a “unified computational scheme” for the ballistic-to-hopping crossover remains unsupported.

minor comments (3)

The abstract quotes numerical values (currents, transmittance, conductance) without accompanying error bars, number of trajectories, or convergence checks typical for stochastic KMC simulations.
No comparison is mentioned to established experimental or theoretical benchmarks for graphene (e.g., minimum conductivity, strain-induced gap opening, or magnetotransport data), which would strengthen the credibility of the reported trends.
Notation for “effective transmittance” and “conductance” extracted from trajectories should be defined explicitly; it is unclear whether these quantities are normalized to the Landauer conductance quantum or to the geometric cross-section.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful and constructive review. The comments correctly identify the need for greater clarity on how a classical hopping model relates to the ballistic regime. We address each point below and describe the revisions that will be incorporated.

read point-by-point responses

Referee: [Abstract] Abstract: the reported effective transmittance of 0.98-1.00 (and conductance (5.8-7.8)×10^{-5} S) for pristine graphene at 0.10 V cannot be reconciled with a purely classical, incoherent KMC hopping model on a fixed lattice. Ballistic graphene transport is dominated by phase-coherent effects (Klein tunneling, Fabry-Pérot resonances, weak localization) that are absent from stochastic hopping trajectories; the manuscript therefore provides no internal justification for recovering the ballistic limit from the same hopping-rate and connectivity rules used for the disordered cases.

Authors: We agree that the model is strictly classical and incoherent and therefore cannot reproduce quantum interference phenomena such as Klein tunneling or weak localization. The near-unity transmittance obtained for the pristine lattice follows directly from the absence of scattering sites: every hop is allowed and the stochastic trajectories sample the fully connected network, yielding high net transmission. This classical high-transmission limit is used only as a reference point against which the suppression caused by vacancies, strain, and magnetic field is quantified. In the revised manuscript we will modify the abstract and add a short clarifying paragraph in the introduction that explicitly distinguishes the classical hopping baseline from the quantum ballistic regime while retaining the claim that the same rate rules are applied uniformly across all disorder levels. revision: partial
Referee: [Abstract] Abstract and method description: the hopping rates themselves are not derived or validated against known ballistic or Landauer-Büttiker results. Without an explicit demonstration that the chosen rates and lattice connectivity reproduce the expected ballistic conductance (or at least the correct scaling with system size) when disorder is removed, the claim of a “unified computational scheme” for the ballistic-to-hopping crossover remains unsupported.

Authors: The hopping rates are constructed from the local electrostatic energy difference produced by the applied bias together with a fixed attempt frequency and a temperature-activated factor; no phenomenological fitting parameters are introduced. In the pristine case this prescription produces length-independent transmission for the simulated device sizes. To meet the referee’s request we will insert in the Methods section a concise derivation of the rate expression and an additional figure (or supplementary panel) that plots effective transmittance versus sample length for the pristine lattice, confirming the expected saturation to near-unity values. This addition will be accompanied by a brief statement that the model is intended as an effective classical scheme rather than a direct replacement for Landauer-Büttiker calculations. revision: yes

Circularity Check

0 steps flagged

No circularity: direct extraction from trajectories

full rationale

The paper's central method extracts current and transmittance directly from stochastic KMC trajectories on a fixed lattice, without any fitted parameters, self-referential definitions, or load-bearing self-citations. Pristine-graphene results (near-unity transmittance) emerge as simulation outputs under the model's assumptions rather than being presupposed or renamed from inputs. No step reduces by construction to its own premises.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

Based solely on the abstract, the central claim rests on the domain assumption that hopping on a fixed lattice captures the essential physics; no explicit free parameters, new entities, or additional axioms are stated.

axioms (1)

domain assumption Charge transport beyond the ballistic regime can be modeled via stochastic hopping on a predefined atomic lattice
This is the foundational premise of the kinetic Monte Carlo trajectory method described in the abstract.

pith-pipeline@v0.9.0 · 5561 in / 1424 out tokens · 81315 ms · 2026-05-10T07:32:59.146185+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

6 extracted references

[1]

G.; Barros, E

(1) Meunier, V.; Souza Filho, A. G.; Barros, E. B.; Dresselhaus, M. S. Physical properties of low-dimensionalsp 2-based carbon nanostructures.Rev. Mod. Phys.2016,88, 025005. (2) Castro Neto, A. H.; Guinea, F.; Peres, N. M.; Novoselov, K. S.; Geim, A. K. The electronic properties of graphene.Reviews of modern physics2009,81, 109–162. (3) Das Sarma, S.; Ada...

2016
[2]

Impurity Conduction at Low Concentrations.Phys

(10) Miller, A.; Abrahams, E. Impurity Conduction at Low Concentrations.Phys. Rev. 1960,120, 745–755. (11) Shklovskii, B. I.; Efros, A. L.Electronic properties of doped semiconductors; Springer, 1984; pp 202–227. (12) de Sousa Araujo Cassiano, T.; Ribeiro Junior, L. A.; Magela e Silva, G.; Henrique de Oliveira Neto, P. Strain-Tuneable Bipolaron Stability ...

1960
[3]

(34) Kramer, B.; MacKinnon, A

(33) Grundmann, M.The Physics of Semiconductors: An Introduction Including Nanophysics and Applications; Springer International Publishing: Cham, 2021; pp 135–178. (34) Kramer, B.; MacKinnon, A. Localization: theory and experiment.Reports on Progress in Physics1993,56,

2021
[4]

Charge transport in disordered organic photoconductors

(35) Bässler, H. Charge transport in disordered organic photoconductors. A Monte Carlo simulation study.Physica Status Solidi B (Basic Research);(Germany)1993,175. 42 (36) Lherbier, A.; Dubois, S. M.-M.; Declerck, X.; Niquet, Y.-M.; Roche, S.; Charlier, J.- C. Transport properties of graphene containing structural defects.Physical Review B—Condensed Matte...

1993
[5]

Localiza- tion landscape theory of disorder in semiconductors

(41) Filoche, M.; Piccardo, M.; Wu, Y.-R.; Li, C.-K.; Weisbuch, C.; Mayboroda, S. Localiza- tion landscape theory of disorder in semiconductors. I. Theory and modeling.Physical Review B2017,95, 144204. (42) Piccardo, M.; Li, C.-K.; Wu, Y.-R.; Speck, J. S.; Bonef, B.; Farrell, R. M.; Filoche, M.; Martinelli, L.; Peretti, J.; Weisbuch, C. Localization lands...

2017
[6]

B.; Yoo, J

44 (53) Skákalová, V.; Kaiser, A. B.; Yoo, J. S.; Obergfell, D.; Roth, S. Correlation between resistance fluctuations and temperature dependence of conductivity in graphene.Phys. Rev. B2009,80, 153404. (54) Raikh, M.; Wessels, G. Single-scattering-path approach to the negative magnetore- sistance in the variable-range-hopping regime for two-dimensional el...

2006

[1] [1]

G.; Barros, E

(1) Meunier, V.; Souza Filho, A. G.; Barros, E. B.; Dresselhaus, M. S. Physical properties of low-dimensionalsp 2-based carbon nanostructures.Rev. Mod. Phys.2016,88, 025005. (2) Castro Neto, A. H.; Guinea, F.; Peres, N. M.; Novoselov, K. S.; Geim, A. K. The electronic properties of graphene.Reviews of modern physics2009,81, 109–162. (3) Das Sarma, S.; Ada...

2016

[2] [2]

Impurity Conduction at Low Concentrations.Phys

(10) Miller, A.; Abrahams, E. Impurity Conduction at Low Concentrations.Phys. Rev. 1960,120, 745–755. (11) Shklovskii, B. I.; Efros, A. L.Electronic properties of doped semiconductors; Springer, 1984; pp 202–227. (12) de Sousa Araujo Cassiano, T.; Ribeiro Junior, L. A.; Magela e Silva, G.; Henrique de Oliveira Neto, P. Strain-Tuneable Bipolaron Stability ...

1960

[3] [3]

(34) Kramer, B.; MacKinnon, A

(33) Grundmann, M.The Physics of Semiconductors: An Introduction Including Nanophysics and Applications; Springer International Publishing: Cham, 2021; pp 135–178. (34) Kramer, B.; MacKinnon, A. Localization: theory and experiment.Reports on Progress in Physics1993,56,

2021

[4] [4]

Charge transport in disordered organic photoconductors

(35) Bässler, H. Charge transport in disordered organic photoconductors. A Monte Carlo simulation study.Physica Status Solidi B (Basic Research);(Germany)1993,175. 42 (36) Lherbier, A.; Dubois, S. M.-M.; Declerck, X.; Niquet, Y.-M.; Roche, S.; Charlier, J.- C. Transport properties of graphene containing structural defects.Physical Review B—Condensed Matte...

1993

[5] [5]

Localiza- tion landscape theory of disorder in semiconductors

(41) Filoche, M.; Piccardo, M.; Wu, Y.-R.; Li, C.-K.; Weisbuch, C.; Mayboroda, S. Localiza- tion landscape theory of disorder in semiconductors. I. Theory and modeling.Physical Review B2017,95, 144204. (42) Piccardo, M.; Li, C.-K.; Wu, Y.-R.; Speck, J. S.; Bonef, B.; Farrell, R. M.; Filoche, M.; Martinelli, L.; Peretti, J.; Weisbuch, C. Localization lands...

2017

[6] [6]

B.; Yoo, J

44 (53) Skákalová, V.; Kaiser, A. B.; Yoo, J. S.; Obergfell, D.; Roth, S. Correlation between resistance fluctuations and temperature dependence of conductivity in graphene.Phys. Rev. B2009,80, 153404. (54) Raikh, M.; Wessels, G. Single-scattering-path approach to the negative magnetore- sistance in the variable-range-hopping regime for two-dimensional el...

2006