at the time of the birth of general relativity ( gr ) , experimental confirmation was almost a side issue .
einstein did calculate observable effects of general relativity , such as the perihelion advance of mercury , which he knew to be an unsolved problem , and the deflection of light , which was subsequently verified , but compared to the inner consistency and elegance of the theory , he regarded such empirical questions as almost peripheral . but
today , experimental gravitation is a major component of the field , characterized by continuing efforts to test the theory s predictions , to search for gravitational imprints of high - energy particle interactions , and to detect gravitational waves from astronomical sources .
the modern history of experimental relativity can be divided roughly into four periods : genesis , hibernation , a golden era , and the quest for strong gravity .
the genesis ( 18871919 ) comprises the period of the two great experiments which were the foundation of relativistic physics  the michelson - morley experiment and the etvs experiment  and the two immediate confirmations of gr  the deflection of light and the perihelion advance of mercury .
following this was a period of hibernation ( 19201960 ) during which theoretical work temporarily outstripped technology and experimental possibilities , and , as a consequence , the field stagnated and was relegated to the backwaters of physics and astronomy . but beginning around 1960 , astronomical discoveries ( quasars , pulsars , cosmic background radiation ) and new experiments pushed gr to the forefront .
experimental gravitation experienced a golden era ( 19601980 ) during which a systematic , world - wide effort took place to understand the observable predictions of gr , to compare and contrast them with the predictions of alternative theories of gravity , and to perform new experiments to test them .
the period began with an experiment to confirm the gravitational frequency shift of light ( 1960 ) and ended with the reported decrease in the orbital period of the hulse - taylor binary pulsar at a rate consistent with the general relativity prediction of gravity wave energy loss ( 1979 ) . the results all supported gr , and most alternative theories of gravity fell by the wayside ( for a popular review , see ) .
since 1980 , the field has entered what might be termed a quest for strong gravity .
many of the remaining interesting weak - field predictions of the theory are extremely small and diffcult to check , in some cases requiring further technological development to bring them into detectable range .
the sense of a systematic assault on the weak - field predictions of gr has been supplanted to some extent by an opportunistic approach in which novel and unexpected ( and sometimes inexpensive ) tests of gravity have arisen from new theoretical ideas or experimental techniques , often from unlikely sources .
examples include the use of laser - cooled atom and ion traps to perform ultra - precise tests of special relativity ; the proposal of a  fifth  force , which led to a host of new tests of the weak equivalence principle ; and recent ideas of large extra dimensions , which have motived new tests of the inverse square law of gravity at sub - millimeter scales .
several major ongoing efforts also continue , principally the stanford gyroscope experiment , known as gravity probe - b . instead ,
much of the focus has shifted to experiments which can probe the effects of strong gravitational fields .
the principal figure of merit that distinguishes strong from weak gravity is the quantity gm / rc , where g is the newtonian gravitational constant , m is the characteristic mass scale of the phenomenon , r is the characteristic distance scale , and c is the speed of light . near the event horizon of a non - rotating black hole , or for the expanding observable universe , 0.5 ; for neutron stars , 0.2 .
these are the regimes of strong gravity . for the solar system <10 ; this is the regime of weak gravity . at one extreme
. will unification of the forces , or quantization of gravity at this scale leave observable effects accessible by experiment ?
dramatically improved tests of the equivalence principle or of the inverse square law are being designed , to search for or bound the imprinted effects of planck - scale phenomena . at the other extreme
are the strong fields associated with compact objects such as black holes or neutron stars .
astrophysical observations and gravitational wave detectors are being planned to explore and test gr in the strong - field , highly - dynamical regime associated with the formation and dynamics of these objects . in this living review
, we shall survey the theoretical frameworks for studying experimental gravitation , summarize the current status of experiments , and attempt to chart the future of the subject .
we shall not provide complete references to early work done in this field but instead will refer the reader to the appropriate review articles and monographs , specifically to theory and experiment in gravitational physics , hereafter referred to as tegp .
additional recent reviews in this subject are [ 142 , 145 , 150 , 139 , 37 , 117 ] . references to tegp will be by chapter or section , e.g.  tegp 8.9  .
the principle of equivalence has historically played an important role in the development of gravitation theory .
newton regarded this principle as such a cornerstone of mechanics that he devoted the opening paragraph of the principia to it . in 1907 , einstein used the principle as a basic element of general relativity .
we now regard the principle of equivalence as the foundation , not of newtonian gravity or of gr , but of the broader idea that spacetime is curved .
one elementary equivalence principle is the kind newton had in mind when he stated that the property of a body called  mass  is proportional to the  weight  , and is known as the weak equivalence principle ( wep ) .
an alternative statement of wep is that the trajectory of a freely falling body ( one not acted upon by such forces as electromagnetism and too small to be affected by tidal gravitational forces ) is independent of its internal structure and composition . in the simplest case of dropping two different bodies in a gravitational field ,
wep states that the bodies fall with the same acceleration ( this is often termed the universality of free fall , or uff ) .
a more powerful and far - reaching equivalence principle is known as the einstein equivalence principle ( eep ) .
it states that : 
 wep is valid.the outcome of any local non - gravitational experiment is independent of the velocity of the freely - falling reference frame in which it is performed.the outcome of any local non - gravitational experiment is independent of where and when in the universe it is performed . 
 
the outcome of any local non - gravitational experiment is independent of the velocity of the freely - falling reference frame in which it is performed .
the outcome of any local non - gravitational experiment is independent of where and when in the universe it is performed .
the second piece of eep is called local lorentz invariance ( lli ) , and the third piece is called local position invariance ( lpi ) . for example , a measurement of the electric force between two charged bodies is a local non - gravitational experiment ; a measurement of the gravitational force between two bodies ( cavendish experiment ) is not .
the einstein equivalence principle is the heart and soul of gravitational theory , for it is possible to argue convincingly that if eep is valid , then gravitation must be a  curved spacetime 
phenomenon , in other words , the effects of gravity must be equivalent to the effects of living in a curved spacetime . as a consequence of this argument , the only theories of gravity that can embody eep are those that satisfy the postulates of  metric theories of gravity  , which are : 
 spacetime is endowed with a symmetric metric.the trajectories of freely falling bodies are geodesics of that metric.in local freely falling reference frames , the non - gravitational laws of physics are those written in the language of special relativity . 
 
the trajectories of freely falling bodies are geodesics of that metric . in local freely falling reference frames ,
the non - gravitational laws of physics are those written in the language of special relativity .
the argument that leads to this conclusion simply notes that , if eep is valid , then in local freely falling frames , the laws governing experiments must be independent of the velocity of the frame ( local lorentz invariance ) , with constant values for the various atomic constants ( in order to be independent of location ) .
the only laws we know of that fulfill this are those that are compatible with special relativity , such as maxwell s equations of electromagnetism .
furthermore , in local freely falling frames , test bodies appear to be unaccelerated , in other words they move on straight lines ; but such  locally straight  lines simply correspond to  geodesics  in a curved spacetime ( tegp 2.3 ) .
general relativity is a metric theory of gravity , but then so are many others , including the brans - dicke theory .
the nonsymmetric gravitation theory ( ngt ) of moffat is not a metric theory . neither , in this narrow sense ,
2.3 ) , which , while based fundamentally on a spacetime metric , introduces additional fields ( dilatons , moduli ) that can couple to material stress - energy in a way that can lead to violations , say , of wep .
so the notion of curved spacetime is a very general and fundamental one , and therefore it is important to test the various aspects of the einstein equivalence principle thoroughly .
a direct test of wep is the comparison of the acceleration of two laboratory - sized bodies of different composition in an external gravitational field .
the simplest way to quantify such possible violations of wep in a form suitable for comparison with experiment is to suppose that for a body with inertial mass mi , the passive gravitational mass mp is no longer equal to mi , so that in a gravitational field g , the acceleration is given by mia = mpg .
now the inertial mass of a typical laboratory body is made up of several types of mass - energy : rest energy , electromagnetic energy , weak - interaction energy , and so on .
if one of these forms of energy contributes to mp differently than it does to mi , a violation of wep would result .
one could then write 
 ( 1)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 { m_{\rm{p } } } = { m_{\rm{i } } } + \sum\limits_a { { \eta ^a}{e^a } } /{c^2 } , 
 $ $ \end{document } 
 where e is the internal energy of the body generated by interaction a , and  is a dimensionless parameter that measures the strength of the violation of wep induced by that interaction , and c is the speed of light .
a measurement or limit on the fractional difference in acceleration between two bodies then yields a quantity called the  etvs ratio  given by 
 ( 2)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 \eta \equiv \frac{{2\left| { { a_1 } - { a_2 } } \right|}}{{\left| { { a_1 } + { a_2 } } \right| } } = \sum\limits_a { { \eta ^a}\left ( { \frac{{e_1^a}}{{{m_1}{c^2 } } } - \frac{{e_2^a}}{{{m_2}{c^2 } } } } \right ) , } 
 $ $ \end{document } 
 where we drop the subscript 
thus , experimental limits on  place limits on the wep - violation parameters . 
  many high - precision etvs - type experiments have been performed , from the pendulum experiments of newton , bessel and potter , to the classic torsionbalance measurements of etvs , dicke , braginsky   and their collaborators . in the modern torsion - balance experiments ,
two objects of different composition are connected by a rod or placed on a tray and suspended in a horizontal orientation by a fine wire .
if the gravitational acceleration of the bodies differs , there will be a torque induced on the suspension wire , related to the angle between the wire and the direction of the gravitational acceleration g. if the entire apparatus is rotated about some direction with angular velocity  , the torque will be modulated with period 2/. in the experiments of etvs and his collaborators , the wire and g were not quite parallel because of the centripetal acceleration on the apparatus due to the earth s rotation ; the apparatus was rotated about the direction of the wire . in the dicke and braginsky experiments ,
g was that of the sun , and the rotation of the earth provided the modulation of the torque at a period of 24 hr ( tegp 2.4 ( a ) ) . beginning in the late 1980s , numerous experiments were carried out primarily to search for a  fifth force  ( see sec .
2.3 ) , but their null results also constituted tests of wep . in the  free - fall galileo experiment  performed at the university of colorado , the relative free - fall acceleration of two bodies made of uranium and copper was measured using a laser interferometric technique .
et - wash  experiments carried out at the university of washington used a sophisticated torsion balance tray to compare the accelerations of various materials toward local topography on earth , movable laboratory masses , the sun and the galaxy [ 121 , 10 ] , and have recently reached levels of 410 .
the resulting upper limits on  are summarized in figure 1 ( tegp 14.1 ; for a bibliography of experiments , see ) . 
 figure 1 
 selected tests of the weak equivalence principle , showing bounds on  , which measures fractional difference in acceleration of different materials or bodies .
the free - fall and et - wash experiments were originally performed to search for a fifth force .
the blue band shows current bounds on  for gravitating bodies from lunar laser ranging ( lure ) . 

 

selected tests of the weak equivalence principle , showing bounds on  , which measures fractional difference in acceleration of different materials or bodies .
the free - fall and et - wash experiments were originally performed to search for a fifth force .
the blue band shows current bounds on  for gravitating bodies from lunar laser ranging ( lure ) . 
 
the second ingredient of eep , local lorentz invariance , has been tested to high - precision in the  mass anisotropy  experiments : the classic versions are the hughes - drever experiments , performed in the period 195960 independently by hughes and collaborators at yale university , and by drever at glasgow university ( tegp 2.4 ( b ) ) .
dramatically improved versions were carried out during the late 1980s using laser - cooled trapped atom techniques ( tegp 14.1 ) . a simple and useful way of interpreting these experiments is to suppose that the electromagnetic interactions suffer a slight violation of lorentz invariance , through a change in the speed of electromagnetic radiation c relative to the limiting speed of material test particles ( c0 , chosen to be unity via a choice of units ) , in other words , c1 ( see sec .
such a violation necessarily selects a preferred universal rest frame , presumably that of the cosmic background radiation , through which we are moving at about 300 km / s .
such a lorentz - non - invariant electromagnetic interaction would cause shifts in the energy levels of atoms and nuclei that depend on the orientation of the quantization axis of the state relative to our universal velocity vector , and on the quantum numbers of the state .
the presence or absence of such energy shifts can be examined by measuring the energy of one such state relative to another state that is either unaffected or is affected differently by the supposed violation . one way is to look for a shifting of the energy levels of states that are ordinarily equally spaced , such as the four j=3/2 ground states of the li nucleus in a magnetic field ( drever experiment ) ; another is to compare the levels of a complex nucleus with the atomic hyperfine levels of a hydrogen maser clock .
these experiments have all yielded extremely accurate results , quoted as limits on the parameter c1 in figure 2 . also included for comparison
is the corresponding limit obtained from michelson - morley type experiments ( for a review , see ) . 

figure 2 
 selected tests of local lorentz invariance showing the bounds on the parameter  , which measures the degree of violation of lorentz invariance in electromagnetism .
the michelson - morley , joos , and brillet - hall experiments test the isotropy of round - trip speed of light , the latter experiment using laser technology .
the centrifuge , two - photon absorption ( tpa ) and jpl experiments test the isotropy of light speed using one - way propagation . the remaining four experiments test isotropy of nuclear energy levels .
the limits assume a speed of earth of 300 km / s relative to the mean rest frame of the universe . 

 

selected tests of local lorentz invariance showing the bounds on the parameter  , which measures the degree of violation of lorentz invariance in electromagnetism .
the michelson - morley , joos , and brillet - hall experiments test the isotropy of round - trip speed of light , the latter experiment using laser technology .
the centrifuge , two - photon absorption ( tpa ) and jpl experiments test the isotropy of light speed using one - way propagation . the remaining four experiments test isotropy of nuclear energy levels .
the limits assume a speed of earth of 300 km / s relative to the mean rest frame of the universe . 
 
recent advances in atomic spectroscopy and atomic timekeeping have made it possible to test lli by checking the isotropy of the speed of light using oneway propagation ( as opposed to round - trip propagation , as in the michelson - morley experiment ) . in one experiment , for example ,
the relative phases of two hydrogen maser clocks at two stations of nasa s deep space tracking network were compared over five rotations of the earth by propagating a light signal one - way along an ultrastable fiberoptic link connecting them ( see sec .
although the bounds from these experiments are not as tight as those from mass - anisotropy experiments , they probe directly the fundamental postulates of special relativity , and thereby of lli ( tegp 14.1 , ) .
the principle of local position invariance , the third part of eep , can be tested by the gravitational redshift experiment , the first experimental test of gravitation proposed by einstein . despite the fact that einstein regarded this as a crucial test of gr , we now realize that it does not distinguish between gr and any other metric theory of gravity , but is only a test of eep .
a typical gravitational redshift experiment measures the frequency or wavelength shift z/=-/ between two identical frequency standards ( clocks ) placed at rest at different heights in a static gravitational field . if the frequency of a given type of atomic clock is the same when measured in a local , momentarily comoving freely falling frame ( lorentz frame ) , independent of the location or velocity of that frame , then the comparison of frequencies of two clocks at rest at different locations boils down to a comparison of the velocities of two local lorentz frames , one at rest with respect to one clock at the moment of emission of its signal , the other at rest with respect to the other clock at the moment of reception of the signal .
the frequency shift is then a consequence of the first - order doppler shift between the frames .
the result is a shift 
 ( 3)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 z = \delta u/{c^2 } , 
 $ $ \end{document } 
 where u is the difference in the newtonian gravitational potential between the receiver and the emitter .
if lpi is not valid , then it turns out that the shift can be written 
 ( 4)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 z = ( 1 + \alpha ) \delta u/{c^2 } , 
 $ $ \end{document } 
 where the parameter  may depend upon the nature of the clock whose shift is being measured ( see tegp 2.4 ( c )   for details ) . 
 
the first successful , high - precision redshift measurement was the series of pound - rebka - snider experiments of 19601965 that measured the frequency shift of gamma - ray photons from fe as they ascended or descended the jefferson physical laboratory tower at harvard university . the high accuracy achieved  one percent  was obtained by making use of the mssbauer effect to produce a narrow resonance line whose shift could be accurately determined . other experiments since 1960 measured the shift of spectral lines in the sun s gravitational field and the change in rate of atomic clocks transported aloft on aircraft , rockets and satellites .
figure 3 summarizes the important redshift experiments that have been performed since 1960 ( tegp 2.4 ( c ) ) . 

figure 3selected tests of local position invariance via gravitational redshift experiments , showing bounds on  , which measures degree of deviation of redshift from the formula /=u / c . 
  selected tests of local position invariance via gravitational redshift experiments , showing bounds on  , which measures degree of deviation of redshift from the formula /=u / c . the most precise standard redshift test to date was the vessot - levine rocket experiment that took place in june 1976 .
a hydrogen - maser clock was own on a rocket to an altitude of about 10,000 km and its frequency compared to a similar clock on the ground .
the experiment took advantage of the masers frequency stability by monitoring the frequency shift as a function of altitude . a sophisticated data acquisition scheme accurately eliminated all effects of the first - order doppler shift due to the rocket s motion , while tracking data were used to determine the payload s location and the velocity ( to evaluate the potential difference u , and the special relativistic time dilation )
redshift experiment performed in 1978 tested whether the relative rates of two different clocks depended upon position .
two hydrogen maser clocks and an ensemble of three superconducting - cavity stabilized oscillator ( scso ) clocks were compared over a 10-day period . during the period of the experiment ,
the solar potential u / c changed sinusoidally with a 24-hour period by 310 because of the earth s rotation , and changed linearly at 310 per day because the earth is 90 degrees from perihelion in april .
however , analysis of the data revealed no variations of either type within experimental errors , leading to a limit on the lpi violation parameter <210 .
this bound has been improved using more stable frequency standards [ 68 , 109 ] .
the varying gravitational redshift of earth - bound clocks relative to the highly stable millisecond pulsar psr 1937 + 21 , caused by the earth s motion in the solar gravitational field around the earth - moon center of mass ( amplitude 4000 km ) , has been measured to about 10 percent , and the redshift of stable oscillator clocks on the voyager spacecraft caused by saturn s gravitational field yielded a one percent test .
the solar gravitational redshift has been tested to about two percent using infrared oxygen triplet lines at the limb of the sun , and to one percent using oscillator clocks on the galileo spacecraft ( tegp 2.4 ( c )   and 14.1 ( a ) ) .
modern advances in navigation using earth - orbiting atomic clocks and accurate time - transfer must routinely take gravitational redshift and time - dilation effects into account .
for example , the global positioning system ( gps ) provides absolute accuracies of around 15 m ( even better in its military mode ) anywhere on earth , which corresponds to 50 nanoseconds in time accuracy at all times . yet
the difference in rate between satellite and ground clocks as a result of special and general relativistic effects is a whopping 39 microseconds per day ( 46 s from the gravitational redshift , and -7 s from time dilation ) .
if these effects were not accurately accounted for , gps would fail to function at its stated accuracy .
this represents a welcome practical application of gr ! ( for the role of gr in gps , see ; for a popular essay , see . )
if lpi is satisfied , the fundamental constants of non - gravitational physics should be constants in time . table 1 shows current bounds on cosmological variations in selected dimensionless constants . for discussion and references to early work , see tegp 2.4 ( c ) . 

table 1bounds on cosmological variation of fundamental constants of nongravitational physics . for references to earlier work , see tegp 2.4 ( c ) .constant klimit on k / k permethodhubble time 1.210 yrfine structure constant =e / c410h - maser vs. hg ion clock 910rb fountain vs. cs clock 610oklo natural reactor 61021-cm vs. molecular absorption at z=0.7 weak interaction
constant =gfmp2c / 1re , k decay rates0.1oklo natural reactor 0.06big bang nucleosynthesis [ 91 , 112]e - p mass ratio1mass shift in quasar spectra at z2proton g - factor ( gp)1021-cm vs. molecular absorption at z=0.7 
  bounds on cosmological variation of fundamental constants of nongravitational physics . for references to earlier work , see tegp 2.4 ( c ) . because the three parts of the einstein equivalence principle discussed above are so very different in their empirical consequences , it is tempting to regard them as independent theoretical principles . on the other hand ,
any complete and self - consistent gravitation theory must possess suffcient mathematical machinery to make predictions for the outcomes of experiments that test each principle , and because there are limits to the number of ways that gravitation can be meshed with the special relativistic laws of physics , one might not be surprised if there were theoretical connections between the three sub - principles .
for instance , the same mathematical formalism that produces equations describing the free fall of a hydrogen atom must also produce equations that determine the energy levels of hydrogen in a gravitational field , and thereby the ticking rate of a hydrogen maser clock .
hence a violation of eep in the fundamental machinery of a theory that manifests itself as a violation of wep might also be expected to show up as a violation of local position invariance . around 1960 ,
schiff conjectured that this kind of connection was a necessary feature of any self - consistent theory of gravity .
more precisely , schiff s conjecture states that any complete , self - consistent theory of gravity that embodies wep necessarily embodies eep . in other words , the validity of wep alone guarantees the validity of local lorentz and position invariance , and thereby of eep . if schiff s conjecture is correct , then etvs experiments may be seen as the direct empirical foundation for eep , hence for the interpretation of gravity as a curved - spacetime phenomenon .
of course , a rigorous proof of such a conjecture is impossible ( indeed , some special counter - examples are known ) , yet a number of powerful  plausibility 
the most general and elegant of these arguments is based upon the assumption of energy conservation .
this assumption allows one to perform very simple cyclic gedanken experiments in which the energy at the end of the cycle must equal that at the beginning of the cycle .
this approach was pioneered by dicke , nordtvedt and haugan ( see , e.g. ) .
a system in a quantum state a decays to state b , emitting a quantum of frequency . the quantum falls a height h in an external gravitational field and is shifted to frequency  , while the system in state b falls with acceleration gb . at the bottom , state a is rebuilt out of state b , the quantum of frequency  , and the kinetic energy mbgbh that state b has gained during its fall . the energy left over must be exactly enough , magah , to raise state a to its original location .
( here an assumption of local lorentz invariance permits the inertial masses ma and mb to be identified with the total energies of the bodies . ) if ga and gb depend on that portion of the internal energy of the states that was involved in the quantum transition from a to b according to 
 ( 5)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 \begin{array}{*{20}{c } } 
 { { g_a } = g(1 + \alpha { e_a}/{m_a}{c^2}),}&{{g_b } = g(1 + \alpha { e_b}/{m_b}{c^2}),}&{{e_a } - { e_b } \equiv h\nu } 
 \end{array } 
 $ $ \end{document } 
 ( violation of wep ) , then by conservation of energy , there must be a corresponding violation of lpi in the frequency shift of the form ( to lowest order in h/mc ) 
 ( 6)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 z = ( \nu ' - \nu ) /\nu '
u/{c^2}. 
 $ $ \end{document } 
 haugan generalized this approach to include violations of lli , ( tegp 2.5 ) . 
 

coordinate system and conventions : 
 x = t = time coordinate associated with the static nature of the static spherically symmetric ( sss ) gravitational field ; x=(x , y , z)=isotropic quasi - cartesian spatial coordinates ; spatial vector and gradient operations as in cartesian space . 

matter and field variables : 
 m0a = rest mass of particle a.ea=charge of particle a.xa(t)=world line of particle a.va=dxa/dt=coordinate velocity of particle a.a=electromagnetic vector potential ; e=a0a/t , b=a 
 gravitational potential : 
 u(x)arbitrary functions : 
 t(u ) , h(u ) , (u ) ,
(u ) ; eep is satisfied if ==(h / t ) for all u. 
 action : 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 i = - \sum\limits_a { { m_{0a}}\int { { { \left ( { t - hv_a^2 }
\right)}^{1/2}}dt + \sum\limits_a { { e_a } } \int { { a_\mu } \left ( { x_a^\nu } \right)v_a^\mu dt + { { \left ( { 8\pi } \right)}^ { - 1}}\int { \left ( { \epsilon { e^2 } - { \mu ^ { - 1}}{b^2 } } \right){d^4}x . } } } } 
 $ $ \end{document } 
 non - metric parameters : 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 \begin{array}{*{20}{l } } 
 { { \gamma _ 0 } = - c_0 ^ 2\left ( { \partial /\partial u } \right)\ln { { \left [ { \epsilon { { \left ( { t / h } \right)}^{1/2 } } } \right]}_0},}\\ 
 { { \lambda _ 0 } = - c_0 ^ 2\left ( { \partial /\partial u } \right)\ln { { \left [ { \mu { { \left ( { t / h } \right)}^{1/2 } } } \right]}_0},}\\ 
 { { \upsilon _ 0 } = 1 - { { \left ( { t{h^ { - 1 } } \epsilon \mu } \right)}_0 } , } 
 \end{array } 
 $ $ \end{document } 
 where c0=(t0/h0 ) and subscript  0  refers to a chosen point in space .
coordinate system and conventions : 
 x = t = time coordinate associated with the static nature of the static spherically symmetric ( sss ) gravitational field ; x=(x , y , z)=isotropic quasi - cartesian spatial coordinates ; spatial vector and gradient operations as in cartesian space . 

matter and field variables : 
 m0a = rest mass of particle a.ea=charge of particle a.xa(t)=world line of particle a.va=dxa/dt=coordinate velocity of particle a.a=electromagnetic vector potential ; e=a0a/t , b=a 
  m0a = rest mass of particle a. ea = charge of particle a. xa(t)=world line of particle a. va=dxa/dt = coordinate velocity of particle a. a=electromagnetic vector potential ; e=a0a/t , b=a gravitational potential : 
 u(x ) arbitrary functions : 
 t(u ) , h(u ) , (u ) , (u ) ; eep is satisfied if ==(h / t ) for all u. 
 action : 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 i = - \sum\limits_a { { m_{0a}}\int { { { \left ( { t - hv_a^2 } \right)}^{1/2}}dt + \sum\limits_a { { e_a } } \int { { a_\mu } \left ( { x_a^\nu } \right)v_a^\mu dt + { { \left ( { 8\pi } \right)}^ { - 1}}\int { \left ( { \epsilon { e^2 } - { \mu ^ { - 1}}{b^2 } } \right){d^4}x . } } } } 
 $ $ \end{document } 
  non - metric parameters : 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 \begin{array}{*{20}{l } } 
 { { \gamma _ 0 } = - c_0 ^ 2\left ( { \partial /\partial u } \right)\ln { { \left [ { \epsilon { { \left ( { t / h } \right)}^{1/2 } } } \right]}_0},}\\ 
 { { \lambda _ 0 } = - c_0 ^ 2\left ( { \partial /\partial u } \right)\ln { { \left [ { \mu { { \left ( { t / h } \right)}^{1/2 } } } \right]}_0},}\\ 
 { { \upsilon _ 0 } = 1 - { { \left ( { t{h^ { - 1 } } \epsilon \mu } \right)}_0 } , } 
 \end{array } 
 $ $ \end{document } 
 where c0=(t0/h0 ) and subscript  0  refers to a chosen point in space . if eep is satisfied , 0000 . 
 
the first successful attempt to prove schiff s conjecture more formally was made by lightman and lee .
they developed a framework called the th formalism that encompasses all metric theories of gravity and many non - metric theories ( box 1 ) .
it restricts attention to the behavior of charged particles ( electromagnetic interactions only ) in an external static spherically symmetric ( sss ) gravitational field , described by a potential u. it characterizes the motion of the charged particles in the external potential by two arbitrary functions t(u ) and h(u ) , and characterizes the response of electromagnetic fields to the external potential ( gravitationally modified maxwell equations ) by two functions (u ) and (u ) .
the forms of t , h ,  and  vary from theory to theory , but every metric theory satisfies 
 ( 7)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 \epsilon = \mu = { ( h / t)^{1/2 } } , 
 $ $ \end{document } 
 for all u. this consequence follows from the action of electrodynamics with a  minimal  or metric coupling : 
 ( 8)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 i = - \sum\limits_a { { m_{0a}}\int { { { \left ( { { g_{\mu \nu } } v_a^\mu v_a^\nu } \right)}^{1/2}}dt + \sum\limits_a { { e_a}\int { { a_\mu } \left ( { x_a^\nu } \right ) } } } } v_a^\mu dt - \frac{1}{{16\pi } } \int { \sqrt { - g } { g^{\mu \alpha } } { g^{\nu \beta } } { f_{\mu \nu } } { f_{\alpha \beta } } { d^4}x , } 
 $ $ \end{document } 
 where the variables are defined in box 1 , and where fa,a,. by identifying g00=t and gij = hij in a sss field , fi0=ei and fij=ijkbk , one obtains eq . ( 7 ) .
metric  form . in a given non - metric theory , the functions t , h
,  and  will depend in general on the full gravitational environment , including the potential of the earth , sun and galaxy , as well as on cosmological boundary conditions . which of these factors has the most influence on a given experiment will depend on the nature of the experiment . 
 
body made up of interacting charged particles , and found that the rate was independent of the internal electromagnetic structure of the body ( wep ) if and only if eq .
in other words , wepeep and schiff s conjecture was verified , at least within the restrictions built into the formalism .
certain combinations of the functions t , h ,  and  reflect different aspects of eep .
for instance , position or u - dependence of either of the combinations (t / h ) and (t / h ) signals violations of lpi , the first combination playing the role of the locally measured electric charge or fine structure constant .
the  non - metric parameters  0 and 0 ( box 1 ) are measures of such violations of eep . similarly , if the parameter 01(th)0 is non - zero anywhere , then violations of lli will occur .
this parameter is related to the difference between the speed of light c , and the limiting speed of material test particles c0 , given by 
 ( 9)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 \begin{array}{*{20}{c } } 
 { c = { { \left ( { { \epsilon _ 0}{\mu _ 0 } } \right)}^ { - 1/2}},}&{{c_0 } = { { \left ( { { t_0}/{h_0 } } \right)}^{1/2}}. } 
 \end{array } 
 $ $ \end{document } 
 in many applications , by suitable definition of units , c0 can be set equal to unity . if eep is valid , 000=0 everywhere . 
 
the rate of fall of a composite spherical test body of electromagnetically interacting particles then has the form 
 ( 10)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 { \rm{a } } = \left ( { { m_{\rm{p}}}/m } \right)\nabla u , 
 $ $ \end{document } 
 ( 11)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 { m_{\rm{p}}}/m = 1 + \left ( { e_{\rm{b}}^{{\rm{es}}}/mc_0 ^ 2 } \right)\left [ { 2{\gamma _ 0 } - \frac{8}{3}{\upsilon _ 0 } } \right ] + \left ( { e_{\rm{b}}^{{\rm{ms}}}/mc_0 ^ 2 } \right)\left [ { 2{\lambda _ 0 } - \frac{4}{3}{\upsilon _ 0 } } \right ] + \ldots , 
 $ $ \end{document } 
 where ebes and ebms are the electrostatic and magnetostatic binding energies of the body , given by 
 ( 12)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 e_{\rm{b}}^{{\rm{es } } } = - \frac{1}{4}t_0^{1/2}h_0^ { - 1 } \epsilon _ 0^ { - 1}\left\langle { \sum\limits_{ab } { { e_a}{e_b}r_{ab}^ { - 1 } } } \right\rangle , 
 $ $ \end{document } 
 ( 13)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 e_{\rm{b}}^{{\rm{ms } } } = - \frac{1}{8}t_0^{1/2}h_0^ { - 1}{\mu _ 0}\left\langle { \sum\limits_{ab } { { e_a}{e_b}r_{ab}^ { - 1}\left [ { { { \rm{v}}_a } \cdot { { \rm{v}}_b } + \left ( { { { \rm{v}}_a } \cdot { { \rm{n}}_{ab } } } \right)\left ( { { { \rm{v}}_b } \cdot { { \rm{n}}_{ab } } } \right ) } \right ] } } \right\rangle , 
 $ $ \end{document } 
 where rab=|xaxb| , nab=(xaxb)/rab , and the angle brackets denote an expectation value of the enclosed operator for the system s internal state
11 ) , and ultimately place limits on the non - metric parameters |0|<210 and |0|<310 .
( we set 0=0 because of very tight constraints on it from tests of lli . )
these limits are sufficiently tight to rule out a number of non - metric theories of gravity thought previously to be viable ( tegp 2.6 ( f ) ) . 
 
the th formalism also yields a gravitationally modified dirac equation that can be used to determine the gravitational redshift experienced by a variety of atomic clocks . for the redshift parameter  ( eq . ( 4 ) )
, the results are ( tegp 2.6 ( c ) ): 
 ( 14)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 \alpha = \left\ { { \begin{array}{*{20}{l } } 
 { - 3{\gamma _ 0 } + { \lambda _ 0}}&{{\rm{hydrogen\ hyperfine\ transition,\ h\ -\ maser\ clock}}}\\ 
 { - \frac{1}{2}\left ( { 3{\gamma _ 0 } + { \lambda _ 0 } } \right)}&{{\rm{electromagnetic\ mode\ in\ cavity,\ scso\ clock,}}}\\ 
 { - 2{\gamma _ 0}}&{{\rm{phonon\ mode\ in\ solid,\ principal\ transition\ in\ hydrogen}}{\rm{. } } } 
 \end{array } } \right . 

$ $ \end{document } 
  the redshift is the standard one ( =0 ) , independently of the nature of the clock if and only if 000 .
thus the vessot - levine rocket redshift experiment sets a limit on the parameter combination |300| ( figure 3 ) ; the null - redshift experiment comparing hydrogen - maser and scso clocks sets a limit on 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 \left| { { \alpha _ { \rm{h } } } - { \alpha _ { { \rm{scso } } } } } \right| = \frac{3}{2}\left| { { \gamma _ 0 } - { \lambda _ 0 } } \right| 
 $ $ \end{document } 
 . alvarez and mann [ 4 , 3 , 5 , 6 , 7 ] extended the th formalism to permit analysis of such effects as the lamb shift , anomalous magnetic moments and non - baryonic effects , and placed interesting bounds on eep violations . 
 
the th formalism can also be applied to tests of local lorentz invariance , but in this context it can be simplified . since most such tests do not concern themselves with the spatial variation of the functions t , h ,  , and  , but rather with observations made in moving frames , we can treat them as spatial constants .
then by rescaling the time and space coordinates , the charges and the electromagnetic fields , we can put the action in box 1 into the form ( tegp 2.6 ( a ) ) 
 ( 15)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 i = - { \sum\limits_a { { m_{0a}}\int { \left ( { 1 - v_a^2 } \right ) } } ^{1/2}}dt + \sum\limits_a { { e_a}\int { { a_\mu } \left ( { x_a^\nu } \right)v_a^\mu dt + { { \left ( { 8\pi } \right)}^ { - 1}}\int { \left ( { { e^2 } - { c^2}{b^2 } } \right){d^4}x , } } } 
 $ $ \end{document } 
 where ch0/t000=(10 )
. this amounts to using units in which the limiting speed c0 of massive test particles is unity , and the speed of light is c. if c1 , lli is violated ; furthermore , the form of the action above must be assumed to be valid only in some preferred universal rest frame .
the natural candidate for such a frame is the rest frame of the microwave background . 
  the electrodynamical equations which follow from eq .
( 15 ) yield the behavior of rods and clocks , just as in the full th formalism .
for example , the length of a rod moving through the rest frame in a direction parallel to its length will be observed by a rest observer to be contracted relative to an identical rod perpendicular to the motion by a factor 1v/2+o(v ) .
the energy and momentum of an electromagnetically bound body which moves with velocity v relative to the rest frame are given by 
 ( 16)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 \begin{array}{*{20}{l } } 
 { e = { m_{\rm{r } } } + { \tfrac{1}{2}}{m_{\rm{r}}}{v^2 } + { \tfrac{1}{2}}\delta m_{\rm{i}}^{ij}{v^i}{v^j},}\\ 
 { { p^i } = { m_{\rm{r}}}{v^i } + \delta m_{\rm{i}}^{ij}{v^j } , } 
 \end{array } 
 $ $ \end{document } 
 where mr = m0ebes , m0 is the sum of the particle rest masses , ebes is the electrostatic binding energy of the system ( eq .
( 12 ) with t01/2h00 - 1=1 ) , and 
 ( 17)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 \delta m_{\rm{i}}^{ij } = - 2\left ( { \frac{1}{{{c^2 } } } - 1 } \right)\left [ { \frac{4}{3}e_{\rm{b}}^{{\rm{es}}}{\delta ^{ij } } + \tilde e_{\rm{b}}^{{\rm{es}}ij } } \right ] , 
 $ $ \end{document } 
 where 
 ( 18)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 \tilde e_{\rm{b}}^{{\rm{es}}ij } = - \frac{1}{4}\left\langle { \sum\limits_{ab } { { e_a}{e_b}r_{ab}^ { - 1}\left ( { ( n_{ab}^in_{ab}^j - \frac{1}{3}{\delta ^{ij } } } \right ) } } \right\rangle . 

$ $ \end{document } 
 note that ( c1 ) corresponds to the parameter  plotted in figure 2 . 
 
( 15 ) can also be quantized , so that we may treat the interaction of photons with atoms via perturbation theory .
the energy of a photon is  times its frequency  , while its momentum is /c . using this approach
, one finds that the difference in round trip travel times of light along the two arms of the interferometer in the michelson - morley experiment is given by l0(v / c)(c1 ) .
the experimental null result then leads to the bound on ( c1 ) shown on figure 2 .
similarly the anisotropy in energy levels is clearly illustrated by the tensorial terms in eqs .
( 16 ) and ( 18 ) ; by evaluating ebesij for each nucleus in the various hughes - drever - type experiments and comparing with the experimental limits on energy differences , one obtains the extremely tight bounds also shown on figure 2 .
the behavior of moving atomic clocks can also be analysed in detail , and bounds on ( c1 ) can be placed using results from tests of time dilation and of the propagation of light . in some cases , it is advantageous to combine the c framework with a  kinematical  viewpoint that treats a general class of boost transformations between moving frames .
such kinematical approaches have been discussed by robertson , mansouri and sexl , and will ( see ) .
for example , in the  jpl  experiment , in which the phases of two hydrogen masers connected by a fiberoptic link were compared as a function of the earth s orientation , the predicted phase difference as a function of direction is , to first order in v , the velocity of the earth through the cosmic background , 
 ( 19)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 \delta \phi /\tilde \phi \approx - \frac{4}{3}\left ( { 1 - { c^2 } } \right)\left ( { { \rm{v } } \cdot { \rm{n } } - { \rm{v } } \cdot { { \rm{n}}_0 } } \right ) , 
 $ $ \end{document } 
 where =2l ,  is the maser frequency , l=21 km is the baseline , and where n and n0 are unit vectors along the direction of propagation of the light at a given time and at the initial time of the experiment , respectively .
the observed limit on a diurnal variation in the relative phase resulted in the bound |c1|<310 .
tighter bounds were obtained from a  two - photon absorption  ( tpa ) experiment , and a 1960s series of  mssbauer - rotor  experiments , which tested the isotropy of time dilation between a gamma ray emitter on the rim of a rotating disk and an absorber placed at the center . 
  in 1986 , as a result of a detailed reanalysis of etvs original data , fischbach et al .   suggested the existence of a fifth force of nature , with a strength of about a percent that of gravity , but with a range ( as defined by the range  of a yukawa potential , e / r ) of a few hundred meters .
this proposal dovetailed with earlier hints of a deviation from the inverse - square law of newtonian gravitation derived from measurements of the gravity profile down deep mines in australia , and with ideas from particle physics suggesting the possible presence of very low - mass particles with gravitational - strength couplings . during the next four years numerous experiments looked for evidence of the fifth force by searching for composition - dependent differences in acceleration , with variants of the etvs experiment or with free - fall galileo - type experiments .
although two early experiments reported positive evidence , the others all yielded null results . over the range between one and 104 meters , the null experiments produced upper limits on the strength of a postulated fifth force between 10 and 10 of the strength of gravity .
interpreted as tests of wep ( corresponding to the limit of infinite - range forces ) , the results of two representative experiments from this period , the free - fall galileo experiment and the early et - wash experiment , are shown in figure 1 . at the same time , tests of the inverse - square law of gravity were carried out by comparing variations in gravity measurements up tall towers or down mines or boreholes with gravity variations predicted using the inverse square law together with earth models and surface gravity data mathematically  continued  up the tower or down the hole . despite early reports of anomalies , independent tower ,
analyses of orbital data from planetary range measurements , lunar laser ranging , and laser tracking of the lageos satellite verified the inverse - square law to parts in 10 over scales of 10 to 10 km , and to parts in 10 over planetary scales of several astronomical units .
the consensus at present is that there is no credible experimental evidence for a fifth force of nature . for reviews and bibliographies , see [ 61 , 63 , 64 , 2 , 143 ]
nevertheless , theoretical evidence continues to mount that eep is likely to be violated at some level , whether by quantum gravity effects , by effects arising from string theory , or by hitherto undetected interactions , albeit at levels well below those that motivated the fifth - force searches . roughly speaking ,
in addition to the pure einsteinian gravitational interaction , which respects eep , theories such as string theory predict other interactions which do not . in string theory , for example , the existence of such eep - violating fields is assured , but the theory is not yet mature enough to enable calculation of their strength ( relative to gravity ) , or their range ( whether they are long range , like gravity , or short range , like the nuclear and weak interactions , and thus too short - range to be detectable ) .
in one simple example , one can write the lagrangian for the low - energy limit of string theory in the so - called  einstein frame  , in which the gravitational lagrangian is purely general relativistic : 
 ( 20)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\begin{array}{*{20}{c } } 
   { \tilde { \mathcal{l } } = \sqrt { - \tilde g } \left ( { { { \tilde g}^{\mu \nu } } \left [ { \frac{1}{{2k}}{{\tilde r}_{\mu \nu } } - \frac{1}{2}\tilde g(\varphi ) { \partial _ \mu } \varphi { \partial _ v}\varphi } \right ] } \right .
- u(\varphi ) { { \tilde g}^{\mu \nu } } { { \tilde g}^{\alpha \beta } } { f_{\mu \alpha } } { f_{\nu \beta } } } \\ 
   { \left .
{ + \bar \tilde \psi \left [ { i\tilde e_a^\mu { \gamma ^a}({\partial _ \mu } + { { \tilde \omega } _ \mu } + q{a_\mu } ) - \tilde m(\varphi ) } \right]\tilde \psi } \right),\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\ ; } 
 \end{array } 
 $ $ \end{document } 
 where g is the non - physical metric , r is the ricci tensor derived from it ,  is a dilaton field , and g , u and m are functions of . the lagrangian includes that for the electromagnetic field f , and that for particles , written in terms of dirac spinors . this is not a metric representation because of the coupling of  to matter via m( ) and u( ) .
a conformal transformation g=f()g , =f() , puts the lagrangian in the form (  jordan  frame ) 
 ( 21)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 { \mathcal l } = \sqrt { - g } \left ( { { g^{\mu \nu } } \left [ { \frac{1}{{2\kappa } } f\left ( \varphi \right){r_{\mu \nu } } - \frac{1}{2}f\left ( \varphi \right)\tilde g\left ( \varphi \right){\partial _ \mu } \varphi { \partial _ \nu } \varphi + \frac{3}{{4\kappa f\left ( \varphi \right)}}{\partial _ \mu } f{\partial _ \nu } f } \right ] - u\left ( \varphi \right){g^{\mu \nu } } { g^{\alpha \beta } } { f_{\mu \alpha } } { f_{\nu \beta } } + \overline \psi \left [ { ie_a^\mu { \gamma ^a}\left ( { { \partial _ \mu } + { \omega _ \mu } + q{a_\mu } } \right ) - \tilde m\left ( \varphi \right){f^{1/2 } } } \right]\psi } \right ) . 

$ $ \end{document } 
 one may choose f()=const./m( ) so that the particle lagrangian takes the metric form ( no explicit coupling to  ) , but the electromagnetic lagrangian will still couple non - metrically to u( ) .
the gravitational lagrangian here takes the form of a scalar - tensor theory ( sec .
but the non - metric electromagnetic term will , in general , produce violations of eep . for examples of specific models , see [ 125 , 50 ] . 
 
thus , eep and related tests are now viewed as ways to discover or place constraints on new physical interactions , or as a branch of  non - accelerator particle physics  , searching for the possible imprints of high - energy particle effects in the low - energy realm of gravity .
whether current or proposed experiments can actually probe these phenomena meaningfully is an open question at the moment , largely because of a dearth of firm theoretical predictions . despite this uncertainty , a number of experimental possibilities are being explored .
the project microscope , designed to test wep to 10 has been approved by the french space agency cnes for a possible 2004 launch .
another , known as satellite test of the equivalence principle ( step ) , is under consideration as a possible joint effort of nasa and the european space agency ( esa ) , with the goal of a 10 test .
the gravitational redshift could be improved to the 10 level using atomic clocks on board a spacecraft which would travel to within four solar radii of the sun .
laboratory tests of the gravitational inverse square law at sub - millimeter scales are being developed as ways to search for new short - range interactions or for the existence of large extra dimensions ; the challenge of these experiments is to distinguish gravitation - like interactions from electromagnetic and quantum mechanical ( casimir ) effects .
the overwhelming empirical evidence supporting the einstein equivalence principle , discussed in the previous section , supports the conclusion that the only theories of gravity that have a hope of being viable are metric theories , or possibly theories that are metric apart from possible weak or short - range non - metric couplings ( as in string theory ) .
therefore for the remainder of this article , we shall turn our attention exclusively to metric theories of gravity , which assume that ( i ) there exists a symmetric metric , ( ii ) test bodies follow geodesics of the metric , and ( iii ) in local lorentz frames , the non - gravitational laws of physics are those of special relativity .
the property that all non - gravitational fields should couple in the same manner to a single gravitational field is sometimes called  universal coupling  .
because of it , one can discuss the metric as a property of spacetime itself rather than as a field over spacetime .
this is because its properties may be measured and studied using a variety of different experimental devices , composed of different non - gravitational fields and particles , and , because of universal coupling , the results will be independent of the device .
thus , for instance , the proper time between two events is a characteristic of spacetime and of the location of the events , not of the clocks used to measure it .
consequently , if eep is valid , the non - gravitational laws of physics may be formulated by taking their special relativistic forms in terms of the minkowski metric  and simply  going over  to new forms in terms of the curved spacetime metric g , using the mathematics of differential geometry .
the details of this  going over  can be found in standard textbooks ( [ 94 , 136 ] , tegp 3.2 . ) . in any metric theory of gravity , matter and
non - gravitational fields respond only to the spacetime metric g. in principle , however , there could exist other gravitational fields besides the metric , such as scalar fields , vector fields , and so on . if , by our strict definition of metric theory , matter does not couple to these fields , what can their role in gravitation theory be ?
their role must be that of mediating the manner in which matter and non - gravitational fields generate gravitational fields and produce the metric ; once determined , however , the metric alone acts back on the matter in the manner prescribed by eep .
what distinguishes one metric theory from another , therefore , is the number and kind of gravitational fields it contains in addition to the metric , and the equations that determine the structure and evolution of these fields . from this viewpoint
, one can divide all metric theories of gravity into two fundamental classes :  purely dynamical  and  prior - geometric  . by  purely dynamical metric theory
 we mean any metric theory whose gravitational fields have their structure and evolution determined by coupled partial differential field equations . in other words , the behavior of each field is influenced to some extent by a coupling to at least one of the other fields in the theory . by 
absolute elements  , fields or equations whose structure and evolution are given a priori , and are independent of the structure and evolution of the other fields of the theory .
these  absolute elements  typically include flat background metrics  , cosmic time coordinates t , and algebraic relationships among otherwise dynamical fields , such as g=h+kk , where h and k may be dynamical fields .
general relativity is a purely dynamical theory since it contains only one gravitational field , the metric itself , and its structure and evolution are governed by partial differential equations ( einstein s equations ) .
brans - dicke theory and its generalizations are purely dynamical theories ; the field equation for the metric involves the scalar field ( as well as the matter as source ) , and that for the scalar field involves the metric .
rosen s bimetric theory is a prior - geometric theory : it has a non - dynamical , riemann - flat background metric  , and the field equations for the physical metric g involve . by discussing metric theories of gravity from this broad point of view , it is possible to draw some general conclusions about the nature of gravity in different metric theories , conclusions that are reminiscent of the einstein equivalence principle , but that are subsumed under the name  strong equivalence principle  . consider a local , freely falling frame in any metric theory of gravity .
let this frame be small enough that inhomogeneities in the external gravitational fields can be neglected throughout its volume .
on the other hand , let the frame be large enough to encompass a system of gravitating matter and its associated gravitational fields . the system could be a star , a black hole , the solar system or a cavendish experiment . call this frame a 
quasi - local lorentz frame  . to determine the behavior of the system we must calculate the metric .
first we determine the external behavior of the metric and gravitational fields , thereby establishing boundary values for the fields generated by the local system , at a boundary of the quasi - local frame  far  from the local system .
but because the metric is coupled directly or indirectly to the other fields of the theory , its structure and evolution will be influenced by those fields , and in particular by the boundary values taken on by those fields far from the local system .
this will be true even if we work in a coordinate system in which the asymptotic form of g in the boundary region between the local system and the external world is that of the minkowski metric .
thus the gravitational environment in which the local gravitating system resides can influence the metric generated by the local system via the boundary values of the auxiliary fields .
consequently , the results of local gravitational experiments may depend on the location and velocity of the frame relative to the external environment .
of course , local non - gravitational experiments are unaffected since the gravitational fields they generate are assumed to be negligible , and since those experiments couple only to the metric , whose form can always be made locally minkowskian at a given spacetime event .
local gravitational experiments might include cavendish experiments , measurement of the acceleration of massive self - gravitating bodies , studies of the structure of stars and planets , or analyses of the periods of  gravitational clocks  .
( i)a theory which contains only the metric g yields local gravitational physics which is independent of the location and velocity of the local system .
this follows from the fact that the only field coupling the local system to the environment is g , and it is always possible to find a coordinate system in which g takes the minkowski form at the boundary between the local system and the external environment . thus the asymptotic values of g are constants independent of location , and are asymptotically lorentz invariant , thus independent of velocity .
general relativity is an example of such a theory.(ii)a theory which contains the metric g and dynamical scalar fields a yields local gravitational physics which may depend on the location of the frame but which is independent of the velocity of the frame .
this follows from the asymptotic lorentz invariance of the minkowski metric and of the scalar fields , but now the asymptotic values of the scalar fields may depend on the location of the frame .
an example is brans - dicke theory , where the asymptotic scalar field determines the effective value of the gravitational constant , which can thus vary as  varies . on the other hand , a form of velocity dependence in local physics
can enter indirectly if the asymptotic values of the scalar field vary with time cosmologically .
then the rate of variation of the gravitational constant could depend on the velocity of the frame.(iii)a theory which contains the metric g and additional dynamical vector or tensor fields or prior - geometric fields yields local gravitational physics which may have both location and velocity - dependent effects . 
  a theory which contains only the metric g yields local gravitational physics which is independent of the location and velocity of the local system .
this follows from the fact that the only field coupling the local system to the environment is g , and it is always possible to find a coordinate system in which g takes the minkowski form at the boundary between the local system and the external environment .
thus the asymptotic values of g are constants independent of location , and are asymptotically lorentz invariant , thus independent of velocity .
a theory which contains the metric g and dynamical scalar fields a yields local gravitational physics which may depend on the location of the frame but which is independent of the velocity of the frame .
this follows from the asymptotic lorentz invariance of the minkowski metric and of the scalar fields , but now the asymptotic values of the scalar fields may depend on the location of the frame .
an example is brans - dicke theory , where the asymptotic scalar field determines the effective value of the gravitational constant , which can thus vary as  varies .
on the other hand , a form of velocity dependence in local physics can enter indirectly if the asymptotic values of the scalar field vary with time cosmologically .
then the rate of variation of the gravitational constant could depend on the velocity of the frame . a theory which contains the metric g and additional dynamical vector or tensor fields or prior - geometric fields yields local gravitational physics which may have both location and velocity - dependent effects .
these ideas can be summarized in the strong equivalence principle ( sep ) , which states that : 
 wep is valid for self - gravitating bodies as well as for test bodies.the outcome of any local test experiment is independent of the velocity of the ( freely falling ) apparatus.the outcome of any local test experiment is independent of where and when in the universe it is performed . 

the distinction between sep and eep is the inclusion of bodies with self - gravitational interactions ( planets , stars ) and of experiments involving gravitational forces ( cavendish experiments , gravimeter measurements ) .
note that sep contains eep as the special case in which local gravitational forces are ignored . 
 
the outcome of any local test experiment is independent of the velocity of the ( freely falling ) apparatus .
the outcome of any local test experiment is independent of where and when in the universe it is performed .
the above discussion of the coupling of auxiliary fields to local gravitating systems indicates that if sep is strictly valid , there must be one and only one gravitational field in the universe , the metric g. these arguments are only suggestive however , and no rigorous proof of this statement is available at present .
empirically it has been found that every metric theory other than gr introduces auxiliary gravitational fields , either dynamical or prior geometric , and thus predicts violations of sep at some level ( here we ignore quantum - theory inspired modifications to gr involving  r 
3.6 , we shall discuss experimental evidence for the validity of sep . despite the possible existence of long - range gravitational fields in addition to the metric in various metric theories of gravity
, the postulates of those theories demand that matter and non - gravitational fields be completely oblivious to them .
the only gravitational field that enters the equations of motion is the metric g. the role of the other fields that a theory may contain can only be that of helping to generate the spacetime curvature associated with the metric .
matter may create these fields , and they plus the matter may generate the metric , but they can not act back directly on the matter .
thus the metric and the equations of motion for matter become the primary entities for calculating observable effects , and all that distinguishes one metric theory from another is the particular way in which matter and possibly other gravitational fields generate the metric .
the comparison of metric theories of gravity with each other and with experiment becomes particularly simple when one takes the slow - motion , weak - field limit .
this approximation , known as the post - newtonian limit , is sufficiently accurate to encompass most solar - system tests that can be performed in the foreseeable future .
it turns out that , in this limit , the spacetime metric g predicted by nearly every metric theory of gravity has the same structure .
it can be written as an expansion about the minkowski metric ( =diag(-1 , 1 , 1 , 1 ) ) in terms of dimensionless gravitational potentials of varying degrees of smallness .
these potentials are constructed from the matter variables ( box 2 ) in imitation of the newtonian gravitational potential 
 ( 22)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 u\left ( { { \rm{x}},t } \right ) \equiv \int { \rho \left ( { { \rm{x'}},t } \right){{\left| { { \rm{x } } - { \rm{x ' } } } \right|}^ { - 1}}{d^3}x ' . } 

$ $ \end{document } 
 the  order of smallness  is determined according to the rules \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 u \sim { v^2 } \sim \pi \sim p/\rho \sim \in , { v^i } \sim \left| { d / dt } \right|/\left| { d / dx } \right| \sim { \in ^{1/2 } } 
 $ $ \end{document } , and so on ( we use units
a consistent post - newtonian limit requires determination of g00 correct through o( ) , g0i through o( ) and gij through o( ) ( for details see tegp 4.1 ) .
the only way that one metric theory differs from another is in the numerical values of the coefficients that appear in front of the metric potentials .
the parametrized post - newtonian ( ppn ) formalism inserts parameters in place of these coefficients , parameters whose values depend on the theory under study . in the current version of the ppn formalism ,
summarized in box 2 , ten parameters are used , chosen in such a manner that they measure or indicate general properties of metric theories of gravity ( table 2 ) . under reasonable assumptions about the kinds of potentials that can be present at post - newtonian order ( basically only poisson - like potentials ) , one finds that ten ppn parameters exhaust the possibilities . 

table 2the ppn parameters and their significance ( note that 3 
 has been shown twice to indicate that it is a measure of two effects).pameterwhat it means relative to grvalue in grvalue in semi - conservative theoriesvalue in fully - conservative theories 
  
 how much space - curvature produced by unit rest mass?1 
  

  

  
 how much  non - linearity  in the superposition law for gravity?1 
  

  

  
 preferred - location effects?0 
  

  

  
 1 
 preferred - frame effects?0 
  
 1 
 0 
  
 2 
 0 
  
 2 
 0 
  
 3 
 000 
  
 3 
 violation of conservation of total momentum?000 
  
 1 
 000 
  
 2 
 000 
  
 3 
 000 
  
 4 
 000 
  the ppn parameters and their significance ( note that 3 
 has been shown twice to indicate that it is a measure of two effects ) . the parameters  and  are the usual eddington - robertson - schiff parameters used to describe the  classical  tests of gr , and are in some sense the most important ; they are the only non - zero parameters in gr and scalar - tensor gravity .
the parameter  is non - zero in any theory of gravity that predicts preferred - location effects such as a galaxy - induced anisotropy in the local gravitational constant gl ( also called 
effects ) ; 1 , 2 , 3 measure whether or not the theory predicts post - newtonian preferred - frame effects ; 3 , 1 , 2 , 3 , 4 measure whether or not the theory predicts violations of global conservation laws for total momentum .
next to  and  , the parameters 1 and 2 occur most frequently with non - trivial null values . in table 2
we show the values these parameters take ( i ) in gr , ( ii ) in any theory of gravity that possesses conservation laws for total momentum , called  semi - conservative  ( any theory that is based on an invariant action principle is semi - conservative ) , and ( iii ) in any theory that in addition possesses six global conservation laws for angular momentum , called  fully conservative  ( such theories automatically predict no post - newtonian preferred - frame effects ) .
semi - conservative theories have five free ppn parameters (  ,  ,  , 1 , 2 ) while fully conservative theories have three (  ,  ,  ) .
the ppn formalism was pioneered by kenneth nordtvedt , who studied the post - newtonian metric of a system of gravitating point masses , extending earlier work by eddington , robertson and schiff ( tegp 4.2 ) .
a general and unified version of the ppn formalism was developed by will and nordtvedt .
the canonical version , with conventions altered to be more in accord with standard textbooks such as , is discussed in detail in tegp 4 .
other versions of the ppn formalism have been developed to deal with point masses with charge , fluid with anisotropic stresses , bodies with strong internal gravity , and post - post - newtonian effects ( tegp 4.2 , 14.2 ) . 

coordinate system : the framework uses a nearly globally lorentz coordinate system in which the coordinates are ( t , x , x , x ) .
all coordinate arbitrariness (  gauge freedom  ) has been removed by specialization of the coordinates to the standard ppn gauge ( tegp 4.2 ) .
units are chosen so that g = c=1 , where g is the physically measured newtonian constant far from the solar system . 

matter variables : 
  : density of rest mass as measured in a local freely falling frame momentarily comoving with the gravitating matter;v=(dxi / dt ) : coordinate velocity of the matter;w : coordinate velocity of the ppn coordinate system relative to the mean rest - frame of the universe;p : pressure as measured in a local freely falling frame momentarily comoving with the matter; internal energy per unit rest mass .
it includes all forms of non - rest - mass , non - gravitational energy , e.g. energy of compression and thermal energy . 

ppn parameters :  ,  ,  , 1 , 2 , 3 , 1 , 2 , 3 , 4 . 

\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 \begin{array}{*{20}{l } } 
 { { g_{00 } } = - 1 + 2u - 2\beta { u^2 } - 2\xi { \phi_w } + \left ( { 2\gamma + 2 + { \alpha _ 3 } + { \zeta _ 1 } - 2\xi } \right){\phi _ 1 } + 2\left ( { 3\gamma - 2\beta + 1 + { \zeta _ 2 } + \xi } \right){\phi _ 2 } + 2\left ( { 1 + { \zeta _ 3 } } \right){\phi _ 3 } + 2\left ( { 3\gamma + 3{\zeta _ 4 } - 2\xi } \right){\phi _ 4 } - \left ( { { \zeta _ 1 } - 2\xi } \right){\mathcal a } - \left ( { { \alpha _ 1 } - { \alpha _ 2 } - { \alpha _ 3 } } \right){w^2}u - { \alpha _ 2}{w^i}{w^j}{u_{ij } } + \left ( { 2{\alpha _ 3 } - { \alpha _ 1 } } \right){w^i}{v_i } + o\left ( { { \epsilon ^3 } } \right),}\\ 
 { { g_{0i } } = - \frac{1}{2}\left ( { 4\gamma + 3 + { \alpha _ 1 } - { \alpha _ 2 } + { \zeta _ 1 } - 2\xi } \right){v_i } - \frac{1}{2}\left ( { 1 + { \alpha _ 2 } - { \zeta _ 1 } + 2\xi } \right){w_i } - \frac{1}{2}\left ( { { \alpha _ 1 } - 2{\alpha _ 2 } } \right){w^i}u - { \alpha _ 2}{w^j}{u_{ij } } + o\left ( { { \epsilon ^{5/2 } } } \right),}\\ 
 { { g_{ij } } = \left ( { 1 + 2\gamma u + o\left ( { { \epsilon ^2 } } \right ) } \right){\delta _ { ij . } } } 
 \end{array } 
 $ $ \end{document } 

 metric potentials : 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 \begin{array}{*{20}{c } } 
 { u = \int { \frac{{\rho ' } } { { \left| { { \rm{x } } - { \rm{x ' } } } \right|}}{d^3}x ' , } } \\ 
 { { u_{ij } } = \int { \frac{{\rho ' { { \left ( { x - x ' } \right)}_i}{{\left ( { x - x ' } \right)}_j}}}{{{{\left| { { \rm{x } } - { \rm{x ' } } } \right|}^3}}}{d^3}x ' , } } \\ 
 { { \phi _ w } = \int { \frac{{\rho ' \rho ' ' \left ( { { \rm{x } } - { \rm{x ' } } } \right)}}{{{{\left| { { \rm{x } } - { \rm{x ' } } } \right|}^3 } } } \cdot \left ( { \frac{{{\rm{x ' } } - { \rm{x''}}}}{{\left| { { \rm{x } } - { \rm{x '' } } } \right| } } - \frac{{{\rm{x } } - { \rm{x''}}}}{{\left| { { \rm{x ' } } - { \rm{x '' } } } \right| } } } \right){d^3}x'{d^3}x '' , } } \\ 
 { { \mathcal a } = \int { \frac{{\rho ' { { \left [ { { \rm{v ' } } \cdot \left ( { { \rm{x } } - { \rm{x ' } } } \right ) } \right]}^2}}}{{{{\left| { { \rm{x } } - { \rm{x ' } } } \right|}^3}}}{d^3}x ' , } } \\ 
 { { \phi _ 1 } = \int { \frac{{\rho ' { { v'}^2}}}{{\left| { { \rm{x } } - { \rm{x ' } } } \right|}}{d^3}x ' , } } \\ 
 { { \phi _ 2 } = \int { \frac{{\rho ' u'}}{{\left| { { \rm{x } } - { \rm{x ' } } } \right|}}{d^3}x ' , } } \\ 
 { { \phi _ 3 } = \int { \frac{{\rho ' \pi ' } } { { \left| { { \rm{x } } - { \rm{x ' } } } \right|}}{d^3}x ' , } } \\ 
 { { \phi _ 4 } = \int { \frac{{p'}}{{\left| { { \rm{x } } - { \rm{x ' } } } \right|}}{d^3}x ' , } } \\ 
 { { v_i } = \int { \frac{{\rho ' { { v'}_i}}}{{\left| { { \rm{x } } - { \rm{x ' } } } \right|}}{d^3}x ' , } } \\ 
 { { w_i } = \int { \frac{{\rho ' \left [ { { \rm{v ' } } \cdot \left ( { { \rm{x } } - { \rm{x ' } } } \right ) } \right]{{\left ( { x - x ' } \right)}_i}}}{{{{\left| { { \rm{x } } - { \rm{x ' } } } \right|}^3}}}{d^3}x ' . } } 
 \end{array } 
 $ $ \end{document } 
 stress - energy tensor ( perfect fluid ) : 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 \begin{array}{*{20}{l } } 
 { { t^{00 } } = \rho \left ( { 1 + \pi + { v^2 } + 2u } \right),}\\ 
 { { t^{0i } } = \rho { v^i}\left ( { 1 + \pi + { v^2 } + 2u + p/\rho } \right),}\\ 
 { { t^{ij } } = \rho { v^i}{v^j}\left ( { 1 + \pi + { v^2 } + 2u + p/\rho } \right ) + p{\delta ^{ij}}\left ( { 1 - 2\gamma u } \right ) . } 
 \end{array } 
 $ $ \end{document } 

 stressed matter : \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 { t^{\mu \nu } } _ { ; \nu } = 0 , 

$ $ \end{document}test bodies : \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 { d^2}{x^\mu } /d{\lambda ^2 } + { \gamma ^\mu } _ { \nu \lambda } \left ( { d{x^\nu } /d\lambda } \right)\left ( { d{x^\lambda } /d\lambda }
\right ) = 0 , 
 $ $ \end{document}maxwell s equations : \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 \begin{array}{*{20}{c } } 
 { { f^{\mu \nu } } _ { ; \nu } = 4\pi { j^\mu } , } & { { f_{\mu \nu } } = { a_{\nu ; \mu } } - { a_{\mu ; \nu } } . } 

\end{array } 
 $ $ \end{document } 

  coordinate system : the framework uses a nearly globally lorentz coordinate system in which the coordinates are ( t , x , x , x ) .
all coordinate arbitrariness (  gauge freedom  ) has been removed by specialization of the coordinates to the standard ppn gauge ( tegp 4.2 ) .
units are chosen so that g = c=1 , where g is the physically measured newtonian constant far from the solar system . 

matter variables : 
  : density of rest mass as measured in a local freely falling frame momentarily comoving with the gravitating matter;v=(dxi / dt ) : coordinate velocity of the matter;w : coordinate velocity of the ppn coordinate system relative to the mean rest - frame of the universe;p : pressure as measured in a local freely falling frame momentarily comoving with the matter; internal energy per unit rest mass .
it includes all forms of non - rest - mass , non - gravitational energy , e.g. energy of compression and thermal energy . 
 
 : density of rest mass as measured in a local freely falling frame momentarily comoving with the gravitating matter ; v=(dxi / dt ) : coordinate velocity of the matter ; w : coordinate velocity of the ppn coordinate system relative to the mean rest - frame of the universe ; p : pressure as measured in a local freely falling frame momentarily comoving with the matter ;  internal energy per unit rest mass .
it includes all forms of non - rest - mass , non - gravitational energy , e.g. energy of compression and thermal energy .
ppn parameters :  ,  ,  , 1 , 2 , 3 , 1 , 2 , 3 , 4 . 

\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 \begin{array}{*{20}{l } } 
 { { g_{00 } } = - 1 + 2u - 2\beta { u^2 } - 2\xi { \phi_w } + \left ( { 2\gamma + 2 + { \alpha _ 3 } + { \zeta _ 1 } - 2\xi } \right){\phi _ 1 } + 2\left ( { 3\gamma - 2\beta + 1 + { \zeta _ 2 } + \xi } \right){\phi _ 2 } + 2\left ( { 1 + { \zeta _ 3 } } \right){\phi _ 3 } + 2\left ( { 3\gamma + 3{\zeta _ 4 } - 2\xi } \right){\phi _ 4 } - \left ( { { \zeta _ 1 } - 2\xi } \right){\mathcal a } - \left ( { { \alpha _ 1 } - { \alpha _ 2 } - { \alpha _ 3 } } \right){w^2}u - { \alpha _ 2}{w^i}{w^j}{u_{ij } } + \left ( { 2{\alpha _ 3 } - { \alpha _ 1 } } \right){w^i}{v_i } + o\left ( { { \epsilon ^3 } } \right),}\\ 
 { { g_{0i } } = - \frac{1}{2}\left ( { 4\gamma + 3 + { \alpha _ 1 } - { \alpha _ 2 } + { \zeta _ 1 } - 2\xi } \right){v_i } - \frac{1}{2}\left ( { 1 + { \alpha _ 2 } - { \zeta _ 1 } + 2\xi } \right){w_i } - \frac{1}{2}\left ( { { \alpha _ 1 } - 2{\alpha _ 2 } } \right){w^i}u - { \alpha _ 2}{w^j}{u_{ij } } + o\left ( { { \epsilon ^{5/2 } } } \right),}\\ 
 { { g_{ij } } = \left ( { 1 + 2\gamma u + o\left ( { { \epsilon ^2 } } \right ) } \right){\delta _ { ij . } } } 
 \end{array } 
 $ $ \end{document } 
 
 metric potentials : 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 \begin{array}{*{20}{c } } 
 { u = \int { \frac{{\rho ' } } { { \left| { { \rm{x } } - { \rm{x ' } } } \right|}}{d^3}x ' , } } \\ 
 { { u_{ij } } = \int { \frac{{\rho ' { { \left ( { x - x ' } \right)}_i}{{\left ( { x - x ' } \right)}_j}}}{{{{\left| { { \rm{x } } - { \rm{x ' } } } \right|}^3}}}{d^3}x ' , } } \\ 
 { { \phi _ w } = \int { \frac{{\rho ' \rho ' ' \left ( { { \rm{x } } - { \rm{x ' } } } \right)}}{{{{\left| { { \rm{x } } - { \rm{x ' } } } \right|}^3 } } } \cdot \left ( { \frac{{{\rm{x ' } } - { \rm{x''}}}}{{\left| { { \rm{x } } - { \rm{x '' } } } \right| } } - \frac{{{\rm{x } } - { \rm{x''}}}}{{\left| { { \rm{x ' } } - { \rm{x '' } } } \right| } } } \right){d^3}x'{d^3}x '' , } } \\ 
 { { \mathcal a } = \int { \frac{{\rho ' { { \left [ { { \rm{v ' } } \cdot \left ( { { \rm{x } } - { \rm{x ' } } } \right ) } \right]}^2}}}{{{{\left| { { \rm{x } } - { \rm{x ' } } } \right|}^3}}}{d^3}x ' , } } \\ 
 { { \phi _ 1 } = \int { \frac{{\rho ' { { v'}^2}}}{{\left| { { \rm{x } } - { \rm{x ' } } } \right|}}{d^3}x ' , } } \\ 
 { { \phi _ 2 } = \int { \frac{{\rho ' u'}}{{\left| { { \rm{x } } - { \rm{x ' } } } \right|}}{d^3}x ' , } } \\ 

{ { \phi _ 3 } = \int { \frac{{\rho ' \pi ' } } { { \left| { { \rm{x } } - { \rm{x ' } } } \right|}}{d^3}x ' , } } \\ 
 { { \phi _ 4 } = \int { \frac{{p'}}{{\left| { { \rm{x } } - { \rm{x ' } } } \right|}}{d^3}x ' , } } \\ 
 { { v_i } = \int { \frac{{\rho ' { { v'}_i}}}{{\left| { { \rm{x } } - { \rm{x ' } } } \right|}}{d^3}x ' , } } \\ 
 { { w_i } = \int { \frac{{\rho ' \left [ { { \rm{v ' } } \cdot \left ( { { \rm{x } } - { \rm{x ' } } } \right ) } \right]{{\left ( { x - x ' } \right)}_i}}}{{{{\left| { { \rm{x } } - { \rm{x ' } } } \right|}^3}}}{d^3}x ' . } } 
 \end{array } 
 $ $ \end{document } 
  stress - energy tensor ( perfect fluid ) : 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 \begin{array}{*{20}{l } } 
 { { t^{00 } } = \rho \left ( { 1 + \pi + { v^2 } + 2u } \right),}\\ 
 { { t^{0i } } = \rho { v^i}\left ( { 1 + \pi + { v^2 } + 2u + p/\rho } \right),}\\ 
 { { t^{ij } } = \rho { v^i}{v^j}\left ( { 1 + \pi + { v^2 } + 2u + p/\rho } \right ) +
\end{array } 
 $ $ \end{document } 
 
 stressed matter : \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 { t^{\mu \nu } } _ { ; \nu } = 0 , 
 $ $ \end{document}test bodies : \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 { d^2}{x^\mu } /d{\lambda ^2 } + { \gamma ^\mu } _ { \nu \lambda } \left ( { d{x^\nu } /d\lambda } \right)\left ( { d{x^\lambda } /d\lambda } \right ) = 0 , 
 $ $ \end{document}maxwell s equations : \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 \begin{array}{*{20}{c } } 
 { { f^{\mu \nu } } _ { ; \nu } = 4\pi { j^\mu } , } & { { f_{\mu \nu } } = { a_{\nu ; \mu } } - { a_{\mu ; \nu } } . } 

\end{array } 
 $ $ \end{document } 
  stressed matter : \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 { t^{\mu \nu } } _ { ; \nu } = 0 , 
 $ $ \end{document } test bodies : \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 { d^2}{x^\mu } /d{\lambda
^2 } + { \gamma ^\mu } _ { \nu \lambda } \left ( { d{x^\nu } /d\lambda } \right)\left ( { d{x^\lambda } /d\lambda }
\right ) = 0 , 
 $ $ \end{document } maxwell s equations : \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 \begin{array}{*{20}{c } } 
 { { f^{\mu \nu } } _ { ; \nu } = 4\pi { j^\mu } , } & { { f_{\mu \nu } } = { a_{\nu ; \mu } } - { a_{\mu ; \nu } } . } 

\end{array } 
 $ $ \end{document } one of the important applications of the ppn formalism is the comparison and classification of alternative metric theories of gravity .
the population of viable theories has fluctuated over the years as new effects and tests have been discovered , largely through the use of the ppn framework , which eliminated many theories thought previously to be viable .
the theory population has also fluctuated as new , potentially viable theories have been invented . in this article
, we shall focus on general relativity and the general class of scalar - tensor modifications of it , of which the jordan - fierz - brans - dicke theory ( brans - dicke , for short ) is the classic example .
the reasons are several - fold : 
 a full compendium of alternative theories is given in tegp 5 .many alternative metric theories developed during the 1970s and 1980s could be viewed as  straw - man  theories , invented to prove that such theories exist or to illustrate particular properties .
few of these could be regarded as well - motivated theories from the point of view , say , of field theory or particle physics . examples are the vector - tensor theories studied by will , nordtvedt and hellings.a number of theories fall into the class of  prior - geometric  theories , with absolute elements such as a flat background metric in addition to the physical metric .
most of these theories predict  preferred - frame  effects , that have been tightly constrained by observations ( see sec .
an example is rosen s bimetric theory.a large number of alternative theories of gravity predict gravitational wave emission substantially different from that of general relativity , in strong disagreement with observations of the binary pulsar ( see sec .
7).scalar - tensor modifications of gr have recently become very popular in unification schemes such as string theory , and in cosmological model building .
because the scalar fields are generally massive , the potentials in the post - newtonian limit will be modified by yukawa - like terms . 
 
many alternative metric theories developed during the 1970s and 1980s could be viewed as  straw - man 
few of these could be regarded as well - motivated theories from the point of view , say , of field theory or particle physics . examples are the vector - tensor theories studied by will , nordtvedt and hellings .
a number of theories fall into the class of  prior - geometric  theories , with absolute elements such as a flat background metric in addition to the physical metric .
most of these theories predict  preferred - frame  effects , that have been tightly constrained by observations ( see sec .
a large number of alternative theories of gravity predict gravitational wave emission substantially different from that of general relativity , in strong disagreement with observations of the binary pulsar ( see sec . 7 ) .
scalar - tensor modifications of gr have recently become very popular in unification schemes such as string theory , and in cosmological model building .
because the scalar fields are generally massive , the potentials in the post - newtonian limit will be modified by yukawa - like terms .
the metric g is the sole dynamical field , and the theory contains no arbitrary functions or parameters , apart from the value of the newtonian coupling constant g , which is measurable in laboratory experiments . throughout this article
, we ignore the cosmological constant . although  has significance for quantum field theory , quantum gravity , and cosmology , on the scale of the solar - system or of stellar systems its effects are negligible , for values of  corresponding to a cosmological closure density .
the field equations of gr are derivable from an invariant action principle i=0 , where 
 ( 23)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 i = { \left ( { 16\pi g } \right)^ { - 1}}\int { r{{\left ( { - g } \right)}^{1/2}}{d^4}x + { i_{\rm{m}}}\left ( { { \psi _ { \rm{m}}},{g_{\mu \nu } } } \right ) , } 
 $ $ \end{document } 
 where r is the ricci scalar , and i m is the matter action , which depends on matter fields m universally coupled to the metric g. by varying the action with respect to g , we obtain the field equations 
 ( 24)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 { g_{\mu \nu } } \equiv { r_{\mu \nu } } - \frac{1}{2}{g_{\mu \nu } } r = 8\pi g{t_{\mu \nu } } , 
 $ $ \end{document } 
 where t is the matter energy - momentum tensor .
general covariance of the matter action implies the equations of motion t;=0 ; varying i m with respect to m yields the matter field equations . by virtue of the absence of prior - geometric elements ,
the equations of motion are also a consequence of the field equations via the bianchi identities g;=0 . 
 
the general procedure for deriving the post - newtonian limit is spelled out in tegp 5.1 , and is described in detail for gr in tegp 5.2 .
table 3metric theories and their ppn parameter values ( 3=i=0 for all cases).theoryarbitrary functions or constantscosmic matching parametersppn parameters 
  

  

  

  
 1 

  
 2 
 general relativitynonenone11000scalar - tensorbrans - dicke 
  

  
 0 

 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\frac{{1 + \omega } } { { 2 + \omega } } $ $ \end{document } 
 1000generala( ) , v( ) 
  
 0 

 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\frac{{1 + \omega } } { { 2 + \omega } } $ $ \end{document } 
 1+000rosen s bimetricnonec0 , c11100 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\frac{{{c_0}}}{{{c_1 } } } - 1$$\end{document } 

  metric theories and their ppn parameter values ( 3=i=0 for all cases ) .
these theories contain the metric g , a scalar field  , a potential function v( ) , and a coupling function a( ) ( generalizations to more than one scalar field have also been carried out ) .
for some purposes , the action is conveniently written in a non - metric representation , sometimes denoted the 
einstein frame  , in which the gravitational action looks exactly like that of gr : 
 ( 25)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 \tilde i = { \left ( { 16\pi g } \right)^ { - 1}}\int { \left [ { \tilde r - 2{{\tilde g}^{\mu \nu } } { \partial _ \mu } \varphi { \partial _ \nu } \varphi - v\left ( \varphi \right ) } \right ] } { \left ( { - \tilde g } \right)^{1/2}}{d^4}x + { i_{\rm{m}}}\left ( { { \psi _ { \rm{m}}},{a^2}\left ( \varphi \right){{\tilde g}_{\mu \nu } } } \right ) , 
 $ $ \end{document } 
 where rgr is the ricci scalar of the  einstein 
metric g. ( apart from the scalar potential term v( ) , this corresponds to eq .
one because the matter fields m couple to a combination of  and g. despite appearances , however , it is a metric theory , because it can be put into a metric representation by identifying the  physical metric  
 ( 26)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 { g_{\mu \nu } } \equiv { a^2}\left ( \varphi \right){\tilde g_{\mu \nu } } . 

$ $ \end{document } 
 the action can then be rewritten in the metric form 
 ( 27)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 i = { \left ( { 16\pi g } \right)^ { - 1}}\int { \left [ { \phi r - { \phi ^ { - 1}}\omega \left ( \phi \right){g^{\mu \nu } } { \partial _ \mu } \phi { \partial _ \nu } \phi - { \phi ^2}v } \right ] } { \left ( { - g } \right)^{1/2}}{d^4}x + { i_{\rm{m}}}\left ( { { \psi _ { \rm{m}}},{g_{\mu \nu } } } \right ) , 
 $ $ \end{document } 
 where 
 ( 28)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 \begin{array}{*{20}{c } } 
 { \phi \equiv a{{\left ( \varphi \right)}^ { - 2}},}\\ 
 { 3 + 2\omega \left ( \phi \right ) \equiv \alpha { { \left ( \varphi \right)}^ { - 2}},}\\ 
 { \alpha \left ( \varphi \right ) \equiv d\left ( { \ln a\left ( \varphi \right ) } \right)/d\varphi . } 
 \end{array } 
 $ $ \end{document } 
 the einstein frame is useful for discussing general characteristics of such theories , and for some cosmological applications , while the metric representation is most useful for calculating observable effects .
the field equations , post - newtonian limit and ppn parameters are discussed in tegp 5.3 , and the values of the ppn parameters are listed in table 3 . 
 
the parameters that enter the post - newtonian limit are 
 ( 29)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 \begin{array}{*{20}{c } } 
 { \omega \equiv \omega \left ( { { \phi _ 0 } } \right),}&{\lambda \equiv { { \left [ { \left ( { d\omega /d\phi } \right){{\left ( { 3 + 2\omega } \right)}^ { - 2}}{{\left ( { 4 + 2\omega } \right)}^ { - 1 } } } \right]}_{\phi 0 } } , } 
 \end{array } 
 $ $ \end{document } 
 where 0 is the value of  today far from the system being studied , as determined by appropriate cosmological boundary conditions .
the following formula is also useful : 1/(2+)=202/(1+02 ) . in brans - dicke theory ( ()=constant ) , the larger the value of  , the smaller the effects of the scalar field , and in the limit (00 ) , the theory becomes indistinguishable from gr in all its predictions . in more general theories , the function ( )
could have the property that , at the present epoch , and in weak - field situations , the value of the scalar field 0 is such that  is very large and  is very small ( theory almost identical to gr today ) , but that for past or future values of  , or in strong - field regions such as the interiors of neutron stars ,  and  could take on values that would lead to significant differences from gr .
indeed , damour and nordtvedt have shown that in such general scalar - tensor theories , gr is a natural  attractor  : regardless of how different the theory may be from gr in the early universe ( apart from special cases ) , cosmological evolution naturally drives the fields toward small values of the function  , thence to large . estimates of the expected relic deviations from gr today in such theories depend on the cosmological model , but range from 10 to a few times 10 for 1 [ 47 , 48 ] . 
 
scalar fields coupled to gravity or matter are also ubiquitous in particle - physics - inspired models of unification , such as string theory . in some models
, the coupling to matter may lead to violations of wep , which are tested by etvs - type experiments . in many models
the scalar field is massive ; if the compton wavelength is of macroscopic scale , its effects are those of a  fifth force  .
only if the theory can be cast as a metric theory with a scalar field of infinite range or of range long compared to the scale of the system in question ( solar system ) can the ppn framework be strictly applied .
if the mass of the scalar field is sufficiently large that its range is microscopic , then , on solar - system scales , the scalar field is suppressed , and the theory is essentially equivalent to general relativity .
this is the case , for example , in the  oscillating - g  models of accetta , steinhardt and will ( see ) , in which the potential function v( ) contains both quadratic ( mass ) and quartic ( self - interaction ) terms , causing the scalar field to oscillate ( the initial amplitude of oscillation is provided by an inflationary epoch ) ; high - frequency oscillations in the  effective  newtonian constant geffg/=ga( ) then result .
the energy density in the oscillating scalar field can be enough to provide a cosmological closure density without resorting to dark matter , yet the value of  today is so large that the theory s local predictions are experimentally indistinguishable from gr . in other models , explored by damour and esposito - farse , non - linear scalar - field couplings can lead to  spontaneous scalarization  inside strong - field objects such as neutron stars , leading to large deviations from gr , even in the limit of very large . with the ppn formalism in hand , we are now ready to confront gravitation theories with the results of solar - system experiments . in this section
we focus on tests of the parameter  , consisting of the deflection of light and the time delay of light . a light ray ( or photon ) which passes the sun at a distance d
is deflected by an angle 
 ( 30)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 \delta \theta = \frac{1}{2}\left ( { 1 + \gamma } \right)\left ( { 4{m _ \odot } /d } \right)\left [ { \left ( { 1 + \cos \phi } \right)/2 } \right ] 
 $ $ \end{document } 
 ( tegp 7.1 ) , where m is the mass of the sun and  is the angle between the earth - sun line and the incoming direction of the photon ( figure 4 ) . for a grazing ray , dd , 0 , and 
 ( 31)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 \delta \theta \approx \frac{1}{2}\left ( { 1 + \gamma } \right)1.''75 , 
 $ $ \end{document } 
 independent of the frequency of light .
another , more useful expression gives the change in the relative angular separation between an observed source of light and a nearby reference source as both rays pass near the sun : 
 ( 32)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 

\delta \theta = \frac{1}{2}\left ( { 1 + \gamma } \right)\left [ { - \frac{{4{m _ \odot } } } { d}\cos { \rm{\chi + } } \frac{{4{m _ \odot } } } { { { d_{\rm{r}}}}}\left ( { \frac{{1 + \cos { \phi _ { \rm{r}}}}}{2 } } \right ) } \right ] , 
 $ $ \end{document } 
 where d and dr are the distances of closest approach of the source and reference rays respectively , r is the angular separation between the sun and the reference source , and  is the angle between the sun - source and the sun - reference directions , projected on the plane of the sky ( figure 4 ) .
thus , for example , the relative angular separation between the two sources may vary if the line of sight of one of them passes near the sun ( dr , drd ,  varying with time ) . 

it is interesting to note that the classic derivations of the deflection of light that use only the principle of equivalence or the corpuscular theory of light yield only the  1/2  part of the coefficient in front of the expression in eq .
but the result of these calculations is the deflection of light relative to local straight lines , as established for example by rigid rods ; however , because of space curvature around the sun , determined by the ppn parameter  , local straight lines are bent relative to asymptotic straight lines far from the sun by just enough to yield the remaining factor  /2  . the first factor  1/2  holds in any metric theory ,
thus , calculations that purport to derive the full deflection using the equivalence principle alone are incorrect .
the prediction of the full bending of light by the sun was one of the great successes of einstein s gr .
eddington s confirmation of the bending of optical starlight observed during a solar eclipse in the first days following world war i helped make einstein famous .
however , the experiments of eddington and his co - workers had only 30 percent accuracy , and succeeding experiments were not much better : the results were scattered between one half and twice the einstein value ( figure 5 ) , and the accuracies were low . 

the arrows denote anomalously large values from early eclipse expeditions . the shapiro time - delay measurements using viking spacecraft yielded an agreement with gr to 0.1 percent , and vlbi light deflection measurements have reached 0.02 percent .
the arrows denote anomalously large values from early eclipse expeditions . the shapiro time - delay measurements using viking spacecraft yielded an agreement with gr to 0.1 percent , and vlbi light deflection measurements have reached 0.02 percent .
however , the development of vlbi , very - long - baseline radio interferometry , produced greatly improved determinations of the deflection of light .
these techniques now have the capability of measuring angular separations and changes in angles as small as 100 microarcseconds .
early measurements took advantage of a series of heavenly coincidences : each year , groups of strong quasistellar radio sources pass very close to the sun ( as seen from the earth ) , including the group 3c273 , 3c279 , and 3c48 , and the group 0111 + 02 , 0119 + 11 and 0116 + 08 .
as the earth moves in its orbit , changing the lines of sight of the quasars relative to the sun , the angular separation  between pairs of quasars varies ( eq .
the time variation in the quantities d , dr ,  and r in eq .
( 32 ) is determined using an accurate ephemeris for the earth and initial directions for the quasars , and the resulting prediction for  as a function of time is used as a basis for a least - squares fit of the measured  , with one of the fitted parameters being the coefficient 1/2(1+ ) .
a number of measurements of this kind over the period 19691975 yielded an accurate determination of the coefficient 1/2(1+ ) which has the value unity in gr . a 1995 vlbi measurement using 3c273 and 3c279 yielded ( 1+)/2=0.99960.0017 .
a recent series of transcontinental and intercontinental vlbi quasar and radio galaxy observations made primarily to monitor the earth s rotation (  vlbi  in figure 5 ) was sensitive to the deflection of light over almost the entire celestial sphere ( at 90 from the sun , the deflection is still 4 milliarcseconds )
analysis of observations made by the hipparcos optical astrometry satellite yielded a test at the level of 0.3 percent .
a vlbi measurement of the deflection of light by jupiter was reported ; the predicted deflection of about 300 microarcseconds was seen with about 50 percent accuracy .
a radar signal sent across the solar system past the sun to a planet or satellite and returned to the earth suffers an additional non - newtonian delay in its round - trip travel time , given by ( see figure 4 ) 
 ( 33)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 \delta t = 2\left ( { 1 + \gamma } \right){m _ \odot } \ln \left [ { \left ( { { r _ \oplus } + { { \rm{x } } _ \oplus } \cdot { \rm{n } } } \right)\left ( { { r_{\rm{e } } } - { { \rm{x}}_{\rm{e } } } \cdot { \rm{n } } } \right)/{d^2 } } \right ] , 
 $ $ \end{document } 
 where xe ( x ) are the vectors , and re ( rx ) are the distances from the sun to the source ( earth ) , respectively ( tegp 7.2 ) . for a ray which passes close to the sun , 
 ( 34)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 \delta t \approx \frac{1}{2}\left ( { 1 + \gamma } \right)\left
[ { 240 - 20\ln \left ( { { d^2}/r } \right ) } \right]\mu { \rm{s } } , 
 $ $ \end{document } 
 where d is the distance of closest approach of the ray in solar radii , and r is the distance of the planet or satellite from the sun , in astronomical units . 
  in the two decades following irwin shapiro s 1964 discovery of this effect as a theoretical consequence of general relativity ,
several high - precision measurements were made using radar ranging to targets passing through superior conjunction . since one does not have access to a 
newtonian  signal against which to compare the round - trip travel time of the observed signal , it is necessary to do a differential measurement of the variations in round - trip travel times as the target passes through superior conjunction , and to look for the logarithmic behavior of eq .
, one must take into account the variations in round - trip travel time due to the orbital motion of the target relative to the earth .
this is done by using radar - ranging ( and possibly other ) data on the target taken when it is far from superior conjunction ( i.e. when the time - delay term is negligible ) to determine an accurate ephemeris for the target , using the ephemeris to predict the ppn coordinate trajectory xe(t ) near superior conjunction , then combining that trajectory with the trajectory of the earth x(t ) to determine the newtonian round - trip time and the logarithmic term in eq .
the resulting predicted round - trip travel times in terms of the unknown coefficient 1/2(1+ ) are then fit to the measured travel times using the method of least - squares , and an estimate obtained for 1/2(1+ ) .
the targets employed included planets , such as mercury or venus , used as passive reflectors of the radar signals (  passive radar  ) , and artificial satellites , such as mariners 6 and 7 , voyager 2 , and the viking mars landers and orbiters , used as active retransmitters of the radar signals (  active radar  ) .
the results for the coefficient 1/2(1+ ) of all radar time - delay measurements performed to date ( including a measurement of the one - way time delay of signals from the millisecond pulsar psr 1937 + 21 ) are shown in figure 5 ( see tegp 7.2   for discussion and references ) .
the viking experiment resulted in a 0.1 percent measurement . from the results of vlbi light - deflection experiments
, we can conclude that the coeficient 1/2(1+ ) must be within at most 0.014 percent of unity .
the explanation of the anomalous perihelion shift of mercury s orbit was another of the triumphs of gr .
this had been an unsolved problem in celestial mechanics for over half a century , since the announcement by le verrier in 1859 that , after the perturbing effects of the planets on mercury s orbit had been accounted for , and after the effect of the precession of the equinoxes on the astronomical coordinate system had been subtracted , there remained in the data an unexplained advance in the perihelion of mercury .
a number of ad hoc proposals were made in an attempt to account for this excess , including , among others , the existence of a new planet vulcan near the sun , a ring of planetoids , a solar quadrupole moment and a deviation from the inverse - square law of gravitation , but none was successful .
general relativity accounted for the anomalous shift in a natural way without disturbing the agreement with other planetary observations .
the predicted advance per orbit  , including both relativistic ppn contributions and the newtonian contribution resulting from a possible solar quadrupole moment , is given by 
 ( 35)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 \delta \tilde \omega = \left ( { 6\pi m / p } \right)\left [ { \frac{1}{3}\left ( { 2 + 2\gamma - \beta } \right ) + \frac{1}{6}\left ( { 2{\alpha _ 1 } - { \alpha _ 2 } + { \alpha _ 3 } + 2{\zeta _ 2 } } \right)\mu /m + { j_2}\left ( { { r^2}/2mp } \right ) } \right ] , 
 $ $ \end{document } 
 where mm1+m2 and m1m2/m are the total mass and reduced mass of the two - body system respectively ; pa(1e ) is the semi - latus rectum of the orbit , with the semi - major axis a and the eccentricity e ; r is the mean radius of the oblate body ; and j2 is a dimensionless measure of its quadrupole moment , given by j2=(ca)/m1r , where c and a are the moments of inertia about the body s rotation and equatorial axes , respectively ( for details of the derivation see tegp 7.3 ) . we have ignored preferred - frame and galaxy - induced contributions to  ; these are discussed in tegp 8.3 . 
  the first term in eq .
( 35 ) is the classical relativistic perihelion shift , which depends upon the ppn parameters  and . the second term depends upon the ratio of the masses of the two bodies ; it is zero in any fully conservative theory of gravity ( 12320 ) ; it is also negligible for mercury , since /mmmerc / m210 .
the third term depends upon the solar quadrupole moment j2 . for a sun that rotates uniformly with its observed surface angular velocity , so that the quadrupole moment is produced by centrifugal flattening
this actually agrees reasonably well with values inferred from rotating solar models that are in accord with observations of the normal modes of solar oscillations ( helioseismology ) . substituting standard orbital elements and physical constants for mercury and the sun we obtain the rate of perihelion shift \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\dot \tilde \omega$$\end{document } , in seconds of arc per century , 
 ( 36)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 \dot \tilde \omega = 42.''98\left [ { \frac{1}{3}\left ( { 2 + 2\gamma - \beta } \right ) + 3 \times { { 10}^ { - 4}}\left ( { { j_2}/{{10}^ { - 7 } } } \right ) } \right ] . 

$ $ \end{document } 
 now , the measured perihelion shift of mercury is known accurately : after the perturbing effects of the other planets have been accounted for , the excess shift is known to about 0.1 percent from radar observations of mercury between 1966 and 1990 .
the solar oblateness effect is smaller than the observational error , so we obtain the ppn bound |21|<310 . 
  the next class of solar - system experiments that test relativistic gravitational effects may be called tests of the strong equivalence principle ( sep ) . in sec .
3.1.2 we pointed out that many metric theories of gravity ( perhaps all except gr ) can be expected to violate one or more aspects of sep . among the testable violations of sep
are a violation of the weak equivalence principle for gravitating bodies that leads to perturbations in the earth - moon orbit ; preferred - location and preferred - frame effects in the locally measured gravitational constant that could produce observable geophysical effects ; and possible variations in the gravitational constant over cosmological timescales . in a pioneering calculation using his early form of the ppn formalism , nordtvedt   showed that many metric theories of gravity predict that massive bodies violate the weak equivalence principle  that is , fall with different accelerations depending on their gravitational self - energy .
dicke   linked such an effect to the possibility of a spatially varying gravitational constant , in theories such as scalar - tensor gravity . for a spherically symmetric body , the acceleration from rest in an external gravitational potential u has the form 
 ( 37)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 \begin{array}{*{20}{c } } 
 { { \rm{a = } } \left ( { { m_{\rm{p}}}/m } \right)\nabla u,}\\ 
 { { m_{\rm{p}}}/m = 1 - \eta \left ( { { e_{\rm{g}}}/m } \right),}\\ 
 { \eta = 4\beta - \gamma - 3 - \frac{{10}}{3}\xi - { \alpha _ 1 } + \frac{2}{3}{\alpha _ 2 } - \frac{2}{3}{\zeta _ 1 } - \frac{1}{3}{\zeta _ 2 } , } 
 \end{array } 
 $ $ \end{document } 
 where eg is the negative of the gravitational self - energy of the body ( eg>0 ) . this violation of the massive - body equivalence principle is known as the  nordtvedt effect  . the effect is absent in gr ( =0 ) but present in scalartensor theory ( =1/(2+)+4 )
. the existence of the nordtvedt effect does not violate the results of laboratory etvs experiments , since for laboratory - sized objects , eg / m10 , far below the sensitivity of current or future experiments .
however , for astronomical bodies , eg / m may be significant ( 10 for the sun , 10 for jupiter , 4.610 for the earth , 0.210 for the moon ) .
if the nordtvedt effect is present ( 0 ) then the earth should fall toward the sun with a slightly different acceleration than the moon .
this perturbation in the earth - moon orbit leads to a polarization of the orbit that is directed toward the sun as it moves around the earth - moon system , as seen from earth .
this polarization represents a perturbation in the earth - moon distance of the form 
 ( 38)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 \delta r = 13.1\eta \cos \left ( { { \omega _ 0 } - { \omega _ { \rm{s } } } } \right)t\left [ { \rm{m } } \right ] , 
 $ $ \end{document } 
 where 0 and s are the angular frequencies of the orbits of the moon and sun around the earth ( see tegp 8.1   for detailed derivations and references ; for improved calculations of the numerical coefficient , see [ 102 , 53 ] ) . 
 
since august 1969 , when the first successful acquisition was made of a laser signal reflected from the apollo 11 retroreflector on the moon , the lunar laser - ranging experiment ( lure ) has made regular measurements of the round - trip travel times of laser pulses between a network of observatories and the lunar retroreflectors , with accuracies that are approaching 50 ps ( 1 cm ) .
these measurements are fit using the method of least - squares to a theoretical model for the lunar motion that takes into account perturbations due to the sun and the other planets , tidal interactions , and post - newtonian gravitational effects .
the predicted round - trip travel times between retroreflector and telescope also take into account the librations of the moon , the orientation of the earth , the location of the observatory , and atmospheric effects on the signal propagation .
the  nordtvedt  parameter  along with several other important parameters of the model are then estimated in the least - squares method .
several independent analyses of the data found no evidence , within experimental uncertainty , for the nordtvedt effect ( for recent results see [ 56 , 154 , 96 ] ) . their results can be summarized by the bound ||<0.001 .
these results represent a limit on a possible violation of wep for massive bodies of 5 parts in 10 ( compare figure 1 ) . for brans - dicke theory
note that , at this level of precision , one can not regard the results of lunar laser ranging as a 
test of sep until one eliminates the possibility of a compensating violation of wep for the two bodies , because the chemical compositions of the earth and moon differ in the relative fractions of iron and silicates . to this end ,
the et - wash group carried out an improved test of wep for laboratory bodies whose chemical compositions mimic that of the earth and moon .
the resulting bound of four parts in 10   reduces the ambiguity in the lunar laser ranging bound , and establishes the firm limit on the universality of acceleration of gravitational binding energy at the level of <1.310 . in gr , the nordtvedt effect vanishes ; at the level of several centimeters and below , a number of non - null general relativistic effects should be present .
some theories of gravity violate sep by predicting that the outcomes of local gravitational experiments may depend on the velocity of the laboratory relative to the mean rest frame of the universe ( preferred - frame effects ) or on the location of the laboratory relative to a nearby gravitating body ( preferred - location effects ) . in the post - newtonian limit
, preferred - frame effects are governed by the values of the ppn parameters 1 , 2 , and 3 , and some preferred - location effects are governed by  ( see table 2 ) .
the most important such effects are variations and anisotropies in the locally - measured value of the gravitational constant , which lead to anomalous earth tides and variations in the earth s rotation rate ; anomalous contributions to the orbital dynamics of planets and the moon ; self - accelerations of pulsars , and anomalous torques on the sun that would cause its spin axis to be randomly oriented relative to the ecliptic ( see tegp 8.2 , 8.3 , 9.3 and 14.3 ( c ) ) .
an improved bound on 3 of 210 from the period derivatives of 20 millisecond pulsars was reported in ; improved bounds on 1 were achieved using lunar laser ranging data , and using observations of the circular binary orbit of the pulsar j2317 + 1439 .
negative searches for these effects have produced strong constraints on the ppn parameters ( table 4 ) . 

here  is a combination of other parameters given by =4310/3122/321/32/3.parametereffectlimitremarks1time delay210viking ranginglight deflection310vlbi1perihelion shift310j2=10 from helioseismologynordtvedt effect610=43 assumed 
  
 earth tides10gravimeter data 
  
 1 
 orbital polarization10lunar laser ranging210psr j2317 + 1439 
  
 2 
 spin precession410solar alignment with ecliptic 
  
 3 
 pulsar acceleration210pulsar  statistics 
  
 nordtvedt effect10lunar laser ranging 
  
 1 
 -210combined ppn bounds 
  
 2 
 binary acceleration410pp for psr 1913 + 16 
  
 3 
 newton s 3rd law10lunar acceleration 
  
 4 
 --not independent 
  current limits on the ppn parameters . here
most theories of gravity that violate sep predict that the locally measured newtonian gravitational constant may vary with time as the universe evolves . for the scalar - tensor theories listed in table 3
, the predictions for  / g can be written in terms of time derivatives of the asymptotic scalar field . where g does change with cosmic evolution ,
its rate of variation should be of the order of the expansion rate of the universe , i.e.  / gh0 , where h0 is the hubble expansion parameter and is given by h0=100 h km s mpc = h10 yr , where current observations of the expansion of the universe give h0.7 .
several observational constraints can be placed on  / g using methods that include studies of the evolution of the sun , observations of lunar occultations ( including analyses of ancient eclipse data ) , lunar laser - ranging measurements , planetary radar - ranging measurements , and pulsar timing data .
laboratory experiments may one day lead to interesting limits ( for review and references to past work see tegp 8.4 and 14.3 ( c ) ) .
table 5 
 constancy of the gravitational constant . for the pulsar data , the bounds are dependent upon the theory of gravity in the strong - field regime and on neutron star equation of state . 

method / g(10 yrlunar laser ranging0835viking radar24 - 210binary pulsar1111pulsar psr 0655 + 64<55 
 
 constancy of the gravitational constant . for the pulsar data , the bounds are dependent upon the theory of gravity in the strong - field regime and on neutron star equation of state . 
 
the best limits on  / g still come from ranging measurements to the viking landers and lunar laser ranging measurements [ 56 , 154 , 96 ] .
it has been suggested that radar observations of a mercury orbiter over a two - year mission ( 30 cm accuracy in range ) could yield ( / g)10 yr .
although bounds on  / g from solar - system measurements can be correctly obtained in a phenomenological manner through the simple expedient of replacing g by g0+0(tt0 ) in newton s equations of motion , the same does not hold true for pulsar and binary pulsar timing measurements .
the reason is that , in theories of gravity that violate sep , such as scalar - tensor theories , the  mass  and moment of inertia of a gravitationally bound body may vary with variation in g. because neutron stars are highly relativistic , the fractional variation in these quantities can be comparable to g / g , the precise variation depending both on the equation of state of neutron star matter and on the theory of gravity in the strong - field regime .
the variation in the moment of inertia affects the spin rate of the pulsar , while the variation in the mass can affect the orbital period in a manner that can subtract from the direct effect of a variation in g , given by b / pb=-2 / g .
thus , the bounds quoted in table 5 for the binary pulsar psr 1913 + 16   and the pulsar psr 0655 + 64   are theory - dependent and must be treated as merely suggestive . of the five  conservation law 
ppn parameters 1 , 2 , 3 , 4 , and 3 , only three , 2 , 3 and 3 , have been constrained directly with any precision ; 1 is constrained indirectly through its appearance in the nordtvedt effect parameter  , eq .
there is strong theoretical evidence that 4 , which is related to the gravity generated by fluid pressure , is not really an independent parameter  in any reasonable theory of gravity there should be a connection between the gravity produced by kinetic energy ( v ) , internal energy (  ) , and pressure ( p ) . from such considerations , there follows   the additional theoretical constraint 
 ( 39)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 6{\zeta _ 4 } = 3{\alpha _ 3 } + 2{\zeta _ 1 } - 3{\zeta _
3}. 
 $ $ \end{document } 
  a non - zero value for any of these parameters would result in a violation of conservation of momentum , or of newton s third law in gravitating systems .
an alternative statement of newton s third law for gravitating systems is that the  active gravitational mass  , that is the mass that determines the gravitational potential exhibited by a body , should equal the  passive gravitational mass  , the mass that determines the force on a body in a gravitational field .
such an equality guarantees the equality of action and reaction and of conservation of momentum , at least in the newtonian limit .
a classic test of newton s third law for gravitating systems was carried out in 1968 by kreuzer , in which the gravitational attraction of fluorine and bromine were compared to a precision of 5 parts in 10 .
they noted that current understanding of the structure of the moon involves an iron - rich , aluminum - poor mantle whose center of mass is offset about 10 km from the center of mass of an aluminum - rich , iron - poor crust .
the direction of offset is toward the earth , about 14 to the east of the earth - moon line .
such a model accounts for the basaltic maria which face the earth , and the aluminum - rich highlands on the moon s far side , and for a 2 km offset between the observed center of mass and center of figure for the moon .
because of this asymmetry , a violation of newton s third law for aluminum and iron would result in a momentum non - conserving self - force on the moon , whose component along the orbital direction would contribute to the secular acceleration of the lunar orbit .
improved knowledge of the lunar orbit through lunar laser ranging , and a better understanding of tidal effects in the earth - moon system ( which also contribute to the secular acceleration ) through satellite data , severely limit any anomalous secular acceleration , with the resulting limit 
 ( 40)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 \left| { \frac{{{{\left ( { { m_{\rm{a}}}/{m_{\rm{p } } } } \right)}_{{\rm{al } } } } - { { \left ( { { m_{\rm{a}}}/{m_{\rm{p } } } } \right)}_{{\rm{fe}}}}}}{{{{\left ( { { m_{\rm{a}}}/{m_{\rm{p } } } } \right)}_{{\rm{fe } } } } } } } \right| < 4 \times { 10^ { - 12}}. 
 $ $ \end{document } 
 according to the ppn formalism , in a theory of gravity that violates conservation of momentum , but that obeys the constraint of eq .
( 39 ) , the electrostatic binding energy ee of an atomic nucleus could make a contribution to the ratio of active to passive mass of the form 
 ( 41)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 { m_{\rm{a } } } = { m_{\rm{p } } } + \frac{1}{2}{\zeta _
3}{e_{\rm{e}}}/{c^2}. 
 $ $ \end{document } 
 the resulting limit on 3 from the lunar experiment is 3<110 ( tegp 9.2 , 14.3 ( d ) ) . 
 
another consequence of a violation of conservation of momentum is a self - acceleration of the center of mass of a binary stellar system , given by 
 ( 42)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 { { \rm{a}}_{{\rm{cm } } } } = \frac{1}{2}\left ( { { \zeta _ 2 } + { \alpha _ 3 } } \right)\frac{m}{{{a^2}}}\frac{\mu } { a}\frac{{\delta m}}{m}\frac{e}{{{{\left ( { 1 - { e^2 } } \right)}^{3/2}}}}{{\rm{n}}_{\rm{p } } } , 
 $ $ \end{document } 
 where m = m1m2 , a is the semi - major axis , and np is a unit vector directed from the center of mass to the point of periastron of m1 ( tegp 9.3 )
. a consequence of this acceleration would be non - vanishing values for dp = dt , where p denotes the period of any intrinsic process in the system ( orbit , spectra , pulsar periods ) .
the observed upper limit on dpp / dt of the binary pulsar psr 1913 + 16 places a strong constraint on such an effect , resulting in the bound |3+2|<410 . since 3 has already been constrained to be much less than this ( table 4 )
a gyroscope moving through curved spacetime suffers a precession of its axis given by 
 ( 43)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 \begin{array}{*{20}{c } } 
 { d{\rm{s / d}}\tau { \rm { = } } { \omega _ { \rm{g } } } \times { \rm{s,}}}&{{\omega _ { \rm{g } } } = \left ( { \gamma + \frac{1}{2 } } \right){\rm{v } } \times \nabla u , } 
 \end{array } 
 $ $ \end{document } 
 where v is the velocity of the gyroscope , and u is the newtonian gravitational potential of the source ( tegp 9.1 ) .
the predicted precession is about 2 arcseconds per century , an effect first calculated by de sitter .
this effect has been measured to about 0.7 percent using lunar laser ranging data [ 56 , 154 ] . 
  for a gyroscope orbiting the earth
the stanford gyroscope experiment has as one of its goals the measurement of this effect to 510 ( see below ) ; if achieved , this would substantially improve the accuracy of the parameter . according to gr , moving or rotating matter should produce a contribution to the gravitational field that is the analogue of the magnetic field of a moving charge or a magnetic dipole . although gravitomagnetism plays a role in a variety of measured relativistic effects , it has not been seen to date , isolated from other post - newtonian effects ( for a discussion of the evidence for gravitomagnetism in solar system measurements and the binary pulsar , see [ 100 , 99 ] ) . the relativity gyroscope experiment ( gravity probe b or gp - b ) at stanford university , in collaboration with nasa and lockheed - martin corporation , is in the advanced stage of developing a space mission to detect this phenomenon directly .
a set of four superconducting - niobium - coated , spherical quartz gyroscopes will be flown in a low polar earth orbit , and the precession of the gyroscopes relative to the distant stars will be measured . in the ppn formalism , the predicted effect of gravitomagnetism is a precession ( also known as the lense - thirring effect , or the dragging of inertial frames ) , given by 
 ( 44)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 \begin{array}{*{20}{c } } 
 { d{\rm{s / d}}\tau { \rm { = } } { \omega _ { { \rm{lt } } } } \times { \rm{s}},}&{{\omega _ { { \rm{lt } } } } = - \frac{1}{2}\left ( { 1 + \gamma + \frac{1}{4}{\alpha _ 1 } } \right)\frac{{{\rm{j } } - 3{\rm{n}}\left ( { { \rm{n } } \cdot { \rm{j } } } \right)}}{{{r^3 } } } , } 
 \end{array } 
 $ $ \end{document } 
 where j is the angular momentum of the earth , n is a unit radial vector , and r is the distance from the center of the earth ( tegp 9.1 ) . for a polar orbit at about 650 km altitude , this leads to a secular angular precession at a rate 1/2(1++1/41)4210 arcsec / yr . the accuracy goal of the experiment is about 0.5 milliarcseconds per year .
the science instrument package and the spacecraft are in the final phases of assembly , with launch scheduled for 2002 . 
 
another proposal to look for an effect of gravitomagnetism is to measure the relative precession of the line of nodes of a pair of laser - ranged geodynamics satellites ( lageos ) , ideally with supplementary inclination angles ; the inclinations must be supplementary in order to cancel the dominant nodal precession caused by the earth s newtonian gravitational multipole moments .
unfortunately , the two existing lageos satellites are not in appropriately inclined orbits , and no plans exist at present to launch a third satellite in a supplementary orbit .
nevertheless , by combing nodal precession data from lageos i and ii with perigee advance data from the slightly eccentric orbit of lageos ii , ciufolini et al . reported a partial cancellation of multipole effects , and a resulting 20 percent confirmation of gr .
a number of advanced space missions have been proposed in which spacecraft orbiters or landers and improved tracking capabilities could lead to significant improvements in values of the ppn parameters , of j2 of the sun , and of  / g .
doppler tracking of the cassini spacecraft ( launched to orbit and study saturn in 1997 ) during its 2003 superior conjunction could measure  to a few parts in 10 , by measuring the time variation of the shapiro delay .
a mercury orbiter , in a two - year experiment , with 3 cm range capability , could yield improvements in the perihelion shift to a part in 10 , in  to 410 , in  / g to 10 yr , and in j2 to a few parts in 10 .
proposals are being developed , primarily in europe , for advanced space missions which will have tests of ppn parameters as key components , including gaia , a high - precision astrometric telescope ( successor to hipparcos ) , which could measure light - deflection and  to the 10 level .
nordtvedt   has argued that  grand fits  of large solar system range data sets , including ranging to mercury , mars and the moon , could yield substantially improved measurements of ppn parameters .
in the solar system , gravity is weak , in the sense that the newtonian gravitational potential and related variables ( u(x , t)vp/ ) are much smaller than unity everywhere .
this is the basis for the post - newtonian expansion and for the  parametrized post - newtonian  framework described in sec .
systems are those for which the simple 1pn approximation of the ppn framework is no longer appropriate .
this can occur in a number of situations : 
 the system may contain strongly relativistic objects , such as neutron stars or black holes , near and inside which 1 , and the post - newtonian approximation breaks down . nevertheless , under some circumstances , the orbital motion may be such that the interbody potential and orbital velocities still satisfy 1 so that a kind of post - newtonian approximation for the orbital motion might work ; however , the strong - field internal gravity of the bodies could ( especially in alternative theories of gravity ) leave imprints on the orbital motion.the evolution of the system may be affected by the emission of gravitational radiation .
the 1pn approximation does not contain the effects of gravitational radiation back - reaction . in the expression for
the metric given in box 2 , radiation back - reaction effects do not occur until o( ) in g00 , o( ) in g0i , and o( ) in gij .
consequently , in order to describe such systems , one must carry out a solution of the equations substantially beyond 1pn order , sufficient to incorporate the leading radiation damping terms at 2.5pn order.the system may be highly relativistic in its orbital motion , so that uv1 even for the interbody field and orbital velocity . systems like this include the late stage of the inspiral of binary systems of neutron stars or black holes , driven by gravitational radiation damping , prior to a merger and collapse to a final stationary state .
binary inspiral is one of the leading candidate sources for detection by a world - wide network of laser interferometric gravitational wave observatories nearing completion .
a proper description of such systems requires not only equations for the motion of the binary carried to extraordinarily high pn orders ( at least 3.5pn ) , but also requires equations for the far - zone gravitational waveform measured at the detector , that are equally accurate to high pn orders beyond the leading  quadrupole  approximation . 
 
the system may contain strongly relativistic objects , such as neutron stars or black holes , near and inside which 1 , and the post - newtonian approximation breaks down .
nevertheless , under some circumstances , the orbital motion may be such that the interbody potential and orbital velocities still satisfy 1 so that a kind of post - newtonian approximation for the orbital motion might work ; however , the strong - field internal gravity of the bodies could ( especially in alternative theories of gravity ) leave imprints on the orbital motion .
the 1pn approximation does not contain the effects of gravitational radiation back - reaction . in the expression for
the metric given in box 2 , radiation back - reaction effects do not occur until o( ) in g00 , o( ) in g0i , and o( ) in gij .
consequently , in order to describe such systems , one must carry out a solution of the equations substantially beyond 1pn order , sufficient to incorporate the leading radiation damping terms at 2.5pn order .
the system may be highly relativistic in its orbital motion , so that uv1 even for the interbody field and orbital velocity .
systems like this include the late stage of the inspiral of binary systems of neutron stars or black holes , driven by gravitational radiation damping , prior to a merger and collapse to a final stationary state .
binary inspiral is one of the leading candidate sources for detection by a world - wide network of laser interferometric gravitational wave observatories nearing completion .
a proper description of such systems requires not only equations for the motion of the binary carried to extraordinarily high pn orders ( at least 3.5pn ) , but also requires equations for the far - zone gravitational waveform measured at the detector , that are equally accurate to high pn orders beyond the leading  quadrupole  approximation .
of course , some systems can not be properly described by any post - newtonian approximation because their behavior is fundamentally controlled by strong gravity .
these include the imploding cores of supernovae , the final merger of two compact objects , the quasinormal - mode vibrations of neutron stars and black holes , the structure of rapidly rotating neutron stars , and so on .
this field of  numerical relativity  is a rapidly growing and maturing branch of gravitational physics , whose description is beyond the scope of this article .
when dealing with the motion and gravitational wave generation by orbiting bodies , one finds a remarkable simplification within general relativity . as long as the bodies are sufficiently well - separated
that one can ignore tidal interactions and other effects that depend upon the finite extent of the bodies ( such as their quadrupole and higher multipole moments ) , then all aspects of their orbital behavior and gravitational wave generation can be characterized by just two parameters : mass and angular momentum . whether their internal structure is highly relativistic , as in black holes or neutron stars , or non - relativistic as in the earth and sun
furthermore , both quantities are measurable in principle by examining the external gravitational field of the bodies , and make no reference whatsoever to their interiors .
it is a consequence of the strong equivalence principle ( sep ) , described in section 3.1.2 .
general relativity satisfies sep because it contains one and only one gravitational field , the spacetime metric g. consider the motion of a body in a binary system , whose size is small compared to the binary separation . surround the body by a region that is large compared to the size of the body , yet small compared to the separation .
because of the general covariance of the theory , one can choose a freely - falling coordinate system which comoves with the body , whose spacetime metric takes the minkowski form at its outer boundary ( ignoring tidal effects generated by the companion ) .
there is thus no evidence of the presence of the companion body , and the structure of the chosen body can be obtained using the field equations of gr in this coordinate system .
far from the chosen body , the metric is characterized by the mass and angular momentum ( assuming that one ignores quadrupole and higher multipole moments of the body ) as measured far from the body using orbiting test particles and gyroscopes .
a black hole of mass m and a planet of mass m would produce identical spacetimes in this outer region .
the geometry of this region surrounding the one body must be matched to the geometry provided by the companion body .
einstein s equations provide consistency conditions for this matching that yield constraints on the motion of the bodies .
the motion of two planets of mass and angular momentum m1 , m2 , j1 and j2 is identical to that of two black holes of the same mass and angular momentum ( again , ignoring tidal effects ) .
this effacement does not occur in an alternative gravitional theory like scalar - tensor gravity .
there , in addition to the spacetime metric , a scalar field  is generated by the masses of the bodies , and controls the local value of the gravitational coupling constant ( i.e. g is a function of  ) .
now , in the local frame surrounding one of the bodies in our binary system , while the metric can still be made minkowskian far away , the scalar field will take on a value 0 determined by the companion body .
this can affect the value of g inside the chosen body , alter its internal structure ( specifically its gravitational binding energy ) and hence alter its mass .
effectively , each mass becomes several functions ma( ) of the value of the scalar field at its location , and several distinct masses come into play , inertial mass , gravitational mass ,  radiation  mass , etc
. the precise nature of the functions will depend on the body , specifically on its gravitational binding energy , and as a result , the motion and gravitational radiation may depend on the internal structure of each body . for compact bodies such as neutron stars , and black holes these internal structure effects could be large ; for example , the gravitational binding energy of a neutron star can be 40 percent of its total mass . at 1pn order ,
this is how the study of orbiting systems containing compact objects provides strong - field tests of general relativity .
even though the strong - field nature of the bodies is effaced in gr , it is not in other theories , thus any result in agreement with the predictions of gr constitutes a kind of  null  test of strong - field gravity .
the motion of bodies and the generation of gravitational radiation are long - standing problems that date back to the first years following the publication of gr , when einstein calculated the gravitational radiation emitted by a laboratory - scale object using the linearized version of gr , and de sitter calculated n - body equations of motion for bodies in the 1pn approximation to gr .
it has at times been controversial , with disputes over such issues as whether einstein s equations alone imply equations of motion for bodies ( einstein , infeld and hoffman demonstrated explicitly that they do , using a matching procedure similar to the one described above ) , whether gravitational waves are real or are artifacts of general covariance ( einstein waffled ; bondi and colleagues proved their reality rigorously in the 1950s ) , and even over algebraic errors ( einstein erred by a factor of 2 in his first radiation calculation ; eddington found the mistake ) .
shortly after the discovery of the binary pulsar psr 1913 + 16 in 1974 , questions were raised about the foundations of the  quadrupole formula  for gravitational radiation damping ( and in some quarters , even about its quantitative validity ) .
these questions were answered in part by theoretical work designed to shore up the foundations of the quadrupole approximation , and in part ( perhaps mostly ) by the agreement between the predictions of the quadrupole formula and the observed rate of damping of the pulsar s orbit ( see section 5.1 ) .
the problem of motion and radiation has received renewed interest since 1990 , with proposals for construction of large - scale laser interferometric gravitational - wave observatories , such as the ligo project in the us , virgo and geo600 in europe , and tama300 in japan , and the realization that a leading candidate source of detectable waves would be the inspiral , driven by gravitational radiation damping , of a binary system of compact objects ( neutron stars or black holes ) [ 1 , 127 ] .
the analysis of signals from such systems will require theoretical predictions from gr that are extremely accurate , well beyond the leading - order prediction of newtonian or even post - newtonian gravity for the orbits , and well beyond the leading - order formulae for gravitational waves .
this presented a major theoretical challenge : to calculate the motion and radiation of systems of compact objects to very high pn order , a formidable algebraic task , while addressing a number of issues of principle that have historically plagued this subject , sufficiently well to ensure that the results were physically meaningful .
this challenge is in the process of being met , so that we may soon see a remarkable convergence between observational data and accurate predictions of gravitational theory that could provide new , strong - field tests of gr .
the einstein equations g=8t are elegant and deceptively simple , showing geometry ( in the form of the einstein tensor g , which is a function of spacetime curvature ) being generated by matter ( in the form of the material stress - energy tensor t ) . however , this is not the most useful form for actual calculations . for post - newtonian calculations ,
a far more useful form is the so - called  relaxed  einstein equations : 
 ( 45)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 \square{h^{\alpha \beta } } = - 16\pi { \tau ^{\alpha \beta } } , 
 $ $ \end{document } 
 where -/t+ is the flat - spacetime wave operator , h is a
 gravitational tensor potential  related to the deviation of the spacetime metric from its minkowski form by the formula h(-g)g , g is the determinant of g , and a particular coordinate system has been specified by the de donder or harmonic gauge condition h/x=0 ( summation on repeated indices is assumed ) .
this form of einstein s equations bears a striking similarity to maxwell s equations for the vector potential a in lorentz gauge : a=-4j , a/x=0 .
there is a key difference , however : the source on the right hand side of eq .
stress - energy pseudotensor 
 ( 46)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 { \tau ^{\alpha \beta } } = \left ( { - g } \right){t^{\alpha \beta } } + { \left ( { 16\pi } \right)^ { - 1}}{\lambda ^{\alpha \beta } } , 
 $ $ \end{document } 
 where  is the non - linear  field  contribution given by terms quadratic ( and higher ) in h and its derivatives ( see , eqs .
, the gravitational field itself generates gravity , a reflection of the nonlinearity of einstein s equations , and in contrast to the linearity of maxwell s equations . 
 
equation ( 45 ) is exact , and depends only on the assumption that spacetime can be covered by harmonic coordinates .
it is called  relaxed  because it can be solved formally as a functional of source variables without specifying the motion of the source , in the form 
 ( 47)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 { h^{\alpha \beta } } \left ( { t,{\rm{x } } } \right ) = 4\int_c { \frac{{{\tau ^{\alpha \beta } } \left ( { t - \left| { { \rm{x } } - { \rm{x ' } } } \right|,{\rm{x ' } } } \right)}}{{\left| { { \rm{x } } - { \rm{x ' } } } \right| } } } { d^3}x ' , 
 $ $ \end{document } 
 where the integration is over the past flat - spacetime null cone c of the field point ( t , x ) . the motion of the source is then determined either by the equation /x=0 ( which follows from the harmonic gauge condition ) , or from the usual covariant equation of motion t;=0 , where the subscript ;  denotes a covariant divergence .
this formal solution can then be iterated in a slow motion ( v<1 ) weak - field ( h<1 ) approximation .
one begins by substituting h0=0 into the source  in eq . ( 47 ) , and solving for the first iterate h1 , and then repeating the procedure sufficiently many times to achieve a solution of the desired accuracy .
for example , to obtain the 1pn equations of motion , two iterations are needed ( i.e. h2 must be calculated ) ; likewise , to obtain the leading gravitational waveform for a binary system , two iterations are needed . 
  at the same time , just as in electromagnetism , the formal integral ( 47 ) must be handled differently , depending on whether the field point is in the far zone or the near zone . for field points in the far zone or radiation zone , |x|>>|x| (  is the gravitational wavelength/2 ) , the field can be expanded in inverse powers of r=|x| in a multipole expansion , evaluated at the  retarded time
the leading term in 1/r is the gravitational waveform . for field points in the near zone or induction zone , |x||x|< , the field is expanded in powers of |xx| about the local time t , yielding instantaneous potentials that go into the equations of motion .
however , because the source  contains h itself , it is not confined to a compact region , but extends over all spacetime . as a result
, there is a danger that the integrals involved in the various expansions will diverge or be ill - defined .
this consequence of the non - linearity of einstein s equations has bedeviled the subject of gravitational radiation for decades .
the  post - minkowskian  method of blanchet , damour and iyer [ 19 , 20 , 21 , 45 , 22 , 15 ] solves einstein s equations by two different techniques , one in the near zone and one in the far zone , and uses the method of singular asymptotic matching to join the solutions in an overlap region .
( dire ) approach of will , wiseman and pati [ 152 , 105 ] retains eq . (
47 ) as the global solution , but splits the integration into one over the near zone and another over the far zone , and uses different integration variables to carry out the explicit integrals over the two zones . in the dire method ,
these methods assume from the outset that gravity is sufficiently weak that h<1 and harmonic coordinates exists everywhere , including inside the bodies .
thus , in order to apply the results to cases where the bodies may be neutron stars or black holes , one relies upon the strong equivalence principle to argue that , if tidal forces are ignored , and equations are expressed in terms of masses and spins , one can simply extrapolate the results unchanged to the situation where the bodies are ultrarelativistic .
while no general proof of this exists , it has been shown to be valid in specific circumstances , such as at 2pn order in the equations of motion , and for black holes moving in a newtonian background field .
methods such as these have resolved most of the issues that led to criticism of the foundations of gravitational radiation theory during the 1970s . among the results of these approaches are formulae for the equations of motion and gravitational waveform of binary systems of compact objects , carried out to high orders in a pn expansion . here
for example , the relative two - body equation of motion has the form 
 ( 48)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 { \rm{a } } = \frac{{d{\rm{v}}}}{{dt } } = \frac{m}{{{r^2}}}\left\ { { - { \rm{\hat n } } + { { \rm{a}}_{1{\rm{pn } } } } + { { \rm{a}}_{{\rm{2pn } } } } + { { \rm{a}}_{2.5{\rm{pn } } } } + { { \rm{a}}_{{\rm{3pn } } } } + { { \rm{a}}_{{\rm{3}}{\rm{.5pn } } } } + \ldots } \right\ } , 
 $ $ \end{document } 
 where m = m1+m2 is the total mass , r=|x1x2| , v = v1v2 , and n=(x1x2)/r .
the notation anpn indicates that the term is o( ) relative to the newtonian term -nn .
explicit formulae for terms through various orders have been calculated by various authors : non - radiative terms through 2pn order [ 40 , 38 , 72 , 39 , 27 ] , radiation reaction terms at 2.5pn and 3.5pn order [ 78 , 79 , 17 ] , and non - radiative 3pn terms [ 80 , 81 , 46 , 26 , 25 ] . here we quote only the first pn corrections and the leading radiation - reaction terms at 2.5pn order : 
 ( 49)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 { { \rm{a}}_{1{\rm{pn } } } } = \left\ { { \left ( { 4 + 2\eta } \right)\frac{m}{r } - \left ( { 1 + 3\eta } \right){v^2 } + \frac{3}{2}\eta { { \dot r}^2 } } \right\}{\rm{\hat n } } + \left ( { 4 - 2\eta } \right)\dot r{\rm{v , } } 
 $ $ \end{document } 
 ( 50)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 { { \rm{a}}_{2.5{\rm{pn } } } } = - \frac{8}{{15}}\eta \frac{m}{r}\left\ { { \left ( { 9{v^2 } + 17\frac{m}{r } } \right)\dot r{\rm{\hat n } } - \left ( { 3{v^2 } + 9\frac{m}{r } } \right){\rm{v } } } \right\ } , 
 $ $ \end{document } 
 where =m1m2/(m1+m2 )
. these terms are sufficient to analyse the orbit and evolution of the binary pulsar ( sec .
5.1 ) . for example , the 1pn terms are responsible for the periastron advance of an eccentric orbit , given by =6fbm / a(1e ) , where a and e are the semi - major axis and eccentricity , respectively , of the orbit , and fb is the orbital frequency , given to the needed order by kepler s third law 2fb=(m / a ) . 
 
another product is a formula for the gravitational field far from the system , written schematically in the form 
 ( 51)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 { h^{ij } } = \frac{{2m}}{r}\left\ { { { q^{ij } } + q_{0.5{\rm{pn}}}^{ij } + q_{1{\rm{pn}}}^{ij } + q_{1.5{\rm{pn}}}^{ij } + q_{2{\rm{pn}}}^{ij } + q_{2.5{\rm{pn}}}^{ij } + \ldots } \right\ } , 
 $ $ \end{document } 
 where r is the distance from the source , and the variables are to be evaluated at retarded time tr .
the leading term is the so - called quadrupole formula 
 ( 52)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 { h^{ij}}\left ( { t,{\rm{x } } } \right ) = \frac{2}{r}{\ddot i^{ij}}\left ( { t - r } \right ) , 
 $ $ \end{document } 
 where i is the quadrupole moment of the source , and overdots denote time derivatives . for a binary system
this leads to 
 ( 53)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 { q^{ij } } = 2\eta \left ( { { v^i}{v^j } - m{{\hat n}^i}{{\hat n}^j}/r } \right ) . 

$ $ \end{document } 
 for binary systems , explicit formulae for all the terms through 2.5pn order have been derived by various authors [ 135 , 155 , 106 , 23 , 24 , 152 , 16 , 18 ] . given the gravitational waveform
, one can compute the rate at which energy is carried off by the radiation ( schematically d , the gravitational analog of the poynting flux ) .
the lowest - order quadrupole formula leads to the gravitational wave energy flux 
 ( 54)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 \dot e = \frac{8}{{15}}{\eta ^2}\frac{{{m^4}}}{{{r^4}}}\left ( { 12{v^2 } - 11{{\dot r}^2 } } \right ) . 

$ $ \end{document } 
 formulae for fluxes of angular and linear momentum can also be derived .
the 2.5pn radiation - reaction terms in the equation of motion ( 48 ) result in a damping of the orbital energy that precisely balances the energy flux ( 54 ) determined from the waveform . averaged over one orbit , this results in a rate of increase of the binary s orbital frequency given by 
 ( 55)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 { \dot f_{\rm{b } } } = \frac{{192\pi } } { 5}f_{\rm{b}}^2{\left ( { 2\pi { \mathcal m}{f_{\rm{b } } } } \right)^{5/3}}f\left ( e \right ) , 
 $ $ \end{document } 
 where \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\mathcal m}$$\end{document }
is the so - called  chirp  mass , given by \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\mathcal m } = { \eta ^{3/5}}m$$\end{document } , and f(e)=(1 + 73e/24 + 37e/96)/(1e ) .
 
 f(t)dt of either the orbit or the gravitational waves ( for which f=2fb for the dominant component ) as a function of the frequency , one in effect measures the  chirp  mass of the system . 
 
these formalisms have also been generalized to include the leading effects of spin - orbit and spin - spin coupling between the bodies [ 83 , 82 ] .
another approach to gravitational radiation is applicable to the special limit in which one mass is much smaller than the other .
one begins with an exact background spacetime of a black hole , either the non - rotating schwarzschild or the rotating kerr solution , and perturbs it according to g=g(0)+h. the particle moves on a geodesic of the background spacetime , and a suitably defined source stress - energy tensor for the particle acts as a source for the gravitational perturbation and wave field h. this method provides numerical results that are exact in v , as well as analytical results expressed as series in powers of v , both for non - rotating and for rotating black holes . for non - rotating holes , the analytical expansions
all results of black hole perturbation agree precisely with the m10 limit of the pn results , up to the highest pn order where they can be compared ( for a detailed review see ) .
a gravitational wave detector can be modelled as a body of mass m at a distance l from a fiducial laboratory point , connected to the point by a spring of resonant frequency 0 and quality factor q. from the equation of geodesic deviation , the infinitesimal displacement  of the mass along the line of separation from its equilibrium position satisfies the equation of motion 
 ( 56)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 \ddot \xi + 2\frac{{{\omega _ 0}}}{q}\dot \xi + \omega _ 0 ^ 2\xi = \frac{l}{2}\left ( { { f _ + } \left ( { \theta , \phi , \psi } \right){{\ddot h } _
+ } \left ( t \right ) + { f _ \times } \left ( { \theta , \phi , \psi } \right){{\ddot h } _ \times } \left ( t \right ) } \right ) , 
 $ $ \end{document } 
 where f+(,, ) and f(,, ) are  beam - pattern  factors , that depend on the direction of the source ( , ) and on a polarization angle  , and h+(t ) and h(t ) are gravitational waveforms corresponding to the two polarizations of the gravitational wave ( for a review , see ) . in a source coordinate system
in which the xy plane is the plane of the sky and the z - direction points toward the detector , these two modes are given by 
 ( 57)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 \begin{array}{*{20}{c } } 
 { { h _ + } \left ( t \right ) = \frac{1}{2}\left ( { { h^{xx}}\left ( t \right ) - { h^{yy}}\left ( t \right ) } \right),}&{{h _ \times } \left ( t \right ) = { h^{xy}}\left ( t \right ) , } 
 \end{array } 
 $ $ \end{document } 
 where h represent transverse - traceless ( tt ) projections of the calculated waveform of eq .
( 51 ) , given by 
 ( 58)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 h_{{\rm{tt}}}^{ij } = { h^{kl}}\left [ { \left ( { \delta ik - { { \hat n}^i}{{\hat n}^k } } \right)\left ( { \delta jl - { { \hat n}^j}{{\hat n}^l } } \right ) - \frac{1}{2}\left ( { \delta ij - { { \hat n}^i}{{\hat n}^j } } \right)\left ( { \delta kl - { { \hat n}^k}{{\hat n}^l } } \right ) } \right ] , 
 $ $ \end{document } 
 where n is a unit vector pointing toward the detector .
the beam pattern factors depend on the orientation and nature of the detector . for a wave approaching along the laboratory z - direction , and for a mass
whose location on the xy plane makes an angle  with the x axis , the beam pattern factors are given by f+=cos2 and f=sin2. for a resonant cylinder oriented along the laboratory z axis , and for source direction (  ,  ) , they are given by f+=sincos2 , f=sinsin2 ( the angle  measures the relative orientation of the laboratory and source x - axes ) . for a laser interferometer with one arm along the laboratory x - axis ,
the other along the y - axis , and with  dedined as the differential displacement along the two arms , the beam pattern functions are f+=1/2(1+cos)cos2cos2cossin2sin2 and f=1/2(1+cos)cos2sin2+cossin2cos2. 
  the waveforms h+(t ) and h(t ) depend on the nature and evolution of the source .
for example , for a binary system in a circular orbit , with an inclination i relative to the plane of the sky , and the x - axis oriented along the major axis of the projected orbit , the quadrupole approximation of eq .
( 53 ) gives 
 ( 59)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 \begin{array}{*{20}{l } } 
 { { h _ + } \left ( t \right ) = - \frac{{2{\mathcal m}}}{r}{{\left ( { 2\pi { \mathcal m}{f_{\rm{b } } } } \right)}^{2/3}}\left ( { 1 + { { \cos } ^2}i } \right)\cos 2{\phi _ { \rm{b}}}\left ( t \right),}\\ 
 { { h _ \times } \left ( t \right ) = - \frac{{2{\mathcal m}}}{r}{{\left ( { 2\pi { \mathcal m}{f_{\rm{b } } } } \right)}^{2/3}}2\cos i\cos 2{\phi _ { \rm{b}}}\left ( t \right ) , } 
 \end{array } 
 $ $ \end{document } 
 where b(t)=2fb(t)dt is the orbital phase . 

the 1974 discovery of the binary pulsar psr 1913 + 16 by joseph taylor and russell hulse during a routine search for new pulsars provided the first possibility of probing new aspects of gravitational theory : the effects of strong relativistic internal gravitational fields on orbital dynamics , and the effects of gravitational radiation reaction . for reviews of the discovery and current status ,
see the published nobel prize lectures by hulse and taylor [ 76 , 123 ] . for
a thorough review of pulsars , including binary and millisecond pulsars , see the living review by dunc lorimer .
the system consists of a pulsar of nominal period 59 ms in a close binary orbit with an as yet unseen companion .
the orbital period is about 7.75 hours , and the eccentricity is 0.617 . from detailed analyses of the arrival times of pulses ( which amounts to an integrated version of the doppler - shift methods used in spectroscopic binary systems ) , extremely accurate orbital and physical parameters for the system
have been obtained ( table 6 ) . because the orbit is so close ( 1r ) and
because there is no evidence of an eclipse of the pulsar signal or of mass transfer from the companion , it is generally believed that the companion is compact : evolutionary arguments suggest that it is most likely a dead pulsar .
thus the orbital motion is very clean , free from tidal or other complicating effects .
clean  in the sense that by exploiting the intrinsic stability of the pulsar clock combined with the ability to maintain and transfer atomic time accurately using such devices as the global positioning system , the observers can keep track of the pulsar phase with an accuracy of 15 s , despite extended gaps between observing sessions ( including a several - year gap during the middle 1990s upgrade of the arecibo radio telescope ) .
the pulsar has shown no evidence of  glitches  in its pulse period . 

data taken from an online catalogue of pulsars maintained by stephen thorsett of the university of california , santa cruz , see [ 
 128 
 ] . 

parametersymbolvalue(units ) 
 ( i )  physical  parameters 
 right ascension 
  
 191528.00018(15)declination 
  
 160627.4043(3)pulsar periodpp ( ms)59.029997929613(7)derivative of period 
  
 p 
 8.62713(8)10 
 ( ii )  keplerian  parameters 
 projected semimajor axisap sin i ( s)2.3417592(19)eccentricity 
 e 
 0.6171308(4)orbital periodpb ( day)0.322997462736(7)longitude of periastron0 ( )226.57528(6)julian date of periastront0 ( mjd)46443.99588319(3 ) 
 ( iii )  post - keplerian 
parameters 
 mean rate of periastron advance (  yr)4.226621(11)redshift / time dilation ( ms)4.295(2)orbital period derivativeb ( 10)-2.422(6 ) 
 
 parameters of the binary pulsar psr 1913 + 16 .
data taken from an online catalogue of pulsars maintained by stephen thorsett of the university of california , santa cruz , see [ 
 128 
 ] . 
 
three factors make this system an arena where relativistic celestial mechanics must be used : the relatively large size of relativistic effects [ vorbit(m / r)10 ] , a factor of 10 larger than the corresponding values for solar - system orbits ; the short orbital period , allowing secular effects to build up rapidly ; and the cleanliness of the system , allowing accurate determinations of small effects . because the orbital separation is large compared to the neutron stars compact size
, tidal effects can be ignored . just as newtonian gravity is used as a tool for measuring astrophysical parameters of ordinary binary systems , so gr is used as a tool for measuring astrophysical parameters in the binary pulsar .
the observational parameters that are obtained from a least - squares solution of the arrival - time data fall into three groups : ( i ) non - orbital parameters , such as the pulsar period and its rate of change ( defined at a given epoch ) , and the position of the pulsar on the sky ; ( ii ) five  keplerian  parameters , most closely related to those appropriate for standard newtonian binary systems , such as the eccentricity e and the orbital period pb ; and ( iii ) five  post - keplerian  parameters .
the five post - keplerian parameters are :  , the average rate of periastron advance ;  , the amplitude of delays in arrival of pulses caused by the varying effects of the gravitational redshift and time dilation as the pulsar moves in its elliptical orbit at varying distances from the companion and with varying speeds ; b , the rate of change of orbital period , caused predominantly by gravitational radiation damping ; and r and s = sin i , respectively the  range  and  shape  of the shapiro time delay of the pulsar signal as it propagates through the curved spacetime region near the companion , where i is the angle of inclination of the orbit relative to the plane of the sky . in gr , these post - keplerian parameters can be related to the masses of the two bodies and to measured keplerian parameters by the equations ( tegp 12.1 , 14.6 ( a ) ) 
 ( 60)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 \begin{array}{*{20}{c } } 
 { \left\langle { \dot \omega } \right\rangle = 6\pi { f_{\rm{b}}}{{\left ( { 2\pi m{f_{\rm{b } } } } \right)}^{2/3}}{{\left ( { 1 - { e^2 } } \right)}^ { - 1}},}\\ 
 { \gamma ' = e{{\left ( { 2\pi { f_{\rm{b } } } } \right)}^ { - 1}}{{\left ( { 2\pi m{f_{\rm{b } } } } \right)}^{2/3}}\left ( { { m_2}/m } \right)\left ( { 1 + { m_2}/m } \right),}\\ 
 { { { \dot p}_{\rm{b } } } = - \left ( { 192\pi /5 } \right){{\left ( { 2\pi { \mathcal m}{f_{\rm{b } } } } \right)}^{5/3}}f\left ( e \right),}\\ 
 { s = \sin i,}\\ 
 { r = { m_{2 , } } } 
 \end{array } 
 $ $ \end{document } 
 where m1 and m2 denote the pulsar and companion masses , respectively .
the formula for  ignores possible non - relativistic contributions to the periastron shift , such as tidally or rotationally induced effects caused by the companion ( for discussion of these effects , see tegp 12.1 ( c ) ) .
the formula for b includes only quadrupole gravitational radiation ; it ignores other sources of energy loss , such as tidal dissipation ( tegp 12.1 ( f ) ) .
notice that , by virtue of kepler s third law , ( 2fb)=m / a , ( 2mfb)m / a , thus the first two post - keplerian parameters can be seen as o( ) , or 1pn corrections to the underlying variable , while the third is an o( ) , or 2.5pn correction .
the current observed values for the keplerian and post - keplerian parameters are shown in table 6 .
the parameters r and s are not separately measurable with interesting accuracy for psr 1913 + 16 because the orbit s 47 inclination does not lead to a substantial shapiro delay . 
  because fb and e are separately measured parameters , the measurement of the three post - keplerian parameters provide three constraints on the two unknown masses .
the periastron shift measures the total mass of the system , b measures the chirp mass , and  measures a complicated function of the masses .
gr passes the test if it provides a consistent solution to these constraints , within the measurement errors . from the intersection of the  and  constraints
( 60 ) then predicts the value b=-2.402430.0000510 . in order to compare the predicted value for b with the observed value of table 6 , it is necessary to take into account the small effect of a relative acceleration between the binary pulsar system and the solar system caused by the differential rotation of the galaxy
this effect was previously considered unimportant when b was known only to 10 percent accuracy .
damour and taylor   carried out a careful estimate of this effect using data on the location and proper motion of the pulsar , combined with the best information available on galactic rotation , and found 
 ( 61)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 \dot p_{\rm{b}}^{{\rm{gal } } } \simeq - \left ( { 1.7 \pm 0.5 } \right ) \times { 10^ { - 14}}. 
 $ $ \end{document } 
 subtracting this from the observed b ( table 6 ) gives the residual 
 ( 62)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 \dot p_{\rm{b}}^{{\rm{corr } } } = - \left ( { 2.408 \pm 0.010\left [ { { \rm{obs } } } \right ] \pm 0.005\left [ { { \rm{gal } } } \right ] } \right ) \times { 10^ { - 12 } } , 
 $ $ \end{document } 
 which agrees with the prediction within the errors .
in other words , 
 ( 63)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 \frac{{\dot p_{\rm{b}}^{{\rm{gr}}}}}{{\dot p_{\rm{b}}^{{\rm{corr } } } } } = 1.0023 \pm 0.0041\left [ { { \rm{obs } } } \right ] \pm 0.0021\left [ { { \rm{gal } } } \right ] . 

$ $ \end{document } 
 the consistency among the measurements is displayed in figure 6 , in which the regions allowed by the three most precise constraints have a single common overlap . 

figure 6 
 constraints on masses of the pulsar and its companion from data on psr 1913 + 16 , assuming gr to be valid .
the width of each strip in the plane reflects observational accuracy , shown as a percentage .
an inset shows the three constraints on the full mass plane ; the intersection region ( a ) has been magnified 400 times for the full figure . 

 

constraints on masses of the pulsar and its companion from data on psr 1913 + 16 , assuming gr to be valid .
the width of each strip in the plane reflects observational accuracy , shown as a percentage .
an inset shows the three constraints on the full mass plane ; the intersection region ( a ) has been magnified 400 times for the full figure . 
 
a third way to display the agreement with general relativity is by comparing the observed phase of the orbit with a theoretical template phase as a function of time . if fb varies slowly in time , then to first order in a taylor expansion , the orbital phase is given by b(t)=2fb0t+b0 t .
the time of periastron passage tp is given by (tp)=2n , where n is an integer , and consequently , the periastron time will not grow linearly with n. thus the cumulative difference between periastron time tp and n / fb0 , the quantities actually measured in practice , should vary according to tpn / fb0=b0n/2fb03-b0-2fb0)t .
figure 7 shows the results : the dots are the data points , while the curve is the predicted difference using the measured masses and the quadrupole formula for b0 . 

the gap during the middle 1990s was caused by a closure of arecibo for upgrading .
[ j. h. taylor and j. m. weisberg , 2000 , private communication . ] 

 

the gap during the middle 1990s was caused by a closure of arecibo for upgrading .
] 
  the consistency among the constraints provides a test of the assumption that the two bodies behave as  point 
masses , without complicated tidal effects , obeying the general relativistic equations of motion including gravitational radiation .
it is also a test of strong gravity , in that the highly relativistic internal structure of the neutron stars does not influence their orbital motion , as predicted by the strong equivalence principle of gr .
recent observations [ 84 , 138 ] indicate variations in the pulse profile , which suggests that the pulsar is undergoing precession as it moves through the curved spacetime generated by its companion , an effect known as geodetic precession .
the amount is consistent with gr , assuming that the pulsar s spin is suitably misaligned with the orbital angular momentum .
unfortunately , the evidence suggests that the pulsar beam may precess out of our line of sight by 2020 .
since 1990 , several new massive binary pulsars similar to psr 1913 + 16 were discovered , leading to the possibility of new or improved tests of gr . 

its pulses are significantly stronger and narrower than those of psr 1913 + 16 , so timing measurements are more precise , reaching 3 s accuracy .
the orbital plane appears to be almost edge - on relative to the line of sight ( i80 ) ; as a result the shapiro delay is substantial , and separate values of the parameters r and s have been obtained with interesting accuracy .
assuming general relativity , one infers that the two masses are m1=1.3350.002m and m2=1.3440.002m. the rate of orbit decay b agrees with gr to about 15 percent , the precision limited by the poorly known distance to the pulsar , which introduces a significant uncertainty into the subtraction of galactic acceleration . independently of b ,
measurement of the four other post - keplerian parameters gives two tests of strong - field gravity in the non - radiative regime . 

table 7parameters of other binary pulsars . from [ 156 , 118 , 119 ] and from an online catalogue of pulsars maintained by stephen thorsett of the university of california , santa cruz , see .parameterb1534 + 12b2127 + 11cb1855 + 09b0655 + 64 
 ( i )  keplerian 
parameters 
 ap sin i ( s)3.729464(3)2.520(3)9.2307802(4)4.125612(5 ) 
 e 
 0.2736775(5)0.68141(2)0.00002168(5)0.0000075(11)pb ( day)0.42073729933(3)0.335282052(6)12.3271711905(6)1.028669703(1 ) 
 ( ii )  post - keplerian  parameters 
  (  yr)1.755794(19)4.457(12) ( ms)2.071(6)4.9(1.1)b ( 10)-0.131(9)<0.5r(s)6.3(1.3)1.27(10)s = sin i0.983(8)0.9992(5 ) 
 psr 2127 + 11c .
this system appears to be a clone of the hulse - taylor binary pulsar , with very similar values for orbital period and eccentricity ( see table 7 ) .
the inferred total mass of the system is 2.7060.011m. because the system is in the globular cluster m15 ( ngc 7078 ) , it suffers doppler shifts resulting from local accelerations , either by the mean cluster gravitational field or by nearby stars , that are more difficult to estimate than was the case with the galactic system psr 1913 + 16 .
this may make a separate , precision measurement of the relativistic contribution to b impossible.psr 1855 + 09 .
this binary pulsar system is not particularly relativistic , with a long period ( 12 days ) and highly circular orbit .
however , because we observe the orbit nearly edge on , the shapiro delay is large and measurable , as reflected in the post - keplerian parameters r and s.psr 0655 + 64 .
this system consists of a pulsar and a white dwarf companion in a nearly circular orbit .
its pulses are significantly stronger and narrower than those of psr 1913 + 16 , so timing measurements are more precise , reaching 3 s accuracy .
the orbital plane appears to be almost edge - on relative to the line of sight ( i80 ) ; as a result the shapiro delay is substantial , and separate values of the parameters r and s have been obtained with interesting accuracy .
assuming general relativity , one infers that the two masses are m1=1.3350.002m and m2=1.3440.002m. the rate of orbit decay b agrees with gr to about 15 percent , the precision limited by the poorly known distance to the pulsar , which introduces a significant uncertainty into the subtraction of galactic acceleration .
independently of b , measurement of the four other post - keplerian parameters gives two tests of strong - field gravity in the non - radiative regime . 

table 7parameters of other binary pulsars . from [ 156 , 118 , 119 ] and from an online catalogue of pulsars maintained by stephen thorsett of the university of california , santa cruz , see .parameterb1534 + 12b2127 + 11cb1855 + 09b0655 + 64 
 ( i )  keplerian 
parameters 
 ap sin i ( s)3.729464(3)2.520(3)9.2307802(4)4.125612(5 ) 
 e 
 0.2736775(5)0.68141(2)0.00002168(5)0.0000075(11)pb ( day)0.42073729933(3)0.335282052(6)12.3271711905(6)1.028669703(1 ) 
 ( ii )  post - keplerian  parameters 
  (  yr)1.755794(19)4.457(12) ( ms)2.071(6)4.9(1.1)b ( 10)-0.131(9)<0.5r(s)6.3(1.3)1.27(10)s = sin i0.983(8)0.9992(5 ) 
  parameters of other binary pulsars . from [ 156 , 118 , 119 ] and from an online catalogue of pulsars maintained by stephen thorsett of the university of california , santa cruz , see .
this system appears to be a clone of the hulse - taylor binary pulsar , with very similar values for orbital period and eccentricity ( see table 7 ) .
the inferred total mass of the system is 2.7060.011m. because the system is in the globular cluster m15 ( ngc 7078 ) , it suffers doppler shifts resulting from local accelerations , either by the mean cluster gravitational field or by nearby stars , that are more difficult to estimate than was the case with the galactic system psr 1913 + 16 .
this may make a separate , precision measurement of the relativistic contribution to b impossible .
this binary pulsar system is not particularly relativistic , with a long period ( 12 days ) and highly circular orbit .
however , because we observe the orbit nearly edge on , the shapiro delay is large and measurable , as reflected in the post - keplerian parameters r and s. psr 0655 + 64 .
this system consists of a pulsar and a white dwarf companion in a nearly circular orbit .
only an upper limit on b has been placed . soon after the discovery of the binary pulsar it was widely hailed as a new testing ground for relativistic gravitational effects . as we have seen in the case of gr , in most respects ,
the system has lived up to , indeed exceeded , the early expectations . in another respect , however , the system has only partially lived up to its promise , namely as a direct testing ground for alternative theories of gravity .
the origin of this promise was the discovery that alternative theories of gravity generically predict the emission of dipole gravitational radiation from binary star systems .
in general relativity , there is no dipole radiation because the  dipole moment  ( center of mass ) of isolated systems is uniform in time ( conservation of momentum ) , and because the  inertial mass  that determines the dipole moment is the same as the mass that generates gravitational waves ( sep ) . in other theories , while the inertial dipole moment may remain uniform , the  gravity wave 
dipole moment need not , because the mass that generates gravitational waves depends differently on the internal gravitational binding energy of each body than does the inertial mass ( violation of sep ) .
schematically , in a coordinate system in which the center of inertial mass is at the origin , so that mi,1x1+mi,2x2=0 , the dipole part of the retarded gravitational field would be given by 
 ( 64)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 h \sim \frac{1}{r}\frac{d}{{dt}}\left ( { { m_{{\rm{gw,1}}}}{{\rm{x}}_{\rm{1 } } } + { m_{{\rm{gw,2}}}}{{\rm{x}}_2 } } \right ) \sim \frac{{\eta m}}{r}{\rm{v}}\left ( { \frac{{{m_{{\rm{gw,1}}}}}}{{{m_{{\rm{i}},1 } } } } - \frac{{{m_{{\rm{gw,2}}}}}}{{{m_{{\rm{i}},2 } } } } } \right ) , 
 $ $ \end{document } 
 where v = v1v2 and  and m are defined using inertial masses . in theories that violate sep , the difference between gravitational wave mass and inertial mass is a function of the internal gravitational binding energy of the bodies .
this additional form of gravitational radiation damping could , at least in principle , be significantly stronger than the usual quadrupole damping , because it depends on fewer powers of the orbital velocity v , and it depends on the gravitational binding energy per unit mass of the bodies , which , for neutron stars , could be as large as 40 percent ( see tegp 10   for further details ) .
as one fulfillment of this promise , will and eardley worked out in detail the effects of dipole gravitational radiation in the bimetric theory of rosen , and , when the first observation of the decrease of the orbital period was announced in 1979 , the rosen theory suffered a terminal blow .
a wide class of alternative theories also fails the binary pulsar test because of dipole gravitational radiation ( tegp 12.3 ) . 
  on the other hand ,
the early observations of psr 1913 + 16 already indicated that , in gr , the masses of the two bodies were nearly equal , so that , in theories of gravity that are in some sense  close  to gr , dipole gravitational radiation would not be a strong effect , because of the apparent symmetry of the system .
the rosen theory , and others like it , are not  close  to general relativity , except in their predictions for the weak - field , slow - motion regime of the solar system .
when relativistic neutron stars are present , theories like these can predict strong effects on the motion of the bodies resulting from their internal highly relativistic gravitational structure ( violations of sep ) . as a consequence
, the masses inferred from observations of the periastron shift and  may be significantly different from those inferred using general relativity , and may be different from each other , leading to strong dipole gravitational radiation damping .
by contrast , the brans - dicke theory is  close  to gr , roughly speaking within 1/bd of the predictions of the latter , for large values of the coupling constant bd ( here we use the subscript bd to distinguish the coupling constant from the periastron advance  ) .
thus , despite the presence of dipole gravitational radiation , the binary pulsar provides at present only a weak test of brans - dicke theory , not yet competitive with solar - system tests . making the usual assumption that both members of the system are neutron stars , and using the methods summarized in tegp 1012 , one can obtain formulas for the periastron shift , the gravitational redshift / second - order doppler shift parameter , and the rate of change of orbital period , analogous to eqs .
these formulas depend on the masses of the two neutron stars , on their self - gravitational binding energy , represented by
 sensitivities  s and  * , and on the brans - dicke coupling constant bd .
first , there is a modification of kepler s third law , given by 
 ( 65)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 2\pi { f_{\rm{b } } } = { \left ( { { \mathcal g}m/{a^3 } } \right)^{1/2}}. 
 $ $ \end{document } 
 then , the predictions for  ,  and b are 
 ( 66)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 \left\langle { \dot \omega } \right\rangle = 6\pi { f_{\rm{b}}}{\left ( { 2\pi m{f_{\rm{b } } } } \right)^{2/3}}{\left ( { 1 - { e^2 } } \right)^ { - 1}}{\mathcal p}{{\mathcal g}^ { - 4/3 } } , 
 $ $ \end{document } 
 ( 67)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 \gamma ' = e{\left ( { 2\pi { f_{\rm{b } } } } \right)^ { - 1}}{\left ( { 2\pi m{f_{\rm{b } } } } \right)^{2/3}}\left ( { { m_2}/m } \right){{\mathcal g}^ { - 1/3}}\left ( { \alpha _ 2^ * + { \mathcal g}{m_2}/m + \kappa _ 1^*\eta _ 2^ * } \right ) , 
 $ $ \end{document } 
 ( 68)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 { \dot p_{\rm{b } } } = - \left ( { 192\pi /5 } \right){\left ( { 2\pi { \mathcal m}{f_{\rm{b } } } } \right)^{5/3}}f'\left ( e \right ) - 4\pi \left ( { 2\pi \mu { f_{\rm{b } } } } \right)\xi { { \mathcal s}^2}g\left ( e \right ) , 
 $ $ \end{document } 
 where \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\mathcal m } \equiv { \chi ^{3/5}}{{\mathcal g}^ { - 4/5}}{\eta ^{3/5}}m$$\end{document } , and , to first order in (2+bd ) , we have 
 ( 69)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 \begin{array}{*{20}{c } } 
 { f'\left ( e \right ) = f\left ( e \right ) + \tfrac{5}{144}\xi { { \left ( { \gamma + 3\gamma ' } \right)}^2}\left ( { \tfrac{1}{2}{e^2 } + \tfrac{1}{8}{e^4 } } \right){{\left ( { 1 - { e^2 } } \right)}^ { - 7/2}},}\\ 
 { g\left ( e \right ) = { { \left ( { 1 - { e^2 } } \right)}^ { - 5/2}}\left ( { 1 + \tfrac{1}{2}{e^2 } } \right),}\\ 
 { { \mathcal s } = { s_1 } - { s_2},}\\ 
 { { \mathcal g } = 1 - \xi \left ( { { s_1 } + { s_2 } - 2{s_1}{s_2 } } \right),}\\ 
 { { \mathcal p } = { \mathcal g}\left [ { 1 - \tfrac{2}{3}\xi + \tfrac{1}{3}\xi \left ( { { s_1 } + { s_2 } - 2{s_1}{s_2 } } \right ) } \right],}\\ 
 { \alpha _ 2^ * = 1 - \xi { s_2},}\\ 
 { \eta _ 2^ * = \left ( { 1 - 2{s_2 } } \right)\xi , } \\ 
 { \chi = { { \mathcal g}^2}\left [ { 1 - \tfrac{1}{2}\xi + \tfrac{1}{12}\xi { \gamma ^2 } } \right],}\\ 
 { \gamma = 1 - 2\left ( { { m_1}{s_2 } + { m_2}{s_1 } } \right)/m,}\\ 
 { \gamma ' = 1 - { s_1 } - { s_2}. } 
 \end{array } 
 $ $ \end{document } 
 the quantities sa and a * a are defined by 
 ( 70)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 \begin{array}{*{20}{c } } 
 { { s_{\rm{a } } } = - { { \left ( { \frac{{\partial \left ( { \ln { m_{\rm{a } } } } \right)}}{{\partial \left ( { \ln g } \right ) } } } \right)}_n},}&{\kappa _ { \rm{a}}^ * = - { { \left ( { \frac{{\partial \left ( { \ln { i_{\rm{a } } } } \right)}}{{\partial \left ( { \ln g } \right ) } } } \right)}_n } , } 
 \end{array } 
 $ $ \end{document } 
 and measure the  sensitivity  of the mass ma and moment of inertia ia of each body to changes in the scalar field ( reflected in changes in g ) for a fixed baryon number n ( see tegp 11 , 12 and 14.6 ( c )   for further details ) . the quantity sa is related to the gravitational binding energy . notice how the violation of sep in brans - dicke theory introduces complex structure - dependent effects in everything from the newtonian limit ( modification of the effective coupling constant in kepler s third law ) to gravitational radiation . in the limit 0 , we recover gr , and all structure dependence disappears .
the first term in b ( eq . ( 68 ) ) is the effect of quadrupole and monopole gravitational radiation , while the second term is the effect of dipole radiation . 
  in order to estimate the sensitivities sa and a * , one must adopt an equation of state for the neutron stars .
it is sufficient to restrict attention to relatively stiff neutron star equations of state in order to guarantee neutron stars of sufficient mass , approximately 1.4m. the lower limit on bd required to give consistency among the constraints on  ,  and b as in figure 6 is several hundred .
the combination of  and  give a constraint on the masses that is relatively weakly dependent on  , thus the constraint on  is dominated by b and is directly proportional to the measurement error in b ; in order to achieve a constraint comparable to the solar system value of 310 , the error in bobs would have to be reduced by more than a factor of ten .
alternatively , a binary pulsar system with dissimilar objects , such as a white dwarf or black hole companion , would provide potentially more promising tests of dipole radiation .
unfortunately , none has been discovered to date ; the dissimilar system b0655 + 64 , with a white dwarf companion is in a highly circular orbit , making measurement of the periastron shift meaningless , and is not as relativistic as 1913 + 16 . from the upper limit on b ( table 7 )
damour and esposito - farse   have generalized these results to a broad class of scalar - tensor theories .
these theories are characterized by a single function ( ) of the scalar field  , which mediates the coupling strength of the scalar field . for application to the solar system or to binary systems ,
one expands this function about a cosmological background field value 0 : 
 ( 71)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 \alpha \left ( \varphi \right ) = { \alpha _ 0}\left ( { \varphi - { \varphi _ 0 } } \right ) + \frac{1}{2}{\beta _ 0}{\left ( { \varphi - { \varphi _ 0 } } \right)^2 } + \ldots
. 
 $ $ \end{document } 
 a purely linear coupling function produces brans - dicke theory , with 02=1/(2bd+3 ) .
the function ( ) acts as a potential function for the scalar field  , and , if 0>0 , during cosmological evolution , the scalar field naturally evolves toward the minimum of the potential , i.e. toward 00 , bd , or toward a theory close to , though not precisely gr [ 47 , 48 ] .
bounds on the parameters 0 and 0 from solar - system , binary - pulsar and gravitational wave observations ( see sec .
negative values of 0 correspond to an unstable scalar potential ; in this case , objects such as neutron stars can experience a  spontaneous scalarization  , whereby the interior values of  can take on values very different from the exterior values , through non - linear interactions between strong gravity and the scalar field , dramatically affecting the stars internal structure and the consequent violations of sep .
on the other hand , 0<0 is of little practical interest , because , with an unstable  potential , cosmological evolution would presumably drive the system away from the peak where 00 , toward parameter values that could easily be excluded by solar system experiments . on the 00 plane shown in figure 8 , the 0 axis corresponds to pure brans - dicke theory , while the origin corresponds to pure gr . as discussed above , solar system bounds ( labelled  1pn  in figure 8) still beat the binary pulsars .
figure 8the region of the scalar - tensor theory 00 
 plane allowed by solarsystem , binary - pulsar , and future gravitational wave observations . a polytropic equation of state for the neutron stars was assumed .
the shaded region is that allowed by all tests . for positive values of 
 0 , solar - system bounds ( labelled 1pn )
( from , 1998 by the american physical society , reproduced by permission . ) 
 
the region of the scalar - tensor theory 00 
 plane allowed by solarsystem , binary - pulsar , and future gravitational wave observations .
the shaded region is that allowed by all tests . for positive values of 
 0 , solar - system bounds ( labelled 1pn )
some time in the next decade , a new opportunity for testing relativistic gravity will be realized , with the commissioning and operation of kilometer - scale , laser interferometric gravitational wave observatories in the u.s .
( ligo project ) , europe ( virgo and geo600 projects ) and japan ( tama300 project ) .
the ligo broad - band antennas will have the capability of detecting and measuring the gravitational waveforms from astronomical sources in a frequency band between about 10 hz ( the seismic noise cutoff ) and 500 hz ( the photon counting noise cutoff ) , with a maximum sensitivity to strain at around 100 hz of hl / l10 ( rms ) .
the most promising source for detection and study of the gravitational wave signal is the  inspiralling compact binary   a binary system of neutron stars or black holes ( or one of each ) in the final minutes of a death dance leading to a violent merger .
such is the fate , for example , of the hulse - taylor binary pulsar psr 1913 + 16 in about 300 million years . given the expected sensitivity of the  advanced ligo 
( around 2007 ) , which could see such sources out to hundreds of megaparsecs , it has been estimated that from 3 to 100 annual inspiral events could be detectable .
other sources , such as supernova core collapse events , instabilities in rapidly rotating nascent neutron stars , signals from non - axisymmetric pulsars , and a stochastic background of waves , may be detectable ( for reviews , see [ 1 , 127 ] ; for updates on the status of various projects , see [ 65 , 32 ] ) .
a similar network of cryogenic resonant - mass gravitational antennas have been in operation for many years , albeit at lower levels of sensitivity ( h10 ) .
while modest improvements in sensitivity may be expected in the future , these resonant detectors are not expected to be competitive with the large interferometers , unless new designs involving bars of spherical , or nearly spherical shape come to fruition .
these systems are primarily sensitive to waves in relatively narrow bands about frequencies in the hundreds to thousands of hz range [ 104 , 73 , 14 , 110 ] . in addition
, plans are being developed for an orbiting laser interferometer space antenna ( lisa for short ) .
such a system , consisting of three spacecraft separated by millions of kilometers , would be sensitive primarily in the very low frequency band between 10 and 10 hz , with peak strain sensitivity of order h10 .
in addition to opening a new astronomical window , the detailed observation of gravitational waves by such observatories may provide the means to test general relativistic predictions for the polarization and speed of the waves , and for gravitational radiation damping . a laser interferometric or resonant bar gravitational wave detector measures the local components of a symmetric 33 tensor which is composed of the  electric  components of the riemann curvature tensor , r0i0j , via the equation of geodesic deviation , given for a pair of freely falling particles by  =- r0i0jx , where x denotes the spatial separation . in general there are six independent components , which can be expressed in terms of polarizations ( modes with specific transformation properties under rotations and boosts ) .
three are transverse to the direction of propagation , with two representing quadrupolar deformations and one representing a monopole  breathing  deformation .
three modes are longitudinal , with one an axially symmetric stretching mode in the propagation direction , and one quadrupolar mode in each of the two orthogonal planes containing the propagation direction .
figure 9 shows the displacements induced on a ring of freely falling test particles by each of these modes .
general relativity predicts only the first two transverse quadrupolar modes ( a ) and ( b ) independently of the source ; these correspond to the waveforms h+ and h discussed earlier ( note the cos 2 and sin 2 dependences of the displacements ) .
scalar - tensor gravitational waves can in addition contain the transverse breathing mode ( c ) .
more general metric theories predict additional longitudinal modes , up to the full complement of six ( tegp 10.2 ) . 

figure 9the six polarization modes for gravitational waves permitted in any metric theory of gravity .
there is no displacement out of the plane of the picture . in ( a ) ,
( b ) and ( c ) , the wave propagates out of the plane ; in ( d ) , ( e ) , and ( f ) , the wave propagates in the plane . in general relativity , only ( a ) and ( b ) are present ; in scalar - tensor gravity , ( c ) may also be present . 
 
there is no displacement out of the plane of the picture . in ( a ) , ( b ) and ( c ) , the wave propagates out of the plane ; in ( d ) , ( e ) , and ( f ) , the wave propagates in the plane . in general relativity , only ( a ) and ( b ) are present ; in scalar - tensor gravity , ( c ) may also be present .
a suitable array of gravitational antennas could delineate or limit the number of modes present in a given wave .
there are eight unknowns ( six polarizations and two direction cosines ) , but only six measurables ( r0i0j ) .
if the direction can be established by either association of the waves with optical or other observations , or by time - of - flight measurements between separated detectors , then six suitably oriented detectors suffice to determine all six components . if the direction can not be established , then the system is underdetermined , and no unique solution can be found .
however , if one assumes that only transverse waves are present , then there are only three unknowns if the source direction is known , or five unknowns otherwise .
then the corresponding number ( three or five ) of detectors can determine the polarization .
if distinct evidence were found of any mode other than the two transverse quadrupolar modes of gr , the result would be disastrous for gr .
on the other hand , the absence of a breathing mode would not necessarily rule out scalar - tensor gravity , because the strength of that mode depends on the nature of the source .
some of the details of implementing such polarization observations have been worked out for arrays of resonant cylindrical , disk - shaped , spherical and truncated icosahedral detectors ( tegp 10.2 , for recent reviews see [ 87 , 133 ] ) ; initial work has been done to assess whether the ground - based or space - based laser interferometers ( or combinations of the two types ) could perform interesting polarization measurements [ 134 , 33 , 90 , 67 ] .
unfortunately for this purpose , the two ligo observatories ( in washington and louisiana states , respectively ) have been constructed to have their respective arms as parallel as possible , apart from the curvature of the earth ; while this maximizes the joint sensitivity of the two detectors to gravitational waves , it minimizes their ability to detect two modes of polarization . in the binary pulsar ,
a test of gr was made possible by measuring at least three relativistic effects that depended upon only two unknown masses .
the evolution of the orbital phase under the damping effect of gravitational radiation played a crucial role .
another situation in which measurement of orbital phase can lead to tests of gr is that of the inspiralling compact binary system .
the key differences are that here gravitational radiation itself is the detected signal , rather than radio pulses , and the phase evolution alone carries all the information . in the binary pulsar , the first derivative of the binary frequency b , was measured ; here the full nonlinear variation of fb as a function of time is measured .
broad - band laser interferometers are especially sensitive to the phase evolution of the gravitational waves , which carry the information about the orbital phase evolution .
the analysis of gravitational wave data from such sources will involve some form of matched filtering of the noisy detector output against an ensemble of theoretical  template  waveforms which depend on the intrinsic parameters of the inspiralling binary , such as the component masses , spins , and so on , and on its inspiral evolution .
how accurate must a template be in order to  match  the waveform from a given source ( where by a match we mean maximizing the cross - correlation or the signal - to - noise ratio ) ? in the total accumulated phase of the wave detected in the sensitive bandwidth , the template must match the signal to a fraction of a cycle . for two inspiralling neutron stars , around
16,000 cycles should be detected during the final few minutes of inspiral ; this implies a phasing accuracy of 10 or better . since v1/10 during the late inspiral , this means that correction terms in the phasing at the level of v or higher are needed .
more formal analyses confirm this intuition [ 35 , 60 , 36 , 108 ] . because it is a slow - motion system ( v10 ) , the binary pulsar is sensitive only to the lowest - order effects of gravitational radiation as predicted by the quadrupole formula .
nevertheless , the first correction terms of order v and v to the quadrupole formula were calculated as early as 1976   see tegp 10.3 ) .
but for laser interferometric observations of gravitational waves , the bottom line is that , in order to measure the astrophysical parameters of the source and to test the properties of the gravitational waves , it is necessary to derive the gravitational waveform and the resulting radiation back - reaction on the orbit phasing at least to 2pn order beyond the quadrupole approximation , and probably to 3pn order .
for the special case of non - spinning bodies moving on quasi - circular orbits ( i.e. circular apart from a slow inspiral ) , the evolution of the gravitational wave frequency f=2fb through 2pn order has the form 
 ( 72)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 \dot f = \frac{{96\pi } } { 5}{f^2}{\left ( { \pi { \mathcal m}f } \right)^{5/3}}\left [ { 1 - \left ( { \frac{{743}}{{336 } } + \frac{{11}}{4}\eta } \right){{\left ( { \pi mf } \right)}^{2/3 } } + 4\pi \left ( { \pi mf } \right ) + \left ( { \frac{{34103}}{{18144 } } + \frac{{13661}}{{2016}}\eta + \frac{{59}}{{18}}{\eta ^2 } } \right){{\left ( { \pi mf } \right)}^{4/3 } } + o\left [ { { { \left ( { \pi mf } \right)}^{5/3 } } } \right ] } \right ] , 
 $ $ \end{document } 
 where =m1m2/m2 . the first term is the quadrupole contribution ( compare eq . ( 55 ) ) , the second term is the 1pn contribution , the third term , with the coefficient 4 , is the  tail  contribution , and the fourth term is the 2pn contribution , first reported jointly by blanchet et al .
these losses react back on the orbit to circularize it and cause it to inspiral .
the result is that the orbital phase ( and consequently the gravitational wave phase ) evolves non - linearly with time .
it is the sensitivity of the broad - band ligo and virgo - type detectors to phase that makes the higher - order contributions to df / dt so observationally relevant .
a ready - to - use set of formulae for the 2pn gravitational waveform template , including the non - linear evolution of the gravitational wave frequency ( not including spin effects ) have been published   and incorporated into the gravitational radiation analysis and simulation package ( grasp ) , a software toolkit used in ligo .
if the coefficients of each of the powers of f in eq . ( 72 ) can be measured , then one again obtains more than two constraints on the two unknowns m1 and m2 , leading to the possibility to test gr . for example , blanchet and sathyaprakash [ 30 , 29 ] have shown that , by observing a source with a sufficiently strong signal , an interesting test of the 4 coefficient of the  tail  term could be performed .
another possibility involves gravitational waves from a small mass orbiting and inspiralling into a ( possibly supermassive ) spinning black hole . a general non - circular , non - equatorial orbit will precess around the hole , both in periastron and in orbital plane , leading to a complex gravitational waveform that carries information about the non - spherical , strong - field spacetime around the hole . according to gr
, this spacetime must be the kerr spacetime of a rotating black hole , uniquely specified by its mass and angular momentum , and consequently , observation of the waves could test this fundamental hypothesis of gr [ 114 , 107 ] .
thirdly , the dipole gravitational radiation predicted by scalar - tensor theories will result in a modification of the gravitational radiation back - reaction , and thereby of the phase evolution .
including only the leading quadrupole and dipole contributions , one obtains , in brans - dicke theory , 
 ( 73)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 \dot f = \frac{{96\pi } } { 5}{f^2}{\left ( { \pi { \mathcal m}f } \right)^{5/3}}\left [ { 1 + b{{\left ( { \pi mf } \right)}^ { - 2/3 } } } \right ] , 
 $ $ \end{document } 
 where \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\mathcal m } = \left ( { { \chi ^{3/5}}{{\mathcal g}^ { - 4/5 } } } \right){\eta ^{3/5}}m$$\end{document } , and b is the coefficient of the dipole term , given by \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$b = \left ( { 5/48 } \right)\left ( { { \chi ^ { - 1}}{{\cal g}^{4/3 } } } \right)\xi { { \cal s}^2}$$\end{document } , where  , g , s are given by eqs . ( 69 ) , and =1/(2+bd ) .
double neutron star systems are less promising because the small range of masses available near 1.4m results in suppression of dipole radiation by symmetry ( the sensitivity s turns out to be a relatively weak function of mass near 1.4m , for typical equations of state ) . for black holes , s=0.5 identically ,
consequently double black - hole systems turn out to be observationally identical in the two theories . 
 
but for a 1.4m neutron star and a 10 m ( 3 m ) black hole at 200 mpc , the bound on bd could be 600 ( 1800 ) ( using advanced ligo noise curves ) .
if one demands that this test be performed annually , thus requiring observation of frequent , and therefore more distant , weaker sources , the bounds on bd will be too weak to compete with existing solar - system bounds ( this corresponds to the  ligo - virgo  bound on the 0 axis in figure 8 , which assumes a signal - to - noise ratio of 10 ) .
however , if one is prepared to wait 10 years for the lucky observation of a nearby , strong source , the resulting bound could exceed the current solar - system bound .
figure 10 assumes a double - neutron - star inspiral rate of 10 per year per galaxy , and a black - hole - neutron - star rate  times that , where  is highly uncertain .
for the general class of scalar - tensor theories , the corresponding bounds are plotted on the 00 plane in figure 8 , under the restrictive assumption of a signal - to - noise ratio ( s / n ) of 10 .
all other factors being equal , the bound achievable on the 00 parameters of scalar - tensor gravity is inversely proportional to s / n . 

figure 10the bounds on the scalar - tensor coupling constant bd 
 from gravitational wave observations of inspiralling black - hole - neutron - star systems .
the solar system bound is around bd=3500 . 
  the bounds on the scalar - tensor coupling constant bd 
 from gravitational wave observations of inspiralling black - hole - neutron - star systems .
in which the wavelength of gravitational waves is small compared to the radius of curvature of the background spacetime , the waves propagate along null geodesics of the background spacetime , i.e .
they have the same speed c as light ( in this section , we do not set c=1 ) .
in other theories , the speed could differ from c because of coupling of gravitation to  background  gravitational fields .
for example , in the rosen bimetric theory with a flat background metric  , gravitational waves follow null geodesics of  , while light follows null geodesics of g ( tegp 10.1 ) .
another way in which the speed of gravitational waves could differ from c is if gravitation were propagated by a massive field ( a massive graviton ) , in which case vg would be given by , in a local inertial frame , 
 ( 74)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 \frac{{v_{\rm{g}}^2}}{{{c^2 } } } = 1 - \frac{{m_{\rm{g}}^2{c^4}}}{{{e^2 } } } , 
 $ $ \end{document } 
 where mg and e are the graviton rest mass and energy , respectively . for a recent review of the idea of a massive graviton along with a model theory , see . 
  the most obvious way to test
this is to compare the arrival times of a gravitational wave and an electromagnetic wave from the same event , e.g. a supernova . for a source at a distance d ,
the resulting value of the difference 1vg / c is 
 ( 75)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 1 - \frac{{{v_{\rm{g}}}}}{c } = 5 \times { 10^ { - 17}}\left ( { \frac{{200{\rm { mpc}}}}{d } } \right)\left ( { \frac{{\delta t}}{{1{\rm { s } } } } } \right ) , 
 $ $ \end{document } 
 where tta(1+z)te is the  time difference  , where ta and te are the differences in arrival time and emission time , respectively , of the two signals , and z is the redshift of the source . in many cases ,
te is unknown , so that the best one can do is employ an upper bound on te based on observation or modelling .
the result will then be a bound on 1vg / c . 
  for a massive graviton , if the frequency of the gravitational waves is such that hfmgc , where h is planck s constant , then vg / c11/2(c/gf ) , where g = h / mgc is the graviton compton wavelength , and the bound on 1vg / c can be converted to a bound on g , given by 
 76\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 { \lambda _ { \rm{g } } } > 3 \times { 10^{12}}{\rm{km}}{\left ( { \frac{d}{{200{\rm { mpc}}}}\frac{{100{\rm { hz}}}}{f } } \right)^{1/2}}{\left ( { \frac{1}{{f\delta t } } } \right)^{1/2}}. 
 $ $ \end{document } 
 
the foregoing discussion assumes that the source emits both gravitational and electromagnetic radiation in detectable amounts , and that the relative time of emission can be established to sufficient accuracy , or can be shown to be sufficiently small . however , there is a situation in which a bound on the graviton mass can be set using gravitational radiation alone .
that is the case of the inspiralling compact binary . because the frequency of the gravitational radiation sweeps from low frequency at the initial moment of observation to higher frequency at the final moment , the speed of the gravitons emitted will vary , from lower speeds initially to higher speeds ( closer to c ) at the end
this will cause a distortion of the observed phasing of the waves and result in a shorter than expected overall time ta of passage of a given number of cycles .
furthermore , through the technique of matched filtering , the parameters of the compact binary can be measured accurately , ( assuming that gr is a good approximation to the orbital evolution , even in the presence of a massive graviton ) , and thereby the emission time te can be determined accurately .
( 76 ) can be measured to an accuracy 1/ , where  is the signal - to - noise ratio .
thus one can estimate the bounds on g achievable for various compact inspiral systems , and for various detectors . for stellar - mass inspiral ( neutron stars or black holes )
observed by the ligo / virgo class of ground - based interferometers , d200 mpc , f100 hz , and ft1/10 .
the result is g>10 km . for supermassive binary black holes ( 10 to 10m )
observed by the proposed laser interferometer space antenna ( lisa ) , d3 gpc , f10 hz , and ft1/1000 .
a full noise analysis using proposed noise curves for the advanced ligo and for lisa weakens these crude bounds by factors between two and 10 .
these potential bounds can be compared with the solid bound g>2.810 km ,   derived from solar system dynamics , which limit the presence of a yukawa modification of newtonian gravity of the form 
 ( 77)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 v\left ( r \right ) = \left ( { gm / r } \right)\exp \left ( { - r/{\lambda _ { \rm{g } } } } \right ) , 
 $ $ \end{document } 
 and with the model - dependent bound g>610 km from consideration of galactic and cluster dynamics . 
 
one of the central difficulties of testing general relativity in the strong - field regime is the possibility of contamination by uncertain or complex physics . in the solar system , weak - field gravitational effects could in most cases be measured cleanly and separately from non - gravitational effects .
the remarkable cleanliness of the binary pulsar permitted precise measurements of gravitational phenomena in a strong - field context .
still , under suitable conditions , qualitative and even quantitative strong - field tests of general relativity can be carried out .
one example is in cosmology . from a few seconds after the big bang until the present
, the underlying physics of the universe is well understood , although significant uncertainties remain ( amount of dark matter , value of the cosmological constant , the number of light neutrino families , etc . ) .
some alternative theories of gravity that are qualitatively different from gr fail to produce cosmologies that meet even the minimum requirements of agreeing qualitatively with big - bang nucleosynthesis ( bbn ) or the properties of the cosmic microwave background ( tegp 13.2 ) .
others , such as brans - dicke theory , are sufficiently close to gr ( for large enough bd ) that they conform to all cosmological observations , given the underlying uncertainties .
the generalized scalar - tensor theories , however , could have small bd at early times , while evolving through the attractor mechanism to large bd today .
one way to test such theories is through big - bang nucleosynthesis , since the abundances of the light elements produced when the temperature of the universe was about 1 mev are sensitive to the rate of expansion at that epoch , which in turn depends on the strength of interaction between geometry and the scalar field . because the universe is radiation - dominated at that epoch , uncertainties in the amount of cold dark matter or of the cosmological constant are unimportant .
the nuclear reaction rates are reasonably well understood from laboratory experiments and theory , and the number of light neutrino families ( 3 ) conforms to evidence from particle accelerators .
thus , within modest uncertainties , one can assess the quantitative difference between the bbn predictions of gr and scalar - tensor gravity under strong - field conditions and compare with observations .
the most sophisticated recent analysis   places bounds on the parameters 0 and 0 of the generalized framework of damour and esposito - farse ( see sec . 5.4 and fig . 8) that are weaker than solar - system bounds for 0<0.3 , but substantially stronger for >0.3 .
observations of low - luminosity binary x - ray sources suggest that a form of accretion known as advection - dominated accretion flow ( adaf ) may be important . in this kind of flow ,
the accreting gas is too thin to radiate its energy efficiently , but instead transports ( advects ) it inward toward the central object .
if the central object is a neutron star , the matter hits the surface and radiates the energy away ; if it is a black hole , the matter and its advected energy disappear .
systems in which the accreting object is believed to be a black hole from estimates of its mass are indeed observed to be underluminous , compared to systems where the object is believe to be a neutron star .
this has been regarded as the first astrophysical evidence for the existence of black hole event horizons ( for a review , see ) . while supporting one of the critical strong - field predictions of gr , the observations and models are not likely any time soon to be able to distinguish one gravitational theory from another ( except for theories that do not predict black holes at all ) .
another example involving accretion purports to explore the strong - field region just outside massive black holes in active galactic nuclei . here ,
iron in the inner region of a thin accretion disk is irradiated by x - ray - emitting material above or below the disk , and fluoresces in the k line .
the spectral shape of the line depends on relativistic doppler and curved - spacetime effects as the iron orbits the black hole near the innermost stable circular orbit , and could be used to determine whether the hole is a non - rotating schwarzschild black hole , or a rotating kerr black hole .
because of uncertainties in the detailed models , the results are inconclusive to date , but the combination of higher - resolution observations and better modelling could lead to striking tests of strong - field predictions of gr .
one reason is that gravity is a fundamental interaction of nature , and as such requires the most solid empirical underpinning we can provide .
another is that all attempts to quantize gravity and to unify it with the other forces suggest that the standard general relativity of einstein is not likely to be the last word .
furthermore , the predictions of general relativity are fixed ; the theory contains no adjustable constants so nothing can be changed .
thus every test of the theory is either a potentially deadly test or a possible probe for new physics .
although it is remarkable that this theory , born 80 years ago out of almost pure thought , has managed to survive every test , the possibility of finding a discrepancy will continue to drive experiments for years to come .
 the status of experimental tests of general relativity and of theoretical frameworks for analysing them are reviewed . 
 einstein s equivalence principle ( eep ) is well supported by experiments such as the etvs experiment , tests of special relativity , and the gravitational redshift experiment . 
 future tests of eep and of the inverse square law will search for new interactions arising from unification or quantum gravity . 
 tests of general relativity at the post - newtonian level have reached high precision , including the light defl 
 ection the shapiro time delay , the perihelion advance of mercury , and the nordtvedt effect in lunar motion . 
 gravitational wave damping has been detected in an amount that agrees with general relativity to half a percent using the hulse - taylor binary pulsar , and new binary pulsar systems may yield further improvements . when direct observation of gravitational radiation from astrophysical sources begins , new tests of general relativity will be possible . 
Introduction
Tests of the Foundations of Gravitation Theory
Tests of Post-Newtonian Gravity
Strong Gravity and Gravitational Waves: A New Testing Ground
Stellar System Tests of Gravitational Theory
Gravitational Wave Tests of Gravitational Theory
Conclusions