Integral Variable Range Hopping (IVRH)

• IVRH is a theoretical framework for electronic transport that models temperature dependence with an integral over hopping probabilities, capturing both low-temperature Mott and high-temperature Arrhenius behaviors.
• It constructs conductivity by integrating over all possible hopping distances weighted by energy and geometric factors, thereby yielding physically meaningful parameters without empirical adjustments.
• Monte Carlo simulations validate IVRH, and its application to monolayer materials like MoS₂ and WS₂ shows robust, gate-tunable transport properties across differing dimensional regimes.

Integral Variable Range Hopping (IVRH) is a theoretical framework for modeling temperature-dependent electrical transport in disordered systems. IVRH refines and extends the Variable Range Hopping (VRH) paradigm by constructing the temperature-dependent conductivity from a physics-driven integral over hopping probabilities, rather than the empirical, single-exponent temperature laws of standard VRH. IVRH reproduces both Mott VRH behavior at low temperatures and Arrhenius (nearest-neighbor hopping) at high temperatures, yielding a smooth, parameter-free crossover that captures the full thermal range of hopping conduction. This model accurately describes low-dimensional, amorphous materials, including monolayer transition metal dichalcogenides, and produces physically meaningful fit parameters with increased robustness and interpretive power (Qin et al., 15 Jan 2026).

1. Formulation of the Hopping Probability and Effective Volume

The core of IVRH is the hopping probability $ω(R,ΔE)$ between localized states separated by distance $R$ and energy difference $ΔE$: $ω(R,ΔE) \propto \exp[-2αR - ΔE/(k_B T)]$ where:

• $α$ = inverse localization length,
• $k_B$ = Boltzmann constant,
• $T$ = temperature.

$ΔE$ is estimated by demanding the presence of one accessible final state within an effective “hopping volume” $V(R)$. With $D_0$ as the constant density of localized states at the Fermi level and Monte Carlo-determined numerical factor $R$0, this gives: $R$1

yielding the IVRH hopping kernel: $R$2 The form of $R$3 encodes geometric dependence:

• Homogeneous $R$4-dimensional solid:

$R$5

with $R$6, $R$7, $R$8.

2. Integral Conductivity Model

The DC conductivity is constructed by integrating $R$9 over all possible hopping distances, weighted by the number of candidate hops: $ΔE$0 with $ΔE$1 an overall prefactor, $ΔE$2 a nearest-neighbor cutoff. This integral framework departs from standard VRH by not assuming a single optimal hop length, but rather summing the contributions of all relevant hops (Qin et al., 15 Jan 2026).

Variable transformation and adaptive quadrature can be employed for numerical stability. All parameters ($ΔE$3, $ΔE$4, $ΔE$5, $ΔE$6, $ΔE$7) have direct and physically meaningful interpretations.

3. Limiting Behaviors: Arrhenius and Mott VRH Regimes

IVRH recovers known analytical limits:

• High-Temperature (Arrhenius) Limit: For $ΔE$8,

$ΔE$9

with $ω(R,ΔE) \propto \exp[-2αR - ΔE/(k_B T)]$0. Conductivity is dominated by nearest-neighbor hopping.

• Low-Temperature (Mott VRH) Limit: For $ω(R,ΔE) \propto \exp[-2αR - ΔE/(k_B T)]$1, the integral is estimated via steepest descent, giving Mott’s law:

$ω(R,ΔE) \propto \exp[-2αR - ΔE/(k_B T)]$2

where $ω(R,ΔE) \propto \exp[-2αR - ΔE/(k_B T)]$3, and the exponent $ω(R,ΔE) \propto \exp[-2αR - ΔE/(k_B T)]$4, e.g., $ω(R,ΔE) \propto \exp[-2αR - ΔE/(k_B T)]$5 in 2D, $ω(R,ΔE) \propto \exp[-2αR - ΔE/(k_B T)]$6 in 3D.

The model automatically captures the crossover between these regimes without artificial partitioning of the temperature range or empirical interpolation. The crossover temperature $ω(R,ΔE) \propto \exp[-2αR - ΔE/(k_B T)]$7 is implicitly defined by equating the two dominant contributions.

4. Dimensionality, Layered Systems, and Universal Scaling

IVRH is generalizable to multilayered and finite-thickness geometries. For an $ω(R,ΔE) \propto \exp[-2αR - ΔE/(k_B T)]$8-layered system,

$ω(R,ΔE) \propto \exp[-2αR - ΔE/(k_B T)]$9

with each $α$0 describing a truncated-sphere volume within each layer. This yields IVRH predictions that interpolate continuously between pure 2D ($α$1) and bulk 3D.

Universal geometric scaling factors emerge from the dependence of $α$2 on the dimensionality and stacking sequence. The applicability extends to nanoribbons and networks by explicitly calculating the geometric $α$3 corresponding to system topology.

5. Monte Carlo Validation and Parameter Inference

IVRH parameters are validated via Monte Carlo simulations based on tight-binding lattice Hamiltonians with uniform disorder. Key simulation features include:

• Lattices: $α$4 (2D), $α$5 (3D), periodic boundaries;
• Hopping: Random spatial proposals weighted by $α$6;
• Accept/reject step governed by energy change $α$7 under a small applied field and the Fermi-Dirac distribution;
• Conductivity extracted from mean displacement in linear response.

Monte Carlo fitting delivers highly stable dimension-specific $α$8: $α$9 The spread in inferred $k_B$0 is an order of magnitude lower than that of the exponent $k_B$1 in classical VRH, permitting robust dimensional assignment and reducing model ambiguity (Qin et al., 15 Jan 2026).

6. Experimental Applications: Monolayer MoS₂ and WS₂

IVRH provides a unified description of transport in monolayer MoS₂ and WS₂, resolving prior ambiguities in fitting regimes:

• MoS₂: Fits the measured $k_B$2 over the full thermal span, with parameters $k_B$3 and fixed $k_B$4 (2D). Extracted localization length $k_B$5–$k_B$6 nm is invariant under carrier density changes, and the Arrhenius–Mott crossover temperature $k_B$7 is sharply defined.
• WS₂: Accurately captures the gate-induced insulator–metal transition by extracting $k_B$8, revealing a universal scaling of $k_B$9 with gate voltage. Fitting errors are reduced compared to traditional two-regime approaches.

Observations support IVRH’s applicability to gate-tunable and chemically modulated two-dimensional systems.

7. Physical Interpretation and Outlook

IVRH supplants empirical fitting of VRH and Arrhenius laws with a physics-based integral whose parameters correspond directly to microscopic mechanisms: localization length ($T$0), density of states ($T$1), and hopping cutoff ($T$2). The only empirical parameter, $T$3, is precisely fixed for each dimension by simulation. No arbitrary partitioning of experimental data or manual adjustment of the Mott exponent is required.

Potential extensions include modeling with energy-dependent density of states (to incorporate Coulomb gap effects), application to exotic network topologies via calculation of non-trivial $T$4, and integration with first-principles disorder models for amorphous two-dimensional materials.

IVRH provides a quantitatively robust, physically interpretable platform for analyzing electrical transport in disordered systems across all temperature regimes and dimensionalities (Qin et al., 15 Jan 2026).