View source for Parameterized post-Newtonian formalism

{{Short description|Calculational tool concerning energy, linking general relativity to classical physics}}
{{General relativity sidebar |equations}}
{{Being merged|spacetype=|Non-relativistic gravitational fields||section=|inactive=|date=April 2023|nocat=|dir=from}}
In [[physics]], precisely in the study of the [[General relativity|theory of general relativity]] and many [[Alternatives to general relativity|alternatives to it]], the [[Post-Newtonian expansion|post-Newtonian formalism]] is a calculational tool that expresses [[Einstein field equations|Einstein's (nonlinear) equations of gravity]] in terms of the lowest-order deviations from [[Newton's law of universal gravitation]]. This allows approximations to [[Albert Einstein|Einstein]]'s equations to be made in the case of weak fields. Higher-order terms can be added to increase accuracy, but for strong fields, it may be preferable to solve the complete equations numerically. Some of these post-Newtonian approximations are expansions in a small parameter, which is the ratio of the velocity of the matter forming the gravitational field to the [[speed of light]], which in this case is better called the [[speed of gravity]]. In the limit, when the fundamental speed of gravity becomes infinite, the post-Newtonian expansion reduces to Newton's law of [[gravity]].

The '''parameterized post-Newtonian formalism''' or '''PPN formalism''', is a version of this formulation that explicitly details the parameters in which a general theory of gravity can differ from Newtonian gravity.  It is used as a tool to compare Newtonian and Einsteinian gravity in the limit in which the [[gravitational field]] is weak and generated by objects moving slowly compared to the speed of light. In general, PPN formalism can be applied to all metric theories of gravitation in which all bodies satisfy the Einstein [[equivalence principle]] (EEP). The speed of light remains constant in PPN formalism and it assumes that the [[Metric tensor (general relativity)|metric tensor]] is always symmetric.

== History ==
The earliest parameterizations of the post-Newtonian approximation were performed by Sir [[Arthur Stanley Eddington]] in 1922. However, they dealt solely with the vacuum gravitational field outside an isolated spherical body. [[Ken Nordtvedt]] (1968, 1969) expanded this to include seven parameters in papers published in 1968 and 1969. [[Clifford Martin Will]] introduced a stressed, continuous matter description of celestial bodies in 1971.

The versions described here are based on [[Wei-Tou Ni]] (1972), Will and Nordtvedt (1972), [[Charles W. Misner]] et al. (1973) (see [[Gravitation (book)|''Gravitation'' (book)]]), and Will (1981, 1993) and have ten parameters.

== Beta-delta notation ==

Ten '''post-Newtonian parameters''' completely characterize the weak-field behavior of the theory. The formalism has been a valuable tool in [[tests of general relativity]]. In the notation of Will (1971), Ni (1972) and Misner et al. (1973) they have the following values:
{| class="wikitable"
|-
|<math>\gamma</math> ||How much space curvature <math>g_{ij}</math> is produced by unit rest mass?
|-
|<math>\beta</math> ||How much [[Nonlinear system|nonlinearity]] is there in the [[Superposition principle|superposition]] law for gravity <math>g_{00}</math>?
|-
|<math>\beta_1</math> ||How much gravity is produced by unit [[kinetic energy]] <math>\textstyle\frac12\rho_0v^2</math>?
|-
|<math>\beta_2</math> ||How much gravity is produced by unit [[gravitational energy|gravitational potential energy]] <math>\rho_0/U</math>?
|-
|<math>\beta_3</math> ||How much gravity is produced by unit internal energy <math>\rho_0\Pi</math>?
|-
|<math>\beta_4</math> ||How much gravity is produced by unit pressure <math>p</math>?
|-
|<math>\zeta</math> ||Difference between radial and transverse kinetic energy on gravity
|-
|<math>\eta</math> ||Difference between radial and transverse stress on gravity
|-
|<math>\Delta_1</math> ||How much dragging of [[Inertial frame of reference|inertial frames]] <math>g_{0j}</math> is produced by unit [[momentum]] <math>\rho_0v</math>?
|-
|<math>\Delta_2</math> ||Difference between radial and transverse momentum on dragging of inertial frames
|}
<math>g_{\mu\nu}</math> is the 4 by 4 symmetric metric tensor with indexes <math>\mu</math> and <math>\nu</math> going from 0 to 3. Below, an index of 0 will indicate the time direction and indices <math>i</math> and <math>j</math> (going from 1 to 3) will indicate spatial directions.

In Einstein's theory, the values of these parameters are chosen (1) to fit Newton's Law of gravity in the limit of velocities and mass approaching zero, (2) to ensure [[conservation of energy]], [[Conservation of mass|mass]], [[Conservation of momentum|momentum]], and [[Conservation of Angular Momentum|angular momentum]], and (3) to make the equations independent of the [[Frame of reference|reference frame]].  In this notation, general relativity has PPN parameters <math>\gamma=\beta=\beta_1=\beta_2=\beta_3=\beta_4=\Delta_1=\Delta_2=1</math> and <math>\zeta=\eta=0.</math>

== Alpha-zeta notation ==

In the more recent notation of Will & Nordtvedt (1972) and Will (1981, 1993, 2006) a different set of ten PPN parameters is used.
: <math>\gamma=\gamma</math>
: <math>\beta=\beta</math>
: <math>\alpha_1=7\Delta_1+\Delta_2-4\gamma-4</math>
: <math>\alpha_2=\Delta_2+\zeta-1</math>
: <math>\alpha_3=4\beta_1-2\gamma-2-\zeta</math>
: <math>\zeta_1=\zeta</math>
: <math>\zeta_2=2\beta+2\beta_2-3\gamma-1</math>
: <math>\zeta_3=\beta_3-1</math>
: <math>\zeta_4=\beta_4-\gamma</math>
: <math>\xi</math> is calculated from <math>3\eta=12\beta-3\gamma-9+10\xi-3\alpha_1+2\alpha_2-2\zeta_1-\zeta_2</math>

The meaning of these is that <math>\alpha_1</math>, <math>\alpha_2</math> and <math>\alpha_3</math> measure the extent of [[preferred frame]] effects. <math>\zeta_1</math>, <math>\zeta_2</math>, <math>\zeta_3</math>, <math>\zeta_4</math> and <math>\alpha_3</math> measure the failure of conservation of energy, momentum and angular momentum.

In this notation, general relativity has PPN parameters
: <math>\gamma=\beta=1</math> and <math>\alpha_1=\alpha_2=\alpha_3=\zeta_1=\zeta_2=\zeta_3=\zeta_4=\xi=0</math>

The mathematical relationship between the metric, metric potentials and PPN parameters for this notation is:
: <math>\begin{matrix}g_{00} = -1+2U-2\beta U^2-2\xi\Phi_W+(2\gamma+2+\alpha_3+\zeta_1-2\xi)\Phi_1 +2(3\gamma-2\beta+1+\zeta_2+\xi)\Phi_2 \\ \ +2(1+\zeta_3)\Phi_3+2(3\gamma+3\zeta_4-2\xi)\Phi_4-(\zeta_1-2\xi)A-(\alpha_1-\alpha_2-\alpha_3)w^2U \\ \ -\alpha_2w^iw^jU_{ij}+(2\alpha_3-\alpha_1)w^iV_i+O(\epsilon^3) \end{matrix}</math>
: <math>g_{0i}=-\textstyle\frac12(4\gamma+3+\alpha_1-\alpha_2+\zeta_1-2\xi)V_i-\textstyle\frac12(1+\alpha_2-\zeta_1+2\xi)W_i
-\textstyle\frac12(\alpha_1-2\alpha_2)w^iU-\alpha_2w^jU_{ij}+O(\epsilon^{\frac52})</math>
: <math>g_{ij}=(1+2\gamma U)\delta_{ij}+O(\epsilon^2)</math>
where repeated indexes are summed.  <math>\epsilon</math> is on the order of potentials such as <math>U</math>, the square magnitude of the coordinate velocities of matter, etc. <math>w^i</math> is the velocity vector of the PPN coordinate system relative to the mean rest-frame of the universe. <math>w^2=\delta_{ij}w^iw^j</math> is the square magnitude of that velocity. <math>\delta_{ij}=1</math> if and only if <math>i=j</math>, <math>0</math> otherwise.

There are ten metric potentials, <math>U</math>, <math>U_{ij}</math>, <math>\Phi_W</math>,  <math>A</math>, <math>\Phi_1</math>, <math>\Phi_2</math>, <math>\Phi_3</math>, <math>\Phi_4</math>,  <math>V_i</math> and  <math>W_i</math>, one for each PPN parameter to ensure a unique solution. 10 linear equations in 10 unknowns are solved by inverting a 10 by 10 matrix. These metric potentials have forms such as:
: <math>U(\mathbf{x},t)=\int{\rho(\mathbf{x}',t)\over|\mathbf{x}-\mathbf{x}'|}d^3x'</math>
which is simply another way of writing the Newtonian gravitational potential,
: <math>U_{ij}=\int{\rho(\mathbf{x}',t)(x-x')_i(x-x')_j\over|\mathbf{x}-\mathbf{x}'|^3}d^3x'</math>
: <math>\Phi_W=\int{\rho(\mathbf{x}',t)\rho(\mathbf{x}'',t)(x-x')_i\over|\mathbf{x}-\mathbf{x}'|^3}\left({(x'-x'')^i\over|\mathbf{x}-\mathbf{x}'|}-{(x-x'')^i\over|\mathbf{x}'-\mathbf{x}''|}\right)d^3x'd^3x''</math>
: <math>A=\int{\rho(\mathbf{x}',t)\left(\mathbf{v}(\mathbf{x}',t)\cdot(\mathbf{x}-\mathbf{x}')\right)^2\over|\mathbf{x}-\mathbf{x}'|^3}d^3x'</math>
: <math>\Phi_1=\int{\rho(\mathbf{x}',t)\mathbf{v}(\mathbf{x}',t)^2\over|\mathbf{x}-\mathbf{x}'|}d^3x'</math>
: <math>\Phi_2=\int{\rho(\mathbf{x}',t)U(\mathbf{x}',t)\over|\mathbf{x}-\mathbf{x}'|}d^3x'</math>
: <math>\Phi_3=\int{\rho(\mathbf{x}',t)\Pi(\mathbf{x}',t)\over|\mathbf{x}-\mathbf{x}'|}d^3x'</math>
: <math>\Phi_4=\int{p(\mathbf{x}',t)\over|\mathbf{x}-\mathbf{x}'|}d^3x'</math>
: <math>V_i=\int{\rho(\mathbf{x}',t)v(\mathbf{x}',t)_i\over|\mathbf{x}-\mathbf{x}'|}d^3x'</math>
: <math>W_i=\int{\rho(\mathbf{x}',t)\left(\mathbf{v}(\mathbf{x}',t)\cdot(\mathbf{x}-\mathbf{x}')\right)(x-x')_i\over|\mathbf{x}-\mathbf{x}'|^3}d^3x'</math>
where <math>\rho</math> is the density of rest mass, <math>\Pi</math> is the internal energy per unit rest mass, <math>p</math> is the pressure as measured in a local freely falling frame momentarily comoving with the matter, and <math>\mathbf{v}</math> is the coordinate velocity of the matter.

Stress-energy tensor for a perfect fluid takes form
: <math>T^{00}=\rho(1+\Pi+\mathbf{v}^2+2U)</math>
: <math>T^{0i}=\rho(1+\Pi+\mathbf{v}^2+2U+p/\rho)v^i</math>
: <math>T^{ij}=\rho(1+\Pi+\mathbf{v}^2+2U+p/\rho)v^iv^j+p\delta^{ij}(1-2\gamma U)</math>

== How to apply PPN ==

Examples of the process of applying PPN formalism to alternative theories of gravity can be found in Will (1981, 1993). It is a nine step process:
* Step 1: Identify the variables, which may include: (a) dynamical gravitational variables such as the metric <math>g_{\mu\nu}</math>, scalar field <math>\phi</math>, vector field <math>K_\mu</math>, tensor field <math>B_{\mu\nu}</math> and so on; (b) prior-geometrical variables such as a flat background metric <math>\eta_{\mu\nu}</math>, cosmic time function <math>t</math>, and so on; (c) matter and non-gravitational field variables.
* Step 2: Set the cosmological boundary conditions. Assume a homogeneous isotropic cosmology, with isotropic coordinates in the rest frame of the universe. A complete cosmological solution may or may not be needed. Call the results <math>g^{(0)}_{\mu\nu} = \operatorname{diag}(-c_0,c_1,c_1,c_1)</math>, <math>\phi_0</math>, <math>K^{(0)}_\mu</math>, <math>B^{(0)}_{\mu\nu}</math>.
* Step 3: Get new variables from <math>h_{\mu\nu}=g_{\mu\nu}-g^{(0)}_{\mu\nu}</math>, with <math>\phi-\phi_0</math>, <math>K_\mu-K^{(0)}_\mu</math> or <math>B_{\mu\nu}-B^{(0)}_{\mu\nu}</math> if needed.
* Step 4: Substitute these forms into the field equations, keeping only such terms as are necessary to obtain a final consistent solution for <math>h_{\mu\nu}</math>. Substitute the perfect fluid stress tensor for the matter sources.
* Step 5: Solve for <math>h_{00}</math> to <math>O(2)</math>. Assuming this tends to zero far from the system, one obtains the form <math>h_{00}=2\alpha U</math> where <math>U</math> is the Newtonian gravitational potential and <math>\alpha</math> may be a complicated function including the gravitational "constant" <math>G</math>. The Newtonian metric has the form <math>g_{00}=-c_0+2\alpha U</math>, <math>g_{0j}=0</math>, <math>g_{ij}=\delta_{ij}c_1</math>. Work in units where the gravitational "constant" measured today far from gravitating matter is unity so set <math>G_{\mathrm{today}} = \alpha/c_0 c_1=1</math>.
* Step 6: From linearized versions of the field equations solve for <math>h_{ij}</math> to <math>O(2)</math> and <math>h_{0j}</math> to <math>O(3)</math>.
* Step 7: Solve for <math>h_{00}</math> to <math>O(4)</math>. This is the messiest step, involving all the nonlinearities in the field equations. The stress–energy tensor must also be expanded to sufficient order.
* Step 8: Convert to local quasi-Cartesian coordinates and to standard PPN gauge.
* Step 9: By comparing the result for <math>g_{\mu\nu}</math> with the equations presented in [[Parameterized post-Newtonian formalism#Alpha-zeta notation|PPN with alpha-zeta parameters]], read off the PPN parameter values.

== Comparisons between theories of gravity ==
{{main | Alternatives to general relativity#Parametric post-Newtonian parameters for a range of theories}}

A table comparing PPN parameters for 23 theories of gravity can be found in [[Alternatives to general relativity#Parametric post-Newtonian parameters for a range of theories]].

Most metric theories of gravity can be lumped into categories. [[Scalar theories of gravitation]] include conformally flat theories and stratified theories with time-orthogonal space slices.

In conformally flat theories such as [[Nordström's theory of gravitation]] the metric is given by <math>\mathbf{g}=f\boldsymbol{\eta}</math> and for this metric <math>\gamma=-1</math>, which drastically disagrees with observations. In stratified theories such as [[Yilmaz theory of gravitation]] the metric is given by <math>\mathbf{g}=f_1\mathbf{d}t \otimes \mathbf{d} t +f_2\boldsymbol{\eta}</math> and for this metric <math>\alpha_1=-4(\gamma+1)</math>, which also disagrees drastically with observations.

Another class of theories is the quasilinear theories such as [[Whitehead's theory of gravitation]]. For these <math>\xi=\beta</math>.  The relative magnitudes of the harmonics of the Earth's tides depend on <math>\xi</math> and <math>\alpha_2</math>, and measurements show that quasilinear theories disagree with observations of Earth's tides.

Another class of metric theories is the [[bimetric theory]]. For all of these <math>\alpha_2</math> is non-zero. From the precession of the solar spin we know that <math>\alpha_2 < 4\times 10^{-7}</math>, and that effectively rules out bimetric theories.

Another class of metric theories is the [[scalar–tensor theory|scalar–tensor theories]], such as [[Brans–Dicke theory]]. For all of these, <math>\gamma=\textstyle\frac{1+\omega}{2+\omega}</math>. The limit of <math>\gamma-1<2.3\times10^{-5}</math> means that <math>\omega</math> would have to be very large, so these theories are looking less and less likely as experimental accuracy improves.

The final main class of metric theories is the vector–tensor theories. For all of these the gravitational "constant" varies with time and <math>\alpha_2</math> is non-zero. Lunar laser ranging experiments tightly constrain the variation of the gravitational "constant" with time and <math>\alpha_2 < 4\times 10^{-7}</math>, so these theories are also looking unlikely.

There are some metric theories of gravity that do not fit into the above categories, but they have similar problems.

== Accuracy from experimental tests ==
Bounds on the PPN parameters from Will (2006) and Will (2014)
{| class="wikitable"
|-
! Parameter !! Bound !! Effects !! Experiment
|-
| <math>\gamma-1</math> || {{val|2.3|e=-5}} || Time delay, light deflection || [[Cassini–Huygens#Tests of general relativity|Cassini tracking]]
|-
| <math>\beta-1</math> || {{val|8|e=-5}} || Perihelion shift || Perihelion shift
|-
| <math>\beta-1</math> || {{val|2.3|e=-4}} || Nordtvedt effect with assumption <math>\eta_N=4\beta-\gamma-3</math> || [[Nordtvedt effect]]
|-
| <math>\xi</math> ||  {{val|4|e=-9}} || Spin precession || Millisecond pulsars
|-
| <math>\alpha_1</math> || {{val|1|e=-4}} || Orbital polarization || [[Lunar Laser Ranging experiment|Lunar laser ranging]]
|-
| <math>\alpha_1</math> || {{val|4|e=-5}} || Orbital polarization || PSR J1738+0333
|-
| <math>\alpha_2</math> || {{val|2|e=-9}} || Spin precession || Millisecond pulsars
|-
| <math>\alpha_3</math> || {{val|4|e=-20}} || Self-acceleration || Pulsar spin-down statistics
|-
| <math>\eta_N</math> || 9{{e|-4}} || Nordtvedt effect || Lunar laser ranging
|-
| <math>\zeta_1</math> || 0.02 || Combined PPN bounds || —
|-
| <math>\zeta_2</math> || 4{{e|-5}}{{px2}}<sup>†</sup> || Binary-pulsar acceleration || [[PSR 1913+16]]
|-
| <math>\zeta_3</math> || 1{{e|-8}} || Newton's 3rd law || Lunar acceleration
|-
| <math>\zeta_4</math> || 0.006{{px2}}<sup>‡</sup> || — || Kreuzer experiment
|}
<sup>†</sup> {{cite journal |last=Will |first=C. M. |title=Is momentum conserved? A test in the binary system PSR 1913 + 16 |journal=Astrophysical Journal Letters |issn=0004-637X |volume=393 |issue=2 |date=10 July 1992 |page=L59–L61 |doi=10.1086/186451 |bibcode=1992ApJ...393L..59W }}{{br}}
<sup>‡</sup> Based on <math>6\zeta_4 = 3\alpha_3 + 2\zeta_1 - 3\zeta_3</math> from Will (1976, 2006). It is theoretically possible{{clarify|reason=Under what conditions?|date=August 2019}} for an alternative model of gravity to bypass this bound, in which case the bound is <math>|\zeta_4| < 0.4</math> from Ni (1972).

== See also ==
{{Portal|Physics}}
* [[Alternatives to general relativity#Parametric post-Newtonian parameters for a range of theories]]
* [[Effective one-body formalism]]
* [[Linearized gravity]]
* [[Peskin–Takeuchi parameter]] The same thing as PPN, but for [[electroweak theory]] instead of [[gravitation]]
* [[Tests of general relativity]]

== References ==

* Eddington, A. S. (1922) The Mathematical Theory of Relativity,  Cambridge University Press.
* Misner, C. W., Thorne, K. S. & Wheeler, J. A.  (1973) Gravitation, W. H. Freeman and Co.
* {{cite journal | last=Nordtvedt | first=Kenneth | title=Equivalence Principle for Massive Bodies. II. Theory | journal=Physical Review | publisher=American Physical Society (APS) | volume=169 | issue=5 | date=1968-05-25 | issn=0031-899X | doi=10.1103/physrev.169.1017 | pages=1017–1025| bibcode=1968PhRv..169.1017N }}
* {{cite journal | last=Nordtvedt | first=K. | title=Equivalence Principle for Massive Bodies Including Rotational Energy and Radiation Pressure | journal=Physical Review | publisher=American Physical Society (APS) | volume=180 | issue=5 | date=1969-04-25 | issn=0031-899X | doi=10.1103/physrev.180.1293 | pages=1293–1298| bibcode=1969PhRv..180.1293N }}
* {{cite journal | last=Will | first=Clifford M. | title=Theoretical Frameworks for Testing Relativistic Gravity. II. Parametrized Post-Newtonian Hydrodynamics, and the Nordtvedt Effect | journal=The Astrophysical Journal | publisher=IOP Publishing | volume=163 | year=1971 | issn=0004-637X | doi=10.1086/150804 | page=611-628| bibcode=1971ApJ...163..611W }}
* {{cite journal | last=Will | first=C. M. | title=Active mass in relativistic gravity - Theoretical interpretation of the Kreuzer experiment | journal=The Astrophysical Journal | publisher=IOP Publishing | volume=204 | year=1976 | issn=0004-637X | doi=10.1086/154164 | page=224-234| bibcode=1976ApJ...204..224W | doi-access=free }}
* Will, C. M. (1981, 1993) Theory and Experiment in Gravitational Physics, Cambridge University Press. {{ISBN|0-521-43973-6}}.
* Will, C. M., (2006) The Confrontation between General Relativity and Experiment, https://web.archive.org/web/20070613073754/http://relativity.livingreviews.org/Articles/lrr-2006-3/
* {{cite journal | last=Will | first=Clifford M. | title=The Confrontation between General Relativity and Experiment | journal=Living Reviews in Relativity | volume=17 | issue=1 | date=2014-06-11 | issn=2367-3613 | doi=10.12942/lrr-2014-4 | page=4| pmid=28179848 | pmc=5255900 | arxiv=1403.7377 | bibcode=2014LRR....17....4W |doi-access=free}}
* {{cite journal | last1=Will | first1=Clifford M. | last2=Nordtvedt | first2=Kenneth Jr. | title=Conservation Laws and Preferred Frames in Relativistic Gravity. I. Preferred-Frame Theories and an Extended PPN Formalism | journal=The Astrophysical Journal | publisher=IOP Publishing | volume=177 | year=1972 | issn=0004-637X | doi=10.1086/151754 | page=757| bibcode=1972ApJ...177..757W | doi-access=free }}

{{theories of gravitation}}
{{Relativity}}
{{Isaac Newton}}

[[Category:Theories of gravity]]
[[Category:Formalism (deductive)]]
[[Category:General relativity]]