Shrinkers, expanders, and the unique continuation beyond generic blowup in the heat flow for harmonic maps between spheres

Abstract: Using mixed analytical and numerical methods we investigate the development of singularities in the heat flow for corotational harmonic maps from the $d$-dimensional sphere to itself for $3\leq d\leq 6$. By gluing together shrinking and expanding asymptotically self-similar solutions we construct global weak solutions which are smooth everywhere except for a sequence of times $T_1<T_2<...<T_k<\infty$ at which there occurs the type I blow-up at one of the poles of the sphere. We show that in the generic case the continuation beyond blow-up is unique, the topological degree of the map changes by one at each blow-up time $T_i$, and eventually the solution comes to rest at the zero energy constant map.

• The paper demonstrates that the unique, linearly stable shrinker f1 acts as the attractor for generic finite-time blow-up.
• It employs rigorous adaptive numerical simulations and matched asymptotic expansions to validate the gluing of shrinkers to expanders.
• The study reveals that the continuation beyond blow-up is achieved through a unique topological degree change, decrementing by one per singularity.

Analysis of Shrinkers, Expanders, and Unique Continuation in the Heat Flow for Harmonic Maps Between Spheres

Problem Context and Statement

The paper investigates the global behavior of the heat flow for corotational harmonic maps from $S^d$ to $S^d$ for $3 \leq d \leq 6$, focusing on singularity formation (finite-time blow-up) and the question of unique continuation of weak solutions beyond blow-up. The work explores the landscape of self-similar solutions—shrinking (“shrinkers”) and expanding (“expanders”) profiles—and demonstrates, by combined analytical and rigorous adaptive numerical methods, that gluing particular shrinkers to suitable expanders enables a unique continuation for generic data. The results have implications for the theory of geometric flows and the study of parabolic PDEs with singularities.

Self-Similar Blow-Up: Shrinkers

The analysis centers on the reduction, via symmetry and scaling arguments, of the harmonic map heat flow to a radial parabolic equation, with forms specialized to corotational symmetry. In the vicinity of a singularity at the pole, the equation reduces to one on $\mathbb{R}^d$ with the blow-up profile described in self-similar variables.

The self-similar (shrinking) solutions, called shrinkers, satisfy a nonlinear second-order ODE. It is established (building on prior shooting arguments by Fan) that for $3 \leq d \leq 6$ there are countably infinite regular shrinkers $f_n(y)$, indexed by their oscillation count. Their analytic structure is characterized via matched asymptotic expansions and the shooting method, supplemented by quantitative scaling laws for parameter dependence.

The crucial point is the linear stability of these shrinkers: only the first (ground-state) shrinker $f_1$ is linearly stable, making it the attractor for the generic blow-up profile. This is confirmed by careful eigenvalue analysis of the associated linearized operator, interpretable via Sturm-Liouville theory and self-adjointness in appropriate weighted Hilbert spaces.

Expanders and Gluing Construction

To achieve a weak solution flowing past the blow-up instant, the paper constructs solutions based on gluing (in self-similar variables) a shrinker to an expander, the latter describing post-singularity evolution. The expander family is parametrized by the initial derivative at the origin.

Unlike shrinkers, expanders form a continuous one-parameter family, regular for any value of the initial slope; their asymptotic behavior is analyzed in detail using both linear and nonlinear techniques. Matching the far-field value of an expander to the singular value attained by the shrinker at blow-up leads to an algebraic equation for the gluing parameter.

A key numerical and analytic finding is that, for the stable shrinker ($n=1$), there exists only a single expander which can be glued to it to continue the flow in a manner compatible with the evolution. Moreover, this continuation necessarily involves a jump in topological degree, as the map value at the origin jumps from $0$ to $\pi$ at the blow-up time, leading to degree decrement by one per singularity.

For higher-index shrinkers, multiple (but unstable) continuation branches exist, reflecting the codimension of nongeneric blow-ups.

Numerical Validation

The paper employs advanced adaptive numerical simulations (moving mesh and Sundman time transformation) to validate the analytic predictions. The results unambiguously show:

• As a generic singularity is approached, the profile locally converges to the stable shrinker.
• Post-blow-up, after the unavoidable loss of resolution across the singularity, the local profile converges to the unique expander dictated by the gluing condition.
• The change in topological degree at each blow-up is clearly visible in simulations involving initial maps of higher degree, with the eventual relaxation to the trivial (zero energy, degree zero) solution.

The numerics are shown to be consistent with the eigenvalue analysis of the linearized problem, affirming the theoretical claim of uniqueness for the generic scenario.

Theoretical and Practical Implications

Theoretical Implications:

• The demonstration of unique continuation past generic (type I) singularities for heat flows of harmonic maps between spheres is in sharp contrast to past results for other geometric/parabolic flows (including mean curvature flow and Yang-Mills flow), where nonuniqueness is generic.
• The explicit connection between the topology (degree change) and singularity formation is clarified: every generic blow-up reduces the degree by one, with the flow eventually relaxing only after eliminating all homotopic (topological) obstructions.
• The mechanism is shown to hinge on the existence of a unique, stable self-similar shrinker and a unique matching expander, a structure not present in higher dimensions.

Practical Implications:

• The work provides a concrete algorithm for regularizing the heat flow numerically and conceptually, by gluing together shrinkers and expanders and tracking degree changes.
• The results generalize to other equivariant geometric flows and highlight how the geometrical/topological constraints manifest in PDE singularity dynamics.
• The analytic techniques—particularly the matched asymptotics, spectral analysis, and adaptive numerical schemes—are broadly applicable to a range of supercritical and critical parabolic PDEs.

Directions for Further Research

• Extension to $\ell$-equivariant harmonic map flows and Yang-Mills heat flows is indicated, where analogous gluing constructions and stability analysis must consider modified symmetry and scalar curvature roles.
• For $d \geq 7$, the scenario changes dramatically as shrinkers disappear; type II blow-up emerges, with as yet unclear continuation and degree change properties.
• Investigation into analogous phenomena in semilinear heat equations and mean curvature flows remains open, particularly for the role of unstable self-similar solutions and uniqueness of continuation.
• Rigorous proof of the conjectured (“Conjecture 1”) uniqueness of continuation remains an outstanding analytic challenge, to be settled beyond the current numerical-analytic evidence.

Conclusion

The paper provides a complete qualitative and quantitative analysis of singularity formation and unique continuation for the heat flow of corotational harmonic maps between spheres for dimensions $3 \leq d \leq 6$. By identifying and leveraging the structure and stability of self-similar shrinkers and expanders, it establishes a scenario where, for generic data, the evolution can be uniquely continued beyond blow-up, with a well-defined topological mechanism for singularity resolution. This offers valuable insights into the interplay between topology, analysis, and numerics in geometric evolution equations.