In classical mechanics, Bertrand's theorem states that only two types of central force potentials produce stable, closed orbits: (1) an inverse-square central force such as the gravitational or electrostatic potential
and (2) the radial harmonic oscillator potential
All attractive central forces can produce circular orbits, which are naturally closed orbits. The only requirement is that the central force exactly equals the centripetal force, which determines the required angular velocity for a given circular radius. Non-central forces (i.e., those that depend on the angular variables as well as the radius) are ignored here, since they do not produce circular orbits in general.
where and the angular momentum L = mr2ω is conserved. For illustration, the first term on the left-hand side is zero for circular orbits, and the applied inwards force equals the centripetal force requirement mrω2, as expected.
The definition of angular momentum allows a change of independent variable from t to θ
giving the new equation of motion that is independent of time
This equation becomes quasilinear on making the change of variables and multiplying both sides by (see also Binet equation)
As noted above, all central forces can produce circular orbits given an appropriate initial velocity. However, if some radial velocity is introduced, these orbits need not be stable (i.e., remain in orbit indefinitely) nor closed (repeatedly returning to exactly the same path). Here we show that stable, exactly closed orbits can be produced only with an inverse-square force or radial harmonic oscillator potential (a necessary condition). In the following sections, we show that those force laws do produce stable, exactly closed orbits (a sufficient condition).
Define J(u) as
where f represents the radial force. The criterion for perfectly circular motion at a radius r0 is that the first term on the left-hand side be zero
Substituting this expansion into the equation for u and subtracting the constant terms yields
which can be written as
where is a constant. β2 must be non-negative; otherwise, the radius of the orbit would vary exponentially away from its initial radius. (The solution β = 0 corresponds to a perfectly circular orbit.) If the right-hand side may be neglected (i.e., for small perturbations), the solutions are
where the amplitude h1 is a constant of integration. For the orbits to be closed, β must be a rational number. What's more, it must be the same rational number for all radii, since β cannot change continuously; the rational numbers are totally disconnected from one another. Using the definition of J along with equation (1),
where is evaluated at . Since this must hold for any value of u0,
which implies that the force must follow a power law
Hence, J must have the general form
For more general deviations from circularity (i.e., when we cannot neglect the higher order terms in the Taylor expansion of J), η may be expanded in a Fourier series, e.g.,
We substitute this into equation (2) and equate the coefficients belonging to the same frequency, keeping only the lowest order terms. As we show below, h0 and h2 are smaller than h1, being of order . h3, and all further coefficients, are at least of order . This makes sense since must all vanish faster than h1 as a circular orbit is approached.
From the cos(βθ) term, we get
where in the last step we substituted in the values of h0 and h2.
Using equations (3) and (1), we can calculate the second and third derivatives of J evaluated at u0,
Substituting these values into the last equation yields the main result of Bertrand's theorem
Hence, the only potentials that can produce stable, closed, non-circular orbits are the inverse-square force law (β = 1) and the radial harmonic oscillator potential (β = 2). The solution β = 0 corresponds to perfectly circular orbits, as noted above.
Inverse-square force (Kepler problem)
The orbit u(θ) can be derived from the general equation
whose solution is the constant plus a simple sinusoid
where e (the eccentricity) and θ0 (the phase offset) are constants of integration.
This is the general formula for a conic section that has one focus at the origin; e = 0 corresponds to a circle, e < 1 corresponds to an ellipse, e = 1 corresponds to a parabola, and e > 1 corresponds to a hyperbola. The eccentricity e is related to the total energy E (cf. the Laplace–Runge–Lenz vector)
Radial harmonic oscillator
To solve for the orbit under a radial harmonic oscillator potential, it's easier to work in components r = (x, y, z). The potential energy can be written
The equation of motion for a particle of mass m is given by three independent Lagrange's equations
where the constant must be positive (i.e., k > 0) to ensure bounded, closed orbits; otherwise, the particle will fly off to infinity. The solutions of these simple harmonic oscillator equations are all similar
where the positive constants Ax, Ay and Az represent the amplitudes of the oscillations and the angles φx, φy and φz represent their phases. The resulting orbit r(t) = [x(t), y(y), z(t)] is closed because it repeats exactly after a period
The system is also stable because small perturbations in the amplitudes and phases cause correspondingly small changes in the overall orbit.