A supersymmetric Sawada-Kotera equation

Similar documents
Recursion Operators of Two Supersymmetric Equations

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

arxiv: v1 [nlin.si] 10 Oct 2011

Nonlocal conservation laws for supersymmetric KdV equations

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

New Integrable Decomposition of Super AKNS Equation

Supersymmetric and Deformed Harry Dym Hierarchies

On construction of recursion operators from Lax representation

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

arxiv: v1 [nlin.si] 23 Aug 2007

INTEGRABILITY AND BIHAMILTONIAN STRUCTURE OF THE EVEN ORDER SKdV HIERARCHIES

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

Solutions of the nonlocal nonlinear Schrödinger hierarchy via reduction

arxiv: v1 [nlin.si] 7 Oct 2013

SYMBOLIC SOFTWARE FOR SOLITON THEORY: INTEGRABILITY, SYMMETRIES CONSERVATION LAWS AND EXACT SOLUTIONS. Willy Hereman

Boussineq-Type Equations and Switching Solitons

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

Introduction to the Hirota bilinear method

New approach for tanh and extended-tanh methods with applications on Hirota-Satsuma equations

Solving the Generalized Kaup Kupershmidt Equation

Generalized bilinear differential equations

Multisolitonic solutions from a Bäcklund transformation for a parametric coupled Korteweg-de Vries system

-Expansion Method For Generalized Fifth Order KdV Equation with Time-Dependent Coefficients

Crum Transformation and Wronskian Type Solutions for Supersymmetric KdV Equation

Exact Solutions of Supersymmetric KdV-a System via Bosonization Approach

Generalized Burgers equations and Miura Map in nonabelian ring. nonabelian rings as integrable systems.

A General Formula of Flow Equations for Harry Dym Hierarchy

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

Universidad del Valle. Equations of Lax type with several brackets. Received: April 30, 2015 Accepted: December 23, 2015

arxiv: v1 [nlin.si] 19 Dec 2017

arxiv: v1 [math-ph] 10 Sep 2015

EXACT SOLUTION TO TIME FRACTIONAL FIFTH-ORDER KORTEWEG-DE VRIES EQUATION BY USING (G /G)-EXPANSION METHOD. A. Neamaty, B. Agheli, R.

Building Generalized Lax Integrable KdV and MKdV Equations with Spatiotemporally Varying Coefficients

arxiv:hep-th/ v2 21 May 1993

The Generalized Harry Dym Equation.

A new integrable system: The interacting soliton of the BO

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

2. Examples of Integrable Equations

Coupled KdV Equations of Hirota-Satsuma Type

The generalized Kupershmidt deformation for integrable bi-hamiltonian systems

REFi INTERNATIONAL CENTRE FOR THEORETICAL PHYSICS IC/92/70 THE BI-HAMILTONIAN STRUCTURES OF THE MANIN-RADUL SUPER KP HIERARCHY. Sudhakar Panda.

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

THE LAX PAIR FOR THE MKDV HIERARCHY. Peter A. Clarkson, Nalini Joshi & Marta Mazzocco

Integration of Bi-Hamiltonian Systems by Using the Dressing Method

Journal of Nonlinear Mathematical Physics, Vol. 19, No. 1 (2012) (11 pages) c S. Carillo and C. Schiebold DOI: /S

Painlevé Test for the Certain (2+1)-Dimensional Nonlinear Evolution Equations. Abstract

A new integrable equation valued on a Cayley-Dickson algebra

arxiv:patt-sol/ v1 25 Sep 1995

Research Article Application of the Cole-Hopf Transformation for Finding Exact Solutions to Several Forms of the Seventh-Order KdV Equation

arxiv: v2 [nlin.si] 27 Jun 2015

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

University of Warwick institutional repository: This paper is made available online in accordance with publisher

Gauge transformations of constrained KP ows: new integrable hierarchies. Anjan Kundu and Walter Strampp. GH{Universitat Kassel. Hollandische Str.

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

Compatible Hamiltonian Operators for the Krichever-Novikov Equation

Wronskians, Generalized Wronskians and Solutions to the Korteweg-de Vries Equation

arxiv: v1 [math-ph] 25 Jul Preliminaries

Two-Componet Coupled KdV Equations and its Connection with the Generalized Harry Dym Equations

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

The supersymmetric Camassa-Holm equation and geodesic flow on the superconformal group

New explicit solitary wave solutions for (2 + 1)-dimensional Boussinesq equation and (3 + 1)-dimensional KP equation

RECURSION OPERATORS AND CLASSIFICATION OF INTEGRABLE NONLINEAR EQUATIONS

Backlund transformation for integrable hierarchies: example - mkdv hierarchy

Freedom in the Expansion and Obstacles to Integrability in Multiple-Soliton Solutions of the Perturbed KdV Equation

Lax Representations and Zero Curvature Representations by Kronecker Product

Prolongation structure for nonlinear integrable couplings of a KdV soliton hierarchy

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

MACSYMA PROGRAM FOR THE PAINLEVÉ TEST FOR NONLINEAR ORDINARY AND PARTIAL DIFFERENTIAL EQUATIONS

HOMOTOPY ANALYSIS METHOD FOR SOLVING KDV EQUATIONS

arxiv:math-ph/ v2 21 Feb 2006

N=2 supersymmetric unconstrained matrix GNLS hierarchies are consistent

arxiv: v1 [math-ph] 4 May 2016

A MULTI-COMPONENT LAX INTEGRABLE HIERARCHY WITH HAMILTONIAN STRUCTURE

Rogue periodic waves for mkdv and NLS equations

Bäcklund transformation and soliton solutions in terms of the Wronskian for the Kadomtsev Petviashvili-based system in fluid dynamics

Exact Solutions for the Nonlinear (2+1)-Dimensional Davey-Stewartson Equation Using the Generalized ( G. )-Expansion Method

Recursion Systems and Recursion Operators for the Soliton Equations Related to Rational Linear Problem with Reductions

Lump-type solutions to nonlinear differential equations derived from generalized bilinear equations

Traveling Wave Solutions For The Fifth-Order Kdv Equation And The BBM Equation By ( G G

Integrable semi-discretizations and full-discretization of the two dimensional Leznov lattice

A NEW VARIABLE-COEFFICIENT BERNOULLI EQUATION-BASED SUB-EQUATION METHOD FOR SOLVING NONLINEAR DIFFERENTIAL EQUATIONS

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

University of Twente. Faculty of Mathematical Sciences. Deformation and recursion for the N =2α =1 supersymmetric KdV hierarchy. Memorandum No.

KdV Equation with Self-consistent Sources in Non-uniform Media

The Construction of Alternative Modified KdV Equation in (2 + 1) Dimensions

Traveling Wave Solutions For Two Non-linear Equations By ( G G. )-expansion method

arxiv:nlin/ v1 [nlin.ps] 12 Jul 2001

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

Soliton and Periodic Solutions to the Generalized Hirota-Satsuma Coupled System Using Trigonometric and Hyperbolic Function Methods.

APPLICATIONS OF MAPLE SOFTWARE TO DERIVE EXACT SOLUTIONS OF GENERALIZED FIFTH-ORDER KORTEWEG-DE VRIES EQUATION WITH TIME-DEPENDENT COEFFICIENTS

arxiv:nlin/ v1 [nlin.si] 7 Sep 2005

Recursion operators of some equations of hydrodynamic type

Elsayed M. E. Zayed 1 + (Received April 4, 2012, accepted December 2, 2012)

Integrable Equations and Motions of Plane Curves

The quasiclassical limit of the symmetry constraint of the KP hierarchy and the dispersionless KP hierarchy with self-consistent sources

Exact Solutions of Matrix Generalizations of Some Integrable Systems

Traveling Wave Solutions For Three Non-linear Equations By ( G G. )-expansion method

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

Algorithmic Computation of Symmetries, Invariants and Recursion Operators for Systems of. Nonlinear Evolution and Differential-Difference.

Transcription:

A supersymmetric Sawada-Kotera equation arxiv:0802.4011v2 [nlin.si] 7 Dec 2008 Kai Tian and Q. P. Liu Department of Mathematics, China University of Mining and Technology, Beijing 100083, P.R. China Abstract A new supersymmetric equation is proposed for the Sawada-Kotera equation. The integrability of this equation is shown by the existence of Lax representation and infinite conserved quantities and a recursion operator. PACS: 02.30.Ik; 05.45.Yv Key Words: integrability; Lax representation; recursion operator; supersymmetry 1 Introduction The following fifth-order evolution equation u t +u xxxxx +5uu xxx +5u x u xx +5u 2 u x = 0 (1) is a well-known system in soliton theory. It was proposed by Sawada and Kotera, also by Caudrey, Dodd and Gibbon independently, more than thirty years ago [1][2], so it is referred as Sawada-Kotera(SK) equation or Caudrey-Dodd-Gibbon-Sawada-Kotera equation in literature. Now there are a large number of papers about it and thus its various properties are established. For example, its Bäcklund transformation and Lax representation were given in [3][4], its bi-hamiltonian structure was worked out by Fuchssteiner and Oevel [5], and a Darboux transformation was derived for this system [6][7], to mention just a few (see also [8][9]). Soliton equations or integrable systems have supersymmetric analogues. Indeed, many equations such as KdV, KP, and NLS equations were embedded into their supersymmetric counterparts and it turns out that these supersymmetric systems have also remarkable properties. Thus, it is interesting to work out supersymmetric extensions for a given integrable equation. Email: qpl@cumtb.edu.cn Tel: 86 10 62339015 1

The aim of the Note is to propose a supersymmetric extension for the SK equation. In this regard, we notice that Carstea [10], based on Hirota bilinear approach, presented the following equation φ t +φ xxxxx + [ 10(Dφ)φ xx +5(Dφ) xx φ+15(dφ) 2 φ ] x = 0 where φ = φ(x,t,θ) is a fermioic super variable depending on usual temporal variable t and super spatial variables x and θ. D = θ + θ x is the super derivative. Rewriting the equation in components, it is easy to see that this system does reduce to the SK equation when the fermionic variable is absent. However, apart from the fact that the system can be put into a Hirota s bilinear form, not much is known for its integrability. We will give an alternative supersymmetric extension for the SK equation and will show the evidence for the integrability of our system. The paper is organized as follows. In section two, by considering a Lax operator and its factorization, we construct the supersymmetric SK (ssk) equation. In section three, we will show that our ssk equation has an interesting property, namely, it does not have the usual bosonic conserved quantities since those, resulted from the super residues of a fractional power for Lax operator, are trivial. Evermore, there are infinite fermionic conserved quantities. In the section four, we construct a recursion operator for our ssk equation. Last section contains a brief summary of our new findings and presents some interesting open problems. 2 Supersymmetric Sawada-Kotera Equation The main purpose of this section is to construct a supersymmetric analogy for the SK equation. To this end, we will work with the algebra of super-pseudodifferential operators on a (1 1) superspace with coordinates (x,θ). We start with the following general Lax operator L = 3 x +Ψ xd +U x +ΦD +V. (2) By the standard fractional power method [12], we have an integrable hierarchy of equations given by L t n = [(L n 3 )+,L] (3) where we are using the standard notations: [A,B] = AB ( 1) A B BA is the supercommutator and the subscript + means taking the projection to the differential part for a given super-pseudodifferential operator. It is remarked that the system (3) is a kind of even order generalized SKdV hierarchies considered in [11]. In the following, we will consider the particular t 5 flow. Our interest here is to find a minimal supersymmetric extension for the SK equation, so we have to do reductions for the general Lax operator (2). To this end, we impose L+L = 0 2

where means taking formal adjoint. Then we find Ψ = 0, V = 1 2 (U x (DΦ)) that is L = 3 x +U x +ΦD + 1 2 (U x (DΦ)) a Lax operator with two field variables. In this case, we take B = 9(L 5 3) +, namely B = 9 5 x +15U 3 x +15ΦD 2 x +15(U x +V) 2 x +15Φ x D x +(10U xx +15V x +5U 2 ) x +10(Φ xx +ΦU)D +10V xx +10UV +5Φ(DU) for convenience. Then, the flow of equations, resulted from L t = [B,L] reads as U t +U xxxxx +5 (UU xx + 34 U2x + 13 U3 +Φ x (DU)+ 12 Φ(DU x)+ 12 ΦΦ x 34 ) (DΦ)2 ( Φ t +Φ xxxxx +5 UΦ xx + 1 2 U xxφ+ 1 2 U xφ x +U 2 Φ+ 1 2 Φ(DΦ x) 1 ) 2 (DΦ)Φ x where we identify t 5 with t for simplicity. Remarks: x x = 0(4a) = 0(4b) 1. It is interesting to note that the above system has an obvious reduction. Indeed, setting Φ = 0, we will have the standard Kaup-Kupershimdt (KK) equation. Therefore, we may consider it as a supersymmetric extension of the KK equation. 2. The coupled system (4a-4b) admits the following simple Hamiltonian structure ( ) ( ) Ut 0 x = δh x 0 Φ t where the Hamiltonian is given by H = [ 5 4 Φ(Dφ)2 (DU x )(DΦ xx ) 5 3 ΦU3 5 4 (DU x)(du)φ + 5 4 (DU)U x(dφ)+5(du)u(dφ x )+ 5 ] 2 (DU)Φ xφ dxdθ. 3

At this point, it is not clear how this system (4a-4b) is related to the SK equation. To find a supersymmetric SK equation from it, we now consider the factorization of the Lax operator in the following way which gives us a Miura-type transformation L = 3 x +U x +ΦD + 1 2 (U x (DΦ)) = (D 3 +WD+Υ)(D 3 DW +Υ), (5) U = 2W x W 2 +(DΥ), Φ = Υ x 2ΥW, and the modified system corresponding to this factorization is given by W t +W xxxxx +5W xxx (DΥ) 5W xxx W x 5W xxx W 2 5W 2 xx +10W xx(dυ x ) 20W xx W x W 5W xx W(DΥ) 5W 3 x 5W 2 x(dυ)+5w x (DΥ xx )+5W x (DΥ) 2 +5W x W 4 5W x W(DΥ x )+10W(DΥ x )(DΥ) 10Υ x ΥW x +5(DW xx )Υ x +5(DW x )Υ xx +5(DW x )ΥW x +10(DW x )ΥW 2 5(DW)Υ x W x +10(DW)Υ x W 2 +10(DW)Υ(DΥ x ) 5(DW)ΥW xx +30(DW)ΥW x W = 0, Υ t +Υ xxxxx +5Υ xxx (DΥ) 5Υ xxx W 2 +5Υ xx (DΥ x )+5Υ xx W xx 25Υ xx W x W +5Υ xx W(DΥ)+5Υ x (DΥ) 2 +5Υ x W xxx 25Υ x W xx W 25Υ x W 2 x +5Υ xw x (DΥ) +10Υ x W x W 2 +5Υ x W 4 10Υ x W 2 (DΥ)+5Υ x W(DΥ x ) 10ΥW xxx W 20ΥW xx W x +10ΥW xx W 2 +30ΥW 2 xw +20ΥW x W 3 30ΥW x W(DΥ) 10ΥW 2 (DΥ x ) 5(DW x )Υ x Υ 5(DW)Υ xx Υ 10(DW)Υ x ΥW = 0. Although this modification does indeed have a complicated form, the remarkable fact is that it allows a simple reduction. What we need to do is simply putting W to zero, namely In this case, we have W = 0, Υ = φ. φ t +φ xxxxx +5φ xxx (Dφ)+5φ xx (Dφ x )+5φ x (Dφ) 2 = 0 (6) this equation is our supersymmetric SK equation. To see the connection with the original SK equation (1), we let φ = θu(x,t)+ξ(x,t) and write the equation (6) out in components u t +u xxxxx +5uu xxx +5u x u xx +5u 2 u x 5ξ xxx ξ x = 0, ξ t +ξ xxxxx +5uξ xxx +5u x ξ xx +5u 2 ξ x = 0. (7a) (7b) It is now obvious that the system reduces to the SK equation when ξ = 0. Therefore, our system (6) does qualify as a supersymmetric SK equation. 4

Our system (6) is integrable in the sense that it has a Lax representation. In fact, the factorization (5) implies that the reduced Lax operator has the following appealing form L = (D 3 +φ)(d 3 +φ) (8) or L = 3 x +(Dφ) x φ x D +(Dφ x ) 3 Conserved Quantities In general, an integrable system has infinite number of conserved quantities. Since the ssk equation hasasimple Laxoperator(8), it isnaturalto takeadvantageofthefractionalpower method of Gel fand and Dickey [12] to find conserved quantities. In the present situation, we have to work with the super residue of a pesudodifferential operator. The obvious choice in this case is to consider the operators L n 3 and their super residues. Then, we have the Proposition 1 sresl n 3 ImD. where sres means taking the super residue of a super pseudodifferential operator. Proof: As observed already in [13], there exists a unique odd operator Λ = D+O(1), whose coefficients are all differential polynomials of φ, such that thus, the Lax operator (8) is written as From it we have (D 3 +φ) = Λ 3, L = (D 3 +φ)(d 3 +φ) = Λ 6. sresl n 3 = sresλ 2n = 1 2 sres{λ2n 1 Λ+ΛΛ 2n 1 } = 1 2 sres[λ2k 1,Λ] ImD. This completes the proof. Remark: The triviality of L n 3 implies that the Lax operator could not generate any Hamiltonian structures for the equation (6). To find nontrivial conserved quantities, we now turn to L n 6 rather than L n 3. It is easy to prove that 5 t L1 6 = [9(L 3 )+,L 1 6 ] thus t Ln 6 = [9(L 5 3 )+,L n 6 ]. 5

Consequently, the super residue of L n 6 is conserved. By direct calculation, we obtain the first two nontrivial conserved quantities sresl 7 1 6 dxdθ = [2(Dφ xx )+(Dφ) 2 6φ x φ]dxdθ 9 sresl 11 1 6 dxdθ = [6(Dφ xxxx )+18(Dφ xx )(Dφ)+9(Dφ x ) 2 81 +4(Dφ) 3 18φ xxx φ+6φ xx φ x 36φ x φ(dφ)]dxdθ Remarks: 1. What is remarkable is that the conserved quantities found in this way, unlike the supersymmetric KdV case [15][14], are local. 2. All those conserved quantities are fermioic. To our knowledge, this is the first supersymmetric integrable system whose only conserved quantities are fermioic. 4 Recursion operator An integrable system often appears as a particular flow of hierarchy equations and an important ingredient in this aspect is the existence of recursion operators. In this section, we deduce the recursion operator for the ssk equation (6) following the method proposed in [16]. We first notice that the ssk hierarchy can be written as t n L = [(L n 3 )+,L] (9) where L is given by (8). It is easy to see that the flow equations are nontrivial only if n is an integer satisfying n 0 mod 3 and n = 1 mod 2. Therefore the next flow which is achieved by applying recursion operator to (9) should be t n+6 L = [(L n+6 3 )+,L]. (10) But [ ] [ (L ] (L n+6 3 )+,L = 2 (L n 3 )+ +L 2 (L n 3 ) )+,L = [ L 2 (L n 3 )+,L ] [ (L ] + 2 (L n 3 ) )+,L = L 2[ (L n 3 )+,L ] +[R n,l] = L 2 t n L+[R n,l] 6

where R n = ( L 2 (L n 3 ) )+ (11) is a differential operator of O( x 5 D), that is, Therefore, R n = (α 5 x +β 4 x +γ 3 x +δ 2 x +ξ x +η)d +a 5 x +b 4 x +c 3 x +d 2 x +e x +f. t n+6 L = L 2 t n L+[R n,l]. (12) Next we may determine the coefficients in R n. Using (12), we obtain a = 1 3 (D 1 φ n ), b = 2(Dφ n ), c = 44 9 (Dφ n,x)+ 5 3 (Dφ)(D 1 φ n )+ 4 9 ( 1 x φ x φ n ), d = 55 9 (Dφ n,xx)+ 19 9 (Dφ)(Dφ n)+ 5 9 φ xφ n + 10 9 (Dφ x)(d 1 φ n ) e = 1 27 {106(Dφ n,xxx)+74(dφ)(dφ n,x ), 14φ x φ n,x +79(Dφ x )(Dφ n ) +27φ xx φ n +[23(Dφ xx )+4(Dφ) 2 ](D 1 φ n )+16(Dφ)( 1 x φ xφ n ) +2D 1 [(φ xxx +φ x (Dφ))(D 1 φ n ) 3(Dφ)(D 1 φ x φ n ) 2φ x ( 1 x φ x φ n )+2D 1 (φ xxx φ n +2φ x (Dφ)φ n )]}, f = 1 27 {28(Dφ n,xxxx)+32(dφ)(dφ n,xx ) 20φ x φ n,xx +54(Dφ x )(Dφ n,x ) +16φ xx φ n,x +[30(Dφ xx )+4(Dφ) 2 ](Dφ n )+[8φ xxx +4φ x (Dφ)]φ n +[10(Dφ xxx )+10(Dφ x )(Dφ)](D 1 φ n ) 8φ x (D 1 φ x φ n ) +12(Dφ x )( 1 x φ x φ n )}, α = 0, β = 1 3 φ n, γ = 5 3 φ n,x, δ = 1 9 {29φ n,xx+5φ n (Dφ)+5φ x (D 1 φ n ) 2(D 1 φ x φ n )}, ξ = 1 9 {26φ n,xxx+16φ x (Dφ n )+3φ n (Dφ x )+14φ n,x (Dφ)+5φ xx (D 1 φ n )}, η = 1 27 {28φ n,xxxx+32(dφ)φ n,xx +28φ x (Dφ n,x )+26(Dφ x )φ n,x +28φ xx (Dφ n )+[2(Dφ xx )+4(Dφ) 2 ]φ n +[10φ xxx +10φ x (Dφ)](D 1 φ n ) 2(Dφ)(D 1 φ x φ n )+12φ x ( 1 x φ xφ n ) 2D 1 [φ xxx φ n +2φ x (Dφ)φ n ]}. where we used the shorthand notation φ n = φ/ t n. 7

Finally, we have the recursion operator R = 6 x +6(Dφ) 4 x +9(Dφ x) 3 x +6φ xx 2 x D+{5(Dφ xx)+9(dφ) 2 } 2 x +{9φ xxx +12φ x (Dφ)} x D +{(Dφ xxx )+9(Dφ x )(DΦ)} x +{5φ xxxx +12φ xx (Dφ)+6φ x (Dφ x )}D +{4(Dφ xx )(Dφ)+4(Dφ) 3 3φ xx φ x } +{φ xxxxx +5φ xxx (Dφ)+5φ xx (Dφ x )+2φ x (Dφ xx )+6φ x (Dφ) 2 }D 1 {2(Dφ xx )+2(Dφ) 2 }D 1 φ x 4φ x (Dφ) 1 x φ x 2(Dφ)D 1 [φ xxx +2φ x (Dφ)] 2φ x D 1 {(φ xxx +φ x (Dφ))D 1 3(Dφ)D 1 φ x 2φ x 1 x φ x +2D 1 [φ xxx +2φ x (Dφ)]} Remark: When calculating the coefficients of R n, one should solve a system of differential equations. Due to nonlocality (those underlined terms), there is certain ambiguity and to avoid it, we used the t 7 -flow 5 Conclusion φ t7 = φ xxxxxxx +7φ xxxxx (Dφ)+14φ xxxx (Dφ x )+14φ xxx (Dφ xx ) +14φ xxx (Dφ) 2 +7φ xx (Dφ xxx )+28φ xx (Dφ x )(Dφ) +14φ x (Dφ xx )(Dφ)+7φ x (Dφ x ) 2 + 28 3 φ x(dφ) 3. Summarizing, we find a supersymmetric SK equation which has Lax representation. We also obtain infinite conserved quantities and a recursion operator for this new proposed system. These imply that the system is integrable. It is interesting to establish other properties for it, such as Bäcklund transformation, Hirota bilinear form, etc.. Acknowledgements The calculations were done with the assistance of SUSY2 package of Popowicz [17]. We would like to thank him for helpful discussion about his package. The comments of anonymous referee has been very useful. The work is supported in part by National Natural Science Foundation of China under the grant numbers 10671206 and 10731080. References [1] K. Sawada and T. Kotera, Prog. Theor. Phys. 51 (1974) 1355. [2] P. J. Caudrey, R. K. Dodd and J. D. Gibbon, Proc. R. Soc. London A 351 (1976) 407. [3] J. Satsuma and D. J. Kaup, J. Phys. Soc. Japan 43 (1977) 692. [4] R. K. Dodd and J. D. Gibbon, Proc. R. Soc. London A 358 (1977) 287. [5] B. Fuchssteiner and W. Oevel, J. Math. Phys. 23 (1982) 358. [6] R. N. Aiyer, B. Fuchssteiner and W. Oevel, J. Phys. A 19 (1986) 3755. 8

[7] D. Levi and O. Ragnisco, Inverse Problems 4 (1988) 815. [8] A. P. Fordy and J. Gibbons, Phys. Lett. A 75 (1980) 325. [9] R. Hirota, J. Phys. Soc. Japan 58 (1989) 2285. [10] A. S. Carstea, Nonlinearity 13 (2000) 1645. [11] J. M. Figueroa-O Farrill, E. Ramos and J. Mas, Rev. Math. Phys. 3 (1991) 479. [12] L. A. Dickey, Soliton Equations and Hamiltonian Systems, 2nd Editon, (World Scientific, Singapore(2003)). [13] Yu I Manin and A. O. Radul, Commun. Math. Phys. 98 (1985) 65. [14] P. H. M. Kersten, Phys. Letts, A 134 (1988) 25. [15] P. Labelle and P. Mathieu, J. Math. Phys. 32 (1991) 923. [16] M. Gürses, A. Karasu and V. Sokolov, J. Math. Phys. 40 (1999) 6473. [17] Z. Popowicz, Compt. Phys. Commun. 100 (1997) 277. 9