Journal of Nonlinear Mathematical Physics

Volume 26, Issue 3, May 2019, Pages 390 - 403

Two-component generalizations of the Novikov equation

Authors
Hongmin Li
Schools of Mathematics, Huaqiao University, Quanzhou, Fujian 362021, P R China,lihongmin@hqu.edu.cn
Received 1 November 2018, Accepted 2 March 2019, Available Online 6 January 2021.
DOI
10.1080/14029251.2019.1613048How to use a DOI?
Keywords
Camassa-Holm type equation; Lax pair; Bi-Hamiltonian structure
Abstract

Some two-component generalizations of the Novikov equation, except the Geng-Xue equation, are presented, as well as their Lax pairs and bi-Hamiltonian structures. Furthermore, we study the Hamiltonians of the Geng-Xue equation and discuss the homogeneous and local properties of them.

Copyright
© 2019 The Authors. Published by Atlantis and Taylor & Francis
Open Access
This is an open access article distributed under the CC BY-NC 4.0 license (http://creativecommons.org/licenses/by-nc/4.0/).

1. Introduction

In recent decades, the Camassa-Holm (CH) type equations raised a lot of interest because of their specific properties, one of which is that they possess peakon solutions (peaked soliton solutions with discontinuous derivatives at the peaks). The most celebrated member of them is the CH equation [1]

mt+umx+2uxm=0,m=uuxx,(1.1)
which was found to be completely integrable with a Lax pair and associated bi-Hamiltonian structure [1, 2], and related by a reciprocal transformation to the first negative flow in the KdV hierarchy [18]. It is worthwhile to note that an inverse scattering approach for the CH equation was developed in the Refs. [5, 11, 14], which showed that after a suitable change of variables, the CH-flow becomes a linear flow at constant speed, this being the strong version of integrability (the infinite-dimensional counterpart of the classical Liouville theorem). Moreover, while an initial m without change of sign leads to solutions of CH that are defined for all times t ≥ 0 [3,8], among the initial m that change sign one encounters solutions of CH that model breaking waves, in the sense that the solution itself remains bounded but its slope becomes unbounded [10].

In 1999, Degasperis and Procesi used an asymptotic integrability approach to isolate integrable third-order equations, and discovered a new CH type equation, i.e., the Degasperis-Procesi (DP) equation

mt+umx+3uxm=0,m=uuxx,(1.2)
which turns out to be integrable with a bi-Hamiltonian structure and a Lax pair [16,17], and is reciprocal linked to a negative flow in the Kaup-Kupershmidt hierarchy [22]. Furthermore, the inverse scattering problem for the DP equation was formulated in Ref. [12]. Especially, for CH and DP equations one can prove that the peakons are orbitally stable [15, 25, 31], in the sense that their shape is stable under small perturbations, and therefore these patterns are detectable. This is an essential aspect since both equations arise as models for water waves [13, 23, 24]. Note that the travelling wave solutions of greatest height of the governing equations for water waves have a peak at their crest [6, 7, 9, 35]. The proof of the orbital stability of the peakons relies on the conservation laws for CH and DP.

Both the CH and DP equations are the third order CH type equations with quadratic nonlinearity. Recently, by the symmetry classification study of nonlocal partial differential equations with quadratic or cubic nonlinearity, Novikov discovered a new CH type equation with cubic nonlinearity [32]

mt+u2mx+3uuxm=0,m=uuxx.(1.3)

Subsequently, Hone and Wang presented a Lax pair for the Novikov equation (1.3) and proved it is related to the negative flow in the Sawada-Kotera hierarchy by a reciprocal transformation. Infinitely many conserved quantities and a bi-Hamiltonian structure of the Novikov equation are also constructed [21].

Later Geng and Xue [19] presented a two-component generalization for the Novikov equation

mt+3uxvm+uvmx=0,nt+3vxun+uvnx=0,m=uuxx,n=vvxx,(1.4)
which admits a Lax representation
φx=M¯φ,φt=N¯φ,(1.5)
where
M¯=(0mλ100nλ100),N¯=(13λ2uxvuxλuvmλuxvxvλ23λ2+uxvuvxuvnλvxλuvuλ13λ2+uvx).(1.6)

They also gave an infinite sequence of conserved quantities of the system (1.4) (also called Geng-Xue equation). Li et al proved it is bi-Hamiltonian [29] and reciprocal linked to a negative flow in the modified Boussinesq hierarchy [30]. Very recently, we constructed a Liouville transformation to connect it with another negative modified Boussinesq equation, and Lax pairs as well as bi-Hamiltonian structures of them are connected [27].

Recently, we make the vector prolongation of the Lax pair (1.5) as follows [28]

Φx=M˜Φ,Φt=N˜Φ,(1.7)
with
M˜=(0λQT10T0λR100),N˜=(UxTVUxTλλUTVQUxTVxVλINλ2+VUxTVxUTVxλλUTVRUTVUTλUTVx),(1.8)
where 0 and 0N are respectively N dimension row vector and N × N zero matrix and T is the transpose of a vector. IN denotes the N × N identity matrix and Q, R, U, V are the N-component column vector potentials, and UxT=UTx.

Then the zero-curvature equation for (1.7) yields the multi-component Novikov equation

Qt=2UxTVQUTVxQUTVQxQTVUx+QTVxU,Rt=2UTVxRUxTVRUTVRxRTUVx+RTUxV,Q=UUxxx,RVVxx,(1.9)
which is bi-Hamiltonian [28] proved by the multi-vector method [34]. This multi-component equation can be reduced to Geng-Xue equation (1.4), DP equation (1.2) and Novikov equation (1.3) under the constraints N = 1, Q = m, R = n, U = u, V = v; N = 1, Q = m, U = u, R = V = 1 and N = 1, Q = R = m, V = U = u respectively. The aim of this paper is to consider some reductions of the Lax pair (1.7) and the corresponding evolution equations.

This paper is organized as follows. In Section 2, we consider two interesting reductions of the multi-component Novikov equation (1.9) and construct the bi-Hamiltonian structures for them. In Section 3, we compute two infinite sequences of conserved quantities for the Geng-Xue equation, and discuss the homogeneous and local properties of the Hamiltonians of the Geng-Xue equation.

2. Two-component generalizations of the Novikov equation

In this section, we will study two reductions of the multi-component Novikov equation (1.9) and obtain two reduced systems, as well as their bi-Hamiltonian structures using Dirac reduction [33].

2.1. The first two-component generalization of the Novikov equation

Setting U = V, Q = R, the multi-component Novikov equation (1.9) is reduced to

Qt=3UxTUQUTUQxQTUUx+QTUxU,Q=UUxx,(2.1)
which can be transformed to a two-component Novikov equation
qt=(u2+v2)qx3(uux+vvx)q+r(uvxuxv),rt=(u2+v2)rx3(uux+vvx)r+q(uxvuvx),q=uuxx,r=vvxx.(2.2)
as Q = R = (q, r)T, U = V = (u,v)T. Especially when u = v, the equation (2.2) is reduced to the Novikov equation (1.3).

As pointed in [28], the equation (1.9) can be written in the bi-Hamiltonian form

(QR)t=𝒦(δH2δQδH2δR)=𝒥(δH1δQδH1δR),(2.3)
using the two compatible Hamiltonian operators
𝒦=(0N(21)IN(12)IN0N),(2.4)
𝒥=𝒥1+𝒥2,(2.5)
with
𝒥1=(32Q+Qx32R+Rx)(34)1(3QT+QxT3RT+RxT),𝒥2=(12Q1QT+(Q1QT)T12Q1RTQT1RIN12R1QTRT1QIN12R1RT+(R1RT)T),
and the two Hamiltonian functionals
1=12QTV+RTUdx,2=12QTVUxTVRTUVxTU+(RTUxQTVx)UTVdx.(2.6)

In what follows we will construct the bi-Hamiltonian structure for the multi-component Novikov equation (2.1) from (2.3). Introducing P=Q+R2, S=QR2, the new system with the variables P, S has a Hamiltonian operator:

(P,S)=14(ININININ)𝒥(ININININ),
then after using the Dirac reduction [33] under the constraint Q = R or S = 0, we have
Π(P,0)=12(Π11(P,0)Π12(P,0)Π21(P,0)Π22(P,0)),
where Π11)=(3Q+2Qx)(34)1(3QT+Q xT)+(Q1QT)TQT1QIN, Π12(P,0) = Π21(P,0) = 0, Π22(P,0) = Q∂−1QT + (Q∂−1QT)T + QT −1QIN.

So one may obtain the reduced Hamiltonian operator for the multi-component equation (2.1) as follows:

Πred(Q)=12(Π11(P,0)Π12(P,0)[Π22(P,0)]1Π21(P,0))=12Π01(P,0).(2.7)

To obtain the second Hamiltonian operator for the multi-component equation (2.1), we consider the Hamiltonian operator 𝒦 𝒥−1𝒦 and its Dirac reduction, then the following equality holds:

Λ(P,S)=14(ININININ)𝒦𝒥1𝒦(ININININ).

Similarly, after the Dirac reduction, one may get

Λ(P,0)=12(21)(Λ11(P,0)Λ12(P,0)Λ21(P,0)Λ22(P,0))1(21),
with Λ11(P,0) = Π22(P,0), Λ12(P,0) = Λ21(P,0) = 0 and Λ22(P,0) = Π11(P,0). So we have
Λ(P,0)=12(21)([Π22(P,0)]100[Π11(P,0)]1)(21),
which means the second Hamiltonian operator for the equation (2.1) is
Λred(Q)=12(21)[ Π22(P,0) ]1(21).(2.8)

Specially, from the constraint Q = R = (q, r)T, U = V = (u,v)T and the Hamiltonian pair (2.72.8), the two-component system (2.2) possesses the bi-Hamiltonian operators

Πred(q,r)=(3q+2qr3r+2rx)(34)1(3q+qx,3r+rx)+(rq)1(r,q),(2.9)
Λred(q,r)=12(21)(3q1q+r1rq1r+r1qq1r+r1qq1q+3r1r)1(21).

Remark 2.1.

The two-component system (2.2) appears in the bi-Hamiltonian form

(qr)t=Πred(q,r)(δHδqδHδr)=Λred(q,r)(δWδqδWδr),(2.10)
with Hamiltonian pair Πred(q, r), Λred(q, r) given by (2.9). The associated Hamiltonian functional is H=12qu+rvdx, and W is nonlocal and looks very complicated, so we omit it.

2.2. The second two-component generalization of the Novikov equation

Assuming N is an even number, we consider the constraint

Q=(Q1Q2),R=(R1R2),U=(U1U2),V=(V1V2),(2.11)
where R1 = Q2, R2 = Q1, V1 = U2, V2 = U1 and all of them are N2 dimension column vectors, then the multi-component Novikov equation (1.9) reduces to
Q1t=3(U1TU2)xQ12U1TU2Q1x+Q1T(U2xU1U2U1x)+Q2T(U1xU1U1U1x),Q2t=3(U1TU2)xQ22U1TU2Q2x+Q1T(U2xU2U2U2x)+Q2T(U1xU2U1U2x),Q1=U1U1xxQ2=U2U2xx.(2.12)

In particular, setting Q = (q, r)T, R = (r, q)T, U = (u, v)T, V = (v, u)T, the multi-component equation (2.12) leads to another two-component Novikov equation

qt=2uvqx2uvxq4uxvq,rt=2uvrx2uxvr4uvxr,q=uuxx,r=vvxx,(2.13)
which may be reduced to the Novikov equation (1.3) under the constraint u = v.

To construct the bi-Hamiltonian structure of the multi-component system (2.12), we apply the same procedure in subsection 2.1 and introduce these symbols P1=Q1+R22, P2=Q2+R12, S1=Q1R22, S2=Q2R12, then the Hamiltonian operator Π(P1, P2, S1, S2) for the new system with the variables P1, P2, S1, S2 can be written as:

(P1,P2,S1,S2)=14𝒥*,
where =(I00I0II0I00I0II0), 0 and I are respectively N2×N2 zero matrix and N2×N2 identity matrix.

After using the Dirac reduction under the constraint Q1 = R2, Q2 = R1 or S1 = S2 = 0, the operator Π(P1, P2, S1, S2) reduces to

Π(P1,P2,0,0)=(Π11(P1,P2,0,0)Π12(P1,P2,0,0)00Π21(P1,P2,0,0)Π22(P1,P2,0,0)0000Π33(P1,P2,0,0)Π34(P1,P2,0,0)00Π43(P1,P2,0,0)Π44(P1,P2,0,0))
with
Π11(P1,P2,0,0)=12(3Q1+2Q1x)(34)1(3Q 1T+Q 1xT)+12(Q11Q 1T)T,Π1212(Q 1T1Q2+Q 2T1Q1)I.Π21(P1,P2,0,0)=12(3Q2+2Q2x)(34)1(3Q 1T+Q 1xT)+12(Q11Q 2T)T12(Q 1T1Q2+Q 2T1Q1)I,Π22(P1,P2,0,0)=12(3Q2+2Q2x)(34)1(3Q 2T+Q 2xT)+12(Q21Q 2T)T,Π33(P1,P2,0,0)=12(Q11Q 1T+(Q11Q 1T)T,Π34(P1,P2,0,0)=12(Q11Q 2T+(Q21Q 1T)T)+12(Q 1T1Q2+Q 2T1Q1)I,Π43(P1,P2,0,0)=12(Q21Q 1T+(Q11Q 2T)T)+12(Q 2T1Q1+Q 1T1Q2)I,Π44(P1,P2,0,0)=12(Q21Q 2T+(Q21Q 2T)T).

And then the reduced Hamiltonian operator for the multi-component Novikov equation (2.12) is

Πred(Q1,Q2)=(Π11(P1,P2,0,0)Π12(P1,P2,0,0)Π21(P1,P2,0,0)Π22(P1,P2,0,0)).

Moreover, by studying the Hamiltonian operator 𝒦 𝒥−1𝒦 and its Dirac reduction, the second Hamiltonian operator for the multi-component system (2.12) may be obtained as

Λred(Q1,Q2)=14(21)(Π44(P1,P2,0,0)Π43(P1,P2,0,0)Π34(P1,P2,0,0)Π33(P1,P2,0,0))1(21).

Especially under the constraint Q = (q, r)T, R = (r, q)T, U = (u, v)T, V = (v, u)T, the two-component system (2.13) has a Hamiltonian pair

Πred(q,r)=12(3q+2qx3r+2rx)(34)1(3q+qx3r+rx)+12(q1qq1rr1qr1r),Λred(q,r)=14(21)(r1rq1+r1qq1r+r1qq1q)(21)=12(21)(1qqr1qr1q1q(2(qr)x1rq)1r1r(2(rr)x1qr)1q1rrq1rq1r)(21).(2.14)

Remark 2.2.

The two-component system (2.13) can be written in bi-Hamiltonian form (2.10) with bi-Hamiltonian operators given by (2.14). The associated Hamiltonian functional for Πred(q, r) is H = ∫qv+rudx. The omission of W here is caused by the nonlocal property and complication of it.

3. Conserved quantities and locality of the Hamiltonians of the Geng-Xue equation (1.4)

In this section, we give two sequences of conserved quantities for the Geng-Xue equation (1.4), some of them seem to be new. Furthermore, we find all of them are homogeneous and all the Hamiltonians Hj, j ≤ −1 are local.

3.1. Conserved quantities

Defining

φ=(φ1,φ2,φ3)T,a=φ1φ2,b=φ3φ2,(3.1)
based on the Lax pair (1.5), we have
ρ=(lnφ2)x=λnb.(3.2)

The infinitely many conserved quantities for the Geng-Xue equation (1.4) can be constructed using the standard algorithm. Next we will utilize the procedure in Ref. [30] to obtain the required Hamiltonians.

Substituting these equalities (3.1) to the spectral problem (1.5) yields

ax=λm+baρ,bx=abρ,(3.3)
after expanding a, b into a series of the spectral parameter λ and equating the coefficients of different powers of λ in these equations (3.3), one can find a series of conserved densities from coefficients of ρ in power of λ.
  • Case 1: Expanding a, b as a = ∑i≥1 aiλi, b = ∑j≥1 bjλ j and substituting them into the system (3.3), we establish that

    a1=ux,b1=u,a2=0,b2=0,(12)b3=(nu2)x+nuux,a3=b3x+nb12,,(12)ak=i,j1i+j=k1[nbibj+(naibj)x],(12)bk=i,j1i+j=k1[naibj+(nbibj)x].

    Hence an infinite sequence of conserved quantities are obtained, and the first two are

    H1=undx,H2=(uxvuvx)undx.(3.4)

    Remark 3.1.

    The above two conserved quantities can be obtained from the Hamiltonian functionals (2.6) 1, 2 under the constraint N = 1, Q = m, R = n, U = u, V = v respectively.

  • Case 2: The second expansion is a=i0aiλ12i3, b=j0biλ2j+13. Substituting it into the system (3.3) yields:

    a0=m23n13,b0=m13n23,a1=13(mn)23mx,b1=13m1n2(mnxmxn),ak+1=13[2(mn)13(bk1akk)+bkk+i,j1i+j=k+1(nbibj2b01aibj)],bk+1=13(mn)23(bk1akx)13(mn)13[bkk+i,j1i+j=k+1(b01ai+nbi)bj].

    In particular, we only consider the conserved quantities associated with b3k, k ∈ Z, and obtain the first two conserved quantities as

    H1=3(mn)13dx,H2=127[5mx2nxn53m833(mn)53mxxnx18nxm23n53]dx.(3.5)

    The negative powers of m and n that are involved in the above expressions for the conserved quantities raise an important structural property, whose solution established their analytic validity. Namely, we expect that if m > 0 and n > 0 hold initially (at time t = 0), then this property is preserved by the Geng-Xue flow. That this might be so is suggested by the granting, in the context of the CH equation, of such a property by means of a pointwise conservation law that has in that context very strong consequences for the long-time behavior of solutions [4, 20]. For the Geng-Xue flow we can proceed as follows. Define the smooth diffeomorphisms ψ(x,t) by ψ(x,0) = x and ψt(x,t) = u(ψ(x,t),t)v(ψ(x,t),t) for t > 0 and consider the expressions Q¯(x,t)=m(ψ(x,t),t)ψx3(x,t) and R¯(x,t)=n(ψ(x,t),t)ψx3(x,t). Since

    ψxt=[ux(ψ,t)v(ψ,t)+u(ψ,t)vx(ψ,t)]ψx,t>0,
    and ψx(x, 0) = 1 yield
    ψx(x,t)=exp(0t[ux(ψ,s)v(ψ,s)+u(ψ,s)vx(ψ,s)]ds)>0,t0,
    using (1.4) we can verify that the relations
    tQ¯(x,t)=3u(ψ(x,t),t)vx(ψ(x,t),t)ψx3(x,t)Q¯(x,t),tR¯(x,t)=3ux(ψ(x,t),t)v(ψ(x,t),t)ψx3(x,t)R¯(x,t),(3.6)
    hold for any fixed x. Therefore
    Q¯(x,t)=Q¯(x,0)exp(30t[u(ψ(x,s),s)vx(ψ(x,s),s)ψx3(x,s)ds),t0,R¯(x,t)=R¯(x,0)exp(30t[ux(ψ(x,s),s)v(ψ(x,s),s)ψx3(x,s)ds),t0.(3.7)

    These considerations prove that m(x, t) > 0 and n(x, t) > 0 at any t ≥ 0, provided that m(x, 0) > 0 and n(x, 0) > 0.

3.2. Homogeneous and local properties of the Hamiltonians

The first nontrivial negative flow in the Geng-Xue hierarchy is the Geng-Xue equation (1.4) [19]. It admits a bi-Hamiltonian structure [29], viz:

θt=KδH2δθ=JδH1δθ,θ=(mn),(3.8)
where H1, H2 are given by (3.4), and K, J are respectively the reductions of the Hamiltonian operators 𝒦, 𝒥 (2.4), (2.5) under the constraint N = 1, Q = m, R = n, U = u, V = v, namely
K=(021120),J=(J11J12J21J22),
with
J11=32m1m+(mx+32m)(34)1(mx+3m),J12=32m1n+(mx+32m)(34)1(nx+3n),J21=32n1m+(nx+32n)(34)1(mx+3m),J22=32n1n+(nx+32n)(34)1(nx+3n),

Furthermore, the first nontrivial positive flow in the hierarchy reads

mt=(m13n23)xxm13n23,nt=n13m23(n13m23)xx,
which can also be reformulated as a bi-Hamiltonian form
θt=JδH2δθ=KδH1δθ,(3.9)
herein H−1, H−2 are given by (3.5).

Hence, the Geng-Xue hierarchy may be constructed as

θtj=JδHjδθ=KδHj+1δθ,jZ,j0,1,(3.10)
with infinitely many Hamiltonian functionals Hj generated by the recursive definition. It is immediately clear that all Hjs are homogeneous and Hjs, j ≥ 3 are nonlocal. Hence we will discuss the local property of Hj, j ≤ −1.

Lemma 3.1 ([26,34]).

If a differential function M^[θ] satisfies

M^[θ]dx=0
for all θ, then there exists a unique differential function N^[θ] up to addition of a constant such that M^[θ] is the total x-derivative M^[θ]=(N^[θ])x.

Theorem 3.1.

Let Xj[θ]=δHj[θ]δθ, then for each j ≤ −1, Xj[θ] and Hj[θ] are local.

Proof.

We only need to demonstrate the local property of Xjs because of Hj[θ]=01Xj[ɛθ]θdɛ. As j = −1, X−1[θ] is local since

X1[θ]=(m23n13,m13n23)T.

Now suppose that Xj[θ] is local for j = −k. When j = −(k + 1), we have

X(k+1)[θ]=J1KXk[θ]=(J1K)kX1[θ].(3.11)

To find an explicit expression for Xk[θ], a direct choice is to calculate the inverse operator of J. However, it seems to be difficult. Hence, we find an indirect way. Suppose h=(mn)13 and define

Xk[θ]=(Ak,Bk)T,Ek=(34)1(3m+mx,3n+nx)Xk[θ],
it follows from the recursive relation (3.11) that
32m1(mAk+1nBk+1)+(32m+mx)Ek+1=(21)Bk,32n1(mAk+1nBk+1)+(32n+nx)Ek+1=(12)Ak.

Eliminating Ek+1 from the above system and combating the definition of Ek+1 yield

(3m+mx3n+nxmn)(Ak+1Bk+1)=13Θ(AkBk),(3.12)
where
Θ=(0002m(21))+(432+2mx3m)1h1(m13n23m23n13)(12).

Then it is not difficult to solve Ak+1 and Bk+1 from the above system, and we find that Ak+1 and Bk+1 are local if there exist two differential functions Mk, Nk such that

Γ1=1h2[m(12)Ak+n(21)Bk]=Mkx,(3.13)
Γ2=1h([34+3(2mx3m)(+2mx3m)]Mkh(2mx3m)6m(21)Bk)=Nkx.(3.14)

Using the Lemma 3.1, the equality (3.13) follows immediately from

Γ1dx=1h2[m(12)Ak+n(21)Bk]dx=X1TKXkdx=X1TK(J1K)k1X1dx=X1T(KJ1)k1KX1dx=X1TK(J1K)k1X1dx=X1TKXkdx,
where we have used the skew-symmetric property of K, J. So ∫Γ1dx = 0. Similarly, the equality (3.14) may be demonstrated from
Γ2dx=1h[34+(32mxm)(2+2mx3m)]Mkh1h(32mxm)2m(21)Bkdx=f(12)Ak+g(21)Bkdx=(g,f)K(Ak,Bk)Tdx=(Ak,Bk)K(g,f)Tdx=J1K(Ak,Bk)TJ(g,f)Tdx=0,
where
f=1h4(23mxx43mnxxn109mx2m+149mnx2n2+49mxnxn2m),g=1h4(23nxx43nmxxm109nx2n+149nmx2m2+49mxnxn2n).

In the above we have eliminated the derivatives of Mk in the integration by the equality (3.13) and used the skew-symmetric property of K, J, as well as the fact that (g, f)T is a kernel of J. So all the Hamiltonians Hj, j ≤ −1 for the Geng-Xue equation (1.4) are local.

Acknowledgements

This work is partially supported by the National Natural Science Foundation of China (Grant Nos. 11805071 and 11747010) and the Initial Founding of Scientific Research for the introduction of talents of Huaqiao University (Project No. 16BS513).

References

[31]Y. Liu and Z. Lin, Stability of peakons for the Degasperis-Procesi equation, Comm. Pure Appl. Math, Vol. 62, 2009, pp. 125-146.
[34]P.J. Olver, Applications of Lie Groups to Differential Equations, Springer, Berlin, 1993.
Journal
Journal of Nonlinear Mathematical Physics
Volume-Issue
26 - 3
Pages
390 - 403
Publication Date
2021/01/06
ISSN (Online)
1776-0852
ISSN (Print)
1402-9251
DOI
10.1080/14029251.2019.1613048How to use a DOI?
Copyright
© 2019 The Authors. Published by Atlantis and Taylor & Francis
Open Access
This is an open access article distributed under the CC BY-NC 4.0 license (http://creativecommons.org/licenses/by-nc/4.0/).

Cite this article

TY  - JOUR
AU  - Hongmin Li
PY  - 2021
DA  - 2021/01/06
TI  - Two-component generalizations of the Novikov equation
JO  - Journal of Nonlinear Mathematical Physics
SP  - 390
EP  - 403
VL  - 26
IS  - 3
SN  - 1776-0852
UR  - https://doi.org/10.1080/14029251.2019.1613048
DO  - 10.1080/14029251.2019.1613048
ID  - Li2021
ER  -