Bill_K said:
This agrees with the result in my
blog post,
Ω = γ
2ω[r - 3M(1 - aω)]/r + γ
2Ma(1 - aω)
2/r
3
γ
2 = (1 - (r
2 + a
2)ω
2 - 2M(1 - aω
)2/r)
-1
in the nonrotating case, ω = 0.
I am happy to report that Bill's equation also agrees with the result by MTW for the ZAMO case, the geodesic orbit case by Rindler and is algebraically equivalent to the
general equatorial equation (35) by Vishveshwara, when the gamma factor is fixed as mentioned
in post 57 of the parallel thread.
With some tutoring on notation from WBN, this is how I obtained the result for the ZAMO case from
the MTW precession equation:
##\Omega = \frac{1}{2}\sqrt{\frac{g_{\phi\phi}}{|g_{t\phi} - \omega^2 g_{\phi\phi}|}} \left[\frac{\omega_{,\theta}}{p}\epsilon_{\hat r}-\frac{\Delta^{1/2}\omega_{,r}}{p}\epsilon_{\hat \theta}\right]##
where ##\Delta = (r^2-2mr+a^2)##. By considering only ZAMOs in the equatorial plane we can set ##\omega_{,\theta}=0##,##p=r## and after inserting the ZAMO orbital velocity ##\omega^2 = (-g_{t\phi}/ g_{\phi\phi})^2## I obtain after some simplification:
##\Omega = \frac{1}{2}\frac{g_{\phi\phi}}{\Delta^{1/2}} \left[-\frac{\Delta^{1/2}\omega_{,r}}{r}\epsilon_{\hat \theta}\right]##
Interpreting ##\omega_{,r}## as ##\frac{\partial{}}{\partial{r}}(-\frac{g_{t\phi}}{ g_{\phi\phi}})## and after carrying out the derivative the result is:
##\Omega = \frac{1}{2}\frac{g_{\phi\phi}}{\Delta^{1/2}} \left[-\frac{\Delta^{1/2}}{r}\left(\frac{-2ma(a^2+3r^2)}{r^2 (g_{\phi\phi})^2}\right)\epsilon_{\hat \theta}\right]##
##\Omega = \frac{ma(a^2+3r^2)}{r^3 (r^2+a^2+2ma^2/r)}\epsilon_{\hat \theta}##
where ##\epsilon_{\hat \theta}## is the unit vector indicating the axis of the resulting precession is orthogonal to the equatorial plane.
The same result (except for the sign) is obtained by inserting the ZAMO velocity into the general equatorial equation (35) by Vishveshwara or Bill's equation. The general equation also agrees with the result of ##\sqrt{m/r^3}## given by Rindler when we enter the geodesic orbital angular velocity ##\omega = 1/(a\pm\sqrt{r^3}{m})## also given by Rindler. None of these results are particularly easy to obtain from the general formulas and required the use of Maple and some manual simplification.
Bill's general equation which is slightly simpler can be expressed as:
##\Omega = \omega\gamma^2 - \omega\frac{3m}{r}(1-a\omega)\gamma^2 +\frac{ma}{r^3}(1-a\omega)^2\gamma^2##
where the first term is the kinematic Thomas precession, the second term is the geodetic precession due the curvature of spacetime around a gravitational body and the third term is the Lense-Thirring precession due the angular velocity (a) of a rotating gravitational body.
##\Omega## in all the above formulas is the rotation rate of a rigid rotating lattice relative to a set of gyroscopes as measured by an observer at rest with (and co-spinning with) the lattice, so the precession rate of the gyroscopes is actually ##-Omega##. The rotation rate ##\omega## of the rigid lattice is measured by the coordinate observer at infinity in the respective metric. The rotating lattice is equivalent to a congruence of observers all with identical angular velocity, (not momentum). To obtain the precession rate of the gyroscopes as measured by the observer at infinity we correct ##\Omega## by the time dilation factor and subtract the result from ##\omega## so that:
##\Omega_{\infty} = \omega -\omega\gamma + \omega\frac{3m}{r}(1-a\omega)\gamma -\frac{ma}{r^3}(1-a\omega)^2\gamma##
Note that setting a=0 gives the expected Schwarzschild precession and setting m=0 gives the expected Minkowski precession. This strongly suggests that the Thomas precession given by the second term is valid in a gravitational field and has the same sign as in flat space, contrary to some claims in the literature.
I have attached a maple worksheet that compares Bill's and Vishveshwara's precession formulas for the ZAMO case.
The hardest part is keeping track of the signs to determine the direction of precession and a lot of the literature is vague on this issue using ##\pm## signs or ##\Omega^2## to duck the issue. The NASA handouts for the direction for the Lense-Thirring precession in the GBP are full of contradictions on the direction. At the moment on balance, based on what is commonaly assumed for the Schwarzschild case, I am going with the L-T precession being retrograde. Also the MTW equation gives the opposite sign for the precession compared to the other equations. If anyone can shed some light on the direction of precession, I would be grateful.