arxiv: v2 [nlin.si] 14 Apr 2012

Transcription

1 Solutions to the modified Korteweg-de Vries equation Da-jun Zhang, Song-lin Zhao, Ying-ying Sun, Jing Zhou Department of Mathematics, Shanghai University, Shanghai 00444, P.R. China arxiv:03.585v [nlin.si] 4 Apr 0 April 7, 0 Abstract This is a continuation of Ref.[]arXiv:nlin.SI/ In the present paper we review solutions to the modified Korteweg-de Vries equation in terms of Wronskians. The Wronskian entry vector needs to satisfy a matrix differential equation set which contains complex operation. This is different from the case of the Korteweg-de Vries equation. We introduce an auxiliary matrix to deal with the complex operation and then we are able to give complete solution expressions for the matrix differential equation set. The obtained solutions to the modified Korteweg-de Vries equation can simply be categorized by two types: solitons and breathers, together with their limit cases. Besides, we give rational solutions to the modified Korteweg-de Vries equation in Wromskian form. This is derived with the help of the Galilean transformed modified Korteweg-de Vries equation. Finally, typical dynamics of the obtained solutions is analyzed and illustrated. We list out the obtained solutions and their corresponding basic Wronskian vectors in the conclusion part. Keywords: The modified Korteweg-de Vries equation, Wronskian, breathers, rational solutions, dynamics MSC 00: 37K40, 35Q35, 35Q35 PACS: 0.30.Ik, 0.30.Jr, Yv Contents Introduction 3 Wronskian solutions of the mkdv equation 4. Preliminary CES of the mkdv equation Simplification of the CES Solutions of the mkdv equation 7 3. Solitons Breathers Rational solutions 4 4. Backgrounds Rational solutions Corresponding author. djzhang@staff.shu.edu.cn

2 5 Dynamics analysis 7 5. Dynamics of solitons Asymptotic behavior of limit solutions Breathers Dynamics of rational solutions Conclusions 7 6. Summation List of solutions A Proof of Theorem 30 B Eigen-polynomial of A = B B 3 C Discussions on the trivial solutions to the CES.8 3 D Proof for Theorem 33

3 Introduction It is well known that the modified Kordeweg-de Vries mkdv equation, v t + 6εv v x + v xxx = 0, ε = ±,. played an important role in constructing infinitely many conservation laws[] and Lax pair for the Korteweg-de Vries KdV equation. The Lax pair led to the breakthrough of the Inverse Scattering TransformIST[3] and then the Soliton Theory. The mkdv equation is also famous for its special soliton behavior, breathers. Actually, the mkdv equation, or its Galilean transformed version, V t + εv 0 V V X + 6εV V X + V XXX = 0, ε = ±,. which is usually referred to as the mixed KdV-mKdV equation, arose in many physics contexts, such as anharmonic lattices[4], Alf ven waves[5], ion acoustic solitons[6, 7, 8], traffic jam[9, 0], Schottky barrier transmission lines[], thin ocean jets[, 3], internal waves[4, 5], heat pulses in solids[6], and so on. With regard to exact solutions, many classical solving methods, such as Hirota s bilinear method[7], the IST[8, 9], commutation methods[0] and Wronskian technique[,, 3] has been used to solve the mkdv equation. For more references on solutions of the mkdv equation, one can refer to Ref.[0] and the references therein. In general, for a soliton equation with bilinear form, its solutions can be expressed through a Wronskian by imposing certain conditions on its Wronskian entry vector[4, 5, 6, ]. For convenience in the following we refer to such conditions as condition equation setces. Usually for an +-dimensional soliton equation the crucial part in its CES is a coefficient matrix and the matrix and its any similar form leads to same solutions for the corresponding soliton equation. Thus it is possible to give a complete classification or structure for the solutions of the soliton equation by considering the canonical form of the coefficient matrix[6, ]. It has been understood that the solutions generated from a Jordan form coefficient matrix are related to the solutions generated from a diagonal form coefficient matrix via some limiting proceure[]. Therefore the latter solutions can be referred to as limit solutions. Actually, from the viewpoint of the IST, N solitons are identified by N distinct eigenvalues of the corresponding spectral problem, or in other words, N distinct simple poles {k j } of transparent coefficient ak. When {k j } are multiple-poles, the related multiple-pole solution can be obtained through a limiting procedure like k k from simple-pole solution. This limiting procedure is easily realized for solutions in Wronskian form[]. Such a procedure is also helpful to understand the dynamics of limit solutions[7, 8, 9]. In Ref.[] we mentioned four topics related to solutions in Wronskian form: to find the CES, to solve the CES, to describe relations between different kinds of solutions, and to discuss dynamics of the solutions. In the present paper, following this line, we will try to review the Wronskian solutions of the mkdvε = equation together with its Galilean transformed version, The main results of the paper are the following. v t + 6v v x + v xxx = 0,.3 V t + v 0 V V X + 6V V X + V XXX =

4 The CES for the mkdv equation can be given by ϕ x =B ϕ, ϕ t = 4ϕ xxx,.5a.5b where ϕ is the N-th order Wronskian entry vector, bar stands for complex conjugate and B is a nontrivial N N constant complex matrix. We note that there is a complex conjugate involved in the CES and this makes difficulties when solving the CES. There are only two kinds of solutions led by the CES.5: Solitons together with their limit case and breathers together with their limit case. No rational solutions arise from.5 due to B 0. Rational solutions to the mkdv equation.3 can be derived by means of the KdV-mKdV equation.4. Dynamics of obtained solutions is analyzed and illustrated. The paper is organized as follows. In Sec., we give a general CES of the mkdv equation and simplify the CES by introducing an auxiliary equation. Then in Sec.3 we solve the CES and classify the solutions as solitons and breathers. In Sec.4 we derive rational solutions to the mkdv equation. This is done with the help of the KdV-mKdV equation.4. Sec.5 consists of dynamic analysis and illustrations. Finally, in the conclusion section we list out the obtained solutions and their corresponding basic Wronskian vectors. Wronskian solutions of the mkdv equation. Preliminary An N N Wronskian is defined as W φ, φ,, φ N = φ, φ,, φ N =,..... φ 0 N φ N φ N N φ 0 φ φ N φ 0 φ φ N where φ l j = l φ j / x l and φ = φ, φ,, φ N T is called the entry vector of the Wronskian. Usually we use the compact form[4] W φ = φ, φ,, φ N = 0,,, N = N,. where N j indicates the set of consecutive columns 0,,, N j. In the paper we also employ the notation Ñ j to indicate the set of consecutive columns,,, N j. A Wronskian can provide simple forms for its derivatives and this advantage admits direct verification of solutions that are expressed in terms of Wronskians. The following matrix properties are usually necessary in Wronskian verification. 4

5 Proposition. [30, ] Suppose that Ξ is an N N matrix with column vector set {Ξ j }; Ω is an N N operator matrix with column vector set {Ω j } and each entry Ω j,s being an operator. Then we have N N Ω j Ξ = Ω T j Ξ T,.3 j= j= where for any N-order column vectors A j and B j we define A j B j = A,j B,j, A,j B,j,, A N,j B N,j T.4 and A j Ξ = Ξ,, Ξ j, A j Ξ j, Ξ j+,, Ξ N..5 Proposition. [4] Suppose that D is an N N matrix and a, b, c, d are Nth-order column vectors, then D, a, b D, c, d D, a, c D, b, d + D, a, d D, b, c = 0..6 Besides, to write solutions of CES in simple forms one may make use of lower triangular Toeplitz matrices. An N th-order lower triangular Toeplitz matrix means a matrix in the following form A = a a a a a a a N a N a N 3 a a 0 N N, a j C..7 All such matrices form a commutative semigroup G N C with identity with respect to matrix multiplication and inverse, and the set G N C = { A A GN C, A 0 } makes an Abelian group. Besides.7, we will also need the following block lower triangular Toeplitz matrix, A B = A A A A A A A N A N A N 3 A A 0 N N,.8 aj 0 where A j = and {a 0 a js } are arbitrary complex numbers. All such block matrices j also compose a commutative semigroup with identity which we denote by G B N C, and the set G B N C = { A B A B G B N C, AB = 0 } makes an Abelian group, too. If all the elements are real, then we correspondingly denote the above mentioned matrix sets by G N R, G N R, G B N R and GB N R. For more properties of such matrices please refer to Ref.[]. 5

6 . CES of the mkdv equation By the transformation v = i ln f f the mkdv equation.3 can be bilinearized as[7] x = i f x f ff x ff,.9 D t + D 3 x f f = 0,.0a D x f f = 0,.0b where i is the imaginary unit, f is the complex conjugate of f, and D is the well-known Hirota s bilinear operator defined by[3, 3] D m t D n xat, x bt, x = m s m n y n at + s, x + ybt s, x y s=0,y=0, m, n = 0,,. The bilinear mkdv equation.0 admits a solution f in Wronskian form. Theorem. A Wronskian solution to the bilinear mkdv equation.0 is given as provided that its entry vector φ satisfies f = W φ = N,. φ x =Bt φ, φ t = 4φ xxx + Ctφ,.a.b where Bt = B ij t N N and Ct = C ij t N N are two N N matrices of t but independent of x, and satisfy Bt t = 0, trct Rt,.3a B t t + Bt Ct = CtBt..3b The proof is given in A..3 Simplification of the CES. To solve the CES. with arbitrary Bt and Ct which satisfy.3 we first introduce a non-singular N N complex matrix Ht C N N t such that [33], also see [] By Ht we then introduce a new Wronskian entry vector which transfers the CES. to the following, H t t = HtCt..4 ψ = Htφ,.5 ψ x = B ψ,.6a ψ t = 4ψ xxx,.6b 6

7 where B = HtBt H t has to be a constant matrix independent of both x and t due to the compatibility condition.3b noting that now Ct = 0. We note that the Wronskians composed by φ and ψ, which we respectively denote by fφ and fψ, yield same solutions to the mkdv equation through the transformation.9 due to fψ = Ht fφ. That means in the following one only needs to focus on the CES.6. However, since there exists a complex operation in.6, solutions can not be classified in terms of the canonical form of B, as done in [6, ]. To overcome the difficulty we introduce an auxiliary equation ψ xx = Ãψ,.7 where à = B B. Our plan now is to first solve the equation set composed by.7 and.6b, which is nothing but the CES of the KdV equation and have been well studied in Refs.[6] and []. Then in the second step we impose condition.6a on the obtained solution ψ and finally get solutions for.6. In addition, noting that à and its any similar form A = P ÃP always generates same solutions to the mkdv equation, in the following let us focus on the CES ϕ xx =Aϕ, ϕ x =B ϕ, ϕ t = 4ϕ xxx,.8a.8b.8c where ϕ = P ψ, B = P B P and A = B B..9 In the next section we will see that solutions to the mkdv equation can be classified in terms of the canonical form of A rather than the canonical form of B. 3 Solutions of the mkdv equation In this section, we list several possible choices of A and derive the related Wronskian entry vectors, one of which is for breathers. We will also discuss the limiting relationship of some solutions. First, for the matrix A defined by.9, we have the following result on its eigenvalues. Proposition 3. The eigenvalues of A defined by.9 are either real or, if there are some complex ones, appear as conjugate pairs. We leave the proof in B. According to the eigenvalues of A, we can categorize solutions to the mkdv equation as solitons and breathers, which correspond to real eigenvalues and complex eigenvalues of conjugate pairs, respectively. 3. Solitons Case I. Solitons: reads When A has N distinct real positive eigenvalues {λ j } its canonical form A = Diagλ, λ,, λ N, 3. 7

8 where, for convenient to discuss, we let λ j = ε j k j 0, in which ε j = ± and k j can either real or complex numbers with distinct absolute values. We consider two subcases.. k j R, i.e. A = Diagk, k,, k N. 3. Following the relation.9, we can take B = Diag±k, ±k,, ± k N. 3.3 We neglect the sign ± because this can be compensated by the arbitrariness of k j. So next we take B = Diagk, k,, k N, k j R. 3.4 For the above matrix B, the solution to the CES.8 can be ϕ = ϕ, ϕ,, ϕ N T, 3.5a where When a + j = j, a j ϕ j = a + j eξ j + ia j e ξ j, ξ j = k j x 4k 3 j t + ξ 0 j, a + j, a j, k j, ξ 0 j R. 3.5b =, 3.5b reads ϕ j = j e ξ j + ie ξ j. 3.6 In this case, the corresponding Wronskian solution.9 with fϕ can be written as[34] N f = j= e ξ j j<l N k j k l µ=0, { N exp µ j η j + π i + j= j<l N where the sum over µ = 0, refers to each of µ j = 0, for j =,,, N, and η j = ξ j 4 N l=,l j A jl, e A kl k jl j =. k l + k j µ j µ l A jl }, This coincides with the N-soliton solution in Hirota s exponential polynomial form[37].. k j = k j + ik j C, i.e. In this case we have A = Diagλ, λ,, λ N, λ j = k j + k j. 3.7 B = Diagk, k,, k N, k j C, 3.8 and the solution to the CES.8 can be given by 3.5a with ϕ j = γ j a + j eξ j + ia j e ξ j, ξ j = λ j x 4λ 3 jt + ξ 0 j, a + j, a j, ξ0 j R, 3.9a where γ j = + iλ j k j k j, λ j = ε j k j. 3.9b 8

9 Now, if we compare 3.5b and 3.9a, we can find that both of them provide same solution to the mkdv equation through the transformation.9. In particular, when k j = 0, i.e., and we take λ j = k j, the entry function 3.9a reduces to B = Diagik, ik,, ik N, 3.0 ϕ j = + ia + j eξ j + ia j e ξ j, ξ j = k j x 4kj 3 t + ξ 0 j, a + j, a j, ξ0 j R. 3. This was first given by Nimmo and Freeman[] as Wronskian entries for soliton solutions. Let us remark this case as follows. Remarks: When B in the CES.8b and.8c is a diagonal matrix 3.8 of which the diagonal elements have different absolute values, no matter these diagonal elements are real or complex, the related Wronskian fϕ generates N-soliton solutions to the mkdv equation. Besides, we specify the non-degenerate condition k i k j, i j, 3. i.e., k i and k j i j can not appear on the same circle with the original point as the center of the circle. Otherwise, the solution degenerates. Case II. Limit solutions of solitons: Corresponding to Case I, we discuss two subcases.. Let k k k A = k k 0 0 0, k R k k In this subcase, general solution to the equation set.8a and.8c is[] ˆφ = AQ BQ 0, A, B G N C, 3.4 N where Q ± 0 = Q± 0,0, Q± 0,,, Q± 0,N T, Q ± 0,s = s! s k e ±ξ, 3.5 ξ is defined in 3.5b, GN C is the commutative set of all the Nth-order lower triangular Toeplitz matrices, see Sec... The above non-degenerate relation can also be described as follows. Define the equivalent relation on the complex plane C: k i k j iff k i = k j. The quotient space C/ denotes the positive half real axis. Then to get a non-degenerate soliton solution one needs k i k j, i j. 9

10 For the matrix B satisfying.9 we take B = Then, substituting into.8b yields Meanwhile, it can be verified that k k k N. 3.6 AQ + 0,x + BQ 0,x = BĀQ+ 0 + BQ Q + 0,x = BQ+ 0, Q 0,x = BQ Noting that B G N R, substituting 3.8 into 3.7, and making use of the commutative property of G N C, we have Then, compared with 3.7 we immediately get BAQ + 0 BQ 0 = BĀQ+ 0 + BQ A = Ā, B = B, 3.0 which means A is real and B pure imaginary. In the end, the solution to the CES.8 can be described as ϕ = A + Q ia Q 0, A± G N R. 3. We note that solutions generated from 3. can also be derived from the solution given in Case I. by a limiting procedurecf.[]. Let us explain this procedure by starting from the following Wronskian fϕ = W ϕ, ϕ,, ϕ N N j= k j k j 3. with ϕ = ϕ k, x, t = a + eξ +ia e ξ as defined in 3.5b and ϕ j = ϕ k j, x, t. 3. gives a Wronskian solution to the bilinear mkdv equation.0. Taking the limit k j k successively for j =, 3,, N and using L Hospital rule, the Wronskian 3. goes to a Wronskian W ϕ with entry vector 3., where the arbitrary coefficient matrices A, B can come from a ± by considering a ± to be some polynomials of k cf.[].. Corresponding to Case I., let us consider A = k + k k k + k k k + k k k + k N, 3.3 0

11 where k, k R and k 0. It then follows from.9 that one can take B = ik ik ik N, k = k + ik. 3.4 For the matrix A defined by 3.3, the general solution to the equation set.8a and.8c can be written as ϕ = AP BP 0, A, B G N C, 3.5 where P 0 ± = P 0,0 ±, P± 0,,, P± 0,N T, P 0,s ± = s! s k e ±ξ, 3.6a ξ = λ x 4λ 3 t + ξ 0, λ = k + k. 3.6b Next, we substitute 3.5 together with 3.6 into.8b so that we identify A, B for.8b. This substitution yields AP + 0,x + BP 0,x = BĀP+ 0 + BP 0 = ĀBP+ 0 + BBP 0, 3.7 where we have made use of the commutative property of G N C. Then, noting that where W = P + 0,x = W P+ 0, P 0,x = W P 0, 3.8 w w w w w w w N w N w N 3 w w 0 and substituting 3.8 into 3.7, we get the relation To have a clearer result, we write N, w j = j! j k λ, 3.9 AW = ĀB, BW = BB A = A + ia, B = B + ib, B = B + ib, 3.3 where A, A, B, B are in G N R and k 0 k B =......, B = k 0 0 k k 0 0 k Here we define P ± 0,s by taking derivative with respect to k rather than k because in this case we always have k 0.

12 Then 3.30 yields A W B A B = 0, A B A B + W = 0, 3.33 B W + B + B B = 0, B B B B W = 0, 3.34 which provides the equation sets for A, A, B, B. We note that since all the elements in 3.33 and 3.34 are in G N R, we can treat 3.33 and 3.34 as ordinary linear equation sets. Let us look at It first indicates A = B W B A, 3.35 and further, by eliminating A, B + B W A = We note that In fact, by calculation we find B + B W = { 0, i < j, W i,j = i j! i j k k + k, i j, from which it is easy to verify Thus, it turns out that A can be an arbitrary element in G N R and then A = A M, M = I N + ib W B, 3.38 where I N is the Nth-order unit matrix. Similarly, we can find B = ib M, B G N R. Thus we can conclude that the solution to the CES. can be given by ϕ = MA + P i A P 0, A± G N R, 3.39 where M is defined in Obviously, the matrix M contributes nothing through the transformation.9 to the mkdv equation. Therefore in practice we may remove M and use the effective part ϕ = A + P i A P 0, A± G N R, 3.40 which is as same as Breathers Case III. Breathers: When A has distinct complex conjugate-pair eigenvalues, we may have breather solutions. Let us consider a N-th order matrix, A = Diagk, k,, k N, k N N, k j 0, 3.4 where k j C, j =,,, N. The matrix B which generates breathers is 0 kj B = DiagΘ, Θ,, Θ N N, Θ j =. 3.4 k j 0

13 For this case, the Wronskian entry vector, i.e., the solution to the CES.8, can be taken as ϕ = ϕ, ϕ, ϕ, ϕ,, ϕ N, ϕ N T, 3.43a where ϕ j = a j e ξ j + b j e ξ j, ϕ j = ā j e ξ j b j e ξ j, ξ j = k j x 4k 3 j t + ξ 0 j, a j, b j, ξ 0 j C. 3.43b With the above ϕ as basic entry vector, the Wronskian fϕ will provide breather solutions for the mkdv equation. For non-trivial solutions we need k j 0, see Sec.5.3. Case IV. Limit solutions of breathers: In this case let us consider the following block matrix, K K K A = I K K , 3.44a I K K where k K = 0 0 k, K = k 0 0 k The matrix B satisfying.9 can be taken as B Ĩ B = B Ĩ B where Ĩ = 0 0, B = 0 k k 0 The Wronskian vector of this case is also in the form of N 0, I = 0 N. 3.44b, 3.45a. 3.45b ϕ = ϕ,, ϕ,, ϕ,, ϕ,,, ϕ N,, ϕ N, T, 3.46a and for convenience we set ϕ + = ϕ,, ϕ,,, ϕ N, T, ϕ = ϕ,, ϕ,,, ϕ N, T. 3.46b Substituting 3.46 into.8 with the above A and B we get ϕ + xx = A N ϕ +, ϕ xx = ĀNϕ, 3.47a ϕ + x = B N ϕ, ϕ x = B N ϕ +, 3.47b ϕ ± t = 4ϕ ± xxx, 3.47c 3

14 where A N = B N = k k k k k k k From 3.47a, 3.47c and 3.48a, one first has k k k N N, 3.48a. 3.48b ϕ + = AQ BQ 0, ϕ = C Q D Q 0, A, B, C, D G N C, 3.49 where Q ± 0 are defined by 3.5 together with ξ defined in 3.43b. Then substituting 3.49 into 3.47b and making using of 3.8, one finds Thus, 3.49 reads C = Ā, D = B ϕ + = AQ BQ 0, ϕ = Ā Q + 0 B Q 0, A, B G N C. 3.5 The Wronskian fϕ with ϕ defined by 3.46 with 3.5 will provide a limit solution of breathers for the mkdv equation. In fact, quite similar to the procedure for the limit solutions of solitons we described in Sec.3., here the Wronskian fϕ with 3.46 is related to the N-breather solution by taking the limit k j k successively for j =, 3,, N. 4 Rational solutions 4. Backgrounds Following solution structures of the KdV equation[], rational solutions should be led from zero eigenvalues of the coefficient matrix A. For the mkdv equation this requires B = 0. However, this is not allowed in the Wronskian verification see Appendix C.4 because a trivial B will lead to a zero Wronskian in the light of A.. Thus, it is clear that to get non-trivial rational solutions to the mkdv equation, we need a non-trivial matrix B. Let us go back to the KdV-mKdV equation.4, i.e., V t + v 0 V V X + 6V V X + V XXX = 0, 4. which is related to the mkdv equation.3 through a Galilean transformation vx, t = v 0 + V X, t, x = X + 6v 0 t, 4. where v 0 is a real parameter. We note that the equation 4. admits non-trivial and nonsingular rational solutions when v 0 0. Then using the transformation 4. rational solutions to the mkdv equation can be obtained. This fact has been realized via Bäcklund transformationbt[35] and Hirota method with a limiting procedure[36, 37], but the presentation for high order rational solutions is complicated. In the following we will derive rational solutions in terms of Wronskian, which provides not only explicit but also impact forms for high order rational solutions. 4

15 4. Rational solutions Still employing the same transformation as.9, i.e., V = i ln f f, 4.3 X the KdV-mKdV equation 4. can be written into the bilinear form[38, 39] D t + D 3 X f f = 0, 4.4a D X iv 0 D X f f = b For the solutions to 4.4 in Wronskian form, we have Theorem. The bilinear equation 4.4 admits Wronskian solution f = W φ = N, 4.5 where the entry vector φ satisfies iφ X = v 0 φ + Bt φ, φ t = 4φ XXX + Ctφ, 4.6a 4.6b in which Bt and Ct are two N N matrices of t but independent of x, and satisfy Bt 0, trct Rt, B t t + Bt Ct CtBt = a 4.7b 4.7c The proof is similar to the one for Theorem, but 4.6a results in complicated expression for f and its derivatives. We leave the proof in D. As in Sec..3, with the help of the auxiliary matrix A and auxiliary equation one can simplify the CES 4.6 to ϕ XX = Aϕ, 4.8 iϕ X = v 0 ϕ + B ϕ, ϕ t = 4ϕ XXX, 4.9a 4.9b where both A and B are N N complex constant matrices and related by A = B B v 0I N, 4.0 in which I N is the Nth-order unit matrix. It might be possible that here we discuss all possible solutions according to the eigenvalues of A, as we have done in the previous section. However, since we have had a clear description for solitons and breathers of the mkdv equation in the previous section, and the parameter v 0 will bring more complexity, in the following we can 5

16 neglect the discussion of A and let us only focus on rational solutions. In fact, rational solutions correspond to the zero eigenvalues of A. We derive rational solutions as the limit solutions of solitons. The N-soliton solution corresponds to B = diag v v 0 + k, 0 + k,, v0 + k N, 4. where we take k j to be N distinct real positive numbers. In this case, a solution to the CES 4.9 is ϕ = ϕ, ϕ,, ϕ N T, 4. with ϕ j = v 0 + ik j e η j + v 0 ik j e η j, η j = k j X 4k 3 j t, 4.3 This provides an N-soliton solution to the KdV-mKdV equation 4. through the transformation 4.3 with fφ. Now in the CES 4.9a we take B to be a lower triangular Toeplitz matrix, with B = j α 0 α j = j! j k In this case the CES 4.9 admits a solution with α α , 4.4 α N α α 0 v 0 + k, j = 0,,, N. 4.5 k =0 ψ = ψ, ψ,, ψ N T, 4.6a ψ j+ = j j! j k ϕ j = 0,,, N, 4.6b k =0, where ϕ is defined by 4.3. Then, with such a ψ as basic column vector, the Wronskian f = W ψ = N 4.7 provides non-singular rational solutions to the KdV-mKdV equation 4.. A simplified form of these solutions is V X, t = F,XF F F,X F + F, F = Re[f], F = Im[f]. 4.8 We list the first three non-trivial f for N =, 3, 4, respectively in the following, f =4v 0 X + i, f = 6v 0 v0 [ X 3 + t 3X 3 4v0 f = v 4 0 4X v 0 4.9a + 3i X + ] v 0 4v0, 4.9b 6 3 tx + X v 070t 60tX 3 X 6 + 4iv 0 4t v0 + X v tX X5. 4.9c For the rational solutions to the mkdv equation.3, we have 6

17 Proposition 4. The non-trivial rational solutions to the mkdv equation is given by vx, t = v 0 F,XF F F,X F + F, X = x 6v0t, 4.0 where f is the Wronskian 4.7 composed with 4.6, F = Re[f], F = Im[f]. The first non-trivial rational solution to the mkdv equation reads and the second one v = v 0 v = v 0 4v 0 4v 0 X +, X = x 6v 0t, 4. v 0 X X 3 4Xt v0 6v0 4 4v0 X3 + t 3X 4v + 9X +, X = x 0 4v 6v 0t We note that there is a bilinear BT[] related to 4.4 and the BT admits multi-soliton solutions[] and rational solutions[40] in Wronskian form. Besides, rational solutions of the mkdv equation for both ε = ± may also be derived through the Galilean transformed version by using the long-wave-limit approach described in Ref.[4] also see [3, 4]. In next section we will discuss dynamics of solutions including these rational solutions. 5 Dynamics analysis In this section we investigate dynamics of two-soliton solutions, limit solutions, breathers and rational solutions. To describe the relationship between two-soliton solution and the simplest limit solution, let us start from the asymptotic behaviors of two-soliton interactions. 5. Dynamics of solitons Soliton solutions to the mkdv equation.3 can be described by vx, t = F,xF F F,x F + F, F = Re[f], F = Im[f], 5. where f = fϕ = N is the Wronskian composed by the basic column vector ϕ which is defined by 3.5, or equivalently, by either 3.9a or 3. in Case I. We note that from the transformation.9 the solution to the mkdv equation.3 can also be written as v = arctan F F = arctan F x F, F = Re[f], F = Im[f], 5. x while 5. gives a more explicit form. In the following we investigate - and -soliton solutions, which are, respectively, corresponding to f = f =fϕ = ϕ, 5.3a f = f =fϕ, ϕ T = ϕ ϕ ϕ,x, ϕ,x 5.3b 7

18 with ϕ j defined by 3.5b, i.e., ϕ j = a + j eξ j + ia j e ξ j, ξ j = k j x 4k 3 j t + ξ 0 j, a ± j, k j, ξ 0 j R. 5.4 Then, one-soliton solution to the mkdv equation.3 reads a ] v = sgn[ a + k sech k x 8kt 3 ln a a +, 5.5 as depicted in Fig., where for convenience, we call a soliton and b anti-soliton due to the signs of their amplitudes x a x b Figure : The shape of one-soliton solution given by 5.5 at t = 0. a. Soliton for a + = a =, k = and ξ 0 = 0. b. Anti-soliton for a + = a =, k = and ξ 0 = 0. Obviously, this soliton is identified by the amplitude and top trace trajectory Amp = sgn [ a ] a + k xt = 4kt + ln a, k a or velocity 4k. Obviously, solitons of the mkdv equation are single-direction waves. Next, let us look at two-soliton solution. The two-soliton solution of the mkdv.3 can be expressed by 5. where from 5.3b F = k k a a e4k3 +k3 t k +k x + a + a+ e 4k3 +k3 t+k +k x, F = k + k a a+ e4k3 k3 t k k x a + a e4 k3 +k3 t+k k x. 5.7a 5.7b We assume that a ± j 0 and set sgn [ a a + ] [ a ] = sgn a, so that + F, F are not zero at same time. Without loss of generality, we also take 0 < k < k. For the analysis of asymptotic behaviors, it is convenient to use the following expression v = arctan F F. 5.8 x 8

19 a b Figure : Two-soliton interactions. a. Soliton-soliton for a + = a = a =, b =, k = 0.8, k = 0.5 and ξ 0 = ξ 0 = 0. b. Soliton-anti-soliton for a + = a = a =, b =, k = 0.8, k = 0.5 and ξ 0 = ξ 0 = 0. There are two types of -soliton interactions, soliton-soliton or anti-soliton-anti-soliton interaction and soliton-anti-soliton interaction, as shown in Fig.. To investigate asymptotic behaviors of the two solitons involved in interaction, we first name them k -soliton and k -soliton, respectively. Then we rewrite the two-soliton solution 5.8 in terms of the following coordinates, X = x 4k t, t, 5.9 which then gives v = arctan a+ a ek X a + a e8k k k t+k X k + k a a a+ a+ e8k k k t+k +k X k k X Noting that for any 0 < k < k it is always valid that k k > 0 and k +k k k > 0, we can keep X to be constant and let t go to infinity. Then we can find there is only k -soliton left along the line X = const. and also find how the k -soliton is asymptotically identified by its top trace and amplitude, for both t ±. As for details, when k > 0, t or k < 0, t +, i.e. sgn[k ] t, the solution 5.0 becomes v = arctan a+ k + k e k X a k k X a ] = sgn[ a + k sech k X ln a a + + ln k + k ; 5. k k and when k > 0, t + or k < 0, t, i.e. sgn[k ] t +, 5.0 becomes v = arctan a k + k a + ekx k k = sgn[ a a + ] k sech X k X ln a a + ln k + k. 5. k k 9

20 We can also rewrite the two-soliton solution 5.8 in terms of the coordinates X = x 4k t, t, 5.3 and do a similar asymptotic analysis for the k -soliton. Finally, we reach to Theorem 3. Suppose that sgn [ a ] [ a ] a = sgn + a, a ± + j 0 and 0 < k < k in 5.0. Then, when sgn[k ] t ±, the k -soliton asymptotically follows top trace : xt = 4kt + ln a k a + ± ln k + k, 5.4a k k k a ] amplitude : sgn[ a + k, 5.4b and when sgn[k ] t ±, the k -soliton asymptotically follows top trace : xt = 4kt + ln k amplitude : sgn[ a a + a a + The phase shift for the k j -soliton after interactions is k j ln k +k k k. ± k ln k + k k k, 5.5a ] k. 5.5b Now it is completely clear how the two-soliton interactions are related to the parameters {k j, a j, b j }. This will be helpful to understand the asymptotic behavior of limit solutions. 5. Asymptotic behavior of limit solutions The simplest limit solution in Case II is v = arctan a a + a 48k3 t + 4k x e 8k 3 t k x + a + e 8k 3 t+k x x. 5.6 This is derived from 5. with f = ϕ ϕ,x k ϕ k ϕ x, 5.7a where ϕ is defined in 3.5b, i.e., ϕ = a + eξ + ia e ξ, ξ = k x 4k 3 t + ξ 0, a±, k, ξ 0 R. 5.7b The solution 5.6 is depicted in Fig.3. We characterize dynamics of the solution by the following two points, which are typically different from the interaction of two normal solitons that we described in the previous subsection. These two points are Soliton-anti-soliton interaction with asymptotically asymmetric wave shape, Top trace of each soliton is asymptotically governed by logarithm and linear functions. 0

21 a b Figure 3: Limit solution given by 5.6 for a + = a =, k = 0.5 and ξ 0 = 0. a. Shape and motion. b. Asymmetric wave shape at t = 60. The first point can be explained as follows. Recall the two-soliton interaction with a ± = a ± 0 and k k > 0. According to Theorem 3, this is a soliton-anti-soliton interaction and the absolute of amplitude of each soliton is k j. Obviously, the asymmetric shape of limit solution 5.6 coincides well with the limit k k. To understand the second point, again we put the limit solution in the coordinate system Y = x 4k t, t, 5.8 and this gives which is described in Fig.4. v = arctan 4a+ a eky k Y 8k t a + a +, 5.9 e 4k Y Y a b Figure 4: Limit solution given by 5.9 for a + = a = and k = 0.5. a. Shape and motion. b. Asymmetric wave shape at t = 60. For convenience here we suppose that k > 0. Then by analyzing the leading terms as t ± in the numerator and denominator in 5.9, we can conclude the asymptotic behaviors of the limit solution 5.9 as follows.

22 Theorem 4. Suppose that a ± 0 and k > 0 in 5.9. When t there are two waves moving with amplitudes sgn [ a ] a + k and top traces governed by the logarithm functions Y ± = ln k a a + ± ln 3kt 3, 5.0a where the subscript ± of Y stands for Y ±. When t +, there are also two waves moving with amplitudes ± sgn [ a ] a + k and top traces Y ± = ln k a a + ± ln3kt b Replacing Y by x by using 5.8 the top traces of the waves in Fig.3 are then asymptotically governed by linear and logarithm functions of t. Let us end up this subsection with the following remark. With regard to the limit solutions the top traces or waves trajectories are governed by linear and logarithm functions should be a typical characteristic, which differs from normal soliton interactions with straight line trajectories. See [30, 7, 8, 9] for more examples. 5.3 Breathers Wronskian entries in Case III provides breather solutions to the mkdv equation.3. The simplest one corresponds to where f = ϕ ϕ,x ϕ ϕ,x = F + if, 5.a ϕ = a e ξ + b e ξ, ϕ = ā e ξ b e ξ, ξ = k x 4k 3 t, k, a, b C, 5.b 5.c and we have written F =4k a b + a b cos4k k t 8k 3 t k x + 4k a b a b sin4k k t 8k 3 t k x, 5.d F = k e k 4k t+k t+x[ b + b e 6k3 t + a + a e 4k k t+x], 5.e k = k + ik, a = a + ia, b = b + ib. Then the breather solution is expressed by 5. with the above F, F. We further assume that sin θ = a b + a b, α = a b + a b α + a b a b,

23 and then rewrite 5. as where v = arctan P Q, x 5.a P = k α sin [ k x 4t3k k θ ], Q = k [ a + a e k x+4t3k k + b + b e k x+4t3k k ]. 5.b 5.c Such a breather is described in Fig.5. a b c Figure 5: Breathers given by 5.. a. A moving breather with a = + i, b = + 0.5i, k = i. b. A stationary breather with a = b =, k = 3 + i. c. A D-plot of b at t = 0. Fig.5a shows an oscillating wave moving along a straight line. The oscillation comes from the sine function and the frequency depends on both x and t. To understand more on the wave we use the coordinates Z = x + 4t3k k, t 5.3 to rewrite the solution 5. as v = arctan k α sin k Z 8tk + k θ k a + a ekz + b + b e k Z, 5.4 Z and then fix Z = 0, i.e., looking at the wave along the straight line Z = 0. Then it is clear that The breather travels along the straight line Z = 0, in other words, the wave speed is 43k k, which means it admits bi-direction traveling. A stationary breather appears when 3k = k, as depicted in Fig.5b. Under the coordinate system 5.3 the frequency only depends on t and the period reads π T = 8k k k. From 5. we can see that k j 0 is necessary for getting nontrivial breather solutions. 3

24 We note that in two dimensions fixing t the breather is in fact a spindle-shape wave. Let us go back to the solution 5.. If we fix time t, then the breather oscillates with frequency k π and its amplitude decays by the rate e xk as x ±. That means if k is small and k is relatively large so that the amplitude decay becomes the dominating factor, we get a normal breather as shown in Fig.6a; while if k is small enough and k >> k so that the oscillation dominates, we will see a spindle-like wave shown in Fig.6b. In the latter case, the wave will travel with high speed and high oscillating frequency. This can lead to overlaps of normal oscillating waves like Fig.6a during their traveling, which makes a spindle shape. a b Figure 6: Shape of breathers given by 5.. a. Normal shape at t = 0 with a = + i, b = + 0.5i, k = i. b. Spindle-shape breather at t = 0 with a = + i, b = + i, k = i. Finally, in this subsection we list and depict two-breather solution and the simplest limit breather solution, without further asymptotic analysis. Both solutions can be given by 5. with f being a 4 by 4 Wronskian where for the -breather solution ϕ = ϕ, ϕ, ϕ, ϕ T, f = ϕ, ϕ x, ϕ xx, ϕ xxx, a ϕ j = a j e ξ j + b j e ξ j, ϕ j = ā j e ξ j b j e ξ j, ξ j = k j x 4k 3 j t + ξ 0 j, 5.7b in which a j, b j, k j, ξ 0 j C, and for the limit breather ϕ = ϕ, ϕ, k ϕ, k ϕ T, 5.8 in which ϕ, ϕ are defined by 5.7b. Fig.7 shows the two-breather interaction where from the density plot b one can clearly see that the two breathers are traveling along straight lines and a phase shift appears after interaction. Fig.8 shows the shape and motion of a limit breather solution, where from the density plot b one can clearly see that the breather trajectories are not any longer straight lines, Here we conjecture that they are governed by logarithm functions. 4

25 a b Figure 7: Shape and motion of two breather solution given by 5. with 5.6 and 5.7. a. 3D-plot for a = b = a = b =, k = + 0.5i, k = i and ξ 0 = ξ 0 = 0. b. Density plot of a for x [, ], t [ 4, 4]. a b Figure 8: Shape and motion of the limit breather solution given by 5. with 5.6 and 5.8. a. 3D-plot for a = b =, k = i and ξ 0 = 0. b. Density plot of a for x [ 0, 0], t [ 8, 8]. 5.4 Dynamics of rational solutions The first non-trivial rational solution to the mkdv equation is 4., i.e., v = v 0 4v 0 4v 0 X +, X = x 6v 0t. 5.9 This is a non-singular traveling wave moving with the constant speed 6v0, constant amplitude 3v 0 and asymptotic line v = v 0. It is depicted in Fig.9. The next rational solution is given by 4., i.e., v = v 0 v 0 X X 3 4Xt v0 6v0 4 4v0 X3 + t 3X 4v + 9X +, X = x 0 4v 6v 0t

26 a b c Figure 9: Shape and motion of the rational solution given by 5.9 for v 0 = 0.8 in a, v 0 = in b and v 0 = in c. It can be viewed as a double-traveling wave solution with v = v 0 v 0 X XY 3 v0 6v0 4 4v0 X3 3Y 4v + 9X +, v 0 X = x 6v 0t, Y = x v 0t. 5.3 However, it does not show interactions of two single rational solutions. Only one wave is left for large x, t see Fig.0. We re-depict Fig.0a in Fig. by a density plot so that we can see the wave top trace clearer. a b Figure 0: Shape and motion of the rational solution given by 5.30 for v 0 = 0.8 in a and v 0 = 0.8 in b. To realize the asymptotic behavior analytically, we rewrite the solution 5.30 in the following coordinates system X, T = X 3 + t + 3 6v0 4X

27 a b Figure : a: Density plot of Fig.0 a with X [ 0, 0], t [ 30, 30]. b: a overlapped by the trajectory curve given by and this gives v = v 0 4v0 T 3 v 0 3X X T X + 3 v0 6v0 4 6v0 4X 3X 4v0 + 9X v 0 Then, by calculation it can be found that for given large X the wave 5.34 has a single stationary point at T = 0 where v gets a local extreme value v = v 0 v 0 X + 4v 0 X, 5.35 which goes to 3v 0 as X ±. Thus we can conclude that for large X, t the wave asymptotically travels along the curve T = X 3 + t + 3 6v 4 0 X = with amplitude 3v 0. Fig.b displays a density plot overlapped by the above wave trajectory curve. More details on the rational solutions to the mkdv equation can be found in [40] where the rational solutions are derived via bilinear Bäcklund transformation. 6 Conclusions 6. Summation In the paper we reviewed the Wronskian solutions to the mkdv equation.3 in terms of Wronskians. When a solution is expressed through the Wronskian f = fϕ = N, 6. one needs to solve the finalized CES.8 together with.9, i.e., ϕ xx =Aϕ, ϕ x =B ϕ, ϕ t = 4ϕ xxx, 6.a 6.b 6.c 7

28 and A = B B. 6.3 A is the auxiliary matrix that we introduced to deal with the complex operation in 6.b and it works in practice. As a result, with the help of A we solved the above CES and then categorized the solutions to the mkdv equation in terms of the canonical form of A rather than canonical form of B. Solutions are categorized by solitons together with their limit case and breathers together with their limit case. There are no rational solutions arising from 6. because no rational solutions correspond to zero eigenvalues of A while we need B 0 to finish Wronskian verification. To derive rational solutions for the mkdv equation.3, we employed the Galilean transformed equation, i.e., the KdV-mKdV equation.4 which admits rational solutions in Wronskian form. Then the rational solutions to the mkdv equation can be recovered through the inverse transformation. Dynamics of some obtained solutions was analyzed and illustrated. Here, particularly, we would like to sum up a typical characteristic of limit solitons: the wave trajectories asymptotically follow logarithm curves combined with linear functions. This point is based on several examples we have examined[43]. Obviously, through the Galilean transformation 4., all these obtained solutions of the mkdv equation.3 can easily be used for the KdV-mKdV equation.4 which often appears in physics contexts. In fact, in the paper we do not differ them from each other. In addition to the KdV-mKdV equation, our treatment to the complex operation in 6. can also be applied to the sine-gordon equation. There are Miura transformations between the KdV equation and the mkdv equation for both ε = ±. For the mkdvε = equation the Miura transformation provides a real map between solutions of the KdV equation and the mkdv equationε =, c.f.[0]. However, when ε =, i.e., for the mkdv equation.3, the Miura transformation has to be complex. In more detail, it maps the real mkdv equation.3 to a complex KdV equation. That means we have had a nice determinant expression for the complex KdV equation. Further investigation about this will be considered separately. 6. List of solutions Let us list out the obtained solutions and their corresponding basic Wronskian vectors. Solutions to the mkdv equation.3 can be given by where v = arctan F F = F,xF F F,x x F + F, 6.4a f = fϕ = N = F + if, F = Re[f], F = Im[f]. 6.4b The available Wronskian vectors are the following. For soliton solutions: with ϕ = ϕ [s] N = ϕ, ϕ,, ϕ N T, 6.5a ϕ j = a + j eξ j + ia j e ξ j, ξ j = k j x 4k 3 j t + ξ 0 j, a + j, a j, k j, ξ 0 j R. 6.5b 8

29 For limit solutions of solitons: ϕ = ϕ [ls] N k = A + Q ia Q 0, A± G N R, 6.6a with Q ± 0 = Q± 0,0, Q± 0,,, Q± 0,N T, Q ± 0,s = s! s k e ±ξ, 6.6b where ξ is defined in 6.5b. For breather solutions: ϕ = ϕ [b] N = ϕ, ϕ, ϕ, ϕ,, ϕ N, ϕ N T, 6.7a with For limit solutions of breathers: ϕ j = a j e ξ j + b j e ξ j, ϕ j = ā j e ξ j b j e ξ j, 6.7b ξ j = k j x 4k 3 j t + ξ 0 j, a j, b j, ξ 0 j C. 6.7c ϕ = ϕ [lb] N k = ϕ +,, ϕ,, ϕ+,, ϕ,,, ϕ+ N,, ϕ N, T, 6.8a and the elements are given through ϕ + = ϕ +,, ϕ+,,, ϕ+ N, T = AQ BQ 0, ϕ = ϕ,, ϕ,,, ϕ N, T = Ā Q + 0 B Q 0, 6.8b 6.8c where A, B G N C, and ξ is defined in 6.7c. Q ± 0 = Q± 0,0, Q± 0,,, Q± 0,N T, Q ± 0,s = s! s k e ±ξ, 6.8d We note that, thanks to the linear property of the CES.8, one may also get mixed solutions by arbitrarily combining the above vectors to be a new Wronskian vector. For example, take ϕ [s] ϕ = N ϕ [ls]. 6.9 N k N + The related solution corresponds to the interaction between N -soliton and a N -order limit-soliton solutions. Finally, for the rational solution, it is given by vx, t = v 0 F,XF F F,X F + F, X = x 6v0t, v 0 0 R, 6.0a where still f = fψ = N = F + if, F = Re[f], F = Im[f], 6.0b 9

30 and the Wronskian is composed by ψ = ψ, ψ,, ψ N T, 6.a with and ψ j+ = j j! j k ϕ j = 0,,, N, 6.b k =0, ϕ = v 0 + ik e η + v 0 ik e η, η = k X 4kt, 3 k R. 6.c Acknowledgments The authors sincerely thank Prof. Gesztesy for kindly providing Refs.[0, 3]. This project is supported by the NSF of China No. 0757, Specialized Research Fund for the Doctoral Program of Higher Education of China No , Shanghai Leading Academic Discipline Project No. J500 and Postgraduate Innovation Foundation of Shanghai University No. SHUCX07. A Proof of Theorem Proof. The compatibility of.a and.b, i.e., φ xt = φ tx, yields.3b. Using.a one gets the complex conjugate form of f as f = Bt, N. A. Then the necessary derivatives of f and f are presented as the following, f x = N, N, A.a f xx = N 3, N, N + N, N +, A.b f xxx = N 4, N, N, N + N, N + + N 3, N, N +, A.c f t = 4 N 4, N, N, N N 3, N, N + + N, N + + trct N, A.d and f x = Bt, N 3, N, f xx = Bt, N 4, N, N +, N 3, N, f xxx = Bt, N 3, N + +, N 4, N, N A.3a A.3b +, N 5, N 3, N, N, A.3c f t = 4 Bt, N 5,, N 3, N, N, N 4, N, N +, N 3, N + + Bt t, N + trct Bt, N. A.3d 30

31 Using the condition.a the complex conjugate of A.d is f t = 4 Bt, N 5, N 3, N, N, N 4, N, N +, N 3, N + + tr Ct Bt, N, which should be same as A.4. This requires Bt t = 0 and trct Rt, i.e., the condition.3a. Noting that φ xx = Bt Btφ and using Proposition with Ω j,s = x, we find trbt Bt, N 3, N =, N 5, N 3, N, N +, N 3, N +, A.4 trbt Bt N, N = N 4, N, N, N + N, N +, trbt Bt, N =, N 4, N, N +, N 3, N, trbt Bt N = N 3, N, N + N, N +. Then, substituting A. and A.3 into.0a and making use of A.7, we have f t f ff t + f xxx f 3 f xx f x + 3 f x f xx ff xxx = 6 Bt, N 3, N + N, N 4, N, N N, N + N 4, N, N, N, N N 3, N, N +, N + N, N +, N 3, N, N 4, N, N N, A.5 A.6 A.7 which is zero in the light of Proposition. Similarly, one can prove.0b. Thus the proof is completed. B Eigen-polynomial of A = B B We prove the Proposition 3 through the following Lemmas. Lemma B.. For two arbitrary Nth-order complex matrices A and B, where I N is the Nth-order unit matrix. detλi N AB = detλi N BA, B. Proof. Assuming that ranka = r, then there exist N th-order non-singular matrices P and Q such that Ir 0 A = P Q. B. 0 0 Thus Ir 0 detλi N AB = detλi N P 0 0 = detp Ir 0 λi N P 0 0 = detλi N P Ir 0 P 0 0 Ir 0 = detλi N QB QBP. QBP QBP

32 In a similar way we have detλi N BA = detλi N QBP Ir B.3 If we rewrite the matrix QBP into the following block matrix form with same structure as Ir 0, 0 0 then we have Ir which further means B B QBP = 0 0 B B QBP = B B Ir 0, QBP 0 0, B.4 = B 0 B 0, B.5 detλi N AB = λ N r detλi r B = detλi N BA. B.6 We complete the proof. Lemma B.. Assuming that B is an arbitrary Nth-order complex matrix and B is its complex conjugate, then detλi BB is a polynomial of λ with real coefficients. Proof. Write fλ = detλi N BB = a N λ N + + a λ + a 0. B.7 Then using Lemma B. with A = B we have fλ = detλi N BB = detλi N B B = det λi N BB = f λ, B.8 which means all the coefficients {a j } are real. C Discussions on the trivial solutions to the CES.8 Let the square matrices A, B follow the relation A = B B. C. We start from the case. Noting that the product of all the eigenvalues of A is nonnegative, in the following we first look at k A = 0 0 k, k k 0, k, k R, C. and suppose a b B = c d C.3 3

33 with undetermined a, b, c, d C. However, in this case it can be found that the matrix relation C. does not have any solutions unless k = k. So next we turn to consider A = k 0 0 k, k R. C.4 In this case, the equation C. admits a non-diagonal matrix solution B as d d +k B = c c d e iθ, c 0, c, d, θ R. C.5 We note that such a B does not lead to any nontrivial solutions to the mkdv equation. In fact, in the CES.8, the general solution to the equation set.8a and.8c is ϕ e iξ ϕ = = H C.6a with arbitrary matrix H C and ϕ e iξ ξ = kx + 4k 3 t + ξ 0, k, ξ 0 R. C.6b However, no matter what condition the matrix H should satisfy under the equation.8b, the Wronskian fϕ = H fe iξ, e iξ T is always a constant, which leads to a trivial solution to the mkdv equation. In the case of the N N matrix N is even A = diag k, k,, k, k R, similar to C.6, the general solution to.8a and.8c is φ = H e iξ, e iξ, 0, 0,, 0 T, H C N N, which leads to a zero Wronskian fϕ. With these discussions we can conclude that for the CES.8 with an N N matrix A = diagα, α,, α, α R, C.7 the possible solution ϕ to.8 composes a trivial Wronskian fϕ. D Proof for Theorem Proof. The parameter v 0 will lead to complicated expressions for f and its derivatives. For simplification let us introduce the notation j where the subscript j indicates the absence of the φ j column[], for example, N j = φ 0,, φ j, φ j+,, φ N. 33

34 Derivatives of f have already given in A.. For f, using the CES 4.6 we can reach f = B t N v 0 j i N j N j, f X = B t N v 0 j+ i N j N N j + v 0 j i N j N, N + j, f XX = B t N v 0 N N, N + + v 0 j i N j N, N + j N + N v 0 j i N j N, N, N + j + N + v 0 j+ i N j N, N + j, v 0 j+ i N j N j f XXX = B t N 3 3 v 0 j+ i N j N, N + j + v 0 N N, N + N + 3 v 0 j+ i N j N, N + j + v 0 N N 3, N, N + N + N v 0 j i N j N, N + 3 j + v 0 j i N j N, N, N + j N v 0 j+ i N j N, N, N + j + v 0 N i N 3, N, N + N 3 + N 3 v 0 j i N j N 3, N, N, N + j + v 0 j+3 i N j 3 N j, f t = 4 B t N N 3 + N 3 v 0 j i N j N, N + 3 j + v 0 j+3 i N j 3 N j v 0 j i N j N 3, N, N, N + j + v 0 N N, N + v 0 j i N j N, N, N + j v 0 N i N 3, N, N + N v 0 N N 3, N, N + + trct B t N v 0 j i N j N j. 34

35 Using the condition 4.6 the complex conjugate of A.d should be same as the above f t. This requires trct Rt, i.e., the condition 4.7b. Besides, noting that φ XX = Bt Bt v 0 I Nφ and using Proposition with Ω j,s = X, we can have the following identities: trbt Bt v 0I N N j = N, N + j N, N, N + j N j, j = 0,,, N, trbt Bt v 0I N N, N + j = N, N + 3 j N 3, N, N, N + j N, N + j, j = 0,,, N 3, trbt Bt v 0I N N = N 3, N, N + N, N +, trbt Bt v 0I N N, N = N 4, N, N, N + N, N +. With these results and A. in hand, for 4.4a we have f t f ff t + f XXX f 3 f XX f X + 3 f X f XX ff XXX =6 v 0 N i N 3, N, N + N N, N N 3, N, N + =0. + N, N + N 3, N, N Similarly, one can prove 4.4b. References [] D.J. Zhang, Notes on solutions in Wronskian form to soliton equations: Korteweg de Vriestype, arxiv:nlin.si/ [] R.M. Miura, C.S. Gardner and M.D. Kruskal, KdV equation and generalizations. II. Existence of conservation laws and constant of motion, J. Math. Phys [3] C.S. Gardner, J.M. Green, M.D. Kruskal and R.M. Miura, Method for solving the Kortewegde Vries equation, Phys. Rev. Lett [4] H. Ono, Soliton fission in anharmonic lattices with reflectionless inhomogeneity, J. Phys. Soc. Jpn [5] T. Kakutani and H. Ono, Weak non-linear hydromagnetic waves in a cold collision-free plasma, J. Phys. Soc. Jpn [6] K. Konno and Y.H. Ichikawa, A modified Korteweg de Vries equation for ion acoustic waves, J. Phys. Soc. Jpn [7] S. Watanabe, Ion acoustic soliton in plasma with negative ion, J. Phys. Soc. Jpn [8] K.E. Lonngren, Ion acoustic soliton experiment in a plasma, Opt. Quant. Electron