# Lagrangian Mechanics Made Simple Part 4 Orbital Motion and Kepler's Equation

Ron Steinke <rsteinke@w-link.net>

Last month we introduced the problem of a system two particles whose potential energy depended only on the distance between them. The dynamics of this system reduced to

 (1)

where the quantities and are the angular momentum and energy of the system, and
 (2)

is the reduced mass.

# Kepler's Equation

For the force of gravity, 1the potential is

 (3)

where is a constant. The dependence of the potential reflects the dependence of the gravitational force. The potential is negative, since the particles have less potential energy when they are close together than when they are far apart. Inserting (3) in (1) gives
 (4)

The solution to this equation is
 (5)

This solution is unique up to a constant shift in the angle , .

Eq. (5) describes a conic section,

 (6)

The quantity
 (7)

is called the eccentricity of the orbit. An eccentricity of 0 gives a circular orbit, between 0 and 1 gives an ellipse, 1 gives a parabola, and greater than 1 gives a hyperbola. Notice that corresponds to . A parabolic orbit has just enough energy to escape the system. Higher (positive) energies give hyperbolic orbits, and lower (negative) energies give circular or elliptical orbits, where the two particles are bound together by the gravitational force. Note also that for an orbit of zero angular momentum, (5) gives a constant value for , and the two bodies fall straight into each other.

# The Lagrange Stability Points

The three body problem is not solvable in the general case. However, some approximate solutions are quite interesting. Consider the case of a small body of mass orbiting two larger bodies of masses and . For simplicity we will assume that the two large bodies orbit each other circularly with angular velocity , and are unaffected by their attraction to the smaller body. If and are the distances of the two bodies from their common center of rotation, then

 (8)

where is the distance between the two bodies,
 (9)

Here, is Newton's gravitational constant.

The two larger bodies lie at positions and in a coordinate system with angular velocity . The Lagrangian for the lighter particle at position is

 (10)

where we have used the Lagrangian for a rotating coordinate system, which we derived in Part 2 of this series. To simplify this a bit, we use the fact that (8) and (9) tell us that
 (11)

This substitution gives the Lagrangian the simpler form
 (12)

Notice that all terms are now proportional to . This constant factor will drop out of the equations of motion, and all dependence on the masses of the particles will be parameterized by the lengths and .

To find metastable points, we look for solutions where is zero. We do this by examining the effective potential,

 (13)

which is just with set to zero. Examining , we see that when is near or , the potential is dominated by gravitational attraction, so the light particle will tend to fall towards one of the heavy ones. Similarly, for large , the potential is dominated by the centripetal term, and the light particle will tend to fly off to infinity.

There exist points for which the centripetal acceleration is balanced by the attraction of the larger bodies. At these points, the force,

 (14)

is zero. The component of has no centripetal term, so it only vanishes when is zero. Setting the other two components of the force equal to zero gives us equations for and ,
 (15)

where we have dropped the common factor of .

It is possible to solve the second equation in (15) by setting to zero. The value of is then found by solving the first equation. There are three of these solutions, one with , between the two larger bodies, and one each with and , off to each side. These three points are saddle points. If you perturb the small mass off one of these points by a small amount along the direction of its orbit, it will tend to fall back into place. If you perturb it into a higher or lower orbit, however, it will tend to move away from the stability point, either towards one of the two heavier bodies, or pulled away by the centripetal component of the potential (13).

There are also two solutions where is not equal to zero. These are

 (16)

These are what are commonly known as the fourth and fifth Lagrange points, and are known to be stable. However, a quick examination of the force (14) shows that they do not derive this stability from the effective potential. Let be the two solutions (16), and examine the force at a point , where is small. The components of the force are
 (17)

The force in the direction is proportional to , indicating that it is pulling against the displacement. Any small perturbation of a particle in the direction will result in a force which tends to pull it back to . The same is not true of perturbations in the and directions, however. The effect of and on the force are determined by a portion of the matrix in (17). The eigenvalues of this submatrix are
 (18)

which are both positive. Therefore, for small perturbations in the and directions, the force (14) tends to push particles away from the Lagrange points. The stability of orbits at these points is due to the Coriolis force, which we have neglected. We'll discuss how that works next month.

# About this document ...

Lagrangian Mechanics Made Simple
Part 4
Orbital Motion and Kepler's Equation

#### Footnotes

... gravity, 1
This section leans heavily on J. Marion and S. Thornton, Classical Dynamics of Particles and Systems, 3rd ed., pp. 256-260

Ron Steinke 2002-08-11