Vortex VRH in Disordered 2D Superconductors

• Vortex variable-range hopping is a framework that describes nonzero resistance in disordered 2D superconducting films due to thermally activated and quantum-tunneling vortex motion over random pinning sites.
• The 2D Mott VVRH law, characterized by a 1/3 exponent, quantitatively links vortex hopping to temperature and magnetic field through a well-defined field-dependent parameter T₀.
• Experimental studies on amorphous InOₓ films confirm the VVRH predictions and reveal duality symmetry at the superconductor–insulator transition, highlighting unresolved theoretical puzzles.

Vortex variable–range hopping (VVRH) is a theoretical and experimental framework describing low-temperature, finite-resistance transport in disordered two-dimensional (2D) superconducting films subjected to a perpendicular magnetic field. In such systems, despite the expectation of zero resistance below the superconducting transition, a nonzero resistance persists due to the thermally activated and quantum-tunneling motion of magnetic vortices among a disordered landscape of pinning sites. VVRH draws a close analogy between vortex dynamics in a superconducting condensate and charge transport via variable-range hopping in Anderson-localized insulators, providing a predictive tool for analyzing the resistive behavior of strongly disordered superconducting films (Percher et al., 2017).

1. Theoretical Foundations and Physical Significance

In ideal, clean superconductors below the upper critical field $H_{c2}$, flux vortices are immobilized by perfect crystalline order, resulting in dissipationless transport. In contrast, disordered 2D superconductors contain randomly distributed weak pinning sites capable of trapping vortices. At finite temperature, thermally activated vortex motion (flux creep) leads to resistance described by an Arrhenius law,

$R(T) \propto \exp\left(\frac{E_p}{k_B T}\right)\,,$

where $E_p$ is a characteristic pinning barrier. At temperatures below $E_p/k_B$, quantum tunneling of vortices through these barriers becomes significant, especially due to the broad distribution of barrier heights and hopping distances. This leads to a regime where the resistance follows a variable-range hopping (VRH) temperature dependence analogous to that found in disordered electronic systems. The physical significance of VVRH is twofold: it predicts the persistence of nonzero resistance down to the lowest measurable temperatures on the superconducting side of the field-tuned superconductor–insulator transition (SIT), and it underlines the analogy between vortex motion and VRH of localized charges in electronic insulators (Percher et al., 2017).

2. Mathematical Formulation: 2D Mott VVRH Law

The mapping of vortex motion to charge hopping underlines the derivation of the VVRH law. For spatial dimension $d$, the Mott law gives the vortex conductance as

$$\sigma_v(T)\;\propto\;\exp\left[-\left(\frac{T_0}{T}\right)^{1/(d+1)}\right]\,,$$

where $T_0$ depends on the density of localized states and the localization length. For $d=2$, this yields the 2D Mott VVRH law,

$R(T) = R_0\,\exp\left[-(T_0/T)^{1/3}\right]\,,$

with

$T_0 = \frac{\beta}{k_B\,g(\mu)\,a^2}\,,$

where $R(T) \propto \exp\left(\frac{E_p}{k_B T}\right)\,,$0, $R(T) \propto \exp\left(\frac{E_p}{k_B T}\right)\,,$1 is the 2D density of states at the vortex chemical potential, $R(T) \propto \exp\left(\frac{E_p}{k_B T}\right)\,,$2 is the vortex localization length (comparable to the microstructural disorder scale), and $R(T) \propto \exp\left(\frac{E_p}{k_B T}\right)\,,$3 is a resistance prefactor approximately equal to the normal-state sheet resistance. This functional form is a direct analogue of Mott’s law for electrons, with the $R(T) \propto \exp\left(\frac{E_p}{k_B T}\right)\,,$4 exponent reflecting the two-dimensionality of vortex hopping (Percher et al., 2017).

3. Granular Modeling and Parameter Estimation

To quantitatively model VVRH transport, disordered InO$R(T) \propto \exp\left(\frac{E_p}{k_B T}\right)\,,$5 films are represented as 2D arrays of superconducting grains (diameter $R(T) \propto \exp\left(\frac{E_p}{k_B T}\right)\,,$6 nm), coupled by Josephson junctions with varying critical currents ($R(T) \propto \exp\left(\frac{E_p}{k_B T}\right)\,,$7). Vortex cores reside in the inter-grain voids, where the self-energy $R(T) \propto \exp\left(\frac{E_p}{k_B T}\right)\,,$8—dominated by the local Josephson coupling—is broadly distributed. At low vortex densities, relevant occupied pinning sites lie in the low-energy tail of this distribution,

$R(T) \propto \exp\left(\frac{E_p}{k_B T}\right)\,,$9

with $E_p$0 as the mean vortex self-energy and $E_p$1 as the dispersion. The vortex concentration is $E_p$2 with $E_p$3, and the density of states at the chemical potential is $E_p$4. This leads to a field-dependent hopping parameter,

$E_p$5

so that $E_p$6 decreases linearly with increasing $E_p$7. For typical parameters ($E_p$8 nm, $E_p$9m, $E_p/k_B$0), $E_p/k_B$1 K·T, in precise agreement with experimental results (Percher et al., 2017).

4. Experimental Validation in Amorphous InOₓ Films

Empirical tests utilized amorphous InO$E_p/k_B$2 films (thickness 55 nm) measured down to $E_p/k_B$3 mK in perpendicular magnetic fields up to 12 T. On the superconducting side ($E_p/k_B$4), the isotherms $E_p/k_B$5 cross near a critical field $E_p/k_B$6 T at resistance $E_p/k_B$7 k$E_p/k_B$8. The resistance $E_p/k_B$9 for fixed $d$0 exhibits $d$1 down to the lowest $d$2. Comprehensive analysis shows that $d$3 versus $d$4, with $d$5 varying from $d$6 down to $d$7, is minimized for $d$8 (using reduced $d$9). Plots of $$\sigma_v(T)\;\propto\;\exp\left[-\left(\frac{T_0}{T}\right)^{1/(d+1)}\right]\,,$$0 versus $$\sigma_v(T)\;\propto\;\exp\left[-\left(\frac{T_0}{T}\right)^{1/(d+1)}\right]\,,$$1 at $$\sigma_v(T)\;\propto\;\exp\left[-\left(\frac{T_0}{T}\right)^{1/(d+1)}\right]\,,$$2–$$\sigma_v(T)\;\propto\;\exp\left[-\left(\frac{T_0}{T}\right)^{1/(d+1)}\right]\,,$$3 T yield linear behavior over 1–2 decades in $$\sigma_v(T)\;\propto\;\exp\left[-\left(\frac{T_0}{T}\right)^{1/(d+1)}\right]\,,$$4, confirming the Mott VRH law. Extracted values of $$\sigma_v(T)\;\propto\;\exp\left[-\left(\frac{T_0}{T}\right)^{1/(d+1)}\right]\,,$$5 are consistent with $$\sigma_v(T)\;\propto\;\exp\left[-\left(\frac{T_0}{T}\right)^{1/(d+1)}\right]\,,$$6, with experimental $$\sigma_v(T)\;\propto\;\exp\left[-\left(\frac{T_0}{T}\right)^{1/(d+1)}\right]\,,$$7 K·T (Percher et al., 2017).

5. Phenomena on the Insulating Side and Duality Symmetry

For $$\sigma_v(T)\;\propto\;\exp\left[-\left(\frac{T_0}{T}\right)^{1/(d+1)}\right]\,,$$8, the system is insulating ($$\sigma_v(T)\;\propto\;\exp\left[-\left(\frac{T_0}{T}\right)^{1/(d+1)}\right]\,,$$9). Here, $T_0$0 is also fit by a VRH-type form—interpreted as the hopping of localized Cooper pairs or single electrons, depending on field strength. At the highest fields, the magnetoresistance changes sign, indicating pair breaking. The data obey a duality symmetry,

$T_0$1

linking the resistance behaviors on both sides of the SIT. Over $T_0$2 mK–$T_0$3 K, $T_0$4. This duality is a hallmark of field-tuned quantum phase transitions in disordered superconductors (Percher et al., 2017).

6. Open Questions and Theoretical Puzzles

Despite theoretical expectations that long-range vortex–vortex interactions should alter the Mott exponent from $T_0$5 to values between $T_0$6 and $T_0$7, experimental data consistently supports the Mott value of $T_0$8 across broad regimes. The reason for this apparent suppression of inter-vortex interactions remains unresolved. The hypothesis advanced is that, near the field-tuned SIT, charge carriers may be weakly interacting composite fermions—namely, a Cooper pair bound to a vortex—effectively screening vortex–vortex interactions and restoring the strict 2D Mott form. A comprehensive microscopic theory addressing this phenomenon is currently lacking (Percher et al., 2017).

7. Summary of Experimental and Model Parameters

Parameter	Typical Value	Physical Meaning
$T_0$9	55 nm (film thickness)	Thickness of amorphous InO$d=2$0 film
$d=2$1	50 nm	Grain diameter/localization length
$d=2$2	2.8 T	Critical magnetic field for SIT
$d=2$3	$d=2$4 k$d=2$5	Critical resistance
$d=2$6	13	Numerical prefactor in Mott VVRH law
$d=2$7	60 K·T	Slope of $d=2$8 vs. $d=2$9

The granular Josephson array model and the VVRH framework provide a coherent quantitative account of low-temperature, finite-resistance transport in disordered 2D superconductors, in close agreement with experimental measurements on InO$R(T) = R_0\,\exp\left[-(T_0/T)^{1/3}\right]\,,$0 films. The persistence of the 1/3 exponent, despite theoretical expectations to the contrary, remains an unresolved question of both experimental and theoretical significance (Percher et al., 2017).