arxiv: v2 [nlin.si] 27 Jun 2015

Similar documents
Canonicity of Bäcklund transformation: r-matrix approach. I. arxiv:solv-int/ v1 25 Mar 1999

New conditional integrable cases of motion of a rigid body with Kovalevskaya's configuration Author(s): Yehia, HM (Yehia, H. M.)[ 1 ] Elmandouh, AA

A MULTI-COMPONENT LAX INTEGRABLE HIERARCHY WITH HAMILTONIAN STRUCTURE

Supersymmetric Sawada-Kotera Equation: Bäcklund-Darboux Transformations and. Applications

Equivalence of superintegrable systems in two dimensions

Models of quadratic quantum algebras and their relation to classical superintegrable systems

Nondegenerate 2D complex Euclidean superintegrable systems and algebraic varieties

Lattice geometry of the Hirota equation

Prolongation structure for nonlinear integrable couplings of a KdV soliton hierarchy

arxiv: v1 [math-ph] 15 Sep 2009

2. Examples of Integrable Equations

Linearization of Mirror Systems

1 Introduction The search for constants of motion in the Lagrangian approach has been traditionally related with the existence of one{parameter subgro

Superintegrable 3D systems in a magnetic field and Cartesian separation of variables

COULOMB SYSTEMS WITH CALOGERO INTERACTION

Relativistic Collisions as Yang Baxter maps

New Analytical Solutions For (3+1) Dimensional Kaup-Kupershmidt Equation

Bäcklund transformation and special solutions for Drinfeld Sokolov Satsuma Hirota system of coupled equations

A supersymmetric Sawada-Kotera equation

From Continuous to Discrete Equations via Transformations

Boussineq-Type Equations and Switching Solitons

Variable separation and second order superintegrability

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

Introduction to the Hirota bilinear method

A Short historical review Our goal The hierarchy and Lax... The Hamiltonian... The Dubrovin-type... Algebro-geometric... Home Page.

arxiv: v1 [nlin.si] 10 Feb 2018

arxiv: v1 [math-ph] 7 Apr 2014

Solutions of the nonlocal nonlinear Schrödinger hierarchy via reduction

arxiv: v2 [math-ph] 24 Feb 2016

Warped product of Hamiltonians and extensions of Hamiltonian systems

arxiv:math-ph/ v1 25 Feb 2002

Nondegenerate three-dimensional complex Euclidean superintegrable systems and algebraic varieties

arxiv: v2 [nlin.si] 22 Dec 2016

arxiv:nlin/ v2 [nlin.si] 15 Sep 2004

A connection between HH3 and KdV with one source

Diagonalization of the Coupled-Mode System.

Models for the 3D singular isotropic oscillator quadratic algebra

arxiv: v1 [nlin.si] 23 Aug 2007

Nonlinear Single-Particle Dynamics in High Energy Accelerators

DUAL HAMILTONIAN STRUCTURES IN AN INTEGRABLE HIERARCHY

Math 3191 Applied Linear Algebra

Exact Solution for Camassa-Holm Equations which Desribe Pseudo-Spherical Surface

Generalized Stäckel Transform and Reciprocal Transformations for Finite-Dimensional Integrable Systems

Symmetries and Group Invariant Reductions of Integrable Partial Difference Equations

Superintegrable three body systems in one dimension and generalizations

A new integrable system: The interacting soliton of the BO

Superintegrability in a non-conformally-at space

Poisson Algebras on Elliptic Curves

Integrable Curves and Surfaces

The Higher Dimensional Bateman Equation and Painlevé Analysis of Nonintegrable Wave Equations

On recursion operators for elliptic models

Math 240 Calculus III

Second Order Lax Pairs of Nonlinear Partial Differential Equations with Schwarzian Forms

From the Kadomtsev-Petviashvili equation halfway to Ward s chiral model 1

arxiv: v1 [math-ph] 31 Jan 2015

arxiv:nlin/ v2 [nlin.si] 14 Sep 2001

Exact Solutions for the Nonlinear (2+1)-Dimensional Davey-Stewartson Equation Using the Generalized ( G. )-Expansion Method

Reductions to Korteweg-de Vries Soliton Hierarchy

Introduction to Integrability

Illustrating Dynamical Symmetries in Classical Mechanics: The Laplace-Runge-Lenz Vector Revisited

QUADRATIC EQUATIONS AND MONODROMY EVOLVING DEFORMATIONS

On universality of critical behaviour in Hamiltonian PDEs

arxiv: v1 [nlin.si] 25 Mar 2009

Building Generalized Lax Integrable KdV and MKdV Equations with Spatiotemporally Varying Coefficients

Poisson Brackets and Lie Operators

Integrable Hamiltonian systems generated by antisymmetric matrices

New Approach of ( Ǵ/G ) Expansion Method. Applications to KdV Equation

Exact solutions through symmetry reductions for a new integrable equation

Physics 106b: Lecture 7 25 January, 2018

Bäcklund transformations for fourth-order Painlevé-type equations

EXACT SOLITARY WAVE AND PERIODIC WAVE SOLUTIONS OF THE KAUP-KUPERSCHMIDT EQUATION

Integrable Discrete Nets in Grassmannians

The (G'/G) - Expansion Method for Finding Traveling Wave Solutions of Some Nonlinear Pdes in Mathematical Physics

arxiv:math-ph/ v1 10 May 2005

Symbolic Computation of Conserved Densities and Symmetries of Nonlinear Evolution and Differential-Difference Equations WILLY HEREMAN

On the numerical integration of variational equations

About Integrable Non-Abelian Laurent ODEs

Geometry and Exact Soliton Solutions of the Integrable su(3)-

An introduction to Birkhoff normal form

Some Elliptic Traveling Wave Solutions to the Novikov-Veselov Equation

Symmetries in Semiclassical Mechanics

A New Integrable Couplings of Classical-Boussinesq Hierarchy with Self-Consistent Sources

Normal form for the non linear Schrödinger equation

MATHEMATICAL STRUCTURES IN CONTINUOUS DYNAMICAL SYSTEMS

GEOMETRY OF DISCRETE INTEGRABILITY. THE CONSISTENCY APPROACH

Holonomy groups. Thomas Leistner. School of Mathematical Sciences Colloquium University of Adelaide, May 7, /15

arxiv: v1 [math-ph] 8 May 2016

Multisolitonic solutions from a Bäcklund transformation for a parametric coupled Korteweg-de Vries system

Representations of Sp(6,R) and SU(3) carried by homogeneous polynomials

Transparent connections

General N-th Order Integrals of Motion

Elementary maths for GMT

Recursion Systems and Recursion Operators for the Soliton Equations Related to Rational Linear Problem with Reductions

arxiv: v2 [math-ph] 13 Feb 2013

What an Effective Criterion of Separability says about the Calogero Type Systems

Lie Theory of Differential Equations and Computer Algebra

THE RICCI BRACKET FLOW FOR SOLVABLE LIE GROUPS

arxiv:solv-int/ v1 2 Dec 1997

Curves in the configuration space Q or in the velocity phase space Ω satisfying the Euler-Lagrange (EL) equations,

The Kepler Problem and the Isotropic Harmonic Oscillator. M. K. Fung

Transcription:

On auto and hetero Bäcklund transformations for the Hénon-Heiles systems A. V. Tsiganov St.Petersburg State University, St.Petersburg, Russia e mail: andrey.tsiganov@gmail.com arxiv:1501.06695v2 [nlin.si] 27 Jun 2015 Abstract We consider a canonical transformation of parabolic coordinates on the plain and suppose that this transformation together with some additional relations may be considered as a counterpart of the auto and hetero Bäcklund transformations associated with the integrable Hénon-Heiles systems. 1 Introduction According to classical definition by Darboux [4], a Bäcklund transformation between the two given PDEs Eu,x,t) = 0 and Ẽv,y,τ) = 0 is a pair of relations F 1,2 u,x,t,v,y,τ) = 0 1.1) and some additional relationsbetween x,t) and y,τ), which allow to get both equations E and Ẽ. The BT is called an auto-bt or a hetero-bt depending whether the two PDEs are the same or not. The hetero-bts describe a correspondence between equations rather than a one-to-one mapping between their solutions [1, 12]. In the modern theory of partial differential equation Bäcklund transformations are seen also as a powerful tool in the discretization of PDEs [2]. A counterpart of the auto Bäcklund transformations for finite dimensional integrable systems can be seen as the canonical transformation u,p u ) v,p v ), {u i,p uj } = {v i,p vj } = δ ij, i,j = 1,...,n, 1.2) preserving the algebraic form of the Hamilton-Jacobi equations H i u, S ) = α i and H i v, S ) = α i u v associated with the Hamiltonians H 1,...,H n [15]. The counterpart of the discretization for finite dimensional systems is also currently accepted: by viewing the new v-variables as the old u-variables, but computed at the next time step; then the Bäcklund transformation1.2) defines an integrable symplectic map or discretization of the continuous model, see discussion in [8, 17]. The counterpart of the hetero Bäcklund transformations for finite dimensional integrable systems has to be a canonical transformation 1.2), which has to relate two different systems of the Hamilton-Jacobi equations H i u, S ) ) = α i and Hi v, S = α i 1.3) u v 1

and has to satisfy some additional conditions. It is necessary to add these conditions to 1.2) and 1.3) in order to get a non-trivial, usable and efficient theory. The question of how to do it remains open. One of the possible additional conditions may be found in the theory of superintegrable systems. For instance, let us consider integrals of motion for the two-dimensional harmonic oscillator H 1 = p 2 x +p 2 y +ax 2 +y 2 ), H 2 = p 2 x p 2 y +ax 2 y 2 ), which yield the Hamilton-Jacobi equations separable in Cartesian coordinates u = x, y) on the plane. Another pair of Hamiltonians for the same harmonic oscillator H 1 = p 2 x +p2 y +ax2 +y 2 ), H2 = xp y yp x is separable in polar coordinates v = r, ϕ) on the plain. Canonical transformation of variables u,p u ) = x,y,p x,p y ) v,p v ) = r,ϕ,p r,p ϕ ) 1.4) defines a correspondence between the two different systems of the Hamilton-Jacobi equations H 1,2 x,y, S x, S ) ) = α 1,2 and H1,2 r,ϕ, S y r, S = α 1,2. ϕ This correspondence may be considered as a hetero-bt defined by the generating function relations between x,y) and r,ϕ) F = p x rcosϕ+p y rsinϕ, x = rcosϕ, y = rsinϕ, and with an additional condition that Hamilton function H 1 = H 1 is simultaneously separable in u and v variables. Canonical transformation 1.4) can be considered as the semi hetero-bt relating different Hamilton-Jacobi equations, which are various faces of the same superintegrable system. We know that theory of such semi hetero-bts is a profound and very useful theory, both in classical and quantum cases [10]. The main aim of this note is to discuss a correspondence between integrable Hénon-Heiles systems proposed in [13]. This correspondence between different Hamiltonians may be considered as a counterpart of the generic hetero-bts relating different but simultaneously separable in v-variables Hamilton-Jacobi equations. 2 The Jacobi method Let us consider some natural Hamilton function on T R n H = p 2 1 + +p 2 n +Vq 1,...,q n ). 2.1) The corresponding Hamilton-Jacobi equation is said to be separable in a set of canonical coordinates u i if it has the additively separated complete integral Su 1,...,u n ;α 1,...,α n ) = n S i u i ;α 1,...,α n ), where S i are found by quadratures as solutions of ordinary differential equations. In order to express initial physical variables q,p) in terms of canonical variables of separation u,p u ) we have to obtain momenta p ui from the second Jacobi equations i=1 p ui = S iu i ;α 1,...,α n ) u i, i = 1,...,n. 2.2) 2

Solving these equations with respect to α i one gets integrals of motion H i = α i as functions on variables of separation u,p u ). According to Jacobi we can use canonical transformation q,p) u,p u ) in order to construct different integrable systems simultaneously separable in the same coordinates, see pp. 198-199 in [6]: The main difficulty in integrating a given differential equation lies in introducing convenient variables, which there is no rule for finding. Therefore, we must travel the reverse path and after finding some notable substitution, look for problems to which it can be successfully applied. For instance, let us consider some natural Hamiltonian on the plane separable in the parabolic coordinates u 1,2 H 1 = p 2 1 +p2 2 +Vq 1,q 2 ) λ 2q 2 q2 1 λ = λ u 1)λ u 2 ). 2.3) λ In this case second Jacobi equations 2.2) may be rewritten in the following form p 2 u i +U i u i ) = H 1 + H 2 u i, i = 1,2, 2.4) where U i u i ) are functions defined by the potential Vq 1,q 2 ) [9]. Adding together the separated relations 2.4) one gets another integrable Hamiltonian H 1 = 1 ) p 2 u 2 1 +U 1 u 1 )+p 2 u 2 +U 2 u 2 ) = H 1 + H 2 1 + 1 ), 2.5) 2 u 1 u 2 which can be considered as an integrable perturbation of H 2.1) because there exists second independent integrals of motion ) 1 H 2 = p 2 u 1 +U 1 u 1 ) p 2 u 2 U 2 u 2 ) = H 2 1 ). u 1 u 2 Of course, in generic case this perturbation has no physical meaning. Auto-BTs of the initial Hamilton-Jacobi equation u i,p ui ) v i,p vi ) preserve an algebraic form of the initial Hamiltonian H 2.1) and change the form of the second Hamiltonians H 1 = H 1 + H 2 1 + 1 ) 1, H2 = H 2 1 ). 2.6) 2 v 1 v 2 v 1 v 2 We can try to pick out a special and maybe unique auto-bt, which gives physical meaning to the second Hamiltonian H, as some of the possible counterparts of the hetero-bts relating different but simultaneously separable Hamilton-Jacobi equations. 2.1 The Hénon-Heiles systems There are three integrable Hénon-Heiles systems on the plane, which can be identified with appropriate finite-dimensional reductions of the integrable fifth order KdV, Kaup-Kupershmidt and Sawada-Kotera equations [5]. An explicit integration for all these cases is discussed in [3]. We can try to get a hetero-bt for the finite-dimensional Hénon-Heiles systems taking the hetero-bt between these integrable PDEs and then applying the Fordy finite-dimensional reduction [5]. We believe that the same information may be directly extracted from the wellknown Lax representation for the Hénon-Heiles system separable in parabolic coordinates. 3

Let us take a Lax matrix for the first Hénon-Heiles system separable in parabolic coordinates p 2 2 + p 1q 1 λ 2q 2 q2 1 2λ λ Lλ) = aλ 2 +2aq 2 λ+aq1 2 +4q2)+ 2 p2 1 Characteristic polynomial of this matrix p 2 2 p 1q 1 2λ ) det Lλ) µ = µ 2 aλ 3 H 1 4 + H 2 λ, a R. 2.7) contains the Hamilton function of the first Hénon-Heiles system associated with a fifth order KdV H 1 = p 2 1 +p 2 2 16aq 2 q 2 1 +2q 2 2) 2.8) and a second integral of motion H 2 = aq 2 1 q2 1 +4q2 2 )+ p 1q 2 p 1 q 1 p 2 ) 2. 2.9) The auto-bts preserve the algebraic form of the Hamiltonians[15]. Since the characteristic polynomial is the generating function of these integrals of motion, their invariance amounts to requiring the existence of a similarity transformation for the Lax matrix ˆL = VLV 1, associated with the given auto-bt. The matrix V needs not to be unique because a dynamical system can have different auto BTs [8, 17]. In contrast with [7] we do not require that the transformed Lax matrix ˆLλ) have the same structure in the spectral parameter λ as the original Lax matrix. This requirement is a property of the particular BTs associated with the special translations on hyperelliptic Jacobians. Let us consider a special, unique similarity transformation associated with matrix ) L 12 0 V = 4 L 11 ˆL 11 λ) ), 2.10) 4L 12 where L ij are entries of the Lax matrix 2.7) and ˆL 11 λ) = p 2 2 + p 1λ 2q 2 ) 2q 1. The Lax matrix ˆLλ) = VLV 1 has the following properties: 1. first off-diagonal element of the Lax matrix ˆL 12 λ) = λ u 1)λ u 2 ) yields initial parabolic coordinates on the plane 2.3); 2. second off-diagonal element ˆL 21 = 4aλ v 1 )λ v 2 ) 2.11) = 4aλ 2 + 8aq2 1q 2 p 2 1)λ q 2 1 +4aq 2 1 +4q 2 2)+ 2p 1p 1 q 2 p 2 q 1 ) q 2 1 has only two commuting and functionally independent zeroes v 1,2 ; 4

3. the conjugated momenta for u and v variables are the values of the diagonal element p ui = ˆL 11 λ = u i ), p vi = ˆL 11 λ = v i ), i = 1,2. In generic case such 2 2Lax matrices ˆLλ) exist only if the genus of hyperelliptic curve defined by equation ) det Lλ) µ = 0 is no more a number of degrees of freedom. The corresponding transformation of the classical r-matrix is discussed in [13]. In u,p u ) and v,p v ) variables entries of the transformed Lax matrix ˆL have the following form: ˆL 11 = λ u 2 u 1 u 2 p u1 + λ u 1 u 2 u 1 p u2 = λ v 2 v 1 v 2 p v1 + λ v 1 v 2 v 1 p v2 ˆL 12 = λ u 1)λ u 2 ) = λ2 +λv 1 +v 2 )+v 2 1 +v 1 v 2 +v 2 2 p v 1 p v2 ) 2 4av 1 v 2 ) 2 p2 v 1 p 2 v 2 v 1 v 2 ) ˆL 21 = 4a λ 2 +λu 1 +u 2 )+u 2 1 +u 1u 2 +u 2 ) p u1 p u2 ) 2 2 u 1 u 2 ) 2 4p2 u 1 p 2 u 2 ) u 1 u 2 = 4aλ v 1 )λ v 2 ) Thus, we have two families of variables of separation for the first Hénon-Heiles system and canonical transformation between them: u 1,2 = v 1 +v 2 2 + p v 1 p v2 ) 2 ± A 2av 1 v 2 ) 2, p u1,2 = p v 1 p v2 ) p v1 p v2 ) 2 ± A ) 2av 1 v 2 ) 3 p v 1 v 1 +3v 2 ) p v2 v 2 +3v 1 ) 2v 1 v 2 ) 2.12) where A = p v1 p v2 ) 4 +2av 1 v 2 ) 2 p v1 p v2 ) p v1 v 1 3v 2 ) p v2 v 2 3v 1 ) ) a 2 3v 2 1 +2v 1v 2 +3v 2 2 )v 1 v 2 ) 4. We can directly prove that Hamiltonains H 1,2 2.8-2.9) have the same algebraic form in u,p u ) and v,p v ) variables using this explicit canonical transformation. Proposition 1 The auto-bt for the first Hénon-Heiles system is a correspondence between two equivalent systems of the Hamilton-Jacobi equations H 1,2 λ, S ) = α 1,2, λ = u,v, λ where variables u,p u ) and v,p v ) are related by canonical transformation 2.12) and Hamiltonians H 1,2 are defined by the following equations Φλ,µ) = µ 2 aλ 3 = H 1 4 H 2 λ, λ = u 1,2,v 1,2, µ = p u1,2,p v1,2. 2.13) 5

This auto Bäcklund transformation changes coordinates on an algebraic invariant manifold defined by H 1,2 without changing the manifold itself [8]. We can convert this special, unique auto-bt to some analogue of the hetero-bt by adding one more relation. Namely, substituting roots of the off-diagonal element ˆL 2,1 2.11) into the definition 2.6) one gets the Hamilton function for the second integrable Hénon-Heiles system associated with the Kaup-Kupershmidt equation up to rescaling p 1 2p 1 and q 1 q 1 / 2. After canonical transformation H 1 = p 2 1 +p 2 2 2aq 2 3q 2 1 +16q 2 2) 2.14) q,p) Q,P), P 1,2 = p v 1 ±p v2, Q 1,2 = v 1 ±v 2 2.15) 2 2 the same Hamiltonian H 1 = P 2 1 +P 2 2 2aQ 2 3Q 2 1 +Q 2 2) 2.16) defines a third integrable Hénon-Heiles system associated with the the Sawada-Kotera equation. According [5] canonical transformation 2.15) is a counterpart of the gauge equivalence of the Sawada-Kotera and Kaup-Kupershmidt equations. So, all the Hénon-Heiles systems on the plane are simultaneously separable in v-variables, and we suppose that this fact allows us to define some natural counterpart of the hetero-bt. Proposition 2 For the Hénon-Heiles systems on the plane 2.8) and 2.14,2.16) an analogue of the hetero-bt is the correspondence between two different systems of the Hamilton-Jacobi equations H 1,2 u, S ) ) = α 1,2 and H1,2 v, S = α 1,2, u v where variables u,p u ) and v,p v ) are related by canonical transformation 2.12) and Hamiltonians are defined by the following equations and Φλ,µ) = µ 2 aλ 3 = H 1 4 H 2 λ, λ = u 1,2,v 1,2, µ = p u1,2,p v1,2, H 1,2 = Φv 1,p v1 )±Φv 2,p v2 ). This analogue of the hetero-bt relates different algebraic invariant manifolds associated with Hamiltonians H 1,2 and H 1,2 similar to the well-studied relations between different invariant manifolds in the theory of superintegrable systems. For the first Hénon-Heiles system on the plane 2.8) we can consider parabolic variables u 1,2,p u1,2 ) as coordinates on the Jacobian variety defined by equations 2.13). In order to get integrals of motion for the second or third Hénon-Heiles systems 2.14,2.16) we have to take linear combinations of these equations and to make simultaneously the special shift of the coordinates u,p u ) v,p v ) on the Jacobian variety. 3 Integrable Hamiltonian with velocity dependent potential It is well-known that Hamilton-Jacobi equation is separable in parabolic coordinates u 1,2 if the Hamilton function has the form H = p 2 1 +p2 2 +V Nq 1,q 2 ), [N/2] V N = 4a 2 1 2k N k k k=0 ) q 2k 1 qn 2k 2, 6

where the positive integer N enumerates the members of the hierarchy. At N = 3 one gets the Hénon-Heiles system 2.8), at N = 4 the next member of hierarchy is a 1:12:16) system with the following Hamiltonian H = p 2 1 +p2 2 4a q 4 1 +12q2 1 q2 2 +16q4 2). 3.1) The corresponding Lax matrix is equal to p 2 2 + p 1q 1 2λ Lλ) = aλ 3 +2aq 2 λ 2 +aq1 2 +4q2 2 )λ+4aq 2q1 2 +2q2 2 )+ p2 1 After similarity transformation of Lλ) with matrix V 2.10), where λ 2q 2 q2 1 λ p 2 2 p 1q 1 2λ. 3.2) ˆL 11 λ) = aλ 2 4 aq 2 q 1 p 1 λ 2 aq 3 1 +2p 1 q 2 p 2 q 1, 2q 1 2q 1 one gets the transformed Lax matrix with two off-diagonal elements ˆL 12 λ) and ˆL 21 λ), which yield two families of variables of separation. As above first coordinates are parabolic coordinates u 1,2, whereas second coordinates v 1,2 are zeroes of the polynomial ˆL 21 = 44aq 1q 2 ap 1 ) λ 2 + 8aq2 1 q2 1 +2q2 2 )+4 aq1 2p 1 q 2 p 2 q 1 ) p 2 1 λ q 1 q 2 1 + 16aq2 1q 2 q 2 1 +2q 2 2)+2p 1 p 1 q 2 p 2 q 1 ) q 2 1 = 44aq 1q 2 ap 1 ) q 1 λ v 1 )λ v 2 ). Substituting roots of this polynomial into the definition 2.6) one gets integrable Hamiltonian with velocity dependent potential H = p2 1 2 +p2 2 +4 ap 1 q 1 q 2 2 ap 2 q 2 1 8aq2 2 5q2 1 +8q2 2 ). Using canonical transformation we can rewrite this Hamiltonian in a more symmetric form H = p 2 1 +p2 2 3 ap 2 q1 2 +2aq4 1 12q2 1 q2 2 32q4 2 ). 3.3) The corresponding second integral of motion is fourth order polynomial in momenta H 2 = p 4 1 +4q4 1 q 4 1 8q1 2 ) q2 2 112q4 2 a 2 +4q1 3 64p1 q2 3 p 2q1 3 12p 2q 1 q2) 2 a 3/2 + q 2 1 4p 2 1 q1 2 48p2 1 q2 2 +32p 1p 2 q 1 q 2 +p 2 ) 2 q2 1 a 6a 1/2 p 2 1 p 2q1 2, which also can be obtained from the Lax matrix ˆLλ). Of course, this integrable system on the plane 3.3) could be obtained in the framework of different theories, see [5, 11, 16] and references within. Canonical transformation 2.15) allows us to identify a Hamilton function with velocity dependent potential 3.3) and Hamilton function Ĥ = P 2 1 +P 2 2 aq 4 1 +6Q 2 1Q 2 2 +Q 4 2) similar to the relation between second and third Hénon-Heiles systems. Canonical transformation u,p u ) v,p v ) is the special auto-bt for the 1:12:16) system, which can be considered as a hetero-bt relating two different Hamilton-Jacobi equations associated with Hamiltonians H 3.1) and H 3.3), respectively. 7

4 Conlusion The problem of finding separation coordinates for the Hamilton-Jacobi equations is highly nontrivial. The problem was originally formulated by Jacobi when he invented elliptic coordinates and successfully applied them to solve several important mechanical problems with quadratic integrals of motion in momenta. We suppose that after suitable Bäcklund transformations standard elliptic, parabolic etc. coordinates turn into variables of separation for physically interesting integrable systems with higher order integrals of motion. For example, in this note we have constructed a canonical transformation of the standard parabolic coordinates, which yields variables of separation for the three integrable Hénon-Heiles systems. Moreover, we believe that information about such suitable Bäcklund transformations and the corresponding integrable systems is incorporated into the Lax matrices associated with these elliptic, parabolic etc. coordinates. In order to prove it we obtained integrals of motion, variables of separation and separated relations for some new integrable system with velocity dependent potential and fourth order integral of motion in momenta. In similar manner we can construct various simultaneously separable integrable systems associated with other curvilinear coordinates on the Riemannian manifolds of constant curvature, see examples in [13, 14]. This work was partially supported by RFBR grant 13-01-00061 and SPbU grant 11.38.664.2013. References [1] Ablowitz, M. J. and Segur, H., Solitons and the Inverse Scattering Transform, SIAM, Philadelphia, 1981. [2] Bobenko, A.I. and Suris,Yu.B., Discrete differential geometry: integrable structure, Graduate Studies in Mathematics, v. 98, AMS, Providence, 2008. [3] Conte, R., Musette, M. and Verhoeven, C., Explicit integration of the Hénon-Heiles Hamiltonians, J. Nonlinear Math. Phys., v.12, Suppl. 1, pp.212-227, 2005. [4] Darboux, G., Leçons sur la théorie générale des surfaces et les applications géométriques du calcul infinitésimal, vol. III Gauthier-Villars, Paris, 1894). [5] Fordy A. P., The Hénon-Heiles system revisited, Phys. D, v.52, no.2-3, pp.204-210, 1991. [6] C. G. J. Jacobi, Vorlesungen über Dynamik, Georg Reimer, Berlin, 1866. [7] Kuznetsov, V.B. and Sklyanin, E.K., On Bäcklund transformations for mahy-body systems, J. Phys. A, v. 31, n.1, pp.2241-2251, 1998. [8] Kuznetsov, V. and Vanhaecke, P., Bäcklund transformations for finite-dimensional integrable systems: a geometric approach, J. Geom. Phys., v. 44, n.1, pp.1-40, 2002. [9] Landau L.D., Lifshitz E.M., Mechanics, v. 1 of A Course of Theoretical Physics, Pergamon Press, 1969. [10] Miller W., Post S. and Winternitz P., Classical and quantum superintegrability with applications, J. Phys. A: Math. and Theor., v.46, n.42, 423001, 2013. [11] Pucacco G., On integrable Hamiltonians with velocity dependent potentials, Celestial Mech. and Dyn. Astronomy, v.90, n.1-2, pp.109-123, 2004. [12] Rogers, C. and Schief, W.K., Bäcklund and Darboux transformations: geometry and modern applications in soliton theory, vol. 30, Cambridge University Press, 2002. [13] Tsiganov, A.V., Simultaneous separation for the Neumann and Chaplygin systems, Reg. Chaotic Dyn., v. 20, n.1, pp.74-93, 2015, 8

[14] Tsiganov, A.V., On the Chaplygin system on the sphere with velocity dependent potential, Journal of Geometry and Physics, v.92, p. 9499, 2015. [15] Wojciechowski S., The analogue of the Bäcklund transformation for integrable many-body systems, J. Phys. A: Math. Gen., v.15, pp.653-657, 1982. [16] Yehia H.M. and Hussein A.M., New families of integrable two-dimensional systems with quartic second integrals, arxiv:1308.1442, 2013. [17] Zullo F., Bäcklund transformations and Hamiltonian flows, J. Phys. A: Math. and Theor., v.46, n.14, 145203, 2013. 9