# Peakon

In the theory of integrable systems, a peakon ("peaked soliton") is a soliton with discontinuous first derivative; the wave profile is shaped like the graph of the function $e^{-|x|}$. Some examples of non-linear partial differential equations with (multi-)peakon solutions are the Camassa–Holm shallow water wave equation, the Degasperis–Procesi equation and the Fornberg–Whitham equation. Since peakon solutions are only piecewise differentiable, they must be interpreted in a suitable weak sense. The concept was introduced in 1993 by Camassa and Holm in the short but much cited paper where they derived their shallow water equation.

## A family of equations with peakon solutions

The primary example of a PDE which supports peakon solutions is
    $u_{t}-u_{xxt}+(b+1)uu_{x}=bu_{x}u_{xx}+uu_{xxx},\,$
where $u(x,t)$ is the unknown function, and b is a parameter. In terms of the auxiliary function $m(x,t)$ defined by the relation $m=u-u_{xx}$, the equation takes the simpler form
    $m_{t}+m_{x}u+bmu_{x}=0.\,$
This equation is integrable for exactly two values of b, namely b = 2 (the Camassa–Holm equation) and b = 3 (the Degasperis–Procesi equation).

## Single peakon solution

The PDE above admits the travelling wave solution $u(x,t)=c\,e^{-|x-ct|}$, which is a peaked solitary wave with amplitude c and speed c. This solution is called a (single) peakon solution, or simply a peakon. If c is negative, the wave moves to the left with the peak pointing downwards, and then it is sometimes called an antipeakon.
It is not immediately obvious in what sense the peakon solution satisfies the PDE. Since the derivative u$_x$ has a jump discontinuity at the peak, the second derivative u$_{xx}$ must be taken in the sense of distributions and will contain a Dirac delta function; in fact, $m=u-u_{xx}=c\,\delta (x-ct)$. Now the product $mu_{x}$ occurring in the PDE seems to be undefined, since the distribution m is supported at the very point where the derivative u$_x$ is undefined. An ad hoc interpretation is to take the value of u$_x$ at that point to equal the average of its left and right limits (zero, in this case). A more satisfactory way to make sense of the solution is to invert the relationship between u and m by writing $m=(G/2)*u$, where $G(x)=\exp(-|x|)$, and use this to rewrite the PDE as a (nonlocal) hyperbolic conservation law:
    $\partial _{t}u+\partial _{x}\left[{\frac {u^{2}}{2}}+{\frac {G}{2}}*\left({\frac {bu^{2}}{2}}+{\frac {(3-b)u_{x}^{2}}{2}}\right)\right]=0.$
(The star denotes convolution with respect to x.) In this formulation the function u can simply be interpreted as a weak solution in the usual sense.

## Multipeakon solutions

Multipeakon solutions are formed by taking a linear combination of several peakons, each with its own time-dependent amplitude and position. (This is a very simple structure compared to the multisoliton solutions of most other integrable PDEs, like the Korteweg–de Vries equation for instance.) The n-peakon solution thus takes the form
    $u(x,t)=\sum _{i=1}^{n}m_{i}(t)\,e^{-|x-x_{i}(t)|},$
where the 2n functions $x_{i}(t)$ and $m_{i}(t)$ must be chosen suitably in order for u to satisfy the PDE. For the "b-family" above it turns out that this ansatz indeed gives a solution, provided that the system of ODEs
    ${\dot {x}}_{k}=\sum _{i=1}^{n}m_{i}e^{-|x_{k}-x_{i}|},\qquad {\dot {m}}_{k}=(b-1)\sum _{i=1}^{n}m_{k}m_{i}\operatorname {sgn}(x_{k}-x_{i})e^{-|x_{k}-x_{i}|}\qquad (k=1,\dots ,n)$
is satisfied. (Here sgn denotes the sign function.) Note that the right-hand side of the equation for $x_{k}$ is obtained by substituting $x=x_{k}$ in the formula for u. Similarly, the equation for $m_{k}$ can be expressed in terms of $u_{x}$, if one interprets the derivative of $\exp(-|x|)$ at x = 0 as being zero. This gives the following convenient shorthand notation for the system:
    ${\dot {x}}_{k}=u(x_{k}),\qquad {\dot {m}}_{k}=-(b-1)m_{k}u_{x}(x_{k})\qquad (k=1,\dots ,n).$
The first equation provides some useful intuition about peakon dynamics: the velocity of each peakon equals the elevation of the wave at that point.

## Explicit solution formulas

In the integrable cases b = 2 and b = 3, the system of ODEs describing the peakon dynamics can be solved explicitly for arbitrary n in terms of elementary functions, using inverse spectral techniques. For example, the solution for n = 3 in the Camassa–Holm case b = 2 is given by
    ${\begin{aligned}x_{1}(t)&=\log {\frac {(\lambda _{1}-\lambda _{2})^{2}(\lambda _{1}-\lambda _{3})^{2}(\lambda _{2}-\lambda _{3})^{2}a_{1}a_{2}a_{3}}{\sum _{j<k}\lambda _{j}^{2}\lambda _{k}^{2}(\lambda _{j}-\lambda _{k})^{2}a_{j}a_{k}}}\\x_{2}(t)&=\log {\frac {\sum _{j<k}(\lambda _{j}-\lambda _{k})^{2}a_{j}a_{k}}{\lambda _{1}^{2}a_{1}+\lambda _{2}^{2}a_{2}+\lambda _{3}^{2}a_{3}}}\\x_{3}(t)&=\log(a_{1}+a_{2}+a_{3})\\m_{1}(t)&={\frac {\sum _{j<k}\lambda _{j}^{2}\lambda _{k}^{2}(\lambda _{j}-\lambda _{k})^{2}a_{j}a_{k}}{\lambda _{1}\lambda _{2}\lambda _{3}\sum _{j<k}\lambda _{j}\lambda _{k}(\lambda _{j}-\lambda _{k})^{2}a_{j}a_{k}}}\\m_{2}(t)&={\frac {\left(\lambda _{1}^{2}a_{1}+\lambda _{2}^{2}a_{2}+\lambda _{3}^{2}a_{3}\right)\sum _{j<k}(\lambda _{j}-\lambda _{k})^{2}a_{j}a_{k}}{\left(\lambda _{1}a_{1}+\lambda _{2}a_{2}+\lambda _{3}a_{3}\right)\sum _{j<k}\lambda _{j}\lambda _{k}(\lambda _{j}-\lambda _{k})^{2}a_{j}a_{k}}}\\m_{3}(t)&={\frac {a_{1}+a_{2}+a_{3}}{\lambda _{1}a_{1}+\lambda _{2}a_{2}+\lambda _{3}a_{3}}}\end{aligned}}$
where $a_{k}(t)=a_{k}(0)e^{t/\lambda _{k}}$, and where the 2n constants $a_{k}(0)$ and $\lambda _{k}$ are determined from initial conditions. The general solution for arbitrary n can be expressed in terms of symmetric functions of $a_{k}$ and $\lambda _{k}$. The general n-peakon solution in the Degasperis–Procesi case b = 3 is similar in flavour, although the detailed structure is more complicated.