# 3-body Coulomb dynamics

We have to introduce some basics before study helium atom.
See [1-3] for details or go directly to the "Collinear helium dynamics"
at the bottom to skip the math.
# Classical hydrogen atom

Consider at first classical dynamics of the hydrogen atom. It is well
known that a bounded electron moves along a Kepler's elliptic orbit with
nucleus in one of its focal points.

If the electron has momentum pointing directly to the nucleus then the
ellipse degenerates in a line segment and electron oscillates along the
line. In this case we get a system of one degree of freedom with the Hamiltonian

*H = p*^{2}/2m - e^{2}/r = E

where *E* is the total energy. We put *m = e = 1* further.
# Scaling

Kepler dynamics remains invariant under a change of energy up to a simple
scaling transformation; a solution of the equations of motion at an arbitrary
energy *E < 0* can be transformed into a solution at a fixed energy
*E*_{o} = -1 by scaling the coordinates as

*r(E) = r / (-E),
p(E) = (-E)*^{1/2} p

together with a time transformation

*t(E) = t / (-E)*^{-3/2} .
# Regularization

*H(p,r)* is singular at *r = 0* therefore whenever the electron comes
close to the nucleus accelerations become infinitely large. A regularization
of the two-body collinear collisions is achieved [1] by means of the Levi-Civita
transformation, which consists of a coordinate dependent time transformation,
which stretches the time scale near the origin, and a canonical transformation
of the phase space coordinates.
A time transformation *dt = f(p,q)ds* for a system described by
*H(q,p) = E* leaves the Hamiltonian structure of the
equations of motion unaltered, if the Hamiltonian itself is transformed into
*H' = f(q,p)[H(p,q) - E]*. We choose *dt = r ds* which lifts the
*|r| -> 0* singularity and leads to a new Hamiltonian

*H' = rp*^{2}/2 - 1 - Er = 0 .

Then a canonical transformation of form

*r = Q*^{2}, p = P/2Q

maps the Kepler problem into that of a harmonic oscillator with
Hamiltonian

*H(Q,P) = P*^{2}/8 - EQ^{2} = 1 .

# Collinear helium

Collinear *eZe* helium with the two electrons situated along a line
on opposite sides of the nucleus [1-3] is a system of two degrees of freedom
with the Hamiltonian

*H = p*_{1}^{2}/2 +
p_{2}^{2}/2 - 2/r_{1} - 2/r_{2} +
1/(r_{1} + r_{2}) = E

We put

*r*_{1} = Q_{1}^{2},
r_{2} = Q_{2}^{2},
p_{1} = P_{1}/2Q_{1},
p_{2} = P_{2}/2Q_{2},
R_{12}^{2} = Q_{1}^{2} +
Q_{2}^{2},

ds = dt/r_{1}r_{2} .

Then the new Hamiltonian is

*H = (Q*_{2}^{2}P_{1}^{2} +
Q_{1}^{2}P_{2}^{2})/8 -
2R_{12}^{2} +
Q_{1}^{2}Q_{2}^{2}(
1/R_{12}^{2} - E) = 0

and equations of motion are

*P*_{1}' = 2Q_{1}[2 -
P_{2}^{2}/8 + Q_{2}^{2}(E -
Q_{2}^{2}/R_{12}^{4})]

P_{2}' = 2Q_{2}[2 -
P_{1}^{2}/8 + Q_{1}^{2}(E -
Q_{1}^{2}/R_{12}^{4})]

Q_{1}' = P_{1}Q_{2}^{2}/4

Q_{2}' = P_{2}Q_{1}^{2}/4

where prime denotes derivative with respect to fictitious time *s*.
Scaling transformations are now

*Q*_{i}(E) = Q_{i} / (-E)^{1/2} ,
P_{i}(E) = P_{i} ,
s(E) = (-E)^{1/2} s .

Note, that *P*_{i} does not depend on *E* and
*|P*_{i}| -> 4 when *Q*_{i} -> 0 .

# Leapfrog algorithm

For the first-order equations

**r**' = **v** , **v**' = **F**(**v**)

the simple leapfrog algorithm is

**r**_{1/2} = **r**_{0} +
**v**_{0} *dt/2*

**v**_{1} = **v**_{0} +
**F**(**r**_{1/2}) *dt*

**r**_{1} = **r**_{1/2} +
**v**_{1} *dt/2*

where subscripts *0* and *1* denote the values at the beginning
and end of a step respectively.
# Collinear helium dynamics

The overall dynamics depends critically on whether *E > 0* or *E < 0* .
If the energy is positive both electrons can escape to infinity. More
interestingly, a single electron can still escape even if *E* is negative,
carrying away an unlimited amount of kinetic energy, as the total energy of the
remaining inner electron has no lower bound. Not only that, but one electron
will escape eventually for almost all starting conditions!
To the left below you can see dynamics of the two electrons on the
phase plane

*x = r*^{1/2} = |Q| ,

y = 2pr^{1/2} = P sign(Q) .

To the right the same orbit is ploted on the
*(r*_{1}^{1/2}, r_{2}^{1/2}) plane.

An animation stops when one of *|P*_{i}| is more then *7* .
Click mouse to set new coordinates of electrons. Press "Enter" to set new
*ds* or "Delay".

[1] *P.Cvitanovic et al.*
Classical and Quantum
Chaos Chap.24: Helium atom

[2] *Gregor Tanner, Klaus Richter and Jan-Michael Rost*

"The theory of two-electron atoms: between ground state and complete
fragmentation"

Rev.Mod.Phys. **72**, 497 (2000)

[3] *Klaus Richter, Gregor Tanner and Dieter Wintgen*

"Classical mechanics of two-electron atoms"
Phys.Rev.A **48**, 4182 (1993)

Contents
Previous: Area-preserving Henon map
Next: Poincare maps

*updated* 9 Apr 2003