Stability of Radial Solutions

• Stability of radial solutions is the persistence of symmetry-invariant profiles under small, admissible perturbations in nonlinear PDEs and physical models.
• Spectral, variational, and energy-based methods rigorously determine stability through linearized operators, critical exponents, and Lyapunov functionals.
• These approaches underpin applications from elliptic and evolution equations to fluid dynamics and gravitational models, setting benchmarks for pattern formation.

The stability of radial solutions refers to the persistence (under small, admissible perturbations) of special, symmetry-invariant solution profiles in nonlinear evolution, elliptic, or stationary PDEs and in self-interacting physical models. Radial solutions—functions depending only on the radial variable $r=|x|$—arise in numerous equations with Euclidean or more general rotational symmetry. The study of their stability encompasses techniques from linear spectral analysis, dynamical systems, variational theory, and energy methods, and has direct implications for singularity formation, pattern robustness, and long-time dynamics in mathematical physics, PDEs, and geometric analysis.

1. Mathematical Setting and Definitions

Stability of radial solutions is typically formulated in one of the following contexts:

• Elliptic equations: A radial solution $u(r)$ to $-\Delta u = f(u)$ or generally $L[u]=0$ is linearly stable if the second variation of the associated energy

$$Q_u(v) = \int \left( |\nabla v|^2 - f'(u) v^2 \right) dx$$

is nonnegative for all $v$ in an admissible class (e.g., compact support, finite energy), and unstable otherwise. For systems, the linearized operator becomes matrix-valued, and stability is characterized by the absence of negative eigenvalues.

• Evolution equations (parabolic/hyperbolic): A steady state is (nonlinearly) stable if small perturbations produce solutions which remain close in a given norm (e.g., Lyapunov stability), or (orbitally) stable if the orbit remains close modulo symmetry (e.g., gauge, translation).
• Special structures: Radially symmetric steady states in fluid or quantum models, e.g., the barotropic Navier–Stokes system (Feireisl et al., 2024), Schrödinger–Poisson (Raynor et al., 2015), or Ginzburg–Landau vortices (Lamy et al., 2021), require adapted stability frameworks that incorporate conservation laws, anisotropy, or nonlocality.

2. Spectral, Variational, and Energy-Based Criteria

Stability is closely linked to the spectral properties of the linearized operator around the radial solution. In scalar equations,

• The principal (lowest) eigenvalue $\lambda_1$ of the linearized operator

$L_u[v] = -\Delta v - f'(u) v$

determines linear stability: $\lambda_1 \ge 0$ (stable), $\lambda_1 < 0$ (unstable).

• For systems or higher-order problems, stability is tied to the sign structure and the Fredholm index of an appropriate quadratic form or self-adjoint matrix (e.g., Ginzburg–Landau, Schrödinger–Poisson, boson stars).

Variational methods are central: many radial solutions are minimizers or critical points of an energy under symmetry or mass constraints. Stability may follow from strict minimization (e.g., ground states) or from coercivity of the second variation modulo symmetries.

Energy- or Lyapunov-function-based approaches, especially in evolutionary problems, allow one to conclude global or unconditional stability: if a properly constructed Lyapunov functional is monotone and vanishes only at the steady state, global convergence follows (Feireisl et al., 2024).

3. Stability Thresholds, Exponents, and Structural Dichotomies

A central finding across models is the existence of critical exponents or indices demarcating changes in the stability of radial solutions:

• Semilinear equations in $\mathbb R^N$: There are large and small stable solutions, separated by explicit dimension-dependent exponents (Villegas, 2014, Miyamoto et al., 18 Jan 2026). In the elliptic Lane–Emden case, the Joseph–Lundgren exponent $p_{JL}$ is critical:
  • For $p < p_{JL}$ (supercritical but sub-Joseph–Lundgren), positive radial solutions are unstable.
  • For $p \ge p_{JL}$, stable (including singular) radial solutions exist provided spectral/structural quantities (involving $f'(u) F(u)$, with $F(u)=\int^\infty_u \frac{dt}{f(t)}$) are in a specified range (Miyamoto et al., 18 Jan 2026).
  • Similar phenomena arise for biharmonic and polyharmonic equations, but with more elaborate spectral root structures and dependence on dimension and order (Karageorgis, 2011, Farina et al., 2014).
• Mode decompositions: For equations with angular Fourier modes, such as Ginzburg–Landau vortices, the classical monotonicity $Q_n \ge Q_1$ of the stability quadratic form breaks down with anisotropy, leading to “mode-dependent” stability loss; higher modes may destabilize even if lower modes are stable (Lamy et al., 2021).
• Nonlocal/nonlinear systems: For (e.g.) repulsive-attractive nonlocal interaction equations, the sign of kernel derivatives at the diagonal (fattening and shift conditions) gives necessary and sufficient conditions for strict exponential stability (or instability) toward radial shell states (Balague et al., 2011).

4. Asymptotic and Quantitative Behavior of Radial Solutions

The asymptotic profile and “quantitative closeness” to radial symmetry are frequently quantified:

• Sharp dichotomies: Any nonconstant stable radial solution outside a compact set either grows or decays at rates $r^{-N/2\pm \sqrt{N-1}+2}$ (or logarithmic rates at critical dimensions) in $\mathbb R^N$ (Villegas, 2014). These rates are proven sharp by optimal Hardy-type examples.
• ODE–autonomous methods: Emden–Fowler transforms convert the radial elliptic PDE into an autonomous ODE or system, enabling spectral analysis of asymptotic decay, e.g., in parabolic Hénon–Lane–Emden systems (Devine et al., 2024). In these cases, regular solutions converge monotonically (as $r\to\infty$) to explicit singular power-law profiles, with subleading corrections governed by roots of a characteristic polynomial (with logarithmic terms at threshold).
• Quantitative symmetry breaking: For semilinear equations with small non-constant coefficients or in non-energy settings, quantitative “deficit” measures control proximity to perfect radial symmetry in strong norms; explicit logarithmic or algebraic rates of “almost radiality” are proved depending on criticality and regularity (Ciraolo et al., 20 Jan 2025).

5. Special Models and Systematic Results

Theoretical developments apply across diverse equations and physical models:

• Anisotropic Ginzburg–Landau: The degree-1 radial vortex is stable for $\delta$ in a sharp subinterval $(\delta_1,0]$ of the anisotropy parameter; outside this range, either lower or higher angular modes render the solution unstable. The breakdown of monotone mode stability is a key novelty (Lamy et al., 2021).
• Hénon–Lane–Emden systems: For the parabolic system, if $(p,q)$ is on or above the Joseph–Lundgren critical curve (defined via an explicit quartic polynomial in exponents and dimension), then every positive radial steady state is Lyapunov stable in weighted norms. Logarithmic corrections to the asymptotics and norm structure are needed at the threshold (Devine et al., 2024).
• Navier–Stokes and fluid models: Radially symmetric, steady compressible flows in domains (with in/outflow boundary conditions) can be unconditionally stable: all solutions (regardless of symmetry) converge to the unique radial stationary state, via the relative-energy Lyapunov method (Feireisl et al., 2024).
• Radial oscillations in stellar and black hole models: For self-gravitating systems (boson stars, anisotropic stars, Skyrmions), the spectrum of radial oscillations and the mass–amplitude curve detect stability: stability holds up to the point of maximal mass, where an explicit mode crosses zero (via the Sturm–Liouville or Hamilton–Jacobi methods) (Alcubierre et al., 2021, Horvat et al., 2010, Doneva et al., 2011).

6. Techniques and Analytical Approaches

A range of robust mathematical tools support the stability analysis:

• Spectral and index theory: Sturm–Liouville, Floquet, and Hamiltonian index methods yield sharp instability/stability criteria.
• Hardy–Rellich and energy inequalities: These underpin optimal pointwise and integral bounds, crucial for “outside a compact set” stability and for polyharmonic/exponential equations (Farina et al., 2014, Karageorgis, 2011).
• Moving planes and monotonicity formulas: Key in classifying uniqueness and structure of positive radial solutions, and in quantifying radial symmetry (Ciraolo et al., 20 Jan 2025, Ciraolo et al., 2013).
• Dynamical systems: Autonomous ODE reductions (Emden–Fowler, shooting methods) underpin the analysis of both steady state profiles and their attractivity.
• Lyapunov functionals and mass transport: Global convergence results leverage the construction of monotonic functionals (relative energy, Wasserstein distance dissipation in nonlocal models) (Feireisl et al., 2024, Balague et al., 2011).

7. Physical and Mathematical Significance

Stability of radial solutions is central in understanding the emergence, persistence, and breakdown of symmetry-driven patterns in nonlinear PDEs and mathematical physics. The phenomena uncovered—mode-dependent instabilities, sharp existence and nonexistence thresholds, connection to spectral curves (e.g., Joseph–Lundgren), and the corresponding attractor structure—establish radial solutions as both mathematically tractable and physically meaningful benchmarks in nonlinear analysis. Advances in classification, quantitative stability rates, and their extension to systems and non-Euclidean geometries reflect a broad unification of ideas across PDE theory, dynamical systems, and geometric analysis (Miyamoto et al., 18 Jan 2026, Berchio et al., 2022, Berchio et al., 2012).