Advanced Topics
Laplace-Runge-Lenz Vector
For an inverse square law force of the form
and the corresponding potential energy , the super useful Laplace-Runge-Lenz vector, or LRL vector is defined as
and are the linear and angular momentum of the particle at any instant, and is the unit vector pointing from the center of mass to the particle. The constant parameter k describes the strength of the central force; it is equal to for gravitational forces. The LRL vector is a scaled version of the eccentricity vector, , and hence is a constant of motion. The magntidue of the LRL vector is
where E is the total energy of the system.
Trajectory of the momentum vector
The conservation of the LRL vector and angular momentum vector is useful in showing that the momentum vector moves on a circle under an inverse-square central force. The trajectory of the momentum vector is called a hodograph.
Choosing along the -axis, and the major semiaxis as the -axis, yields the locus equation for
Thus, the momentum vector is confined to a circle of radius centered on . For unbounded orbits, and hence the circle does not intersect the -axis.
Precession of LRL vector under a perturbed potential
Consider a small perturbation to the potential energy , such that
Since the force is no longer inverse-square, the LRL vector is no longer conserved, and precesses, leading to precessiong of apsis.
Since and are still conserved, the magnitude of remains constant. Writing ,
where . The average precession of the apsis
The rate of precession is therefore given by
If the perturbation is small, the orbit during one revolution can be approximated as keplerian
and hence the integral can be evaluated.
For example, the potential due to a point mass in general relativity is given by
Evaluating the integral gives
Orbit Equation
The energy of a particle of mass moving in a central potential is
Differentiating with respect to time, and substituting the , we get the equaiton of motion of the particle
Substituting and simplifying, we get
This is the orbit equation. Also referred to as Binet’s equation, it describes the shape of the orbit of a particle moving under a central force. If the shape of the trajectory is known, the force law can be determined by solving for .
For an inverse square force having potential , the orbit equation is given by
This is, again, the equation of a conic section. The solution is
where (the eccentricity) and (the phase offset) are constants of integration.
Let us clarify what bound, closed and stable orbits mean
- A bound orbit is one where the total energy . The particle cannot escape to infinity, and will always remain within a certain distance from the center of force.
- A closed orbit is one where the particle returns to its initial position after a finite period of time. Closed orbits are a subset of bound orbits, but not all bound orbits are closed. For example, in an inverse-square law force, elliptical orbits are closed, but parabolic and hyperbolic orbits are not.
- A stable orbit is one where small perturbations to the particle’s position or velocity do not lead to large deviations from the original orbit. In other words, if the particle is slightly displaced from its orbit, it will oscillate around the original orbit rather than diverging away from it. Stability is a more general concept that can apply to both bound and unbound orbits.
Now, as it turns out, there exist only two types of central force potentials that produce closed orbits for all bound particles: the inverse-square law potential () and the radial harmonic oscillator potential (). This result is known as Bertrand’s theorem.
As far as stability is concerned, all closed orbits must be stable to begin with. Hence, the inverse square law and radial harmonic oscillator potentials produce stable closed orbits.
Aside from these two, there are also closed circular orbits in any central force potential, but these are not stable in general. (in a potential , circular orbits are stable only if )
Daniel, a strange guy, finds a strange force as well. He finds it to be a central force proportional to , where is the distance from the source of the force to the body affected by it (perhaps, think of it as a strange type of “gravity”). What value of could make it such that closed, stable, nearly circular orbits can occur by bodies under the influence of this force (that is, orbits without precession each period)?
A central force proportional to would give a potential of the form
Using the binet equation,
Let the radius of a circular orbit be . For a circular orbit, , which gives
To check for stability, we perturb the orbit a little. Substituting and linearizing (ignoring quadratic and higher order terms), we get
This gives a linear equation for , which can be analyzed for stability.
This is an equation for a simple harmonic motion. Thus for small perturbations, the orbit oscillates about the circular orbit with frequency . Therefore, for stability, we need .
For the orbit to be closed, we must have that the frequency of oscillation is a rational multiple of the orbital frequency. The orbital frequency is 1 (since increases by in one orbit), and the oscillation frequency is . Therefore, for closed orbits, must be a rational number.
The only option which satisfies both these conditions is .
Precession of apasis in General Relativity
The orbit equation in general relativity is given by
Defining a dimensionless and ,
Since is very small, we can use perturbation methods to find an approximate solution for . In the first order, can be expanded in terms of a power series as
just gives the newtonian solution
For , we get
This gives the solution for as
Putting everything together and neglecting the small terms, we get
The period of the ellipse is not , and hence precesses at a rate of
which matches the result which we obtained via the precession of the LRL vector. This apsidal precession is prominent in the orbit of Mercury.

Orbits in General Relativity
As shown earlier, elliptical orbits in general relativity tend to precess. The energy conservation equation is given by
where the central body of mass is much more massive than the orbiting body of mass and is essentially at rest. Here is a term related to the energy of the system, and is the specific angular momentum of the orbiting body. For a circular orbit, , which gives the radius of the orbit as
We see that for a given angular momentum , two circular orbits are possible. Introducing and the schwarschild radius ,
Therefore for a circular orbit, . This puts a limit on the minimum angular momentum of the body.
The inner orbit
is unstable, while the outer orbit
is stable. Therefore the minimum radius of a stable circular orbit is when
In the limit ,
Circular orbits with radius smaller than are not possible. For massless particles (photons), the only circular orbit possible is at .
Gravitational Waves
Consider two massive bodies and orbiting each other, with position vectors and relative to the center of mass. We define and .
According to general theory of relativity, accelerated masses with non zero quadrouple moments radiate gravitataional waves (GWs) and lose energy. For small enough velocities, the emitted GWs
- have a frequency twice as large as the orbital frequency
- can be characterized by a luminosity, which is dominated by the expression
where is the orbital separation and is the angular velocity of each mass.
The energy of the system is . Differentiating, we get
where we define the chirp mass .
Now, gives
where is a constant of integration.
Poynting Robertson Effect
The force of radiation acting on a particle of mass density and radius is , where is the cross section of the particle, is the radiation pressure coefficient, and is the radiation pressure. The radiation pressure is given by
where is the luminosity of the star and is the distance from the star. For sun, the peak wavelength is . For , where sun’s gravity dominates over radiation pressure, the Poynting-Robertson effect acts as a brake on orbiting particles, decreasing their angular momenta so that they slowly spiral into the sun. In the frame of the orbiting particle, the phenomena of abberation causes slight displacement in the direction of the motion for photons striking the particle. This applies a torque on the particle, causing it to lose angular momentum.
The angular momenta of the orbiting particle is
Hence we get that the radius of the particle’s orbit decays as
The Poynting Robertson time scale is the time it takes for the particle to spiral into the sun
where is the initial radius of the particle’s orbit.