arxiv: v1 [nlin.si] 7 Oct 2013

Similar documents
Journal of Geometry and Physics

arxiv:nlin/ v1 [nlin.si] 25 Sep 2006

arxiv: v1 [nlin.si] 19 Aug 2008

arxiv: v2 [nlin.si] 3 Aug 2017

Bifurcations of Traveling Wave Solutions for a Generalized Camassa-Holm Equation

A NONLINEAR GENERALIZATION OF THE CAMASSA-HOLM EQUATION WITH PEAKON SOLUTIONS

Derivation of Generalized Camassa-Holm Equations from Boussinesq-type Equations

arxiv: v1 [nlin.si] 10 Oct 2011

The generalized Kupershmidt deformation for integrable bi-hamiltonian systems

A MULTI-COMPONENT LAX INTEGRABLE HIERARCHY WITH HAMILTONIAN STRUCTURE

A supersymmetric Sawada-Kotera equation

arxiv:nlin/ v1 [nlin.si] 23 Jul 2003

arxiv: v1 [math-ph] 17 Sep 2008

OSCILLATION-INDUCED BLOW-UP TO THE MODIFIED CAMASSA HOLM EQUATION WITH LINEAR DISPERSION

On Camassa Holm Equation with Self-Consistent Sources and Its Solutions

Supersymmetric and Deformed Harry Dym Hierarchies

Exact solutions through symmetry reductions for a new integrable equation

Symmetry reductions and travelling wave solutions for a new integrable equation

arxiv:nlin/ v1 [nlin.si] 3 Sep 2002

An integrable shallow water equation with peaked solitons

A CAUCHY-KOVALEVSKY THEOREM FOR NONLINEAR AND NONLOCAL EQUATIONS. In memory of M. Salah Baouendi

Lax Representations for Matrix Short Pulse Equations

Supersymmetric Sawada-Kotera Equation: Bäcklund-Darboux Transformations and. Applications

Reciprocal transformations and reductions for a 2+1 integrable equation

Breaking soliton equations and negative-order breaking soliton equations of typical and higher orders

Symmetry Reductions of (2+1) dimensional Equal Width. Wave Equation

arxiv:nlin/ v1 [nlin.si] 17 Jun 2006

New Integrable Decomposition of Super AKNS Equation

HAMILTONIAN STRUCTURE OF PEAKONS AS WEAK SOLUTIONS FOR THE MODIFIED CAMASSA-HOLM EQUATION

Recursion Operators of Two Supersymmetric Equations

A New Integrable Couplings of Classical-Boussinesq Hierarchy with Self-Consistent Sources

About Integrable Non-Abelian Laurent ODEs

Prolongation structure for nonlinear integrable couplings of a KdV soliton hierarchy

2. Examples of Integrable Equations

OF THE PEAKON DYNAMICS

Does the Supersymmetric Integrability Imply the Integrability of Bosonic Sector?

Bäcklund transformation and special solutions for Drinfeld Sokolov Satsuma Hirota system of coupled equations

A new integrable system: The interacting soliton of the BO

Continuous limits and integrability for a semidiscrete system Zuo-nong Zhu Department of Mathematics, Shanghai Jiao Tong University, P R China

The Cauchy problem for nonlocal evolution equations

Second Order Lax Pairs of Nonlinear Partial Differential Equations with Schwarzian Forms

Proceedings of WASCOM 2005, pages World Sci. Publishing, 2006

arxiv: v1 [math-ph] 13 Feb 2008

Evolutionary Hirota Type (2+1)-Dimensional Equations: Lax Pairs, Recursion Operators and Bi-Hamiltonian Structures

On construction of recursion operators from Lax representation

Symbolic Computation and New Soliton-Like Solutions of the 1+2D Calogero-Bogoyavlenskii-Schif Equation

Commutator representations of nonlinear evolution equations: Harry-Dym and Kaup-Newell cases

14th Panhellenic Conference of Mathematical Analysis University of Patras May 18-19, 2012

Some Elliptic Traveling Wave Solutions to the Novikov-Veselov Equation

An initial-boundary value problem for the coupled focusing-defocusing complex short pulse equation with a 4 4 Lax pair

An Improved F-Expansion Method and Its Application to Coupled Drinfel d Sokolov Wilson Equation

Keywords: Exp-function method; solitary wave solutions; modified Camassa-Holm

GLASGOW Paolo Lorenzoni

ON TRAVELING WAVE SOLUTIONS OF THE θ-equation OF DISPERSIVE TYPE

arxiv: v3 [nlin.si] 27 Jan 2009

On bosonic limits of two recent supersymmetric extensions of the Harry Dym hierarchy

Multipeakons and a theorem of Stieltjes

Integrable Rosochatius Deformations for an Integrable Couplings of CKdV Equation Hierarchy

Complete integrability of a modified vector derivative nonlinear Schrödinger equation

ANALYSIS OF A NONLINEAR SURFACE WIND WAVES MODEL VIA LIE GROUP METHOD

On recursion operators for elliptic models

Travelling Wave Solutions for the Gilson-Pickering Equation by Using the Simplified G /G-expansion Method

A PERIODICITY PROBLEM FOR THE KORTEWEG DE VRIES AND STURM LIOUVILLE EQUATIONS. THEIR CONNECTION WITH ALGEBRAIC GEOMETRY

Surfaces of Arbitrary Constant Negative Gaussian Curvature and Related Sine-Gordon Equations

A Block Matrix Loop Algebra and Bi-Integrable Couplings of the Dirac Equations

New Exact Travelling Wave Solutions for Regularized Long-wave, Phi-Four and Drinfeld-Sokolov Equations

Yong Chen a,b,c,qiwang c,d, and Biao Li c,d

Introduction to the Hirota bilinear method

Hamiltonian and quasi-hamiltonian structures associated with semi-direct sums of Lie algebras

Integration of Bi-Hamiltonian Systems by Using the Dressing Method

Multi-component bi-hamiltonian Dirac integrable equations

Coupled KdV Equations of Hirota-Satsuma Type

Global conservative solutions of the Camassa-Holm equation

On universality of critical behaviour in Hamiltonian PDEs

B.7 Lie Groups and Differential Equations

New families of non-travelling wave solutions to a new (3+1)-dimensional potential-ytsf equation

arxiv: v1 [nlin.ps] 3 Sep 2009

EXACT SOLITARY WAVE AND PERIODIC WAVE SOLUTIONS OF THE KAUP-KUPERSCHMIDT EQUATION

Compacton-like solutions in some nonlocal hydrodynamic-type models

Three types of generalized Kadomtsev Petviashvili equations arising from baroclinic potential vorticity equation

A Short historical review Our goal The hierarchy and Lax... The Hamiltonian... The Dubrovin-type... Algebro-geometric... Home Page.

Nonlinear Analysis. Variational identities and applications to Hamiltonian structures of soliton equations

The dual modified Korteweg-de Vries Fokas Qiao equation: Geometry and local analysis

Rational Form Solitary Wave Solutions and Doubly Periodic Wave Solutions to (1+1)-Dimensional Dispersive Long Wave Equation

arxiv: v1 [math.ap] 24 Oct 2014

New Formal Solutions of Davey Stewartson Equation via Combined tanh Function Method with Symmetry Method

Approximate Similarity Reduction for Perturbed Kaup Kupershmidt Equation via Lie Symmetry Method and Direct Method

Lax Representations and Zero Curvature Representations by Kronecker Product

Integrable evolution equations on spaces of tensor densities

Soliton surfaces and generalized symmetries of integrable equations

The Exact Solitary Wave Solutions for a Family of BBM Equation

On the bi-hamiltonian structures of the Camassa-Holm and Harry Dym equations

No. 11 A series of new double periodic solutions metry constraint. For more details about the results of this system, the reader can find the

arxiv: v5 [math-ph] 23 Oct 2018

HOMOTOPY ANALYSIS METHOD FOR SOLVING KDV EQUATIONS

EXISTENCE OF NONTRIVIAL SOLUTIONS FOR A QUASILINEAR SCHRÖDINGER EQUATIONS WITH SIGN-CHANGING POTENTIAL

A Method for Obtaining Darboux Transformations

arxiv:nlin/ v2 [nlin.si] 14 Sep 2001

CONVERGENCE OF SOLITARY-WAVE SOLUTIONS IN A PERTURBED BI-HAMILTONIAN DYNAMICAL SYSTEM. I. COMPACTONS AND PEAKONS.

Department of Applied Mathematics, Dalian University of Technology, Dalian , China

Transcription:

A four-component Camassa-Holm type hierarchy arxiv:1310.1781v1 [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 100083, P R China 2 Institute of Theoretical Physics University of Wroc law pl. M. Borna 9, 50-205 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 0 0 1 ϕ x λm 0 0 ϕ. (4) 1 1 0 Preprint submitted to Elsevier June 12, 2018

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 + 1 2 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 1 0 0 λn 1 0 0 ϕ, (11) 2

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

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 0 1 0 1 0 0 0 0 1 0 1 0, 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 1 2 1 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 1 1 + 1 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

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

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 3 + 3 4 ) 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 3 + 1 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 δ

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

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

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. 4.1. a three component reduction Assuming m 1 u 1 0, we have φ 1 φ 2 φ 3 x 0 0 1 λ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 0 1 0 1+λ 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 2 1 2. m 2 m 2 ( m 2 + m 2 ) +S 0 ( m 2 + m 2 ) S S +S + ( m 2 + m 2 ) 1 2 + 1 2 0 2 1 0 ( 3 4 ) 1 m 2 +2m 2x m 2 2 +3n 1 +2n 1x 3 +2x, m 2 +2m 2x m 2 2 +3n 1 +2n 1x 3 +2x 3m 2 1 m 2 m 2 2 +3m 2 1 n 1 3m 2 1 m 2 2 +3n 1 1 m 2 m 2 m 2 +3n 1 1 n 1 m 2 3n 1 1 3 1 m 2 m 2 3 1 n 1 3 1 where S m 2 (1 2 ). We remark that it is not clear how to obtain above Hamilitonian pair from (19)., 9

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 ) 1. 2 2 10

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

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

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: 10971222 and 11271366) 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) 1661-1664. 13

[2] R. Camassa, D. D. Holm and J. M. Hyman, A new integrable shallow water equation, Adv. Appl. Mech. 31 (1994) 1-33. [3] A. Degasperis and M. Procesi, Asymptotic integrability, in: Symmetry and Perturbation Theory A. Degasperis and G. Gaeta (eds.), Singapore: World Scientific (1999), pp 23-37. [4] A. Degasperis, D. D. Holm and A. N. W. Hone, A new integrable equation with peakon solutions, Theor. Math. Phys. 133 (2002) 1463-74. [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), pp93-101. [6] A. S. Fokas, On a class of physically important integrable equations, Physica D 87 (1995) 145-150. [7] B. Fuchssteiner, Some tricks from the symmetry-toolbox for nonlinear equations: generalizations of the Camassa-Holm equation, Phys. D 95 (1996) 229-243. [8] X. G. Geng and B. Xue, An extension of integrable peakon equations with cubic nonlinearity, Nonlinearity 22 (2009) 1847-1856. [9] X. G. Geng and B. Xue, A three-component generalization of Camassa-Holm equation with N-peakon solutions, Adv. Math. 226 (2011) 827-839. [10] A. N. W. Hone and J. P. Wang, Prolongation algebras and Hamiltonian operators for peakon equations, Inverse Problems 19 (2003) 129-145. [11] A. N. W. Hone and J. P. Wang, Integrable peakon equations with cubic nonlinearity, J. Phys. A: Math. Theor. 41 (2008) 372002. [12] N. H. Li and Q. P. Liu, On Bi-Hamiltonian structure of two-component Novikov equation, Phys. Lett. A 377 (2012) 257-261. [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) 126-134. [14] H. Lundmark and J. Szmigielski, Multi-peakon solutions of the Degasperis-Procesi equation, Inverse problems 19 (1994) 1241-5. [15] V. S. Novikov, Generalisations of the Camassa-Holm equation, J. Phys. A: Math. Theor. 42 (2009) 342002 (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) 1900-1906. 14

[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) 112701 (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) 013503 (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) 043513 (10 pages). [23] B. Q. Xia and Z. J. Qiao, Integrable multi-component Camassa-Holm system, arxiv: 1310.0268v1. 15