In a Yahoo discussion group, a question was raised as to whether the GEM unified field proposal was consistent with the precession of the perihelion of Mercury results. My stock reply is that the coefficients of the GEM metric are identical to those of the Schwarzschild metric only to first order PPN accuracy, the level used in the calculations, ergo the results must be identical. I like a short, solid reply.
I also like a long-winded one, because it shows all the nuts and bolts. I have read in many places about the precession of the perihelion of Mercury, yet didn't get how they actually did the darn calculation. There were always a few steps that I did not follow. While reading through the Sean Carroll's Lecture notes on GR, I decided to try and figure out the details. Here I write it all out. This is not easy or short, but for those willing to work at it, might be a unique information source.
Start with the GEM field equations with the gauge choice of a constant 4-potential, then for a spherically symmetric, non-rotating, uncharged mass, the metric which solves the field equations is:
A couple of procedural things. I like to baby step through derivations: if there is one division to make, I write it out, just like I would do with a pencil, but this time with TeX. Because I was initially trained to get the units right, I still keep that tradition going as an internal consistency check, instead of adopting natural units.
Let's make the metric simpler by making some reasonable assumptions about the
angle , namely that it is always in a plane at angle
so that
and
. Write out the metric with these
substitutions:
What one is suppose to do is work with the equation of motion, looking for a
constant of the equation of motion. Sometime though, it is quicker to cheat a
little, and just divide the above expression by to get here:
Let's get rid of the exponential in front of the
.
For our solar system, the exponent is really tiny, so we can use the Taylor series expansion to one term of the exponent. One slight trick is that the -4 one is like a -2 exponent squared. Approximate away:
Notice that this expression is not a function of either time or angle
. This means there is a conserved quantity associated with change in
the time (energy
) and a change in angle (angular momentum
). Come back
at a later time, the expression stays the same. Spin around a few degrees,
the metric stays the same. Something like this happens for Lagrange
densities. The same logic applies to metrics. This time the idea goes by a
different name: Killing vectors. I do not understand enough of the details,
sorry, but it turns out that a Killing vector works like this:
Do this for energy and angular momentum, normalizing to a test mass:
Calculate the squares of these:
Plug these back into the metric expression:
Let's pause and catch our breath, there still is a long way to go. All that
has happened is to introduce two constant quantities, the energy and
angular momentum
into the metric expression. The equation still displays
the ``bending factor'', the
, which is equal to
one for Newtonian gravity. We could stop here because this expression is
identical to the one for GR, but let's continue.
The next task is to do classical celestial mechanics. Talk about an area of
weakness for me! The way this is done is to work with the two invariants:
energy because Newtonian gravity is conservative, and angular momentum because
the system is isolated. The cause of the planets whirling around is a
potential, gravity. For that reason, we need to change the variable.
Plug these into the equation-that-used-to-be-the-metric:
Multiply this out for , and bring all the terms to one side of the
equation:
Take the derivative with respect to which will drop a few constant
terms.
Divide all the terms by
:
The first three terms are classical Newtonian gravitational physics. The fourth term is the correction required by GR and GEM. First we need to understand the Newtonian gravity solution, ignoring the correction. Here is the Newtonian expression to solve:
Guess a solution. Fortunately, this differential equation is a variation on a
common equation with a cosine solution. Thing is, the circle will be a tad
eccentric, so include the eccentricity :
Guess:
Put it all together:
At the perihelion,
. Now will can get an
expression for
at the perihelion with no corrections:
the perihelion distance is
for an ellipse with a
semi-major axis of
. Plug this in, and divide by
:
This will be useful, because in a short while we will see a
, and we can drop the
right in.
How do we now deal with the correction? Well the Newtonian solution is pretty
darn good, so we can use that as a place to start. The correction term has
, so calculate that:
Put this back into the differential equation:
Which of these three terms matters? The first one is a constant term, much
smaller than the constant term we already have since it has a factor of
. A small number cubed can be safely ignored. The second term
oscillates just like the Newtonian solution, so it will always be adding into
the results. This ``on resonance'' terms will continually add up so the term
matters. The squared cosine term will not be in tune with the dominant
Newtonian solution, so it will add a periodic perturbation to the path, a form
of noise. Keep only the middle term of
:
What is the solution for this differential equation? If we can find an expression that drops the cosine term, then that expression with the previous solution would do the trick. We will need a sine or cosine, and something to make terms drop.
Guess:
The right coefficients in front of this correction term will cancel the term. Combine the correction term with the Newtonian solution:
Cosine is an even function, sine is an odd function, and
is
an even function. Because
is so small, only the
first term of its Taylor series expansion is going to matter, and it will
behave like it subtracts a bit off of a cosine. Bring the correction into the
cosine, noticing a
is already shared:
)
How much does change? In one revolution it would go
. Calculate
the ratio of no correction to correction:
The advance is
. Now we can use the
result found earlier for the major axis and the eccentricity[
]:
This is what appears in all the books. I think this is tough to get all the details clear and correct. I hope I achieved that goal.
Home Page |
Quaternion Physics |
Pop Science
Java |
The Bike |
Lindy Hop |
Contact Doug
Copyright © 2005, doug <sweetser@alum.mit.edu> All rights reserved worldwide.