Peakons in Nonlinear Wave Equations

Updated 27 January 2026

Peakons are non-smooth weak solutions of nonlinear dispersive PDEs characterized by a continuous profile and a jump in the first derivative at the crest.
They emerge in integrable models such as the Camassa-Holm and Degasperis-Procesi equations, with dynamics governed by finite-dimensional Hamiltonian systems.
Their stability and interaction properties underpin modern research in nonlinear wave phenomena, offering insights into energy localization and wave breaking.

Peakons are a distinguished class of non-smooth, weak solutions for several integrable and non-integrable nonlinear dispersive partial differential equations (PDEs). Characterized by continuous profiles with a discontinuity in their first spatial derivative at the crest, peakons have deep structural connections to integrable Hamiltonian systems, nonlinear wave dynamics, and stability theory. Their prototypical example is the traveling wave solution $u(x,t) = a e^{-|x-ct|}$ for the Camassa-Holm (CH), Degasperis-Procesi (DP), and Novikov equations and their generalizations.

1. Mathematical Definition and General Properties

A peakon for a scalar PDE has the explicit form $u(x,t) = a\,e^{-|x - ct|}$ where $a > 0$ is the amplitude and $c$ is the wave speed. At $x=ct$ , $u$ is continuous but the derivative $u_x$ has a jump discontinuity. In weak-solution theory, such profiles satisfy the original nonlinear PDE under integration against suitable test functions, typically exploiting the distributional properties of derivatives and convolution operators.

For generalized equations—such as the Fornberg-Whitham ( $u_t-u_{xxt} + u_x + uu_x = u u_{xxx} + 3 u_x u_{xx}$ ), Novikov's cubic CH ( $u_t-u_{xxt} + 4u^2 u_x = 3u u_x u_{xx} + u^2 u_{xxx}$ ), and higher-order $b$ -family systems—peakons and pseudo-peakons (with smoothed higher-order derivative jumps) form natural exact weak solutions (Zhou et al., 2009, Lafortune, 2023, Zhu et al., 24 Feb 2025).

Key generic features:

Peak amplitude and speed are typically related by $a = c$ (CH/DP) or $a = \sqrt{c}$ (Novikov).
Exponential decay away from the peak with a slope parameter fixed by the equation's linear terms.
Cusp at the crest: the jump in $u_x$ at $x=ct$ is $2a$ (CH/DP/Novikov variants).
Nonzero linear momentum for individual peakons and momentum maps for multi-peakon configurations.

2. Integrable Structures and Hamiltonian Formulation

In canonical integrable models—Camassa-Holm, Degasperis-Procesi, Novikov, modified CH (mCH), and multi-component extensions—peakons emerge from Lax integrability, infinite hierarchies of conservation laws, and bi-Hamiltonian structures (Kabakouala, 2016, Anco et al., 2017, Chang et al., 2017, Xia et al., 2013, Lis, 2011).

Finite peakon trains are described by ODEs for the positions and amplitudes: $m(x,t) = u(x,t) - u_{xx}(x,t) = 2 \sum_{k=1}^N m_k(t)\,\delta(x-x_k(t))$ The evolution equations (obtained by matching distributional coefficients) are Hamiltonian with respect to non-canonical brackets, e.g.,

$\dot{x}_j = \frac{\partial H}{\partial m_j}, \quad \dot{m}_j = -\frac{\partial H}{\partial x_j}$

where $H$ is typically $H = \sum_{i,k} m_i m_k e^{-|x_i - x_k|}$ (Chang et al., 2017). In the multi-component setting (Popowicz, Geng–Xue, Novikov 2-component), the ODE structure generalizes to include cross-component peakon dynamics and additional invariants (Barnes et al., 2018, Lundmark et al., 2016, He et al., 2021).

Inverse spectral methods allow for explicit integration of multi-peakon systems via the corresponding Lax pairs and Weyl functions, resulting in soliton scattering, bound states, and breathers (Anco et al., 2017, Lundmark et al., 2016).

3. Stability Theory: Orbital and Asymptotic Results

Extensive research establishes orbital stability (Lyapunov stability modulo translation) for single and multi-peakon configurations in natural energy spaces ( $H^1$ , $L^2$ , or equation-specific weighted Sobolev norms). Key techniques include:

Modulation by the peak's position using maximum principles.
Construction of Lyapunov functionals from invariants (energy, Hamiltonians, higher degree moments).
Piecewise auxiliary functions, yielding polynomial inequalities connecting invariants and localized value at the crest.
Rigorous control on maxima, coercivity estimates, and sharp inequalities.

Notable stability results:

CH peakons, DP peakons, Novikov cubic peakons, and periodic peakons are orbitally stable in $H^1$ or analogous spaces (Kabakouala, 2016, Kabakouala, 2016, Wang et al., 2018, Chen et al., 2010, Zhou et al., 2013, Palacios, 2020).
Orbital stability for $N$ -peakon trains, provided sufficient separation and positivity of momentum density (Kabakouala, 2016, He et al., 2021).
Asymptotic stability under $H^1$ -small perturbations for Novikov peakons, with strong convergence to a shifted peakon (Palacios, 2020).
Precise error estimates for deviation in maximum, norm, and profile shape governed by conserved quantities (Zhou et al., 2013, Kabakouala, 2016).

The function space is critical: spectral and linear stability in $H^1$ ; instability may arise in weaker ( $L^2$ , $W^{1,\infty}$ ) spaces (Lafortune, 2023). Stability for higher-order and pseudo-peakons remains an open direction (Zhu et al., 24 Feb 2025).

4. Multi-Component, Higher-Order, and Non-Integrable Generalizations

Peakons are not confined to scalar integrable equations. Recent advances include:

Multi-component systems: Popowicz's coupling of CH and DP, Geng–Xue's two-field extension of Novikov, and three-component CH systems (Barnes et al., 2018, Lundmark et al., 2016, He et al., 2021, Xia et al., 2013).
Higher-order $b$ -family equations: explicit classification of $b$ -independent and $b$ -dependent peakons and pseudo-peakon families, validated up to order $J=14$ via symbolic computation (Zhu et al., 24 Feb 2025).
U(1)-invariant models from the NLS hierarchy admit conservative peakons with time-invariant $H^1$ norm and Hamiltonian Poisson geometry (Anco et al., 2017).
Non-integrable and stochastic variants: formation and persistence of peakons under transport noise in the stochastic Camassa-Holm equation, handled via variational finite element discretization and probabilistic diagnostics (Bendall et al., 2019).

Dynamical peakons—solutions with time-dependent amplitude and speed—appear in non-integrable generalizations, exhibiting behaviors such as extinction, blow-up, reversal, and breathers (Anco et al., 2019).

5. Analytical Techniques and Weak-Solution Machinery

Peakons require weak (distributional) interpretation due to their non-smoothness. Standard approaches:

Utilization of convolution with Green's functions: $(1-\partial_x^2)^{-1}f$ .
Jump conditions for derivatives, e.g., $u_x(x^+) - u_x(x^-) = 2a$ at the crest.
Averaging prescriptions for products (conservative sector), e.g., $(u^2 m)(x_j) = \tfrac12 [u^2(x_j^+) + u^2(x_j^-)] m_j$ .
Testing in the space of distributions $\mathcal D'$ for verifying weak solutions (Bendall et al., 2019, Anco et al., 2017, Chang et al., 2017).

Energy localization, monotonicity, and modulation arguments are central for controlling multi-peakon trains, especially in the presence of interactions and nonlocal effects (Kabakouala, 2016, He et al., 2021).

6. Physical and Geometric Significance

Peakons model sharply localized pulsating waves in shallow-water and related nonlinear dispersive media. Their exponential decay and crest cusp echo physical requirements for water-wave propagation and critical phenomena (wave breaking). The connection to Lie–Poisson structures in the dual of diffeomorphism groups or extensions thereof illuminates the geometric underpinnings of integrable peakon dynamics (Bendall et al., 2019, Anco et al., 2017).

Parameter dependence and bifurcation phenomena—e.g., transition from smooth solitons to peakons, emergence due to loss of well-posedness (vanishing coefficient of the highest derivative)—are especially pronounced in baby Skyrme models and higher-order nonlinear equations (Lis, 2011, Zhu et al., 24 Feb 2025).

7. Open Problems and Research Directions

Current research targets include:

Rigorous analysis of interactions between peakons and pseudo-peakons in higher-order, multi-field, and non-integrable settings (Zhu et al., 24 Feb 2025).
Stability and asymptotic theory for peakon solutions in $L^2$ , $W^{1,\infty}$ , and other non-energy function spaces (Lafortune, 2023).
Classification of geometric structures for conservative sectors in multi-component and NLS-derived systems (Anco et al., 2017).
Extension to stochastic PDEs: persistence, formation, and classification of peakons under random perturbations (Bendall et al., 2019).
Analytical proofs of symbolic and computational results for high-order peakon and pseudo-peakon families (Zhu et al., 24 Feb 2025).

The continued exploration of peakon theory thus represents a nexus of integrable systems, nonlinear wave phenomena, weak solution theory, and modern geometric analysis.

Selected Key References: