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

Transcription

1 Integrable dispersionless KdV hierarchy with sources arxiv:nlin/ v [nlin.si] 7 Jun 2006 Zhihua Yang, Ting Xiao and Yunbo Zeng Department of Mathematical Sciences, Tsinghua University, Beijing 00084, P. R. China Abstract. An integrable dispersionless KdV hierarchy with sources dkdvhws) is derived. Lax pair equations and bi-hamiltonian formulation for dkdvhws are formulated. Hodograph solution for the dispersionless KdV equation with sources dkdvws) is obtained via hodograph transformation. Furthermore, the dispersionless Gelfand-Dicey hierarchy with sources dgdhws) is presented. PACS number: I Introduction In recent years, research in the dispersionless hierarchies has become quite active see, for example, []-[] and references therein). Dispersionless hierarchies arise as the quasiclassical limit of the original dispersionfull hierarchies [2]. The operators in the Lax equations for dispersionfull hierarchy are replaced by phase functions for dispersionless hierarchy, commutators are replaced by Poisson bracets and the role of Lax pair equations by dispersionless Lax pair equations. The dispersionless hierarchies have Hamiltonian formulation [, 3] and many other aspects [4, 5, 6], and several methods of solutions of dispersionless hierarchies have been formulated [7, 8, 9, 0, ]. The soliton equations with self-consistent sources SESCS) are another type of integrable models and have important physical applications [2]-[28]. There are some ways to derive the SESCS, for example, the Mel niov way [2, 8, 9, 2, 30] and the Leon approach [3, 4, 5]. In recent years, SESCS were studied based on the constrained flows of soliton equations which are just the stationary equations of SESCSs[6, 7]. There are several methods for solving the SESCSs, for example, the inverse scattering method [8, 9, 20], the matrix theory [2], the method and gauge transformation [22, 23], and yangzhihua@tsinghua.org.cn

2 the Darboux transformation [7, 24, 25, 26]. By treating the variable x as the evolution parameter and t as the spatial variable respectively, and introducing Jacobi-Ostrogradisi coordinates, the SESCS has a t-type Hamiltonian formulation [27, 28]. This paper is devoted to the integrable dispersionless KdV hierarchy with sources dkdvhws). Considering the asymptotic expansion of the wave function, and taing the dispersionless limit of the KdV hierarchy with sources KdVHWS), we can deduce the dkdvhws. The Lax pair equations of the dkdvhws can be deduced by the dispersionless limit of the Lax pair equations of the KdVHWS. Similar to the dkdv hierarchy, the dkdvhws has bi-hamiltonian formulation and can be solved via hodograph transformation. The Gelfand-Dicey hierarchy with sources GDHWS) is the integrable generalization of the Gelfand-Dicey hierarchy, and the corresponding integrable dispersionless hierarchy, i.e. the dgdhws, can be deduced. This paper is organized as follows: in Section 2 we review some definitions and results about the KdVHWS. In Section 3 we derive the dkdvhws as well as its Lax pair equations by taing the dispersionless limit of the KdVHWS. In Section 4, we construct the bi-hamiltonian formulation of dkdvhws. Then we derive the hodograph solution for dispersionless KdV equation with sources dkdvws) in Section 5. In Section 6, we deduce the integrable dispersionless Gelfand-Dicey hierarchy with sources dgdhws). Some conclusion is made in Section 7. 2 The KdV hierarchy with sources We first briefly review some definitions and results about the KdV hierarchy with sources in the framewor of Sato theory. Given a pseudo-differential operator PDO) of the form [29] L = 2 + u 2.) where =, u = ux, t), t = t x 3, t 5,...), and the wave function ψ = ψx, t), consider the Lax pair Lψ = λψ, 2.2a) ψ t2m = B 2m ψ, 2.2b) where λ is a parameter, B 2m = L 2m 2 ) +, m =, 2, 3,..., and A) + here standing for the differential part of A. The compatibility condition of 2.2a) and 2.2b) gives rise to L t 2m = [B 2m, L], 2.3) 2

3 which is the well-nown KdV hierarchy[29]. As was shown in [29], the KdV hierarchy could be written as bi-hamiltonian systems u t2m = B 0 δh 2m+ δu δh 2m = B, m =, 2,..., 2.4) δu where H 2m = h 2m dx is a functional of u, B 0 =, B = u + u 2 x are Hamiltonian operators, and δh 2m is the Euler-Lagrange derivative of the Hamiltonian δu H 2m defined as Let us consider [6] δh 2m δu = δ δu B 2m = B 2m + h 2m dx = =0 ) h 2m u ). 2.5) ψ φ, m =, 2, 3,..., 2.6) where ψ = ψ x, t) and φ = φ x, t) satisfying Lψ = λ ψ, L φ = λ φ, =,..., n, 2.7) here L = ) 2 + u = L is the adjoint operator of L, and λ is a constant, =,..., n. Then the KdV hierarchy with sources KdVHWS) [6, 20, 24, 30] can be defined as L t 2m = [ B 2m, L], 2.8a) Lψ = λ ψ, 2.8b) with the Lax pair given by L φ = λ φ Lψ = λψ, 2.8c) 2.9a) ψ t2m = B 2m ψ, 2.9b) namely, under 2.8b) and 2.8c), the compatibility condition of 2.9a) and 2.9b) gives rise to 2.8a). 3

4 3 Dispersionless limit Following the procedure introduced in [2, 3, 9, 0], we could derive the dispersionless hierarchy by taing the dispersionless limit of the initial system. Taing T = ǫt, X = ǫx, and thining of ux, t) = u X, T ) = UX, T) + Oǫ) as ǫ 0, L in 2.) changes into ǫ ǫ L ǫ = ǫ 2 X 2 + ux ǫ, T ǫ ) = ǫ2 X 2 + UX, T) + Oǫ), 3.) where X =. It can be proved [3] that X satisfies L = σ ǫ L ǫ ) = p 2 + U 3.2) L T2m = {B 2m, L}, 3.3) where σ ǫ denotes the principal symbol [2], the bracet {, } is the Poison bracet defined in 2D phase space p, X) as {A, B} = A B p X A B X p, 3.4) and B 2m = L 2m 2 ) + now refers to nonnegative powers of p. We define 3.3) as the dispersionless KdV dkdv) hierarchy [3], and the first few equations are expressed as U T = U X, 3.5a) U T3 = 3 2 UU X, 3.5b) U T5 = 5 8 U2 U X, 3.5c) U T7 = 35 6 U3 U X, d) As was shown in [3], equation 3.5b) has tri-hamiltonian formulation as where δh 5 U T3 = D = D δh 3 2 = 3 4 D δh 3, 3.6) D = 2 X, D 2 = U X + X U, D 3 = U 2 X + X U 2, 4

5 H = UdX, H 3 = 4 U 2 dx, H 5 = 8 U 3 dx. For the general case, define H 2m = 2 2m TrL2m 2, 3.7) where TrA = ResAdX, and ResA is the residue of the general Laurent polynomial of the form A = + i= a ix)p i, i.e. the coefficient of the p term, then the Hamiltonians 3.7) are in involution with respect to any of the three Poisson bracets {H 2m, H 2l } i = dx δh 2m D δh 2l i = 0, i =, 2, 3, 3.8) and dkdv hierarchy 3.3) have tri-hamiltonian formulation U T2m = D δh 2m+ = D 2 δh 2m = 2m )2m 3) δh 2m 3 D 2m 2) 2 3, m = 2, 3, ) It was also shown in [3] that the solution of 3.5b) can be described through the implicit form where f is an arbitrary function. U = fx UT 3), 3.0) In what follows, we derive the dispersionless KdV hierarchy with sources dkdvhws). By taing T = ǫt, X = ǫx, 2.8) change into ǫl ǫt2m = [B ǫ2m ) + ψ X ǫ, T ǫ )ǫ X) φ X ǫ, T ǫ ), L ǫ], L ǫ ψ X ǫ, T ǫ ) = λ ψ X ǫ, T ǫ ), 3.a) 3.b) L ǫ φ X ǫ, T ǫ ) = λ φ X ǫ, T ǫ ), 3.c) where B ǫ2m ) = L 2m 2 ǫ ) +. Similar to the dispersionless KP case in [2, 9, 0], we consider the following WKB asymptotic expansion of ψ X, T ) and φ ǫ ǫ X, T ), =, 2,..., n, ǫ ǫ ψ X ǫ, T ǫ ) exp{sx, T, λ = λ ) ǫ 5 + β + Oǫ)}, ǫ 0, 3.2a)

6 φ X ǫ, T ǫ ) exp{ SX, T, λ = λ ) + β 2 + Oǫ)}, ǫ 0, 3.2b) ǫ it can be calculated that as ǫ 0. Therefore, we have ψ X ǫ, T ǫ )ǫ X) φ X ǫ, T ǫ ) = e β +β 2 [ǫ X ) + X SX, T, λ = λ ) + Oǫ))ǫ X ) 2 + X SX, T, λ = λ )) 2 + Oǫ))ǫ X ) ], σ ǫ ψ X ǫ, T ǫ )ǫ X) φ X ǫ, T ǫ ) ) = v p p, 3.3) where v = e β +β 2, p = X SX, T, λ = λ ). 3.4) Taing the principal symbol of both sides of 3.a), we have L T2m = {B 2m + v p p, L} = {B 2m, L} + { and the dispersionless limit of 3.b) and 3.c) lead to v p p, L}, 3.5) p 2 + U = λ, X p v ) = ) Under 3.6), it can be found that { v p p, L} = 2v,X, 3.7) therefore, the dispersionless limit of 3.), i.e. dkdvhws, reads U T2m = {B 2m, L} 2 v,x, m =, 2, 3,..., p 2 + U = λ, X p v ) = a) 3.8b) 3.8c) 6

7 Integrating 3.8c) and taing λ m as the integral constants, then we can eliminate v in 3.8a) and rewrite dkdvhws in another form U T2m = {B 2m, L} 2 λ m ) λ U X. 3.9) If we tae the dispersionless limit of 2.9), we will obtain the Lax pair equations of dkdvhws 3.8) as p 2 + U = λ, 3.20a) p T2m = B 2m + v, 3.20b) p p )X namely, under 3.8b) and 3.8c), the compatibility condition of 3.20a) and 3.20b) gives rise to 3.8a). We can eliminate v and p in 3.20b) and rewrite 3.20) in another form p 2 + U = λ, 3.2a) p T2m = B 2m + which are the Lax pair equations of 3.9). λ m ) p, 3.2b) λ U λ U) X We give two examples in the following, the first one is the dispersionless KdV equation with sources dkdvws) U T3 = 3 2 UU X 2 λ 2 ) λ U X, 3.22) with Lax pair equations given by p 2 + U = λ, 3.23a) p T3 = p Up + λ 2 ) p. 3.23b) λ U λ U) X And the second example is the dispersionless KdV5) equation with sources dkdv5)ws) U T5 = 5 8 U2 U X 2 λ 3 ) λ U X, 3.24) 7

8 with Lax pair equations given by p 2 + U = λ, 3.25a) p T5 = p Up U2 p + λ 3 ) p. 3.25b) λ U λ U) X 4 Hamiltonian formulation of dkdvhws It is well nown that the KdV hierarchy has bi-hamiltonian formulation [29], the dkdv hierarchy has tri-hamiltonian formulation [3], or even further, quasi-hamiltonian formulation [], and the KdVHWS has bi-hamiltonian formulation [27]. Motivated by the Hamiltonian formulation of the dkdv case [3], we would construct the bi-hamiltonian formulation of the dkdvhws 3.9). Let us firstly consider dkdvws 3.22). Set A = 2λ 2 λ U dx, B = 2λ λ U dx, 4.) then by direct computation we have δa D = 2 λ 2 ) λ 2 ) X = 2, λ U λ U X δb D 2 = U λ ) X + X U) = λ2 U X λ U λ U) = 2 λ 2 ) 3/2 λ U Therefore, if we denote H 3 = H 3 + B = dx 4 U2 + 2 λ λ U ), H 5 = H 5 + A = dx 8 U3 + 2 then equation 3.22) can be written in two Hamiltonian forms For the dkdvhws 3.9), denote H 2m = H 2m + 2 X. 4.2a) λ 2 λ U ), 4.2b) δ U T3 = D H 5 = D δ H ) 8 λ m λ U dx, 4.4)

9 we can directly prove see Appendix) that the Hamiltonians H 2m, m =, 2,... satisfy { H 2m, H 2l } i = dx δ H 2m D δ H 2l i = 0, i =, 2, 4.5) therefore, the dkdvhws 3.9) have bi-hamiltonian formulation U T2m = D δ H 2m+ δ = D H 2m 2, m =, 2, ) 5 Hodograph solution for dkdvws In this section, using the hodograph transformation[7, 9, 0], we will derive the hodograph solution for the the dkdvws 3.22). Following [7] and let U T3 = BU)U X, we can find from 3.22) that BU) = 3 2 U λ 2 λ U) 3/2. 5.) Maing the hodograph transformation with the change of variables X, T 3 ) U, T 3 ) and let X = XU, T 3 ), we have which implies that It can be integrated as 0 = dx = X U + X = X dt 3 U T 3 T 3 U BU X + X, 5.2) T 3 X = B = 3 T 3 2 U + λ 2 λ U) 3/2. 5.3) 3 X + 2 U λ 2 ) T λ U) 3/2 3 = FU), 5.4) where FU) is an arbitrary function of U, 5.4) gives an implicit solution of 3.22). When F = 0 and λ = λ 2 = = λ n = 0, dkdvws 3.22) degenerates to dkdv equation 3.5b), and 5.4) degenerates to the rational solution of 3.5b), U = 2X 3T. When F 0 and is convertible, then solution of dkdvws 3.22) can be written through the implicit form which is similar to 3.0). 3 U = F X + 2 U λ 2 ) T λ U) 3/2 3 ), 5.5) 9

10 6 The dispersionless Gelfand-Dicey hierarchy with sources The well-nown Gelfand-Dicey hierarchy with sources GDHWS) [6] is defined as L = [ t B m, L] = [B m + ψ φ, L], 6.a) m Lψ = λ ψ, 6.b) where L φ = λ φ, 6.c) L = N + u N 2 N u + u 0, 6.2) u = u N 2,...,u 0 ) T, u i = u i x, t), i = 0,,..., N 2, t = t 2, t 3,...), B m = [L N ) m ] +, L = ) N + ) N 2 u N )u + u 0 is the adjoint operator of L, λ is a constant, =,..., n, ψ = ψ x, t) and φ = φ x, t), and the Lax pair is given by Lψ = λψ, 6.3a) ψ tm = B m ψ, 6.3b) namely, under 6.b) and 6.c), the compatibility condition of 6.3a) and 6.3b) gives rise to 6.a). Following the procedure given above, we can derive the dispersionless Gelfand-Dicey hierarchy with sources dgdhws). Taing T = ǫt, X = ǫx, and let u X, T ) = U ǫ ǫ X, T)+ Oǫ) as ǫ 0, then 6.) change into ǫl ǫtm = [B ǫm + ψ X ǫ, T ǫ )ǫ X) φ X ǫ, T ǫ ), L ǫ], 6.4a) L ǫ ψ X ǫ, T ǫ ) = λ ψ X ǫ, T ǫ ), 6.4b) where L ǫφ X ǫ, T ǫ ) = λ φ X ǫ, T ǫ ), L ǫ = ǫ X ) N + u N 2 X ǫ, T ǫ )ǫ X) N u X ǫ, T ǫ )ǫ X + u 0 X ǫ, T ǫ ) 6.4c) = ǫ X ) N + U N 2 X, T) + Oǫ))ǫ X ) N U X, T) + Oǫ))ǫ X + U 0 X, T) + Oǫ), 6.5) 0

11 and B ǫm = [L N ǫ ) m ] +. Consider the following WKB asymptotic expansion of ψ X, T ) ǫ ǫ and φ X, T ), =,..., n, ǫ ǫ ψ X ǫ, T ǫ ) exp{sx, T, λ = λ ) ǫ + β + Oǫ)}, ǫ 0, 6.6a) φ X ǫ, T ǫ ) exp{ SX, T, λ = λ ) + β 2 + Oǫ)}, ǫ 0, 6.6b) ǫ then the principal symbol [2] of 6.4a) arises where L T m = {B m + v p p, L} 6.7) L = σ ǫ L ǫ ) = p N + U N 2 p N U p + U 0, 6.8) B m = σ ǫ B ǫm ), and v = e β +β 2, p = X SX, T, λ = λ ) are obtained by σ ǫ ψ X ǫ, T ǫ )ǫ X) φ X ǫ, T ǫ ) ) = v p p. 6.9) The dispersionless limit of 6.4b) and 6.4c) give rise to p N + U N 2p N U p + U 0 = λ, 6.0a) v L ) = 0, 6.0b) X p p=p 6.7) together with 6.0a) and 6.0b) give rise to the dispersionless Gelfand-Dicey hierarchy with sources dgdhws) L T m = {B m + v p p, L}, 6.a) whose Lax pair equations are given by p N + U N 2 p N U p + U 0 = λ, 6.b) v L ) = 0, 6.c) X p p=p p N + U N 2 p N U p + U 0 = λ, 6.2a)

12 p Tm = B m + v p p ) X, 6.2b) namely, under 6.b) and 6.c), the compatibility condition of 6.2a) and 6.2b) gives rise to 6.a). Under 6.b) and 6.c), it can be found by a tedious computation that { v p p, L} = a N 2 p N 2 + a N 3 p N a 0, 6.3) where a N 2 = Nv,X, a N 3 = Nv p ) X, and for i = 4,...,N, i 3 a N i = v j= i 2 jp j U N+ i+j,x l=2 N l)v p i 2 l When N = 2, we have L = p 2 + U, and 6.3) is the same as 3.7). ) X U N l Np i 2 v ) X. 6.4) Similar to the dkdvhws, the dgdhws possess bi-hamiltonian formation and their solutions can be obtained via hodograph transformation. 7 Conclusion We derive dkdvhws by taing the dispersionless limit of KdVHWS, meanwhile, Lax pair equations of the dkdvhws can be obtained by taing the dispersionless limit of the corresponding dispersionfull equations. We have constructed the bi-hamiltonian formulation of the dkdvhws, and have obtained the implicit solutions of the dkdvws via the hodograph transformation. For the generalization case, we have deduced the dgdhws which also possess bi-hamiltonian formulation and can be solved via the hodograph transformation. Acnowledgment This wor was supported by the Chinese Basic Research Project Nonlinear Science. 2

13 Appendix Here we give the proof of involution relation of the Hamiltonians H 2m 4.4). Since { H 2m, H 2l } = { H 2m, H 2l 3 } 2 = { H 2l 3, H 2m } 2 = { H 2l 3, H 2m+ } = { H 2m+, H 2l 3 } = = { H 2m+2l 3, H }, it suffices to prove that for any m, { H 2m, H } = 0. We can directly calculate that { H 2m, H } δ H2m = D δ H dx δh2m = δh2m = λ m λ m λ U )2 X ) λ U l= ) dx λl U U X λl Uλ l U) dx. Since all δh 2m, m =, 2,... are of the form cu s, where c are constants, and s N, it suffices to prove U s U X dx = 0, l =,..., n, I) λl Uλ l U) U X λ U dx = 0,, l =,...,n. II) λ l Uλ l U) For I), let F l = λ l U, then U = λ l Fl 2, U X = 2F l F l,x and U s U X λl Uλ l U) dx = 2 λl F 2 which are all zero for the reason that λ l F 2 l )s F 2 l F 2 l l )s df l, are rational polynomials of F 2 l and so λ l Fl 2)s df Fl 2 l are total differentials. For II), when = l, then II) obviously holds; when l, let G l = λ l U, then 3

14 U X = G l,x, and = = = U X λ U λ l Uλ l U) dx G l,x λ λ l + G l Gl G l dx G l,x G 2 λ λ l l λ λ l dx G l + λ λ l G l + d λ λ l G l ), set λ λ l G l = tan 2 H l, then + λ λ l G l = sec 2 H l, d λ λ l G l ) = 2tanH l sec 2 H l dh l, and d λ λ l ) = 2tanH lsec 2 H l dh l = 2dsecH l ), λ λ l G l + G l sech l which are total differentials, and this proves II). References [] Brunelli J C 2000 Braz. J. Phys [2] Taasai T and Taebe T 995 Rev. Math. Phys [3] Brunelli J C 996 Rev. Math. Phys [4] Chen Y T and Tu M H 2003 Lett. Math. Phys [5] Zheltuhin K 2002 Phys. Lett. A [6] Alagesan T Chung Y and Naeeran K 2004 Chaos, Solitons and Fractals 2 63 [7] Kodama Y 988 Phys. Lett. A [8] Kodama Y 989 Phys. Lett. A

15 [9] Xiao T and Zeng Y B 2006 Phys. Lett. A [0] Xiao T and Zeng Y B 2006 J. Nlin. Math. Phys. at press [] Chen Y T and Tu M H 2003 J. Phys. A: Math. Gen [2] Mel niov V K 990 J. Math. Phys [3] Leon J988 J. Math. Phys [4] Leon J and Latifi A 990 J. Phys. A [5] Leon J 990 Phys. Lett. A [6] Antonowicz M 992 Phys. Lett. A [7] Zeng Y B 994 Phys. D 73 7 [8] Mel niov V K 988 Phys. Lett. A [9] Mel nilov V K 992 Phys. Lett. A 8 33 [20] Lin R L, Zeng Y B and Ma W X 200 Phys. A [2] Mel niov V K 989 Comm. Math. Phys [22] Dotorov E V and Shchesnovich V S 995 Phys. Lett. A [23] Shchesnovich V S and Dotorov E V 996 Phys. Lett. A [24] Zeng Y B, Ma W X and Shao Y S 200 J. Math. Phys [25] Zeng Y B, Shao Y J and Xue W M 2003 J. Phys. A: Math. Gen [26] Zeng Y B, Shao Y J 2005 J. Phys. A: Math. Gen [27] Blasza M 995 J. Math. Phys [28] Zeng Y B 999 Phys. A [29] Dicey L A 2002 Soliton equations and Hamiltonian systems 2nd ed.)world Scientific: Singapore) [30] Mel niov V K 990 Nonlinear evolution equations and dynamical systems Berlin: Springer) 5