---
title: Variable Range Hopping in Disordered Materials
url: https://www.emergentmind.com/topics/variable-range-hopping-vrh
type: topic
---

# Variable Range Hopping in Disordered Materials

Variable Range Hopping (VRH) describes electronic transport in disordered materials where carriers are localized and conduction proceeds via phonon-assisted tunneling between spatially separated states. Unlike nearest-neighbor hopping, VRH involves a trade-off between tunneling distance and energy difference, resulting in temperature-dependent conductivity characterized by stretched exponential laws. The VRH framework provides a unifying description of transport in a wide range of amorphous, molecular, granular, and low-dimensional systems, and elucidates effects due to disorder, electron-electron interactions, and dimensionality.

## 1. Fundamental Principles and Universal Law

VRH is distinguished from activated hopping by the variable distance over which carriers hop. The rate for an electron (or hole) hopping from site $i$ to $j$, separated by $r_{ij}$ with energy mismatch $\epsilon_{ij}$, is
\[
\Gamma_{ij}\;\sim\;\exp\!\bigl[ -2 r_{ij}/\xi\;-\;|\epsilon_{ij}|/T \bigr]
\]
where $\xi$ is the localization length. The DC conductivity is dominated by the "optimal" hops that minimize the total exponent, subject to the condition that available final states exist within the hop range and energy window.

The resulting Mott law for conductivity in $d$ dimensions is
\[
\sigma(T) \;\simeq\; B(T) \exp\left[ - (T_0/T)^{1/(d+1)} \right]
\]
where $T_0 = \beta_d / [\rho \xi^d ]$, $\rho$ is the single-particle density of states (DOS) at the Fermi level, and $\beta_d$ is a numerical constant. The exponent $1/(d+1)$ reflects optimization over both spatial and energy variables [1404.3981], [1407.2478].

### Parameters and Scaling

- $\xi$: Localization length (nm scale in strong localization)
- $\rho$: Density of localized states at $E_F$
- $T_0$: Characteristic temperature, inversely proportional to $\rho$ and $\xi^d$
- $B(T)$: Weakly temperature-dependent prefactor, often sub-leading (power-law)

## 2. Extensions: Coulomb Gap and Efros–Shklovskii VRH

In the presence of long-range unscreened Coulomb interactions, a soft gap opens in the single-particle DOS at $E_F$ ("Coulomb gap"), leading to the Efros–Shklovskii (ES) law:
\[
\sigma(T) \propto \exp\left[ - (T_{ES}/T)^{1/2} \right]
\]
where $T_{ES} = C e^2 / (k_B \varepsilon \xi)$ and $C \sim 2.8$ [1702.06729], [1210.1876], [2106.08401]. This $1/2$ exponent is universal in any dimension under ES conditions. The VRH regime often exhibits a crossover from Mott ($p=1/(d+1)$) to ES ($p=1/2$) behavior as temperature decreases and interaction effects become dominant [1009.1795], [1702.06729].

## 3. Influence of Disorder, Dimensionality, and Energy-Dependent DOS

The VRH exponent $\alpha$ can differ from canonical values if the DOS near $E_F$ is energy-dependent or if the transport is constrained (e.g., quasi-1D structures). For a DOS $g(E)\propto |E-E_F|^p$, the hopping exponent generalizes to
\[
\alpha = \frac{p+1}{p+d+1}
\]
in the VRH expression $\sigma \sim \exp[ - (T_0/T)^{\alpha} ]$ [2203.07823]. As the dimensionality $d$ is reduced, or as $p$ increases (stronger pseudogap), the hopping law smoothly evolves towards nearest-neighbor hopping ($\alpha \to 1$). In disordered carbon aerogels, variation from $\alpha \approx 0.4$ to $0.9$ is observed as the material transitions from 3D to quasi-1D hopping [2203.07823].

In graphene with covalently attached impurities, the impurity wavefunction decays as a power-law rather than exponentially. There, VRH conductivity follows a power law of temperature, $\sigma \propto T^{\eta}$, with $\eta = 2s/d + s'$ and $s$ the localization exponent [1208.5026].

## 4. Advanced Theoretical Refinements: Memory Effects and Percolation

A rigorous treatment reveals that occupation-number "memory" effects strongly suppress cross-Fermi-energy hops: transitions that create an electron-hole pair rarely contribute to DC transport due to near-immediate recombination. This splits the VRH percolation network into two weakly coupled subnetworks ("two-color percolation"). The main transport proceeds via (a) electron-like hops between normally empty sites and (b) hole-like hops between normally full sites; type (c) cross-chemical-potential hops are exponentially suppressed [1404.3981]. The result is an additional subleading exponential correction to Mott's law,
\[
\sigma(T) = B(T)\,\exp\left[ - (T_0/T)^{1/(d+1)} + \alpha_d (T_0/T)^{\mu_d} \right]
\]
where $\mu_d$ is a universal exponent $(\mu_2 \approx 0.20,\, \mu_3 \approx 0.10)$. This refinement can explain unexpectedly large extracted prefactors in experimental data [1404.3981].

## 5. Thermoelectric Response and Seebeck Coefficient in the VRH Regime

Standard Mott VRH predicts a Seebeck coefficient (thermopower) scaling with temperature as
\[
S(T) \propto T^{(d-1)/(d+1)}
\]
from the Mott-Cutler formula. However, when localization length is energy-dependent, the refined scaling is
\[
S(T) \propto T^{d/(d+1)}
\]
as predicted by Kubo-Luttinger linear response with Anderson localization scaling [2204.02764]. Experimental data for thiospinel CuCrTiS$_4$ confirm this $T^{3/4}$ behavior for $d=3$ [2204.02764]. In CVD single-layer MoS$_2$, simultaneous measurement of $\sigma(T)$ and $S(T)$ verifies $S \propto T^{1/3}$ for $d=2$, with room-temperature $S$ up to $30$ mV/K [1407.2478].

## 6. Field-Driven and Nonlinear Regimes

In sufficiently large electric fields, the VRH process may become field-driven. The conductance crosses over from thermally-activated (Pollak–Riess) to field-driven (Shklovskii-type):
\[
G(T,E) \sim \exp\left[ - (T_0/T)^{1/(d+1)} + E/E_a \right] \;\; \text{(intermediate field)}
\]
\[
G(E) \sim \exp\left[ - (E_0/E)^{1/(d+1)} \right] \;\; \text{(high field)}
\]
The crossover field $E_c$ separating these regimes is analytically determined from the two characteristic field scales $E_0$ and $E_a$; e.g., in 2D, $E_c^{(2D)} = \left(E_a E_0^{1/3}/3\right)^{3/4}$ [1305.0315].

Non-Ohmic (field-driven) VRH is also directly observed in reduced graphene oxide sheets, with field scaling consistent with ES physics: $R(E) \propto \exp[(E_0/E)^{1/2}]$ [1210.1876].

## 7. Generalizations: Integral VRH and Transport Modeling in Disordered Systems

Conventional VRH theories rely on "optimal-hop" approximations, but integral VRH (IVRH) models provide a physically motivated formula for the full conductivity:
\[
\sigma(T) = A \int_{R_0}^\infty \exp\left[ -2\alpha R - \frac{\beta}{D_0 V(R) k_B T} \right] dR
\]
where $V(R)$ encodes system geometry (2D, 3D, multi-layer), and $\beta$ is a dimensionless factor obtained from Monte Carlo validation [2601.10226]. IVRH unifies the low-$T$ (Mott) and high-$T$ (Arrhenius) regimes without ad hoc regime separation, reproducing a smooth crossover and reducing fitting ambiguity in experimental data, with robust application demonstrated for monolayer MoS$_2$ and WS$_2$ [2601.10226].

## 8. Materials, Experimental Validation, and Defect Physics

VRH is ubiquitous in a wide range of systems, including:

- Polycrystalline ZnO, with a clear Mott$\to$ES crossover at cryogenic temperatures controlled by oxygen defects [1702.06729].
- Disordered organic semiconductors: energetically disordered systems require VRH models accounting for non-nearest-neighbor transitions; a critical ratio of inter-site distance to localization radius, $a_{NN}/a \approx 3$, determines the regime [2103.03856].
- Graphene and reduced graphene oxide: ES VRH dominates due to structural disorder and interaction-driven Coulomb gaps; measured localization lengths and extracted bandgaps validate theory [1210.1876].
- Pyrite (FeS$_2$) thin films: p-type conduction proceeds via defect-induced VRH associated with sulfur vacancy clusters [2106.08401].
- Semiconductor nanocrystal arrays: quantum-confinement-induced donor-number fluctuations generate the disordered Coulomb landscape essential for ES-VRH, with marked crossover to activated transport as disorder is reduced [1203.3889].

Negative magnetoresistance and weak-localization signatures co-exist with a variable $\alpha$ as the system transitions from VRH to nearest-neighbor hopping in quasi-one-dimensional disordered carbon aerogels [2203.07823].

## 9. Non-Standard Mechanisms and Environmental Coupling

Environmental effects, such as coupling to marginally localized phonons or proximity to metallic gates, modify VRH transport:

- In 1D, phonons with only a delocalized zero mode ("marginally localized phonon bath") suppress low-order phonon processes; VRH is dominated by exponentially rare many-phonon processes, resulting in a singular prefactor to the Mott law [1511.03676].
- Proximity to a metallic plate statically screens the Coulomb gap (restoring a flat DOS at low energy) and dynamically induces polaronic suppression of phonon-assisted hopping via coupling to an Ohmic bath of metal electrons. This introduces a multiplicative suppression factor in $\sigma(T)$, notable over a wide intermediate temperature regime, and can reduce VRH conductivity by more than an order of magnitude without altering the underlying exponent [2007.12137].

---

Variable Range Hopping remains a powerful conceptual and quantitative framework for understanding electronic conduction in disordered and low-dimensional systems. Recent theoretical advances, rigorous experimental validation, and development of integral models continue to extend its applicability and resolve longstanding questions pertaining to interaction effects, scaling of thermoelectric properties, and the roles of geometry, dimensionality, and disorder.

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