# Detection of the relativistic corrections to the gravitational potential using a Sagnac interferometer.

General Relativity predicts the existence of relativistic corrections to the static Newtonian potential which can be calculated and verified experimentally. The idea leading to quantum corrections at large distances is that of the interactions of massless particles which only involve their coupling energies at low energies. In this short paper we attempt to propose the Sagnac intrerferometric technique as a way of detecting the relativistic correction suggested for the Newtonian potential, and thus obtaining an estimate for phase difference using a satellite orbiting at an altitude of 250 km above the surface of the Earth.1 Introduction

The potential acting between to masses M and m that separated from their centers by a distance r is:

V(r) = -GMm/r, (1)

where s the Newton's constant of gravitation. This potential is of course only approximately valid [1]. For large masses and or large velocities the theory of General Relativity predicts that there exist relativistic corrections which can be calculated and also verified experimentally [2]. In the microscopic distance domain, we could expect that quantum mechanics, would predict a modification in the gravitational potential in the same way that the radiative corrections of quantum electrodynamics leads to a similar modification of the Coulombic interaction [3].

Even though the theory of General Relativity constitutes a very well defined classical theory, it is still not possible to combine it with quantum mechanics in order to create a satisfied theory of quantum gravity. One of the basic obstacles that prevent this from happening is that General Relativity does not actually fit the present paradigm for a fundamental theory that of a renormalizable quantum field theory. Gravitational fields can be successfully quantized on smooth-enough space-times [4], but the form of gravitational interactions is such that they induce unwanted divergences which can not be absorbed by the renormalization of the parameters of the minimal General Relativity [5]. Somebody can introduce new coupling constants and absorb the divergences then, one is unfortunately led to an infinite number of free parameters. In spite the difficulty above quantum gravity calculations can predict long distance quantum corrections.

The main idea leading to quantum corrections at large distances is due to the interactions of massless particles which only involve their coupling energies at low energies, something that it is known from the theory of General Relativity, even though at short distances the theory of quantum gravity differs resulting to finite correction of the order, O (G[??]/[c.sup.3][r.sup.e]). The existence of a universal long distance quantum correction to the Newtonian potential should be relevant for a wide class of gravity theories. It is well known that the ultraviolet behaviour of Einstein's pure gravity can be improved, if higher derivative contributions to the action are added, which in four dimensions take the form:

[alpha][R.sup.[kappa][lambda]][R.sub.[kappa][lambda]] + [beta][R.sup.2] (2)

where [alpha] and [beta] are dimensionless coupling constants. What makes the difference is that the resulting classical and quantum corrections to gravity are expected to significantly alter the gravitational potential at short distances comparable to that of Planck length [l.sub.p] = [square root of G[??]/[c.sup.3]] = [10.sup.-35] m, but it should not really affect its behaviour at long distances. At long distances is the structure of the Einstein-Hilbert action that actually determines that. At this point we should mentioned that some of the calculation to the corrections of the Newtonian gravitational potential result in the absence of a cosmological constant [LAMBDA] which usually complicates the perturbative treatment to a significant degree due to the need to expand about a non-flat background.

In one loop amplitude computation one needs to calculate all first order corrections in G, which will include both the relativistic O([G.sup.2][m.sup.2]/[c.sup.2]) and the quantum mechanical O ([G[??]/[c.sup.3]) corrections to the classical Newtonian potential [6].

2 The corrections to the potential

Our goal is not to present the details of the one loop treatment that leads to the corrections of the Newtonian gravitational potential but rather sate the result and then use it in our calculations. Valid in order of [G.sup.2] we have that the corrected potential now becomes [6]:

V(r) = - GMm/r [1 - G(M + m)/2[c.sup.2]r - 122G[??]/15[pi][c.sup.3][r.sup.2]. (3)

Observing (3) we see that in the correction of the static Newtonian potential two different length scales are involved. First, the Planck length [l.sub.p] = [square root of G[??]/[c.sup.3] = [10.sup.-35] m and second the Schwarzschild radii of the heavy sources [r.sub.sch] = 2G[M.sub.n]/[c.sup.2]. Furthermore there are two independent dimensionless parameters which appear in the correction term, and involve the ratio of these two scales wit respect to the distance r. Presumably for meaningful results the two length scales are much smaller than r.

3 Perturbations due to oblateness [J.sub.2]

Because the Earth's gravitational potential is not that of a perfect spherical body, we can approximate its potential as a spherical harmonic expansion of the following form:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (4)

where:

r = geocentric distance,

[phi] = geocentric latitude.

[R.sub.e] = means equatorial radius of the Earth,

[P.sub.n] = Legendre polynomial of degree n and order zero,

[J.sub.n] = [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] jonal harmonics of order zero, that depend on the latitude [phi] only,

and the first term GMm/r now describes the potential of a homogeneous sphere and thus refers to Keplerian motion, the remaining part represents the Earth's oblateness via the zonal harmonic coefficients and [7]

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (5)

Therefore equation (4) can be further written:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (6a)

Since [J.sub.2] is 400 larger that any other [J.sub.n] coefficients, we can disregard them and write the following expression for the Earth's potential function including only the relativistic correction and omitting the quantum corrections as being very small we have:

V (r, [phi]) = - G[M.sub.e]m/r + G[M.sub.e]m[R.sup.2.sub.e][J.sub.2]/[r.sup.3] (3/2 [sin.sup.2][phi] - 1/2) +

+ [G.sup.2][M.sub.e]m ([M.sub.e] + m)/[c.sup.2 - [r.sup.2] (7)

Since we propose a satellite in orbit that carries the Sagnac instrument it will of a help to express equation (7) for the potential in terms of the orbital elements. We know that sin [phi] = sin i sin(f + w) where i is the inclination of the orbit, f is the true anomaly and w is the argument of the perigee. Ignoring long and short periodic terms (those containing w and f) we write (7) in terms of the inclination as follows:

V (r, [phi]) = - G[M.sub.e]m/r + 3G[M.sub.e]m[R.sup.2.sub.e][J.sub.2]/2[r.sup.3] ([sin.sup.2]i/2 - 1/3) +

+ [G.sup.2][M.sub.e]m ([M.sub.e] + m)/[c.sup.2 - [r.sup.2] (8)

therefore the corresponding total acceleration that a mass m at r > [R.sub.e] has becomes:

[g.sub.tot] = -1/m [partial derivative]/[partial derivative]r [-G[M.sub.e]m/r + 3G[M.sub.e]m[R.sup.2.sub.e] [J.sub.2]/2[r.sup.3] x

x ([sin.sup.2]i/2 - 1/3) + [G.sup.2][M.sub.e]m(M + m)/[r.sup.2][c.sup.2]] (9)

so that:

[g.sub.tot] = -G[M.sub.e]/[r.sup.2] + 9G[M.sub.e][R.sup.2.sub.e][J.sub.2]/2[r.sup.4] ([sin.sup.2]i/2 - 1/3) +

+ [G.sup.2][M.sub.e] ([M.sub.e] + m)/[c.sup.2][r.sup.3]. (10)

4 Basic Sagnac interferometric theory

The Sagnac interferometer is based on the Sagnac effect, reported by G. Sagnac in 1913 [8]. Two beams are sent in opposite directions around the interferometer until they meet

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (15)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (16)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (17)

again to create a phase pattern. By rotating the interferometer in the direction of either the clockwise (CW) or counter-clockwise (CCW) beam, a phase difference results between the two beams that its given by:

[DELTA][[PHI].sub.rs] = 8[[pi].sup.2][R.sup.2.sub.sag]N[OMEGA][upsilon]/([c.sup.2] - [a.sup.2][[OMEGA].sup.2]), (11)

where [OMEGA] is the angular velocity of the interferometer, [R.sub.sag] is the radius of the interferometer, N is the number of turns of fiber around the radius and [upsilon] is the frequency of light in the fiber.

Let us now assume that the Sagnac interferometer and its light laser beams are in the region of space around the Earth where the gravitational potential is given by equation (3) and let us further assume that the quantum correction to the potential is really negligible. If the Sagnac light loop area has a unit vector that is perpendicular to the acceleration of gravity vector, then the motion of the interferometer will exhibit a red-shift that will be given by:

[f.sub.rs] = f/1 - [DELTA]V/[c.sup.2] = f/1 - gcorz/[c.sup.2], (12)

where [DELTA]V is the difference in the potential between to different points [P.sub.1] and [P.sub.2], and [g.sub.cor] is the corrected or total acceleration of gravity and z is the difference in vertical distance between the two beams as the interferometer coil rotates. This distance z that the laser beams see is given by:

z = [R.sub.sag] {1 - cos [2[pi][OMEGA][R.sub.sag] / c(1 + [R.sub.sag][OMEGA]/c)]}. (13)

This Sagnac effect can also be amplified by an interferometer that is in orbit, where the orbital velocity of the interferometer with respect to the Earth's surface produces an increased phase shift. Both terms involved in the acceleration of gravity in the first one:

[DELTA][[PHI].sub.rs] = 8[[phi].sup.2][R.sup.2.sub.sag]N[OMEGA][upsilon]([R.sub.sag][OMEGA] + [v.sub.orb]/[[pi].sup.2])/([c.sup.2] - [R.sup.2.sub.sag] [[OMEGA].sup.2]) [1 - [g.sub.tot] z / [c.sup.2]] (14)

using (14) and taking into account that M [much greater than] m we further obtain (15), where M is the source of the gravitational field = the mass of the Earth in our case [M.sub.e], and R is the radius of the massive body = [R.sub.e], and r = [R.sub.e] + [z.sub.orb] it's orbital height plus Earth radius for an Earth-based satellite.

This Sagnac effect can also be amplified by an interferometer that is in orbit, where the orbital velocity of the interferometer with respect to the Earth's surface produces an increased phase shift. Both terms involved in the acceleration of gravity in the first one:

5 Sagnac in circular orbit of known inclination

Let now a Sagnac interferometer be aboard a satellite in a circular polar orbit of inclination i = 90 degrees. If the inclination is 90 degrees the term [sin.sup.2] i/2 - 1/3 = 1/6 and the orbital velocity at some height z above the surface of the Earth is [v.sub.orb] = [square root of G[M.sub.e]/([R.sub.e] + [Z.sub.orb]) and (6) takes the form (16) can be finally written as (17).

6 Sagnac in elliptical orbit of known inclination

If now a satellite is carrying a Sagnac device is in an elliptical orbit of eccentricity e and semi-major axis a we have that the radial orbital vector and the orbital velocity are given by:

r(f)= a(1 - [e.sup.2])/1 + e cos f, (18)

[v.sup.2] = G[M.sub.e](2/r - 1/a) = G[M.sub.e]/a [2(1 + e cos f)/(1 - [e.sup.2]) - 1], (19)

where f is the true anomaly of the orbit. Substituting now in (8) we obtain (20).

If we use the fact that G[M.sub.e] = [n.sup.2][a.sup.3] where n is the mean motion of the satellite, equation (20) can be further written as (21).

When the satellite approaches perigee its orbital velocity will increase, so we will expect to see a higher phase difference than any other point of the orbit, and similarly the effect

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (20)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (21)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (24)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (25)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (26)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (27)

will be minimum at the point of apogee because the satellite's velocity is minimal. The distance at perigee and apogee are given by the equations below:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (22)

also the corresponding velocities are:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (23)

therefore the phase difference detected by the Sagnac due to the contribution of the Earth's oblateness plus relativistic correction to the potential at perigee and apogee can be written as (24) or again (25).

Similarly the phase difference at apogee can be written as (26) or again (27).

For this last case of the elliptical orbit in (25) and (26) where the Sagnac interferometer is on the satellite and we assume [R.sub.S] = 1 m, [upsilon] = 2 x [10.sup.14] Hz, N = [10.sup.6], [[OMEGA].sub.s] = 400 rad/sec, a = 8 x [10.sup.6] m, e = 0.2, [R.sub.e] = 6.378 x [10.sup.6] meters we arrive at the following values for [DELTA][phi]:

[DELTA][phi] (perigee) = 3.57 x [10.sup.-16] radians, [DELTA][phi] (apogee) = 2.44 x [10.sup.-16] radians.

These values are based on the dominant potential correction in (11) of section 3 which is the first term in (11) or the Newtonian correction:

Newtonian correction = 2.17 x [10.sup.-16] radians.

In comparison, the second and third terms in (11) are the oblateness and relativistic corrections respectively and they produce the following values based on the given parameters:

Oblateness correction = 8.52 x [10.sup.-20], Relativistic correction = 7.91 x [10.sup.-26].

So by comparison of the values above, the Newtonian correction is much easier to measure.

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (30)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (31)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (32)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (33)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (34)

The [DELTA][phi] values given above may be more easily measured using a QPSK-modulator inserted in the CCW or CW beam path to improve phase resolution. Also, the use of higher wavelengths (factor of 10 higher in frequency) will increase resolution.

7 We suggest a Sagnac with an elliptic fiber loop

To attempt increasing the resolution of the phase difference of the Sagnac interferometer let us now propose a Sagnac loop, that has the shape of an ellipse that rotates with an angular velocity [OMEGA]. In this case it can be shown that the height difference between two points on the ellipse can be given by:

z = a [1 + ([e.sup.2] + e - 1) cos [theta]/1 + e cos [theta]]. (28)

To check the validity of the formula we derived we can set e=0 which is the case of a circular Sagnac fiber optical path we can see that the (13) in now retrieved since

[R.sub.sag] = [a.sub.loop(sag)] = [a.sub.s] is the semi major axis of the elliptical fiber loop. When the ellipse spins with angular velocity [OMEGA] that would force it to trace out a circle whose radius r, will be that of the semi-major axis a of the ellipse, and therefore we can finally write for (13):

z = [a.sub.s] [1 + ([e.sup.2] + e - 1) cos{2[pi][a.sub.s][[OMEGA].sub.s]]/c (1 + [a.sub.s][[OMEGA].sub.s]/c).sup.-1]}]/ (1 + e cos {2[pi][a.sub.s][[OMEGA].sub.s]/c (1 + [a.sub.s] [[OMEGA].sub.s]/c).sup.-1}) (29)

8 Circular orbit formula for the phase difference of the Sagnac

Let now as before have a Sagnac interferometer be aboard a satellite in a circular polar orbit of inclination i = 90 degrees. If the inclination is 90 degrees the term [sin.sup.2] i/2 - 1/3 = 1/6 and the orbital velocity at some height z above the surface of the Earth is [v.sub.orb(circ)] = [square root of G[M.sub.e]/([R.sub.e] + [z.sub.orb]) and (6) takes the form (30) that can be finally written as (31).

9 Sagnac in elliptical orbit of known inclination

If now a satellite is carrying a Sagnac device is in an elliptical orbit of eccentricity e and semi-major axis a we have that the radial orbital vector and the orbital velocity are given by (32).

At perigee the equation (32) becomes (33) and also (34).

For (33) and (34) above the following values are computed assuming e = 0.2, [upsilon] = 2 x [10.sup.14] Hz, a = 8 x [10.sup.6] meters, N = 1 (because the orbit is the Sagnac loop), [R.sub.sag] = [R.sub.perigee] or [R.sub.apogee] as determined by (22), [[OMEGA].sub.perigee] = 0.001 rad/sec, and [[OMEGA].sub.apogee] = 6 x [10.sup.-4] rad/sec we find,

[DELTA][phi] (perigee) = 6.05 x [10.sup.10] radians, [DELTA][phi] (apogee) = 2.36 x [10.sub.10] radians.

These values are for measuring the dominant Newtonian contribution as described in Section 6. To detect relativistic contribution which is 3.64 x [10.sup.-10] smaller than the Newtonian contribution the corresponding phase-shifts from (33) and (34) are:

[DELTA][phi] (perigee) = 22 radians, [DELTA][phi] (apogee) = 8.59 radians.

Thus, the relativistic contribution in (11) of Section 3 is easily measurable using a Sagnac interferometer where the satellites in orbit are the Sagnac loop. In this scenario, the light path can be implemented by transmitting laser beams from one satellite to the next satellite in orbit ahead of it. Also, by using the maximum spacing possible between satellites in orbit this will allow line of site transmission while reducing the number of satellites required for the Sagnac loop. With the potential to measure such small relativistic corrections, the merit of using satellites to implement a large Sagnac loop of radius [R.sub.s] = [R.sub.ap] or [R.sub.per] is well worth considering.

Submitted on March 27, 2008 Accepted on April 04, 2008

References

[1.] Donoghue J. F. arXiv: gr-gc/9310024.

[2.] Bjorken J. D. and Drell S. Relativistic quantum mechanics and relativistic quantum field theory. Mc Graw-Hill, New York 1964.

[3.] t' Hooft G. and Veltman M. Ann. Inst. H. Poincare, 1974, v. 20, 69.

[4.] Capper D. M., Leibrand G. and Medrano R. M. Phys. Rev. D, 1973, v. 8, 4320.

[5.] Goroff M. and Sagnotti A. Nucl. Phys. B, 1986, v. 22, 709.

[6.] Hamber H. W. and Liu S. arXiv: hep-th/9505128.

[7.] Blanco V M. and Mc Cuskey S. W. Basic physics of the Solar system. Addison-Wesley Publishing Inc., 1961.

[8.] Sagnac G. R. Acad. Sci., 1913, v.157, 708.

[9.] Camacho A. General Relativity and Gravitation, 2004, v.36, no. 5, 1207-1211.

[10.] Roy A. E. Orbital motion. Adam Hilger, Third Edition, 1988.

Ioannis Iraklis Haranas * and Michael Harney ([dagger])

* Department of Physics and Astronomy, York University, 314A Petrie Building, North York, Ontario, M3J-1P3, Canada E-mail: ioannis@yorku.ca

([dagger]) 841 North 700 West, Pleasant Grove, Utah, 84062, USA E-mail: michael.harney@signaldisplay.com

Printer friendly Cite/link Email Feedback | |

Author: | Haranas, Ioannis Iraklis; Harney, Michael |
---|---|

Publication: | Progress in Physics |

Geographic Code: | 1USA |

Date: | Jul 1, 2008 |

Words: | 3323 |

Previous Article: | A blind pilot: who is a super-luminal observer? |

Next Article: | Resolving spacecraft earth-flyby anomalies with measured light speed anisotropy. |

Topics: |