arxiv: v1 [math-ph] 10 Sep 2015

Transcription

1 A N = 2 etension of the Hirota bilinear formalism and the supersymmetric KdV equation arxiv: v [math-ph] 0 Sep 205 Laurent Delisle,2 September, 205 Abstract We present a bilinear Hirota representation of the N = 2 supersymmetric etension of the Korteweg-de Vries equation. This representation is deduced using binary Bell polynomials, hierarchies and fermionic limits. We, also, propose a new approach for the generalisation of the Hirota bilinear formalism in the N = 2 supersymmetric contet. Introduction The study of eact solutions of completely integrable supersymmetric systems is of current interest in modern mathematical physics research. In particular, the N = 2 supersymmetric etension of the Korteweg-de Vries(KdV) equation[] has been largely studied in terms of integrability conditions, eact solutions and symmetry group structures[2, 3, 4, 5, 6, 7, 8, 9, 0]. The equation is described by a bosonic superfield A defined on the superspace R 2 2 [] of local coordinates (,t,θ,θ 2 ). The variables (,t) are the usual Euclidean space-time coordinates and (θ,θ 2 ) are real Grassmann coordinates satisfying the usual anti-commutation relations θ θ 2 +θ 2 θ = 0, θ 2 = θ 2 2 = 0. () The superfield A satisfies the Labelle-Mathieu N = 2 supersymmetric etension of the KdV equation [] A t = ( A +(a+2)ad D 2 A+(a )(D A)(D 2 A)+aA 3 ), (2) where a R is a real parameter, D, D 2 are the covariant derivatives defined as D i = θi +θ i, i =,2 (3) and satisfy D 2 = D 2 2 =. The bosonic superfield A can be decomposed [] using a Taylor epansion as A(,t,θ,θ 2 ) = u(,t)+θ ξ (,t)+θ 2 ξ 2 (,t) θ θ 2 v(,t), (4) where u, v are comple-valued even functions and ξ, ξ 2 are comple-valued odd functions. Labelle and Mathieu [] showed that equation (2) is completely integrable for the special choices a = 2,, 4. This fact suggests that, for these special values, the supersymmetric KdV equation possesses travelling wave and multi-soliton solutions [2, 3, 4, 5, 6, 7]. An algebraic direct method to find such solutions is described by the Hirota bilinear formalism. This method was used numerous time for non-supersymmetric integrable evolution equation to construct soliton and similarity solutions, Bäcklund and Darbou transformations and Institut de mathématiques de Jussieu-Paris Rive Gauche, UP7D-Campus des Grands Moulins, Bâtiment Sophie Germain, Cases 702, Paris Cede laurent.delisle@imj-prg.fr

2 to obtain integrability conditions [2]. Carstea has adapted this formalism to N = supersymmetric etensions [3] such as the KdV, modified KdV and Sine-Gordon equations. The generalisation of this formalism to N = 2 etensions has been confronted to numerous difficulties. Zhang and al. [5] used the strategy of decomposing equation (2) into two N = equations for which the Hirota formalism is well adapted. The way of achieving this is to re-write the bosonic superfield A given in (4) as A(,t,θ,θ 2 ) = A 0 (,t,θ )+θ 2 Ξ(,t,θ ), (5) where A 0 (,t,θ ) = u(,t)+θ ξ (,t) and Ξ(,t,θ ) = ξ 2 (,t)+θ v(,t) are, respectively, bosonic and fermionic superfields of (,t,θ ) R 2. In order to use the bilinear Hirota formalism [2], we have to re-write the superfields A 0 and Ξ in terms of dimensionless bosonic superfields. This is done by dimensional analysis following the fact that equation (2) is invariant under the dilatation vector field [2] +3t t + 2 θ θ + 2 θ 2 θ2 A A. (6) This vector field shows that, under the transformation ( ) (,t,θ,θ 2,A) λ,λ 3 t,λ 2 θ,λ 2 θ2,λ A, (7) equation (2) is invariant, where λ is a free parameter. Under these transformations, we deduce the dimension of these quantities [A] =, [ ] =, [D ] = 2 (8) and this allows us to re-write the superfields A 0 and Ξ as A 0 = B, Ξ = D p, (9) where B = B(,t,θ ) and p(,t,θ ) are bosonic superfields such that [B] = [p] = 0. Introducing the superfield A, given by (5) together with the relations (9), in equation (2), we get two N = supersymmetric equations B t = B +(a+2)b p +(a )(D B )(D p )+ab 3, (0) D p t = D p +3p D p (a+2)b D B +( a)b D B +3aBD 2 p (). The study of these two equations is the main object of the paper for the special cases a = 2,,4, i.e. for the cases for which equation (2) is completely integrable []. Indeed, we will give a bilinear Hirota representation of these equations using different approaches using the binary Bell polynomials [4], hierarchies [8, 9, 0] and fermionic limits [4]. The later approach is simple, it relies on taking ξ = ξ 2 = 0 in the representation of the bosonic superfield A given in (4). The binary Bell polynomials have, recently, found a connection with the Hirota bilinear formalism [4] and this connection will be used throughout this paper. The one-variable Bell polynomials Y are defined as Y k,k tt,k θ (f) = e f k kt t Dk ef, (2) where k µ are integer constants and f = f(,t,θ ) is a bosonic superfield. Using the polynomials Y, we define the binary Bell polynomials Y as Y k,k tt,k θ (w,w 2 ) = Y k,k tt,k θ (f k k tt k θ ), (3) where the different derivatives of f are replaced by the superfields w and w 2 following the procedure { w f k k tt k θ =, k k tt k θ if k + k t + k is odd, if k + k (4) t + k is even. w 2, k k tt k θ 2

3 Note here that we are using the notation f k k tt k θ = k k t Hirota bilinear formalism is given by Y k,k tt,k θ (w = ln where the Hirota derivative is defined as ( ) ) f,w 2 = ln(fg) g t D k f. The link with the = (fg) S k Dk Dt kt (f g), (5) S k Dk Dkt t (f g) = (D D )k ( ) k ( t t ) kt f(,t,θ )g(,t,θ ) (,t,θ )=(,t,θ). (6) These relations will be used to transform equations (0) and () into binary Bell polynomial equations. To achieve this, we will need auiliary tools such as hierarchies of the supersymmetric KdV equation (2) and fermionic limits. This paper is divided as follows. In the following three sections, we give a bilinear representation of the supersymmetric KdV equation (2) for, respectively, a =, a = 4 and a = 2. In section 2, we use directly the binary Bell polynomials to get a general Hirota formulation. In the a = 4 case, we use the binary Bell polynomials and the Two-Boson supersymmetric equation [5, 8], which are members of the same hierarchy, to obtain a Hirota equation. In section 4, we get the bilinear representation using fermionic limits and retrieve the bosonic Miura transformation [2] linking a solution of the KdV equation with the modified KdV equation. The last section addresses the open problem of generalizing the Hirota bilinear formalism to the N = 2 supersymmetric contet. We ehibit this generalisation threw the supersymmetric KdV equation with a =. We conclude the paper with some future outlooks and remarks. The novelty of this paper is based on the use of the binary Bell polynomials (3) to get a bilinear representation of the supersymmetric KdV equation (2). As of today, the Hirota formulation of the supersymmetric KdV equation with a = 2 was an open problem and here we propose a partial answer to this question. We also, for the first time, generalize the Hirota bilinear formulation to N = 2 etensions of certain integrable systems. 2 The supersymmetric KdV equation with a = In this section, we directly use the binary Bell polynomials to get a bilinear representation of the supersymmetric KdV equation (2) with a =. In this case, the two N = equations (0) and () reduces to B t = B +3B p +B 3, (7) D p t = D p +3p D p 3B D B +3B 2 D p. (8) We can notice from these equations that the superfield B is associated to odd numbers of derivatives while the superfield p to even numbers. This fact is compatible with the binary Bell polynomials [4]. Indeed, we have the epressions Y t (cb,dp) = cb t, Y 3 (cb,dp) = cb +3cdB p +c 3 B 3, (9) from which we easily deduce an equivalent representation of equation (7) given as Y t (ib, p)+y 3 (ib, p) = 0. (20) For the second equation (8), we have the following binary Bell polynomials, assuming equation (7) is satisfied, Y tθ (ib, p) = D p t +B D B 3B p D B B 3 D B, (2) Y θ (ib, p) = D p +B 3 D B +3B p D B B D B +3B 2 D p + 3p D p 3B D B (22) 3

4 and it is direct to show that equation (8) is equivalent to the binary Bell polynomials equation Y tθ (ib, p)+y θ (ib, p) = 0. (23) Using the link between the binary Bell polynomials and the Hirota derivative(5), we obtain, casting the change of variables ib = ln(f/g) and p = ln(fg), the bilinear representation of the supersymmetric KdV equation with a = [5] given as (D t +D 3 )(f g) = 0, S (D t +D 3 )(f g) = 0. (24) 3 The supersymmetric KdV equation with a = 4 In this case, we give a Bell polynomial perspective to the supersymmetric KdV equation with a = 4 using its integrable hierarchy [8, 0]. Unlike the a = case, the Bell polynomial approach may not be directly applied to equation (0) and () with a = 4. Indeed, these equations are eplicitly given as B t = B +6B p +3(D B )(D p )+4B 3, (25) D p t = D p +3p D p 6B D B 3B D B +2B 2 D p, (26) and the terms (D B )(D p ), B D B are incompatible with the definition of the binary Bell polynomials. These terms are incompatible in the sense that the superfield B is associated to odd number of derivatives in the binary Bell polynomial perspective. So, we have to find an other way of writing these problematic terms. One way of achieving this is by considering the integrable hierarchy associated to the supersymmetric KdV equation with a = 4. One member of this hierarchy is the Two-Boson supersymmetric equation [5, 8] for which its flow commutes with the flow associated to equation (2) with a = 4. Regarding this matter, the Two-Boson supersymmetric equation is given by w t2 = ( w +w 2 +2D ρ), (27) ρ t2 = (ρ +2wρ), (28) where w and ρ are, respectively, bosonic and fermionic superfields. Making use of the dilatation invariant vector field, we cast the change of variables w = C and ρ = D q, where C and q are dimensionless bosonic superfields. In this case, after integration with respect to the variable, equations (27) and (28) reduce to C t2 = C +C 2 +2q, (29) D q t2 = D q +2C D q. (30) At first sight, this system is again incompatible with the binary Bell polynomials, but, taking m = 2q C and n = C, these equations are given as n t2 = m +n 2, (3) D m t2 = D n +2n D m. (32) Wethusobservethatmisassociatedtoevennumbersofderivatives, whilentoodd numbers. From the binary Bell polynomials Y t2 (n,m) = n t2, Y θ (n,m) = m D n+n 2 D n+d n +2n D m (33) Y (n,m) = m +n 2, Y t2θ (n,m) = D m t2 +m D n+n 2 D n, (34) it is direct to see that the system of equations (3) and (32) is equivalent to the system of binary Bell polynomials Y t2 (n,m) Y (n,m) = 0, Y t2θ (n,m) Y θ (n,m) = 0. (35) 4

5 Making the change of variables n = ln(f/g), m = ln(fg) and using the link between the Bell polynomials and the Hirota derivative (5), we get the bilinear representation of the Two-Boson equation [5] given as (D t2 D 2 )(f g) = 0, S (D t2 D 2 )(f g) = 0. (36) Let us now focus our attention on the first equation(25) of the supersymmetric KdV equation for a = 4. We have, taking n = αb and m = βp, the following binary Bell polynomial epression: Y t2(n,m) = Y t2(αb,βp) = αb +3αβB p +α 3 B 3 +2αβ(D B )(D p ), (37) where we have retrieve the previously problematic term (D B )(D p ). Supposing that equation (25) has the binary Bell polynomial representation α Y t(αb,βp) = γy (αb,βp)+δy t2(αb,βp), (38) for α,β,γ,δ C constants, we directly find that α = 2i, β = 2, γ = i 8, δ = 3i 8 (39) and, thus, equation (25) as the following Hirota bilinear formulation, taking 2iB = ln(f/g) and 2p = ln(fg), (D t + 4 D D D t2 ) (f g) = 0. (40) The Bell polynomial analysis of the second equation (26) is similar as for the first one. Indeed, it lies on the following binary Bell polynomial epressions: β Y tθ (αb,βp) = D p t 2B D B +2B p D B +6(D B )(D p )D B + 8B 3 D B, β Y t 2θ (αb,βp) = 2B D B 8(D B )(D p )D B 2B p D B 8BD 3 B + D p +4B D B +6B D B 2BD 2 p 2p D p, β Y θ (αb,βp) = 2B D B 2B p D B 8B 3 D B +D p 6p D p + 6B D B 2B 2 D p. From these epressions, it is elementary algebra to show that equation (26) is equivalent to the binary Bell polynomials equation Y tθ (2iB, 2p)+ 4 Y θ (2iB, 2p)+ 3 4 Y t 2θ (2iB, p) = 0 (4) and, using the same change of variables for B and p as for equation (25), we get the bilinear representation S ( D t + 4 D D D t2 )(f g) = 0. (42) Hence, equations (40) and (42) represent the Hirota formulation of the supersymmetric KdV equation with a = 4 [5], where the functions B, p, f and g are related as B = i ( ) f 2 ln, p = ln(fg). (43) g 2 5

6 4 The supersymmetric KdV equation with a = 2 As of today, a Hirota bilinear representation of the supersymmetric KdV equation with a = 2 was an open problem [5]. In this section, we propose a partial answer to this problem by considering the fermionic limit [4] of the superfield A given in (4). This means that we take ξ = ξ 2 = 0 in (4) and we consider bosonic superfield A of the form A(,t,θ,θ 2 ) = u(,t) θ θ 2 v(,t) (44) solution of equation (2) with a = 2. In this case, the bosonic comple-valued functions u and v satisfy the system of partial differential equations u t = u 6u 2 u, (45) v t = v +6vv +6u u 6u 2 v 2uu v. (46) We can make some observations on this system: equation (45) is the bosonic modified KdV equation for which a Hirota formulation is known and, as a second observation, we have that, taking u = 0, this system reduces to the bosonic KdV equation [2]. Again, at first sight, this system is incompatible with the definition of the binary Bell polynomials. To solve this problem, we use a Miura-type transformation relating the supersymmetric KdV equation with a = 2 to the equation [0] Q t = Q + 2 Q (D 2Q )(D 2 Q)Q 3 4 (D 2Q )(D Q)(D D 2 Q) 3 4 (D 2Q)(D Q )(D D 2 Q)+ 3 4 (D Q )(D Q)Q, (47) where Q is a bosonic superfield. The equation (47) is the first non-trivial flow of a N = 2 supersymmetric hierarchy, as shown by Tian and Liu[0], and the Miura-type transformation is eplicitly given as A = 2 D D 2 Q+ 4 (D 2Q)(D Q). (48) Using the Taylor epansion Q(,t,θ,θ 2 ) = q 0 (,t) θ θ 2 q 2 (,t) (fermionic limit), the Miura-type transformation (48) is equivalent, in components, to u = 2 q 2, v = 4 (q2 2 +q 2 0,) 2 q 0,. (49) Note that setting u = 0 in these transformations leads to the classical Miura transformation v = 4 q2 0, 2 q 0, relating a solution of the modified KdV equation to a solution of the KdV equation [2]. This observation is compatible with the fact that the bosonic comple-valued functions q 0 and q 2 satisfy the decoupled system of partial differential equations q 0,t = q 0, + 2 q3 0,, (50) q 2,t = q 2, 3 2 q2 2 q 2,. (5) These two equations are all of the modified KdV-type and we know that they possess a Hirota bilinear representation[2]. Indeed, they are equivalent to the binary Bell polynomial equations Y t (/2q 0,β q 0 )+Y (/2q 0,β q 0 ) 3Y (/2q 0,β q 0 )Y (/2q 0,β q 0 ) = 0, (52) Y t (i/2p 2,γ p 2 )+Y (i/2p 2,γ p 2 ) 3Y (i/2p 2,γ p 2 )Y (i/2p 2,γ p 2 ) = 0, (53) where q 2 = p 2,, β and γ are arbitraryconstants and q 0 and p 2 are auiliaryand arbitrary bosonic functions. Using the change of variables ( ) ( ) f 2 q i f 0 = ln, β q 0 = ln(fg), g 2 p 2 = ln, γ p 2 = ln( f g), (54) g 6

7 we get the Hirota bilinear representation (D t +D 3 )(f g) = 3λD (f g), D 2 (f g) = λfg, (55) (D t +D 3 )( f g) = 3 λd ( f g), D 2 ( f g) = λ f g, (56) where λ and λ are arbitrary parameters and the functions u and v take the eplicit forms ( ) f u = i ln, v = ( lnf/g) 2 ( ln g f/ g) 2 lnf/g. 2 (57) In the reminder of this section, we give plots of different soliton solutions. For a -soliton profile, the functions f, g, f and g, solution of the bilinear equations (55) and (56) with λ = λ = 0, may be chosen as f = +e η, g = e η, f = +ie η, g = ie η, (58) where η = κ κ 3 t and η = κ κ 3 t. In the case of a 2-soliton profile, the functions f 2, g 2, f 2 and g 2 given as f 2 = +e η +e η2 +A 2 e η+η2, g 2 = e η e η2 +A 2 e η+η2, (59) f 2 = +ie η +ie η2 Ã2e η+ η2, g 2 = ie η ie η2 Ã2e η+ η2 (60) solvesthe bilinearequations (55) and (56) forλ = λ = 0, whereη i = κ i κ 3 i t, η i = κ i κ ( ) 3 i t 2 ( ) κ for i =,2, A 2 = κ 2 κ +κ 2 and κ 2. Ã 2 = κ 2 κ + κ 2 In the figures below, we use the notation ( ) fn u (m,n) = i ln, v (m,n) = ( lnf m /g m ) 2 ( ln g f n / g n ) 2 2 lnf m/g m (6) n and we have made the choices κ =, κ = 4 5, κ = 3 5, κ 2 = 2, κ = 3 4 and κ 2 = 2 3. Figure represents the functions u (,) and u (2,2), figure 2 gives the behavior of the functions v (,) and v (2,2), while figure 3 plots the functions v (,2) and v (2,). 0.8 u, 0.7 u 2, Figure : The plot of the functions u (,) and u (2,2). The dashed curve correspond to time t = 2, the thick curve to time t = 0 and the other curve to time t = 2. 5 The N = 2 formalism: a first approach In this section, we make an attempt to formalize the Hirota bilinear approach in the N = 2 supersymmetric contet using the Bell polynomials. We will illustrate this new approach on the supersymmetric KdV equation (2) with a =. Using the change of variable A = M and integrating once with respect to the variable, the supersymmetric KdV equation for a = reads as M t = M +3M D D 2 M +M 3, (62) 7

8 0.3 v, 0.4 v 2, Figure 2: The plot of the functions v (,) and v (2,2). The dashed curve correspond to time t = 2, the thick curve to time t = 0 and the other curve to time t = 2. v,2 0.6 v 2, Figure 3: The plot of the functions v (,2) and v (2,). The dashed curve correspond to time t = 2, the thick curve to time t = 0 and the other curve to time t = 2. where M = M(,t,θ,θ 2 ) is a dimensionless bosonic superfield. In orderto adapt the Hirota bilinear formalism in the N = 2 contet, we make the following observation D i = D i = θ i +θ i = D D 2 M = 2 (D D 2 M ) = C, (63) where C = D D 2 M and is the integration operator with respect to. In this setting, the equation (62) can be re-written as M t = M +3M C +M 3 (64) and we observe that the Bell polynomial approach can be directly used. Indeed, the above equation is equivalent to the binary Bell polynomial equation Y t (im, C)+Y (im, C) = 0 (65) which, using the identification im = ln(f/g) and C = ln(fg), can be written as (D t +D 3 )(f g) = 0, (66) where f = f(,t,θ,θ 2 ) and g = g(,t,θ,θ 2 ) are bosonic superfields. The identification imposes a further constraint given as ( ) f D 2 D ln(fg) = i ln (67) g and it can be shown that this relation is equivalent to the bilinear Hirota equations S 2 S (f f) = 2iff, S 2 S (g g) = 2igg. (68) In order to study this additional constraint, we consider the following Taylor epansions ln(f) = f 0 +θ f +θ 2 f 2 +θ θ 2 f 2, ln(g) = g 0 +θ g +θ 2 g 2 +θ θ 2 g 2 (69) 8

9 and, once introduce in (67), we get ln(f) = f 0 +(θ iθ 2 )f +iθ θ 2 f 0,, ln(g) = g 0 +(θ +iθ 2 )g iθ θ 2 g 0,. (70) For the -soliton solution, we can choose f = +ae η+θζ+θ2ζ2+θθ2m2, g = +be η+θν+θ2ν2+θθ2n2 (7) as a solution of the bilinear equation (66) for a,b,m 2,n 2 bosonic quantities, ζ, ζ 2, ν, ν 2 fermionic quantities and η = κ κ 3 t. These free parameters have to be determined in order that the constraint (67) be satisfied. This can be done by substituting the above epressions for f and g in (70). We get the eplicit forms e η f 0 = ln(+ae η ), f = aζ +ae η, ζ 2 = iζ, m 2 = iκ, (72) e η g 0 = ln(+be η ), g = bν +be η, ν 2 = iν, n 2 = iκ (73) and, as a consequence, the components of the superfield A defined in (4) are given as u = i(g 0, f 0, ), ξ = i(g, f, ), ξ 2 = (f, +g, ), v = (f 0, +g 0, ). (74) In this section, we have produced, for the first time, a N = 2 supersymmetric Hirota representation of the KdV equation with a =. Indeed, the Hirota formulation is described by equations (66) and (68). We have thus partially solved the open problem of finding a supersymmetric N = 2 generalisation of the Hirota bilinear formalism [5]. 6 Future outlooks and remarks In this paper, we have presented a systematic way of finding the Hirota bilinear representation of the N = 2 supersymmetric KdV equation using its decomposition into two N = supersymmetric equations and the Bell polynomials. For the completely integrable cases a =,4, we have obtain a complete representation using the integrable hierarchy and the binary Bell polynomials. For the a = 2 case, it was an open problem to find an Hirota formulation. Here, we have proposed a partial answer to this question using fermionic limits and have retrieved the well known Miura transformation relating the KdV and modified KdV equations. It still remains to find its general representation. An other open problem was to find a N = 2 generalisation of the Hirota bilinear formalism. We have, for the first time, succeeded in given such a representation for the supersymmetric KdV equation with a =. The main idea in this construction was to re-write the operator D D 2 as a second derivative with respect to of a given quantity. This had the effect of transforming the N = 2 equation into a new equation involving only derivatives of the bosonic variables (, t) for which the binary Bell polynomials could be directly applied. This as led to its N = 2 Hirota representation. Our future goals is to generalize the results of section 5. This new approach avoids treating a N = 2 equation as two N = supersymmetric equations and this allows us to obtain a N = 2 bilinear representation of the equation. As a final eample to illustrate the efficiency of this new proposed procedure, we consider the N = 2 supersymmetric etension of the potential Burgers equation M t = D D 2 M + 2 M2. (75) As in section 5, we re-write the quantity D D 2 M as C where C = D D 2 M and, as a consequence, the Burgers equations [4] reads as M t = C + 2 M2. (76) 9

10 This equation may be cast into a binary Bell polynomials equation as 2λY t (λm,2λ 2 C) Y (λm,2λ 2 C) = 0, (77) where λ is a free comple parameter. Making the change of variables λm = ln(f/g) and 2λ 2 C = ln(fg), we get the Hirota bilinear equation (2λD t D 2 )(f g) = 0 (78) together with the constraint D 2 D ln(fg) = 2λ ln ( ) f, (79) g which can be re-write into a Hirota bilinear form as S 2 S (f f) = 4λff, S 2 S (g g) = 4λgg. (80) In conclusion, equations (76) and (80) constitute the Hirota bilinear representation of the supersymmetric potential Burgers equation (76). Acknowledgement The author acknowledges a Natural Sciences and Engineering Research Council of Canada (NSERC) postdoctoral fellowship. The author would like to thank its Ph. D. advisor Véronique Hussin for helpful discussions. References [] P. Labelle and P. Mathieu (99) A new supersymmetric Korteweg-de Vries equation, J. Math. Phys. 32, [2] M. A. Ayari, V. Hussin and P. Winternitz (999) Group invariant solutions for the N = 2 super Korteweg-de Vries equation, J. Math. Phys. 40, [3] S. Ghosh and D. Sarma (200) Soliton solutions for the N = 2 supersymmetric KdV equation, Phys. Lett. B 522, [4] V. Hussin and A. V. Kiselev (2009) Virtual Hirota s multi-soliton solutions of N = 2 supersymmetric Korteweg-de Vries equations, Theor. Math. Phys. 59, [5] M.-X. Zhang, Q. P. Liu, Y.-L. Shen and K. Wu (2008) Bilinear approach to N = 2 supersymmetric KdV equations, China Series A: Mathematics, 52, [6] L. Delisle and V. Hussin (202) New solution of the N = 2 supersymmetric KdV equation via Hirota methods, J. Phys.: Conf. Ser. 343, [7] L. Delisle and V. Hussin (202) Soliton and similarity solutions of N = 2,4 supersymmetric equations, Symmetry 4(3), [8] V. Hussin, A. V. Kiselev, A. O. Krutov and T. de Wolf (200) N = 2 supersymmetric a = 4- Korteweg-de Vries hierarchy derived via Gardner s deformation of Kaup- Boussinesq equation, J. Math. Phys. 5, [9] Q. P. Liu (997) On the integrable hierarchies associated with the N = 2 super W n algebra, Phys. Lett. A 4, [0] K. Tian and Q. P. Liu (202) The transformations between N = 2 supersymmetric Korteweg-de Vries and Harry Dym equations, J. Math. Phys. 53,

11 [] J. F. Cornwell, Group theory in Physics: Supersymmetries and Infinite-Dimensional Algebras, Techniques of Physics Vol. 3, Academic, New York, 989. [2] M. J. Ablowitz and H. Segur, Solitons and the Inverse Scattering Transform, SIAM, 98. [3] A. S. Carstea (2000) Etension of the bilinear formalism to supersymmetric KdV-type equations, Nonlinearity 3, [4] E. Fan and Y. C. Hon (202) Super etension of Bell polynomials with applications to supersymmetric equations, J. Math. Phys. 53, [5] Q. P. Liu and X.-X. Yang (2006) Supersymmetric two-boson equation: Bilinearization and solutions, Phys. Lett. A 35, 3 35.