arxiv:math-ph/ v2 21 Feb 2006

Transcription

1 SOLVING BI-DIRECTIONAL SOLITON EQUATIONS IN THE KP HIERARCHY BY GAUGE TRANSFORMATION arxiv:math-ph/050308v Feb 006 JINGSONG HE, YI CHENG AND RUDOLF A. RÖMER Department of Mathematics, University of Science and Technology of China, Hefei, 3006 Anhui, P.R. China Department of Physics and Centre for Scientific Computing, University of Warwic, Coventry CV4 7AL, United Kingdom Abstract. We present a systematic way to construct solutions of the n = 5-reduction of the BKP and CKP hierarchies from the general τ function τ n of the KP hierarchy. We obtain the one-soliton, two-soliton, and periodic solution for the bi-directional Sawada-Kotera bsk, the bi-directional Kaup-Kupershmidt and also the bi-directional Satsuma-Hirota equation. Different solutions such as left- and right-going solitons are classified according to the symmetries of the 5th roots of e iε. Furthermore, we show that the soliton solutions of the n- reduction of the BKP and CKP hierarchies with n = j, j =,, 3,..., can propagate along j directions in the space-time domain. Each such direction corresponds to one symmetric distribution of the nth roots of e iε. Based on this classification, we detail the existence of two-pea solitons of the n-reduction from the Grammian τ function of the sub-hierarchies BKP and CKP. If n is even, we again find two-pea solitons. Last, we obtain the stationary soliton for the higher-order KP hierarchy.. Introduction The Kadomtsev-Petviashvili KP hierarchy is of central interest for integrable systems and includes several well-nown partial differential equations such as the Korteweg-de Vries KdV and the KP equation. With pseudo-differential Lax operator L given as 3] the corresponding generalized Lax equation L = u u 3, L t n = B n, L], n =,, 3,, gives rise to the infinite number of partial differential equations PDEs of the KP hierarchy with dynamical variables {u i t, t, t 3, } with i =, 3, 4,. Here B n = n b n,i i L n denotes the differential part of L n and in following we will use L n L n B n to denote the integral part. The simplest nontrivial PDE constructed from is the KP equation given as x 4 u u u t 3 x 3 u x 3 i=0 3 u t = 0. 3 In Table we show the Lax operator and corresponding -dimensional examples of subhierarchies of the KP hierarchy. An alternative way to express the KP hierarchy is given by Date: Revision :., compiled February 5, 008.

2 JINGSONG HE, YI CHENG AND RUDOLF A. RÖMER the Zaharov-Shabat ZS equation ], B n t m B m t n B n, B m ] = 0, m, n =, 3, 4,. 4 The eigenfunction φ and the adjoint eigenfunction ψ of the KP hierarchy associated with equation 4 are defined by φ ψ = B n φ, = B t n t nψ, 5 n where φ = φλ; t and ψ = ψλ; t and t = t, t,. The n-reduction of the KP hierarchy corresponds to the situation L n = 0 such that Ln = B n = n v n n v v 0. Then the v i, i = 0,,, n, are independent of t n, t n, t 3n,. In this way the Lax pair of the -dimensional integrable system can be found. Well-nown examples of such n-reductions include the 4-reduction of the KP hierarchy 3] with Lax pair 4 x 4u x 4u x x u xx 4u vφ = λφ, 6 t φ = 3 x 3u x 3 u xφ, t = x, t 3 = t, 7 corresponding to the Satsuma-Hirota SH equation 3] 4u t uu x u xxx 3v x = 0, v t 6uv x v xxx = 0. 8 Furthermore, eliminating v in the above equations, we can obtain a 6th order equation u = z x 8z tt z xxxxxx z xxxt 8z x z xxxx 36z xx z xxx 7z xz xx = 0, 9 which has been called bi-directional Satsuma-Hirota equation 3]. Naturally, there also exist n-reductions of the BKP and CKP hierarchies. For example, the 5-reduction of the BKP hierarchy with u = u is given as 5x 5u 3x 5u x x 5u 03 u xx 53 z t x ] φ = λφ, 0 t φ = 3 x 3u x φ, u = z x, t 3 = t, t = x, which is the Lax pair corresponding to bi-directional Sawada-Kotera bsk equation 0,] zxxxxx 5z x z xxx 5z 3 x 5z x z t 5z xxt x 5z tt = 0. The 5-reduction of the CKP hierarchy u = u with Lax pair 0,] x 5 5u x 3 5 u x x 5u 35 6 u xx 5 3 z t x 5uu x 5 3 u xxx 5 ] 6 u t φ = λφ, 3 t φ = 3 x 3u x 3 u xφ, u = z x, t 3 = t, t = x, 4 gives the bi-directional Kaup-Kupershmidt equation z xxxxx 5z x z xxx 5zx 3 5z x z t 5z xxt 45 4 z xx 5z tt = 0. 5 x

3 SOLVING BI-DIRECTIONAL SOLITON EQUATIONS 3 An essential characteristic of the KP hierarchy is the existence of the τ-function and all dynamical variables {u i }, i =, 3,..., can be constructed from it,], e.g., u = log τ, 6 x u 3 = 3 x log τ, 7 3 x t. So it is a central tas to construct the τ-function in order to solve the nonlinear PDEs associated with the KP hierarchy. In the following, we will show that φ and ψ play a ey role in this construction. Gauge transformations 4, 5] offer an efficient route towards the construction of the τ function of the KP hierarchy. In Ref. 6] two inds of such a gauge transformation have been proposed, namely, T D φ = φ φ, T I ψ = ψ ψ 8 resulting in a very general and universal τ function see equation 3.7 of 6] and also IW,n in 7]. The determinant representation of the gauge transformation operators with n steps is given in Ref. 7]. In particular, the Grammian τ function 8] of the KP hierarchy can be generated by an iteration of the transformation 6, 9, 30]. This is straightforwardly understood from Chau s τ function and the determinant representation 7] if we impose a restriction on the generating functions of the gauge transformation. Grammian τ function have also been used to solve the reduction of the constrained BKP and CKP hierarchies 3 33,36]. There are two issues that arise when one wants to study the solutions of the - dimensional solitons equations given by the n-reduction of the BKP and CKP hierarchies. The first is how it retain the restrictions, i.e. L = L for BKP and L = L for CKP, for the transformed Lax operators L = TLT. In other words, the problem is how to obtain the τ-functions τ n BKP and τn CKP from the general τ-function τn = IW,n τ 0 with the gauge transformation T n of the KP hierarchy. Here τ 0 is the initial value of the τ-function of the KP hierarchy. Also, the generating functions φ i, ψ i of the gauge transformation will be complex-valued and related to the n-th roots of e iε. The second issue therefore is how to choose generating functions φ i = φλ i ; x, t and ψ i = ψµ i ; x, t such that τ n BKP to a physical τ-function ˆτ n Eq, which is real and positive on the full x, t plane. and τn CKP correspond In fact, the and bsk equations have been introduced recently by Dye and Parer 0,] when looing for the bidirectional soliton analogues of the Sawada-Kotera SK 6,7] and Kaup- Kupershmidt KK 8,9] equations. The Lax pairs of and bsk related similarly as the Lax pairs of KdV and Boussineq equation, thus ensuring their integrability. Both and bsk equation have a bidirectional soliton solution 0, ] which have been obtained by the Hirota bilinear method 37]. The profile of the solitons depend on their direction of propagation. The right-going solitons of are standard one-pea solitons, but the left-going solitons have two peas. Very recently, Verhoeven and Musette 3] have plotted the bi-directional solitons for the and equation based on the Grammian τ function. In this paper, we want to study why the 5-reduction of the BKP and CKP hierarchies have bidirectional soliton solutions, whereas their 3-reduction does not. As a first step, we will therefore exhibit the relationship between the periodic, left-going and right-going solitons of the 5-reduction and the 5-th roots of e iε. In order to do so, we derive the τ functions of the BKP

4 4 JINGSONG HE, YI CHENG AND RUDOLF A. RÖMER and CKP hierarchies in sections -4. The explicit formulas of the corresponding τ-functions for solitons as well as for the periodic solutions of bsk and are given and the two-pea soliton is discussed in detail. In section 5, we will prove that no two-pea solitons exist for the equation. The one-pea soliton has bi-directional motion and we also obtain the periodic and two-soliton solutions. In section 6, we will discuss the lower and higher-order reductions of BKP and CKP hierarchies and also the n =even-reductions of the KP hierarchy. We will show that the soliton of the j -reduction of BKP and CKP hierarchies can move along j directions j =,,..., investigate the relationship with the symmetric distribution of the j -th roots of e iε. In particular, we will obtain the stationary soliton for the higher reduction of the KP hierarchy. For the higher-order equation and even-reduction of KP hierarchy, we can again find a two-pea soliton.. τ functions for BKP and CKP hierarchies Let us first define the generalized Wronsian determinant IW,n IW,n g 0, g0,, g0 ; f0, f0,, f0 n 0 g f 0 0 g f 0 0 g f g f 0 n 0 g f0 0 g f0 0 g f0 0 3 g f0 n g f 0 0 g f 0 0 g f g f 0 n = f 0 f 0 f 0 3 f n 0 f 0,x f 0,x f 0 3,x f n,x f 0 n f 0 n f 0 3 n f n 0 n 9 Table. Examples of sub-hierarchies of the KP hierarchy, Lax operators used to construct them and resulting equations. The symbol indicates the conjugation, for example, =. There are some abbreviations used in Table: Sawada-Kotera SK, bi-directional Sawada-Kotera bsk, Kaup-Kupershmidt KK, bi-directional Kaup-Kupershmidt, Satsuma-Hirota SH, Yajima-Oiawa YO, Modified KdV MKdV, Non-linear Schrödinger NLS,constrained KPcKP. sub-hierarchy Lax operator example equation BKP,4] L = L SK 6,7], bsk 0,] CKP 4] L = L KK 8,9], 0,] n-thkdv 5] L n = 0 KdV ], Boussinesq-type ], SH 3] n =, 3, 4 ckp 7,] L = φ ψ YO 4], MKdV 5], NLS 6]

5 SOLVING BI-DIRECTIONAL SOLITON EQUATIONS 5 In particular, IW 0,n = W n f 0, f0,, f0 n with f 0 i = f 0 i x is the usual Wronsian determinant of functions {f 0, f 0,, f 0 n }. We shall also use the abbreviation f = fdx with integration constant equal to zero. Lemma. 6,7]. The τ function of the KP hierarchy generated by the gauge transformation T n is given as τ n = IW,n ψ 0, ψ0,, ψ0 ; φ0, φ0,, φ0 n τ0, 0 where φ 0 i, ψ 0 j = φλ i ; t, ψµ j ; t are solutions of equation 5 with initial value τ 0 for the τ-function and the initial values of the {u i } are {u 0 i }. Here {φ 0 i, ψ 0 j } are called generating functionsgfs of gauge transformation. Let us now discuss how to reduce the τ n in 0 to the τ function of the BKP hierarchy. The ey problem is how to eep the restriction L n = L n under the gauge transformation T n 7]. It should be noted that t = t, t 3, t 5, in BKP hierarchy. Proposition. see Refs. 30,39]. The Lax operator transforms as L n = T n LT n under the gauge transformation T n with n = and generating functions ψ 0 i = φ 0 i,x for i =,,, n. The τ function τ nn BKP of the BKP hierarchy is τ nn BKP = IW n,n φ 0 n,x, φ0 n,x,, φ0,x ; φ0, φ0,, φ0 n = 0 φ n,x φ 0 0 φ n,x φ 0. 0 φ,x φ 0 0 φ n,x φ 0 0 φ n,x φ 0. φ0 φ0 0 φ,x φ 0 0 φ n,x φ φ n,x φ 0 n 0 φ n,x φ 0 3 φ0.. 0 φ,x φ φ,x φ 0 3 φ0 n n 0 φ n,x φ 0 n 0 φ,x φ 0 0 φ,x φ 0 n n. 0 φ,x φ 0 n 0 φ,x φ 0 n τ 0 BKP Proof. It is clear that a single step of the gauge transformations T D or T I can not eep the restriction. So we use T T = T I ψ T D φ 0 3 such that the lax operator is L = TLT. Let us chec whether it satisfies the required restriction L = L 4 which means in terms of T that T D ψ T I φ 0 = T I ψ T D φ 0. 5

6 6 JINGSONG HE, YI CHENG AND RUDOLF A. RÖMER Based on the determinant representation of T 7] we see from 5 that φ 0 r.h.s = ψ 0 0 φ0 ψ0 0 φ φ l.h.s = ψ0 ψ 0 0 φ ψ 0 x φ 0 ψ0,x φ0 ψ 0 0 φ ψ 0, 6. 7 This implies ψ 0 = φ 0,x. So we have seen that in order to eep the restriction of the Lax operator, we have to regard T = T as basic building bloc in iteration of the gauge transformations T n. In particular, T = T I T D φ T I φ T D φ 0, T 33 = T I φ 3,x φ 5 3,x T D φ 4 3 T I,x φ 3,x T D φ T I φ,x T D and so on such that = n and ψ 0 i = φ 0 i,x for i =,,..., n. According to the determinant of T n 7] and τ n 6,7] with = n and ψ 0 i = φ 0 i,x, i =,,, n, τ nn BKP can be obtained directly from τn as in Lemma.. For the CKP hierarchy, we have again t = t, t 3, t 5,... and the restriction is L n = L n. Proposition. see Refs. 30,39]. The appropriate gauge transformation T n is given by n = and generating functions ψ 0 i = φ 0 i for i =,,..., n. The τ function τ nn CKP of the CKP hierarchy has the form φ 0, 8 τ nn CKP = IW n,n φ 0 n, φ0 n,, φ 0 ; φ 0, φ 0,, φ 0 n 9 = 0 φ n φ 0 0 φ n φ 0. 0 φ φ 0 0 φ φ 0 0 φ n φ 0 0 φ n φ 0. 0 φ φ 0 0 φ φ 0 0 φ n φ φ n φ φ φ φ φ φ n φ 0 n 0 φ n φ 0 n. 0 φ φ 0 0 φ φ 0 n n 0 φ n φ 0 n 0 φ n φ 0 n τ 0 CKP 30 0 φ φ 0 n 0 φ φ 0 n Proof. Similar to the BKP hierarchy we have to try the two-step gauge transformation T T = T I ψ T D. 3 φ 0 With L = TLT, the restriction L = L then implies T D ψ T I = T I T D φ 0 ψ φ 0. 3

7 SOLVING BI-DIRECTIONAL SOLITON EQUATIONS 7 Based on the determinant representation of T 7], we find from 3 that r.h.s = l.h.s = Then ψ 0 = φ 0. Again, we have to regard T = T I bloc such that T = T I T 33 = T I φ 3 φ 5 3 T D T D φ φ 4 3 T I T I φ0 ψ 0 0 φ ψ 0, 33 ψ0 φ 0 0 φ ψ0. 34 T D φ φ 3 so = n and ψ 0 i = φ 0 i for i =,,..., n. T D T D φ 0 φ φ, T I φ φ 0 T D φ 0 as basic building According to the determinant of T n 7] and τ n 6,7] with = n and ψ 0 i = φ 0 i, i =,,, n, τ nn CKP is obtained directly from τn in Lemma.. In fact, we can let ψ 0 i = c i φ 0 i,x or ψ0 i = c i φ 0 i with constants c i. However, the new τ nn BKP or τ nn CKP associated with ψ 0 i = c i φ 0 i,x or ψ0 i = c i φ 0 i is equivalent to the ones in Proposition. or Proposition.. Although Refs. 30, 39] have results similar to our Propositions. and., our approach is more direct and simpler for the construction τ nn BKP and τ nn CKP. If the initial values of dynamical variables {u i} of BKPCKP hierarchy are zero, then the equations in 5 of φ 0 i, ψ 0 j = φλ i ; t, ψµ j ; t become more simpler as, 35 φλ; t t n = n x φλ; t], t= t, t 3, t 5,, 36 ψµ; t t n = n n x ψµ; t], t= t, t 3, t 5,, 37 and τ 0 BKP = τ0 CKP =. Last, we note that for the generalized KP gkp hierarchy with Lax operator ˆL = L n, n =, 4, 6, 8,..., and ˆL = ˆL, the τ function τ n gkp generated by gauge transformations T n has the same form as for the CKP hierarchy. This result will afford a simple way to construct the τ function of equation in Section Soliton solutions of the bsk equation As pointed out in the introduction, there are two steps en route from a τ function τ n generated by the gauge transformations T n of the KP hierarchy to the τ function of equations as the n-reduction of BKP or CKP hierarchies. The second step is to build physical τ functions from the complex-valued τ nn BKP and τ nn CKP constructed in the last section. In the following Sections, we will illustrate our approach by computing the τ function for the 5-reduction of BKP and CKP, i.e., for the bsk and equations.

8 8 JINGSONG HE, YI CHENG AND RUDOLF A. RÖMER The 5-reduction of the BKP hierarchy is the bsk equation. Assume for the initial value u = 0 in equations 0 and, then φ 0 i = φλ i ; x, t are solutions of 5 x φλ i; x, t = λ i φλ i ; x, t, φλ i ; x, t t = 3 x φλ i; x, t. 38 So proposition. with τ 0 BKP = implies that the τ function of bsk is given as follows. Proposition 3.. The τ function of the bsk equation generated by T nn from initial value is = IW n,n φ 0 n,x, φ0 n,x,, φ0,x ; φ0, φ0,, φ0 n 39 τ nn bsk and the solution of the bsk equation generated by T nn from initial value 0 is u = x log τ nn bsk Here φ 0 i = φλ i ; x, t are solutions of equation 38. In general, this τ function τ nn bsk for bsk is complex and related to 5-th roots of e iε. We have to find the real and non-zero τ function from it such that u in equation 40 is a real and smooth solution of bsk. This is main tas of this section. We start by analysing the solution φλ; x, t of equation 38and mae the universal ansatz φλ; x, t = 40 5 A j e xp jtp 3 j, with p 5 j = λ. 4 j= Here p j = exp επj i, 5 = λ, R, 0 ε < π and j = 0,,, 3, 4. There are two 5 important ingredients which we can use to find the desired solution. The first is that the 5- th roots ε j = exp επj i of e iε are distributed uniformly on the unit circle in C. So for a 5 suitable value of ε there exist combinations of p j s which are symmetric upon reflection on the x-axes; similarly for the y-axes for other values of ε. The second ingredient is that τ bsk and exp αx βtτ bsk will imply the same solution u since u = x log τ bsk. Here, α and β are arbitrary, complex constants. Therefore we can obtain the desired real and smooth solutions of the bsk if τ bsk can be expressed as τ bsk = e αxβtˆτ bsk ˆτ bsk, in which ˆτ bsk is a real and nonzero function although τ bsk is complex. We call ˆτ bsk the physical τ function for the bsk equation. Based on the above arguments, let us mae the refined ansatz or φλ ; x, t = A e p xp 3 t B e q xq 3 t, p = e iε, q = e iε, 5 = λ, R, 4 φλ ; x, t = A e p xp 3 t B e q xq 3 t, p = e iε, q = e iε, 5 = λ, R, 43 and in the next step we need to fix the ratio B A. We stress that the above analysis is also true for the derivation of the equation. Proposition 3.. Define ξ = x cosε t 3 cos 3ε. Then the physical τ function of bsk generated by T is ˆτ bsk = B eξ e ξ B 44 A A

9 SOLVING BI-DIRECTIONAL SOLITON EQUATIONS 9 and the corresponding one soliton u = x log ˆτ bsk is given as u = B 6 A cos ε e ]. 45 ξ B A e ξ B A Here B A > 0. The velocity of the moving soliton is v = cos 3ε cos ε and can be both positive and negative depending on the choice of ε. Specifically, we have v = v ε = π < 0 and v = 0 v ε = 3π > 0. 0 Proof. τ bsk = φ 0 = A e ix sin ε t 3 sin3ε e ξ B A e ξ B A ] 46 and φ 0 = φλ ; x, t defined by equation 4. Let us point out a relation between the distribution of the 5-th roots of e iε and the direction of movement for the soliton. e iε ε =π/0, e iε ε =π/0 one distribution of 5-th roots of e iε p = e iε ε =π/0, q = e iε ε =π/0 in equation 4 v ε =π/0 < 0, left-going soliton u in equation 45 ; e iε ε =3π/0, e iε ε =3π/0 another distribution of 5-th roots of e iε p = e iε ε =3π/0, q = e iε ε =3π/0 in equation 4 v ε =3π/0 > 0, right-going soliton u in equation 45. We can see from equation 45 that the one-soliton of bsk has only one pea in its profile. The process of generating a two-soliton by T is more complicated.

10 0 JINGSONG HE, YI CHENG AND RUDOLF A. RÖMER Lemma 3.. With φλ ; x, t and φλ ; x, t as in equation 4, and using Proposition 3., is given by τ bsk τ bsk = A A eix sin ε sinε t 3 sin3ε 3 sin 3ε ] { 4z f e ξ ξ 4 cosε ε ] 4z f e ξξ B B 4 cosε ε ] A A 4z 3 f 3 e ξ ξ B 4z3 f 3 e ξ ξ 4 cosε ε ] A 4 cosε ε ] iz 5 f 5 e ξ cosε ε ] cosε ε ] iz 5 f 5e ξ cosε ε ] cosε ε ] iz 7 f 5 e ξ cosε ε ] cosε ε ] iz 7 f 5e ξ cosε ε ] cosε ε ] cosε ε ] cosε ε ] B A B A B A B A B A B A B A B A B } 47 Here the z i, i =, 3, 5, 7 are given in A and f = cosε ε ], f 3 = cosε ε ] and f 5 = f f 3, as well as ξ i = x i cosε i t 3 i cos 3ε i for i = and. We now need to find a suitable solution of B i A i, i =, such that the summation of terms in the {} bracet of equation 47 is a positive function on the whole x,t plane. The following two lemmas is useful. Let z = 4z f, z 3 = 4z 3 f 3, z 5 = iz 5 f 5, z 7 = iz 7 f 5, and z i, i =, 3, 5, 7 given in A. Lemma 3.. For z i, i =, 3, 5, 7 there exist relations Lemma 3.3. Let B hold. z 3 z z z z 5 z A = z, B z 7 A = z z 5 B A B A B A z z 3 = z 5 ; z z 3 = z 7. 48, and g = z, g z 3 6 = g 8 = z, g z 5 9 = =, z 3 B B = g, A B = g 6, A z 7 z z B z A B A A A then z 3 =, 49 B = g 9, 50 A B = g 8, 5 The B i A i, i =, in Lemma 3.3 are what we are looing for, and then the physical τ function is give by following proposition. ˆτ bsk A

11 SOLVING BI-DIRECTIONAL SOLITON EQUATIONS Proposition 3.3. ˆτ bsk = { e ξ ξ 4 cosε ε g e ξξ 4 cosε ε e ξ ξ 4 cosε ε e ξ ξ 4 cosε ε e ξ cosε ε cosε ε g 6 e ξ cosε ε cosε ε e ξ cosε ε cosε ε g 8 e ξ cosε ε cosε ε g } 9 cosε ε cosε ε, 5 the two soliton solution is u = x log ˆτ bsk. In particular, ε = ε = π 0 overtaing solitons and moving in negative direction; ε = ε = 3π 0 solitons and moving in positive direction; ε = π 0, ε = 3π 0 results in two produces two overtaing results in two head-on solitons. We have plotted the one solitons in figure associated with parameters A = B =, = of equation 45. The two solitons in Proposition 3.3 are shown in figure associated with parameters = and =.3. Figure. Single left- and right-going solitons for the bsk equation : ε = π 0 left, ε = 3π 0 right.

12 JINGSONG HE, YI CHENG AND RUDOLF A. RÖMER Figure. Two left- and right-going as well as head-on colliding solitons for the bsk equation : ε = ε = π left, ε 0 = ε = 3π right and ε 0 = π, ε 0 = 3π 0 collision. 4. Periodic and soliton solutions of equation The 5-reduction of the CKP hierarchy yields the equation equation 5. Let the initial value be u = 0 in Eqs. 3 and 4, then φ 0 i = φλ i ; x, t are solutions of 5 x φλ i; x, t = λ i φλ i ; x, t, φλ i ; x, t t = 3 x φλ i; x, t. 53 So the Proposition. implies the τ function of equation.

13 SOLVING BI-DIRECTIONAL SOLITON EQUATIONS 3 Proposition 4.. The τ function of the equation generated by T nn from initial value is τ nn = IW n,nφ 0 n, φ 0 n,, φ 0 ; φ 0, φ 0,, φ 0 n 0 φ n φ 0 0 φ n φ 0 0 φ n φ φ n φ 0 n 0 φ n φ 0 0 φ n φ 0 0 φ n φ φ n φ 0 = φ φ 0 0 φ φ 0 0 φ φ φ φ 0 n 0 φ φ 0 0 φ φ 0 0 φ φ φ φ 0 n n 0 φ n φ 0 n 0 φ n φ 0 n 0 φ φ 0 n 0 φ φ 0 n and the solution u of the from initial value zero is u = x log τ nn. 55 Here φ 0 i = φλ i ; x, t are solutions of equation 53. As before, τ nn is complex and related to the 5-th roots of eiε and again we have to find a physical τ function ˆτ nn such that u in equation 55 is real and smooth solution includes solitons and periodic solutions. The case of n = and n = will be discussed in detail. Similar to the bsk equation, we should assume the solutions of equation 53 as or 54 φλ ; x, t = A e p xp 3 t B e q xq 3 t, p = e iε, q = e iε, 5 = λ, R, 56 φλ ; x, t = A e p xp 3 t B e q xq 3 t, p = e iε, q = e iε, 5 = λ, R, 57 to extract physical τ function ˆτ nn from τnn. At first, we would lie to give the two simple cases which are generated by the gauge transformations T. Proposition 4.. Let ξ = x cosε t 3 cos 3ε, B A = ie iε and φ 0 = φλ ; x, t defined by equation 56, then the physical τ function of extracted from τ nn is n= ˆτ = eξ e ξ 58 sin ε and the corresponding one soliton u = x log ˆτ is u = 4 cos ε cosh ξ sin ε 59 cosh ξ sinε with ε = π 3π or. The velocity of the soliton is v = 0 0 cos 3ε cos ε. In particular, the left-going soliton have two peas in its profile and the negative speed v = v ε = the right-going soliton 0; π have only one pea and positive speed v = v ε =. 3π 0

14 4 JINGSONG HE, YI CHENG AND RUDOLF A. RÖMER Proof. Taing φ 0 = φλ ; x, t of equation 56 and n = bac into Proposition 4., the straightforward calculation leads to τ = φ 0 A = eix sinε t 3 sin3ε p Here ξ = x cosε t 3 cos 3ε. e ξ e ξ ]. 60 sin ε If φ 0 = φλ ; x, t defined by equation 57, then we can get a periodic solution as following proposition. Proposition 4.3. Let η = x sin ε t 3 sin 3ε, and φ 0 = φλ ; x, t defined by equation 57, A = B = in φ 0, then the physical τ function of extracted from τnn is n= ˆτ = cosε cosη ε 6 and the corresponding solution 4 u = x log ˆτ sin ε cosη ε cos ε = 6 cos ε cosη ε is periodic. Here ε = π 4π or. The velocity for the solution is v = 0 0 sin3ε sinε. If ε = π, u 0 in equation 6 is a left-going periodic wave. If ε = 4π, u in equation 6 is a right-going 0 periodic wave. Proof. τ = φ 0 e x cos ε ] t 3 cos 3ε = cosη ε cosε Here η = x sin ε t 3 sin 3ε, and φ 0 = φλ ; x, t is defined by equation 57. Based on Propositions 4. and 4.3, we can find following corresponding relationship between symmetrical distributions of 5-th roots of e iε and moving direction of solutions. e iε, e iε ε = π the first distribution of 5-th roots of eiε p 0 = e iε, q = e iε ε = π in equation 56 left-going two-pea soliton in equation 59; 0 e iε, e iε ε = 3π the second distribution of 5-th roots of e iε p = e iε, q = 0 e iε ε = 3π in equation 56 right-going one-pea soliton in equation 59; 0 3 e iε, e iε ε = π the third distribution of 5-th roots of e iε p = e iε, q = 0 e iε ε = π in equation 57 left-going periodic wave in equation 6; 0 4 e iε, e iε ε = 4π the fourth distribution of 5-th roots of e iε p = e iε, q = 0 e iε ε = 4π in equation 57 right-going periodic wave in equation 6. 0 There are only four distributions of 5-th roots of e iε, which are symmetric respect with x-axes or y-axes. However, there exist several other pairs of roots in the above four distributions which will result in divergent solutions of through the above procedure. For example, p = e i3π 0, q = e i3π 0 or p = e i8π 0, q = e i8π 0. Let us now concentrate on the two-pea soliton solution in equation 59.

15 SOLVING BI-DIRECTIONAL SOLITON EQUATIONS 5 Lemma 4.. Let x, constant a > 0 and function y = yx = x a x, 64 a then if a > /, y < 0; x if a = /, then y x x= = 0; 3 if > a > 0, then there exists one point x > such that y x x=x = 0, x is one extreme maximum point of y, and y x= > 0. x Proof. We have y x = y x = a x a a x a Firstly, y x < 0 if a > /. Secondly, if a = /, y x = 0 when x =. At last, if / > a > 0, there exist x > such that y x = 0. Note that y x > 0 if x, x, y x < 0 if x > x. So x is one extreme maximum point of y. Proposition 4.4. Let a, b, be positive constants, ξ = x ct, c R, for following ind of solution u = b cosh ξ a cosh ξ, 66 a if a /, u has one pea in its profile defined by ξ = 0; if 0 < a <, then there exist two peas in profile; 3 There exist no more than two peas in a soliton give by equation 66. Proof. By calculation, we have u x = b sinh ξ a cosh ξ a a cosh ξ a According to the Lemma 4., we have a > /, there exist ξ = 0 such that u x = 0 because sinh ξ ξ=0 = 0. Note a cosh ξ a < 0. a = /, there exist ξ = 0 such that u x = 0 because sinh ξ ξ=0 = 0 and a cosh ξ a ξ=0 = 0. However, let ξ be sufficiently small, we have u x < 0 if ξ > 0 and u x > 0 if ξ < 0. So ξ = x ct = 0 defines one extreme maximum line of ux, t on x,t plane. 3 / > a > 0, there exist ξ = 0 and ξ > 0 and ξ = ξ < 0 such that u x = 0. But u x > 0 if ξ < ξ ; u x < 0 if ξ ξ, 0; u x > 0 if ξ 0, ξ ; u x < 0 if ξ > ξ. So ξ = x ct = 0 defines one extreme minimal line on x,t plane; 0 < ξ = x ct and 0 > ξ = x ct define two extreme maximum lines on the x,t plane. Using u 0 if ξ, conclusions are proven. Comparing equation 66 with equation 59 we get a = sin ε, and then can understand why sin ε ε =π/0 will lead to two peas in one soliton of but sinε ε =3π/0 will lead only to one pea in one soliton of. On the other hand, one soliton solution of u in equation 59 have one pea or two peasmaximum case in its profile. According to analysis above, we can claim from the point of view of reduction in KP hierarchy that the existence of two peas in the soliton is traced to three facts:

16 6 JINGSONG HE, YI CHENG AND RUDOLF A. RÖMER The Grammian τ function in Proposition 4. which determines the form of soliton in equation 59; The order of n-reduction, i.e. n 5 can produce two peas soliton in KP hierarchy; 3 The phase ε of n-th rootn 5 of e iε, such that 0 < a = sin ε < /. Now we turn to the more complicated τ from Proposition 4., which generates the two soliton and periodic solution with two spectral parameters of equation. The first case is the two soliton solution. Lemma 4.. Let φ 0 i = φλ i ; x, t, i =,, defined by Eq.56, ξ i = x i cosε i ti 3 cos 3ε i, η i = x i sin ε i ti 3 sin 3ε i, i =,, then τ nn n= gives out τ = A A eiη η { z eξ ξ 4 cosε ε z e ξξ 4 cosε ε B B A A z 3e ξ ξ 4i cosε ε B z 3 e ξ ξ A 4i cosε ε z e ξ i sin ε cosε ε cosε ε z e ξ i sin ε cosε ε cosε ε z 4 eξ i sin ε cosε ε cosε ε z 4 e ξ i sin ε cosε ε cosε ε B A B A B A B 4 cos ε cos ε sin ε sin ε sin ε sin ε cosε ε cosε ε Here z i, i =,, 3, 4, are given in B. z i means the complex conjugation of z i. B A B A A B A B A B } In order to extract physical τ function ˆτ, we need following two Lemmas for suitable B i A i, i =,. Lemma 4.3. For z i, i =,, 3, 4, as given in B, the following identities A 68 z = z z 3, z 4 = z z 3 69 hold.

17 Lemma 4.4. Let B A = i z SOLVING BI-DIRECTIONAL SOLITON EQUATIONS 7, B z4 A = i z z, and g 5 = z 3, g 6 = g 8 = z z, g 9 = z 4 z 4, then z 3 B z A z B z A z B A iz =, z 3 B = g 9, A B = g 8, A z z 4 iz B A B A z B A =, 70 B = g 5, 7 A B = g 6 7 A hold. Taing B i A i, i =,, and relations in Lemma 4.4 bac into Lemma 4., the physical τ function is obtained. ˆτ Proposition 4.5. ˆτ = { e ξ ξ 4 cosε ε g 9 e ξξ 4 cosε ε e ξ ξ 4 cosε ε e ξ ξ 4 cosε ε e ξ sin ε cosε ε cosε ε g 8 e ξ sin ε cosε ε cosε ε e ξ sin ε cosε ε cosε ε g 6 e ξ sin ε cosε ε cosε ε g 5 4 cos ε cos ε sin ε sin ε sin ε sin ε cosε ε cosε ε } 73 The two solitons solution is u = x log ˆτ. In particular, ε = ε = π results in two 0 overtaing soltions moving in negative direction; ε = ε = 3π results in two overtaing soltions 0 moving in positive direction; ε = π, ε 0 = 3π results in head-on colliding two soltions. 0 The second case is a periodic solution with two spectral parameters of equation from Proposition 4..

18 8 JINGSONG HE, YI CHENG AND RUDOLF A. RÖMER Lemma 4.5. Let φ 0 i = φλ i ; x, t, i =,, defined by Eq. 57, ξ i = x i cosε i ti 3 cos 3ε i, η i = x i sin ε i ti 3 sin 3ε ii =,, then τ nn gives n= τ = A A e iξ ξ { ze iη η 4 cosε ε z e iηη 4 cosε ε B B A A z3 eiη η 4 cosε ε B z 3 e iη η A 4 cosε ε z eiη cosε cosε ε cosε ε z e iη cosε cosε ε cosε ε z 4 eiη cosε cosε ε cosε ε z 4 e iη cosε cosε ε cosε ε 4 cos ε cos ε sin ε sin ε B A B A B A B cosε cosε cosε ε cosε ε Here z i, i =,, 3, 4, are given in C. z i indicates the complex conjugation of z i. B A B A A B A B A B } Similar to the two solitons solution of equation, we need following two Lemmas to find suitable B i A i, i =,, to extract physical ˆτ from equation 74 for periodic solution. Lemma 4.6. For z i, i =,, 3, 4, as given in C, the following identities hold. A 74 z = z z 3, z 4 = z z 3 75 Lemma 4.7. Let z = z e iθ, =,, 3, 4, be as given in C, and B A = e iθ, B A = e iθ 4, and g = z, g z 3 = z 3, g z 4 = z 4, g z 5 =, then z z 3 B = z 3 B = g 3, 76 hold. z A z B z A z z B A = z z z A B =, A B A z B A B = g, A We can get the physical τ function ˆτ bac into Lemma 4.5. B = g 5, 77 z4 z A B A = z 4 z B A B = g 4 78 A by taing B i A i, i =, and relations in Lemma 4.7

19 Proposition 4.6. SOLVING BI-DIRECTIONAL SOLITON EQUATIONS 9 ˆτ = { cos η η 4 cosε ε g 3 cos η η 4 cosε ε g cos η cosε cosε ε cosε ε g 4 cos η cosε cosε ε cosε ε g 5 4 cos ε cos ε sin ε sin ε cosε cosε cosε ε cosε ε } The periodic solution with two parameters and is u = x log ˆτ 79. Furthermore, ε = ε = π results in two overtaing waves moving in negative direction ; ε 0 = ε = 4π results 0 in two overtaing waves moving in positive direction; ε = π, ε 0 = 4π results in two head-on 0 colliding waves. We have plotted soliton solutions of in figure 3, and there periodic solutions with two spectral parameters in Fig Periodic and soliton solutions of equation The τ function of the equation is still in the form of a Grammian although the equation does not belong to the CKP hierarchy, which is obtained in 38] through the Bäclund transformation. Similar to the equation, its τ function is in the form of Grammian, we can find τ function τ and τ u = 0 in equations 6 and 7, then φ 0 i 4 x φλ i; x, t = λ i φλ i ; x, t, of from Grammian τ function. Let the initial value be = φλ i ; x, t are solutions of φλ i ; x, t t = 3 x φλ i; x, t. 80 Proposition 5. see Ref. 38]. The τ function of equation generated by Bäclund transformation from initial value u = 0 is τ nn = IW n,n φ 0 n, φ0 n,, φ0 ; φ0, φ0,, φ0 n 0 φ n φ 0 0 φ n φ 0 0 φ n φ φ n φ 0 n 0 φ n φ 0 0 φ n φ 0 0 φ n φ φ n φ 0 = φ φ 0 0 φ φ 0 0 φ φ φ φ 0 n 0 φ φ 0 0 φ φ 0 0 φ φ φ φ 0 n and the solution u of from initial value zero is u = x log τnn Here φ 0 i = φλ i ; x, t are solutions of equation 80. n 0 φ n φ 0 n 0 φ n φ 0 n 0 φ φ 0 n 0 φ φ 0 n 8 8

20 0 JINGSONG HE, YI CHENG AND RUDOLF A. RÖMER Figure 3. Soliton solutions of the equation 5. Top left: one left-going two-pea soliton when ε = π and 0 =.. Top right: two left-going soliton when =, =.3, ε = ε = π. Bottom: Head-on collision of left- and 0 right-going solitons when =.8, =.3, ε = 3π, ε 0 = π. 0 In fact, τ nn can be generated by gauge transformation T n n=. The Lax pair of is L = 4 x 4u x 4u x x u xx 4u v, M = 3 x 3u x 3 u x, and satisfy L = L, M = M. Similar to the CKP hierarchy, let T = T = T I ψ T D φ 0, and do gauge transformation L = TL T. So L = L requires T D ψ T Iφ 0 = T Iψ T Dφ 0 as we have seen in CKP hierarchy. The remaining procedure is the same as the gauge transformation of the CKP hierarchy as well as the equation. Of course, the generating functions φ 0 i, ψ 0 i = φλ i ; x, t, ψλ i ; x, t satisfy equation 6 and equation 7 if the initial values are u 0, v 0, or equation 80 if the initial values are u = 0, v = 0. Remar 5.. We should note that L v=0 = x 4 4u x 4u x x u xx 4u = x u = L KdV. The Lax pair of the KdV equation is

21 SOLVING BI-DIRECTIONAL SOLITON EQUATIONS Figure 4. Periodic solutions with two spectral parameters of equation 5. Top left: left-going periodic solution with = 0., = 0.3, ε = ε = π. 0 Top right: right-going periodic solution when = 0.4, = 0.5, ε = ε = 4π. 0 Bottom: Collision of left- and right-going periodic solution when =, =.5, ε = π, ε 0 = 4π. 0 L KdV = x u, M KdV = 3 x 3u x 3 u x. T D φ 0 generates a single soliton solution of the KdV from zero initial value. Here φ 0 = φλ ; x, t satisfy L KdV φλ ; x, t = λ φλ ; x, t and φλ ;x,t = M t KdV φλ ; x, t simultaneously. The left-going multi-soliton can be produced by using repeated iteration of T D. In order to get real and smooth solutions, such as soliton and periodic solution, we should construct physical τ function ˆτ from τ nn which is complex and related to 4-th roots of e iε. The case of n = and n = will be discussed in detail. Let us start to discuss the single soliton with two directional propagation. To do this, similar to the above two sections, we should assume the solution of equation 80 as φλ ; x, t = A e p xp 3 t B e q xq 3 t, p = e iε, q = e iε, 4 = λ, R, 83

22 JINGSONG HE, YI CHENG AND RUDOLF A. RÖMER or φλ ; x, t = A e p xp 3 t B e q xq 3 t, p = e iε, q = e iε, 4 = λ, R. 84 For φλ i ; x, t, i =,, the difference between here and above two sections is the i 4 = λ i, i =,, instead of i 5 = λ i, i =,. From Proposition 5. we can extract physical τ function from τnn. n= ˆτ Proposition 5.. Let ξ = x cosε t 3 cos 3ε, B A = ie ε, and φ 0 = φλ ; x, t as defined by equation 83, then the physical τ function of extracted from τ nn n= is ˆτ = eξ e ξ 85 sin ε and the corresponding single soliton u = u = x log ˆτ is 4cosε cosh ξ sinε. 86 cosh ξ sin ε Here ε = π 4. The velocity of the soliton is v = cos 3ε cos ε ε = π 4 > 0. Proof. τ = φ 0 = A e ix sin ε t3 p sin 3ε e ξ e ξ sin ε As we discussed in Remar 5., the left-going soliton can also be generated by T D. Proposition 5.3. Let ξ = x cosε t 3 cos 3ε ε =0, and φ 0 = φλ ; x, t ε =0 as defined by equation 83, then the physcial τ function of generated by T D φ 0 is and the corresponding single soliton u = ˆτ = A u = x log ˆτ Here A B > 0. The velocity of the soliton is v = < 0. Proof. 87 e ξ 88 B is 4 A B e ξ A B e ξ. 89 τ = φ0 = B e ξ A e ξ B It can be clarified by u = 0, v = 0 T Dφ 0 u 0, v = 0, and then τ = φ0. On the other hand, if φ 0 = φλ ; x, t as defined by equation 84, then we can get periodic solution from Proposition

23 SOLVING BI-DIRECTIONAL SOLITON EQUATIONS 3 Proposition 5.4. Let η = x sin ε t 3 sin 3ε, A = B = in φ 0, then the physical τ function of equation for periodic solution extracted from τ nn n= is ˆτ = cosη ε 9 cosε and the corresponding peroidic solution u = x log ˆτ is 4 sin ε cosη ε cos ε u =. 9 cos ε cosη ε Here ε = π 4. The velocity of the solution is v = sin 3ε sinε ε = π 4 < 0. Proof. τ = φ 0 = ex cos ε t3 4 cos 3ε cosη ε. 93 cosε There are some relationship between the distributions of 4-th roots of e iε and moving direction of solutions. e iε, e iε ε =0 the first distribution of 4-th roots of e iε p = e iε, q = e iε ε =0 in equation 83 left-going soliton in equation 89; e iε, e iε ε = π the second distribution of 4-th roots of eiε p 4 = e iε, q = e iε ε = π in equation 83 right-going soliton in equation 86; 4 3 e iε, e iε ε = π the third distribution of 4-th roots of e iε p 4 = e iε, q = e iε ε = π in equation 84 left-going periodic wave in equation 9. 4 In the above discussion, we now the right-going soliton and left-going periodic wave of of have the completely same form with the equation, except ε = π/4 instead of ε = π/0 and ε = 3π/0. The reason is that the τ function of two equations is in the same Grammian of generating functions φ 0 i, and generating functions φ 0 i for two equations satisfy analogous linear partial differential equations with constant coefficients, i.e. equation 53 for equation, equation 80 for equation. These relations between and are still true for their two soliton and two parameters periodic solutions. Proposition 5.5. The two right-going solitons are given u = x log ˆτ in which ˆτ is ˆτ = ˆτ ε =ε =π/4 95 and ˆτ is given by Proposition 4.5. Proposition 5.6. The right-going periodic wave with two spectral parameters and is given by u = x 96 log ˆτ 94

24 4 JINGSONG HE, YI CHENG AND RUDOLF A. RÖMER in which ˆτ is ˆτ = ˆτ ε =ε =π/4 97 and ˆτ is given by Proposition 4.6. According to the analysis in Remar 5., two left-going solitons of equation can be generated by a chain of gauge transformations u = 0, v = 0 T Dφ 0 u 0, v = 0 T D φ u 0, v = 0 using the notation of 7], φ 0 i = φλ i ; x, t ε =0, i =,, are defined by equation 83. Their τ function of generated by T = T D φ T D φ 0 is τ = φ 0 φ 0 φ 0,x φ 0,x. 98 From τ we can obtain the physical τ function ˆτ and two soliton solution. Proposition 5.7. Let φ 0 i = φλ i ; x, t ε =0 are defined by equation 83, ξ i = i xi 3 t, i =,. If B A > 0, B A < 0, >, then the physical τ function ˆτ is given by The two soliton solution is u = ˆτ = e ξ ξ B A B A e ξ ξ B A e ξ ξ B A e ξ ξ. 99 x log ˆτ, which is left-going. The collision of two soliton is generated by gauge transformation chain u = 0, v = 0 T Dφ 0 u 0, v = 0 T I ψ T Dφ u 0, v 0, φ 0 = φλ ; x, t ε =0 is defined by equation 83, ψ 0 = φ 0, φ 0 = φλ ; x, t is defined by equation 83. The corresponding τ function of is τ = 0 ψ φ0 ψ 0 φ0 φ 0 φ 0 ψ 0 =φ0 === 0 φ φ0 φ 0 φ0 φ 0 φ Taing φ 0 i, i =,, bac into equation 00, we have its explicit expression as following Lemma.

25 SOLVING BI-DIRECTIONAL SOLITON EQUATIONS 5 Lemma 5.. Let ξ = x t 3, ξ = x cosε t 3 cos 3ε, η = x sin ε t 3 sin 3ε, z i = c i d i, i =, 3, 5, as given in D. τ = eiη A A { z eξ ξ cosε z 3 e ξ ξ cosε z e ξ ξ B B cosε A A z 3 e ξ ξ B A cosε A B z 5 e ξ i sin ε cosε cosε z 5 e ξ i sin ε cosε cosε Lemma 5.. For z i = z i e θ i i =, 3, 5, the following identities are true. Lemma 5.3. Let z i, i =, 3, 5, be as given by D, if B then g = z B B, g z 4 = z 3 B A A z A hold. B A B B A A }. 0 z z 5 = z 3 z 5, e iθ 3 θ = e iθ 3 0 A = z z3, B A = i z z5, g = g 4 = z, g 6 = z 5 iz B B A A z 5, g 6 =,, 03 With the help of Lemmata 5. and 5.3, we deduce the physical τ function of colliding two soliton of equation from Lemma 5.. Proposition 5.8. Let g, g 4 be given as in Lemma 5.3, then { ˆτ = e ξ ξ g e ξ ξ cosε cosε e ξ ξ g 4 e ξ ξ cosε cosε e ξ sin ε cosε cosε e ξ sin ε cosε cosε } 04 We have plotted the two soliton solutions of equation in figure 5, and periodic solutions with one spectral parameter and with two spectral parameters of the same equation in figure Lower and Higher order reductions In this section, we want to discuss the general character of soliton equation from lower order to higher order in one same sub-hierarchy. The purpose is to show the relation between propagation of soliton on x,t plane and the order of Lax pair, and show the difference between

26 6 JINGSONG HE, YI CHENG AND RUDOLF A. RÖMER Figure 5. Two left- and right-going as well as head-on colliding solitons for the equation 9. Parameters are chosen as: A = B, A =, B =, =.5, = left; =.5, =.3, ε = ε = π 4 right; = 0.8, = 0.9, ε = π 4 collision. the lower reduction and higher reduction. Let Lax pair of soliton equation is L, M, which defines φλ; x, t by φλ; x, t Lφλ; x, t = λφλ; x, t, = Mφλ; x, t. 05 t There are some examples of n-reduction of the KP hierarchy. For the BKP hierarchy, Here Lax pair 3-reduction 5-reduction 7-reduction 9-reduction L B 3 B 5 B 7 B 9 M B5 B 3 B 3 B 3 Equation SK bsk higher order higher order. 06 B 5 = x 5 5u 3 x 5u x x 5u 0 3 u xx x, 07 SK 6,7] : 9u t 45u u x u xxxxx 5uu xxx 5u x u xx = 0. 08

27 SOLVING BI-DIRECTIONAL SOLITON EQUATIONS 7 Figure 6. Left-going periodic solutions with one left and two right spectral parameters for the equation 9. Parameters: =, ε = π 4 left; = 0., = 0.3, ε = ε = π 4 right. B 5 and B 3 are given by equations 0 and. For the CKP hierarchy, Lax pair 3-reduction 5-reduction 7-reduction 9-reduction L B 3 B 5 B 7 B 9 M B5 B 3 B 3 B 3 Equation KK higher order higher order Lax pair -reduction 4-reduction 6-reduction 8-reduction L B B 4 B 6 B 8 M B 3 B 3 B 3 B 3 Equation KdV higher order higher order. 09 Here B 5 = x 5 5u 3 x 5 u x x 5u 35 6 u xx x 5uu x 5 3 u xxx, 0 KK 8,9] : 9u t 45u u x u xxxxx 5uu xxx 75 u xu xx = 0, B 5 and B 3 are given in equations 3 and 4. There are several even-reductions of the KP hierarchy as following, Now we start to discuss the BKP hierarchy. Lemma 6.. Let ξ = x cosε t 5 cos 5ε, then ˆτ SK ˆτ and the corresponding single soliton is u = 4 cos 5ε cos ε ε = π 6 = 4 > 0. Here ˆτ bsk is expressed by. SK = ˆτ bsk ξ > ξ 3 x. The velocity of soliton v = log ˆτ SK is given by Proposition 3.. Proof. Because the SK equation and bsk equation belong to the same sub-hierarchy BKP, so the results of bsk are also hold by SK equation only if we replace ξ in bsk by ξ =

28 8 JINGSONG HE, YI CHENG AND RUDOLF A. RÖMER x cosε t 5 cos 5ε. For SK equation, the generating functions φ 0 i transformation satisfy 3 x φλ; x, t = λφλ; x, t, φλ; x, t t = φλ i ; x, t of gauge = x 5 φλ; x, t, 4 which are different with equation 38 for bsk equation. So 3 = λ. This difference determines replacement in equation 3. Of course, similar to the bsk, we also should assume the solutions of equation 4 be the form of or φλ ; x, t = A e p xp 5 t B e q xq 5 t, p = e iε, q = e iε, 3 = λ, R, 5 φλ ; x, t = A e p xp 5 t B e q xq 5 t, p = e iε, q = e iε, 3 = λ, R. 6 Taing the generating functions φ 0 i in Eq. 5 bac into the Proposition., then we can extract ˆτ SK from τ SK. The relation ˆτ SK = ˆτ bsk ξ > ξ is given by comparison. In particular, there are two distributions of roots of third-order of e iε on circle, which is symmetric with respect to y-axes. However, they are corresponding to same single soliton solution. e i π 6, e i π 6 one distribution of 3-order root of e iε on unit circle p = e i π 6, q = e i π 6 in equation 5 a single soliton in Lemma 6.. e iπ 6, e i π 6 one distribution of 3-order root of e iε on unit circle p = e iπ 6, q = e iπ 6 in equation 5 one soliton as. Lemma 6.. The higher order equations of the BKP hierarchy are defined by equation 06. For the n-reduction equation of the BKP hierarchy nbkp, n = j, j = 3, 4, 5,, and let ξ mp = x m cosε p tm 3 cos 3ε p, m n = j m = λ m, then the physical τ function of the nbkp generated by T is ˆτ nbkp = ˆτ bsk ξ > ξ p, 7 and the corresponding single soliton of the nbkp is u = x log ˆτ nbkp. Here ε p = p π = n p π, p =,, 3,, j, ˆτ 4j bsk is given by Proposition 3.. So the single soliton can move along j directions in x,t plane, which are given by ξ p = 0 associated with j-value of ε p given before. Proof. Comparing the nbkp with the bsk equation, the main change here is the Lax pair L,M. The Lax pair of the nbkp defines the generating functions φ 0 i = φλ i ; x, t are slight different as x n φλ; x, t = λφλ; x, t, φλ; x, t = x 3 φλ; x, t 8 t and then we assume or φλ ; x, t = A e p xp 3 t B e q xq 3 t, p = e iεp, q = e iεp, n = λ, R, 9 φλ ; x, t = A e p xp 3 t B e q xq 3 t, p = e iεp, q = e iεp, n = λ, R. 0

29 SOLVING BI-DIRECTIONAL SOLITON EQUATIONS 9 In order to avoid the divergence of u, we only tae 0 < ε p < π, and then ε p = p π = p π, p = n 4j,, 3,, j. This change results to the emergence of ξ mp = x m cosε p tm 3 cos 3ε p, m n = m j = λ m. The ˆτ nbkp and single soliton solution u = x log ˆτ nbkp can be derived directly from the Proposition. and the generating functions φ 0 in equation 9 associated with λ for the gauge transformation. Further, for a given p, ξ p = 0 determines one moving direction of the single soltion on x,t plane, then the single soliton solution have j directions for propagation because p =,,, j. From Lemmata 6., 6. and the results of the bsk equation, we have Proposition 6.. The single soliton u = x of the nbkp equation, n = log ˆτ nbkp j, j =, 3, 4,, can move along a direction defined by ξ p = 0 on x,t plane for a given p. e iεp, e iεp one distribution of n-th order roots of e iε on circle p = e iεp,q = e iεp in equation 9 The single soliton moves along a line ξp { } = 0 on x, t π plane. Here ε p, 3π, 5π,, j π. 4j 4j 4j 4j 3 For a given n = j, the single soliton of the nbkp have j directions to propagate on x,t plane, which are defined ξ p = 0, p =,, 3, j. Note that the result of j = in above Proposition is given by Lemma 6.. Now we turn to the lower and higher reductions of the CKP hierarchy. Similar to the discussion of the BKP hierarchy in this section, we can obtain parallel results in the CKP hierarchy, so we write out the results without proof in the following to save space. Lemma 6.3. Let ξ = x cosε t 5 cos 5ε, then ˆτ KK ˆτ and the corresponding single soliton is u = 4 cos 5ε cos ε π 6 = 4 > 0. Here ˆτ can be expressed by KK = ˆτ ξ > ξ x. The velocity of soliton is ˆv = log ˆτ KK is given by Proposition 4.. Lemma 6.4. The higher order equation of CKP defined by equation 09. For n-reduction of CKP hierarchy nckp, n = j, j = 3, 4, 5,. Let ξ mp = x m cosε p tm 3 cos 3ε p, m n = = λ m, then the τ function of the nckp generated by T is j m ˆτ nckp = ˆτ and the corresponding single soliton of the nckp equation is u = p ξ > ξ p, x. Here ε p = log ˆτ nckp π = p π, p =,, 3,, j, and ˆτ n 4j is given by Proposition 4.. So the single soliton can move along j directions on x, t plane, which are given by ξ p = 0 associated with j-value of ε p given before. Using the Lemmata 6.3, 6.4 and results for the equation, we get Proposition 6.. The single soliton u = x log ˆτ nckp of the nckp, n = j, j =, 3, 4,, can move along a direction defined by ξ p = 0 on x, t plane for a given p.

30 30 JINGSONG HE, YI CHENG AND RUDOLF A. RÖMER e iεp, e iεp one distribution of n-th order roots of e iε on circle p = e iεp, q = e iεp in equation { 9 the single soliton } moves along a line ξ p = 0 on x, t π plane. Here ε p, 3π, 5π,, j π. 4j 4j 4j 4j 3 For a given n = j, the single soliton of the nckp can move along j directions on x,t plane, which are defined by ξ p = 0, p =,, 3,, j. 4 In particular, if 0 < ε p < π/6, u = x is a two-pea soliton. log ˆτ nckp In above Proposition, the case of j = is given by Lemma 6.3. This Proposition shows there exist several single two-pea solitons for nckp if n. Corollary 6.. There are two single two-pea solitons for -reduction of CKP hierarchy, i. e. CKP equation, u CKP = x log ˆτ, in which ε = π/ and ε = 3π/ respectively. Here ˆτ is given in equation 58 We have plotted it out in figure 7 with = 0.8. Figure 7. Left-going two-pea soliton with dashed line ε = π/ is faster, left-going two-pea soliton with full line ε = 3π/. The left is plotted when t = 0, the right is plotted when t = 0 We have nown that a = / in Lemma 4. and Proposition 4.4 is one crucial point to exist one-pea soliton or two-pea soliton. It is more interesting that a = / will lead to stationary soliton of higher reductions of the BKP and the CKP hierarchy, which is not moving on x,t plane. When ξ ε =π/6 = x cosε t 3 cos 3ε ε =π/6 = x cos ε, ξ is independent with t. So u is independent with t by taing this ξ into Proposition 3. and Proposition 4..

31 SOLVING BI-DIRECTIONAL SOLITON EQUATIONS 3 Corollary 6.. There exists stationary single soliton for the 9-reduction of BKP hierarchy, which is u 9BKP = x log ˆτ bsk ε =3π/8. Here ˆτ bsk is given by Proposition 3.; There exists stationary single soliton for the 9-reduction of CKP hierarchy, which is u 9CKP = x log ˆτ ε =3π/8. Here ˆτ is given by Proposition 4.. We have plotted out stationary soliton for the 9-reduction of CKP in figure 8 when =. Figure 8. Left: stationary soliton for 9-reduction of CKP, Right: Single two-pea soliton for 8-reduction of KP. Corollary 6.3. There is single two-pea soliton u = x log ˆτ the KP hierarchy; there is stationary single one-pea soliton u = ε =π/8 x for 8-reduction of ε =π/6 for log ˆτ the 6-reduction of the KP hierarchy. Here ˆτ is given by Proposition 5.. The two-pea soliton for the 8-reduction of KP hierarchy is plotted in figure 8 when =. For our best nowledge, this is first time to report the even-reduction of the KP hierarchy also has two-pea soliton solution. The possession of two-pea soliton solution is not sole property of CKP hierarchy. 7. Conclusions and Discussions We have presented a systematic way in which to obtain the solution of the n-reduction n = 4, 5 from the general τ function of the KP hierarchy. Our approach is based on the determinant representation of gauge transformations T n 7] and τ n 6]. It may be summarized as follows: τ n constraints of GFs and =n τ n BKP τn CKP 5 reduction τ n =n Eq = bsk, assume the form of φ i and find suitable B i A i n efficient τ function ˆτ Eq =n=,. We have applied this approach to various equations. The one soliton, two soliton and periodic solution are constructed for bsk, and. We show the corresponding relation between the distribution of 5-th Eq

32 3 JINGSONG HE, YI CHENG AND RUDOLF A. RÖMER or 4th roots of e iε on the unit circle and several types of solutions left-going one soliton, right-going one soliton, left/right-going periodic solutions. We also show the reason for the existence of the two-pea soliton. Furthermore, the lower reduction and higher reduction of BKP, CKP, and the even-reductions are explored by this method. Our results show that the soliton of the n-reduction with n = j, j =,, 3,... of BKP and CKP can move alone j directions, which are defined by ξ p = 0. Each direction corresponds to one symmetry distribution of n-th roots of e iε on the unit circle. This supplies a very natural explanation why the 5-reduction BKP or CKP has bi-directional solitons whereas the 3-reduction of BKP or CKP has only single-directional solitons. At last, the two-pea soliton is not a monopolizing phenomena of only the CKP hierarchy. Rather, we find that the higher-order even-reduction of KP also exhibits two-pea solitons and we elucidate the criterion for its existence from the Grammian τ function. At the same time, we show there is not three and more pea soliton from Grammian τ function. The stationary soliton for higher order reduction of KP hierarchy is also obtained. We thin that it is possible to construct an N-soliton solution of the bsk, and equations by this approach. Namely, there exist suitable B i A i i =,,, N such that we can find a physical τ function ˆτ NN Eq for these equations from a complex-valued τ NN Eq, which is symmetric because we have assumed generating functions φ i in equation 4 and equation 43 with symmetric form. Here Eq = bsk,,. Additionally, it is worthy to discuss the phase shift in the collision of one-pea soliton and two-pea soliton. Furthermore, it is possible to construct solutions for bsk, and from constant initial value u = constant 0, which is parallel to present results. Upon completion of this wor, Prof. V. Soolov indly pointed out Ref. 40] where equations 8, 08, and their Lax operators as well as the Lax operator L = u u for KdV equation have been obtained for the first time. Acnowledgement This wor was supported partly by project 973 Nonlinear Science and the National Natural Science Foundation of China We than Prof. Musette for his warm and valuable answering of questions about Ref. 3]. We than Prof. Soolov for information about references and useful discussion on the gauge transformation technique during his visit at Hefei P. R. China. Support from the Warwic Chinese Fellowship Funds is grateful acnowledged. A short version of this paper has been delivered at 005 International Symposium on Nonlinear Dynamics, Celebration of M.S. El-Naschie s 60 Anniversary, December 0-, 005, Shanghai, China