Thermally Activated Flux Flow (TAFF)

• TAFF is the mixed-state transport regime in type-II superconductors where thermally activated vortex motion leads to finite resistance below the superconducting transition.
• TAFF analyses typically use Arrhenius and modified Arrhenius plots of ln(ρ) versus 1/T to extract a field-dependent activation energy U₀(H) that characterizes vortex dynamics.
• TAFF studies reveal material-dependent factors such as anisotropy, disorder, and microstructure, which impact vortex pinning, collective creep, and critical current behavior.

Searching arXiv for recent and foundational TAFF papers to support the article. arxiv_search(query="thermally activated flux flow superconductor TAFF", max_results=10) Thermally activated flux flow (TAFF) is the mixed-state transport regime of a type-II superconductor in which vortices, or vortex bundles, overcome pinning barriers by thermal activation and thereby generate a finite resistance below the superconducting transition. In the literature, TAFF is identified primarily from field-broadened resistive tails and from Arrhenius or modified Arrhenius analyses of $\ln \rho$ or $\ln R$ versus $1/T$. The central output of such analyses is a field-dependent activation-energy scale $U_0(H)$, whose magnitude, temperature dependence, anisotropy, and field scaling are used to diagnose vortex pinning, collective creep, plastic flow, dimensionality, and the influence of disorder or microstructure (Lei et al., 2011, Lei et al., 2010, Wang et al., 19 Jul 2025).

1. Phenomenological basis and canonical equations

In the standard transport formulation, TAFF is described by an activated resistivity or resistance,

$$\rho(T,H)=\rho_0(H)\exp\!\left[-\frac{U_0(H)}{k_B T}\right], \qquad R(T,H)=R_0\exp\!\left(-\frac{U_0}{k_B T}\right),$$

or, in the notation used by several iron-based-superconductor papers, by

$\ln \rho(T,H)=\ln \rho_0(H)-\frac{U_0(H)}{T},$

with $U_0$ then expressed effectively in kelvin units rather than with an explicit $k_B$ factor (Vinod et al., 2011, Lei et al., 2011). The practical distinction is not conceptual but metrological: some studies report $U_0$ in meV or eV and keep $k_B$ explicit, whereas others absorb it into the definition of the apparent activation energy.

A more microscopic route starts from the low-current limit of stochastic vortex hopping. In Tl$\ln R$0Rb$\ln R$1Fe$\ln R$2Se$\ln R$3 and in $\ln R$4-NbS$\ln R$5, the transport analysis is written in terms of a flux-bundle hopping expression that reduces, for $\ln R$6, to

$\ln R$7

after which the prefactor is approximated as weakly temperature dependent and an Arrhenius form is recovered (Jiao et al., 2011, Wang et al., 19 Jul 2025). In this sense, TAFF is the low-current, ohmic limit of thermally assisted vortex motion in a pinned mixed state.

Several works further parametrize the barrier itself as

$\ln R$8

which implies

$\ln R$9

In $1/T$0-FeSe, the good linearity of $1/T$1 versus $1/T$2 is taken as evidence that this approximately linear softening of the barrier toward zero at $1/T$3 is adequate in the TAFF window (Lei et al., 2011).

2. Experimental identification and extraction procedures

Experimentally, TAFF is usually inferred from the field-broadened low-resistivity tail of the superconducting transition. In single-crystalline Tl$1/T$4Rb$1/T$5Fe$1/T$6Se$1/T$7, the transition broadens strongly under field and develops a pronounced low-temperature tail; in polycrystalline Bi$1/T$8Sr$1/T$9CaCu$U_0(H)$0O$U_0(H)$1, the broadening depends strongly on grain morphology; and in $U_0(H)$2-FeSe the transition shifts to lower temperature without obvious broadening, yet the resistive tail below $U_0(H)$3 still follows TAFF behavior (Jiao et al., 2011, Sharma et al., 2012, Lei et al., 2011). The essential signature is therefore not merely a shift of the transition, but the persistence of a finite-resistivity tail in the mixed state.

The standard extraction protocol is to plot $U_0(H)$4 or $U_0(H)$5 versus $U_0(H)$6 at fixed field and fit the linear segment. The slope gives $U_0(H)$7 or $U_0(H)$8, depending on convention, and the intercept gives the prefactor. Different studies choose the fitting interval differently. In Ba(Fe,Co)$U_0(H)$9As$$\rho(T,H)=\rho_0(H)\exp\!\left[-\frac{U_0(H)}{k_B T}\right], \qquad R(T,H)=R_0\exp\!\left(-\frac{U_0}{k_B T}\right),$$0, the linear TAFF regime is identified for $$\rho(T,H)=\rho_0(H)\exp\!\left[-\frac{U_0(H)}{k_B T}\right], \qquad R(T,H)=R_0\exp\!\left(-\frac{U_0}{k_B T}\right),$$1; in Bi-2212, the linear fitting is taken in the low-resistivity region around $$\rho(T,H)=\rho_0(H)\exp\!\left[-\frac{U_0(H)}{k_B T}\right], \qquad R(T,H)=R_0\exp\!\left(-\frac{U_0}{k_B T}\right),$$2 to $$\rho(T,H)=\rho_0(H)\exp\!\left[-\frac{U_0(H)}{k_B T}\right], \qquad R(T,H)=R_0\exp\!\left(-\frac{U_0}{k_B T}\right),$$3; and several studies simply refer to the linear part of the Arrhenius plot without a fixed numerical cutoff (Vinod et al., 2011, Sharma et al., 2012, Shekhar et al., 2012).

A useful internal consistency check is the common-crossing construction. In $$\rho(T,H)=\rho_0(H)\exp\!\left[-\frac{U_0(H)}{k_B T}\right], \qquad R(T,H)=R_0\exp\!\left(-\frac{U_0}{k_B T}\right),$$4-FeSe, the extrapolated Arrhenius lines converge near $$\rho(T,H)=\rho_0(H)\exp\!\left[-\frac{U_0(H)}{k_B T}\right], \qquad R(T,H)=R_0\exp\!\left(-\frac{U_0}{k_B T}\right),$$5, and the independently fitted relation $$\rho(T,H)=\rho_0(H)\exp\!\left[-\frac{U_0(H)}{k_B T}\right], \qquad R(T,H)=R_0\exp\!\left(-\frac{U_0}{k_B T}\right),$$6 gives $$\rho(T,H)=\rho_0(H)\exp\!\left[-\frac{U_0(H)}{k_B T}\right], \qquad R(T,H)=R_0\exp\!\left(-\frac{U_0}{k_B T}\right),$$7 for both $$\rho(T,H)=\rho_0(H)\exp\!\left[-\frac{U_0(H)}{k_B T}\right], \qquad R(T,H)=R_0\exp\!\left(-\frac{U_0}{k_B T}\right),$$8 and $$\rho(T,H)=\rho_0(H)\exp\!\left[-\frac{U_0(H)}{k_B T}\right], \qquad R(T,H)=R_0\exp\!\left(-\frac{U_0}{k_B T}\right),$$9 (Lei et al., 2011). Similar self-consistency appears in Fe$\ln \rho(T,H)=\ln \rho_0(H)-\frac{U_0(H)}{T},$0(Te$\ln \rho(T,H)=\ln \rho_0(H)-\frac{U_0(H)}{T},$1S$\ln \rho(T,H)=\ln \rho_0(H)-\frac{U_0(H)}{T},$2)$\ln \rho(T,H)=\ln \rho_0(H)-\frac{U_0(H)}{T},$3, where $\ln \rho(T,H)=\ln \rho_0(H)-\frac{U_0(H)}{T},$4 and fitted $\ln \rho(T,H)=\ln \rho_0(H)-\frac{U_0(H)}{T},$5 agree within uncertainty, in flux-free FeSe where the extrapolated lines meet near a bulk $\ln \rho(T,H)=\ln \rho_0(H)-\frac{U_0(H)}{T},$6, and in $\ln \rho(T,H)=\ln \rho_0(H)-\frac{U_0(H)}{T},$7-NbS$\ln \rho(T,H)=\ln \rho_0(H)-\frac{U_0(H)}{T},$8 for $\ln \rho(T,H)=\ln \rho_0(H)-\frac{U_0(H)}{T},$9, where $U_0$0 is close to $U_0$1 and the reciprocal slope of the $U_0$2–$U_0$3 plot is $U_0$4 (Lei et al., 2010, Maheshwari et al., 2016, Wang et al., 19 Jul 2025).

Not all TAFF-related work is based on transport. In weak-field YBCO single crystals near $U_0$5, nonlogarithmic magnetization relaxation in fields of order $U_0$6 was interpreted as thermoactivated creep with dynamics “similar to thermally assisted flux flow,” but this was explicitly presented as TAFF-like magnetization relaxation rather than canonical resistive TAFF transport (Monarkha et al., 2011).

3. Field scaling of activation energy and vortex-dynamics crossovers

The most common quantitative output of TAFF analysis is a field-dependent activation energy of the form

$U_0$7

or equivalently $U_0$8 or $U_0$9. Small exponents correspond to weak field dependence and are typically associated with single-vortex or individual pinning, whereas larger exponents indicate stronger many-vortex effects such as collective creep, collective pinning, or plastic motion. This interpretation is explicit in several iron-based systems and recurs, with material-specific variations, in cuprates and layered dichalcogenides (Vinod et al., 2011, Lei et al., 2011, Choi et al., 2017).

System	Field dependence of $k_B$0	Interpretation
$k_B$1-FeSe	$k_B$2 below $k_B$3, $k_B$4 above for $k_B$5; similar exponents for $k_B$6	single-vortex pinning at low field; collective creep above $k_B$7 (Lei et al., 2011)
Ba(Fe,Co)$k_B$8As$k_B$9	For $U_0$0, $U_0$1–$U_0$2 up to $U_0$3; for $U_0$4, $U_0$5–$U_0$6 below $U_0$7 and $U_0$8–$U_0$9 above	individual pinning over broad range for $k_B$0; crossover to collective pinning for $k_B$1 above $k_B$2 (Vinod et al., 2011)
NaFe$k_B$3Co$k_B$4As	$k_B$5 with crossover near $k_B$6; high-field exponents around $k_B$7–$k_B$8 depending on doping and orientation	coexistence of collective and plastic pinning below $k_B$9; plastic pinning above (Choi et al., 2017)
quenched K$\ln R$00Fe$\ln R$01Se$\ln R$02	$\ln R$03 for $\ln R$04 and $\ln R$05 for $\ln R$06 above $\ln R$07	collective flux creep in high fields; point defects as main pinning source (Lei et al., 2011)
$\ln R$08-NbS$\ln R$09 for $\ln R$10	$\ln R$11 below $\ln R$12 and $\ln R$13 above	plastic strong pinning at low field; steeper high-field suppression attributed to an entangled vortex liquid (Wang et al., 19 Jul 2025)
CaFFe$\ln R$14Co$\ln R$15As	$\ln R$16 for $\ln R$17 and $\ln R$18 for $\ln R$19	weakly pinned intergranular Josephson vortices at low field; crossover toward collective creep (Shekhar et al., 2012)

The same language is used in more weakly anisotropic systems. In Tl$\ln R$20Rb$\ln R$21Fe$\ln R$22Se$\ln R$23, the activation energies at $\ln R$24 are about $\ln R$25 for $\ln R$26 and $\ln R$27 for $\ln R$28, while the field exponents are only $\ln R$29 and $\ln R$30, respectively, which was taken as evidence of moderate pinning with only weak orientation dependence of the pinning-force scaling (Jiao et al., 2011). By contrast, polycrystalline CeFeAsO$\ln R$31F$\ln R$32 shows $\ln R$33 with $\ln R$34 to $\ln R$35 from $\ln R$36 to $\ln R$37 and little doping dependence of the exponent, while the magnitude of $\ln R$38 varies strongly with fluorine content (Chong et al., 2014).

4. Temperature dependence beyond the simple Arrhenius form

Although many TAFF analyses use $\ln R$39, several systems show that this is not universally adequate. A widely used generalization is

$\ln R$40

which leads to

$\ln R$41

This form retains the temperature dependence of the prefactor and allows non-Arrhenius curvature in $\ln R$42 versus $\ln R$43 (Lei et al., 2010, Choi et al., 2017, Wang et al., 19 Jul 2025).

Fe$\ln R$44(Te$\ln R$45S$\ln R$46)$\ln R$47 is a particularly clear case. There, conventional Arrhenius analysis gives only an approximate average barrier over a narrow window, whereas the modified TAFF theory fits the data more faithfully. The fitted exponent is strongly orientation dependent, with $\ln R$48 for $\ln R$49 and $\ln R$50 for $\ln R$51, while $\ln R$52 still shows a low-field to high-field crossover near $\ln R$53 from single-vortex pinning to collective creep (Lei et al., 2010).

NaFe$\ln R$54Co$\ln R$55As provides a different pattern. Across underdoped, optimally doped, and overdoped single crystals, and for both $\ln R$56 and $\ln R$57, the TAFF region is non-Arrhenius but is well described by

$\ln R$58

The universal value $\ln R$59 is interpreted in that work as evidence of 3D vortex behavior, and the same study further argues that the low-temperature phase evolves into a 3D vortex glass (Choi et al., 2017).

The strongest anisotropic breakdown of simple Arrhenius TAFF in the dataset appears in $\ln R$60-NbS$\ln R$61. For $\ln R$62, $\ln R$63 versus $\ln R$64 is linear in the TAFF region and the standard Arrhenius construction works. For $\ln R$65, however, the lines do not converge to a common crossing point, the standard Arrhenius form fails, and the modified TAFF theory gives the best fit with $\ln R$66. The authors interpret this as potentially revealing 2D characteristics of vortices in the TAFF region, and the field dependence changes qualitatively from a power law to a Kramer-like parabolic form (Wang et al., 19 Jul 2025).

By contrast, several materials remain close to the linear-barrier limit. In $\ln R$67-FeSe, the good Arrhenius linearity supports $\ln R$68; in Bi-2212, the fitted temperature exponent is $\ln R$69 for all samples and fields; and in quenched K$\ln R$70Fe$\ln R$71Se$\ln R$72, the Arrhenius linearity is again taken to indicate an approximately linear barrier softening (Lei et al., 2011, Sharma et al., 2012, Lei et al., 2011).

5. Anisotropy, disorder, microstructure, and correlation with critical current

TAFF is highly sensitive to anisotropy, but the form of that sensitivity is material dependent. In tetragonal $\ln R$73-FeSe single crystals, the Arrhenius fits are equally good for $\ln R$74 and $\ln R$75, the extracted $\ln R$76 values are similar, and the exponents on either side of the $\ln R$77 crossover are essentially the same within uncertainty. This near-isotropy of TAFF tracks the nearly isotropic upper critical fields and the broader picture of FeSe as a Pauli-limited type-II superconductor with relatively modest pinning strength and $\ln R$78 much smaller than in several other iron-based systems (Lei et al., 2011).

Other compounds are markedly anisotropic. In Ba(Fe,Co)$\ln R$79As$\ln R$80, $\ln R$81 is larger for $\ln R$82 than for $\ln R$83 in the optimally doped and overdoped crystals, and the field dependence is much weaker for $\ln R$84. The optimally doped sample has the highest $\ln R$85, while increasing Co content reduces the high-field $\ln R$86 exponent from $\ln R$87 to $\ln R$88, which the authors interpret as evidence that Co-induced disorder strengthens individual pinning and competes with collective bundle pinning (Vinod et al., 2011). In Tl$\ln R$89Rb$\ln R$90Fe$\ln R$91Se$\ln R$92, by contrast, the magnitude of $\ln R$93 differs between $\ln R$94 and $\ln R$95, but the field exponents remain close and the pinning-force scaling is only weakly anisotropic (Jiao et al., 2011).

Microstructure is equally consequential. In sol-gel Bi-2212, increasing sintering temperature from $\ln R$96C to $\ln R$97C improves grain connectivity, raises $\ln R$98 from $\ln R$99 to $1/T$00, narrows $1/T$01 from $1/T$02 to $1/T$03, and increases the field exponent from about $1/T$04–$1/T$05 to $1/T$06. The study interprets this as stronger flux pinning in better-coupled samples and suggests that the highest-temperature sample, with $1/T$07, may approach planar-defect-like collective pinning (Sharma et al., 2012). In polycrystalline CaFFe$1/T$08Co$1/T$09As, TAFF is instead attributed primarily to weakly linked granules, amorphous grain boundaries, FeAs impurity regions, and weakly pinned intergranular Josephson vortices; the same defects are invoked to explain the low $1/T$10 and broad field-broadened resistive transition (Shekhar et al., 2012).

A particularly direct connection between TAFF and current-carrying capability appears in CeFeAsO$1/T$11F$1/T$12. There, $1/T$13 and $1/T$14 both peak in the mid-doping region $1/T$15–$1/T$16, and $1/T$17 shows an excellent linear correlation with self-field $1/T$18 at $1/T$19. The authors interpret this within a two-fluid flux-creep model for granular samples, arguing that the doping dependence of the depinning barrier governs the doping dependence of $1/T$20 (Chong et al., 2014).

6. Extensions, special geometries, and methodological caveats

TAFF is not confined to the standard geometry of transport in a perpendicular magnetic field through a bulk or single-crystalline mixed state. In a thin amorphous MoGe film with the magnetic field applied parallel to the film plane, a resistive state appears above a characteristic temperature $1/T$21 and is attributed to thermally activated hopping of spontaneous perpendicular vortices tuned by the pairbreaking effect of the parallel field. The low-current response is Ohmic, the higher-current response follows $1/T$22, and the resistance is exponentially dependent on $1/T$23 rather than on a field-induced vortex density (Kunchur et al., 2015). This suggests that the TAFF concept extends naturally to thin-film, edge-controlled spontaneous-vortex transport, provided the activated object remains a vortex crossing a barrier.

Theoretical work on two-dimensional vortex motion in an asymmetric washboard planar pinning potential reaches a complementary conclusion. In that exactly solvable Fokker–Planck treatment, the linear regimes of motion are the thermally activated flux-flow regime and the ohmic flux-flow regime, and the odd longitudinal and transverse voltages produced by ratchet asymmetry vanish in both of those limits. They appear only in the highly nonlinear crossover from TAFF to FF, where the barrier is sampled asymmetrically (Shklovskij et al., 2010). This places TAFF within a broader nonlinear transport framework rather than treating it as an isolated fitting ansatz.

Several studies also emphasize limitations of the simplest Arrhenius analysis. Tl$1/T$24Rb$1/T$25Fe$1/T$26Se$1/T$27 explicitly notes that resistive broadening can arise from both TAFF and superconducting critical fluctuations, even though the Arrhenius fits work well over a broad temperature range (Jiao et al., 2011). Fe$1/T$28(Te$1/T$29S$1/T$30)$1/T$31 shows that a constant prefactor and a strictly linear barrier can be quantitatively inadequate, so the extracted $1/T$32 from a straight-line fit is better regarded as an apparent activation energy (Lei et al., 2010). Flux-free FeSe adds a materials caveat: the transport-based TAFF analysis is internally consistent, but the crystal contains both majority tetragonal $1/T$33-FeSe and minority magnetic hexagonal $1/T$34-FeSe, so the microscopic origin of the pinning landscape cannot be attributed to phase-pure FeSe alone (Maheshwari et al., 2016).

Taken together, these studies define TAFF not as a single universal formula but as a family of activated vortex-transport regimes. The unifying feature is thermally assisted vortex motion below the superconducting transition; the nonuniversal features are the barrier form, prefactor, dimensionality exponent, field scaling, and the relative importance of single-vortex pinning, collective creep, plastic deformation, weak-link transport, or edge-controlled spontaneous-vortex hopping.