---
title: Thermally Activated Flux Flow (TAFF)
url: https://www.emergentmind.com/topics/thermally-activated-flux-flow-taff
type: topic
---

# Thermally Activated Flux Flow (TAFF)

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 [1108.0407] [1010.0263] [2507.14540].

## 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 [1107.4170] [1108.0407]. 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$_{0.58}$Rb$_{0.42}$Fe$_{1.72}$Se$_2$ and in $2H$-NbS$_2$, the transport analysis is written in terms of a flux-bundle hopping expression that reduces, for $JHVL/T\ll 1$, to
\[
\rho(T,H)=(2\rho_cU/T)e^{-U/T},
\qquad
R(T,H)=\left(\frac{2R_cU}{T}\right)e^{-U/T},
\]
after which the prefactor is approximated as weakly temperature dependent and an Arrhenius form is recovered [1106.2283] [2507.14540]. 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
\[
U(T,H)=U_0(H)\left(1-\frac{T}{T_c}\right),
\]
which implies
\[
\ln \rho_0(H)=\ln \rho_{0f}+\frac{U_0(H)}{T_c}.
\]
In $\beta$-FeSe, the good linearity of $\ln \rho$ versus $1/T$ is taken as evidence that this approximately linear softening of the barrier toward zero at $T_c$ is adequate in the TAFF window [1108.0407].

## 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$_{0.58}$Rb$_{0.42}$Fe$_{1.72}$Se$_2$, the transition broadens strongly under field and develops a pronounced low-temperature tail; in polycrystalline Bi$_2$Sr$_2$CaCu$_2$O$_{8+\delta}$, the broadening depends strongly on grain morphology; and in $\beta$-FeSe the transition shifts to lower temperature without obvious broadening, yet the resistive tail below $T_c$ still follows TAFF behavior [1106.2283] [1203.2057] [1108.0407]. 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 $\ln \rho$ or $\ln R$ versus $1/T$ at fixed field and fit the linear segment. The slope gives $-U_0(H)$ or $-U_0(H)/k_B$, depending on convention, and the intercept gives the prefactor. Different studies choose the fitting interval differently. In Ba(Fe,Co)$_2$As$_2$, the linear TAFF regime is identified for $R(T)<0.1R_n$; in Bi-2212, the linear fitting is taken in the low-resistivity region around $10^{-4}$ to $10^{-5}\ \Omega\,\mathrm{cm}$; and several studies simply refer to the linear part of the Arrhenius plot without a fixed numerical cutoff [1107.4170] [1203.2057] [1202.1473].

A useful internal consistency check is the common-crossing construction. In $\beta$-FeSe, the extrapolated Arrhenius lines converge near $T_{\mathrm{cross}}\approx 11.1\ \mathrm{K}$, and the independently fitted relation $\ln\rho_0(H)=\ln\rho_{0f}+U_0(H)/T_c$ gives $T_c=11.2(1)\ \mathrm{K}$ for both $H\parallel(101)$ and $H\perp(101)$ [1108.0407]. Similar self-consistency appears in Fe$_{1+y}$(Te$_{1-x}$S$_x$)$_z$, where $T_{\mathrm{cross}}$ and fitted $T_c$ agree within uncertainty, in flux-free FeSe where the extrapolated lines meet near a bulk $T_c\approx 11.62\ \mathrm{K}$, and in $2H$-NbS$_2$ for $H\perp ab$, where $T_{\mathrm{cross}}=5.95\ \mathrm{K}$ is close to $T_{c,\mathrm{onset}}=6.17\ \mathrm{K}$ and the reciprocal slope of the $\ln R_0$–$U_0$ plot is $5.90\ \mathrm{K}$ [1010.0263] [1606.08201] [2507.14540].

Not all TAFF-related work is based on transport. In weak-field YBCO single crystals near $T_c$, nonlogarithmic magnetization relaxation in fields of order $1\ \mathrm{Oe}$ 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 [1105.5896].

## 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(H)\propto H^{-\alpha},
\]
or equivalently $H^{-n}$ or $B^{-\gamma}$. 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 [1107.4170] [1108.0407] [1706.08696].

| System | Field dependence of $U_0$ | Interpretation |
|---|---|---|
| $\beta$-FeSe | $H^{-0.25(6)}$ below $30\ \mathrm{kOe}$, $H^{-0.68(6)}$ above for $H\parallel(101)$; similar exponents for $H\perp(101)$ | single-vortex pinning at low field; collective creep above $30\ \mathrm{kOe}$ [1108.0407] |
| Ba(Fe,Co)$_2$As$_2$ | For $H\parallel ab$, $n=0.22$–$0.28$ up to $12\ \mathrm{T}$; for $H\parallel c$, $n\approx 0.29$–$0.36$ below $7\ \mathrm{T}$ and $0.65$–$1.19$ above | individual pinning over broad range for $H\parallel ab$; crossover to collective pinning for $H\parallel c$ above $7\ \mathrm{T}$ [1107.4170] |
| NaFe$_{1-x}$Co$_x$As | $U_0(B)\sim B^{-\alpha}$ with crossover near $2\ \mathrm{T}$; high-field exponents around $0.67$–$1.16$ depending on doping and orientation | coexistence of collective and plastic pinning below $2\ \mathrm{T}$; plastic pinning above [1706.08696] |
| quenched K$_x$Fe$_{2-y}$Se$_2$ | $U_0(H)\sim H^{-0.84(2)}$ for $H\parallel c$ and $H^{-0.78(2)}$ for $H\parallel ab$ above $3\ \mathrm{T}$ | collective flux creep in high fields; point defects as main pinning source [1110.5316] |
| $2H$-NbS$_2$ for $H\perp ab$ | $H^{-0.71}$ below $\mu_0H_{\mathrm{cr}}\sim0.7\ \mathrm{T}$ and $H^{-1.67}$ above | plastic strong pinning at low field; steeper high-field suppression attributed to an entangled vortex liquid [2507.14540] |
| CaFFe$_{0.9}$Co$_{0.1}$As | $B^{-0.05}$ for $B\le 3\ \mathrm{T}$ and $B^{-0.15}$ for $B>3\ \mathrm{T}$ | weakly pinned intergranular Josephson vortices at low field; crossover toward collective creep [1202.1473] |

The same language is used in more weakly anisotropic systems. In Tl$_{0.58}$Rb$_{0.42}$Fe$_{1.72}$Se$_2$, the activation energies at $2\ \mathrm{T}$ are about $4900\ \mathrm{K}$ for $H\perp c$ and $3607\ \mathrm{K}$ for $H\parallel c$, while the field exponents are only $0.27\pm0.02$ and $0.22\pm0.02$, respectively, which was taken as evidence of moderate pinning with only weak orientation dependence of the pinning-force scaling [1106.2283]. By contrast, polycrystalline CeFeAsO$_{1-x}$F$_x$ shows $U_0(B)\propto B^{-\gamma}$ with $\gamma\approx0.41$ to $0.50$ from $1$ to $8\ \mathrm{T}$ and little doping dependence of the exponent, while the magnitude of $U_0$ varies strongly with fluorine content [1410.6546].

## 4. Temperature dependence beyond the simple Arrhenius form

Although many TAFF analyses use $U(T,H)=U_0(H)(1-T/T_c)$, several systems show that this is not universally adequate. A widely used generalization is
\[
U(T,H)=U_0(H)\left(1-\frac{T}{T_c}\right)^q,
\]
which leads to
\[
\ln R = \ln(2R_cU_0) + q\ln\left(1-\frac{T}{T_c}\right)-\ln T-\frac{U_0\left(1-\frac{T}{T_c}\right)^q}{T}.
\]
This form retains the temperature dependence of the prefactor and allows non-Arrhenius curvature in $\ln R$ versus $1/T$ [1010.0263] [1706.08696] [2507.14540].

Fe$_{1+y}$(Te$_{1-x}$S$_x$)$_z$ 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 $q\approx1$ for $H\parallel c$ and $q\approx2$ for $H\parallel ab$, while $U_0(H)$ still shows a low-field to high-field crossover near $5\ \mathrm{T}$ from single-vortex pinning to collective creep [1010.0263].

NaFe$_{1-x}$Co$_x$As provides a different pattern. Across underdoped, optimally doped, and overdoped single crystals, and for both $B\parallel c$ and $B\parallel ab$, the TAFF region is non-Arrhenius but is well described by
\[
U(T,B)=U_0(B)\left(1-\frac{T}{T_c}\right)^{1.5}.
\]
The universal value $q=1.5$ 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 [1706.08696].

The strongest anisotropic breakdown of simple Arrhenius TAFF in the dataset appears in $2H$-NbS$_2$. For $H\perp ab$, $\ln R$ versus $1/T$ is linear in the TAFF region and the standard Arrhenius construction works. For $H\parallel ab$, 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 $q=2$. 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 [2507.14540].

By contrast, several materials remain close to the linear-barrier limit. In $\beta$-FeSe, the good Arrhenius linearity supports $U(T,H)=U_0(H)(1-T/T_c)$; in Bi-2212, the fitted temperature exponent is $\beta=1$ for all samples and fields; and in quenched K$_x$Fe$_{2-y}$Se$_2$, the Arrhenius linearity is again taken to indicate an approximately linear barrier softening [1108.0407] [1203.2057] [1110.5316].

## 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 $\beta$-FeSe single crystals, the Arrhenius fits are equally good for $H\parallel(101)$ and $H\perp(101)$, the extracted $U_0(H)$ values are similar, and the exponents on either side of the $30\ \mathrm{kOe}$ 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 $J_c$ much smaller than in several other iron-based systems [1108.0407].

Other compounds are markedly anisotropic. In Ba(Fe,Co)$_2$As$_2$, $U_0$ is larger for $H\parallel ab$ than for $H\parallel c$ in the optimally doped and overdoped crystals, and the field dependence is much weaker for $H\parallel ab$. The optimally doped sample has the highest $U_0$, while increasing Co content reduces the high-field $H\parallel c$ exponent from $1.19$ to $0.65$, which the authors interpret as evidence that Co-induced disorder strengthens individual pinning and competes with collective bundle pinning [1107.4170]. In Tl$_{0.58}$Rb$_{0.42}$Fe$_{1.72}$Se$_2$, by contrast, the magnitude of $U_0$ differs between $H\perp c$ and $H\parallel c$, but the field exponents remain close and the pinning-force scaling is only weakly anisotropic [1106.2283].

Microstructure is equally consequential. In sol-gel Bi-2212, increasing sintering temperature from $760^\circ$C to $840^\circ$C improves grain connectivity, raises $T_c^{R=0}$ from $31.74\ \mathrm{K}$ to $81.79\ \mathrm{K}$, narrows $\Delta T_c$ from $58.01\ \mathrm{K}$ to $14.61\ \mathrm{K}$, and increases the field exponent from about $0.18$–$0.19$ to $0.45$. The study interprets this as stronger flux pinning in better-coupled samples and suggests that the highest-temperature sample, with $\alpha\approx0.45$, may approach planar-defect-like collective pinning [1203.2057]. In polycrystalline CaFFe$_{0.9}$Co$_{0.1}$As, 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 $J_c$ and broad field-broadened resistive transition [1202.1473].

A particularly direct connection between TAFF and current-carrying capability appears in CeFeAsO$_{1-x}$F$_x$. There, $U_0(1\,\mathrm{T})$ and $J_c(1\,\mathrm{T})$ both peak in the mid-doping region $x=0.15$–$0.20$, and $U_0(1\,\mathrm{T})$ shows an excellent linear correlation with self-field $J_c$ at $5\ \mathrm{K}$. 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 $J_c$ [1410.6546].

## 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 $T_v$ 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 $V\propto I^\beta$, and the resistance is exponentially dependent on $B^2$ rather than on a field-induced vortex density [1510.03259]. 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 [1006.2964]. 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$_{0.58}$Rb$_{0.42}$Fe$_{1.72}$Se$_2$ 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 [1106.2283]. Fe$_{1+y}$(Te$_{1-x}$S$_x$)$_z$ shows that a constant prefactor and a strictly linear barrier can be quantitatively inadequate, so the extracted $U_0$ from a straight-line fit is better regarded as an apparent activation energy [1010.0263]. Flux-free FeSe adds a materials caveat: the transport-based TAFF analysis is internally consistent, but the crystal contains both majority tetragonal $\beta$-FeSe and minority magnetic hexagonal $\delta$-FeSe, so the microscopic origin of the pinning landscape cannot be attributed to phase-pure FeSe alone [1606.08201].

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.

Source: https://www.emergentmind.com/topics/thermally-activated-flux-flow-taff