Celestial Mechanics: Newton’s Equation for Orbital Motion

CM01 introduced Newton’s law of universal gravitation as a force law. CM02 then expressed the same interaction through gravitational potential and potential energy. CM03 now turns that physics into a dynamical equation: given an object’s position and velocity now, how does gravity determine its future motion?

For a particle moving in the gravity of a spherical body of mass M, the central result is

(1)

This equation is the Newtonian point-mass equation of orbital motion. It is a second-order vector differential equation. Nearly everything developed later in elementary celestial mechanics—orbital energy, angular momentum, conic sections, Kepler’s laws, orbital elements, and time-of-flight relations—will follow from this equation plus suitable initial conditions [1, 2, 3, 4, 5].

A second important result is that for two finite masses m₁ and m₂, the exact relative coordinate obeys

(2)

Thus the same mathematical form appears again, but with gravitational parameter

(3)

The derivation of these equations is the main goal of CM03.

1 From gravitational force to acceleration

Place a spherical gravitating body of mass M at the origin. Let a particle of mass m have position vector

(4)

with magnitude

(5)

Newton’s gravitational force on the particle is

(6)

Newton’s second law is

(7)

Because acceleration is the second time derivative of position,

(8)

Therefore

(9)

Cancel the particle mass m:

(10)

Define the gravitational parameter

(11)

Then the orbital equation becomes

(12)

2 Why the test mass cancels

The gravitational force is proportional to the orbiting mass:

(13)

But the inertial resistance to acceleration in Newton’s second law is also proportional to that same mass:

(14)

Therefore

(15)

contains no remaining factor of the test mass. In Newtonian theory, two sufficiently small test bodies released from the same position with the same initial velocity follow the same trajectory in a prescribed gravitational field, independent of their individual masses.

This cancellation is the orbital-motion version of the result already encountered in CM01:

(16)

The source mass determines the field; the test particle responds to that field.

3 Why the vector equation contains r∕r³

The scalar inverse-square magnitude is

(17)

But an equation of motion must specify direction as well as magnitude. The inward radial unit vector is

(18)

Since

(19)

we can write

a	= −r	(20)
	= −	(21)
	= −μ.	(22)

Thus the apparently unusual denominator r³ is simply the combination of

(23)

for the force magnitude and

(24)

for the radial direction.

4 The minus sign means central attraction

The vector r points from the source at the origin toward the orbiting particle. Gravity points in the opposite direction, back toward the source. Hence

(25)

The acceleration is therefore always antiparallel to r:

(26)

This does not mean the velocity must point inward. In orbital motion the velocity often has a large transverse component. Gravity continuously bends the trajectory because the velocity and acceleration generally point in different directions.

5 Cartesian component equations

Because

(27)

we have

(28)

Substitute into the orbital equation:

(29)

Matching components gives

(30)

(31)

and

(32)

These three coupled scalar differential equations are exactly equivalent to the compact vector equation.

6 A dimensional check

The gravitational parameter has units

(33)

The factor

(34)

has dimensions

(35)

Therefore

(36)

which is the correct dimension for acceleration.

7 Connection with gravitational potential

CM02 derived

(37)

for a point or spherically symmetric source outside its mass distribution, together with

(38)

Since a freely moving test body has

(39)

we can also write the equation of motion as

(40)

For

(41)

this becomes

(42)

Thus the force-law, field, and potential descriptions all lead to the same dynamics.

8 The equation is a differential equation, not an orbit formula

The expression

(43)

specifies acceleration at each position. It does not yet tell us directly whether the resulting trajectory will be a circle, ellipse, parabola, or hyperbola.

To obtain a unique trajectory, we must also specify an initial position and initial velocity:

(44)

The pair (r₀,v₀) determines the particular solution of the orbital differential equation.

Different initial velocities at the same position can produce completely different trajectories even though the gravitational law is identical.

9 First-order state form

Numerical mechanics usually rewrites the second-order equation as two coupled first-order vector equations. Define

(45)

Then

(46)

and

(47)

Define the six-component state

(48)

Then

(49)

This is the form that will later be integrated numerically.

10 Example 1: gravitational acceleration at 400 km altitude

Take Earth’s gravitational parameter as

(50)

Use mean Earth radius

(51)

At altitude

(52)

the geocentric distance is

(53)

The acceleration magnitude is

a	=	(54)
	=	(55)
	≈ 8.69 m/s2.	(56)

Thus

(57)

A spacecraft in low Earth orbit is therefore not in a region where gravity is nearly absent. Its orbital motion is continuous free fall under a gravitational acceleration still close to terrestrial surface gravity.

11 Example 2: evaluate the acceleration vector

Suppose a spacecraft has geocentric position

(58)

Convert to meters:

(59)

The radius is

r	=	(60)
	≈ 7.280 × 106 m.	(61)

The acceleration is

(62)

Using the numerical values gives approximately

(63)

Its magnitude is about

(64)

consistent with

(65)

The component signs are not arbitrary: both are negative because the acceleration points from the spacecraft back toward Earth’s center.

12 Circular motion as one special solution

A circular orbit is not an additional force law. It is one special solution of Newton’s orbital equation.

For uniform circular motion of radius r, the inward acceleration magnitude is

(66)

Gravity provides exactly that acceleration:

(67)

Multiply by r:

(68)

Thus the circular speed is

(69)

This formula is only a circular-orbit special case. The general equation of motion remains

(70)

At 400 km altitude around Earth,

(71)

Thus a spacecraft can experience roughly 8.69 m/s² of inward acceleration while maintaining altitude because its transverse velocity continually carries it around Earth as gravity bends its path.

13 Radial motion as another special case

If the velocity is exactly radial and no transverse component exists, the motion remains on a radial line. Let r now denote a signed radial coordinate along that line outside the source. The equation reduces to

(72)

This is a one-dimensional nonlinear free-fall equation. It is not the general orbital problem because ordinary orbits possess transverse motion.

The contrast is important:

(73)

whereas in a circular orbit

(74)

General Keplerian motion lies between these special cases.

14 Two finite masses: inertial equations

The fixed-source equation is extremely useful when one body overwhelmingly dominates the mass. However, Newtonian gravity acts on both bodies. Let two masses have inertial positions

(75)

Define the relative vector from body 1 to body 2:

(76)

The force on body 2 due to body 1 is

(77)

Therefore

(78)

and cancellation of m₂ gives

(79)

By Newton’s third law, the force on body 1 is opposite:

(80)

Hence

(81)

Equal and opposite forces therefore do not imply equal and opposite accelerations unless the masses themselves are equal.

15 Deriving the exact relative equation of motion

Differentiate the relative position twice:

(82)

Substitute the two body accelerations:

r	= −Gm1 − Gm2	(83)
	= −G(m1 + m2).	(84)

Therefore

(85)

Define the two-body gravitational parameter

(86)

Then the exact relative equation becomes

(87)

This has exactly the same mathematical form as the fixed-central-mass equation.

The distinction is hidden in the definition of μ:

(88)

for the test-particle approximation, but

(89)

for exact Newtonian relative two-body motion.

16 The center of mass moves uniformly

Define the center-of-mass coordinate

(90)

Differentiate twice:

(91)

Using Newton’s equations,

(92)

and

(93)

The internal forces cancel:

(94)

Therefore

(95)

So the center of mass moves at constant velocity when no external force acts on the pair.

CM08 will return to this reduction in greater depth and introduce the reduced mass. For CM03, the essential point is that the internal two-body dynamics can already be represented by one relative vector obeying the familiar central-force equation.

17 Recovering the individual positions from R and r

The definitions

(96)

and

(97)

can be solved for the two body positions:

(98)

(99)

Thus the relative orbit does not imply that either body literally sits at the origin. Both bodies orbit their common center of mass unless one mass is negligible compared with the other.

This point becomes especially important in binary-star and exoplanet applications, where the reflex motion of the more massive body is directly observable.

18 The test-particle limit

Suppose

(100)

Then

(101)

so

(102)

The center of mass also lies very close to m₁. In that limit, treating the massive body as fixed at the origin becomes an excellent approximation.

For a planet orbiting the Sun, or a spacecraft orbiting Earth, this approximation is often extremely accurate. For comparable-mass binary stars, it is not.

19 First structural consequence: gravity is a central acceleration

The orbital equation has the form

(103)

with

(104)

Therefore

(105)

This fact will be central in CM04. It implies zero torque about the force center and leads directly to conservation of specific orbital angular momentum:

(106)

CM03 does not yet need that conservation law to establish the dynamics, but the structure of Newton’s orbital equation already contains it.

20 What information is needed to predict an orbit?

For a known gravitational parameter μ, the Newtonian point-mass orbit is determined by six scalar initial conditions:

(107)

and

(108)

Equivalently, we specify

(109)

Those six numbers are the Cartesian state of the body at an epoch.

Later, CM15 will show that the same ideal two-body orbit can be described geometrically by six classical orbital elements. The state-vector description and orbital-element description are two different coordinate descriptions of the same physical trajectory.

21 Astrophysical interpretation

The equation

(110)

is not merely a spacecraft equation. It also describes the idealized relative motion of

• two stars in a binary system,
• a star and an exoplanet,
• a planet and a moon,
• a comet and the Sun,
• or a star orbiting a compact central mass when other perturbations are negligible.

The same equation therefore underlies both celestial mechanics and much of observational astrophysics. Once an orbit can be measured, its geometry and timing can be used to infer otherwise inaccessible masses.

For example, in an exact Newtonian binary, the relative equation depends on

(111)

That is why measured binary orbits can reveal the total system mass. Later astrophysical application lessons will develop this connection quantitatively.

22 Where the ideal equation stops being sufficient

The equation

(112)

assumes an isolated Newtonian two-body system with spherical or point-like masses. Real systems may require additional accelerations:

(113)

where a_p represents perturbations such as

• gravity from additional bodies,
• nonspherical gravity fields,
• atmospheric drag,
• radiation pressure,
• tidal interactions,
• or relativistic corrections.

The ideal equation is nevertheless the correct starting point because perturbation theory is built by understanding this solvable core first.

23 Common mistakes

• Writing r = −μ∕r² without a direction. The left-hand side is a vector, so the right-hand side must also be a vector.
• Confusing −μr∕r³ with an inverse-cube force. Its magnitude is still μ∕r² because |r| = r.
• Forgetting that the position vector must be measured from the chosen force center in the fixed-central-body approximation.
• Treating velocity as if it always points toward the gravitating body. Gravity constrains acceleration, not velocity direction.
• Assuming low Earth orbit has negligible gravity. Orbital free fall occurs precisely because gravity remains strong.
• Using μ = Gm₁ for a comparable-mass binary. Exact relative motion uses μ = G(m₁ + m₂).
• Assuming equal and opposite gravitational forces produce equal accelerations. Accelerations differ when masses differ.
• Thinking the differential equation alone selects one orbit. Initial position and velocity are also required.
• Confusing a state vector (r,v) with orbital elements. They are different representations of the same ideal orbit.

24 What CM03 adds to the series

CM01 supplied the force law

(114)

CM02 supplied the potential description

(115)

CM03 now converts those descriptions into dynamics:

(116)

For a fixed dominant source,

(117)

while exact two-body relative motion uses

(118)

Together with

(119)

this defines the Newtonian orbital initial-value problem.

The next lessons can now ask what deeper structure is hidden inside this differential equation. The first answers are conservation of angular momentum and conservation of orbital energy, which will eventually lead to the conic-section orbit equation.

References