arxiv:nlin/ v1 [nlin.si] 7 Jul 2005

Transcription

1 Completeness of the cubic and quartic Hénon-Heiles Hamiltonians Robert CONTE Micheline MUSETTE and Caroline VERHOEVEN Service de physique de l état condensé URA 464 CEA Saclay F 99 Gif-sur-Yvette Cedex France Conte@drecam.saclay.cea.fr arxiv:nlin/05070v [nlin.si] 7 Jul 005 Dienst Theoretische Natuurkunde Vrije Universiteit Brussel and International Solvay Institutes for Physics and Chemistry Pleinlaan B 050 Brussels Belgium MMusette@vub.ac.be CVerhoev@vub.ac.be November 004 Abstract The quartic Hénon-Heiles Hamiltonian H = P +P /+Ω Q +Ω Q /+CQ 4 +BQ Q +AQ 4 + /α/q +β/q γq passes the Painlevé test for only four sets of values of the constants. Only one of these identical to the traveling wave reduction of the Manakov system has been explicitly integrated Wojciechowski 985 while the three others are not yet integrated in the generic case αβγ 000. We integrate them by building a birational transformation to two fourth order first degree equations in the classification Cosgrove 000 of such polynomial equations which possess the Painlevé property. This transformation involves the stationary reduction of various partial differential equations PDEs. The result is the same as for the three cubic Hénon-Heiles Hamiltonians namely in all four quartic cases a general solution which is meromorphic and hyperelliptic with genus two. As a consequence no additional autonomous term can be added to either the cubic or the quartic Hamiltonians without destroying the Painlevé integrability completeness property. Keywords: Hénon-Heiles Hamiltonian Painlevé property hyperelliptic. separation of variables PACS 995: 0.30.Hq Introduction The considered Hamiltonian originates from celestial mechanics as a system describing the motion of a star in the axisymmetric potential of the galaxy. Denoting q the radius and q the altitude this Hénon- Heiles Hamiltonian HH [8] is the sum of a kinetic energy and a potential energy in which the potential is a cubic polynomial in the position variables q q H = p +p +q +q +q q 3 q3 it is nonintegrable and displays a strange attractor. However if one changes the numerical coefficients in the potential the system may become integrable and this question to find all the integrable cases and to integrate them has attracted a lot of activity in the last three decennia. A prerequisite is to define the word integrability and in section we briefly recall its three main acceptations in the context of Hamiltonian systems. In section 3 we recall all the cases three cubic plus four quartic for which the most general two-degree of freedom classical time-independent Hamiltonian may have a single valued general solution. Corresponding author RC. Preprint S004/047. nlin.si/05070

2 Then discarding the integrated cases see [0] for a review of the current state of this problem we focus on the three cases all quartic for which the general solution is still missing with the aim of finding this general solution. In section 4 we build an equivalent fourth order ordinary differential equation ODE for q t in the hope of finding it listed in one of the classical tables of explicitly integrated ODEs. This hope is deceived because these tables are not yet finished. This is why in the lasttwo sections we adopt adifferent strategy. In front ofthe difficulty to perform the separation of variables in the sense of Arnol d and Liouville we establish a birational transformation between the two second order Hamilton equations and a fourth order ODE listed in a classical table established by Cosgrove [] whose general solution is single valued. Integrability for Hamiltonian systems Given a Hamiltonian system with a finite number N of degrees of freedom three main definitions of integrability are known. the one in the sense of Liouville that is the existence of N independent invariants K j whose pairwise Poisson brackets vanish {K j K l } = 0. the one in the sense of Arnol d-liouville[ chap. 9] which is to find explicitly some canonical variables s j r j j = N which separate the Hamilton-Jacobi equation for the action S which for two degrees of freedom writes as Hq q p p E = 0 p = S q p = S q 3. the one in the sense of Painlevé [8] i.e. the representation of the general solution q j t by an explicit closed form single valued expression of the time t. 3 The seven Hénon-Heiles Hamiltonians Given the most general two-degree of freedom classical time-independent Hamiltonian H = p +p +Vq q = E 3 the requirement that the system made of the two Hamilton equations passes the Painlevé test [8] for at least some integer powers q n qn selects seven and only seven potentials V depending on a finite number of constants namely. three cubic potentials HH3 case [7 5 9] H = p +p +ω q +ω q +αq q 3 βq3 + γq α 0 4 in which the constants αβω ω γ can only take three sets of values SK : β/α = ω = ω 5 KdV5 : β/α = 6 6 KK : β/α = 6ω = 6ω. 7

3 . four quartic potentials HH4 case [4 7] H = P +P +Ω Q +Ω Q +CQ 4 +BQ Q +AQ 4 + α Q + β Q γq B 0 8 in which the constants ABCαβγΩ Ω can only take the four values the notation A : B : C = p : q : r stands for A/p = B/q = C/r = arbitrary A : B : C = : : γ = 0 A : B : C = : 6 : γ = 0 Ω = Ω A : B : C = : 6 : 8 α = 0 Ω = 4Ω A : B : C = : : 6 γ = 0 Ω = 4Ω. All seven cases are integrable in the sense of Liouville with a second constant of the motion K [ 4 9] [0 3 4] either quadratic or quartic in the momenta p p. In the sense of Arnol d-liouville the separation of variables has been performed [ ] except in three cases. HH4 :6: α β. HH4 :6:8 βγ 0 3. HH4 ::6 αβ 0. What is remarkable is the fact that in all cases when the separation of variables is achieved the equations of Hamilton have the Painlevé property the general solution being a hyperelliptic function of genus two. The purpose of this work is to prove equally the Painlevé property in the three remaining cases where the separation of variables is not yet performed. 4 Equivalent fourth order ODEs In the cubic case the two Hamilton equations q q βq +αq = 0 0 q q +αq q γq 3 = 0 together with the Hamiltonian 4 are equivalent [5] to a single fourth order ODE for q t q +8α βq q α+βq 0 3 αβq3 +ω +4ω q +6αω 4βω q +4ω ω q +4αE = 0 independent of the coefficient γ of the nonpolynomial term q and depending on the constant value E of the Hamiltonian H. In the three HH3 cases 5 7 this ODE belongs to a list [] classification of equations enjoying the Painlevé property whose general solution is hyperelliptic with genus two. In the quartic case the similar fourth order equation is built by eliminating Q and Q between the two Hamilton equations Q +Ω Q +4CQ 3 +BQ Q αq 3 +γ = 0 3 Q +Ω Q +4AQ 3 +BQ Q βq 3 =

4 and the Hamiltonian 8 which results in Q + Q Q + +6 A Q B +8 6 AC B C B [ + A B ω 4ω + Q Q Q Q Q Q Q +4B CQ Q +4C 4 AC γ α B B + A 4 +3 A ] B Q B Q 4 Q +6 A α B Q 7 +0 α Q 5 Q A γα B Q AB ω ω γ α Q 3 A +6 B γ +Bα 8 AC B α + 6 A Q B ω 4ω ω 8BE Q +48 AC B γq +4 AC B 4C B Q 5 γ Q Q ω Q 3. 5 This ODE depends on E but not on β and as opposed to the cubic case it does not belong to a classified set of equations because Q is not polynomial in Q. In the three remaining cases since one is yet unable either to perform the separation of variables or to establish a direct link to a classified ODE let us build an indirect link to such a classified ODE. This link which involves soliton equations is the following. For each of the seven cases the two Hamilton equations are equivalent [5 6 3] to the traveling wave reduction of a soliton system made either of a single PDE HH3 or of two coupled PDEs HH4 most of them appearing in lists established from group theory [3]. Among the various soliton equations which are equivalent to them via a Bäcklund transformation some of them admit a traveling wave reduction to a classified ODE.This propertydefines apath [ 30] which starts from oneofthe three remaininghh4 cases goes up to a soliton system of two coupled +-dimensional PDEs admitting a reduction to the considered case then goes to another +-dim PDE system equivalent under a Bäcklund transformation finally goes down by reduction to an already integrated ODE or system of ODEs. 5 General solution of the quartic :6: and :6:8 cases and Let us denote the two constants of the motion of the :6: and :6:8 cases as H = P +P + Ω Q +Q 3 Q4 +6Q Q +Q4 κ Q + κ Q = E : 6 : K = P P +Q Q Q +Q +Ω 8 κ κ Q + κ 4 Q +κ Q κ + κ Q Q P Q P H = p +p + ω 4q +q 6 8q4 +6qq +q 4 γq + β = E q : 6 : 8 K = p q 6 q +4q +ω+ β q +γ 4 q q p q p γq 4q p p + q q 4 +q 3 q +4q p 4ωq q +4q β q

5 There is a canonical transformation [3] between the :6: and :6:8 cases mapping the constants as follows E :6:8 = E :6: K :6:8 = K :6: ω = Ω γ = κ +κ β = κ κ 8 therefore one only needs to integrate either case. The path to an integrated ODE comprises the following three segments. The coordinate q t of the :6:8 case can be identified [4 3] to the component F of the traveling wave reduction fxτ = Fx cτgxτ = Gx cτ of a soliton system of two coupled KdV-like equations c-kdv system denoted c-kdv [4 3] f τ + f xx + 3 ff x f3 +3fg = 0 x g τ +g xxx +6gg x +3fg xx +6gf xx +9f x g x 3f g x f xxxx + 3 ff xxx+9f x f xx 3f f xx 3ffx = 0 with the identification { q = F q = F +F +G ω c = ω K = γ K = E 0 in which K and K are two constants of integration. There exists a Bäcklund transformation between this soliton system and another one of the c-kdv type denoted bi-sh system [3] { uτ + u xx +u +6v x = 0 v τ +v xxx +uv x = 0. This BT is defined by the Miura transformation [] u = 3 g fx f v = 3 fxxx +4ff xx +8gf x +4fg x +3fx 4 f f x f 4 +4gf. Finally the traveling wave reduction uxτ = Ux cτvxτ = Vx cτ can be identified [30] to the autonomous F-VI equation a-f-vi in the classification of Cosgrove [] a-f-vi : y = 8yy +9y 4y 3 +α VI y + α VI 9 y +κ VIt+β VI κ VI = 0 3 an ODE whose general solution is meromorphic expressed with genus two hyperelliptic functions [ Eq. 7.6]. The identification is U = 6 y + c 8 V = y 6y cy c K A 4 α VI = 4c β VI = K B ck A c3 in which K A K B are two constants of integration. In order to perform the integration of both the :6: and the :6:8 cases it is sufficient to express FG 5

6 rationally in terms of UVU V. The result is F = W W + K [ 3U U 3c V +U +3c 4W +36K B 54K] G = U 3 + [ V +3K V +K U 3K 8W U 3c K V +K U +3c] W = V + 3 K + 3 K U 3c K A = K. 5 Making the product of the successive transformations0 5 4 one obtains a meromorphic general solution for Q Q q q [ 9j 3 y + 49 h+e ω 94 ] γ q = W W + γ W q = 6 y 5 9 ω + W [ y + γ 48y 3 6ωy + 4E ωe β 4γ W = h+e 9γ y 5 9 ω 9 ω y ω3 y 59 ω h 44γ y 5 9 ω ] α VI = 4ω β VI = 3 4 γ +ωe 3 6 β 5 43 ω3 K VI = 3 3 K E K VI = 3 3 EK 3 E βγ K = γ K = E K A = E K B = 3 6 β γ. in which h and j are convenient auxiliary variables [ Eqs ] y = Qs s + Qs Qs s C VI s C 36 α VI VI h = 3 4 E s s +C VI + s C VIs C VI VI F VI s +s j = { 6 h+f VI y + α } VI 9 E VI 4s +s. 6 7 In the above the variables s s are defined by the hyperelliptic system [] s s s = Ps s s s = Ps Ps = s C VI Qs Qst = s C VI t C VI α VI s +t C VI + E VI s+t+f VI Qs = Qss. 8 Theexpressions7cannotbewrittenasrationalfunctionsofs s s s andareneverthelessmeromorphic [4 ]. The coefficients αcef VI of the hyperelliptic curve depend algebraically on the parameters of the 6

7 Hamiltonians β γ ω E K [ Eqs ] A VI = 4ω EVI = 6 3 ωf VI E β +4γ C VI EVI = 4 3 F VI 4E +K F VI E F VI +4E+ 9K 4 F VI E 7 4 βγ = 0 and this algebraic dependence could explain the difficulty to separate the variables in the Hamilton-Jacobi equation. Note that in the particular case βγ = 0 i.e. κ = κ these coefficients become rational see [9]. Remark. The F-VI ODE can be written in Hamiltonian form H = P +Q P Q Q Q 3 96 α3 VIQ + 6 α VIQ 8 α VIQ 6β VI Q 6κ VI tq α4 VI + 9 α VIβ VI 30 Q = 6 y α VI Q = 6y P = 6y 08yy P = 6y. 7 In the autonomous case κ VI = 0 the Hamiltonian H is a first integral equal to 36K VI and the other constant of the motion is cubic in the momenta. However because of the nonlinear link between K VI and the two first integrals of the :6:8 case see 4 there exists no canonical transformation between the variables q j p j of :6:8 and the above canonical variables of a-f-vi. 6 General solution of the quartic ::6 case 9 Let us denote the two constants in involution as H = P +P + Ω 8 4Q +Q 3 6Q4 +Q Q +Q4 κ Q + 4κ Q = E : : 6 K = 8Q P Q P P Q Q 4 Q 3 6 Q +ΩQ Q κ +3Q Q +κ Q 4 4 Q P Q. 3 Similarly to the :6:-:6:8 couple there exists a canonical transformation between the ::6 Hamiltonian and another Hamiltonian [3 4] which is however not the sum of a kinetic energy and a potential energy which we denote similarly as 5:9:4 H = p + p 3 q q 8 4q4 +9qq +5q+ 4 ω q +q κq + ζ = E 5 : 9 : 4 K = q q p +q q q q p q 4 4κq q p +q q +p 4q 4 q p q +ωq +4κ +4ζ +4κ q p + q q +4q +q κ +4q q q p 3 q q + q +q q p q 4 ζ +q q 4 +8κq 4κ p q qq +q 4 +κq E 5:9:4 = E ::6 K 5:9:4 = K ::6 ω = Ω κ = κ +κ ζ = κ κ. 7 q 3

8 ThepathtoanintegratedODEisalsoquitesimilarandismadeofthefollowingthreesegments[37]. Firstly the coordinate q t of 3 is identified [4 3] to the component F of the traveling wave reduction fxτ = Fx cτgxτ = Gx cτ of a soliton system of two coupled KdV-like equations denoted c-kdvaf g [4 3] { q = F q = F F G+ω 5 33 c = ω. Secondly there exists a Bäcklund transformation between this soliton system and another one of the c-kdv type denoted bi-sk system u v [3] transformation defined by the Miura map u = 3 3fx f +g 0 v = 9 fxxx +g xx +f x g fg x ff xx +g Finally the traveling wave reduction uxτ = Ux cτvxτ = Vx cτ is identified [30] U = 3 y ω V = 6y +8y ωy + 0 ω 3 E 35 5 to the F-IV equation or to the F-III as well in the classification of Cosgrove [] y = 30yy 60y 3 +α IV y +β IV y = s +s +s +s s +s +A s s s F-IV = Ps s s s = Ps Ps = s +A 3 α IV 3 s +A+Bs+ β IV 3 3B K IV = K IV = 9AB in which K IV K IV denote two polynomial first integrals of F-IV. The general solution of this ODE is meromorphic expressed with genus two hyperelliptic functions []. In order to perform the integration of both Hamiltonians 3 and 3 it is sufficient to express FG rationally in terms of UVU V. The result is F = W W +K ax G = F X X +K a 54U W = X +08K a X 54K a U + 3ω 0 X = V +U 3ωU ω 7 X = 9 5 E 4U UV 8 5 ωu + 48 ωv EU 4 5 ω U 9 κ +κ 9 K a K a = κ κ. U + 3ω W + U + 9ω 0 WX ωe 5 ω3 From the point of view of the separation of variables one should first exhibit a Hamiltonian representation of F-IV. One such structure is that of the cubic SK case. However since the constant value of the Hamiltonian of the cubic SK case when expressed only in terms of the parameters EKωκ κ of the ::6 is not an affine function of E there exists no canonical transformation between the cubic SK case and the ::6 case. 37 8

9 7 Conclusion remaining work The explicit integration of all the seven cases is now achieved in the Painlevé sense finding a closed form single valued expression for the general solution and the common features are the following.. In all cases the general solution is hyperelliptic with genus two and therefore meromorphic.. Each case is birationally equivalent to a fourth order ODE which is complete in the Painlevé sense i.e. which accepts no additional term under penalty of losing its Painlevé property. Consequently for each of the seven Hamiltonians it is impossible to add any term to the Hamiltonian without destroying the Painlevé property and the seven Hénon-Heiles Hamiltonians are complete. About the integration in the Arnol d-liouville sense finding the separating variables of the Hamilton- Jacobi equation two problems remain open.. In the :6:-:6:8 case the hyperelliptic curve y = Ps of F-VI see 8 reduces in the separated cases βγ = 0 to the hyperelliptic curve of the separating variables. Therefore F-VI is the good ODE to consider and the only missing item is to find a Hamiltonian structure of F-VI necessarily distinct from 30 admitting a canonical transformation to :6:-:6:8.. In the ::6-5:9:4 case the hyperelliptic curve y = Ps of F-IV see 36 does not reduce in the separated cases κ κ = 0 [7] to the hyperelliptic curve of the separating variables which is κ κ = 0 : Ps = s 6 ωs 3 +Es + K 0 s+κ +κ = Therefore F-IV as well as its birationally equivalent ODE F-III is not the good ODE to consider and it should be quite instructive to integrate the fourth order equivalent ODE 5 in that case. Acknowledgments The authors acknowledge the financial support of the Tournesol grant no. T CV is a postdoctoral fellow at the FWO-Vlaanderen. References [] S. Abenda and Yu. Fedorov Acta Appl. Math [] V.I. Arnol d Les méthodes mathématiques de la mécanique classique Nauka Moscou 974 Mir Moscou 976. [3] S. Baker PhD Thesis University of Leeds 995. [4] S. Baker V. Z. Enol skii and A. P. Fordy Phys. Lett. A [5] M. B laszak and S. Rauch-Wojciechowski J. Math. Phys [6] T. Bountis H. Segur and F. Vivaldi Phys. Rev. A [7] Chang Y. F. M. Tabor and J. Weiss J. Math. Phys [8] R. Conte The Painlevé property one century later ed. R. Conte CRM series in mathematical physics Springer New York 999. solv-int/ [9] R. Conte A. P. Fordy and A. Pickering Physica D [0] R. Conte M. Musette and C. Verhoeven J. Nonlinear Math. Phys. Supp nlin.si/

10 [] C. M. Cosgrove Stud. Appl. Math [] J. Drach C. R. Acad. Sc. Paris [3] V. G. Drinfel d and V. V. Sokolov Soviet Math. Dokl [4] N. Farkas and E. Kra Riemann surfaces Springer Berlin 980. [5] A. P. Fordy Physica D [6] A.P. Fordy and P. P. Kulish Commun. Math. Phys [7] B. Grammaticos B. Dorizzi and A. Ramani J. Math. Phys [8] M. Hénon and C. Heiles Astron. J [9] J. Hietarinta J. Math. Phys [0] J. Hietarinta Phys. Rep [] D. Mumford Tata lectures on theta II Birkhäuser Basel 983. [] M. Musette and C. Verhoeven Theor. Math. Phys [3] A. Ramani in Fourth international conference on collective phenomena ed. J. L. Lebowitz Annals of the New York Academy of Sciences NY Ac. Sc. New York 98. [4] A. Ramani B. Dorizzi and B. Grammaticos Phys. Rev. Lett [5] V. Ravoson L. Gavrilov and R. Caboz J. Math. Phys [6] V. Ravoson A. Ramani and B. Grammaticos Phys. Letters A [7] C. Verhoeven PhD thesis Vrije Universiteit Brussel May 003. [8] C. Verhoeven M. Musette and R. Conte J. Math. Phys nlin.si/0030. [9] C. Verhoeven M. Musette and R. Conte Theor. Math. Phys nlin.si/0300. [30] C. Verhoeven M. Musette and R. Conte pages Bilinear integrable systems - from classical to quantum continuous to discrete ed. P. van Moerbeke Kluwer Dordrecht 004. nlin.si/ [3] S. Wojciechowski Physica Scripta