Instability of anchored spirals in geometric flows

We investigate existence, stability, and instability of anchored rotating spiral waves in a model for geometric curve evolution. We find existence in a parameter regime limiting on a purely eikonal curve evolution. We study stability and instability theoretically, in the aforementioned limiting regime, and numerically. We find convective and absolute oscillatory instability, as well as saddle-node bifurcations. Our results in particular shed light on the instability of spiral waves in reaction-diffusion systems caused by an instability of wave trains against transverse modulations.

💡 Research Summary

The paper investigates the existence, stability, and instability of anchored rotating spiral waves within a geometric curve‑evolution framework. The authors consider a planar curve γ(t) whose normal velocity is prescribed by a linear combination of curvature κ and its second arclength derivative κ_{ss}:

 c ∂ₜγ = D₂ κ + D₄ κ_{ss},

where V (the base normal speed) is hidden in the constant c, while D₂ and D₄ are coefficients representing line‑tension (or surface tension) and a higher‑order regularizing term, respectively. The curve is anchored at an inner circle of radius R_i and may be bounded externally by an outer circle of radius R_o. Boundary conditions fix contact angles at the inner and outer boundaries and impose zero curvature there.

Main analytical strategy
The authors focus on the regime where the core radius R_i is large (R_i ≫ 1) and the coefficients D₂, D₄ are small, scaling as D₂ = O(ε) and D₄ = O(ε³) for a small parameter ε. In the singular limit ε→0 (i.e., D₂ = D₄ = 0) the evolution reduces to the eikonal equation, and an explicit rigidly rotating Archimedean spiral solution exists. Its angular frequency is given by

 ω = V sin θ_i / R_i,

where θ_i is the prescribed contact angle at the inner boundary. This solution serves as the base state for a singular‑perturbation analysis.

To treat the small curvature terms, the authors recast the fourth‑order ordinary differential equation for the spiral profile into a first‑order system by introducing variables for φ_r, φ_{rr}, and φ_{rrr}. They then apply Fenichel’s geometric singular perturbation theory. A key step is the compactification of the radial coordinate via α = 1/r, which brings the point at infinity (r→∞) to α = 0. After a time‑rescaling τ = ε⁻³ r, the system becomes regular at α = 0, revealing a slow manifold of equilibria that corresponds to the eikonal spiral. The dynamics on and near this manifold capture the influence of D₂ and D₄.

Stability results
Linearizing about the spiral solution yields a non‑self‑adjoint operator whose spectrum consists of an “essential” part (associated with the far‑field planar wave trains) and a discrete set of eigenvalues stemming from the core region. The authors prove:

1. Positive D₂ (line‑tension): All discrete eigenvalues have negative real parts; the spiral is linearly stable. The regularizing D₄ term only damps short‑wavelength perturbations.

2. Negative D₂: As D₂ decreases past a critical value D₂^{crit}(R_i, D₄, V), a pair of complex conjugate eigenvalues crosses the imaginary axis, producing a Hopf bifurcation. The associated eigenfunction grows like e^{α r} with α > 0, i.e., it is exponentially localized near the core but becomes super‑exponential as r→∞. This signals both convective instability (perturbations are advected outward) and absolute instability (growth at a fixed spatial point). The Hopf bifurcation gives rise to oscillatory spiral cores.

3. Saddle‑node bifurcations: Varying parameters can also cause the disappearance or creation of spiral solutions via saddle‑node points on the slow manifold. These are detected numerically and are consistent with the singular‑perturbation picture.

Numerical verification
Two complementary numerical approaches are employed:

• Direct time integration of the full fourth‑order PDE (2.2) with boundary conditions (2.3). Simulations starting from a straight anchored curve evolve into outward‑rotating Archimedean spirals, confirming the predicted frequency ω = V sin θ_i / R_i. When D₂ < 0, the simulations display growing transverse modulations of the far‑field wave trains, eventually destroying the spiral.

• Continuation (AUTO) calculations of the boundary‑value problem for the spiral profile and its linear eigenvalue problem. By tracking eigenvalues as D₂ varies, the Hopf and saddle‑node points are located precisely, matching the analytical critical values derived from the reduced slow‑fast system.

The numerical results also illustrate that the “essential spectrum” at the far field lies on the imaginary axis (as in the eikonal limit), but a cluster of discrete eigenvalues approaches it from the left and crosses into the right half‑plane when D₂ becomes sufficiently negative. This mirrors the phenomenon observed in reaction‑diffusion systems where transverse instabilities of planar wave trains destabilize spiral waves.

Implications for reaction‑diffusion media
The authors argue that the geometric model captures the core mechanism behind the transverse instability of spiral waves in excitable or oscillatory reaction‑diffusion media. In those systems, the spiral’s far field is essentially a train of planar waves; if the wave train is unstable to transverse perturbations, the spiral inherits this instability. The present analysis provides a clean, analytically tractable setting where the same mechanism appears as a Hopf bifurcation driven by a negative line‑tension term (D₂ < 0) and regularized by a higher‑order diffusion term (D₄ > 0).

Conclusions
The paper establishes, via rigorous singular‑perturbation theory and extensive numerics, that anchored rotating spirals exist in a broad parameter regime of the geometric curve‑evolution model. Their stability hinges on the sign and magnitude of the curvature‑dependent line‑tension D₂. Positive D₂ yields stable spirals, while sufficiently negative D₂ triggers a Hopf bifurcation leading to convective and absolute oscillatory instabilities, as well as saddle‑node annihilation of solutions. The work bridges geometric flow theory with the well‑studied instability of spiral waves in reaction‑diffusion systems, offering a new perspective on how transverse modulations of wave trains can destabilize rotating patterns.