Interaction of two peakons — which are sharp-crested soliton solutions to the Camassa–Holm equation. The wave profile (solid curve) is formed by the simple linear addition of two peakons (dashed curves):
The evolution of the individual peakon positions and , as well as the evolution of the peakon amplitudes and is however less trivial: this is determined in a non-linear fashion by the interaction.
The Camassa–Holm equation can be written as the system of equations:
with p the (dimensionless) pressure or surface elevation. This shows that the Camassa–Holm equation is a model for shallow water waves with non-hydrostatic pressure and a water layer on a horizontal bed.
The linear dispersion characteristics of the Camassa–Holm equation are:
with ω the angular frequency and k the wavenumber. Not surprisingly, this is of similar form as the one for the Korteweg–de Vries equation, provided κ is non-zero. For κ equal to zero, the Camassa–Holm equation has no frequency dispersion — moreover, the linear phase speed is zero for this case. As a result, κ is the phase speed for the long-wave limit of k approaching zero, and the Camassa–Holm equation is (if κ is non-zero) a model for one-directional wave propagation like the Korteweg–de Vries equation.
The Camassa–Holm equation is an integrable system. Integrability means that there is a change of variables (action-angle variables) such that the evolution equation in the new variables is equivalent to a linear flow at constant speed. This change of variables is achieved by studying an associated isospectral/scattering problem, and is reminiscent of the fact that integrable classical Hamiltonian systems are equivalent to linear flows at constant speed on tori. The Camassa–Holm equation is integrable provided that the momentum
is positive — see  and  for a detailed description of the spectrum associated to the isospectral problem, for the inverse spectral problem in the case of spatially periodic smooth solutions, and  for the inverse scattering approach in the case of smooth solutions that decay at infinity.
representing waves of permanent shape f that propagate at constant speed c. These waves are called solitary waves if they are localized disturbances, that is, if the wave profile f decays at infinity. If the solitary waves retain their shape and speed after interacting with other waves of the same type, we say that the solitary waves are solitons. There is a close connection between integrability and solitons.  In the limiting case when κ = 0 the solitons become peaked (shaped like the graph of the function f(x) = e-|x|), and they are then called peakons. It is possible to provide explicit formulas for the peakon interactions, visualizing thus the fact that they are solitons. For the smooth solitons the soliton interactions are less elegant. This is due in part to the fact that, unlike the peakons, the smooth solitons are relatively easy to describe qualitatively — they are smooth, decaying exponentially fast at infinity, symmetric with respect to the crest, and with two inflection points — but explicit formulas are not available. Notice also that the solitary waves are orbitally stable i.e. their shape is stable under small perturbations, both for the smooth solitons and for the peakons.
The Camassa–Holm equation models breaking waves: a smooth initial profile with sufficient decay at infinity develops into either a wave that exists for all times or into a breaking wave (wave breaking being characterized by the fact that the solution remains bounded but its slope becomes unbounded in finite time). The fact that the equations admits solutions of this type was discovered by Camassa and Holm and these considerations were subsequently put on a firm mathematical basis. It is known that the only way singularities can occur in solutions is in the form of breaking waves. Moreover, from the knowledge of a smooth initial profile it is possible to predict (via a necessary and sufficient condition) whether wave breaking occurs or not. As for the continuation of solutions after wave breaking, two scenarios are possible: the conservative case and the dissipative case (with the first characterized by conservation of the energy, while the dissipative scenario accounts for loss of energy due to breaking).
It can be shown that for sufficiently fast decaying smooth initial conditions with positive momentum splits into a finite number and solitons plus a decaying dispersive part. More precisely, one can show the following for : Abbreviate . In the soliton region the solutions splits into a finite linear combination solitons. In the region the solution is asymptotically given by a modulated sine function whose amplitude decays like . In the region the solution is asymptotically given by a sum of two modulated sine function as in the previous case. In the region the solution decays rapidly. In the case the solution splits into an infinite linear combination of peakons (as previously conjectured).
In 2003, a unified wave model (UWM) for progressive gravity waves in finite water depth was proposed by Liao. Based on the symmetry and the exact wave equations, the UWM admits not only all traditional smooth periodic/solitary waves but also the peaked solitary waves including the famous peaked solitary waves of Camassa–Holm equation mentioned above. Thus, the UWM unifies the smooth and peaked waves in finite water depth. In other words, the peaked solitary waves are consistent with the traditional, smooth ones, and thus are as acceptable as the smooth ones.
It is found that the peaked solitary waves in finite water depth have some unusual characteristics. First of all, it has a peaked wave elevation with a discontinuous vertical velocity at crest. Secondly, unlike the smooth waves whose horizontal velocity decays exponentially from free surface to the bottom, the horizontal velocity of the peaked solitary waves always increases from free surface to the bottom. Especially, different from the smooth waves whose phase speed is dependent upon wave height, the phase speed of the peaked solitary waves in finite water depth have nothing to do with the wave height! In other words, the peaked solitary waves in finite water depth are non-dispersive.
The above usual characteristics of the peaked solitary waves in finite water depth are quite different from those of the traditional, smooth waves, and thus might challenge some traditional viewpoints. Even so, they could enrich and deepen our understandings about the peaked solitary waves and the Camassa–Holm equation.
Alber, Mark S.; Camassa, Roberto; Holm, Darryl D.; Marsden, Jerrold E. (1994), "The geometry of peaked solitons and billiard solutions of a class of integrable PDEs", Lett. Math. Phys.32 (2): 137–151, Bibcode:1994LMaPh..32..137A, doi:10.1007/BF00739423
Alber, Mark S.; Camassa, Roberto; Holm, Darryl D.; Fedorov, Yuri N.; Marsden, Jerrold E. (2001), "The complex geometry of weak piecewise smooth solutions of integrable nonlinear PDE's of shallow water and Dym type", Comm. Math. Phys.221 (1): 197–227, arXiv:nlin/0105025, Bibcode:2001CMaPh.221..197A, doi:10.1007/PL00005573
Artebrant, Robert; Schroll, Hans Joachim (2006), "Numerical simulation of Camassa–Holm peakons by adaptive upwinding", Applied Numerical Mathematics56 (5): 695–711, doi:10.1016/j.apnum.2005.06.002
Beals, Richard; Sattinger, David H.; Szmigielski, Jacek (2005), "Periodic peakons and Calogero–Françoise flows", J. Inst. Math. Jussieu4 (1): 1–27, doi:10.1017/S1474748005000010
Boutet de Monvel, Anne; Shepelsky, Dmitry (2005), "The Camassa–Holm equation on the half-line", C. R. Math. Acad. Sci. Paris341 (10): 611–616, doi:10.1016/j.crma.2005.09.035
Boutet de Monvel, Anne; Shepelsky, Dmitry (2006), "Riemann–Hilbert approach for the Camassa–Holm equation on the line", C. R. Math. Acad. Sci. Paris343 (10): 627–632, doi:10.1016/j.crma.2006.10.014
Byers, Peter (2006), "Existence time for the Camassa–Holm equation and the critical Sobolev index", Indiana Univ. Math. J.55 (3): 941–954, doi:10.1512/iumj.2006.55.2710
Camassa, Roberto (2003), "Characteristics and the initial value problem of a completely integrable shallow water equation", Discrete Contin. Dyn. Syst. Ser. B3 (1): 115–139, doi:10.3934/dcdsb.2003.3.115
Golovko, V.; Kersten, P.; Krasilʹshchik, I.; Verbovetsky, A. (2008), "On integrability of the Camassa–Holm equation and its invariants: a geometrical approach", Acta Appl. Math.101 (1–3): 59–83, arXiv:0812.4681, doi:10.1007/s10440-008-9200-z
Himonas, A. Alexandrou; Misiołek, Gerard (2001), "The Cauchy problem for an integrable shallow-water equation", Differential and Integral Equations14 (7): 821–831
Himonas, A. Alexandrou; Misiołek, Gerard (2005), "High-frequency smooth solutions and well-posedness of the Camassa–Holm equation", Int. Math. Res. Not.2005 (51): 3135–3151, doi:10.1155/IMRN.2005.3135
Lenells, Jonatan (2007), "Infinite propagation speed of the Camassa–Holm equation", J. Math. Anal. Appl.325 (2): 1468–1478, doi:10.1016/j.jmaa.2006.02.045
Li, Luen-Chau (2008), "Factorization problem on the Hilbert–Schmidt group and the Camassa–Holm equation", Comm. Pure Appl. Math.61 (2): 186–209, doi:10.1002/cpa.20207
Lombardo, Maria Carmela; Sammartino, Marco; Sciacca, Vincenzo (2005), "A note on the analytic solutions of the Camassa–Holm equation", C. R. Math. Acad. Sci. Paris341 (11): 659–664, doi:10.1016/j.crma.2005.10.006
Loubet, Enrique (2006), "Genesis of solitons arising from individual flows of the Camassa–Holm hierarchy", Comm. Pure Appl. Math.59 (3): 408–465, doi:10.1002/cpa.20109
Misiołek, Gerard (2005), "Classical solutions of the periodic Camassa–Holm equation", Geom. Funct. Anal.12 (5): 1080–1104, doi:10.1007/PL00012648
Xin, Zhouping; Zhang, Ping (2002), "On the uniqueness and large time behavior of the weak solutions to a shallow water equation", Comm. Partial Differential Equations27 (9–10): 1815–1844, doi:10.1081/PDE-120016129
Zampogni, Luca (2007), "On algebro-geometric solutions of the Camassa–Holm hierarchy", Adv. Nonlinear Stud.7 (3): 345–380