Astronomy:Hénon–Heiles system

From HandWiki
Contour plot of the Hénon–Heiles potential

While at Princeton in 1962, Michel Hénon and Carl Heiles worked on the non-linear motion of a star around a galactic center with the motion restricted to a plane. In 1964 they published an article titled "The applicability of the third integral of motion: Some numerical experiments".[1] Their original idea was to find a third integral of motion in a galactic dynamics. For that purpose they took a simplified two-dimensional nonlinear axi-symmetric potential and found that the third integral existed only for a limited number of initial conditions. In the modern perspective the initial conditions that do not have the third integral of motion are called chaotic orbits.

Introduction

The Hénon–Heiles potential can be expressed as[2]

[math]\displaystyle{ V(x,y) = \frac{1}{2}(x^2 + y^2) + \lambda \left(x^2y - \frac{y^3}{3}\right). }[/math]

The Hénon–Heiles Hamiltonian can be written as

[math]\displaystyle{ H = \frac12 (p_x^2 + p_y^2) + \frac{1}{2}(x^2 + y^2) + \lambda \left(x^2y - \frac{y^3}{3}\right). }[/math]

The Hénon–Heiles system (HHS) is defined by the following four equations:

[math]\displaystyle{ \dot{x} = p_x, }[/math]
[math]\displaystyle{ \dot{p_x} = -x - 2\lambda xy, }[/math]
[math]\displaystyle{ \dot{y} = p_y, }[/math]
[math]\displaystyle{ \dot{p_y} = -y - \lambda(x^2 - y^2). }[/math]

In the classical chaos community, the value of the parameter [math]\displaystyle{ \lambda }[/math] is usually taken as unity. Since HHS is specified in [math]\displaystyle{ \R^2 }[/math], we need a Hamiltonian with 2 degrees of freedom to model it. It can be solved for some cases using Painlevé analysis.

Quantum Hénon–Heiles Hamiltonian

In the quantum case the Hénon–Heiles Hamiltonian can be written as a two-dimensional Schrödinger equation.

The corresponding two-dimensional Schrödinger equation is given by

[math]\displaystyle{ i\hbar\frac{\partial}{\partial t} \Psi(x,y) = \left[ \frac{-\hbar^2}{2m}\nabla^2 + \frac12 (x^2 + y^2 ) + \lambda \left(x^2 y - \frac13 y^3\right) \right] \Psi(x,y). }[/math]

Wada property of the exit basins

Hénon–Heiles system shows rich dynamical behavior. Usually the Wada property cannot be seen in the Hamiltonian system, but Hénon–Heiles exit basin shows an interesting Wada property. It can be seen that when the energy is greater than the critical energy, the Hénon–Heiles system has three exit basins. In 2001 M. A. F. Sanjuán et al.[3] had shown that in the Hénon–Heiles system the exit basins have the Wada property.

References

  1. Hénon, M.; Heiles, C. (1964). "The applicability of the third integral of motion: Some numerical experiments". The Astronomical Journal 69: 73–79. doi:10.1086/109234. Bibcode1964AJ.....69...73H. 
  2. Hénon, Michel (1983), "Numerical exploration of Hamiltonian Systems", in Iooss, G., Chaotic Behaviour of Deterministic Systems, Elsevier Science Ltd, pp. 53–170, ISBN 044486542X 
  3. Aguirre, Jacobo; Vallejo, Juan C.; Sanjuán, Miguel A. F. (2001-11-27). "Wada basins and chaotic invariant sets in the Hénon-Heiles system". Physical Review E (American Physical Society (APS)) 64 (6): 066208. doi:10.1103/physreve.64.066208. ISSN 1063-651X. 

External links