Orbital stability of the trains of peaked solitary waves for the modi ﬁ ed Camassa-Holm-Novikov equation

: Consideration herein is the stability issue of peaked solitary wave solution for the modi ﬁ ed Camassa-Holm-Novikov equation, which is derived from the shallow water theory. This wave con ﬁ guration accommodates the ordered trains of the modi ﬁ ed Camassa-Holm-Novikov-peaked solitary solution. With the application of conservation laws and the monotonicity property of the localized energy functionals, we prove the orbital stability of this wave pro ﬁ le in the (cid:2) ( ) H 1 energy space according to the modulation argument.


Introduction
In this article, we consider the modified Camassa-Holm-Novikov (mCH-Novikov) equation arising from the water wave theory with constant vorticity [5], namely, (( ) ) ( ) (1.1) where = − m u u xx , ( ) u t x , is related to the horizontal velocity of the fluid: k 1 and k 2 are two arbitrary constants.This model was derived by Chen et al. [5] from the Euler equation for an incompressible fluid with simple bottom and surface conditions.Indeed, it is obtained with the application of suitable scaling, the truncation of the asymptotic expansions of the unknowns to an appropriate order, and the introduction of a special Kodama transformation.With the property of transport equation, the local well-posedness of the mCH-Novikov equation is given in in the study by Mi et al. [20].It is shown in the study by Chen et al. [4] that this model possesses a peaked traveling wave solution, which is orbital stable in both ( ) . Moreover, a blowup criterion is given in in the study by Chen et al. [5] by using the sign preservation of the momentum density and a Riccati-type inequality.
It is observed that the mCH-Novikov equation (1.1) becomes the modified Camassa-Holm (mCH) equation (1. (1.5) are two conserved quantities.Similar to the CH equation, some new features are captured by the mCH equation, which has attracted the attention of numerous mathematicians and physicians in the past decades, including the wave breaking phenomena, blow-up criteria, and peaked solitary wave [7,8,27,28] that differed from the traditional smooth wave pattern.Indeed, the mCH equation possesses the single peakon [12] ( (1.6) and multipeakons where ( ) p t i and ( ) q t i satisfy the corresponding Hamiltonian system [14].It also possesses a train of infinite many peaked solitary waves [3].In the study by Qu et al. [25], with the construction of certain Lyapunov functionals, the single peakon soliton of the mCH equation is orbital stable under small perturbations in H 1 energy space.Li and Liu also establish the orbital stability of peakons in the Sobolev space ∩ H W 1 1 , 4 [17] using the variational approach without the assumption on the positive momentum density initially.Inspired by the work of Qu et al. [25], the stability of ordered trains of peakons is proved by constructing two suitable functionals related to the two conserved quantities [18].On the other hand, taking = k 0 which is another well-known CH-type equation with cubic nonlinearity.The Novikov equation was originally proposed from the classification of generalized CH-type equations, which possess infinite hierarchies of higher symmetries in [21].It can also be derived from the shallow water wave theory [5].Note that the Novikov equation also carries a bi-Hamiltonian structure, which can be written as follows: (1.9) where (1. (1.11) are two conserved quantities.The integrability, well-posedness, blow-up criteria, and existence of peakons to this model have been investigated in [13,26,[29][30][31].More specifically, the single peaked solution to the Novikov equation was derived in [11,15,16], i.e., (1.12) due to the interaction between the nonlinear and dispersion terms.The multipeakon solution is given by Hone et al. with the inverse scattering approach in [15].Using the conservation laws, the orbital stability of peakons for the Novikov equation in the sense of the energy space H 1 norm is given in the study by Liu et al. [19].Moreover, by using the modulation argument, the exponential decay property of H 1 -localized solutions, and the rigidity property, the asymptotic stability for the Novikov peakon is given in the studies by Chen et al. [6] and Palacios [23].
As a successor to the mCH equation and the Novikov equation, equation (1.1) also enjoys the peaked traveling wave solution c x ct (1.13) where a is related to c, k 1 , and k 2 .A direct computation shows that the mCH-Novikov equation retains a linear superposition of N -peaked solitary wave solutions with restrictions on the position functions and the ampli- tude functions [4].Thus, inspired by the aforementioned work, we are wondering the dynamical properties of the mCH-Novikov solitary wave solution has how much in common with the mCH equation and the Novikov equation.More precisely, we are interested in the orbital stability issue on the wave profiles under small perturbation for the mCH-Novikov (equation (1.1), which has significant physical implication.In the present study, the idea for investigating the mCH-Novikov-peaked solitary wave solution developed from [9] and [18], and hence we will recall the approach briefly and reveal the main challenges.
The general framework for the proof of the stability of the N -peaked solitary waves has three principal ingredients: the first one is the modulation argument, which provides the orbits of each traveling wave following the implicit function theorem; the second one is the almost monotonicity of the functionals relating to the conservation laws that describe the energy at the right of ith bump, for = i N 1, 2,…, ; and the third one is the local and global estimates that connect the local maximum value of u with localized energy functionals i and i .
The case for the mCH-Novikov equation is considerably similar to the mCH equation and the Novikov equation.However, note that the expression of the mCH-Novikov equation contains interaction between those cubic nonlinearities.This presents a series of obstacles when applying Dika and Molinet's method.The conservation laws of the mCH-Novikov equation include which are exactly the same as two conserved quantities in mCH equation.Using these two conservation laws may still be able to give the local and global estimates.However, it is not clear whether the appearance of the interaction between those cubic nonlinearities (( and ( ) with coefficients k 1 and k 2 would fit the framework of the modulation argument that generates the orbits, and carry the same almost monotonicity of the localized energy functionals or not.Moreover, the existence of a peaked traveling wave solution depends on the parameters k 1 and k 2 , so prior to the stability of the N -peaked solitary waves, one shall place restrictions on these two parameters.
We now explain the idea and the strategy adopted here to overcome the above difficulties.To prove the orbital stability of the trains of ordered N -peaked solitary waves to the mCH-Novikov equation has the flavor of the Lyapunov sense.In order to apply the modulation argument, we have to estimate the ( ) . Combining all the like terms, using the properties of the kernel function and a subtle analysis, we fulfill the first gap mentioned above.On the other hand, with the help of weighted functions and the tricks adopted in both the mCH equation, and Novikov equation, we achieve the almost monotonicity for the mCH-Novikov equation.This combining with local and global estimates implies the stability result according to the idea originated in [9] and [18].Now, we are ready to state the orbital stability result on the ordered N -peaked solitary waves , where > ∕ s 5 2 be the initial data of the mCH-Novikov equation and be the corresponding solution of the Cauchy problem for the mCH-Novikov equation , and > 0 such that for some < < κ κ 0 0 and then there exist ( ) where ( ) , = i N 2,…, and T is the maximum existence time depending only on initial data u 0 .
The remainder of the article is organized as follows: in Section 2, the local well-posedness, the existence of peaked solitary waves, and the orbital stability of single peaked solitary wave solution for the mCH-Novikov equation are recalled; and then, efforts are devoted to the modulation arguments, the almost monotonicity, and the local and global estimates in Section 3, which prove the orbital stability of ordered trains of N -peaked solitary wave solution to the mCH-Novikov equation.

Preliminary
This section gives a brief review on the existence and stability of peaked solitary waves solution for the mCH-Novikov equation (1.1).Prior to a discussion of the stability issue should be the well-posedness result.

Local well-posedness and conservation laws
Consider the following Cauchy problem, namely, , . where is the kernel function for the operator ( ) − ∂ − 1 (2.3) Now, we recall the following local well-posedness result without proof, which is rudiment of the stability theory established here.For details, refer to the study by Mi et al. [20].
. Furthermore, the mapping ↦ u u .
As aforementioned, the conservation laws play an essential role while proving the stability result.Hence, we also recall the following result.Proposition 2.2.[4] For the strong solutions u obtained in Proposition 2.1, the following functionals (2.4) are conserved with respect to time t.
Moreover, we recall the existence of the peaked solution to equation (1.1), which discovers the relationship between wave amplitude and wave speed.
is a peaked solution to equation (2.3), provided that , and = ± ≠ + a 0 ; or and ≠ a 0.
A direct computation reveals that (2.6) Note that, due to the emergence of singularity, the peaked solitary wave solution should be considered as a weak solution in the distribution sense.Moreover, the H 1 -orbital stability of the single peaked wave is also given in the study by Chen et al. [4].
x ct be the peaked solutions given in Proposition 2.3.Assume that the initial data Then, φ c is H 1 -orbital stable in the following sense: there exists then the corresponding solution ( ) where ( ) * r t is the point at which the solution ( ) u t x , attains its maximum and the constant Note that with a refined continuity argument, the H 1 -orbital stability of the single peaked solitary wave solution of the mCH-Novikov equation removes the sign preservation condition on the initial data.Moreover, the proof of this proposition is essentially based on the following four lemmas.x x (2.11) 4 . ( (2.12) where Lemma 2.4.[4] Assume that ( ) , which satisfies equations (2.12) and (2.13) with ≤ ≪ δ 0 1.Then, we have the following: (2.14) (2.15) These lemmas give a clue for proving the orbital stability of trains of N -ordered peaked solitary wave solutions.
3 Orbital stability of the sum of ordered trains of peaked solitary waves On the one hand, solitons traveling with their own speed will finally be arranged in order.On the other hand, the soliton resolution conjecture says solutions to the dispersive equation with certain initial data should eventually resolve into a finite number of soliton plus a radiative term that goes to zero.Therefore, we mainly concern the orbital stability of the trains of N -ordered peaked solitary wave solution in the present article.
Inspired by the strategy developed in the study by El Dika and Molinet [9], we carry out the framework of proof to Theorem 1.1 in the following.
For > κ 0 and > R 0, the neighborhood of ( ) with radius κ and spatial shifts r i satisfying − ≥ , where = i N 2,…, , is denoted as follows: To prove the orbital stability of the ordered trains is equivalent to showing that there exist > 0, > R 0 0 , and > κ 0 , then for all [ ) ∈ t T 0, , where T is the maximal existence time, the corresponding solution ( ) . We would like to point out that is independent of time t.
Due to the fact that the mapping and if for some then Therefore, assume that equation (3.1) holds for some < < κ κ 0 0 and > R R 0 , where , κ 0 , and R 0 will be determined later, we are going to demonstrate equation (3.2) in the sequel.

Modulation argument
In the following lemma, with the application of orthogonal condition and conservation laws, we will show that the solution u of the mCH-Novikov equation satisfying equation (3.1) has the following properties: ( ) u t remains close to the sum of N -modulated peaked waves in ( ) H 1 , where each peaked wave has no overlap for [ ] ∈ t t 0, * .More precise, solution u can be decomposed into the following form: and the corresponding solution ( ) u t satisfies where < < t T 0 * , then there exist unique C 1 spatial functions ( ) one has the orthogonal condition and the following estimates , with Stability trains peaked solitary waves mCH-Novikov equations  7 Let ( ) * r t i be any point in the spatial line such that then there holds Proof.The standard modulation argument enable us to discover the translations of N -peaked solitary waves. where To determine the spatial shift functions f i , one shall apply the implicit function theorem, which requires two preconditions satisfied.The first one is that function should be C 1 mapping.This can be proved by the dominated convergence theorem.The second condition is that the determinant of the matrix of the first partial derivatives at ( ) 0, …,0, Φ y should be nonzero.For d , where .
Thus, for sufficiently large Hence, with the application of the implicit function theorem we know, there exist < < β κ 0 0 0 and uniquely determined where K 0 and β 0 depend on k 1 , k 2 , c 1 , and , 2 , which can be covered by In addition, the uniqueness in the implicit function theorem carries out the uniqueness of the modulation of u.Therefore, one shall be able to define the modulation of the solution ( ) , which gives rise to equation (3.4).
Inspired by the translation ( ) r t ˜i defined above and the estimates of ( ( )) Attention is now turn to the rate of change for ( ) x c c j j 2 j j (3.13)For the orthogonality condition, taking derivative with respect to time t, we derive that which, when combined with equation (3.13), implies that On the other hand, substituting u by Stability trains peaked solitary waves mCH-Novikov equations  9 For each c i , φ ˜ci is also a solution to the mCH-Novikov equation (2.3), which implies that Then, the aforementioned two identities infer that v satisfies the following condition on [ ] t 0, * , , Then, taking L 2 -inner product with ∂ φ x cj , the above identity deduces where the technique of integration by parts and the properties of exponential decay for both φ ˜ci and its first order derivative are applied.Now, the estimate of Then, V 1 has the following form: With the identity (3.13), V 1 can be rewritten as follows: which, when combined with the embedding formula and the exponential decay of φ ˜cj , gives rise to Second, for the estimate of V 2 , a direct expansion infers With the sign-persistence property of solution u and equation (3.18), we know Stability trains peaked solitary waves mCH-Novikov equations  11 In the similar fashion, with Holder's inequality, we obtain Moreover, using equation (3.12) and the exponential decay of ∂ φ x ci , there holds Hence, we obtain Third, with the expression of Q 1 and integration by parts, V 3 can be written in the following form: The application of equations (3.19) and (3.20) infers that For the estimation V 4 , ( ) With the fact that , we have Then, with equation (3.12) and the exponential decay of φ ˜ci , we deduce , which then gives rise to For the estimate of V 5 , we know i i .Similar to the process in the estimation of V 2 , we know Hence, this deduces that Note that the constants involving in V V V V , , , 1 2 3 4 , and V 5 depend on k k , 1 2 , and ( ) = c i i N 1 .On the other hand, note that by the inequality (3.12), the construction of the modulation guarantees ‖ ( )‖ ( ) . With the application of equations (3.12) and (3.14) and the decay of ∂ φ x ci , we obtain that Taking κ 0 small enough and R 0 large enough depending on k 1 , k 2 , and 2 , j j j j j j j j j j j j j j j j which yields the estimate (3.7).Inspired by equations (3.5) and (3.7), taking .

Monotonicity property
In this subsection, with the application of two conservation laws and the weighted function, we will show the almost monotonicity of functionals that describe the energy at the right side of ith peaked wave of u, and , we carry out the functionals that describe the energy at the right side of ith peaked wave and is close to ‖ ( )‖ ( ( )) .Let > T 0 be the maximal existence time of the corresponding strong solution ( ) u t x , .Then, for arbitrary smooth function ( ) x Θ , the following identity holds: Using integrating by parts gives Adopting equation (1.1), one has Moreover, it can be deduced from equation (3.30) that Similarly, substituting = − m u u xx into the expressions of 1 and 2 , we have Therefore, combining 1 and 2 , one has which proves equation (3.29).□ Consider the weighted function ϒ j K , defined in equations (3.27) and (3.28).For = j N 2,…, , we introduce In the following lemma, we shall show that the functional representing the energy at the right of the − j 1th bump of the solution ( ) u t x , is almost decreasing with respect to time.In the sequel, denote min , , …, .
, is the corresponding strong solution of the mCH-Novikov equation.There exist Proof.The proof of this lemma is carried out by the approach developed in the study by El Dika and Molinet [9].Fix the index , one obtains Stability trains peaked solitary waves mCH-Novikov equations  17 Recall the bound of ( ) r t ˙as equation (3.6) in Lemma 3.1.Then, for ( ) ( ( ) , we have the following estimate: Now, we are going to estimate ( ) . With the sign preservation of momentum density function ( ) m t x , and the fact that (3.32) For further discussion, we divide by = ∪ D D j j c , where which, when combined with equation (3.7) and the definition of χ 0 , infers Moreover, with the fact that which together with the definition of ϒ gives rise to, for ∈ x D j c , ( ) ( ) Therefore, with restriction on D j c , inspired by the conservation law On the other hand, for any ∈ x D j and the index Due to the decay property of ( ( )) ⋅− φ r t c i i and the estimate (3.5), one has .
Hence, combining the estimates on two subregions together, i.e., (3.33) and (3.35), for Moreover, by the Sobolev embedding theorem and the properties of the kernel function, we know where σ 1 and σ 2 are arbitrary constants.According to the definition of the weighted function

, we know
Stability trains peaked solitary waves mCH-Novikov equations  Picking ≥ K 4, the above inequality infers . Therefore, on the subregion D j , we have where σ 1 and σ 2 are arbitrary constants.With the estimate for v in equation (3.5) and the application of conservation law ( ) u , one has which along with equation ( Thus, for , 0 4 Hence, we have in equation (3.37), respectively.With the same pattern, discussing the integrals in two subregions, one has , 0 4 Moreover, by taking advantage of integration by parts and the properties of weighted function, we obtain and any [ ] ∈ t t 0, * with a positive constant C. Finally, with the application of the Gronwall's inequality, we deduce which completes the proof of Lemma 3.3.□

Local and global energy estimates
With the definition and properties of the weighted functions given by equations (3.27) and (3.28), we shall define the bump functions as follows: , we derive the following properties for the bump functions: With the aforementioned properties of bump functions ( ) t x Ψ , i , we shall introduce the following localized conservation laws of i and i , where First, we shall present a estimate describing the distance between ( ) u x and The proof of the global estimate is almost identical to that in ( [18], Lemma 3.1), and hence we omit it. .

Assume that ( ) u t x
, is the solution of the mCH-Novikov equation given by Proposition 2.1.Then, for fixed time t, there holds (3.43) where the constants involving in In the next lemma, we provide the relationship between the localized conservation laws of i and i and the localized maximum value ( ) . For the sake of completeness, its proof will be given in Appendix A.1. . 8).Assume that ( ) u t x , is the solution of the mCH-Novikov equation given by Proposition 2.1.For = i N 1,…, , denote the local maximum value of ( )

Consider the interval
. Then, there holds .
1 2 (3.45) Furthermore, we shall estimate the difference between the local maximum of the solution ( ) u t x , and the maximum of each single peaked solitary wave.The approach is inspired by the study of Liu et al. [18], and the detailed proof is given in Appendix A5.
where the constants in ( ) where , where ( ) * r t* i is defined in equation (3.9), Then, taking advantage of the global estimate (3.43) and the local estimate (3.46), it is deduced that In consequence, for where the depends only on , and ‖ ‖ u H 0 s , not on .Hence, we may choose = ∕2, which then completes the proof of Theorem 1.1.
x i x i Then, one has With an application of integration by parts and the localized energy function, one has which, when combined with the exponential decay property of the bump function Ψ i , gives rise to Furthermore, to connect the two localized conservation laws, we introduce another auxiliary function ( ) ζ x in the following form: , .
It satisfies Using integrating by parts and exponential decay Ψ i , one has In the same fashion, we derive Stability trains peaked solitary waves mCH-Novikov equations  25 which, when combined with the expression of II 1 , gives rise to On the other hand, with the sign preservation of momentum density, we know Therefore, together with the definition of local maximum value, one shall obtain which carries out the following estimate with the identity (A2) Considering both equations (A3) and (A4), we obtain that Based on the construction of the bump function Ψ i , one has Moreover, with the estimation of equation (3.10) From the fact that ( ) it is deduced from the Minkowski inequality that 1 which together with equation (3.3) implies that By the relationship between the localized conservation laws of i and i and the localized maximum value ( ) Taking summation with respect to i from 1 to N for the above identity, one obtains Second, for III 3 , following with a same fashion developed in Lemma 3.5 in [25], we obtain   Finally, according to the estimates on III 1 , III 2 , and III 3 , we conclude that ( ) ( ) which then infers that This completes the proof of Lemma 3.6. □

1 and = k 1 2 ,
the mCH-Novikov equation reduces to the Novikov equation

x 2 1 1 x 2 1
and " * " represents the convolution.With the application of the operator ( ) − ∂ − , the first equation in equation (2.1) can be rewritten as follows:

Remark 2 . 1 .
Even though the mCH-Novikov equation does not admit conserved total mass, with the observation of the corresponding characteristic equation, it still enjoys the sign-persistence property.That is, if the initial data ), then the corresponding solution satisfies (

2 .depending on k 1 ,
Let u 0 be the initial value to the mCH-Novikov equation.There exist > k 2 , and

jLemma 3 . 2 .
defined in equation(3.8).Now, we claim the following fundamental identity holds based on the property of the weighted functional.Assume that ( ) u x 0 is the initial data of the mCH-Novikov equation in ( ) H s , with > s 5 2 3.36) and (3.38)-(3.41)into equation (3.32), one has

FirstN 1 .
, with the application of Sobolev embedding theorem, for < κ 0 and R 0 both depending on k 1 , k 2 , and ( ) = c i i Then, it follows from equation (A5), the exponential decay of peaked solitary wave φ ci , and bump function Ψ i , as well as the definition of localized energy function i , that

4
Third, based on the definition of the bump function Ψ i and the Abel transform, we have, equations (3.5), (3.7), and (3.10) and the exponential decay of φ ci , there holds Stability trains peaked solitary waves mCH-Novikov equations  13 Stability trains peaked solitary waves mCH-Novikov equations  27 i, then Stability trains peaked solitary waves mCH-Novikov equations  29 − III O R .