arxiv: v1 [nlin.si] 7 Oct 2013

Transcription

1 A four-component Camassa-Holm type hierarchy arxiv: v1 [nlin.si] 7 Oct 2013 Abstract Nianhua Li 1, Q. P. Liu 1, Z. Popowicz 2 1 Department of Mathematics China University of Mining and Technology Beijing , P R China 2 Institute of Theoretical Physics University of Wroc law pl. M. Borna 9, Wroc law, Poland, ziemek@ift.uni.wroc.pl We consider a 3 3 spectral problem which generates four-component CH type systems. The bi-hamiltonian structure and infinitely many conserved quantities are constructed for the associated hierarchy. Some possible reductions are also studied. Mathematical Subject Classification: 37K10, 37K05 Keywords: Bi-Hamiltonian structure, Conserved quantity, Lax pair, Integrability. 1. Introduction The Camassa-Holm (CH) equation u t u xxt +3uu x 2u x u xx +uu xxx, (1) has been a subject of steadily growing literature since it was derived from the incompressible Euler equation to model long waves in shallow water by Camassa and Holm in 1993 [1]. It is a completely integrable system, which possesses Lax representation and is a bi-hamiltonian system and especially admits peakon solutions [1][2]. CH equation is a system with quadratic nonlinearity, and another such system was proposed in 1999 by Degasperis and Procesi [3], which reads as m t +um x +3u x m 0, m u u xx, (2) now known as called DP equation whose peakon solutions were studied in [4][14]. The following third order spectral problem for the DP equation (2) was found by Degasperis, Holm and Hone [4] ψ xxx ψ x λψ (3) or equivalently it may be rewritten in the matrix form as ϕ x λm 0 0 ϕ. (4) Preprint submitted to Elsevier June 12, 2018

2 Olver and Rosenau suggested so-called tri-hamiltoninan duality approach to construct the CH type equations and many examples were worked out [17] (see also [5][6][7]). In particular, they obtained a CH type equation with cubic nonlinearity m t +[m(u 2 u 2 x )] x 0, m u u xx (5) which, sometimes is referred as Qiao equation, was studied by Qiao [19]. Recently, Song, Qu and Qiao [21] proposed a two-component generalization of (5) { mt [m(u x v x uv +uv x u x v)] x 0, m u u xx, n t [n(u x v x uv +uv x u x v)] x 0, n v v xx (6) based on the following spectral problem ( 1 λm ϕ x 2 λn 1 2 ) ϕ. (7) This system (7) is shown to be bi-hamiltonian and is associated with WKI equation [22]. Moreover Xia,Qiao and Zhou[23] generalized previous system to the integrable two-component CH type equation m t F +F x 1 2 m(uv u xv x +uv x u x v), n t G+G x n(uv u xv x +uv x u x v), (8) m u u xx, n v v xx where F and G are two arbitrary functions. Working on symmetry classification of nonlocal partial differential equations, Novikov [15] found another CH type equation with cubic nonlinearity. It reads as m t +u 2 m x +3uu x m 0, m u u xx. (9) Subsequently, Hone and Wang proposed a Lax representation for (9) and showed that it is associated to a negative flow in Sawada-Kotera hierarchy [10]. Furthermore, infinitely many conserved quantities and bi-hamiltonian structure are also presented [11]. A two-component generalization of the Novikov equation was constructed by Geng and Xue [8]. { mt +3u x vm+uvm x 0, m u u xx n t +3v x un+uvn x 0, n v v xx, (10) with a spectral problem ϕ x 0 λm λn ϕ, (11) 2

3 which reduces to the spectral Novikov s system as m n. They also calculated the N-peakons and conserved quantities and found a Hamiltonian structure. In this paper we will discuss the properties of the equations which follow from the following generalized spectral problem Φ x UΦ, U 0 λm 1 1 λn 1 0 λm 2 1 λ 0. (12) This spectral problem was proposed by one of us (ZP) recently [18] and interestingly the resulted nonlinear systems can involve an arbitrary function. This freedom allows us to recover many known CH type equations from reductions. As we will show this spectral problem takes the spectral problem considered by Geng and Xue for three-component system [9], one and two-component Novikov s equation [11][8], and one or two component Song-Qu-Qiao equation [21] as special cases. In this sense, almost all known 3 3 spectral problems for the CH type equations are contained in this case, so it is interesting to study this spectral problem. The first negative flow corresponding to the spectral problem (12) is m 1t + g 1 g 2 +m 1 (f 2 g 2 +2f 1 g 1 ) 0, where m 2t n 1 g 1 g 2 m 2 (f 1 g 1 +2f 2 g 2 ) 0, n 1t m 2 f 1 f 2 n 1 (f 2 g 2 +2f 1 g 1 ) 0, t +m 1 f 1 f 2 + (f 1 g 1 +2f 2 g 2 ) 0, m i u i u ixx,n i v i v ixx,i 1,2, f 1 u 2 v 1x, f 2 u 1 +v 2x, g 1 v 2 +u 1x, g 2 v 1 u 2x. which will be shown to be a bi-hamiltonian system. The paper is organized as follows. In the sectio, we will construct bi-hamiltonian operators related to the spectral problem (12) and present a bi-hamilotonian representation for the system (13). In the section 3, we construct infinitely many conserved quantities for the integrable hierarchy of (12). In the section 4, we consider the special reductions of our spectral problem. The last section contains concluding remarks. 2. Bi-Hamiltonian structures in general Let us consider the following Lax pair (13) Φ x UΦ, Φ t VΦ, (14) where U is defined by (12) and V (V ij ) 3 3 at the moment is an arbitrary matrix. The compatibility condition of (14) or the zero-curvature representation U t V x +[U,V] 0, (15) 3

4 is equivalent to with λm 1t V 12x V 32 +λ(m 1 V 11 + V 13 m 1 V 22 ), λm 2t V 23x +V 21 +λ(m 2 V 22 m 2 V 33 n 1 V 13 ), λn 1t V 21x +V 23 +λ(n 1 V 22 m 2 V 31 n 1 V 11 ), λt V 32x V 12 +λ( V 33 +m 1 V 31 V 22 ), V 11 V 31x +V 33 λ V 21 +λn 1 V 32, V 13 V 33x +V 31 +λm 2 V 32 λ V 23, (16) V 22x λ(n 1 V 12 +m 2 V 32 m 1 V 21 V 23 ), (17) 2V 31x +V 33xx λ(( +m 1 )V 23 ( m 2 +n 1 )V 32 m 2 V 12 + V 21 ), 2V 33x +V 31xx λ(( +m 1 )V 21 ( n 1 +m 2 )V 32 n 1 V 12 + V 23 ). Taking account of (17) and through a tedious calculation, the system (16) yields m 1 V 21 m 2 n 1 (λ 1 K+λL) V 32 V 12, (18) V 23 with and where J K t , L J +F, (19) 2m 1 1 m 1 m 1 1 m 2 J 13 J 14 m 2 1 m 1 2m 2 1 m 2 J 23 J 24 J13 J23 2n 1 1 n 1 n 1 1 J 14 J 24 1 n F (2P +S )( 3 4 ) 1 P T (2S +P )( 3 4 ) 1 S T, J 13 2m 1 1 n 1 1 m 2, J 14 m m 1,, J 23 m 2 1 n 1 +n 1 1 m 2, J 24 2m 2 1 n 1 1 m 1, P (m 1,m 2, n 1, ) T, S (,n 1, m 2,m 1 ) T. By specifying the matrix V properly, we can obtain the hierarchy of equations associated with the Lax pair (14). For example, if we assume g f 1 g 1 1 g λ 1 g 2 f V 1 1 g +f λ λ 2 1 g 1 +f 2 g 2 2 λ, f 1 f 2 f 2 λ f 2 g 2 4

5 then the resulted system, the first negative flow of the hierarchy, is exactly the equation (13). On the other hand taking the V matrix in time part of the Lax pair as 0 m 1 Γ 0 Ṽ λ n 1 Γ 0 m 2 Γ. (20) 0 Γ 0 we obtain { m1t +(Γm 1 ) x Γ 0, m 2t +(Γm 2 ) x +n 1 Γ 0, n 1t +(Γn 1 ) x +m 2 Γ 0, t +(Γ ) x m 1 Γ 0, (21) where Γ is an arbitrary function. Now the combinations of V and Ṽ leads us to the following system of equations m 1t +(Γm 1 ) x + (g 1 g 2 Γ)+m 1 (f 2 g 2 +2f 1 g 1 ) 0, m 2t +(Γm 2 ) x n 1 (g 1 g 2 Γ) m 2 (f 1 g 1 +2f 2 g 2 ) 0, n 1t +(Γn 1 ) x m 2 (f 1 f 2 Γ) n 1 (f 2 g 2 +2f 1 g 1 ) 0, t +(Γ ) x +m 1 (f 1 f 2 Γ)+ (f 1 g 1 +2f 2 g 2 ) 0, m i u i u ixx, n i v i v ixx, i 1,2, (22) which has the Lax pair Φ x UΦ, Φ t (V +Ṽ)Φ. It is obvious that the operator K given by (19) is a Hamiltonian operator. Indeed we have the following theorem Theorem 1. The operators K and L defined by (19) constitute a pair of compatible Hamiltonian operators. In particular, the four-component system (13) is a bi-hamiltonian system, namely it can be written as where H 0 m 1 m 2 n 1 t K δh 0 δm 1 δh 0 δm 2 δh 0 δn 1 δh 0 δ L (f 1 g 1 +f 2 g 2 )(m 2 f 2 +n 1 g 1 )dx, H 1 δh 1 δm 1 δh 1 δm 2 δh 1 δn 1 δh 1 δ (23) (m 2 f 2 +n 1 g 1 )dx. (24) Proof: It is easy to check that L is a skew-symmetric operator. Thus, what we need to do is to verify the Jacobi identity for L and the compatibility of two operators K and L. To this end, we follow Olver and use his multivector approach [16]. Let us introduce Θ K 1 θ Kθdx, Θ L 1 θ Lθdx (25) 2 2 5

6 with θ (θ 1,θ 2,θ 3,θ 4 ). A tedious but direct calculation shows that pr v Lθ (Θ L ) 0, pr v Kθ (Θ L )+pr v Lθ (Θ K ) 0 (26) hold. The details of this calculation are postponed to the Appendix. Since we have a bi-hamiltonian pair K and L, we may formulate an integrable hierarchy of the corresponding nonlinear evolutions equations through recursion operator. While the system (21) possesses a Lax pair, we could not expect that it is a bi-hamiltonian system for any Γ and the constraint of Γ will appear if it is requested to commute with the negative flow. To understand the appearance of the arbitrary function occurred in the present case, we now calculate the Casimir functions of the Hamiltonian operator L. Let and define L(A,B,C,D) T 0, (27) K 1 m 1 A n 1 C, K 2 m 2 B D, (28) K 3 m 2 C m 1 D, K 4 A n 1 B. (29) The system (27) consists of four equations of similar type. For example the one of them is m 1 ( 1 (2K 1 K 2 )+( 3 4 ) 1 (2(K 1 +K 2 )+(K 3 +K 4 ) x ) Solving these equations we found ( 1 K ) 1 (2(K 3 +K 4 )+(K 1 +K 2 ) x )). K 1 (m 2 Λ) x +(n 1 m 1 m 2 )Λ, K 2 (m 1 n 1 Λ) x +(n 1 m 1 m 2 )Λ, K 3 (m 1 m 2 Λ) x +(m 1 n 1 m 2 )Λ, K 4 (n 1 Λ) x +(m 1 n 1 m 2 )Λ, k where Λ m 1 n 1 +m 2 and k is an arbitrary number. Substituting above expressions for K i into (28)-(29) and solving the resulted linear equations leads to A n 1 Γ+ n 1 m 1 m 2 K m 1 K 1, B Γ+ 1 m 2 K 2, C m 1 Γ+ 1 m 2 K 3, D m 2 Γ, where Γ is an arbitrary function. This implies that L is a degenerate Hamiltonian operator. For the special case k 0 and Γ m 1 n 1 + m 2, the Casimir is H c 1 2 Γ 2 dx and hence the system (21) in this case is a Hamiltonian system m 1 m 2 n 1 K. t 6 δh c δm 1 δh c δm 2 δh c δn 1 δh c δ

7 3. Conserved quantities An integrable system normally possesses infinity number of concerned quantities and such property has been taken as one of the defining properties for integrability. In this section, we show that infinitely many conserved quantities can be constructed for the nonlinear evolution equations related with the spectral problem (12). Indeed, we may derive two sequences of conserved quantities utilizing the projective coordinates in the spectral problem. We can introduce these coordinates in three different manners as I.) a ϕ 1 ϕ 2, b ϕ 3 ϕ 2, II.) σ ϕ 2 ϕ 1, τ ϕ 3 ϕ 1, III.) α ϕ 1 ϕ 3, β ϕ 2 ϕ 3 Case 1: In these coordinates we obtain that is conserved quantity with a, b satisfy ρ (lnϕ 2 ) x λn 1 a+λm 2 b, (30) a x λm 1 +b aρ, b x a+λ bρ. (31) Substituting the Laurent series expansions in λ of a and b into (31) a i 0 a i λ i, b j 0 b j λ j, then we find a 0 0, a 1 v 2 u 1x g 1, a 2 0, b 0 0, b 1 u 1 v 2x f 2, b 2 0, and a k,x b k (n 1 a i a j +m 2 a i b j ), i+jk 1 b k,x a k (n 1 a i b j +m 2 b i b j ), (k 3). i+jk 1 With the aid of a 1,b 1, we obtain a simple conserved quantity ρ 1 (n 1 g 1 +m 2 f 2 )dx. (32) Also, due to a 3 a 3xx n 1 f 2 g 1 +m 2 f2 2 +(n 1 g1 2 +m 2 f 2 g 1 ) x, b 3 b 3xx n 1 g1 2 +m 2f 2 g 1 +(n 1 f 2 g 1 +m 2 f2 2 ) x, 7

8 we obtain the next conserved quantity by (v1 ρ 3 n 1 a 3 +m 2 b 3 dx (a 3 a 3xx)+u 2 (b 3 b ) 3xx) dx (n 1 g 1 +m 2 f 2 )(f 1 g 1 +f 2 g 2 )dx. (33) In addition, we may consider alternative expansions of a, b in negative powers of λ, namely a Σ i 0ã i λ i, b Σ j 0 b j λ j. As above, inserting these expansions into (31) we may find recursive relations for ã i, b j. The first two conserved quantities are ρ 0 m1 n 1 +m 2 dx, (34) ρ 1 Case 2: The quantity ρ defined as with σ,τ satisfy 2m1 m 2 +2n 1 +m 1 n 1x m 1x n 1 +m 2x m 2 x dx. 4(m 1 n 1 +m 2 ) ρ (lnϕ 1 ) x λm 1 σ +τ, (35) σ x λn 1 +λm 2 τ σ ρ, τ x 1+λ σ τ ρ. (36) is conserved quantity. Expanding σ and τ in Laurent series of λ then once again we may find the corresponding conserved quantities. For instance, in the case k 0, we get ρ 2 1 (m 1 + )(f 1 +g 2 )dx, 2 while in the case k 0, we obtain 2m2 2 ρ 1 +2m 1n 1 m 2x m 1 +4m 1x m 2 3x m 1 m 2 +m 1x m 1 n 1 n 1x m 2 1 dx. 4m 1 (m 1 n 1 +m 2 ) Case 3: For this case the conserved quantity is defined as with α,β satisfy ˆρ (lnϕ 3 ) x α+λ β, (37) α x λm 1 β +1 αˆρ, β x λn 1 α+λm 2 βˆρ. (38) Expanding α and β in Laurent series of λ and substituting them into (38), we may obtain the conserved quantities and apart from those found in last two cases, we have 2n1 2 ˆρ 1 +2m 1m 2 m 2x 2 +4xm 1 n 1 3m 1x n 1 +x m 2 n 1x m 1 dx. 4 (m 1 n 1 +m 2 ) Let us remark that these conserved quantities have been obtained from the x-part of the Lax pair representation only hence they are valid for the whole hierarchy. As we checked they are conserved for the system (13) as well as for the (22). 8

9 4. Reductions We now consider the possible reductions of our four component spectral problem (14) and relate them to the spectral problems known in literatures a three component reduction Assuming m 1 u 1 0, we have φ 1 φ 2 φ 3 x λn 1 0 λm 2 1 λ 0 which, by identifying the variables as follows u, m 2 v ( v u, n 1 w+ u) may be reformulated as ϕ 1 ϕ 2 ϕ 3 x λ 2 v 0 u λ 2 w 0 0 φ 1 φ 2 φ 3 x ϕ 1 ϕ 2 ϕ 3, (39), (40) This spectral problem is nothing but the one proposed by Geng and Xue [9] and the associated integrable flows are also bi-hamiltonian [13]. It is easy to see that (40) is an invertible transformation, so to find the bi-hamiltonian structures of the flows resulted from the spectral problem (39) we may convert those for the Geng-Xue spectral problem into the present case. Direct calculations yield L 1 L m 2 m 2 ( m 2 + m 2 ) +S 0 ( m 2 + m 2 ) S S +S + ( m 2 + m 2 ) ( 3 4 ) 1 m 2 +2m 2x m n 1 +2n 1x 3 +2x, m 2 +2m 2x m n 1 +2n 1x 3 +2x 3m 2 1 m 2 m m 2 1 n 1 3m 2 1 m n 1 1 m 2 m 2 m 2 +3n 1 1 n 1 m 2 3n m 2 m n where S m 2 (1 2 ). We remark that it is not clear how to obtain above Hamilitonian pair from (19)., 9

10 4.2. a two-component reduction or In this case, we assume φ 1 φ 2 φ 3 x n 1 m 2, m 1 0 λm 1 1 λm 2 0 λm 2 1 λm 1 0 which yields ( ) ( φ1 +φ 3 1 2λm1 φ 2 λm x 2 0 By the following change of variables (41) gives ( ψ1 ψ 2 ) φ 1 φ 2 φ 3 )( φ1 +φ 3 φ 2, φ 1 +φ 3 e 1 2 x ψ 1, φ 2 e 1 2 x ψ 2, 2m 1 m, m 2 n, x ( 1 λm 2 λn 1 2 )( φ1 ψ 2 ), ). (41) a spectral problem considered by Song, Qu and Qiao [21]. Therefore the bi-hamiltonian structure of Song-Qu-Qiao system (see [22]) is hidden in our bi-hamiltonian structure. 5. Concluding remarks In this paper, started from a general 3 3 problem, we considered the related four-component CH type systems. We obtained the bi-hamiltonian structure and suggested the way to construction of infinitely many conserved quantities for the integrable equations. Different reductions were also considered. As noticed above, the positive flows allow for an arbitrary function Γ involved and such systems are interesting since different specifications of Γ lead to different CH type equations. Although the flow equations with arbitrary Γ do possess infinitely many conserved quantities, we do not expect they are (bi-) Hamiltonian systems in general case. Also, we explained the appearance of this arbitrary function by studying the kernel of one of the Hamiltonian operators and it seems that further study of such systems is needed. A remarkable property of CH type equations is that it possess peakon solutions. One may find that the first negative flow, which does not depend on Γ, only possess stationary peakons. The systems such as (22) may admit non-stationary peakon solutions. This and other related issues may be considered in further publication. Appendix A. We first prove that L is also a Hamiltonian operator. For this purpose, we first define Θ J 1 θ Jθdx, Θ F 1 θ Fθdx, A ( 3 4 )

11 By direct calculation, we have ((n2 Θ J θ 1 n 1 θ 2 ) 1 (m 1 θ 4 m 2 θ 3 )+(m 1 θ 1 n 1 θ 3 ) 1 ( θ 4 m 2 θ 2 ) + ( θ 4 m 2 θ 2 ) 1 ( θ 4 m 2 θ 2 )+(m 1 θ 1 n 1 θ 3 ) 1 (m 1 θ 1 n 1 θ 3 ) ) dx, and Θ F 1 ((2Q R ) AQ+(Q ) 2R) AR dx 2 (Q AQ R AR+Q AR ) dx, where Q m 1 θ 1 +m 2 θ 2 n 1 θ 3 θ 4, R θ 1 n 1 θ 2 +m 2 θ 3 m 1 θ 4. Since prv Lθ (Θ L ) prv Lθ (Θ J )+prv Lθ (Θ F ), (A1) we calculate prv Lθ (Θ J ) and pr v Lθ (Θ F ). Indeed, we have [( prv Lθ (Θ J ) 1 ( θ 1 n 1 θ 2 ) 1 (m 1 θ 4 m 2 θ 3 ) 1 Q) x +(2m 1 θ 1 +2n 1 θ 3 θ 4 m 2 θ 2 ) (AR x +2AQ) 1 (m 1 θ 1 n 1 θ 3 ) +(2m 2 θ 3 2 θ 1 +m 1 θ 4 n 1 θ 2 ) (AQ x +2AR) 1 (m 1 θ 1 n 1 θ 3 ) +(2m 1 θ 4 2n 1 θ 2 +m 2 θ 3 θ 1 ) (AQ x +2AR) 1 ( θ 4 m 2 θ 2 ) +(m 1 θ 1 +n 1 θ 3 2m 2 θ 2 2 θ 4 ) (AR x +2AQ) 1 ( θ 4 m 2 θ 2 ) +(m 1 θ 4 m 2 θ 3 ) (AR x +2AQ) 1 ( θ 1 n 1 θ 2 ) (n 1 θ 3 + θ 4 ) (AQ x +2AR) 1 ( θ 1 n 1 θ 2 ) +(m 1 θ 1 +m 2 θ 2 ) (AQ x +2AR) 1 (m 1 θ 4 m 2 θ 3 ) ( θ 1 n 1 θ 2 ) (AR x +2AQ) 1 (m 1 θ 4 m 2 θ 3 ) ] dx. Next we consider pr v Lθ (Θ F ). For simplicity, we denote Υ 2m 1 θ 1 m 2 θ 2 2n 1 θ 3 + θ 4, Ω m 1 θ 1 2m 2 θ 2 n 1 θ 3 +2 θ 4, 11

12 then a direct calculation shows that prv Lθ (Θ F ) can be expressed as prv Lθ (Θ F ) [(m1 θ 1 +n 1 θ 3 ) 1 Υ (2AQ+AR x ) +2( θ 1 n 1 θ 2 m 2 θ 3 +m 1 θ 4 ) (AR x +2AQ) (AQ x +2AR) (m 2 θ 2 + θ 4 ) 1 Ω (2AQ+AR x ) +( θ 1 n 1 θ 2 ) 1 (m 1 θ 4 m 2 θ 3 ) (2AQ+AR x ) (m 1 θ 4 m 2 θ 3 ) 1 ( θ 1 n 1 θ 2 ) (2AQ+AR x ) +(m 1 θ 4 n 1 θ 2 ) 1 Υ (2AR+AQ x ) ( θ 1 m 2 θ 3 ) 1 Ω (2AR+AQ x ) +( θ 4 +n 1 θ 3 ) 1 (m 1 θ 4 m 2 θ 3 ) (2AR+AQ x ) (m 2 θ 2 +m 1 θ 1 ) 1 ( θ 1 n 1 θ 2 ) (2AR+AQ x ) ] dx. Letting f θ 1 n 1 θ 2,g m 1 θ 4 m 2 θ 3 and substituting above expansions into (??) lead to prv Lθ (Θ L ) 2( θ 1 n 1 θ 2 m 2 θ 3 +m 1 θ 4 ) (AR x +2AQ) (AQ x +2AR) +2(( θ 1 n 1 θ 2 ) 1 (m 1 θ 4 m 2 θ 3 ) +(m 2 θ 3 m 1 θ 4 ) 1 ( θ 1 n 1 θ 2 )) (2AQ+AR x ) +(( θ 1 n 1 θ 2 m 2 θ 3 +m 1 θ 4 ) 1 Q + 1 ( θ 1 n 1 θ 2 m 2 θ 3 +m 1 θ 4 ) Q) (2AR+AQ x )dx (f 1 2g + 1 f 2g) (2AQ+AR x ) +2(f +g) (AR x +2AQ) (AQ x +2AR) +((f +g) 1 Q+ 1 (f +g) Q) (2AR+AQ x )dx 1 (2g +R) (Q (AQ x +2AR) R (AR x +2AQ)) +(2g +R) ( 1 R (AR x +2AQ)+ 1 Q (AQ x +2AR) +2(AR x +2AQ)(AQ x +2AR))dx 1 (2g +R) ( 1 R (AR xx +2AQ x ) 1 Q (AQ xx +2AR x ) 2(AR xx +2AQ x )(AQ x +2AR) 2(AR x +2AQ)(AQ xx +2AR x ))dx 1 (2g +R)( 1 R (2AQ x +4AR) 1 Q (2AR x +4AQ) 0, 2 1 R (AQ x +2AR) 2(AR x +2AQ) 1 Q)dx. where we use f g R for short. Thus, L given by (??) is a Hamiltonian operator. Finally we prove the compatibility of K and L, which is equivalent to prv Lθ (Θ K )+prv Kθ (Θ L ) prv Kθ (Θ L ) prv Kθ (Θ J )+prv Kθ (Θ F ) 0. (A2) 12

13 To this end, we notice prv Kθ (Θ J ) (2(θ 4 θ 2 ) x +(θ 1 θ 3 ) x +θ 1 θ 2 +θ 3 θ 4 ) 1 ( θ 4 m 2 θ 2 ) +((θ 4 θ 2 ) x +2(θ 1 θ 3 ) x θ 1 θ 2 θ 3 θ 4 ) 1 (m 1 θ 1 n 1 θ 3 ) +((θ 1 θ 2 ) x θ 1 θ 3 θ 2 θ 4 ) 1 (m 1 θ 4 m 2 θ 3 ) +((θ 4 θ 3 ) x +θ 1 θ 3 +θ 2 θ 4 ) 1 ( θ 1 n 1 θ 2 )dx ( (θ 1 θ 2 +θ 3 θ 4 ) 1 Q+(θ 1 θ 3 +θ 2 θ 4 ) 1 R (2θ 4 θ 2 +θ 1 θ 3 ) ( θ 4 m 2 θ 2 ) (θ 1 θ 2 ) (m 1 θ 4 m 2 θ 3 ) (θ 4 θ 2 +2θ 1 θ 3 ) (m 1 θ 1 n 1 θ 3 )+(θ 3 θ 4 ) ( θ 1 n 1 θ 2 )dx (θ 1 θ 2 +θ 3 θ 4 ) 1 Q+(θ 1 θ 3 +θ 2 θ 4 ) 1 Rdx. and prv Kθ (Θ F ) (θ 1 (θ 3x θ 2 )+θ 2 (θ 1 +θ 4x ) θ 3 (θ 1x +θ 4 ) θ 4 (θ 2x θ 3 )) (2AQ+ AR)+(AQ x +2AR) ((θ 2x θ 3 ) θ 1 (θ 1x +θ 4 ) θ 2 +(θ 1 +θ 4x ) θ 3 (θ 3x θ 2 ) θ 4 )dx. (( 2 4)(θ 1 θ 2 +θ 3 θ 4 )) AQ (( 2 4)(θ 1 θ 3 +θ 2 θ 4 )) ARdx (θ 1 θ 2 +θ 3 θ 4 ) 1 Q (θ 1 θ 3 +θ 2 θ 4 ) 1 Rdx. Therefore we arrive at prv Kθ (Θ L ) prv Kθ (Θ J +Θ F ) 0, so K and L are two compatible Hamiltonian operators. Acknowledgments We would like to thank Decio Levi for interesting discussions. NHL and QPL are supported by the National Natural Science Foundation of China (grant numbers: and ) and the Fundamental Research Funds for Central Universities. ZP thanks the China University of Mining and Technology for its support during his visit to Beijing where this work has been carried out. References [1] R. Camassa and D. D. Holm, An integrable shallow water equation with peaked solitons, Phys. Rev. Lett. 71 (1993)

14 [2] R. Camassa, D. D. Holm and J. M. Hyman, A new integrable shallow water equation, Adv. Appl. Mech. 31 (1994) [3] A. Degasperis and M. Procesi, Asymptotic integrability, in: Symmetry and Perturbation Theory A. Degasperis and G. Gaeta (eds.), Singapore: World Scientific (1999), pp [4] A. Degasperis, D. D. Holm and A. N. W. Hone, A new integrable equation with peakon solutions, Theor. Math. Phys. 133 (2002) [5] A. S. Fokas, P. J. Olver and P. Rosenau, A plethora of integrable bi-hamiltonian equations, in: Algebraic Aspects of Integrable Systems: In Memory of Irene Dorfman, A. S. Fokas and I. M. Gelfand (eds.), Birkhauser: Boston (1996), pp [6] A. S. Fokas, On a class of physically important integrable equations, Physica D 87 (1995) [7] B. Fuchssteiner, Some tricks from the symmetry-toolbox for nonlinear equations: generalizations of the Camassa-Holm equation, Phys. D 95 (1996) [8] X. G. Geng and B. Xue, An extension of integrable peakon equations with cubic nonlinearity, Nonlinearity 22 (2009) [9] X. G. Geng and B. Xue, A three-component generalization of Camassa-Holm equation with N-peakon solutions, Adv. Math. 226 (2011) [10] A. N. W. Hone and J. P. Wang, Prolongation algebras and Hamiltonian operators for peakon equations, Inverse Problems 19 (2003) [11] A. N. W. Hone and J. P. Wang, Integrable peakon equations with cubic nonlinearity, J. Phys. A: Math. Theor. 41 (2008) [12] N. H. Li and Q. P. Liu, On Bi-Hamiltonian structure of two-component Novikov equation, Phys. Lett. A 377 (2012) [13] N. H. Li and Q. P. Liu, Bi-Hamiltonian structure of a three-component Camassa-Holm type equation, J. Nonlinear Math. Phys. 20 (2013) [14] H. Lundmark and J. Szmigielski, Multi-peakon solutions of the Degasperis-Procesi equation, Inverse problems 19 (1994) [15] V. S. Novikov, Generalisations of the Camassa-Holm equation, J. Phys. A: Math. Theor. 42 (2009) (14 pages). [16] P. J. Olver, Applications of Lie Groups to Differential Equations, 2nd edn, Berlin: Springer (1993). [17] P. J. Olver and P. Rosenau, Tri-Hamiltonian duality between solitons and solitary-waves solutions having compact support, Phys. Rev. E 53 (1996)

15 [18] Z. Popowicz, The generalizations of the peakon s systems, talk given at the third international conference on Nonlinear Waves: theory and applications, June 12-15, 2013, Beijing. [19] Z. J. Qiao, A new integrable equation with cuspons and W/M-shape-peaks solitons, J. Math. Phys. 47 (2006) (9 pages). [20] C. Z. Qu, J. F. Song and R. X. Yao, Multi-component integrable systems and invariant curve flows in certian geometries, SIGMA 9 (2013) 001 (19 pages). [21] J. F. Song, C. Z. Qu and Z. J. Qiao, A new integrable two-component system with cubic nonlinearity, J. Math. Phys. 52 (2011) (9 pages). [22] K. Tian and Q. P. Liu, Tri-Hamiltonian duality between the Wadati-Konno-Ichikawa hierarchy and Song-Qu-Qiao hierarchy, J. Math. Phys. 54 (2013) (10 pages). [23] B. Q. Xia and Z. J. Qiao, Integrable multi-component Camassa-Holm system, arxiv: v1. 15