Do SR account for one sixth or one third of the anomalous perihelion shift?

Question

In answers to two other questions on this site:

Mercury's Orbital Precession in Special Relativity

and

Do 2-body elliptic orbits precess in special relativity?

it is stated that special relativity could account for one sixth of the "anomalous perihelion shift" of Mercury and other planets.

If I insert the "relativistic mass" (I know people hate that term) on both sides of Newtons equation for gravitational acceleration you first have:

$$\frac{d(m\gamma\bar{v})}{dt}=-\frac{GMm\gamma}{r^2}\hat{r}\tag1$$

and you end up with:

$$\frac{d\bar{v}}{dt}=-\frac{GM}{r^2}(\hat{r}-\frac{v^2}{c^2}(\hat{r}\cdot\hat{v})\hat{v})\tag2$$

($\hat{r}=\bar{r}/r, \hat{v}=\bar{v}/v, \gamma=1/\sqrt{1-v^2/c2}$)

This will give you more or less precisely one third of the experimentally found anomalous perihelion shift. (I checked this out by numerical integration, I do not have an analytical solution)

I guess it can be argued whether adding a "relativistic mass" in the form used in special relativity to the expression for gravitational acceleration by Newton constitutes "special relativity" or not.

Still using this method you get one third of the anomalous perihelion shift and not one sixth.

Questions:

1. Do SR account for one third or one sixth of the so called anomalous "perihelion shift"?

2. Can you say that one third of the weak field perihelion shift is explained by using the "Lorentz factor" in this way and two thirds are due to effects accounted for in general relativity but not in special relativity?

---

In the strong field limit you get orbits using the expression above that are obviously wrong, but still in the ballpark. The green circle represents the Schwarzschild radius and the red circle the innermost stable circular orbit. In this special case one gets three revolutions between consecutive aphelions.

---

According to both current answers by Elio Fabri and Pulsar assuming:

$$\frac{d(m\gamma\bar{v})}{dt}=-\frac{GMm}{r^2}\hat{r}\tag3$$

which results in:

$$\frac{d\bar{v}}{dt}=-\frac{1}{\gamma}\frac{GM}{r^2}(\hat{r}-\frac{v^2}{c^2}(\hat{r}\cdot\hat{v})\hat{v})\tag4$$

results in one sixth of the experimentally found perihelion shift and is the most SR way to do it based on arguments that you can get (4) but probably not (2) from a Lagrangian approach and that (4) but not (2) is consistent with a "conservative force". It could perhaps be argued that (3) violates the equivalence principle as you have "$m\gamma$" on the left side of the equation but only "$m$" on the right side. Note that using (4) we no longer have $v=\sqrt{GM/r}$ for circular orbits as is true in both classical mechanics and general relativity (in Schwarzschild coordinates using coordinate time at least) and using (2). Instead, for circular orbits, we end up with:

$v\sqrt{\gamma}=\sqrt{GM/r} \tag5$

It is a bit strange if you get the correct orbital velocity in both GR and classically but not in SR, if we assume (4) is the correct SR expression.

I plotted an example of a strong field orbit using (4) below. The extra factor of $1/\gamma$ appearing in (4) but not in (2) weakens the gravitational acceleration at high velocities and thus gives less perihelion shift. The kinds of orbits one gets using (4) instead of (2) looks less like what is expected from GR.

Answer 1

1. One sixth; see below.
2. Sort of, but not in a rigorous way. Special Relativity is only valid in flat spacetime, so it's not equipped to deal with gravitational systems. For a rigorous treatment you need the full machinery of GR, which has SR 'baked in', and you can't really disentangle the effects of curved spacetime and 'pure' SR in a clean way. Nonetheless, a purely SR calculation of Mercury's perihelion shift is a fun exercise.

With these caveats in mind, we can try to formulate some sort of SR version of Newton's laws. Your initial equation is almost correct, but if we want to keep the property that the gravitational force is conservative, the equation should have the form $$ \frac{\text{d}(m\gamma\vec{v})}{\text{d}t} = - \frac{GMm}{r^2}{\vec{e}_r}.\tag{1} $$ (again, SR is not able to describe gravity, but this is the best we can do). This way, the force can be written as the gradient of a gravitational potential $\vec{F} = -\vec\nabla V(r)$. A force that only depends on position and not on velocity also satisfies Newton's Third Law and conservation of momentum: keep in mind that the Sun feels an equal and opposite force from Mercury, but contrary to Mercury the velocity of the Sun is nearly zero. So if we want to keep the Newtonian form of gravity (i.e. no curvature of spacetime), the gravitational force can only depend on the distance between Mercury and the Sun.

We can verify this equation by deriving it from a relativistic energy that resembles the traditional Newtonian energy: $$ E = m\gamma c^2 - \frac{k}{r} = \sqrt{m^2c^4 + \vec{p}^2c^2} - \frac{k}{r},\tag{2} $$ with $k = GMm$ and $\vec{p} = m\gamma\vec{v}$. Conservation of energy $\dot{E}= 0$ then leads to $$ \frac{\vec{p}\cdot\dot{\vec{p}}c^2}{m\gamma c^2} = -\frac{k}{r^2}\dot{r}, $$ which indeed reduces to (1). We can derive the SR precession starting from (2), in a very similar way as the Newtonian Kepler problem (see e.g. Goldstein). First, we write the momentum in polar coordinates: $$ E = \sqrt{m^2c^4 + (m\gamma\dot{r})^2c^2 + (m\gamma r\dot{\varphi})^2c^2} - \frac{k}{r}.\tag{3} $$ Next, we introduce the SR angular momentum $$ L = m\gamma r^2\dot{\varphi}, $$ and we can use this to write $$ \dot{r} = \frac{\text{d}r}{\text{d}\varphi}\dot{\varphi} = \frac{\text{d}r}{\text{d}\varphi}\frac{L}{m\gamma r^2}. $$ Plugging this into (3), we obtain $$ \left(E + \frac{k}{r}\right)^2 = m^2c^4 + \frac{L^2c^2}{r^4}\left(\frac{\text{d}r}{\text{d}\varphi}\right)^2 + \frac{L^2c^2}{r^2}, $$ or $$ \frac{\text{d}r}{\text{d}\varphi} = \frac{r^2}{Lc}\sqrt{\left(E + \frac{k}{r}\right)^2 - m^2c^4 - \frac{L^2c^2}{r^2}}, $$ which leads to $$ \varphi = \varphi_0 + \int_{r_0}^{r}\frac{Lc\,\text{d}r}{r^2\sqrt{\left(E + \frac{k}{r}\right)^2 - m^2c^4 - \frac{L^2c^2}{r^2}}}. $$ Without loss of generality, we can define $\varphi_0=0$ at perihelion $r_0$. If we change the variable of integration to $u=1/r$, the integral takes the form $$ \varphi = - \int_{u_0}^{u}\frac{Lc\,\text{d}u}{\sqrt{E^2 - m^2c^4 + 2kEu + (k^2-L^2c^2)u^2}}. $$ Just like in the Newtonian case, this integral has an analytical solution (see Goldstein 3rd Ed, eq. (3.51)) $$ \int\frac{\text{d}x}{\sqrt{\alpha + \beta x + \gamma x^2}} = \frac{1}{\sqrt{-\gamma}}\arccos\left(-\frac{\beta + 2\gamma x}{\sqrt{\beta^2 -4\alpha\gamma}}\right). $$ We find $$ \varphi = -\frac{Lc}{\sqrt{L^2c^2-k^2}}\arccos\left(\frac{(L^2c^2-k^2)u - kE}{\sqrt{\Delta}}\right), $$ with $$ \Delta = k^2E^2 + (L^2c^2-k^2)(E^2 - m^2c^4), $$ so that finally $$ r = \frac{L^2c^2-k^2}{kE + \sqrt{\Delta}\cos\left(\sqrt{1-\frac{k^2}{L^2c^2}}\varphi\right)}. $$ This is indeed a precessing ellipse, with perihelion shift $$ \delta\varphi = \frac{2\pi}{\sqrt{1-\frac{k^2}{L^2c^2}}} - 2\pi \approx \frac{\pi k^2}{L^2c^2} $$ per orbit, using $(1-x)^{-1/2} \approx 1 + x/2$. The General Relativistic precession is (see wiki) $$ \delta\varphi_\text{GR} \approx \frac{6\pi G^2M^2m^2}{L^2c^2} = \frac{6\pi k^2}{L^2c^2}, $$ which is indeed about 6 times larger. With the values (see wiki and Mercury Fact Sheet ) $$ L/m = \frac{r_0v_0}{\sqrt{1-v_0^2/c^2}},\\ r_0 = 46\,001\,200\ \text{km},\\ v_0 = 58.98\ \text{km}\,\text{s}^{-1},\\ k/m = GM = 1.32712440042\times 10^{11}\ \text{km}^3 \text{s}^{-2}, $$ I obtain a perihelion shift of $0.0172''$ per orbit. With an orbital period of $87.969$ days, that amounts to $7.16''$ per century.

---

Suppose we started instead with $$ \frac{\text{d}(m\gamma\vec{v})}{\text{d}t} = - \frac{GMm\gamma}{r^2}{\vec{e}_r}. $$ It's easy to verify that the energy of this system is $$ E = m\gamma c^2\exp\left(-\frac{GM}{c^2r}\right) = \sqrt{m^2c^4 + \vec{p}^2c^2}\exp\left(-\frac{GM}{c^2r}\right), $$ which looks weird; it doesn't resemble either the Newtonian or the GR form. It does reduce to Newtonian physics for weak gravitational fields and low velocities, but this is not what Goldstein in the original questions had in mind. Ultimately though, it's a discussion about definitions, since there is no correct SR version of gravity.

Answer 2

$\let\g=\gamma \def\bp{{\bf p}} \def\br{{\bf r}} \def\bv{{\bf v}} \def\bF{{\bf F}} \def\D#1#2{{d#1 \over d#2}}$ Let me try to put some order in the matter. There are different problems involved which shouldn't be confused:

1. A general problem: how to deal with central forces in SR mechanics.
2. Some special examples, among them Sommerfeld's relativistic treatment of Hydrogen atom.
3. SR approach to the gravitational two body problem.

As to the first problem, it's a special case of a matter discussed in almost all books on classical mechanics, e.g. Goldstein or Landau.

Generally a variational (Lagrangian) approach is followed. The Lagrangian for a mass point subjected to a force is written (with proper time as independent variable): $$L = m\,c^2 - A_\mu u^\mu$$ where $A_\mu$ is a 4-potential and $u^\mu$ the 4-velocity. If in some inertial frame only $A_0$ doesn't vanish the action integral is written $$\int_{\tau_1}^{\tau_2}\!\![m\,c^2 - A_0 u^0]\,d\tau = \int_{t_1}^{t_2} \!\![m\,c^2 \g(\bv) - A_0(\br)]\,dt.\tag1$$ Lagrange equations are $$m\,\D{}t\,[\bv\,\g(\bv)] = -\nabla A_0(\br)$$ $$\D \bp t = \bF(\br).\tag2$$ From eq. (2) a differential equation for the trajectory (a plane curve) is derived following standard approach.

If $A_0=-k/r$ ($k>0$) the equation for trajectory can be solved in close terms and gives a "rotating" ellipse (1/6 the GR result). You find detailed calculations in your first link.

We only have to discuss to what physical situations the above treatment applies. There are no approximations and the only condition is that "potential" is the time component of a 4-vector. This ends discussion of item 1.

---

2) I don't remember how Sommerfeld approached the problem of relativistic Hydrogen atom, but I'd bet he closely followed the general treatment examined above. In Hydrogen atom the only interaction is between electron and proton. Given the high mass ratio (about 1 to 2000) the problem can be safely reduced to that of an electron moving in an electrostatic field of magnitude $e/r^2$. Potential energy is then $-e^2/r^2$.

It's well known that electrostatic potential is the time component of a 4-vector, so we are precisely in the domain where the above treatment can be applied with no further approximation. The rotating ellipse has a consequence: whereas angular momentum is still conserved, this is no longer true for Lenz vector. In QM terms, $l$-degeneracy of levels is broken. This is a satisfying feature, but the resulting "fine structure" doesn't agree with experimental results. We had to wait for Dirac's theory to obtain a good agreement with experiment.

---

3) IMO there's no hope to get the GR result for precession in SR. The main reason is that spacetime in GR is curved and only part of the effects due to curvature can be reproduced in SR. This is dramatically shown in your figure. There is a trajectory merrily entering and exiting the horizon - a thing perfectly understandable since in SR there is no horizon: spacetime is flat.

A comment about your introduction of a $\g$ in force expression. The issue isn't a more or less "hated" relativistic mass. Actually there could be some argument in favour of your $\g$, as in GR what causes spacetime curvature is not plainly mass, but the energy-momentum tensor. Then we could think that the way a body feels an existing curvature may depend on its energy, not on rest mass.

Some time ago I made a calculation (unpublished) aimed at seeing if acceleration of a falling body, as measured in a local rest frame in Schw. geometry depends on $m$ or on $m\g$. I tested two cases:

• vertical fall
• transverse motion (a projectile shot horizontally).

In both cases I found that the measured acceleration if interpreted within SR should be viewed as due to a force proportional to $m\g$.

Should that result be taken as a proof that your $\g$ factor is right? I'd not say so, as the context is entirely different. You're assuming SR and building a gravitational theory where a force exists, acting on the test body, depending on $m\g$. You take that expression as holding in the whole flat spacetime of SR.

My approach was different. I assumed GR and Schw. spacetime. I computed geodesics in that spacetime and asked myself how the corresponding motion of a test body would look if measured in a small neighborhood of a fixed point. Measured, I mean, with rulers and clocks for that neighborhood. Then I searched how a physicist in that neighborhood would interpret the observed motion if he used SR, included the modified second law: which force did cause the observed motion? I already gave the answer I found.

Answer 3

Newton’s gravitational force is $$F_G=-\frac{GMm}{r^3} \vec{r} \tag{1}$$ The mathematical definition of force (in all concepts, classically or not) is $$ \vec{F}=\frac{d \vec{p}}{dt} \tag{2}$$ Therefore, (1) can be written as $$\frac{d \vec{p}}{dt}=-\frac{GMm}{r^3} \vec{r} \tag{3}$$ The SR definition of momentum is $$\vec{p}=mγ \vec{v} \tag{4}$$ Combining (1) and (2), one gets $$\frac{d(mγ \vec{v})}{dt} =-\frac{GMm}{r^3} \vec{r} \tag{5}$$ From the answers of Pulsar and others it gives you one sixth.

Your confusion came from introducing the so called “relativistic mass”, which seems to give you the wrong result. The mass $m$ is invariant in all contexts and doesn’t transform as $m→mγ$. However, (4) is the best definition of momentum and can be considered valid even in classical physics as γ→1.