arxiv: v1 [physics.gen-ph] 29 Feb 2016

Similar documents
arxiv: v1 [physics.gen-ph] 10 May 2017

The Precession of Mercury s Perihelion

KEPLER S LAWS OF PLANETARY MOTION

Newton s Laws, Kepler s Laws, and Planetary Orbits

A New Approach to General Relativity

Geometry of the homogeneous and isotropic spaces

d 2 x 0a d d =0. Relative to an arbitrary (accelerating frame) specified by x a = x a (x 0b ), the latter becomes: d 2 x a d 2 + a dx b dx c

Classical Mechanics Homework set 7, due Nov 8th: Solutions

The Schwartzchild Geometry

Chapter 13 Gravitation

OSCILLATIONS AND GRAVITATION

Tutorial Exercises: Central Forces

GENERAL RELATIVITY: THE GEODESICS OF THE SCHWARZSCHILD METRIC

MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department. Problem Set 10 Solutions. r s

Paths of planet Mars in sky

Appendix B The Relativistic Transformation of Forces

DEVIL PHYSICS THE BADDEST CLASS ON CAMPUS IB PHYSICS

Solving Problems of Advance of Mercury s Perihelion and Deflection of. Photon Around the Sun with New Newton s Formula of Gravity

= 4 3 π( m) 3 (5480 kg m 3 ) = kg.

Physics 235 Chapter 5. Chapter 5 Gravitation

AST 121S: The origin and evolution of the Universe. Introduction to Mathematical Handout 1

Determining solar characteristics using planetary data

Physics: Work & Energy Beyond Earth Guided Inquiry

MODULE 5 ADVANCED MECHANICS GRAVITATIONAL FIELD: MOTION OF PLANETS AND SATELLITES VISUAL PHYSICS ONLINE

Earth and Moon orbital anomalies

University Physics Volume I Unit 1: Mechanics Chapter 13: Gravitation Conceptual Questions

Physics 161: Black Holes: Lecture 5: 22 Jan 2013

KEPLER S LAWS AND PLANETARY ORBITS

Physics 506 Winter 2006 Homework Assignment #9 Solutions

Homework 7 Solutions

Central Force Problem. Central Force Motion. Two Body Problem: Center of Mass Coordinates. Reduction of Two Body Problem 8.01 W14D1. + m 2. m 2.

Mechanics Physics 151

Central Force Motion

Homework # 3 Solution Key

Between any two masses, there exists a mutual attractive force.

Projection Gravitation, a Projection Force from 5-dimensional Space-time into 4-dimensional Space-time

Physics 161 Fall 2011 Extra Credit 2 Investigating Black Holes - Solutions The Following is Worth 50 Points!!!

Universal Gravitation

AY 7A - Fall 2010 Section Worksheet 2 - Solutions Energy and Kepler s Law

Experiment 09: Angular momentum

Light Time Delay and Apparent Position

F(r) = r f (r) 4.8. Central forces The most interesting problems in classical mechanics are about central forces.

Physics 181. Assignment 4

AH Mechanics Checklist (Unit 2) AH Mechanics Checklist (Unit 2) Circular Motion

1 Basics of two-body orbits

From Gravitational Collapse to Black Holes

arxiv: v2 [gr-qc] 18 Aug 2014

Reversed Gravitational Acceleration for High-speed Particles

F 12. = G m m 1 2 F 21 = F 12. = G m 1m 2. Review. Physics 201, Lecture 22. Newton s Law Of Universal Gravitation

Chapter 12. Kinetics of Particles: Newton s Second Law

DEVIL PHYSICS THE BADDEST CLASS ON CAMPUS IB PHYSICS

Gravitation. AP/Honors Physics 1 Mr. Velazquez

Chapter 5 Force and Motion

Chapter 5 Force and Motion

Introduction to General Relativity 2

Scattering in Three Dimensions

DEVIL PHYSICS THE BADDEST CLASS ON CAMPUS IB PHYSICS

The Concept of the Effective Mass Tensor in GR. Clocks and Rods

SIO 229 Gravity and Geomagnetism. Lecture 6. J 2 for Earth. J 2 in the solar system. A first look at the geoid.

= e2. = 2e2. = 3e2. V = Ze2. where Z is the atomic numnber. Thus, we take as the Hamiltonian for a hydrogenic. H = p2 r. (19.4)

HW Solutions # MIT - Prof. Please study example 12.5 "from the earth to the moon". 2GmA v esc

Chapter. s r. check whether your calculator is in all other parts of the body. When a rigid body rotates through a given angle, all

Math Notes on Kepler s first law 1. r(t) kp(t)

Galilean Transformation vs E&M y. Historical Perspective. Chapter 2 Lecture 2 PHYS Special Relativity. Sep. 1, y K K O.

10. Universal Gravitation

arxiv: v1 [physics.pop-ph] 3 Jun 2013

MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department. Problem Set 9

ASTR415: Problem Set #6

The Schwarzschild Solution

B. Spherical Wave Propagation

Lecture 22. PE = GMm r TE = GMm 2a. T 2 = 4π 2 GM. Main points of today s lecture: Gravitational potential energy: Total energy of orbit:

History of Astronomy - Part II. Tycho Brahe - An Observer. Johannes Kepler - A Theorist

Hilbert s forgotten equation of velocity dependent acceleration in a gravitational field

Derivation of the Gravitational Red Shift from the Theorem of Orbits

10. Force is inversely proportional to distance between the centers squared. R 4 = F 16 E 11.

Circular Orbits. and g =

15 Solving the Laplace equation by Fourier method

Physics 207 Lecture 5. Lecture 5

Physics 121 Hour Exam #5 Solution

Problems with Mannheim s conformal gravity program

Conformal transformations + Schwarzschild

arxiv:hep-th/ v11 4 Feb 2016

Is there a magnification paradox in gravitational lensing?

Physics 107 TUTORIAL ASSIGNMENT #8

Chap 5. Circular Motion: Gravitation

Radial Inflow Experiment:GFD III

Revision Guide for Chapter 11

Physics 111. Ch 12: Gravity. Newton s Universal Gravity. R - hat. the equation. = Gm 1 m 2. F g 2 1. ˆr 2 1. Gravity G =

Physics 2A Chapter 10 - Moment of Inertia Fall 2018

Ch 13 Universal Gravitation

Lecture 8 - Gauss s Law

Chem 453/544 Fall /08/03. Exam #1 Solutions

Problems with Mannheim s conformal gravity program

Three dimensional flow analysis in Axial Flow Compressors

arxiv:gr-qc/ v2 8 Jun 2006

Chapter 7. Rotational Motion Angles, Angular Velocity and Angular Acceleration Universal Law of Gravitation Kepler s Laws

20th Century Atomic Theory - Hydrogen Atom

- 5 - TEST 1R. This is the repeat version of TEST 1, which was held during Session.

working pages for Paul Richards class notes; do not copy or circulate without permission from PGR 2004/11/3 10:50

Gravitation. Chapter 12. PowerPoint Lectures for University Physics, Twelfth Edition Hugh D. Young and Roger A. Freedman. Lectures by James Pazun

Transcription:

epl daft Pedicting Mecuy s Pecession using Simple Relativistic Newtonian Dynamics axiv:603.0560v [physics.gen-ph] 9 Feb 06 Y. Fiedman and J. M. Steine Jeusalem College of Technology Jeusalem, Isael PACS 95.30.Sf Relativity and gavitation PACS 95.0.Eg Obit detemination and impovement Abstact We pesent a new simple elativistic model fo planetay motion descibing accuately the anomalous pecession of the peihelion of Mecuy and its oigin. The model is based on tansfoming Newton s classical equation fo planetay motion fom absolute to eal spacetime influenced by the gavitational potential and intoducing the concept of influenced diection. Intoduction. The centennial of the discovey of Geneal Relativity GR) by Albet Einstein is an oppotune moment to stive to fully undestand his ideas on gavitation. Although Einstein s GR oiginates in the logical incompleteness of Newton s Gavitation NG), that incompleteness itself has not been undestood completely, as yet. In ou quest to unavel this poblem, we incopoated the basic idea of tansition fom absolute space and time to influenceable spacetime fo the poblem of planetay motion in NG. In 687, Si Isaac Newton descibed gavity as an instantaneous and invisible foce between two objects. At the end of the 9th centuy, thee seemed to be no dynamic poblem which could not be addessed with classical Newtonian mechanics except the anomalous pecession of Mecuy s obit. In 95, exactly 00 yeas ago, Albet Einstein devised a completely new desciption of gavity. In GR, the gavitational foce is a fictitious foce due to the cuvatue of spacetime. This new theoy enabled Einstein to calculate the obseved anomalous pecession of Mecuy fist ecognized in 859 by Ubain Le Veie. In this lette we pesent a new simple elativistic model fo planetay motion pedicting Mecuy s pecession without GR. The enegy consevation equation fo planetay motion in NG is ewitten in tems of dimensionless enegy and, then, into a nom equation fo the 4-velocity in absolute spacetime. This nom equation is then tansfomed into the coesponding equation in a spacetime influenced by the gavitational potential - the eal spacetime. Intoducing the concept of influenced diection, the esulting equation yields immediately the known equation fo planetay motion pedicted by GR. This model pedicts the obseved value and also povides an intepetation of Mecuy s anomalous pecession fom the point of view of elativistic Newtonian dynamics. Finally, we show how to ecove the Schwazschild metic fom this equation. Planetay motion in NG using the concept of dimensionless enegy. We begin with a bief eview of planetay motion in NG pesented in the liteatuesee fo example [] - [4]) with the intoduction of the concept of dimensionless enegy. To descibe the motion of a planet of mass m unde the gavity of the Sun of mass M, we assume that m M and choose a coodinate system with the oigin at the cente of the Sun. Denote by the position vecto of the cente of the planet elative to the Sun and by G the Newtonian gavity constant. The motion of the planet is govened by Newton s second law m d t) dt = Ut)), ) whee the potential enegy of the gavitational field U = GMm and U denotes the gadient of U. This yields the enegy consevation equation in NG m ) d GMm = E ) dt expessing that the total enegy E the sum of the kinetic and potential enegies) of the obit is conseved. Fo any enegy we define a dimensionless enegy as its atio to the maximum kinetic enegy MKE) of the planet p-

Y. Fiedman J. M. Steine mc, whee c is the speed of light. The dimensionless kinetic enegy DKE) is d c dt) = β, whee β is the known beta-facto. Fo any point on the obit, the potential enegy is negative and we denote by u the absolute value of the dimensionless potential enegy DPE) u = GM c = s, s = GM c, 3) whee s is the Schwazschildadius ofthe Sun the adius ofaspheesuchthat, ifallthe massweetobe compessed within a sphee with this adius, the escape velocity fom the suface of the sphee would equal the speed of light). Finally, we denote by E = E mc the dimensionless total enegy DTE) of the obit. Dividing equation ) by the MKE of the planet we obtain the consevation equation of the DTE on the obit c ) d s dt = E. 4) We intoduce pola coodinates,ϕ in the plane of the obit, whee is the distance of the planet fom the Sun and ϕ is the dimensionless pola angle, measued in adians. Consevation of angula momentum pe unit mass J, allowsus toexpessthe angulavelocityas dϕ dt = J and to decompose the squae of the velocity of the planet as the sum of the squaes of its othogonal adial and tansvese components ) d = dt ) d + J dt. 5) Substituting this into 4) we obtain the classical dimensionless enegy consevation equation ) d c + J dt c s = E. 6) Using the definition 3) of u, and denoting its deivative with espect to ϕ by u, it can be shown that Hence, equation 6) becomes d dt = J s u. 7) J sc u ) +u ) = u+e. 8) Multiplying this equation by µ, whee µ is a unit-fee obit paamete µ = sc J 9) we obtain u ) = u +µu+µe. 0) Diffeentiating this equation with espect to ϕ and dividing by u we obtain a linea diffeential equation with constant coefficients u +u = µ. ) Its solution is uϕ) = µ+εcosϕ ϕ 0 )), ) whee ε - the eccenticity of the obit, and ϕ 0 - the pola angle of the peihelion. This implies that ϕ) = s /µ +εcosϕ ϕ 0 ) 3) and the obit is a non-pecessing ellipse. Since the minima of ϕ), coesponding to the peihelion, occu when ϕ = ϕ 0 +πn, n = 0,,, the position of the peihelion will not change with the evolution of the planet. Moeove, fom equation ) we obtain that µ = π π 0 uϕ)dϕ 4) is the absolute value of the angula aveage DPE on the obit. The peihelion p and aphelion a of the obit coespond to d dt = 0 which by 7) equies u = 0 fo thei coesponding DPE. Hence, the u values of the peihelion u p and aphelion u a ae the two positive) oots of the quadatic on the ight hand side of 0). This equies E, the DTE, to be negative fo a bound obit, and fom ) µ = u p +u a = s L, 5) whee L is the semi-latus ectum of the obit. The Schwazschild adius of the Sun is s = 953.5m, fo Mecuy a = 6.9868 0 0 m and p = 4.600 0 0 m implying that µ = 5.3497 0 8 and ε = 0.0563 as obseved). Note that µ. A new simple elativistic model fo planetay motion. In the pevious section we have seen that NG pedicts a non-pecessing obit fo planetay motion. Special Relativity SR) pedicts a pecession of Mecuy, but much lowethen the obsevedone see, fo example [5] p.033). The coect pediction of the obseved anomalous pecession is povided by Einstein s GR theoy of gavitation based on cuved spacetime. In this lette we pesent a much simple model which also pedicts the coect pecession. Themotionofaplanetcanbedecomposedintotwopeiodic motions: the adial motion and the angula otation. In NG, the peiods of these motions ae equal, esulting in a non-pecessing obit. Since in eality thee is a pecession, these two peiods ae not equal. The eason fo this lies in the inaccuate desciption of the espective velocities adial and tansvese) of these motions by the NG model. Einsten s SR and GR assume that in the classical limit when the velocities DKE ) o gavitational fields o acceleation) DP E ) ae small espectively, the laws of dynamics educe to the classical Newton s laws. Thus, these theoies diffe only in the way how events ae tansfomed fom a local fame attached to the moving p-

Pedicting Mecuy s Pecession object) to a fixed inetial efeence fame. Fo instance, see, fo example[6]) the 3D elativistic dynamics equation in SR can be deived by tansfoming Newton s dynamic law fom an inetial comoving fame to the efeence fame by the use of the Loentz tansfomations. These spacetime tansfomations ae not tivial. Fom Planck s fomula see, fo example [] p.0) it is evident that time is influenced by enegy. SR consides the tansfomation between inetial fames with no potential enegy). Any est object in the moving fame has a non-zeo kinetic enegy in the efeence fame. The Relativity Pinciple [6] implies that in this tansfomation space as well as time ae both influenced by this kinetic enegy. The influence on space is in some given diection, which we call the influenced diection. By the Loentz tansfomations the influenced diection in SR is the diection of the velocity, which entails that time and the component of the spatial) displacement vecto in this diection ae alteed by the Loentz time dilation γ factodepending on DKE), wheeas the spatial components tansvese to this diection, ae not. The same is tue fo the tansfomations of the components of the 4-velocity. Since in planetay motion both adial and tansvese components of the velocity ae alteed, the SR model pedicts a pecession, but with a value significantly lowe than the obseved one. Anothe eason fo this low pecession is that the influence of the potential enegy on spacetime is ignoed by the SR model. We pesent a new simple elativistic model fo planetay motion by tansfoming the NG solution fom an absolute flat) spacetime to the eal spacetime influenced by the gavitational potential, hence by the DPE. The gavitational time dilation facto descibing the effect of the DPE on time in such spacetime is see fo example []) s /, 6) the analogof the Loentz γ facto of SR). Note that γ has the same fom as γ with β eplaced by s /, which also eflects the fact that spacetime in SR is influenced by the DKE, while in ou model it is influenced by the DPE. It is knownthat such γ yields the tue time dilation in the classical limit, and its coectness was veified expeimentally in gavitational ed-shift expeiments. The gavitational time dilation facto γ can be undestood communicated by Z. Weinbege) in tems of the escape velocity v e, defined as the minimum speed needed fo an object to beak fee fom the gavitational attaction of a massive body. Moe paticulaly, it is the velocity speed tavelled away fom the stating point) at which the sum of the object s kinetic and its gavitational potential enegies is equal to zeo. This is also the velocity of the object pulled by the massive body fom infinity. The object keeps ecod of gavity in the fom of the escape velocity, which gets stoed in it in the fom of negative) potential enegy. Fo a spheically symmetic massive body, the escape velocity at a given distance is given by In tems of v e GM v e =. 7) v e /c, which is analogous to the γ in SR. Altenatively, we can view γ as epesenting the influence of the acceleation on time. Fom the definition of the Schwazschild adius it follows that any eal tajectoy must satisfy > s in ode to pevent it to be absobed into the Sun. The magnitude of the classical) acceleation of a fee moving object in the gavitational field of the Sun is a = GM/. Hence, by use of 3), s < s a a = GM/a GM/s = c / s ) <, 8) implying the existence of a maximal acceleation a m in ou sola system and hence a m = c s =.564 0 3 m/s 9) a/a m. 0) The influence of the acceleation on time is discussed in [8] and efeences theein. We ewite the classical NG dimensionless enegy consevation equation 6) as s = c ) d J dt c +E +. ) Dividing this equation by the left-hand side we obtain +E s d ) J c dt c s s =. ) This equation is just anothe desciption fo the 3D Newtonian obit with the DPE instead appeaing as a sepaate tem, it modifies all the othe tems. If we lift this obit in 4D absolute spacetime, this equation esembles a nom equation fo the time, adial and tansvese components of a 4-velocity, all multiplied by γ. In this case thee ae only 3 components, because the plana obit has only two non-zeo space components. The multiplication by γ does not effect the obit, it meely defines an ac-length paametization of the classical obit. Since the dynamic equation ) involves only the fist deivative of the potential U) at any given point, it is sufficient to conside the tansfomation fom an absolute to a eal spacetime influenced by a linea gavitational potential, the linea pat in the expansion of U at this point. By the Equivalence Pinciple, such a spacetime is p-3

Y. Fiedman J. M. Steine equivalent to a unifomly acceleated system with acceleation a = U/m. Using the Genealised Pinciple of Relativity, it was shown in [7] and [8] that, unde this tansfomation, both time and the component of velocity in the diection of a ae alteed by the same γ. Thus, in ou poposed model the influenced diection is the adial diection. Hence, the tansfomation fom the absolute to the eal spacetime should alte by the time dilation facto γ) the time and only the adial component of the 4-velocity and not the tansvese ones. With these ideas we tansfom the 4-velocity nom equation ) descibing the motion in absolute spacetime to the coesponding nom equation in eal spacetime. The fist tem in equation is γ, the squae of the time dilation facto scaled by the DTE of the obit. The second tem is the squae of the 4-velocity component in the influenced adial) diection multiplied by γ with the aclength paametization equiing that the unit-fee adial ) beeplacedwith d ds). Thethidtemis velocity d c dt the squae of the tansvese component of the 4-velocity also multiplied by γ. Howeve, this tem epesents a component tansvese to the influenced diection not influenced by the DPE), which should not be affected by ou tansfomation. Hence,the coefficient γ should be omitted fom this tem. This yields ou modified equation fo planetay motion implying +E d ) s / ds s / J c = 3) ) d + J ds c s ) s = E. 4) A close inspection eveals that this equation is simply a mino modification of the NG dimensionless enegy consevation equation 6) in which only the tansvese component of the DKE is multiplied by γ. This equation is analogous to the known GR equation fo planetay motion, as the geodesic of the Schwazschild metic. Ou model also eveals the souce of the pecession of the planetay obit. As mentioned above, in NG, the adial and the tansvese peiods ae identical, esulting in a non-pecessing obit. In SR, both the adial and tansvese components of the velocity ae alteed, esulting in unequal peiods with elatively small diffeence between them and hence a small pecession. In ou model, only the adial component of the velocity is influenced, while the tansvese angula) component is not. This, in tun, accentuates the diffeence between these peiods, esulting in the obseved pecession, as follows. To define the pecise value of the pecession, we ewite equation 4) as J c s u ) = J c s u 3 J c s u +u+e. Multiplying this equation by µ, with µ defined by 9) we obtain u ) = u 3 u +µu+µe. 5) This equation is identical to 0) in NG, except that it has a vey small since u ) additional tem u 3 on the ight-hand side. This esult is the same esult as that of GR. We seek a solution of this equation in the fom genealizing ), uϕ) = µ+εcosαϕ)) 6) fo some function αϕ). As befoe, two oots of the cubic on the ight-hand side of 5), ae the u values u p and u a of the the peihelion and aphelion, espectively. Moeove, since the coefficients of this cubic ae constant fo a given obit, these values will not change fom one evolution to the next. We denote the thid oot of this cubic by u e. Thus, equation 5) can be factoized as u ) = u u p )u u a )u e u). 7) Fom equation 6), u ) = α ) µ ε sin αφ), u p = µ+µε and u a = µ µε. Moeove, since the sum of the oots of this cubic is, u e = u p +u a ) = µ. 8) Substituting these into 7), yields afte simplification α = dα dϕ = 3µ µεcosαϕ))/. This allows us to obtain the dependence of ϕ on α as ϕα) = ϕ 0 + α 0 3µ µεcos α) / d α. 9) As mentioned above, fo Mecuy, µ is vey small. Expanding the integand into a powe seies in µ, we obtain ϕα) = ϕ 0 +α+ 3 µα+ µ εsinα+. The pola angles of the peihelion coespond to α = πn, n = 0,,,. Thus, the pecession of the peihelion, up to the fist ode in µ, is given by ϕπ) ϕ0) π 3πµ. 30) Fom this equation, the pecession of the peihelion of Mecuy is 5.0866 0 7 adians pe evolution, which is exactly the cuently obseved one. Finally, we show how we can ecove the Schwazschild metic fom equation 3). Since the influenced diection in ou model is the adial diection, ou metic will diffe fom the Minkowski metic only in the time and adial components. Thus, ou metic in spheical coodinates is ds = f)dx 0 g)d) dθ) +cos θdϕ) ). 3) p-4

Pedicting Mecuy s Pecession The obit is the geodesic of this metic optimizing the Lagangian Lx 0,,ϕ,θ,ẋ 0,ṙ, ϕ, θ), L = f)ẋ 0 ) g)ṙ) θ) +cos θ ϕ) ) whee the dot denotes diffeentiation with espect to s. Assuming θ = π/, the Eule-Lagange equations yield f)ẋ 0 = a fo some constant a and ϕ = J/c. The nom of the 4-velocity on the obit is f)ẋ 0 ) g)ṙ) ϕ) = a f) g)ṙ) J c. Compaing this with 3) we obtain f) = s / and g) = s/, implying that the metic is the Schwazschild metic. Discussion. Ou model fo planetay motion pesented hee is a special case of a moe geneal Newtonian elativistic dynamics, which will be pesented in [9]. In this dynamics, the enegy consevation equation fo motion unde a consevative foce in a spacetime influenced by the potential enegy Ux) vanishing at infinity) is m u)ẋ +uẋ n) ) +Ux) = E, 3) [3] Ø. Gøn and H. Sigbjøn,Einstein s Geneal Theoy of Relativity: With Moden Applications in Cosmology Spinge) 007. [4] K. Bown, Reflections on Relativity lulu.com) 05. [5] C. Sigismondi, Poceedings of the Tenth Macel Gossmann Meeting on Geneal Relativity Wold Sci.) 005. [6] Y. Fiedman, Homogeneous Balls and thei Physical Applications Bikhäuse, Boston) 004. [7] Y. Fiedman and Yu. Gofman, Phys. Sc., 8 00) 05004. [8] Y. Fiedman, Poceedings of IARD 0, Jou. of Phys. Confeence Seies, 437 03) 007. [9] Y. Fiedman, Relativistic Newtonian Dynamics unde a Consevative Foce, 06, in pepaation. whee u = U mc, n = U U and dot denotes diffeentiation with espect to time. In the case of a adial potential U), n is in the adial diection and the consevation of angula momentum defines explicitly the velocity decomposition into its adial and tansvese components. Equation 3) yields ) d = J dt u) U m + E m, 33) which educes to 4) fo planetay motion. Ou model consides only the influence of the gavitational potential and ignoes the influence of the kinetic enegy. A complete model should conside both these effects. Nevetheless, ou model coectly pedicts Mecuy s anomalous pecession eithe because the influence of the kinetic enegy is below the expeimental accuacy o, its effect is cancelled out by the inteaction with the neighbouing planets. We wish to acknowledge Pof. Lawence Paul Howitz and M. Zvi Weinbege fo thei constuctive comments. Refeences [] W. Rindle, Relativity, Special, Geneal and Cosmological Oxfod) 00. [] S. Kopeikin, M. Efoimsky and G. Kaplan,Relativistic Celestial Mechanics of the Sola System Wiley-VCH, Belin) 0. p-5