The perihelion of Mercury advance and the light bending calculated in (enhanced) Newton’s theory

We show that results of a simple dynamical gedanken experiment interpreted according to standard Newton’s gravitational theory, may reveal that three-dimensional space is curved. The experiment may be used to reconstruct the curved geometry of space, i.e. its non-Euclidean metric 3gik\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^3g_{ik}$$\end{document}. The perihelion of Mercury advance and the light bending calculated from the Poisson equation 3gik∇i∇kΦ=-4πGρ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^3g^{ik} \nabla _i \nabla _k \varPhi = -4\pi G \rho $$\end{document} and the equation of motion Fi=mai\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$F^i = ma^i$$\end{document} in the curved geometry 3gik\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^3g_{ik}$$\end{document} have the correct (observed) values. Independently, we also show that Newtonian gravity theory may be enhanced to incorporate the curvature of three dimensional space by adding an extra equation which links the Ricci scalar 3R\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^3R$$\end{document} with the density of matter ρ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\rho $$\end{document}. Like in Einstein’s general relativity, matter is the source of curvature. In the spherically symmetric (vacuum) case, the metric of space 3gik\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^3g_{ik}$$\end{document} that follows from this extra equation agrees, to the expected accuracy, with the metric measured by the Newtonian gedanken experiment mentioned above.

in Einstein's general relativity, matter is the source of curvature. In the spherically symmetric (vacuum) case, the metric of space 3 g ik that follows from this extra equation agrees, to the expected accuracy, with the metric measured by the Newtonian gedanken experiment mentioned above.

Introduction
Newton's theory of gravity was formulated in a flat, Euclidean 3-D space but its basic laws, i.e. the Poisson equation and the equation of motion, make perfect sense in a 3-D space with an arbitrary geometry 3 g ik . Indeed, the curvature of space is (potentially) present in Newton's theory. It is easy to argue that the "centrifugal" acceleration of a particle moving with velocity V on a circular orbit equals a C = V 2 /R, and the "gravitational" acceleration in the gravitational field of a spherically symmetric body with the mass M equals a G = −GM/r 2 , with R being the curvature radius of the circle, andr being its circumferential radius. In flat, i.e. Euclidean, 3-D space these two radii are equal, R =r , but in a space with a non-zero gaussian curvature G, they are different, R =r . Therefore, by measuring centrifugal and gravitational accelerations one may independently measure R andr , and thus experimentally find whether the space is flat (Euclidean) or it has a non-zero Gaussian curvature G = 0. Based on that, Abramowicz has recently suggested in [3] that a Newtonian physicist could experimentally determine the metric 3 g ik of the real physical 3-D space and calculate, according to (1) and (2), the perihelion of Mercury advance and the light bending effects. In this paper we follow this suggestion and calculate both effects within Newton's theory. Surprisingly, the values of the perihelion advance and the light bending agree (to the expected order of M/r ) with predictions of Einstein's theory. Here M is the "geometrical" mass of the spherical gravitating body expressed in the convenient "geometrical" units G = 1 = c. It is connected to the mass M expressed in the standard units by M = GM/c 2 and has the dimension of length. Another point discussed in this paper is based on the following two remarks: (i) Obviously, Newton's gravity theory is a limit of Einstein's general theory of relativity. Should the limit necessarily correspond to G = 0? Perhaps not, because Newtonian physicists could discover within Newton's theory that G = 0. (ii) They could also discover that the curvature of space depends on the distance from the gravity center. This would suggest to them, again within the framework of Newton's theory, that gravity and curvature are not independent, but instead they are somehow linked. Here we suggest that it is possible to establish the link within an "enhanced" version of Newton's theory, by adding to its standard version defined by (1) and (2) an extra equation, Fig. 1 For a circle placed in a curved space (here on a curved 2-D surface), its geodesic radius r * , circumferential radiusr , and curvature radius R are all different, r * =r = R where 3 R is the Ricci scalar corresponding to 3 g ik , ρ is the density of matter, and k is a constant. Equations (1), (2) and (3) define our enhanced version of Newtonian gravitational theory. In the special case of a spherically symmetric, vacuum (ρ = 0) space, they uniquely lead to the 3-D metric of the form, where r 0 is a constant. A choice r 0 = 4 M leads to correct values for both the perihelion advance and the light bending effects. 1 Several authors have discussed the idea of Newton's gravity and dynamics in a curved 3-D space, for example most recently Naresh Dadhich in "Einstein is Newton with space curved" [1] or previously, in the context of the optical geometry, [2].

The three radii of a circle
Consider a two dimensional curved, axisymmetric surface. In general, it has the metric, Consider a family of concentric circles r * = const in it. One of them is shown in Fig. 1.
Let τ i =r −1 δ i φ be a unit vector tangent to the circle. From the Frenet formula, one deduces that the curvature radius R may be defined by, Two useful formulae for the curvature of the circle, K = 1/R, and for the Gaussian curvature G of the surface with the metric (5) read, Formula (10) follows from (9). For derivation of (11) see e.g. [4], Section 3.4.

Equations of motion
Let us consider a curve in space given by a parametric equation, where x i are coordinates in space, and s is the length along the curve. If a body moves along this curve, its velocity equals, Here V = ds/dt is the speed of the body and τ i = dx i /ds is a unit vector tangent to the curve (12), i.e. the direction of motion. The acceleration may be calculated as follow, Assuming circular motion with constant velocity, V = const, and applying (9) to calculate the term in brackets,we arrive at which is the well known formula for the centripetal acceleration. Therefore the radial component of the centrifugal acceleration is Consider now circular motion around a spherically symmetric center of gravity. The Newtonian equation of motion, F i = ma i , takes the form, where F i = −m∇ i Φ is the gravitational force, and Φ is the gravitational potential. Three quantities characterize motion on a particular circular orbit: the angular velocity Ω, the angular speed V , and the specific angular momentum L. They are related by, Multiplying left hand side of the equation of motion (17) by λ i and using (19), we transform (17) into a form which will be convenient later, In this expression, λ i is a unit, outside pointing, vector. Here "outside" has the absolute meaning-outside the center, in the direction towards infinity. We will calculate the left-hand side of this equation in the next Section.

Newton's gravity and Kepler's law
In an empty space, the gravitational potential Φ obeys the Laplace equation, Let us integrate (21) over the volume V that is contained between two spheres, concentric with the gravity center, with sphere S 1 being inside sphere S 2 . We transform the volume integral into a surface integral, using the Gauss theorem The oriented surface elements on the spherical surfaces S 1 and S 2 may be written, respectively, as therefore, This means that the value of the integral is the same, say S 0 , for all spheres around the gravity center. In addition, because of the spherical symmetry of the potential, the quantity (∇ i Φ)λ i is constant over the sphere of integration. Thus, From the above expression it follows that the radial component of the gravitational acceleration equals, Combining (25) with (20), we may finally write, This is the Kepler Third Law. Using natural units for radius and frequency, we may write the formula for the Keplerian angular velocity as,

A Newtonian experiment
Imagine a Newtonian physicist who measures the gravitational a G and the centrifugal a C accelerations for a circular motion (with the orbital velocity V ) around the Sun. In the c.g.s. (or any other unit system he may use), a G and a C are just numbers. They are independent on the measuring setup and on the theory used. As we have already explained, in Newton's theory these measurements are equivalent to measurements of the circumferential and curvature radii of the circular orbits, because in Newton's theory it follows from Eqs. (16) and (26) that, In the flat Euclidean space it must ber = R, and therefore the sufficient experimental Newtonian condition for the 3-D space to be curved is, The Newtonian physicist may form a table of pairs [r , R], which follow from his measurements of [a G , a C ] at many different orbits. Thus, he may experimentally establish the function, The geometrical information contained in this function is sufficient to reconstruct the metric of any spherically symmetric 3-D space (spherical symmetry with the Sun at the center would be a natural "working hypothesis" adopted by the Newtonian physicist). The Newtonian physicist may achieve this reconstruction without referring to anything out of the framework of Newton's physics. Obviously, however, he must allow for a non-zero curvature of the 3-D space. The Solar System dynamics is given by the Einstein geodesic motion in the Schwarzschild metric; this metric is the "physical reality" here. Thus, when we calculate results of the gedanken experiment considered now, we must use the Schwarzschild metric. We will now calculate the results of the measurements of the circumferential and curvature orbital radii, expressing them by the Schwarzschild standard radial coordinate r , To avoid a possible confusion here, we stress that Eq. (33) should be understood as a parametric representation of the experimentally established, invariant, function (32). This function does not refer to r explicitly. No explicit knowledge of r is needed for any argument of the Newtonian physicist -he is only usingr and R, which he directly measures. Similarly, the final formulae for the perihelion shift and light bending depend neither on r , nor on any non-Newtonian concepts (spacetime, timelike Killing vectors, etc.) that we have used to calculate the results of the Newtonian gedanken experiment.
In the Schwarzschild metric, the acceleration of a particle (a "planet") moving with the orbital velocity v along a circular orbit equals, Here V = v/(1 − v 2 ) 1/2 , and the scalars Ψ and R are expressed in terms of the time-symmetry Killing vector η i , and the axial-symmetry Killing vector ξ i , In Schwarzschild coordinates this is, at the "equatorial plane" θ = π/2, Equation (34) allows to predict the acceleration measurement results. The acceleration measured with V = 0, interpreted by the Newtonian physicist as gravitational acceleration a G , corresponds to the first term on the right hand side of (34). Velocity dependent part of the acceleration, on the other hand, corresponding to the second term on the right hand side of (34), would be interpreted as the centrifugal acceleration, i.e., By comparing (30) and (37), one concludes that, The above formula allows one to write the metric of the 2-D space geometry of the equatorial plane θ = π/2, measured in this Newtonian experiment, in the form The second line of this Eq. (41), which in an obvious tautology for any axially symmetric 2-D metric, shows nevertheless that the metric "reconstructed" from the Newtonian experiment described in this Section depends only onr and R(r ), i.e. quantities that are measured by the Newtonial physicist. It needs no reference to r . However, one may use r as a convenient parameter, as in the first line of this Eq. (40).

Epicyclic oscillations, the perihelion advance
Suppose that we slightly perturb a test-body on a circular orbit. This means that its angular momentum will not correspond to the Keplerian one, L 2 , given by (27), but will be slightly different L 2 + δL 2 . There will be also a small radial motion with velocity(δr * ) and acceleration(δr * ). From (20) it follows that Keeping the first order term in Eq. (42), and using we arrive at the simple harmonic oscillator equation, where ω is the radial epicyclic frequency, Using Eqs. (27) and (28), we may write the expression for the epicyclic frequency in the form, or comparing this with (29), In a flat space, r * =r = R, and therefore ω = Ω, which implies that the slightly non-circular orbit is a closed curve, indeed an ellipse. In a curved space with G = 0, one has r * =r = R, and consequently ω = Ω. The slightly non-circular orbit would not be a closed curve. It could be represented by a precessing ellipse, with two consecutive perihelia shifted by where T = 2π/Ω is the orbital period. Inserting (38) into the Newtonian perihelion advance formula (48) one gets, Expanding this to the desired accuracy O 2 (M/r ), and using Eq. (38), one finally gets the same value for the perihelion advance as calculated in Einstein's theory,

Light bending
Knowing the space geometry, given by (40), we may calculate the effect of light bending assuming that light travels along geodesic lines in space. In Newton's theory this assumption is equivalent to the Fermat principle, i.e. that light travels (with a constant speed) between two points A, B in space, minimizing the time travel T AB . The equation of motion for the φ coordinate is, in these circumstances, from which we find dφ ds to be equal to .
The integration constant can be evaluated at the perihelion location r = R 0 (i.e. where dφ/ds = 0), yielding Using Eqs. (40) and (53) we find also After dividing (53) by (54) and substituting x = R 0 /r the dφ/dr equation can be integrated from R 0 to ∞ (or x from 0 to 1), which will give us the half of π + δ. Let us also define μ = M/R 0 , then This integration can be expanded in a Taylor series for μ: As the first component on the right hand side is equal to π/2, we conclude that Similar calculations in the Schwarzschild spacetime geometry give the same result 2 Once again, the prediction of the Newtonian theory in the non-flat space is found to be consistent with observations (and with Einstein's general relativity).
In the real curved space of the Solar System (described by the Schwarzschild metric) a Newtonian physicist may discover the curvature by the two specific measurements that we have described, and then correctly estimate the value of the perihelion of Mercury advance and the light bending. Of course, he may perform different measurements, for example he may directly measure the circumferential and curvature radii. If he would use light signals to measure the distance (radar defined distances leading to the optical geometry), then the result will be the same as from measuring these radii "dynamically" (see e.g. Abramowicz [5]). However, if he would measure these radii using the "rigid rods", then he would find a different result. To get the right answers for the light bending and orbital precession, he would need to make the right choice of the measure to use. He may eventually find that the optical geometry agrees with dynamics, but the geometry based on rigid rods does not.

Enhanced Newtonian gravitational theory
Jürgen Ehlers pointed out in 1961 that in Einstein's theory the curvature of the restspace of irrotational matter is determined by its distribution and relative motion (see 1221 in his article [6]). The equations governing such 3-space curvature for arbitrary irrotational flows are given in [7]; see their Eq. (54). Consequently it makes sense to consider gravitational dynamics in the context of 3-D curved Riemannian spaces. As Newtonian theory is an approximation to General Relativity Theory, it is therefore interesting to see what happens in the case of Newtonian theory in a curved 3-D background space 3 .
In the case of isometric flows, θ = σ ab = 0 and there is a potential such thaṫ u a = U ,a where the gravitational potential U relates the Killing vector ξ to the unit 4velocity u a by ξ a = e U u a (see 1234 in Ehlers [6]). Then the relevant equation becomes where∇ a is the 3-D covariant derivative, ρ is the energy density of matter, and we have assumed anisotropic stress is zero (π ab = 0) and a vanishing cosmological constant.
Here h ab = g ab − u a u b is the metric of the three-spaces orthogonal to u a . This case will include static spherically symmetric spacetimes. Taking the trace of this equation gives (see equation (55) in [7]) where the potential terms have gone because of the relation between the 3-D and 4-D covariant derivatives. Together with the Poisson equation and equation of motion it defines the Enhanced Newtonian Gravitational Theory, For spherically symmetric spaces, the most general metric has the form, and in the vacuum case, ρ = 0 = 3 R one has Here the prime denotes a derivative with respect to r . Equation (66) has a unique solution, with r 0 being an integration constant. Its value cannot be determined by Eq. (61), but instead must be chosen by correspondence with experiment 4 . Using the same procedure as in Sects. 6 and 7, one proves that the choice r 0 = 4M gives the correct values for the perihelion advance and light bending (with accuracy O(r 0 /r )).

The two metrics
We have shown that "experimentally" established and the "theoretically" postulated Newtonian metrics of the curved 3-D space corresponding to a spherically symmetric body are, respectively, We have also shown that any spherically symmetric metric that obeys 3 R = 0 must be isometric with (69). The Ricci scalar for the "experimental" metric may be calculated to be On the other hand a metric, has the Ricci tensor, Thus, the experimental metric (68) and the theoretical metric (69) describe, with accuracy O 2 (M/r ), the same geometry of space.

Conclusions
We demonstrated that a Newtonian physicist may experimentally determine the geometry of the 3-D space 3 g E ik by measuring gravitational and centrifugal accelerations. He may then predict by calculations the perihelion advance and the light bending as effects of the curvature of space. The predicted values agree with the ones measured. We also demonstrated that one may extend Newton's theory of gravitation by adding an equation that links Ricci curvature of space with the density of matter. We calculated the resulting theoretical metric of space 3 g T ik assuming spherical symmetry. In this metric, the values of perihelion advance and light bending also agree with those observed. The two metrics represent the same geometry, 3 g E ik = 3 g T ik with accuracy O 2 (M/r ).
Abramowicz [3] has shown that for spaces with constant Gaussian curvature Newton's theory predicts no perihelion advance. We speculate that this is why Gauss (and other XIX century mathematicians) who might have calculated Newtonian orbits in curved spaces, would have missed the effect of perihelion advance. Most probably, they would calculate orbits in spaces with a constant Gaussian curvature first. Gauss almost certainly made this calculation. He was a master in calculating orbits. He made himself famous at the age of 23 by calculating the orbit of Ceres, discovered in 1801 by Piazzi. He seriously considered the possibility that our space is curved. He even attempted to determine the curvature of space by measuring angles in a big triangle (69, 84, 106 km) made by the summits of Brocken, Hoher Hagen and Großer Inselsberg. Gauss was not quick in publishing his results concerning curved spaces. It is known that he discovered most of Bolyai's results, but never published them. Gauss died in February 1885, four years before Le Verrier discovered the effect of the perihelion of Mercury advance.