arxiv:math-ph/ v2 29 Mar 2003

Transcription

1 THE ATIYAH-HITCHIN BRACKET AND 1D INTEGRABLE SYSTEMS arxiv:math-ph/ v2 29 Mar 2003 KL Vaninsky Michigan State University, East Lansing, MI Abstract All fashionable integrable equations analyzed by the inverse spectral transform are Hamiltonian systems We demonstrate that the Hamiltonian formalism is intrinsically build into the spectral theory The Poisson bracket on the phase space is an image of the Atiyah Hitchin bracket on Weyl functions under the inverse spectral transform 1 Introduction All 1 D partial differential equations like Korteweg-de-Vriez, Camassa Holm, sin/sinh Gordon, cubic nonlinear Schrödinger equation, analyzed by the inverse spectral transform are Hamiltonian systems We consider these problems on the entire line, ie x R 1 We do not assume anything from the initial data except continuity The KdV equation q t = 3 q x q x 3 is a Hamiltonian flow with the Gardner Faddeev Zakharov bracket (1) {A, B} = δa δq(x) D δb δq(x) dx, Recently discovered the Camassa Holm equation [ v t + v v x + x G v ( ) ] 2 v = 0 2 x D = x is Hamiltonian with the bracket (2) {A, B} = δa δb (md + Dm) δm(x) The work is partially supported by NSF grant DMS δm(x) dx, 1 Typeset by AMS-TEX

2 2 KL VANINSKY where m = v 2 v x 2 The sin/sinh Gordon equation 2 q t 2 q + sin / sinhq = 0 x2 is Hamiltonian with the classical bracket δa (3) {A, B} = δq(x) δb δp(x) where p = q t Finally, the cubic NLS equation δa δp(x) i t ψ = 2 x ψ + 2 ψ 2 ψ, where ψ(x, t) is a complex function can be written as with the Hamiltonian H = 1 2 (4) {A, B} = 2i ψ = {ψ, H}, δb δq(x) dx, ψ 2 + ψ 4 dx and the bracket δa δψ(x) δb δψ(x) δa δψ(x) δb δψ(x) dx The NLS flow will serve as our model example The equation is a compatibility condition for the zero curvature representation with and V 2 = iλ 2 σ 3 + Y 0 = [ t V 3, x V 2 ] = 0, ( ) iλ 2 0 iλ ( ) 0 ψ ψ 0 V 3 = λ2 2 iσ 3 λy 0 + ψ 2 iσ 3 iσ 3 Y 0 The auxiliary spectral problem f = V 2 (x, λ)f for the vector function f T (x, λ) = (f 1, f 2 ) is, in fact, the eigenvalue problem for the Dirac operator (5) Df = [( ) ( )] iψ i 0 1 x + f = λ iψ 0 2 f We will demonstrate that the Poisson bracket (4) is intimately related with the spectral problem (5) for the Dirac operator Namely, the Weyl functions of the

3 THE ATIYAH-HITCHIN BRACKET 3 Dirac operator cary natural Poisson structure We call it the Atiyah Hitchin bracket It is given by the formula (6) {w(λ), w(µ)} = (w(λ) w(µ))2 λ µ The bracket (4) on the phase space is an image under the inverse spectral transform of the Atiyah Hitchin bracket The main result of this note is that all other brackets (1), (2) and (3) arise in the same way from the single formula (6) This can be proved along the same lines as it will be done for NLS The question of construction of canonical coordinates can not be resolved is such generality It requires assumptions on the spectrum of the operator in question At the moment the only equations which were analyzed are the finite Toda lattice, [V1], and the Camassa Holm equation with rapidly decaying initial data of one sign, [V2] We conclude the introduction with two historical remarks Various approaches to relate the Hamiltonian theory and spectral theory were developed in the previous years One is the r matrix approach developed by Leningrad group, [FT] It allows to put computation of the bracket between the entries of the monodromy matrix in compact form Unfortunately the r matrix is different for different integrable equations It requires significant efforts to find the r matrix for each integrable system Within our approach the Atiah Hitchin bracket gives a universal form of the bracket at least for all fashionable integrable PDE associated with differential operator of the second order In the recent paper Krichever and Phong [KP] proposed a new approach for construction of the symplectic forms for integrable equations with periodic initial data The approach of [KP] works for 2D systems like KP and 2D-Toda lattice, but at the same time requires periodicity of the potential It can be extended to the scattering case, [V3] The regularity of the potential is needed to prove that Jost solutions can be expanded in power series Our approach does not require anything from the potential except for the operator to be in the limit point case 2 Statement of the Theorem Pick any real number y and consider the Dirac operator acting on smooth functions with support on the positive half line [y, + ) The case of negative half line can be considered in the same way The Weyl solution H T (x, λ) = (H 1, H 2 ), [W], is the solution of (5) which belongs to the L 2 ([y, + )) It is well known, [LS], that for a continuous potential ψ(x) and λ with Iλ 0 there exists a unique (up to multiplication) solution with this property Pick some α To construct H consider the fundamental system of solutions ϕ α (x, y, λ) and θ α (x, y, λ) of (5) normalized as ( ) e iα (7) ϕ α (x, y, λ) x=y = e iα ( ie iα θ α (x, y, λ) x=y = ie iα )

4 4 KL VANINSKY The Weyl solution H is a linear combination of ϕ and θ: H = bϕ α + aθ α Evidently the Weyl function w α (λ, y) = a/b is defined correctly and h α (x, y, λ) = ϕ α + w α θ α is square integrable for λ with Iλ 0 As we will show any two functions w α and w β are related by the equation (8) w α = w β cos(α β) sin(α β) w β sin(α β) + cos(α β) Therefore, if w α is known for some value of the parameter α, then it is known for all other values of α The Weyl function w α takes valus in the upper hulf plane (Iλ > 0) when its argument λ lies in the upper half plane It is represented by the formula, [KK], with w α (λ) = bλ + a + + b 0, a R, [ 1 t λ + t ] t 2 dσ α (t), + 1 dσ α (t) t < Note that the measure dσ α is the measure corresponding to the self adjoint extension of the operator D specified by the boundary conditions f 1 (y)e iα = f 2 (y)e iα, α [0, π) Thus we defined the direct spectral transform ψ(x), x [y, + ) w α (λ), Iλ 0 According to the theorem of Marchenko, [MA], this map is injective The inverse map which recovers the potential ψ(x) from the spectral data w α (λ) is called an inverse spectral transform For example, if ψ(x) 0 then ( e iα iλ 2 ϕ α (x, y, λ) = (x y) e iα+ iλ 2 (x y) ) ( ) e iα iλ 2, θ α (x, y, λ) = (x y) ie iα+ iλ 2 (x y) The Weyl function w α (λ, y) = i if Iλ > 0 and i for Iλ < 0 The only potential which have such Weyl function vanishes identically

5 THE ATIYAH-HITCHIN BRACKET 5 Now we are going to change subject Consider m N, a N monopole solution of the Bogomolny equation, see [AH] To any such solution one can associate the scattering function S m (λ), which is the rational function of degree N and such that S m ( ) = 0 Thus we have an analog of the direct spectral transform m N S m (λ) The theorem of S Donaldson [AH] implies that this map is injective The function S m (λ) has the form S m (λ) = N 1 a i λ i i=1 q(λ) = λ N + N 1 p(λ) b j λ j j=0 The space of all such functions we denote by Rat N The monic polynomial p(λ) is determined by its roots λ 0,, λ N 1 The polynomial q(λ) of degree N 1 can be determined from its values at the roots of denominator Therefore, λ 0,, λ N 1, q(λ 0 ),, q(λ N 1 ) are global coordinates on Rat N Atiyah and Hitchin defined ω, the 2N form on Rat N by the formula ω = N 1 i=0 The corresponding Poisson bracket is δq(λ i ) q(λ i ) δλ i (9) {q(λ k ), λ n } = q(λ k ) δ k n, {λ k, λ n } = {q(λ k ), q(λ n )} = 0 In the recent paper Faybusovich and Gehtman, [FG], wrote the Atiyah-Hitchin structure on rational functions in compact coordinate free form (10) {w(λ), w(µ)} = (w(λ) w(µ))2 λ µ The meaning of this formula is the following Fix some λ and µ and consider w(λ) and w(µ) as a functions of coordinates The bracket (9) between these two functions can be computed in terms of w(λ) and w(µ) itself In [V1] for rational functions we gave a direct proof that (10) implies (9) In fact, the formula (10) itself can be a starting point for construction of the

6 6 KL VANINSKY bracket As we will see the Atiyah Hitchin bracket can be extended to much wider class then rational functions Here we establish some of its remarkable properties We think about w(λ) as an element of some commutative complex algebra which depends holomorphicaly on the parameter λ Evidently, (6) is skew symmetric with respect to λ and µ It is natural to require linearity of the bracket (11) {aw(λ) + bw(λ), w(ν)} = a{w(λ), w(ν)} + b{w(λ), w(ν)}, where a and b are constants The symbol w(λ) for λ inside the contour C is given by the Cauchy formula w(λ) = 1 2πi C w(ζ) ζ λ dζ Whence due to (11) the values of the bracket in different points are related {w(λ), w(µ)} = 1 2πi C {w(ζ), w(µ)} dζ ζ λ It can be verified that the AH bracket satisfies this compatibility condition Also it is naturally to require the Leibnitz rule to hold (12) {w(λ)w(µ), w(ν)} = w(λ){w(µ), w(ν)} + w(µ){w(λ), w(ν)} It can be verified in a long but simple calculation that (11) and (12) imply the Jacobi identity {w(λ), {w(µ), w(ν)}} + {w(µ), {w(ν), w(λ)}} + {w(ν), {w(λ), w(µ)}} = 0 The particular useful to us Lemma 1 The bracket (6) is invariant under linear fractional transformations (13) w w = aw + b cw + d, where a, b, c, d are constants Proof Consider transformation of the form (14) w w = 1 cw + d

7 THE ATIYAH-HITCHIN BRACKET 7 Then {w (λ), w (µ)} = = 1 1 (cw(λ) + d) 2 {cw(λ) + d, cw(µ) + d} (cw(µ) + d) ((cw(λ) + d) (cw(µ) + d)) 2 (cw(λ) + d) 2 (cw(µ) + d) 2 λ µ = (w (λ) w (µ)) 2 λ µ To finish the proof we note that two consecutive transformations of the form (14) produce the whole group (3) We are done Now we can formulate our main result Theorem The Poisson bracket (4) on the phase space corresponds to the Atiyah-Hitchin bracket (6) on the Weyl functions w α (y, λ) We would like to remark that if the theorem is proved for some value of α then the formula (4) will produce the bracket for all other w α There is no preference in the choice of the boundary condition Therefore, the AH bracket must be invariant under the one parameter subgroup of rotations of the upper half plane defined by the formula (8) In fact, it is invariant under the whole group of linear fractional transformations 3 The Weyl functions of the Dirac operator Proof First we will obtain the expression for the function w α (λ, y) in terms of any square integrable solution H(x, λ) This solution is proportional to h α (x, y, λ) with some constant c Thus ϕ α + w α θ α = ch and normalization conditions (7) imply the system e iα + w α ie iα = ch 1 (y), e iα w α ie iα = ch 2 (y) Solving for w α (15) w α (λ, y) = H 1(y)e iα H 2 (y)e iα H 1 (y)ie iα + H 2 (y)ie iα This identity can be written as w α = (H 1e iβ H 2 e iβ ) cos(α β) (H 1 ie iβ + H 2 ie iβ ) sin(α β) (H 1 ie iβ + H 2 ie iβ ) cos(α β) + (H 1 e iβ H 2 e iβ ) sin(α β) Dividing the numerator and denominator by H 1 ie iβ + H 2 ie iβ we obtain (8) The standard monodromy matrix M(x, y, λ) = ( M (1), M (2)) is 2 2 matrix solution of the auxiliary linear problem M = V 2 (xλ)m which satisfies the

8 8 KL VANINSKY boundary condition M(x, y, λ) x=y = I In [V4] we considered the Weyl solution ĥ T (x, y, λ) = (ĥ1, ĥ2) defined by the formula ĥ = M(1) ŵ + M (2) with the new Weyl function ŵ(y, λ) Arguing as before (16) ŵ = i w α i + w α e i2α This implies, in particular, that the function ŵ(λ) for Iλ > 0 takes values in the unit circle The transformation z = i w i + w establishes 1:1 correspondence between the upper half plane and the enterior of the unit circle The shift α β described on the w plane by the formula (8) corresponds to the rotation on the angle 2(α β) of the z plane The Weyl function w α (λ, y) and ŵ(λ, y) for fixed λ and α are the functions of y and the potential ψ(x), x [y, + ] Formula (15) and the spectral problem imply the Ricatti type equation for w α (λ, y) dw α (λ, y) dy = λ 2 (w2 α (λ, y)+1) w α(λ, y)(ψ α +ψ α )+ i 2 (w2 a (λ, y)+1)(ψ α ψ α ), where ψ α = ψe i2α Similar, for the function ŵ(y, λ) dŵ(y, λ) dy = iλŵ(y, λ) + ψ ψŵ 2 (y, λ) This is the general situation: the function ŵ is easier to deal with then w α We compute formal variational derivatives for ŵ only Lemma 2 The following identities hold, δŵ(y, λ) δψ(x) δŵ(y, λ) δψ(x) = ] 2 [ĥ1 (x, y, λ), = [ĥ2 (x, y, λ)] 2 for y x The derivatives vanish for y x For the first time the lemma was proved in [V4] Here we present simplified proof of the first formula The argument for the second formula is the same We split the proof into 3 steps

9 THE ATIYAH-HITCHIN BRACKET 9 Step 1 Let M = δm be a variation of M(x, 0) in response to the variation of ψ(y) and ψ(y), 0 y x Then M = V M + V M The solution of this nonhomogenious equation is Therefore, x M (x) = M(x) δm(x, 0) δψ(y) δm(x, 0) δψ(y) 0 M 1 (ξ)v (ξ)m(ξ) dξ ( 0 0 = M(x, y) 1 0 ( 0 1 = M(x, y) 0 0 ) M(y, 0), ) M(y, 0) Step 2 The purpose of this step is to prove the formula [ δŵ(0, λ) m12 (y, 0) = ŵ(0, λ, )ŵ(y, λ) δψ(y) A where A = m 12 (y, 0) ŵ(y, λ)m 22 (y, 0), B = m 11 (y, 0) + ŵ(y, λ)m 21 (y, 0) + m 11(y, 0) B Consider the eigenvalue problem: Df = λ f on the finite interval [y, b] The 2 solution f(x, λ) = M (1) (x, y, λ)ŵ(b, y, λ) + M (2) (x, y, λ) with some ŵ(b, y, λ) satisfies the boundary condition f 1 (b, λ) = f 2 (b, λ) if The limit ŵ(b, y, λ) = m 22 m 12 m 11 m 21 (b, y, λ) ŵ(y, λ) = lim y, λ) b + ŵ(b, exists because the spectral problem is in the limit point case Therefore, To compute the derivative ŵ(y, λ) = lim ŵ(b, y, λ) b + ŵ(b, 0, λ) = m 22 m 12 m 11 m 21 m 22 m 12 m 22 m 12 (b, 0, λ) m 22 m 12 m 11 m 21 m 11 m 21 m 11 m 21 (b, 0, λ) ],

10 10 KL VANINSKY we use the formulae of Step 1 δŵ(b, 0, λ) δψ(y) = ŵ(b, 0, λ) [m 22 m 12 ](b, y)m 12 (y, 0) [m 22 m 12 ](b, 0) ŵ(b, 0, λ) [m 12 m 22 ](b, y)m 11 (y, 0) [m 11 m 21 ](b, 0) = ŵ(b, 0, λ)[m 22 m 12 ](b, y) [ m 12 (y, 0) [m 22 m 12 ](b, 0) + m 11 (y, 0) [m 11 m 21 ](b, 0) Using the identity M(b, 0) = M(b, y)m(y, 0), and simple algebra, one finds ŵ(b, 0, λ)ŵ(b, y, λ) [ m 12 (y, 0) m 12 (y, 0) ŵ(b, y, λ)m 22 (y, 0) + m 11 (y, 0) m 11 (y, 0) + ŵ(b, y, λ)m 21 (y, 0) Now pass to the limit when b and Step 3 Consider ĥ(x, y, λ) and proportional to it ĥ(x, 0, λ), then [ŵ(y, ] ĥ(y, y, λ) = M (1) (y, y)ŵ(y, λ) + M (2) λ) (y, y) =, 1 (17) ŵ(y, λ) = ĥ1(y, y, λ) ĥ 2 (y, y, λ) = ĥ1(y, 0, λ) ĥ 2 (y, 0, λ) Therefore, ŵ(y, λ)[m 21 (y, 0)ŵ(0, λ) + m 22 (y, 0)] = m 11 (y, 0)ŵ(0, λ) + m 12 (y, 0) After simple algebra, ŵ(0, λ) (18) = 1 A B Therefore, using (17) and (18) we obtain [ŵ(0, δŵ(0, λ) λ)m12 (y, 0) = ŵ(y, λ) + ŵ(0, λ)m ] 11(y, 0) δψ(y) A B [ŵ(0, ] λ)m11 (y, 0) + m 12 (y, 0) = ŵ(y, λ) B = ĥ 1 (y, λ) m 11 (y, 0)ĥ2(y, 0) ĥ 1 (y, 0) m 21(y, 0) [ĥ1 (y, λ)] 2 = m 11 (y, 0)ĥ2(y, 0) m 21 (y, 0)ĥ1(y, 0) ] ]

11 THE ATIYAH-HITCHIN BRACKET 11 The denominator does not depend on y and can be computed for y = 0, where it is equal to 1 We are done The next identity for the quartic products of the solution is standard It can be verified by direct calculation Lemma 3 Let f +, f are two solution of (5) corresponding to the spectral parameter λ, and g +, g two solutions corresponding to the spectral parameter µ Then, f 1 + f 1 g+ 2 g 2 f 2 f+ 2 g+ 1 g 1 = 1 i(µ λ) d [( )( )] j T dx + Jg + j T Jg Now we can easily prove the Theorem Due to translation invariance we can assume y = 0 Then for ŵ(λ) = ŵ(0, λ) by the Lemma 2 {ŵ(λ), ŵ(µ)} = 2i = 2i Using the identity of Lemma 3 = 0 δŵ(λ) δψ(x) δŵ(µ) δψ(x) δŵ(λ) δŵ(µ) δψ(x) δψ(x) dx ĥ2 2 (x, λ)ĥ2 1 (x, µ) + ĥ2 1 (x, λ)ĥ2 2 (x, µ) dx 2i ] 2 [ĥ(λ) T i(µ λ) Jĥ(µ) + (ŵ(λ) ŵ(µ))2 0 = 2 µ λ Thus the result is proved for the Weyl functions ŵ The formula (16) and invariance of the AH bracket under linear fractional transformations imply that the same formula holds for the Weyl function w α (λ, y) Theorem is proved References [AH] [FT] [FG] [KK] [KP] [LS] M Atiyah and N Hitchin, Geometry and dynamics of magnetic monopoles, Princeton University press, Princeton, NJ, 1988 LD Faddeev and LA Takhtadzian, Hamiltonian Methods in the Theory of Solitons, Springer Verlag, Berlin, 1988 L Faybusovich and M Gehtman, Poisson brackets on rational functions and multi Hamiltonian structures for integrable lattices, Phys Lett A 272 (2000), IS Kac and MG Krein, r-functions analytic functions mapping the upper half-plane into itself, Transl Amer Math Soc 103 (1974), 1-18 I Krichever and D Phong, Symplectic forms in the theory of solitons, Surveys in differential geometry IV (1998), Int Press, Boston, BM Levitan and Sargsian IS, Sturm-Liouwille and Dirac operators, Kluwer Academic, Dordrecht, 1991

12 12 KL VANINSKY [MA] [V1] [V2] [V3] [V4] [W] VA Marchenko, Sturm-Liouville operator and applications, Birkauser Verlag, Basel, 1986 K Vaninsky, The Atiyah Hitchin bracket and open Toda Lattice, to appear in, Journal of Geometry and Physics K Vaninsky, The Camassa Holm equation on the line and Jacobi elliptic coordinates, submitted to CMP K Vaninsky, Symplectic structure for the nonlinear Schrödinger operator in the periodic and scattering case, submitted to Mathematical Physics, Analysis, Geometry K Vaninsky, A convexity theorem in the scattering theory for the Dirac operator, Trans Amer Math Soc 350 (1998), H Weyl, Uber gewohnliche Differentialgleichungen mit Singulariten und die zugehorigen Entwicklungen willkurlichen Funktionen, Math Ann 68 (1910),