A discussion by John Baez prompts the question of whether Sundman transformations can give insight into systems other than the Kepler potential. A Sundman transformation is a reparametrization of the temporal variable that was first used as part of an analytic solution to the three-body problem. Such transformations provide an intriguing perspective on the two-body problem for any spherically symmetric potential.
Begin with the Hamiltonian for an n-dimensional spherically symmetric system,
where the variables have their usual definitions as vector squares:
The dot product of these two vectors
is the third rotationally invariant quantity available in a spherically symmetric system. The square of the angular momentum tensor is given by
which can be used to replace the momentum variables conjugate to spatial angles in the Hamiltonian in the usual manner:
The equation of motion for the radial variable is
To determine completely the motion in the invariant plane one needs an angle with respect to a fixed direction. This is given by angular momentum conservation,
which is easily integrated for
With a known functional dependence of the radial variable on the temporal variable this in principle provides a complete solution to the orbital motion for any spherically symmetric potential:
The problem of course is that the equation of motion for the radial variable, the energy equation
is rarely integrable and invertible. The reparametrization of the temporal variable via a Sundman transformation allows one to choose a simple form for the dependence of the radial variable on the new parameter, obviating the issue of integration and inversion.
The solution presented here will be determined in parametric form up to the two quadratures directly above. There will remain the issue of inverting the dependence of the temporal variable on the new parameter, but this same issue occurs in the Kepler potential in the form of the Kepler equation linking time with the eccentric anomaly.
The Sundman transformation will be written
where following Baez a prime will be used to denote differentiation with respect to the transformed temporal variable s. A dot will be used to denote differentiation with respect to the original temporal variable t and derivatives with respect to the radial variable will be written explicitly.
The first derivative of the radial vector with respect to the transformed temporal variable is
The first derivative of the radial variable with respect to the transformed temporal variable is
which in terms of the original dynamic variables can be variously written
The second derivative of the radial vector with respect to the transformed temporal variable is
The second derivative of the radial vector with respect to the original temporal variable is given by the equations of motion
so that the second derivative of the radial vector with respect to the transformed temporal variable becomes
The second derivative of the radial variable with respect to the transformed temporal variable can be written
but will be more useful in the form
These two forms are easily shown to be identical using the original dynamic variables.
The third derivative of the radial vector with respect to the transformed temporal variable is
where the second derivative of the reparametrization function is
One way to determine the reparametrization function is to remove the last term in the third derivative of the radial vector with respect to the transformed temporal variable by choosing
so that the third derivative of the radial vector becomes
Since the reparametrization function is here a function of the first derivative of the potential, the quantity in brackets is in general complicated, but there are two very simple particular solutions available. If the the reparametrization function is constant, so that both derivatives in the brackets disappear, then
which corresponds to the potential
This is the potential for the simple harmonic oscillator, and the Sundman transformation in this case is merely a rescaling of the temporal variable. If the reparametrization function is linear in the radial variable, so that the second derivative in the brackets disappears, one has
The potential in this case is
which is the Kepler potential. The second derivative of the radial variable with respect to the transformed temporal variable is now
and the third derivative of the radial vector with respect to the transformed temporal variable is
For both the simple harmonic oscillator and the Kepler potential, the Sundman transformation here will produce harmonic motion for the radial vector. This will occur as well for the potential determined by the general solution of the nonlinear differential equation
which is third order with respect to derivatives of the potential. It will not, however, be true for an arbitrary spherically symmetric potential, so that this Sundman transformation is inadequate for a general solution.
To find a Sundman transformation applicable to all spherically symmetric potentials, note that for the Kepler potential the previous transformation led to the behavior
for the radial variable. Since energy is negative for bound states of the Kepler potential, this represents harmonic motion with a constant offset. Taking a cue from this, consider a Sundman transformation that produces harmonic motion in the radial variable for all potentials. This can be done by setting the second derivative of the radial variable with respect to the transformed temporal variable in its more useful form equal to a linear with constant offset:
With the expression above for the first derivative of the radial variable with respect to the transformed temporal variable, the integration is simple:
Since the radial variable is harmonic by construction, one can immediately write
where parameters analogous to the semimajor axis and eccentricity of the Kepler problem have been introduced. The turning points of the motion are
and correspond to zeroes of the momentum conjugate to the radial variable. Evaluating the energy at both end points
one can easily determine both the angular momentum
and the energy
in terms of the orbit parameters introduced. The transformation as it now stands determines the original temporal variable as
with dynamic constants of the motion as given directly above. While this holds for any constant of integration , this constant can be chosen to simplify the integral. Evaluating part of the numerator under the radical,
the choice cancels the constant term, with the remainder being the square of the derivative of the radial variable with respect to the transformed temporal variable. The integral for the temporal variable becomes
which surprisingly is the original energy equation, but now parametrically dependent on the transformed temporal variable. While at first sight it appears that nothing has really been accomplished, first write out the second derivative of the radial variable with respect to the original temporal variable:
The final right-hand side is a fully determined function of the radial variable, with no leeway to assign a convenient parametric form for the radial variable. The introduction of the Sundman transformation provides this leeway in a consistent manner, allowing any dependence of the radial variable to be chosen. The simple harmonic motion has a significant advantage in that it characterizes the turning points efficiently for evaluation of the energy and angular momentum in terms of parameters introduced.
The complete solution includes the integral for the angle in the invariant plane
as a function parametrically dependent on the transformed temporal variable as well. The integrals can be evaluated with respect to either the radial variable or the transformed temporal variable to determine the motion parametrically. This can be done explicitly for only a handful of potentials, but represents a complete solution in integral form.
The parametric determination of the motion is an extension of the classical eccentric anomaly, which is normally described by circumscribing the elliptical orbit of the Kepler potential. It should more accurately be described with a circular orbit halfway between the inner and outer limits of the actual, and generally precessing, orbit. The normal description holds because the eccentric anomaly sweeps on the same angle on all circles centered with the elliptical orbit. The extension of the eccentric anomaly is labeled “Keplerization” in a conference paper by Vladimir Martinusi and Pini Gurfil.
In illustration of the method, orbits in terms of the transformed temporal variable will be evaluated for two simple cases, the Kepler potential and the simple harmonic oscillator. Using expressions above, the energy and angular momentum for the Kepler potential are
The radical for denominators of integrals is
The integral for the temporal variable is
which is the Kepler equation as expected. The integral for the angle in the invariant plane is
which can be verified by differentiation. Circular functions in the invariant plane are
Multiplying by the parametrization of the radial variable, the radial vector for the Kepler potential is
For the two independent coordinates in the invariant plane, the orbit is
which is an ellipse displaced along the major axis as expected. When thinking about this system in terms of a fourth dimension, the eccentric circle is rotated into the fourth dimension at an angle determined by the eccentricity.
The energy and angular momentum for the simple harmonic oscillator are
Integrals for this systems are more easily evaluated in terms of the radial variable to avoid elliptic functions. The radical for denominators of integrals needs to be written in two different ways:
Using the first form of the radical, the integral for the temporal variable is
where a constant of integration has been included for later use. Using the second form of the radical, the integral for the angle in the invariant plane is
Circular functions in the invariant plane are
Multiplying by the parametrization of the radial variable, the radial vector for the simple harmonic oscillator is
For the two independent coordinates in the invariant plane, the orbit is
which is an ellipse centered at the origin. As the parameter s progresses from zero to π only one quadrant of the ellipse is traversed, so that the entire ellipse requires the parameter to vary from zero to 4π. This differs from a period of 2π for the Kepler orbit and explains the difference between them in the location of the center of force.
Unlike the Kepler equation, the transformed temporal variable for the simple harmonic oscillator is easily inverted. Adjusting the constant of integration so that zeroes of both variables coincide, one has
so that the radial vector becomes
and the moderately complicated intermediate parametrization leads directly back to the expected simple circular functions.
As an aside, it is easy to show from the expression for energy that the Kepler potential is the only power potential for which energy is independent of eccentricity.
To recap, the Sundman transformation
along with expressions for the energy
and the angular momentum
determines harmonic motion in the radial variable
for any spherically symmetric potential in conjunction with the two quadratures
The motion is determined parametrically in an extension of the classical eccentric anomaly.
Uploaded 2015.08.15 — Updated 2015.08.17 analyticphysics.com