Testing Newton’s Law of Gravitation with Modern Solar-System Data:
From Tycho Brahe to JPL

Newton’s law of universal gravitation is often presented in a textbook as a finished formula,

(1)

A more interesting question is: how could observations tell us that gravity has this form?

This article reconstructs that logic twice. First, it follows the historical chain from Tycho Brahe’s precise observations, through Kepler’s empirical orbital rules, to Newton’s dynamical interpretation. Second, it repeats several of the key tests with modern Solar-System data from NASA and the Jet Propulsion Laboratory (JPL).

The word “prove” must be used carefully in experimental physics. Observations do not prove a physical law in the deductive sense of a mathematical theorem. They test the law against nature over a specified domain and precision. Modern gravity models are also not purely Newtonian: high-precision planetary ephemerides include many-body perturbations, nonspherical gravity, relativistic corrections, observational biases, and other effects [8, 7]. What we can test very strongly is that Newton’s central inverse-square law is the dominant Solar-System gravitational law and that modern observations recover its characteristic signatures.

The main empirical chain studied here is

(2)

1. The historical path: Tycho, Kepler, and Newton

1.1 Tycho Brahe supplied precision observations, not Newton’s force law

Tycho Brahe assembled an unusually precise naked-eye observational record of planetary positions during the late sixteenth century. Mars was especially important because its relatively large orbital eccentricity makes departures from uniform circular motion easier to detect than for nearly circular planets.

After working with Tycho in Prague and succeeding him as Imperial Mathematician, Johannes Kepler used Tycho’s observations extensively in the Astronomia Nova of 1609. Kepler tested a traditional circular-orbit construction against Mars observations and found an approximately eight-arcminute discrepancy that he regarded as too large to dismiss given the quality of Tycho’s observations. This failure became part of the route to a new orbital description [2, 1].

The historical chain is therefore more precise than the common statement that “Newton used Tycho Brahe’s data.” Newton did not simply perform a force-law fit to Tycho’s raw Mars measurements. Kepler first transformed observations into empirical regularities about planetary motion. Newton later supplied a dynamical framework in which those regularities, together with observations of planetary satellites and other phenomena, could be understood as consequences of centripetal gravitational forces [4, 3].

1.2 Kepler’s observational regularities

In modern language, three Keplerian statements became particularly important:

1. planets approximately follow ellipses with the Sun at one focus;
2. the radius vector from the Sun sweeps equal areas in equal times;
3. the orbital periods and semimajor axes approximately satisfy

   (3)

The first two emerged from Kepler’s Mars work; the harmonic relation was published later in Harmonices Mundi. Historically, Newton treated such astronomical “phenomena” as high-quality approximations rather than exact statements about an isolated two-body universe [4].

1.3 From the area law to a central force

CM04 derived the differential-area relation

(4)

Constant areal velocity means

(5)

Because

(6)

we obtain

(7)

Thus the acceleration must be parallel or antiparallel to the radius vector. The area law is therefore evidence that the dominant acceleration is central: it points along the line joining the planet and the Sun.

This result does not yet tell us how the magnitude varies with distance. For that, the period-distance relation supplies the next clue.

1.4 Kepler’s harmonic law plus circular dynamics gives the inverse square

For uniform circular motion of radius r and period T,

(8)

The centripetal acceleration is

ac	=	(9)
	= 2	(10)
	= .	(11)

If the orbital data obey

(12)

then

(13)

Therefore

(14)

This is the elementary observational route to an inverse-square acceleration. Newton’s actual reasoning was richer than this nearly-circular approximation. In the Principia, he connected the area rule to centripetal force, considered satellite systems as well as planets, and related the detailed absence or presence of apsidal motion to the force law [4, 3].

2. Turn Newton’s inverse square into a parameter the data can estimate

Instead of assuming an inverse-square law from the start, suppose the radial acceleration magnitude follows a generalized power law

(15)

where K and n are unknown. Newtonian point-mass gravity corresponds to

(16)

For a circular orbit,

(17)

Rearranging,

(18)

Take logarithms:

(19)

Equivalently,

(20)

with

(21)

and hence

(22)

This is useful experimentally. We can fit the observed slope β in log–log space and let the orbital data estimate the force-law exponent n.

For elliptical Keplerian orbits the correct characteristic radius is the semimajor axis a, and Newtonian two-body dynamics gives exactly

(23)

For a planet with m ≪ M_⊙, μ ≃ GM_⊙.

3. Modern data used in this article

The first data set contains the eight planets from Mercury through Neptune. The accompanying CSV records J2000-style semimajor axes and sidereal periods from the NASA planetary fact sheets, together with perihelion/aphelion distances and the corresponding maximum/minimum orbital speeds used in the area-law test. NASA’s fact-sheet notes explain that the compiled values draw on sources including JPL ephemerides, astronomical almanacs, IAU/IAG reports, and current literature; the values are research products rather than a single set of immutable “official” constants [5].

The second data set contains eight regular inner/Galilean moons of Jupiter from JPL’s Planetary Satellite Mean Elements page: Metis, Adrastea, Amalthea, Thebe, Io, Europa, Ganymede, and Callisto. JPL states that these mean elements are precessing ellipses fitted in a least-squares sense to numerically integrated satellite orbits, and warns that accurate ephemeris computation should use horizons rather than the mean-element table [6].

4. Modern test I: does the motion look central?

For a central force,

(24)

At perihelion and aphelion the radial velocity is zero, so the velocity is transverse. Therefore

(25)

A central two-body orbit predicts

(26)

Using the rounded perihelion/aphelion distances and maximum/minimum speeds published in the NASA fact sheets, define the diagnostic

(27)

The resulting values are

The root-mean-square mismatch is about

(28)

and the largest absolute mismatch is about

(29)

The result is not a precision test of angular-momentum conservation, because the tabulated values are rounded and correlated products of orbit determination. It is nonetheless a striking classroom-scale demonstration that modern published orbital quantities obey the central-force signature rv = constant at the apsides to roughly one part in 10³ or better.

5. Modern test II: let the planets estimate the force exponent

For the eight planets, perform an ordinary least-squares fit

(30)

The fitted slope from the values in Table 1 is

(31)

The generalized power-law exponent is therefore

n	= 2β − 1	(32)
	= 2(1.4998265333) − 1,	(33)

so

(34)

The fit has

(35)

The formal ordinary-least-squares standard error propagated to n is approximately

(36)

That number is a formal regression error only. It is not a complete experimental uncertainty on the exponent of gravity. The orbital quantities are model-derived, highly correlated, and subject to systematic assumptions that are not represented by an elementary two-column regression.

The scale range is itself notable: the fit spans semimajor axes from about 0.39 au for Mercury to about 30.07 au for Neptune, nearly two orders of magnitude in radius. The leading-order Solar-System relation has the power expected from inverse-square gravity across that entire interval.

6. Modern test III: repeat the exponent test around Jupiter

A strong empirical law should not depend on choosing the Sun as the central body. Newton emphasized satellite systems precisely because the same dynamical logic could be applied around planets [4, 3].

Apply the same regression to JPL’s eight regular Jovian moons in Table 2. The result is

(37)

and hence

(38)

The fit gives

(39)

The formal OLS standard error of n_J is approximately

(40)

again with the important warning that this is not a complete physical uncertainty budget.

The Sun and Jupiter systems therefore tell the same leading-order story:

(41)

This is much closer in spirit to Newton’s claim of universality than a fit to only one planet or one central body.

7. Recover the scale of the law: estimate GM_⊙

The exponent tells us how the acceleration changes with distance. The constant of proportionality tells us the strength of the central gravitational field.

For a Newtonian two-body orbit,

(42)

Therefore

(43)

As a numerical example, use Venus:

(44)

The International Astronomical union defines

(45)

exactly, and recommends that the solar mass parameter be determined observationally in SI units [9].

Converting units and inserting the Venus values gives

(46)

JPL’s DE440 astrodynamic-parameter value is

(47)

The simple Venus two-body reconstruction differs by approximately

(48)

This comparison is educational rather than an independent precision experiment. The published semimajor axis and period are themselves products of astronomical orbit determination and modern dynamical modeling. The agreement shows internal consistency of the two-body scale at an impressive level; it should not be advertised as an independent 0.32-ppm measurement of Newtonian gravity.

8. Why modern ephemerides are a much stronger test than these tables

The preceding calculations deliberately use compact published orbital summaries because they expose the physics. Modern Solar-System gravity tests operate at a much deeper level.

JPL’s DE440 and DE441 planetary/lunar ephemerides were produced by fitting numerically integrated trajectories to ground- and space-based observations. The data include spacecraft radio ranging, Very Long Baseline Array measurements, optical astrometry, lunar laser ranging, occultation data, and other observations. For example, DE440 incorporated Juno ranging and VLBA data to improve Jupiter’s orbit and Cassini ranging and VLBA data to improve Saturn’s orbit [8].

Conceptually, modern orbit determination looks like Figure 5.

At leading order, the acceleration of body i due to point-mass Newtonian interactions is

(49)

The correction term may contain relativistic contributions, nonspherical gravity, asteroid perturbations, tidal effects, and other modeled accelerations appropriate to the data set. The power of the modern test is not that one draws a straight line through eight points; it is that one integrated dynamical model must account simultaneously for enormous collections of heterogeneous observations over long time spans.

9. Newton’s law is extraordinarily successful, but not exact

The inverse-square form

(50)

is the dominant weak-field, low-speed gravitational acceleration in the Solar System. It is not the complete modern theory of gravity. General relativity predicts corrections that become measurable in precise planetary and spacecraft data, including relativistic perihelion advance and propagation-time effects. Modern ephemerides therefore do not test “Newton versus data” in isolation; they test a hierarchy of dynamical models in which Newtonian gravity provides the leading term [8].

This distinction illustrates how physical science works. A law can be spectacularly successful within a domain while being recognized as the limiting form of a more general theory.

10. Reproducible calculation with Julia

The accompanying file

reproduces the numerical calculations using only the Julia standard library. It reads

and

performs the two log–log linear regressions, computes

(51)

checks the apsidal quantity rv, and reconstructs GM_⊙ from the Venus pair (a,T).

For a set of points (x_i,y_i), the code uses the ordinary least-squares slope

(52)

where

(53)

The numerical values printed by the supplied script are recorded in analysis_results.txt. The article’s results were independently cross-checked during preparation. A Julia installation is sufficient to reproduce the calculation; no astronomy package is required because the source data have already been frozen into the supplied CSV files.

10.1 File downloads from the PhysicsLibrary filebox

The reproducibility files used in this article can be attached directly to the PhysicsLibrary object’s filebox. PhysicsLibrary’s manual file-link command is

Once the files have been uploaded to this object’s filebox with the exact filenames shown below, the following links provide direct downloads:

• Download the Julia analysis script
• Download the planetary orbital data CSV
• Download the Jupiter regular-moons data CSV
• Download the saved numerical results

The Julia script and the two CSV files are sufficient to reproduce the numerical tests in this article. The text results file is included as a convenient reference output. The filebox names are intentionally kept simple and should be uploaded without renaming so that the manual links remain valid.

For a higher-fidelity extension, one could query JPL Horizons or use Astropy/Astroquery to obtain epoch-specific state vectors, numerically differentiate or fit trajectories, and compare observation-based residuals with Newtonian and post-Newtonian models. That would move the exercise from an orbital-summary test toward a simplified form of real orbit determination.

11. What has actually been demonstrated?

The modern data used here support four distinct statements:

1. The small apsidal rv mismatches are consistent with the central-force/constant-areal-velocity structure expected from angular-momentum conservation.
2. The eight planetary periods and semimajor axes give

   (54)

   when a generalized central force a ∝ r⁻ⁿ is inferred from the period–distance scaling.

3. An independent fit to eight regular moons of Jupiter gives

   (55)

   Thus the same inverse-square behavior appears around a second central body.

4. The scale relation

   (56)

   recovers the solar gravitational parameter from a simple planetary orbit at a value extremely close to the modern JPL value.

These are not logically independent, raw-observation proofs of Newton’s law: the orbital elements themselves are derived using modern ephemeris models. They are transparent consistency tests that allow a student to reproduce, with modern numbers, the essential observational logic behind Newtonian celestial mechanics.

The most important result is therefore not merely

(57)

It is the methodological chain

(58)

That chain connects Tycho, Kepler, Newton, and modern spacecraft navigation remarkably directly.

12. Suggested extensions

Several natural computational projects follow from this article:

1. Replace the tabulated mean elements by JPL Horizons state vectors sampled over time and numerically verify constant areal velocity.
2. Fit a generalized acceleration law directly to numerically estimated accelerations rather than inferring the exponent through Kepler’s third-law scaling.
3. Repeat the satellite regression for Saturn, Uranus, or Neptune.
4. Construct Newtonian N-body and relativistic propagators and compare their residuals against modern ephemeris states.
5. Use Mercury to investigate why a pure Newtonian ellipse does not remain exactly fixed in inertial space.

These extensions naturally connect this application article to later PhysicsLibrary lessons on perturbation theory, numerical integration, N-body dynamics, and relativistic celestial mechanics.

References

References