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 rotational 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.
The Hénon–Heiles potential can be expressed as
V(x,y)=
1 | |
2 |
(x2+y2)+λ\left(x2y-
y3 | |
3 |
\right).
H=
12 | |
(p |
2 | |
x |
+
2) | |
p | |
y |
+
1 | |
2 |
(x2+y2)+λ\left(x2y-
y3 | |
3 |
\right).
x |
=px,
px |
=-x-2λxy,
y |
=py,
py |
=-y-λ(x2-y2).
λ
\R2
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
i\hbar | \partial |
\partialt |
\Psi(x,y)=\left[
-\hbar2 | |
2m |
\nabla2+
12 | |
(x |
2+y2)+λ\left(x2y-
13 | |
y |
3\right)\right]\Psi(x,y).
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.[2] had shown that in the Hénon–Heiles system the exit basins have the Wada property.