Lax Representations for Matrix Short Pulse Equations

Lax Representations for Matrix Short Pulse Equations Z. Popowicz arxiv:1705.04030v1 [nlin.si] 11 May 017 May 1, 017 Institute of Theoretical Physics, University of Wrocław, Wrocław pl. M. Borna 9, 50-05 Wrocław Poland, ziemek@ift.uni.wroc.pl Abstract The Lax representation for different matrix generalizations of Short Pulse Equations SPE is considered. The four-dimensional Lax representations of four-component Matsuno, Feng and Dimakis-Müller-Hoissen-Matsuno equations is obtained. The four-component Feng system is defined by generalization of the two-dimensional Lax representation to the four-component case. This system reduces to the original Feng equation or to the two-component Matsuno equation or to the Yao-Zang equation. The three component version of Feng equation is presented. The four-component version of Matsuno equation with its Lax representation is given. This equation reduces the new two-component Feng system. The two-component Dimakis-Müller-Hoissen-Matsuno equations are generalized to the four parameter family of the four-component SPE. The bi-hamiltonian structure of this generalization, for special values of parameters, is defined. This four-component SPE in special case reduces to the new two-component SPE. 1 Introduction The short pulse equation SPE u x,t = u+ 1 6 u3 xx, 1 derived by Schäfer and Wayne [1] as a model of ultra-short optical pulses in nonlinear media has attracted considerable attention. For the first time the SPE has appeared in an attempt to construct integrable differential equations associated with pseudo spherical surfaces []. The associated linear scattering problem appeared for the first time in differential geometry [, 3]. 1

The integrability of the SPE have been studied from various points of view [4, 5, 6]. More specifically it was shown that SPE equation admits a Lax pair and possess a bi-hamiltonian formulation. The problem of mathematical description of the propagation ultra-short pulses has been considered from different points of view. For example, the SPE was generalized in different manners to the -component systems. The first generalization have been considered by Pietrzyk Kanattsikov and Bandelow in the form [7]. u xt = u+ 1 6 u3 +3uv xx, v xt = v + 1 6 v3 +3u v xx, u xt = u+ 1 6 u3 3uv xx, v xt = v 1 6 v3 3u v xx, u xt = u+ 1 6 u3 xx, v xt = v + 1 u v xx. Sakovich [8] presented another -component generalization u xt = u+ 1 6 u3 +uv xx, v xt = v + 1 6 v3 +u v xx, 3 u xt = u+ 1 6 u3 xx, v xt = v + 1 6 u v xx. As was shown, these generalizations are integrable and possess Lax representation and bi-hamiltonian formulation. Recently Matsuno [9, 10] generalized, in two different manners, the SPE to the n-component case as u i,xt = u i + 1 Fu i,x x, F = 1 u i,xt = u i +Fu i,x x 1 Gu i, G = 1 j<k n 1 j,k n c j,k u j u k, 4 c j,k u j,x u k,x. 5 where c j,k are arbitrary constants such as c j,k = c k,j. For these equations Matsuno applied the bilinear method and presented multisolitons solutions in the parametric form. Moreover, Matsuno found the Lax representation, local and nonlocal conservation laws for the equation 5, for special case n =,u 1 = u,u = v and c 11 = c = 0,c 1 = 1 u x,t = u+ 1 vu xx, v x,t = v + 1 uv xx. 6 Matsuno defined the zero curvature condition for the equation 6 as X t Z x + [X,Z] = 0 where 1 ux v X = x u x 7 v x 1+u x v x uv1 ux v Z = x + 1 u 4 x uv u u x uv+v uv1 u x v x 1. 4

Dimakis and Müller-Hoissen [11] studied matrix generalization of the SPE in the form U xt = U + 1 U U x x, 8 where U is matrix valued function. Assuming that U has to be a scalar times the identity matrix, Dimakis and Müller-Hoissen were able to construct the Lax pair. When U = u +v I, where I is an identity matrix, the equation 8 reduces, after the transformations of variables u, v u + v/,u v/i, to u t = 1 u+ 1 uvu x, v t = 1 v + 1 uvv x. 9 This equation, as it was shown by Matsuno [9], is a special case of the system 4 when n =,c 1,1 = c, = 0,c 1, = 1,u 1 = u,u = v. For this equation, Brunelli and Sakovich found [1] the bi-hamiltonian formulation where H 1 = K = u v t = K H1,u H 1,v = J H0,u H 0,v, 10 dx[v 1 u+ 1 8 u v x ], H 0 = 1 dx uv, 11 0 1 u, J = x 1 u x 1 +u x 1 v x 1 0 1 +v x 1 u x v x 1, v x and Matsuno [9] discovered the Lax representation ψ1 ψ x 1 ux = v x 1, ψ1 ψ t = uv+ 1 uvu 4 x 1u uvv x + 1v uv 1 4. 1 Quite different generalization of two-component SPE was proposed by Feng [13] u x,t = u+uu x + 1 u +v u xx, v x,t = v +vv x + 1 u +v v xx. 13 Brunelli and Sakovich in [1] obtained the zero curvature condition X t T t + [X,T] = 0 and the bi-hamiltonian structure for this equation, where 1+ux v X = x u x v x u x v x 1 u x v x T = u +v 1+u x v x + 1 4 u +v u x v x 1 u v u +v u x v x + 1 u v u +v 1+u x v x 1 The last two-component generalizations of SPE were given by Yao-Zeng [14] 14. u x,t = u+ 1 6 u3 xx, v x,t = v + 1 u v x x 15 for which Brunelli and Sakovich found, using the four-dimensional matrices, the Lax representation in [1]. 3

In this paper we would like to study several different generalizations of the twocomponent SPE to the matrix and then to the four-component case. Therefore we investigate the equations obtained from the matrix generalization of the Lax representations eqs. 7,1,14. The paper is organized as follows. In the following two chapters we study the matrix version of the Lax representation of Matsuno and Feng equations. The fourcomponent Matsuno equation is discussed in the first chapter too, and its reduction to the new version of the two-component Feng equation is given. In the third section, we study the four-component version of Feng equation and its different reductions. This generalization describes the interaction between the Feng, Matsuno and Yao- Zeng equations and reduces to the original Feng equation or to the two-component Matsuno equation or to the Yao-Zeng system. Also, the three-component version of Feng equation is obtained as a result of reduction. In the last section we study the matrix Lax representation of the two-component Dimakis-Müller-Hoissen - Matsuno equation. This representation in the case of four-dimensional matrices produces the four parameters family of four-component SPE. For special values of free parameters we defined the bi-hamiltonian structure for this equation. As a result of reduction of this four parameter equations, we obtained a new two-component version of SPE. The last section is a conclusion. The paper contains two appendixes. Lax Representation of Matrix Matsuno Equation Let us consider the Lax representation for matrix Matsuno equation ψ1 Π x = XΠ, Π t = TΠ, Π =, 16 ψ I Ux V X = x U x, V x I +U x V x I Ux V T = x UV +I 1 UVU 4 x U UVV x +V I U x V x UV I 1, 4 where I is identity matrix, U and V are arbitrary dimensional matrices such that UV and U x V x are scalar functions multiplied by identity matrix. From the integrability condition Π x,t = Π t,x we obtained the following matrix equations U t,x = U +VUU x x +V x UU x U x U, 17 V t,x = V +VUV x x +V x UV x U x V. These equations reduce to the Matsuno eq.6 when U,V are scalar functions. The four-component version of the Matsuno equation are concluded from eq.17 assuming that U = 1 v0 +v 1 I iv +v 3 σ 1 v v 3 σ v 1 v 0 σ 3, 18 V = 1 v0 +v 1 I +iv +v 3 σ 1 +v v 3 σ +v 1 v 0 σ 3. 4

where σ i are Pauli matrices and I is identity matrix. v 0,t,x = v 0 + 1 v 1v 0 xx + v v 3 v 0,xx +v 0,x v v 3 x v 0 v,x v 3,x, 19 v 1,t,x = v 1 + 1 v 0v 1 xx + v v 3 v 1,xx +v 1,x v v 3 x v 1 v,x v 3,x, v,t,x = v + 1 v 3v xx + v 0 v 1 v,xx +v,x v 0 v 1 x v v 0,x v 1,x, v 3,t,x = v 3 + 1 v v 3 xx + v 0 v 1 v 3,xx +v 3,x v 0 v 1 x v 3 v 0,x v 1,x. The system eq.19 allows interesting reduction to the two-component case when v 1 = v 0 = 1f +g, v 3 = v = 1 f g f t,x = f +f x f + 1 f +g f xx g t,x = g +g x g + 1 f +g g xx and describes the Feng fields plus additional interaction. 1 f x +g x f +f xg x g, 0 1 f x +g x g +f xg x f, 3 Lax Representation of Matrix Feng Equation We postulate the Lax representation for the matrix Feng equation as Π x = XΠ, Π t = TΠ, Π = I +Ux V X = x U x V x, U x V x I U x V x 1 T = U +V I +U x V x +I 1 4 ψ1 ψ, 1 U +V U x V x + 1 V U U +V U x V x + 1 U V 1 U +V I +U x V x I 1 4 where I is identity matrix, U and V are arbitrary matrices such that UV and U x V x are scalar functions multiplied by the identity matrix. The matrix Feng equation is obtained from the integrability conditionπ x,t = Π t,x and reads U t,x = U + 1 U +V U x x +UU x U x U U x V VU x V x, V t,x = V + 1 Vxx U +V +V UV x +V xv UV x. When U and V are scalar function, then eq. reduces to the Feng equation 13. In the next section we show that this system, in the case where we consider the two dimensional matrices U and V, contains the Feng, Matsuno and Yao-Zheng equations., 5

4 Unification of Feng, Matsuno and Yao - Zeng equations. The four-component Feng equation is concluded from eq. assuming that U = 1 v0 +v 1 I v +v 3 σ 1 +iv v 3 σ v 1 v 0 σ 3, 3 V = 1 v0 +v 1 I +v +v 3 σ 1 iv v 3 σ +v 1 v 0 σ 3, where σ i are Pauli matrices and I is identity matrix. v 0,t,x = v 0 +v 0,xv 0 + 1 v 0 +v 1v 0,xx +v 0,x v v 3 x +v 1 v,x v 3,x, 4 v 1,t,x = v 1 +v1,xv 1 + 1 v 0 +v1v 1,xx +v 1,x v v 3 x +v 0 v,x v 3,x, v,t,x = v +v 3 v 1 xx + v,x v0 +v 1 +v 0,x v 1,x v, x v 3,t,x = v 3 +v v3 1 xx + v3,x v0 +v 1 +v 0,x v 1,x v 3, x We see that this system of equations describes the interaction between Feng fields {v 0,v 1 } and Matsuno fields {v,v 3 }. For this reasons the system of equations could be considered as the interacting system of Feng - Matsuno type. The Matrices X and T in the Lax representation eq.1 are rewritten in terms of the 16 dimensional Lie algebra, spanned by the generators {e 1,e,...,e 16 } as X = 1+v 0,x v 1,x v,x v 3,x e 1 +e + 5 v 0,x v 1,x e 3 +e 4 +v,x e 5 +v 3,x e 6, T = γ1+v 0,x v 1,x v,x v 3,x e 1 +e + 1 4 e 1 +e +γv 0,x v 1,x e 3 +e 4 + γv,x e 5 +γv 3,x e 6 + 1 v 0 v 1 e 7 +e 8 +v e 9 +v 3 e 10, where γ = 1 v 0 +v 1 +v v 3, The explicit representation of the Lie algebra e e,e,...e 16 is given in the appendix A. The system of equations 4 allows three different reductions to the two-component case. When v = v 0 = 1a + b, v 3 = v 1 = 1 a b then eq.4 reduces to the Yao-Zeng eq.15 in the variables a,b. For this type of reduction the matrices X,T are X = e 1 +e +a x e 5 +e 6 +b x e 3 +e 4 +e 5 e 6, 6 T = 1 4 a +1e 1 +e + a a x e 5 +e 6 + 1 ae 9 +e 10 a b x e 3 +e 4 +e 5 e 6 1 be 7 +e 8 +e 9 e 10. 6

This representation is different than the one given in [1]. When v 1 = v 0 = 1a + b,v 3 = v = 1 a b then eq.4 reduces to the Feng eq.13. When v 3 = v 0 = 1a+b,v = v 1 = 1 a b then eq.4 reduces to the Yao-Zeng eq.15 again but now we have different representations of the matrices X,T X = e 1 +e +a x e 5 +e 6 +b x e 3 +e 4 e 5 +e 6, 7 T = a e 1 +e + 1 4 e 1 +e + a a x e 5 +e 6 + 1 ae 10 +e 9 a b x e 3 +e 4 e 5 +e 6 + 1 be 7 +e 8 e 9 +e 10. The system 4 allows the reduction to the three component case as well. For example, assuming that v 0 = a+b+c, v 1 = a b+c 8 v = a+b c, v 3 = v 0 v 1 v then equations 4 reads as a t,x = a+a xx ac+b +ca x +b x c x+4a x b x b+4a x c x a, 9 b t,x = b+b xx ac+b +ba x +b x +c x +4a xb x c+4b x c x a, c t,x = c+c xx ac+b +ab x +c x a x +4a xc x c+4b x c x b. For this type of reduction the matrices X,T are X = 1+a x b x +c x e 1 +e +a x c x e 5 e 6 30 +b x e 3 +e 4 +e 5 +e 6, T = ac+b 1+a x b x +c x e 1 +e + 1 4 e 1 +e c ae 10 e 9 +be 7 +e 8 +e 9 +e 10 +4ac+b b x e 3 +e 4 +e 5 +e 6 +4c x a x ac+b e 6 e 5. 5 Lax representation of Matrix Dimakis-Müller-Hoissen- Matsuno Equation Now let us consider the matrix generalization of Lax pair representation eq.1 I Ux ψ1 Ψ x = ΩΨ = ξ Ψ, Ψ =, 31 V x I ψ Ψ t = PΨ = ξ UV +I 1 4ξ 1 U + ξ 4 U xvu +UVU x 1 V + ξ 4 V xuv +VUV x 1 VU I 1 4ξ Ψ. 7

where I is an N dimensional identity matrix and Ux,t,Vx,t are at the moment arbitrary N and dimensional matrices. The zero-curvature condition Ω t P x + [ Ω,P ] = 0 produces the following equations on U,V U tx = U + 1 4 U xvu +UVU x x, 3 V tx = V + 1 4 V xuv +VUV x x, with the constraints on U,V [ Ux V x,uv ] = 0, [ Vx U x,vu ] = 0 33 The equations 3 with the constraints 33 are our matrix generalization of twocomponent Dimakis-Müller-Hoissen-Matsuno equation different than the one proposed in [11]. Now we consider the case where U and V are two-dimensional matrices. We parametrize U by two functions u 1,u and V by two functions u 3,u 4 as u U = 1 u u, V = 3 u 4 34 s 1 u 1 +s u s 3 u 1 +s 4 u z 1 u 3 +z u 4 z 3 u 3 +z 4 u 4 where s i,z i,i = 1,,3,4 are an arbitrary constants. 1 Substituting U, V defined by eq.34 to the equations 33, we find the connection between z i and s i z 1 = s 1 s 3, z = s s 3 s 1 s 4 s 3, z 3 = 1 s 3, z 4 = s 4 s 3 35 Now the Lax representation generates four parameters family of equations u 1,t = 1 u 1 + 1 u3 s 3 u 1 +s u 4s x +u 4 s 1 s 3 u 1 +s s 4 u +s s 3 u 1 u x, 3 u,t = 1 u + 1 u 4 s 4s 3 u 1 +s s 3 +s 4 s 1 s 4 u +s 3 s 4 u 1 u x 3 u 3,t = 1 u 3 + 1 4 + 1 4s 3 u 4,t = 1 u 4 + 1 4 +u 3 s4 s 1 u +s 3 u 1 u x u 1 u 3 +s s 3 s 1 s 4 u 4 x, 36 u s1 u 3 +s s 3 s 1 s 4 s 4 u 4 +u 3 u 4, x u 1 s1 +s 4 u 4 +u 3 u 4 x + 1 u u 3 +s s 3 +s 4u 4 +s 4 u 3 u 4 s 3 x. 1 In general it is possible to parametrize the two dimensional matrix by two arbitrary functions in six different manners. However, all these parametrizations could be obtained from eq.34 using the linear transformations of these functions. 8

These equations are integrable in the sense that they possess the Lax representation and allow several different reduction to the one or two-component SPE equation case. When u 4 = u 3 = u = u 1 = u,s 1 = s,s 4 = s 3 1 the equations 36 reduce to the original SPE equation 1. For u 1 = u = u,u 3 = u 4 = v,s 1 = 1 s,s 4 = s 3 1 the equations 36 reduce to the two-component Dimakis-Müller-Hoissen-Matsuno 9. When u 1 = u 3 = u,u = u 4 = v,s 1 = s 4 = 0,s 3 = 1 and s = ±1,0 we obtained three generalizations of the SPE equation considered by Pietrzyk at.al. On the other hand for u 1 = u 3 = u,u = u 4 = v,s 1 = 0,s = 1,s 3 = 1,s 4 = we obtained new two-component generalization of the SPE equation u t,x = u+ 1 ux u +v v x u v 37 v t,x = v +u x vu v + 1 v xu +5v 4uv This equation was missing in the classification of Pietrzyk at.all [7]. For this equation the Lax representation is exactly as in the paper [7] I Ux X =, T = U + 1 I 4 6 U3 x 1U U x I 6 U3 x + 1 U 1 U 1 I, 38 4 but now U = u v v u v. 39 We found the bi-hamiltonian formulations for the equations 36 for s 1 = 0,s = 1. u 1,t = 1 u 1 + 1 u3 s 3 u 1 +u 4s x +u 4 s 4 u +s 3 u 1 u x, 40 3 u,t = 1 u + 1 u 4 s 4s 3 u 1 +s 3 +s 4u +s 3 s 4 u 1 u x +u 3 s4 u +s 3 u 1 u, x 3 u 3,t = 1 u 3 + 1 4 u 1 u 3 +s 3 u 4 +u x s 4 u 4 +u 3u 4 x, u 4,t = 1 u 4 + 1 u 1 s4 u 4 4 +u 3u 4 x + 1 u u 3 +s 3 +s 4 s u4 +s 4u 3 u 4 3 x. u 1 u u 3 u 4 t = J 0 H 1,u1 H 1,u H 1,u3 H 1,u4 = J 1 H 0,u1 H 0,u H 0,u3 H 0,u4 41 9

H 1 = dx [ u 1 1 u 3 +u 1 u 4 ] + 1 u +s 3 u 8s 1u 3 x + 4 3 u s 3 +s 4 +s 3u 1 +s 3s 4 u 1 u u 4 x +s 4 u +4s 3u 1 u u 3 u 4 x H 0 = dx u 1 u 3 +u u 4 J 1 = J 0 = 0 0 0 1 1 0 0 0 0 1 0 0 0 0 1 0 0 0 0 1 1 0 0 0 0 1 0 0 43 + 44 s 4 s 3 0 u,x 1 u,x 0 u,x 1 u,x u,x 1 u,x u,x 1 s 4 u,x +s 3 u 1,x +s 3 u 1,x 1 u,x s 3 u,x 1 u 4,x u,x 1 s 4 u 4,x +u 3,x +s 3 u 1,x 1 u 4,x 0 s 3 u 4,x 1 u,x 0 s 3 u 4,x 1 u 4,x u 4,x 1 u,x u 4,x 1 s 4 u,x +s 3 u 1,x +u 3,x 1 u,x s 3 u 4,x 1 u 4,x u 4,x 1 s 4 u 4,x +u 3,x +u 3,x 1 u 3,x + 1 s 3 u,x 1 u,x s 3 u,x 1 u 1,x s 3 u,x 1 u 4,x u,x 1 u,x +s 3 u 1,x 1 u 1,x +s 3 u 1,x 1 u 1,x +s 3 u 1,x 1 u,x +s 3 u 1,x 1 u 4,x s 3 u 1,x 1 u,x +s 3 u 1,x 1 u,x s 3 u,x 1 u,x +s 3u 1,x 1 u+1,x s 3 u,x 1 u 3,x +u 1,x 1 u 4,x s 3 u,x 1 u 4,x +s 3 u 1,x 1 u 3,x s 3 u 4,x 1 u,x s 3 u 3,x 1 u,x s 3 u 3,x 1 u 3,x s 3 u 4,x 1 u 3,x +s 3 u 3,x 1 u 1,x +s 3 u 4,x 1 u 1,x +s 3 u 4,x 1 u 4,x +s 3 u 3,x 1 u 4,x s 3 u 4,x 1 u 1,x s 3 u 4,x 1 u,x s 3 u 4,x 1 u 3,x s 3 u 4,x 1 u 3,x +s 3 u 3,x 1 u,x +s 3 u 3,x 1 u 1,x +s 3 u 3,x 1 u,x +s 3 u 3,x 1 u 3,x The proof that this operator satisfies the Jacobi identity is postponed to the appendix A. Moreover we showed that the Hamiltonians operators J 0 and J 1 are compatible. It means that µj 0 + J 1 is also the Hamiltonian operator, where µ is an arbitrary constant. 10

u0 u U = 1 45 u 1 u 0 The equation 46 for s 4 =,s 3 = 1 possesses very nice bi-hamiltonian structure u 1,t = 1 u 1 + 1 4 u 3u 1 +u x + 1 u 4 u +u 1 u x, u,t = 1 u + 1 4 u 4u 1 +5u 4u 1 u x + 1 u 3 u +u 1 u x, u 3,t = 1 u 3 + 1 4 u 1u 3 +u 4 x + 1 u u 4 +u 3 u 4 x, u 4,t = 1 u 4 + 1 4 u u 3 +5u 4 4u 3u 4 x + 1 u 1 u 4 +u 3u 4 x. u 1 u u 3 u 4 t = H 1 = H 0 = 0 0 1 0 0 0 0 1 1 0 0 0 0 1 0 0 H 1,u1 H 1,u H 1,u3 H 1,u4 = J H 0,u1 H 0,u H 0,u3 H 0,u4, 46 dx [ u 3 1 u 1 +u 4 1 u ] + 1 u 1 8 xu 3 +u 4 + 47 u x u 3 +5u 4 4u 3 u 4 +u 1 u x u 3 u 4 u 4 dx u 1 u 3 +u u 4, 48, U J = 1 U 1 I + U 1 V 1 I + V 1 U V 1 V u1 u U = u3 u, V = 4 1 0,, I = u u 1 u u 4 u 3 u 4 0 1 49. 6 Conclusion In this paper we studied several different generalizations of the two-component SPE to the matrix case and then to the four-component case. In particular we studied four-component version of Feng, Matsuno and Dimakis-Müller-Hoissen-Matsuno equations with its different reductions. The four-component Feng equation is very reticular, because it contains the Feng, Matsuno and Yao-Zeng equations and thus unifies these equations. The four parameters of four-component Dimakis-Müller-Hoissen-Matsuno equations were discussed. For special values of free parameters, we obtained the bi- Hamiltonian structure and presented new two-component SPE 37. We obtained several new integrable equations in the sense that these possess the Lax representation and this open perspective for further investigations. For example the problem of finding the bi-hamiltonian structures for other generalized SPE is still an open issue, similarly as the problem of presentation of their solutions. 11

7 Appendix A. Representation of e 1,e,...e 16 e 1 = e 4 = e 7 = e 10 = e 13 = e 16 = 1 0 0 0 0 0 1 0 0 0 0 1 0 1 0 0 0 0 1 0 1 0 0 0 0 0 0 1 0 1 0 0 1 0 0 0 0 1 0 0 e = e 5 = e 8 = e 11 = e 14 = 0 1 0 0 0 0 0 1 1 0 0 0 0 0 0 1 0 1 0 0 1 0 0 0 1 0 0 0 e 3 = e 6 = e 9 = e 1 = e 15 = 1 0 0 0 0 0 0 1 0 1 0 0 1 0 0 0 0 1 0 0 0 0 0 1 0 0 1 0 0 0 0 1 8 Appendix B. Hamiltonian operators and verification of the Jacobi identity. We assume the most general form on the second Hamiltonian operator as J 1 = J 0 +Ĵ where Ĵ j,s = 4 c j,s,k,l u k,x 1 u l,x, j s 50 k,l=1 Ĵ s,j = Ĵj,s j > s, J 0 = 1 0 0 0 1 1 0 0 0 51 0 1 0 0 1

where c j,k,s,l are arbitrary constants j,k,s,l = 1...4 for j k and c j,j,s,k = c j,j,k,s. Now we will investigate the equations obtained from the second Hamiltonian structure u 1 H 0,u1 u u 3 = J H 0,u H 0,u3 5 u 4 H 0,u4 t and compare the solutions with the solutions obtained from Lax representation eq. 31. This fixes all coefficients c i,j,k,l and moreover restricts s 1 = 0,s = 1. We use traditional manner to verify the Jacobi identity [15]. In order to prove that the operator J 1, defined in 44, satisfies the Jacobi identity, we utilize the standard form of the Jacobi identity Jacobi = dxaj 1 J 1 B C +cyclica,b,c = 0, 53 where A,B and C are the test vector functions, for example A = a 1,a,a 3.a 4 while denotes the Gato derivative along the vector LB. We check this identity utilizing the computer algebra Reduce and package Susy [16]. We will, briefly explain our procedures applied for the verification of the Jacobi identity 53. In the first stage, we remove the derivatives from the test functions a i,x,b i,x and c i,x in the Jacobi identity, using the rule f x = f f 54 where f is an arbitrary function. Because of this, the Jacobi identity can be split into three segments. The first and second segments contain terms in which the integral operator 1 appears twice and once respectively. The last segment does not contain any integral operators. We consider each segment separately. The first segment is a combination of the following expressions: dx n c 1 n a 1 n b + dx ñ a 1 ñ c 1 ñ b +cyclica,b,c where n a denotes n j a i or n j,x a i or n j,xx a i, i,j = 1,,3,4 and similarly for n b,n c, ñ a,ñ b,ñ c. Here 1 is an integral operator, and, therefore, each ingredient could be rewritten as, for example, dx n c 1 n a 1 n b = dx n a 1 n c 1 n b, 55 Now, we replace n a in the last formula by n a = 1 n a 1 n a. Hence, the expression 55 transforms to dx n c 1 n a 1 n b = dx, n c 1 n a 1 n b + dx n b 1 n a 1 n c 13

Now repeating this procedure for n a and ñ a in the first segment, it appears that this segment reduces to zero. The second segment is constructed from the combinations of the following terms: dx Λ a Λ c 1 Λ b + dx Λ b 1 Λa Λc +cyclica,b,c. where Λ a takes values in {n j a i, n j,x a i, n j,xx a i, n j a i,x, n j a i,xx },i,j = 1,,3,4. In a similar manner the Λ b,λ c, Λ a, Λ b, Λ c are defined. These terms are rewritten as dx Λ a Λ c 1 Λ b dx Λ a Λc 1 Λb +cyclica,b,c. Next, using rule 54, we replace once more the derivatives in a k,x and b k,x in the second segment. After this replacement, it appears that the second segment contains no term with an the integral operator. Therefore, we add this segment to the third segment. Now, it is easy to check that this last segment vanishes. Indeed, it is enough to use the rule 54 in order to remove the derivatives from a k,x in the last segment. This finishes the proof. References [1] T. Schäfer and C. E. Wayne, propagation of ultra-short optical pulses in cubic nonlinear media, Physics D 196 004 90-105. [] M.L.Rabelo, On equations which describe pseudospherical surfaces Stud. Appl.Math. 81 1-48 1989. [3] R. Beals, M.L.Rabelo and K. Tenenblat, Bäcklund transformation and inverse scattering solutions for some pseudospherical surfaces equations, Stud. Appl.Math. 81 15-151 1989. [4] A. Sakovich and S. Sakovich, The short pulse equations is integrable, J. Phys. Soc. Jpn. 74 30-41 005. [5] J. C. Brunelli The short pulse hierarchy, J. Math. Phys. 46 13507 005. [6] J. C. Brunelli The bi-hamiltonian structure of the short pulse equation, Phys. Letta Z 353 475-478 006. [7] M. Pietrzyk, I. Kanatsikov and U. Bandelow On the propagation of vector ultrashort pulses, J. Nonlinear Math. Phys. 15, 16 008. [8] S. Sakovich On integrability of the vector short pulse equation J. Phys. Soc. Jpn. 77 13001 008. [9] Y. Matsuno A novel multi-component generalization of the short pulse equation and its multisoliton solutions, J. Math. Phys. 5 1370 0111. 14

[10] Y. Matsuno Integrable multi-com-ponent generalization of a modified short pules equations, J. Math. Phys. 57 111507 0116. [11] A. Dimakis and F. Müller-Hoissen Bidifferentional Calculus Approach to AKNS Hirarchies and Their Solutions Symmetry, Intgrability and Geometry: Methods and Applications 6 055 010. [1] J.C. Brunelli and S. Sakovich Hamiltonian Integrability of Two-Component Short Pulse Equations J. Math. Phys. 54 013 01701. [13] B.F. Feng An integrable coupled short pulse equation J. Phys.A 45 0850 01. [14] Y. Yao and Y.Zeng Coupled Short Pulse Hierarchy and Its Hamiltonian Structure J. Phys.Soc.Jpn. 80, 064004 011. [15] M. Blaszak Multi-Hamiltonian Theory of Dynamical Systems Springer 1998. [16] Z. Popowicz SUSY, Computer Physics Communications 100, 77-96 1997. 15