--- title: 'Novikov Equations: Integrability & Peakons' url: https://www.emergentmind.com/topics/novikov-equations type: topic --- # Novikov Equations: Integrability & Peakons “Novikov equations” denotes a small but non-uniform family of mathematical objects centered on the scalar Novikov equation \[ u_t-u_{txx}+4u^2u_x-3uu_xu_{xx}-u^2u_{xxx}=0, \] equivalently \[ m_t+u^2m_x+3uu_xm=0,\qquad m=u-u_{xx}, \] an integrable Camassa–Holm type equation with cubic nonlinearities, Lax pair representations, and peakon solutions. In current usage, the term also covers weakly dissipative, multicomponent, \(b\)-family, and \(\lambda\)-family variants, the associated Novikov equation, and several distinct Novikov-named equations that are not equivalent to the \(1+1\)-dimensional cubic shallow-water model, notably the Novikov–Veselov equation and the algebraic-geometric “Novikov equations” arising from commuting differential operators [1009.1820][1207.0968][1105.3903][2507.14978]. ## 1. Terminological scope and canonical scalar form In the PDE literature, the canonical Novikov equation is the cubic analogue of the Camassa–Holm and Degasperis–Procesi equations. Besides the momentum form \[ m_t+u^2m_x+3uu_xm=0,\qquad m=u-u_{xx}, \] it appears in periodic, nonperiodic, weakly dissipative, and nonzero-background formulations. In the periodic Cauchy problem it is also written in nonlocal form \[ u_t+u^2u_x=-A^{-2}\bigl(3u^2u_x+2u_x^3+3uu_xu_{xx}\bigr),\qquad A^2=1-\partial_x^2, \] while the nonzero-background inverse-scattering treatment uses \(\hat m=m+1=u-u_{xx}+1\) after reduction to zero background [1009.1820][1603.08842]. | Usage | Defining equation | Relation | |---|---|---| | Scalar Novikov equation | \(u_t-u_{txx}+4u^2u_x-3uu_xu_{xx}-u^2u_{xxx}=0\) | Central \(1+1\)-dimensional cubic CH-type equation [1009.1820] | | Associated Novikov equation | \(f_{xxt}-3(f_xf_t-1)=0\) | Related to Novikov by a chain of transformations [1905.02179] | | Novikov–Veselov equation | \(\partial_\tau q_\tau = - \partial_z^3 q_\tau - \bar\partial_z^3 q_\tau + 3\,\partial_z(q_\tau v_\tau) + 3\,\bar\partial_z(q_\tau \overline{v_\tau})\) | Unrelated \(2+1\)-dimensional KdV analogue [1105.3903] | | Novikov equations for commuting operators | Conditions equivalent to \([L,M]=[L,N]=0\) for \(L=D^3+pD+q\) | Stationary algebraic-geometric usage [2507.14978] | This multiplicity of usage is substantive rather than merely terminological. The scalar Novikov equation belongs to the peakon and shallow-water literature; the associated Novikov equation is a transform-related scalar representative with a third-order scalar spectral problem; the Novikov–Veselov equation belongs to \(2\)-dimensional inverse scattering; and the commuting-operator Novikov equations arise as stationary commutativity conditions in the theory of rank-one commutative rings of differential operators [1905.02179][1105.3903][2507.14978]. ## 2. Symmetry structure, self-adjointness, and local conservation laws The Lie-point symmetry analysis of the scalar Novikov equation yields a five-dimensional algebra with basis \[ X_1=\frac{\partial}{\partial t},\qquad X_2=\frac{\partial}{\partial x}, \] \[ X_3=e^{2x}\left(\frac{\partial}{\partial x}+u\frac{\partial}{\partial u}\right),\qquad X_4=e^{-2x}\left(\frac{\partial}{\partial x}-u\frac{\partial}{\partial u}\right), \] \[ X_5=-2t\frac{\partial}{\partial t}+u\frac{\partial}{\partial u}, \] together with the discrete symmetry \(x\mapsto -x\), which interchanges \(X_3\) and \(X_4\). Using Ibragimov’s formal Lagrangian \[ \mathcal L=v\left(u_t-u_{txx}+4u^2u_x-3uu_xu_{xx}-u^2u_{xxx}\right), \] the adjoint equation satisfies \[ F^*|_{v=u}=-F, \] so the equation is strictly self-adjoint. In this framework, among all Lie point symmetries only the dilation symmetry \(X_5\) yields a nontrivial local conserved vector, namely \[ D_t(u^2+u_x^2)+D_x\bigl(2u^4-2u^3u_{xx}-2uu_{tx}\bigr)=0 \] on solutions, whereas \(X_1,X_2,X_3,X_4\) give only trivial conservation laws. The same analysis produces explicit invariant solutions, including the infinite families \[ u(t,x)=p(t)e^x,\qquad u(t,x)=p(t)e^{-x}, \] dilation-invariant ansätze \(u(t,x)=t^{-1/2}\psi(x)\), and travelling waves \(u(t,x)=\phi(x-ct)\), with \(u(t,x)=e^{x-ct}\) as an explicit example [1202.3954]. This selective symmetry–conservation correspondence is characteristic of the scalar Novikov equation. A broader polynomial family \[ u_t-u_{txx}+au^p u_x-bu^{p-1}u_xu_{xx}-cu^p u_{xxx}=0 \] contains Camassa–Holm at \((p,a,b,c)=(1,3,2,1)\) and Novikov at \((2,4,3,1)\). Within that family, the joint requirements of \(N\)-peakons and the low-order conservation law \[ \int_{-\infty}^{\infty}(u^2+u_x^2)\,dx \] single out the generalized Camassa–Holm–Novikov equation \[ m_t+(p+1)u^{p-1}u_xm+u^pm_x=0,\qquad m=u-u_{xx}, \] which reduces to Novikov at \(p=2\) [1412.4415]. ## 3. Cauchy theory, global continuation, inverse scattering, and weak dissipation For the periodic Cauchy problem on \(\mathbb T=\mathbb R/\mathbb Z\), the Novikov equation is locally well posed in Sobolev spaces \(H^s(\mathbb T)\) for \[ s>\frac52, \] with solutions \[ u\in C\bigl([0,T),H^s(\mathbb T)\bigr)\cap C^1\bigl([0,T),H^{s-1}(\mathbb T)\bigr). \] A continuation criterion of Beale–Kato–Majda type is obtained from \[ \frac{d}{dt}\|u\|_{H^s}^2\lesssim \|u\|_{C^1}\|u\|_{H^s}^2. \] If \(s\ge 3\) and the initial momentum satisfies \[ A^2u_0=u_0-u_{0,xx}\ge 0 \qquad (\text{or } \le 0), \] then the solution exists globally and uniquely. In the analytic category, a Cauchy–Kowalevski type theorem yields short-time existence and uniqueness of solutions analytic in both \(x\) and \(t\) [1009.1820]. A different global theory is available in the \(\lambda\)-family \[ u_t-u_{txx}+(\lambda+2)u^\lambda u_x=(\lambda+1)u^{\lambda-1}u_xu_{xx}+u^\lambda u_{xxx}, \] whose case \(\lambda=2\) is the Novikov equation. For \(\lambda=2\), absolutely continuous initial data \[ u_0\in H^1(\mathbb R)\cap W^{1,4}(\mathbb R) \] generate global energy-conservative weak solutions that are Hölder continuous with exponent \[ 1-\frac{1}{2\lambda}=\frac34. \] These solutions preserve \[ \int_{\mathbb R}(u^2+u_x^2)\,dx \] for every \(t\ge 0\) and continue past wave breaking in a conservative weak-solution framework [2111.01030]. On the line with nonzero constant background, the inverse scattering transform has been formulated through a \(3\times 3\) matrix Riemann–Hilbert problem. After reduction to zero background and introduction of \[ q(x,t)=\hat m(x,t)^{1/3},\qquad \hat m=u-u_{xx}+1>0, \] the Cauchy problem is encoded in a singular \(3\times 3\) RH problem on a six-ray contour in the \(k\)-plane. The solution is reconstructed parametrically from the RH solution \(\hat M\) at the distinguished point \(k=e^{i\pi/6}\) through formulas for \(x(y,t)\) and \(\hat u(y,t)\). This framework exhibits the Novikov equation as a “modified DP equation” in the sense that its RH geometry is closely parallel to that of Degasperis–Procesi, but with a different nonlinear reconstruction [1603.08842]. The weakly dissipative Novikov equation, \[ v_t-v_{txx}+4v^2v_x-3vv_xv_{xx}-v^2v_{xxx}+\lambda(v-v_{xx})=0,\qquad \lambda>0, \] is exactly reducible to the non-dissipative equation by \[ v(t,x)=e^{-\lambda t}u\!\left(\frac{1-e^{-2\lambda t}}{2\lambda},x\right). \] Hence the dissipative flow is the non-dissipative flow viewed through exponential amplitude damping and the compressed time variable \[ s=\frac{1-e^{-2\lambda t}}{2\lambda}. \] This makes the weakly dissipative model structurally equivalent to the ordinary Novikov equation rather than an independent cubic evolution law [1207.0968]. ## 4. Peakons, characteristics, wave breaking, and stability The multipeakon ansatz \[ u(x,t)=\sum_{i=1}^N m_i(t)e^{-|x-x_i(t)|} \] plays the same role for Novikov that it does for Camassa–Holm and Degasperis–Procesi, but the cubic nonlinearity changes both characteristic transport and finite-dimensional dynamics. In the Novikov case, \[ \dot \xi(t)=u(\xi(t),t)^2, \] so the particle speed is the square of the wave elevation, and for peakons \[ \dot x_k=u(x_k)^2,\qquad \dot m_k=-m_k\,u(x_k)\,u_x(x_k). \] This implies \(\dot x_k\ge 0\), so both peakons and antipeakons move to the right. The explicit inverse-spectral formulas for ordinary \(N\)-peakons can be extended to all characteristic curves by introducing “ghostpeakons,” namely zero-amplitude peakons. In each interval between adjacent peakons, the characteristics are \[ \xi(t)=\frac12\ln\frac{Z_{k+1}+\theta Z_k}{W_k+\theta W_{k-1}},\qquad \theta>0, \] and the signed field along such curves is recovered by a separate determinant formula because the characteristic speed only gives \(u^2\), not the sign of \(u\) [1807.01910]. The Novikov equation also admits explicit unbounded solutions with finite-time peakon creation. A Lie symmetry analysis produces nonlinear transformations that send the one-peakon \(ce^{-|x-c^2t|}\) to \[ u_5(x,t)=c\sqrt{1+2\varepsilon e^{-2x}\; e^{-\left|\frac12\ln(e^{2x}+2\varepsilon)-c^2t\right|}}, \] which is smooth up to \[ t_0=\frac{1}{2c^2}\ln(2\varepsilon), \] when a peak is created at \(x=-\infty\). For \(t>t_0\), the solution becomes piecewise smooth, \[ u_5(x,t)= \begin{cases} c\,e^{-x+c^2 t}, & x\ge B(t),\\[1mm] c\,(e^x+2\varepsilon e^{-x})e^{-c^2 t}, & x< B(t), \end{cases} \qquad B(t)=\frac12\ln(e^{2c^2 t}-2\varepsilon), \] and remains a weak solution. This provides an explicit smooth-to-peaked transition inside the Novikov dynamics [1310.4927]. Stability theory bifurcates according to the class of solitary waves considered. For smooth solitary waves on nonzero background, there is a one-parameter family \(u(x,t)=\phi(x-ct)\) with asymptotic endstate \[ 0