---
title: Integral Variable Range Hopping (IVRH)
url: https://www.emergentmind.com/topics/integral-variable-range-hopping-ivrh
type: topic
---

# Integral Variable Range Hopping (IVRH)

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 [2601.10226].

## 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 $β$, this gives:
\[
ΔE(R) = β / [D_0 V(R)]
\]
yielding the IVRH hopping kernel:
\[
ω(R) = ω_0 \exp\left[-2αR - \frac{β}{D_0 V(R) k_B T}\right]
\]
The form of $V(R)$ encodes geometric dependence:
- Homogeneous $d$-dimensional solid:
  \[
  V_d(R) = S_d \frac{R^d}{d}
  \]
  with $S_1=2$, $S_2=2π$, $S_3=4π$.

## 2. Integral Conductivity Model

The DC conductivity is constructed by integrating $ω(R)$ over all possible hopping distances, weighted by the number of candidate hops:
\[
σ(T) = A\int_{R_0}^{∞} \exp\left[-2αR - \frac{β}{D_0 V(R) k_B T}\right] dR
\]
with $A$ an overall prefactor, $R_0$ 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 [2601.10226].

Variable transformation and adaptive quadrature can be employed for numerical stability. All parameters ($α$, $D_0$, $R_0$, $A$, $β$) have direct and physically meaningful interpretations.

## 3. Limiting Behaviors: Arrhenius and Mott VRH Regimes

IVRH recovers known analytical limits:

- **High-Temperature (Arrhenius) Limit:** For $k_B T \gg β/(D_0 V(R_0))$,
  \[
  σ(T) \sim σ_∞ \exp\left[-\frac{E_a}{k_B T}\right]
  \]
  with $E_a = β/(D_0 V(R_0))$. Conductivity is dominated by nearest-neighbor hopping.

- **Low-Temperature (Mott VRH) Limit:** For $T \ll T_c$, the integral is estimated via steepest descent, giving Mott’s law:
  \[
  σ(T) \sim σ_0 \exp\left[-(T_0/T)^{1/(d+1)}\right]
  \]
  where $T_0 = β d/(2α D_0 S_d k_B)$, and the exponent $p = 1/(d+1)$, e.g., $p = 1/3$ in 2D, $p = 1/4$ in 3D.

The model automatically captures the crossover between these regimes without artificial partitioning of the temperature range or empirical interpolation. The crossover temperature $T_c$ 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 $H$-layered system,
\[
\overline V(R) = \frac{1}{H} \sum_{i=1}^{H} V_i(R)
\]
with each $V_i(R)$ describing a truncated-sphere volume within each layer. This yields IVRH predictions that interpolate continuously between pure 2D ($H=1$) and bulk 3D.

Universal geometric scaling factors emerge from the dependence of $V(R)$ on the dimensionality and stacking sequence. The applicability extends to nanoribbons and networks by explicitly calculating the geometric $V(R)$ 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: $20\times20\times1$ (2D), $20\times20\times20$ (3D), periodic boundaries;
- Hopping: Random spatial proposals weighted by $\exp(-2αR_{ij})$;
- Accept/reject step governed by energy change $\Delta E$ 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 $β$:
\[
β_{2D} = 1.06 \pm 0.02,\quad β_{3D} = 6.68 \pm 0.21
\]
The spread in inferred $β$ is an order of magnitude lower than that of the exponent $p$ in classical VRH, permitting robust dimensional assignment and reducing model ambiguity [2601.10226].

## 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 $σ(T)$ over the full thermal span, with parameters $\{A, α, D_0, R_0\}$ and fixed $β=1.06$ (2D). Extracted localization length $\ell = α^{-1} \approx 2$–$3$ nm is invariant under carrier density changes, and the Arrhenius–Mott crossover temperature $T_c$ is sharply defined.
- **WS₂:** Accurately captures the gate-induced insulator–metal transition by extracting $α^{-1}(V_g)$, revealing a universal scaling of $\ell$ 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 ($α^{-1}$), density of states ($D_0$), and hopping cutoff ($R_0$). The only empirical parameter, $β(d)$, 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 $V(R)$, 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 [2601.10226].

Source: https://www.emergentmind.com/topics/integral-variable-range-hopping-ivrh