- #71
- 10,338
- 1,516
Nice reformulation! - it's less abstract and more physical.
I'll take this opportunity to provide what I think is a simpler derivation of the original formula.
First, a very quick review of 4-velocities (I know George and Chris know all this, but I'm not sure of how much our other readers know).
If we have coordinates T, r, [itex]\theta[/itex], and [itex]\phi[/itex], we can write out the components of the 4-velocity as
[tex]
u^a = \left[ \frac{dT}{d\tau}, \frac{dr}{d\tau}, \frac{d\theta}{d\tau}, \frac{d\phi}{d\tau} \right] [/tex]
using geometric units where c=1. [tex]\tau[/itex] should be interpreted as the proper time of the observer, i.e. we can parameterize the wordline of the observer by four functions [tex]T(\tau),r(\tau),\theta(\tau),\phi(\tau)[/tex] where [itex]\tau[/itex] represents the proper time or 'age' of the observer, and the above set of derivatives of the coordinates with respect to the proper time is the 4-velocity.
It should also be apparent that
[tex]g_{ab} u^a u^b [/tex] has a value of plus or minus 1, +1 if we have a +--- metric signature, -1 if we have a -+++ metric signature.
Now let's consider how to find the magnitude of the relative 3-velocity of two observers if one observer has a 4-velocity of [itex]u^a[/itex] and the other has a 4-velocity of [itex]v^b[/itex].
The product [tex]\vec{u} \cdot \vec{v} = g_{ab} u^a v^b[/tex] is going to be independent of the choice of frame or coordinate system. It will determine the angle between the two 4-vectors on a space-time diagram, and hence it will determine the magnitude of the 3-velocity.
By considering the case in flat space-time where u = (1,0,0,0) and v = ( [itex]1/\sqrt{1-\beta^2},\beta/\sqrt{1-\beta^2},0,0)[/itex], we can see that the magnitude of the 3-velocity is just
[tex]\beta = \sqrt{1 - \left(\frac{1}{\vec{u}\cdot\vec{v}}\right)^2}}[/tex]
To apply this formula to the original problem, we need only the 4-velocity of a stationary observer at radius r, and the 4-velocity of an infalling observer with energy E at radius r.
In Schwarzschild coordinates, using a +--- metric, the 4-velocity of a stationary observer will just be
[tex]\left[ \frac{1}{\sqrt{1-2M/r}},0,0,0 \right] [/tex]
This results from the normalization requirement that [itex]g_{00} u^a u^b = 1[/itex]
We only need the first component, [itex]v^0[/itex] of the 4-velocity of an infalling observer to compute the dot product.
http://www.fourmilab.ch/gravitation/orbits/, or the argument I presented about finding conserved quantites from Killing vectors, gives the value of this quantity.
[tex]v^0 = \frac{E}{1-2M/r}[/tex]
Therfore the dot product [itex]\vec{u} \cdot \vec{v} [/itex] is just
[tex]g_{00} \left( \frac{1}{\sqrt{1-2M/r}} \right) \left( \frac{E}{1-2M/r} \right) = \frac{E}{\sqrt{1-2M/r}} [/tex]
since [itex]g_{00}[/itex] = 1-2M/r using the +--- sign convention.
So the relative velocity is just
[tex]\beta = \sqrt{1 - \frac{1-2M/r}{E^2}}[/tex]
We can also re-write this in the form suggested by George by solving for the energy E in terms of the maximum height.
I'll take this opportunity to provide what I think is a simpler derivation of the original formula.
First, a very quick review of 4-velocities (I know George and Chris know all this, but I'm not sure of how much our other readers know).
If we have coordinates T, r, [itex]\theta[/itex], and [itex]\phi[/itex], we can write out the components of the 4-velocity as
[tex]
u^a = \left[ \frac{dT}{d\tau}, \frac{dr}{d\tau}, \frac{d\theta}{d\tau}, \frac{d\phi}{d\tau} \right] [/tex]
using geometric units where c=1. [tex]\tau[/itex] should be interpreted as the proper time of the observer, i.e. we can parameterize the wordline of the observer by four functions [tex]T(\tau),r(\tau),\theta(\tau),\phi(\tau)[/tex] where [itex]\tau[/itex] represents the proper time or 'age' of the observer, and the above set of derivatives of the coordinates with respect to the proper time is the 4-velocity.
It should also be apparent that
[tex]g_{ab} u^a u^b [/tex] has a value of plus or minus 1, +1 if we have a +--- metric signature, -1 if we have a -+++ metric signature.
Now let's consider how to find the magnitude of the relative 3-velocity of two observers if one observer has a 4-velocity of [itex]u^a[/itex] and the other has a 4-velocity of [itex]v^b[/itex].
The product [tex]\vec{u} \cdot \vec{v} = g_{ab} u^a v^b[/tex] is going to be independent of the choice of frame or coordinate system. It will determine the angle between the two 4-vectors on a space-time diagram, and hence it will determine the magnitude of the 3-velocity.
By considering the case in flat space-time where u = (1,0,0,0) and v = ( [itex]1/\sqrt{1-\beta^2},\beta/\sqrt{1-\beta^2},0,0)[/itex], we can see that the magnitude of the 3-velocity is just
[tex]\beta = \sqrt{1 - \left(\frac{1}{\vec{u}\cdot\vec{v}}\right)^2}}[/tex]
To apply this formula to the original problem, we need only the 4-velocity of a stationary observer at radius r, and the 4-velocity of an infalling observer with energy E at radius r.
In Schwarzschild coordinates, using a +--- metric, the 4-velocity of a stationary observer will just be
[tex]\left[ \frac{1}{\sqrt{1-2M/r}},0,0,0 \right] [/tex]
This results from the normalization requirement that [itex]g_{00} u^a u^b = 1[/itex]
We only need the first component, [itex]v^0[/itex] of the 4-velocity of an infalling observer to compute the dot product.
http://www.fourmilab.ch/gravitation/orbits/, or the argument I presented about finding conserved quantites from Killing vectors, gives the value of this quantity.
[tex]v^0 = \frac{E}{1-2M/r}[/tex]
Therfore the dot product [itex]\vec{u} \cdot \vec{v} [/itex] is just
[tex]g_{00} \left( \frac{1}{\sqrt{1-2M/r}} \right) \left( \frac{E}{1-2M/r} \right) = \frac{E}{\sqrt{1-2M/r}} [/tex]
since [itex]g_{00}[/itex] = 1-2M/r using the +--- sign convention.
So the relative velocity is just
[tex]\beta = \sqrt{1 - \frac{1-2M/r}{E^2}}[/tex]
We can also re-write this in the form suggested by George by solving for the energy E in terms of the maximum height.
Last edited: