Orbits and Celestial Dynamics

Scrollable formula. Use Left and Right arrow keys to reveal hidden content.

Scrollable table. Use Left and Right arrow keys to reveal additional columns.

Scrollable figure. Use Left and Right arrow keys to inspect hidden content.

Orbital dynamics is where symmetry, conservation laws, turning points, and stability first become a single calculation. A central force confines every nonradial trajectory to a plane; rotational symmetry fixes the angular momentum; and the remaining radial motion is governed by a one-dimensional effective potential. The chapter develops that structure before moving from exact Kepler conics to encounters, weak perturbations, resonances, and restricted few-body dynamics.

The Newtonian theory is more than a historical approximation. It supplies the language of periapsis, binding energy, impact parameter, orbital elements, and stability that later reappears in relativistic calculations. It also provides discriminating weak-field checks. The transfer must nevertheless be made carefully: relativistic geodesic motion derives its parameter, constants of motion, and normalization from a spacetime metric, not from a Newtonian force law.

Helpful background. Calculus and Linear Algebra for Relativity and Differential Equations for Relativity provide vector calculus and second-order initial-value problems; Variational Mechanics and Symmetries supplies the Euler–Lagrange equations and Noether reasoning; The Newtonian Gravitational Field and Potential fixes the inverse-square field and sign conventions.

Central forces as one-dimensional dynamics

For two masses, introduce the relative separation $\mathbf r=\mathbf x_1-\mathbf x_2$ , total mass $M=m_1+m_2$ , and reduced mass $\mu=m_1m_2/M$ . Their Newtonian gravitational interaction reduces to a particle of mass $\mu$ moving in $V(r)=-Gm_1m_2/r$ . If $\ell$ is the conserved angular-momentum magnitude, the radial energy is

E=\frac{\mu}{2}\dot r^2+V_{\mathrm{eff}}(r), \qquad V_{\mathrm{eff}}(r)=V(r)+\frac{\ell^2}{2\mu r^2}.

The angular term is kinetic energy rewritten at fixed $\ell$ , not an additional interaction. Allowed radii satisfy $E\ge V_{\mathrm{eff}}$ ; simple intersections are turning points, while a circular orbit also requires $V'_{\mathrm{eff}}=0$ . The sign of $V''_{\mathrm{eff}}$ then decides linear radial stability, with a vanishing second derivative requiring a higher-order analysis. This reduction and the conic orbit equation are developed in Tong 2013, §§4.1–4.4, pp. 48–66, PDF; the exact reduced-mass separation is given in Tong 2013, §5.1.5, pp. 71–72, PDF.

Keep two mass combinations distinct. The gravitational parameter $GM$ determines the relative Kepler trajectory, whereas $\mu$ converts specific energy and angular momentum into the binary’s total mechanical quantities. Confusing them gives correct-looking orbit shapes with wrong energies and periods.

A route through orbital phenomena

Begin with Central-Force Reduction and Effective Potentials to derive planarity, turning points, and the circular-orbit test. Kepler Orbits and Orbital Elements then obtains the conic

r(\nu)=\frac{p}{1+e\cos\nu}

and connects $p$ , eccentricity, semimajor axis, anomaly variables, and period to conserved quantities. The Newtonian Two-Body Problem reconstructs both barycentric trajectories and makes the roles of $M$ and $\mu$ explicit.

For unbound motion, Newtonian Parabolic and Hyperbolic Encounters relates asymptotic speed and impact parameter to deflection and gravitational focusing. Perturbed Orbits and Apsidal Precession treats orbital elements as slowly varying coordinates on an osculating Kepler orbit and distinguishes periodic changes from secular drift. Finally, Restricted Few-Body Dynamics, Resonances, and Stability introduces the Jacobi integral, Lagrange equilibria, resonant phases, and the boundary between linear and nonlinear stability.

The sequence is deliberately cumulative. Exact two-body integrability establishes the baseline; scattering changes the boundary data; perturbation theory changes constants into slowly varying quantities; and the restricted three-body problem shows how conserved quantities can coexist with resonance and chaos.

Assumptions and transfer to relativity

The chapter assumes Newtonian absolute time, point masses unless stated otherwise, and instantaneous conservative gravity. Tides, finite size, nonspherical multipoles, dissipation, radiation reaction, close collisions, and relativistic corrections all add structure. An effective-potential minimum proves only radial stability within the chosen reduced model; it does not establish stability against arbitrary three-dimensional, dissipative, or many-body perturbations.

Orbital elements are coordinates on initial-data space. For a perturbed trajectory they normally describe the Kepler conic tangent to the true motion at one epoch, not a path that will repeat forever. Numerical propagation must evolve the full force model and test convergence, especially near close encounters, separatrices, and chaotic regions.

Within this volume, continue to Special Relativity and Minkowski Spacetime to replace absolute time with invariant causal geometry before relativistic particle and field dynamics are introduced.

The geometric continuation is Geodesics, Exponential Maps, and Jacobi Fields. There, freely falling curves follow a connection, conserved quantities arise from spacetime symmetries, and geodesic deviation replaces Newtonian relative-acceleration language.

Review the chapter

For a bound Kepler orbit with semimajor axis $a$ and eccentricity $e$ , express the periapsis and apoapsis radii and verify that their sum is independent of $e$ .

Solution

At periapsis and apoapsis the true anomaly is respectively $0$ and $\pi$ . Using $p=a(1-e^2)$ in the conic equation gives

r_{\mathrm p}=\frac{a(1-e^2)}{1+e}=a(1-e), \qquad r_{\mathrm a}=\frac{a(1-e^2)}{1-e}=a(1+e).

Therefore

r_{\mathrm p}+r_{\mathrm a}=2a.

The eccentricity redistributes the fixed major-axis length between the two apsidal radii. The circular limit $e\to0$ gives $r_{\mathrm p}=r_{\mathrm a}=a$ , while $e\to1^-$ drives periapsis toward zero and apoapsis toward $2a$ within the point-particle Kepler model.

References

Tong, David. Dynamics and Relativity. University of Cambridge lecture notes, Lent 2013. §§4.1–4.4, pp. 48–66; §5.1.5, pp. 71–72, PDF.