Peakon

The primary example of a PDE which supports peakon solutions is

${\displaystyle u_{t}-u_{xxt}+(b+1)uu_{x}=bu_{x}u_{xx}+uu_{xxx},\,}$

where ${\displaystyle u(x,t)}$ is the unknown function, and b is a parameter.[2] In terms of the auxiliary function ${\displaystyle m(x,t)}$ defined by the relation ${\displaystyle m=u-u_{xx}}$, the equation takes the simpler form

${\displaystyle 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).

The PDE above admits the travelling wave solution ${\displaystyle 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, ${\displaystyle m=u-u_{xx}=c\,\delta (x-ct)}$. Now the product ${\displaystyle 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 ${\displaystyle m=(G/2)*u}$, where ${\displaystyle G(x)=\exp(-|x|)}$, and use this to rewrite the PDE as a (nonlocal) hyperbolic conservation law:

${\displaystyle \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.[3]

Two-peakon wave profile (solid curve) formed by adding two peakons (dashed curves): ${\displaystyle u=m_{1}\,e^{-|x-x_{1}|}+m_{2}\,e^{-|x-x_{2}|}}$

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

${\displaystyle u(x,t)=\sum _{i=1}^{n}m_{i}(t)\,e^{-|x-x_{i}(t)|},}$

where the 2n functions ${\displaystyle x_{i}(t)}$ and ${\displaystyle 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

${\displaystyle {\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 ${\displaystyle x_{k}}$ is obtained by substituting ${\displaystyle x=x_{k}}$ in the formula for u. Similarly, the equation for ${\displaystyle m_{k}}$ can be expressed in terms of ${\displaystyle u_{x}}$, if one interprets the derivative of ${\displaystyle \exp(-|x|)}$ at x = 0 as being zero. This gives the following convenient shorthand notation for the system:

${\displaystyle {\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.

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[4]

${\displaystyle {\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 ${\displaystyle a_{k}(t)=a_{k}(0)e^{t/\lambda _{k}}}$, and where the 2n constants ${\displaystyle a_{k}(0)}$ and ${\displaystyle \lambda _{k}}$ are determined from initial conditions. The general solution for arbitrary n can be expressed in terms of symmetric functions of ${\displaystyle a_{k}}$ and ${\displaystyle \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.[5]