Flux-Flow Instability in Superconductors

• Flux-flow instability (FFI) is the loss of stability in the moving vortex state of superconductors, marked by abrupt jumps in the current–voltage characteristic at a critical vortex velocity.
• The underlying mechanisms include the Larkin–Ovchinnikov effect, thermal runaway, and edge-barrier influences, with each contributing to changes in vortex viscosity and effective current margins.
• FFI is pivotal both as a practical quench mechanism in superconducting devices and as a diagnostic tool for probing nonequilibrium quasiparticle relaxation and heat flow dynamics.

Flux-flow instability (FFI) is the loss of stability of the moving-vortex state in a superconductor. In the mixed state of a type-II superconductor, a transport current drives vortices into a dissipative flux-flow regime; at sufficiently large drive, however, the vortex system can no longer sustain steady motion and the current–voltage characteristic exhibits an abrupt jump to a highly resistive or normal-state branch at a critical vortex velocity $v^\ast$ (Leo et al., 2017). Across conventional thin films, iron-based superconductors, disordered two-dimensional superconductors, and layered Josephson systems, FFI has become both a practical quench mechanism and a diagnostic of nonequilibrium quasiparticle relaxation, self-heating, pinning, edge barriers, and local vortex dynamics (Dobrovolskiy, 2023).

1. Definition and phenomenology

In the standard Abrikosov-vortex setting, the Lorentz force per unit length is

$\mathbf{f}_L=\mathbf{J}\times \Phi_0,$

and vortex motion generates an electric field

$\mathbf{E}=\mathbf{B}\times \mathbf{v},$

so that $E=Bv$ in the usual orthogonal geometry of current, field, and vortex velocity (Leo et al., 2017). The resulting dissipative state is the flux-flow regime, often described by a roughly linear or quasi-linear branch of the current–voltage characteristic.

Experimentally, FFI is identified by three observables. The first is the onset current of dissipation, denoted either the critical current $I_c$ or, in depinning language, the depinning current density $j_d$. The second is the instability or quenching current $I^\ast$ (or $j^\ast$), defined by the voltage jump. The third is the instability voltage $V^\ast$, from which the critical vortex velocity is inferred as

$v^\ast=\frac{V^\ast}{LB},$

with $\mathbf{f}_L=\mathbf{J}\times \Phi_0,$0 the voltage-contact separation (Leo et al., 2017). Between $\mathbf{f}_L=\mathbf{J}\times \Phi_0,$1 and $\mathbf{f}_L=\mathbf{J}\times \Phi_0,$2, the superconductor is already dissipative but has not yet quenched.

A useful operational quantity is the interval $\mathbf{f}_L=\mathbf{J}\times \Phi_0,$3. In Fe(Se,Te) microbridges, this difference was explicitly interpreted as a “safe range” before complete quench, and normalized as

$\mathbf{f}_L=\mathbf{J}\times \Phi_0,$4

That definition makes FFI directly relevant to device margins rather than only to vortex kinetics (Leo et al., 2017).

2. Principal physical mechanisms

The canonical electronic mechanism is the Larkin–Ovchinnikov (LO) instability. In that picture, fast vortex motion drives quasiparticles in vortex cores far from equilibrium; the vortex viscosity becomes velocity dependent,

$\mathbf{f}_L=\mathbf{J}\times \Phi_0,$5

so the viscous drag reaches a maximum and then decreases, producing a runaway instability once $\mathbf{f}_L=\mathbf{J}\times \Phi_0,$6 approaches $\mathbf{f}_L=\mathbf{J}\times \Phi_0,$7 (Bezuglyj et al., 2019). Near $\mathbf{f}_L=\mathbf{J}\times \Phi_0,$8, this is associated with nonequilibrium quasiparticles escaping the core and an effective core shrinkage; in the dirty-limit LO formulation, $\mathbf{f}_L=\mathbf{J}\times \Phi_0,$9 is controlled by the quasiparticle energy relaxation time $\mathbf{E}=\mathbf{B}\times \mathbf{v},$0 (Dobrovolskiy, 2023).

At lower temperatures, a distinct hot-electron regime becomes important. In the Kunchur framework, strong electron–electron scattering drives the quasiparticle subsystem toward a thermal-like distribution at an elevated electron temperature $\mathbf{E}=\mathbf{B}\times \mathbf{v},$1, and the vortex core broadens rather than shrinks. The instability then reflects the inability of electron–phonon cooling to remove the injected power rapidly enough (Baeva et al., 18 Sep 2025). In epitaxial TiN, comparison with independent relaxation measurements indicates that the dominant microscopic process is an increase in quasiparticle temperature relative to phonons, not quasiparticle escape from the core (Baeva et al., 18 Sep 2025).

A third class of mechanisms is extrinsic thermal runaway. Joule heating in the flux-flow branch can raise the local temperature, depress the critical current, and trigger an abrupt transition that may resemble an electronic FFI in the current–voltage curve. The Bezuglyj–Shklovskij (BS) framework was developed precisely to describe how finite heat removal modifies LO-type instability and introduces a crossover between electronic and heating-dominated regimes (Shklovskij et al., 2017).

Magnetothermal instability can also be formulated at a more macroscopic level. Taylanov and Samadov analyzed the coupled evolution of magnetic induction, temperature, and vortex velocity in the flux-flow regime, including an effective vortex mass. In their linearized treatment,

$\mathbf{E}=\mathbf{B}\times \mathbf{v},$2

with $\mathbf{E}=\mathbf{B}\times \mathbf{v},$3 a dimensionless inertia parameter and $\mathbf{E}=\mathbf{B}\times \mathbf{v},$4 a magnetothermal parameter. They found instability for $\mathbf{E}=\mathbf{B}\times \mathbf{v},$5, and for the special case $\mathbf{E}=\mathbf{B}\times \mathbf{v},$6 a thermal criterion $\mathbf{E}=\mathbf{B}\times \mathbf{v},$7 (Taylanov et al., 2014). This is a different FFI limit from LO: it emphasizes avalanche-like flux redistribution and oscillatory dynamics rather than quasiparticle-core nonequilibrium.

3. Theoretical descriptions and scaling relations

The LO and BS theories remain the standard starting point for FFI near $\mathbf{E}=\mathbf{B}\times \mathbf{v},$8. In the LO regime, the critical velocity scales as

$\mathbf{E}=\mathbf{B}\times \mathbf{v},$9

while the characteristic electric field is $E=Bv$0 (Bezuglyj et al., 2019). Doettinger-type corrections account for the field dependence that appears once the quasiparticle diffusion length $E=Bv$1 becomes comparable to the intervortex spacing; this produces the frequently observed low-field trend $E=Bv$2 (Dobrovolskiy, 2023).

The BS generalization introduces a characteristic crossover field $E=Bv$3. For $E=Bv$4, instability is predominantly non-thermal and LO-like; for $E=Bv$5, self-heating dominates. In the near-$E=Bv$6 theory with finite heat removal, the critical power density $E=Bv$7 becomes a central quantity because it directly links the instability to substrate cooling and film thickness (Shklovskij et al., 2017). In Fe(Se,Te) microbridges, analysis of $E=Bv$8 gave $E=Bv$9 for sample W10 at 15 K, far above the experimental field range up to 2.5 T, indicating an electronic-dominated regime throughout those measurements (Leo et al., 2017).

A major refinement of recent years is the distinction between global and local FFI. Global models assume homogeneous vortex motion and a uniform nonequilibrium state. Local models assume that instability nucleates only in a narrow strip or channel where vortices move faster than elsewhere. Bezuglyj and co-workers showed that such a local instability can explain abrupt jumps terminating nearly linear global $I_c$0–$I_c$1 curves: the nonlinearity is confined to a small region, while the rest of the sample remains in ordinary flux flow (Bezuglyj et al., 2019). This local formulation has been especially important for samples with strong edge barriers or inhomogeneous pinning (Dobrovolskiy, 2023).

TiN provides a further extension of theory into geometry-dependent nonequilibrium. In 12 nm TiN strips with negligible volume pinning, the measured $I_c$2 could not be reconciled with a pure LO interpretation of $I_c$3. A modified Kunchur relation, including a width-dependent factor $I_c$4 that parameterizes the crossover from effectively one-dimensional to two-dimensional quasiparticle cooling, yielded $I_c$5 consistent with independently measured $I_c$6 (Baeva et al., 18 Sep 2025). This established FFI as a quantitative relaxation probe, provided the transport geometry and cooling dimensionality are explicitly included.

4. Pinning, geometry, and edge barriers

Pinning alters FFI in a highly structured way. In epitaxial Nb films at $I_c$7, the depinning current density followed

$I_c$8

with $I_c$9 for weak random pinning, $j_d$0 for strong random pinning, and $j_d$1 for strong periodic pinning induced by nanogrooves (Dobrovolskiy et al., 2017). The instability current $j_d$2 changed comparatively little, whereas the instability velocity $j_d$3 was extremely sensitive to the pinning landscape: monotonic $j_d$4-like decrease in weakly pinned films, a broad maximum at low fields in strongly pinned films, and an additional sharp maximum at the matching field $j_d$5 mT in nanopatterned films due to commensurability (Dobrovolskiy et al., 2017).

Near $j_d$6, the effect of pinning can be incorporated into the BS framework by replacing the free-flux-flow branch with a nonlinear current–voltage characteristic appropriate for a washboard pinning potential. In that theory, stronger pinning increases $j_d$7 but decreases $j_d$8 and $j_d$9, while $I^\ast$0 and the effective instability temperature $I^\ast$1 remain practically constant (Shklovskij et al., 2017). This result is conceptually important: pinning changes how the instability point is reached in current–field space without strongly changing the critical dissipated power.

Geometry controls the thermal channel as strongly as it controls pinning. In Fe(Se,Te) thin films, a 1 mm-wide bridge displayed no clear FFI signature, whereas narrower 20 $I^\ast$2m and 10 $I^\ast$3m bridges showed distinct voltage jumps. The interpretation was that reducing width increases the film–substrate heat-transfer coefficient $I^\ast$4, suppresses self-heating, and allows the intrinsic electronic instability to manifest; the 10 $I^\ast$5m bridge displayed the steepest jump and the clearest LO-like behavior (Leo et al., 2017).

Edge quality is a further control parameter. In 15 nm MoSi strips, improving the edges by focused ion beam milling produced a factor of 3 larger $I^\ast$6, a factor of 20 higher maximal vortex velocities of 20 km/s, and a factor of 40 shorter apparent $I^\ast$7 than in otherwise identical laser-etched strips with rough, oxygen-enriched edges (Budinska et al., 2021). The accompanying message was methodological: intrinsic $I^\ast$8 can be inferred from FFI only when the field dependence of $I^\ast$9 demonstrates dominant edge-barrier control rather than disorder-dominated local vortex entry (Budinska et al., 2021).

In directly written Nb–C microstrips, a close-to-perfect edge barrier organized the vortex dynamics up to $j^\ast$0–$j^\ast$1 km/s. The spatial development of FFI there was described by an edge-controlled model in which a chain of nucleation points forms along the edge and evolves into self-organized Josephson-like junctions or “vortex rivers” (Dobrovolskiy et al., 2020). This explicitly local picture generalizes the earlier narrow-strip instability theory and connects FFI to nonuniform current crowding at the sample boundary.

5. Manifestations across materials and dimensionalities

In iron chalcogenides, FFI displays a distinctive relation between $j^\ast$2 and $j^\ast$3. In Fe(Se,Te) microbridges, $j^\ast$4 varied only weakly with field—always by less than 24% over 0–2.5 T and by about 3% in the W10 bridge at 15 K—whereas $j^\ast$5 changed by more than 60% over the same range. As a result, $j^\ast$6 increased at low field and then saturated at a crossover field around 1 T (Leo et al., 2017). Comparison with YBa$j^\ast$7Cu$j^\ast$8O$j^\ast$9, where FFI is dominated by thermal runaway, showed the opposite low-field trend of $V^\ast$0, while Nb films with pure electronic instability resembled the Fe(Se,Te) behavior (Leo et al., 2017). This suggests that $V^\ast$1 can serve as an empirical discriminator between predominantly intrinsic and predominantly thermal FFI, although the same study explicitly identified that proposition as an open question rather than a settled universal rule.

Two-dimensional superconductors add BKT-specific phenomenology. In AlO$V^\ast$2/KTaO$V^\ast$3 (111), the high-current dissipation channels below $V^\ast$4 were hot-spots and LO-type FFI, and vortex velocities inferred from the Gor’kov–Josephson relation reached $V^\ast$5 m/s near the instability (Ojha et al., 2022). The field dependence of

$V^\ast$6

followed the expected $V^\ast$7 trend at higher fields, consistent with generalized LO theory including heating. A particularly unusual observation was clockwise hysteresis in the interval $V^\ast$8, interpreted as a signature of electronic LO-type instability rather than of hot-spot-dominated switching (Ojha et al., 2022).

Weak-pinning TiN illustrates the relaxation-spectroscopy role of FFI. In 12 nm strips, narrow devices showed no extended flux-flow branch because the critical current density was so high that vortices effectively entered directly into the unstable regime, while wider strips exhibited a broad interval $V^\ast$9 from which $v^\ast=\frac{V^\ast}{LB},$0 could be extracted reliably (Baeva et al., 18 Sep 2025). The width dependence itself carried physical information, because it reflected the crossover in quasiparticle cooling geometry through the large electron–phonon relaxation length $v^\ast=\frac{V^\ast}{LB},$1 (Baeva et al., 18 Sep 2025).

Layered superconductors introduce a qualitatively different FFI mechanism. In stacks of intrinsic Josephson junctions, a moving Josephson vortex can emit Cherenkov radiation once its velocity exceeds the minimum phase velocity of Josephson plasma waves. Numerical simulations showed a terminal velocity $v^\ast=\frac{V^\ast}{LB},$2, above which steady flux flow breaks down and the moving vortex begins generating expanding vortex–antivortex pairs. The resulting instability yields macrovortex structures, branching flux patterns, large-amplitude standing waves of magnetic flux, oscillations of the total magnetic moment, and a radiated power $v^\ast=\frac{V^\ast}{LB},$3 that increases strongly with the number of layers (Sheikhzada et al., 2019). This is still an FFI, but one driven by radiation drag and nonlinear phase dynamics rather than by LO core shrinkage or ordinary thermal runaway.

6. Applications, controversies, and open problems

FFI sets a hard operational limit for superconducting conductors, microbridges, and devices that carry large current densities in magnetic field. In Fe(Se,Te), the coexistence of high $v^\ast=\frac{V^\ast}{LB},$4, moderately high $v^\ast=\frac{V^\ast}{LB},$5, and a relatively robust $v^\ast=\frac{V^\ast}{LB},$6 led to the proposal that iron-based superconductors should be regarded as high-field superconductors with performance comparable to, or even better than, high-temperature superconductors in applications where stability under high current bias matters (Leo et al., 2017). The same study linked this to high-field magnets, fault-current limiters, and high-current transmission lines.

At the same time, FFI has become an experimental probe of nonequilibrium quasiparticle dynamics. In TiN, the agreement between FFI-derived $v^\ast=\frac{V^\ast}{LB},$7 and independently measured $v^\ast=\frac{V^\ast}{LB},$8 demonstrated that, under negligible volume pinning and controlled geometry, FFI can determine the quasiparticle energy relaxation time accurately (Baeva et al., 18 Sep 2025). The converse warning is equally clear from MoSi: if edge disorder or local pinning dominates the instability, FFI still measures something real, but not necessarily the intrinsic $v^\ast=\frac{V^\ast}{LB},$9 of the material (Budinska et al., 2021).

Single-photon detection provides a major technological motivation. Reviews of fast vortex dynamics emphasize that FFI-derived relaxation times, maximal vortex velocities, and edge-controlled local models are directly relevant for micrometer-wide strips operated near the depairing current, where the same nonequilibrium processes influence detection latency, reset, and timing performance (Dobrovolskiy, 2023). The Nb–C and MoSi results both suggest that weak volume pinning, high-quality edges, and efficient heat removal are the materials conditions most favorable for combining high $\mathbf{f}_L=\mathbf{J}\times \Phi_0,$00 with delayed instability onset (Dobrovolskiy et al., 2020).

Several controversies remain active. One concerns the boundary between intrinsic electronic instability and extrinsic thermal runaway; another concerns when global models are adequate and when local, edge-controlled descriptions are mandatory. A further open question is whether empirical indicators such as the field dependence of $\mathbf{f}_L=\mathbf{J}\times \Phi_0,$01 in iron-based superconductors can be generalized across materials or must always be interpreted through the pinning landscape, anisotropy, and cooling architecture (Leo et al., 2017). More broadly, recent reviews argue that a unified microscopic theory across clean and dirty limits, low and high temperatures, and global and local geometries is still lacking, particularly once microwave driving, magnonic coupling, or three-dimensional nanoarchitectures are introduced (Dobrovolskiy, 2023).

FFI therefore occupies a dual position in superconductivity. It is a quench mechanism that limits the usable mixed-state transport window, and it is a high-resolution probe of the nonequilibrium physics of vortices, quasiparticles, heat flow, and geometry. That duality is the reason it remains central to both fundamental vortex dynamics and superconducting-device engineering.