__On
the Spinning Sun-Engendered Mercury Orbital Precession of .43 Arc Seconds per
Year
__

(2) THE FEYNMAN LECTURES ON PHYSICS, volume 2, Section 14-5.

In this article we compute the orbit of Mercury as it travels around the Sun.
Two cases are examined: (1) with the Sun not spinning (Newton 2 and his force
law of Universal Gravitation), and (2) with the Sun spinning and (according to
GMT) engendering a gravitomagnetic field in addition to the omnipresent
gravitational field.

In GMT the spinning sun is theorized to engender a dipolar gravitomagnetic field
that (a) points in the positive z-direction at internal points of the Sun’s
equatorial plane, and (b) points in the negative z-direction at points in the
external equatorial plane. The orbit of Mercury is also assumed to lie in the xy
plane and, from the perspective of the positive z axis, also to be CCW. In Case
1 Mercury is theorized to be acted upon solely by a gravitational force. In
Case 2 Mercury is theorized to be acted upon by a gravitational force *and* by an opposing gravitomagnetic force.

*open*
ellipse.

It is expected that the computed result may not be precisely what is observed
owing to certain realities. First the spinning Sun is modeled as a spinning
square line mass, but in reality it is a spinning, oblate sphere of mass with a
complicated internal density and with zones that rotate at different rates.
Secondly, the Sun's spin equatorial plane does not actually coincide with the
plane of Mercury's orbit. (Mercury's orbital plane is tilted slightly, relative
to the Sun's spin equatorial plane.) Finally there may be small numeric computer
errors that cannot be ignored, since they are comparable in size to the very
faint gravitomagnetic effects.

Feynman(Reference 2) models a square electric current loop as a square line
charge q, with sides a, and line charge density ^{2},
or m=qwa^{2}/2p.
For a point r>a in the loop's plane he finds that B points in the negative z
direction and has the value B_{z}=-m/(4pe_{o}c^{2}r^{3})
=^{2}/2p(4pe_{o})c^{2}r^{3}.

We *model* the spinning Sun as (a) a
square line mass M, where M is the mass of the Sun, and (b) with side a equal to
R, the radius of the Sun, and (c) with a line mass density of l=M/4R). We assume that this line mass circulates CCW with a mass current
of I=lv,
where v=2wR/p is the line mass speed around the loop.

At points in the equatorial plane, and at distances r>R, our model indicates
that O_{z}=-mG/c^{2}r^{3}
= -MwR^{2}G/2pc^{2}r^{3}.
Given the observed complexities of the real Sun, this model-based value for O_{z}
does not produce an orbit whose precession matches the observed one of .43 arc
seconds per year. Since the outer regions of the Sun are less dense than the
more central ones, it might be guessed that the value of O_{z} in the
real case must be somewhat less than what our model indicates. Indeed one can
experiment with different adjustment values for O_{z} to produce a
precession that approximates the observed .43 arc seconds. It turns out that an
adjustment factor of .959 results in an O_{z} that closely produces the
required .43 arc seconds per year.

__Appendix. Python Program.__

FY=gY*MercuryMass