A pendulum, no matter the gravitational strength (*g*) or the length of the string (*l*), finds its north star in the dimensionless oscillation period: $T/l/g$ is always $2\pi $. In understanding this statement, students seem to pass through similar stages as I did. First, I didn't even understand the claim. Second, I understood it only symbolically, as I would on hearing that my child was running a fever of 40 °C: “Yes, and?” Third and finally, I felt its force, as I would on hearing that the fever was 104 °F: “Oh, damn, that's dangerous!” I became amazed that I could swing my key chain from a string, measure the period, compute $T/l/g$, and get a number valid for pendula everywhere and of all lengths.

This lengthy process of understanding happened for a dimensionless constant, a zero-parameter universal or dimensionless function. Still harder to grasp is the next-simplest kind of universal function, a one-parameter dimensionless function. A famous example, featured in many treatments of dimensional analysis, is Taylor's analysis of the Trinity atomic-bomb test.^{1}

A less warlike example hides among the motion of the planets. Take the orbital data in Table I that include perihelion ($rmin$) and aphelion ($rmax$) distances. The data could serve for a course that I've long wanted to teach: “Patterns in nature as seen in the *CRC Handbook of Chemistry and Physics*.” If the mask …^{2} … lifts, I might even teach such a course in person. Each lecture would begin with CRC data on one topic—perhaps the boiling points of long-chain hydrocarbons, the temperature versus height in the atmosphere, or the thermal conductivity of several materials—and explain the trends qualitatively and quantitatively, approximating as needed.

. | $rmin$ (m) . | $rmax$ (m) . | V ($m2/s2$)
. |
---|---|---|---|

Mercury | $4.6001\xd71010$ | $6.9818\xd71010$ | $\u22121.1462\xd7109$ |

Venus | $1.0748\xd71011$ | $1.0894\xd71011$ | $\u22126.1339\xd7108$ |

Earth | $1.4710\xd71011$ | $1.5210\xd71011$ | $\u22124.4369\xd7108$ |

Mars | $2.0666\xd71011$ | $2.4923\xd71011$ | $\u22122.9119\xd7108$ |

Jupiter | $7.4067\xd71011$ | $8.1601\xd71011$ | $\u22128.5277\xd7107$ |

Saturn | $1.3498\xd71012$ | $1.5036\xd71012$ | $\u22124.6523\xd7107$ |

Uranus | $2.7350\xd71012$ | $3.0063\xd71012$ | $\u22122.3122\xd7107$ |

Neptune | $4.4598\xd71012$ | $4.5370\xd71012$ | $\u22121.4755\xd7107$ |

. | $rmin$ (m) . | $rmax$ (m) . | V ($m2/s2$)
. |
---|---|---|---|

Mercury | $4.6001\xd71010$ | $6.9818\xd71010$ | $\u22121.1462\xd7109$ |

Venus | $1.0748\xd71011$ | $1.0894\xd71011$ | $\u22126.1339\xd7108$ |

Earth | $1.4710\xd71011$ | $1.5210\xd71011$ | $\u22124.4369\xd7108$ |

Mars | $2.0666\xd71011$ | $2.4923\xd71011$ | $\u22122.9119\xd7108$ |

Jupiter | $7.4067\xd71011$ | $8.1601\xd71011$ | $\u22128.5277\xd7107$ |

Saturn | $1.3498\xd71012$ | $1.5036\xd71012$ | $\u22124.6523\xd7107$ |

Uranus | $2.7350\xd71012$ | $3.0063\xd71012$ | $\u22122.3122\xd7107$ |

Neptune | $4.4598\xd71012$ | $4.5370\xd71012$ | $\u22121.4755\xd7107$ |

In these orbital data, *V* is the planet's effective potential: its effective potential energy, $Ueff$, divided by its mass, *m*. $Ueff$ itself is most evident in the total energy written in polar coordinates:

where *M* is the solar mass. Then,

If the orbital data are plotted with *V* on the vertical axis and distances ($rmin$ and $rmax$) on the horizontal axis, the result is a mess with no interesting pattern (Fig. 1). Sure, the outer planets are nearer the horizontal axis than are the inner planets, and all the points lie near a curve that might be $V\u221d\u22121/r$. But this pattern arises from the form of the gravitational potential energy, $\u2212GMm/r$, so it's not news.

But let's replot the data in dimensionless form. To rescale the horizontal axis, we need a length scale. The best scale is

the harmonic mean of $rmin$ and $rmax$. (We'll return to why we should form that particular length.) The rescaled horizontal axis is then $r\xaf\u2261r/r0$.

For rescaling the vertical axis, we need a specific energy (energy per mass). A quantity with these dimensions and based on the chosen length scale is $GM/r0$. Thus, the rescaled vertical axis will be $V\xaf\u2261V/(GM/r0)$.

The resulting dimensionless plot (Fig. 2) is cleaner than the preceding unscaled graph. The data now lie on the decent curve

As promised, let's return to the choice of length scale: Why is the harmonic mean, *r*_{0}, the correct length scale? Other natural length scales include the semi-major axis, *a*, and the semi-minor axis, *b*; see Fig. 3, taken with minor changes from my column on zeroth powers.^{3}

One reason not to choose the semi-major axis as the length scale is that the resulting plot is *too* regular. The total energy in a Newtonian-gravitational or Kepler orbit turns out to be

At the orbital extremes (where $r\u0307=0$), Eq. (1) shows us that $E=Ueff$. Then,

Thus, at the extremes, any planet has $V/(GM/a)=\u22121/2$, placing all the data on one horizontal line (with this choice of length scale).

Meanwhile, the semi-minor axis and indeed all length scales except *r*_{0} are disfavored for a physical reason. The orbital energy, Eq. (1), depends on *r* and $\theta \u0307$. In using energy conservation to find the orbital shape, we would get a two-variable differential equation (from $dE/dt=0$), far more than twice as difficult as a single-variable differential equation. Fortunately, the energy's joker term, $mr2\theta \u03072/2$, which ropes in the second variable, can be rewritten using the other conserved quantity in a Kepler orbit, the angular momentum

Thus,

and voilà: The $\theta \u0307$ has gone, making the effective potential a function now only of *r* (and ostensibly of *L*, but it's a constant).

This rewriting helps us because *r*_{0}, unlike other length scales, is a function just of *L* (and physical constants). In a Kepler orbit, it turns out that

If we wisely choose *r*_{0} as our length scale, the dimensionless plot in Fig. 2 becomes an analog computer that implicitly incorporates the conservation of angular momentum.

Let's put the plot and the curve Eq. (4) to a moderate test by using an orbital extreme case—Halley's comet.^{4} Here are its orbital data:

Do its two points lie on the dimensionless curve of Eq. (4)? Numerically,

For $r\xaf=rmax/r0$, Eq. (4) predicts

and the $V\xaf$ calculated from Eq. (11) is almost identical, −0.032 456. (The reader is invited to confirm that the closest approach, where $r\xaf=0.508\u2009388$, also works.) In Newtonian gravity, any planet orbiting around any sun is described by one universal function.

Having passed through the stages of understanding a zero-parameter universal function and having thought for a while about the orbital one-parameter function, I feel myself nearing the final, third stage of understanding and my amazement growing.

## AUTHOR DECLARATIONS

### Conflict of Interest

The author has no conflicts to disclose.

*Sanjoy Mahajan is interested in the art of approximation and physics education and has taught varying subsets of physics, mathematics, electrical engineering, and mechanical engineering at MIT, the African Institute for Mathematical Sciences, and the University of Cambridge. He is the author of Street-Fighting Mathematics (MIT Press, 2010), The Art of Insight in Science and Engineering (MIT Press, 2014), and A Student's Guide to Newton's Laws of Motion (Cambridge University Press, 2020).*