Skip to content

Parameterized Post-Newtonian Gravity

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.

The parameterized post-Newtonian (PPN) framework is a common language for slowly moving sources in weak, quasistatic gravitational fields. It replaces theory-specific field equations by coefficients in a standardized metric, allowing many Solar-System and laboratory observables to test the same combinations. PPN is a framework for a regime—not a complete theory and not a description of relativistic merger or cosmological evolution.

Required background. Weak Metric Perturbations and Gauge Freedom supplies weak-field gravity, perturbative metric variables, and gauge control. Probability and Inference for Gravity supplies likelihoods, covariance, and parameter inference.

Helpful background. Relativistic Two-Body Dynamics supplies the orbital expansions whose weak-field observables are encoded by PPN coefficients.

Choose the standard PPN gauge and count a characteristic velocity as v=O(ϵ)v=O(\epsilon), the Newtonian potential as U=O(ϵ2)U=O(\epsilon^2), and pressure-to-density and specific internal energy as O(ϵ2)O(\epsilon^2). With the site signature (−,+,+,+)(-,+,+,+) and

U(x,t)=G∫ρ(x′,t)∣x−x′∣ d3x′,U(\boldsymbol x,t)=G\int\frac{\rho(\boldsymbol x',t)} {\lvert\boldsymbol x-\boldsymbol x'\rvert}\,d^3x',

the simplest part of the PPN metric is

g00=−1+2U−2βU2+⋯ ,gij=(1+2γU)δij+⋯ .g_{00}=-1+2U-2\beta U^2+\cdots, \qquad g_{ij}=(1+2\gamma U)\delta_{ij}+\cdots.

The omitted terms include standardized scalar and vector potentials formed from density, velocity, pressure, and internal energy. Their coefficients define ten conventional PPN parameters. General relativity has γ=β=1\gamma=\beta=1 and the remaining eight equal to zero in this convention Will 2014, §3.2, Box 2, and Table 2.

The coefficient γ\gamma measures how much spatial curvature is produced per unit rest mass in this weak-field limit; β\beta measures nonlinear superposition in g00g_{00}. The preferred-frame parameters α1,α2,α3\alpha_1,\alpha_2,\alpha_3, preferred-location parameter ξ\xi, and nonconservation parameters ζ1,…,ζ4\zeta_1,\ldots,\zeta_4 organize other ways a metric theory can differ. These labels describe phenomenology. Their values are derived by solving a specified theory and matching it to the PPN potentials.

Light propagation at leading post-Newtonian order depends on 1+γ1+\gamma. Both deflection and Shapiro delay therefore compare null trajectories with an independently modeled mass distribution, reference frame, and propagation environment. Planetary perihelion advance contains the combination

2+2γ−β3\frac{2+2\gamma-\beta}{3}

relative to the GR post-Newtonian contribution, in addition to multipoles, asteroid masses, solar mass loss, and other ephemeris parameters. When preferred-frame, preferred-location, and nonconservation parameters all vanish, the Nordtvedt combination controlling differential free fall of self-gravitating bodies is

ηN=4β−γ−3.\eta_{\rm N}=4\beta-\gamma-3.

In the full standard PPN framework it is instead

ηN=4β−γ−3−103ξ−α1+23α2−23ζ1−13ζ2.\eta_{\rm N}=4\beta-\gamma-3-\frac{10}{3}\xi -\alpha_1+\frac23\alpha_2-\frac23\zeta_1-\frac13\zeta_2.

Energy-momentum conservation alone does not remove the preferred-location term ξ\xi. Quoting the shorter expression without its assumptions silently imposes a theory prior Will 2014, §4.3.1, Eq. (66).

An observable rarely measures one coefficient alone. A global fit may constrain correlated combinations together with masses, initial conditions, reference-frame orientation, ephemeris parameters, and calibration quantities. The likelihood should use the predicted observations—not post-fit residuals treated as new independent data. Modern planetary-ephemeris analyses emphasize that changing gravity changes the fitted ephemeris itself Fienga and Minazzoli 2024, §§2–5.

To compute PPN parameters, expand the theory’s fields and matter equations consistently to the required order, solve the standardized potentials with declared boundary conditions, and transform the metric to the standard PPN gauge. The matter coupling determines which metric governs clocks and freely falling test bodies. A conformal or disformal field redefinition cannot be interpreted by reading coefficients before the operational matter frame is fixed.

The mapping can depend on a cosmological scalar value or screening environment. If the local solution has a thin shell, a Vainshtein region, or strong-field body sensitivities, constant universal PPN coefficients may not describe all bodies and locations. In that case one must solve the environmental problem or use a framework designed for body-dependent charges.

The independent algebraic check is the GR limit: substitute γ=β=1\gamma=\beta=1 and all other PPN parameters zero, then recover the standard post-Newtonian metric and observables. A second check is gauge consistency: a coordinate transformation may rearrange intermediate potentials but cannot change the timed, ranged, or angular observable.

PPN assumes weak fields, slow motion, an appropriate metric description, and sufficiently smooth near-zone behavior. It does not describe gravitational radiation through merger, a black-hole ringdown spectrum, or arbitrary horizon-scale cosmological perturbations. It also does not by itself impose stability, a high-energy cutoff, or correlations predicted by a complete action.

A bound on γ−1\gamma-1 should consequently be reported with the dataset, global model, assumed values or priors for other coefficients, and the theory map. Combining “one parameter at a time” bounds as if they were a joint likelihood can create a parameter region that no analysis actually tested.

In the conservative PPN subset, an observation constrains 2+2γ−β2+2\gamma-\beta while another constrains 1+γ1+\gamma. Show why the pair can determine both γ\gamma and β\beta, and identify one reason the algebra may not equal a valid joint inference.

Solution

Let the first measured combination be A=2+2γ−βA=2+2\gamma-\beta and the second B=1+γB=1+\gamma. Then γ=B−1\gamma=B-1 and substitution gives β=2+2(B−1)−A=2B−A\beta=2+2(B-1)-A=2B-A. Thus independent exact measurements would determine both coefficients. In real analyses the estimates may share ephemerides, masses, calibration, or priors, so their covariance is not generally zero. Therefore a valid combination uses the joint likelihood or a documented covariance matrix; multiplying two marginalized one-dimensional curves can double-count shared information and impose inconsistent assumptions about the remaining PPN parameters.