- #36
DrGreg
Science Advisor
Gold Member
- 2,474
- 2,083
Proper time along a geodesic
According to Woodhouse(1)(2), the equations of free fall (i.e. geodesics) in the Schwarzschild metric are:
together with the definition of [itex]d\tau[/itex] via the metric. E is energy and J is angular momentum, and both are determined by the initial conditions of motion.
The above is from a reputable textbook so must be true. What follows is my own working, so might not.
In the case of radial motion, [itex]d\theta/d\tau = d\phi/d\tau = 0[/itex], and the metric simplifies to
Eliminating [itex]dt/d\tau[/itex] between (1) and (3) gives
Feed that into your favourite online integrator. (Good luck, I had problems with one of sites I tried! )
If the particle is released from "rest" ([itex]dr/d\tau = 0[/itex]) at radius r3, then (5) tells you the value of E, and so (6) becomes
In particular if [itex]r_3 = \infty[/itex], E = 1 (which makes the integral easy enough to do without computer assistance).
Coordinate time t can be calculated by combining (1) and (5):
As for the distance, the simple answer is zero! The particle passes through both events so the distance between events in the particle's frame is zero. But I guess what you really mean is, what is the distance between two hovering buoys at the start and end points? Off the top of my head, I don't know, I'll have to think about that.
References
1. Woodhouse, N.M.J. (2007), General Relativity, Springer, London, ISBN 978-1-84628-486-1, pages 100, 107, 124.
2. Woodhouse, N.M.J. (2003), http://people.maths.ox.ac.uk/~nwoodh/gr/index.html - online lecture notes on which the book (1) was closely based, pages 54-55 (Lecture 12), p59 (L13), p73 (L15).
Online Integration results
Integrate[1/Sqrt[a^2 - 1 + (2*m)/x], x] ==
(Sqrt[-1 + a^2 + (2*m)/x]*x)/(-1 + a^2) - (m*Log[m + (-1 + a^2 + Sqrt[-1 + a^2]* Sqrt[-1 + a^2 + (2*m)/x])*x])/(-1 + a^2)^(3/2)
Integrate[1/Sqrt[(-2*m)/r3 + (2*m)/x], x] ==
-(r3*(2*Sqrt[m*(-r3^(-1) + x^(-1))]*x + Sqrt[m]*Sqrt[r3]*ArcTan[ (Sqrt[r3]*Sqrt[m*(-r3^(-1) + x^(-1))]* (-r3 + 2*x))/(2*Sqrt[m]*(-r3 + x))]))/ (2*Sqrt[2]*m)
I don't claim to be an expert in this subject.kev said:I wonder if DrGreg or George can provide us with an equation for the time measured by the clock of an observer falling from r2 to r1 and what she would consider the distance r2-r1 to be? For simplicity, consider only the motion of a observer falling from infinity with an inital velocity of zero.
According to Woodhouse(1)(2), the equations of free fall (i.e. geodesics) in the Schwarzschild metric are:
[tex] E = \left(1 - \frac{r_s}{r} \right) \frac{dt}{d\tau} = \mbox{constant} [/tex] ...(1)
[tex] J = r^2 \sin^2 \theta \frac{d\phi}{d\tau} = \mbox{constant} [/tex] ...(2)
[tex] J = r^2 \sin^2 \theta \frac{d\phi}{d\tau} = \mbox{constant} [/tex] ...(2)
together with the definition of [itex]d\tau[/itex] via the metric. E is energy and J is angular momentum, and both are determined by the initial conditions of motion.
The above is from a reputable textbook so must be true. What follows is my own working, so might not.
In the case of radial motion, [itex]d\theta/d\tau = d\phi/d\tau = 0[/itex], and the metric simplifies to
[tex] \left(1 - \frac{r_s}{r} \right) \left(\frac{dt}{d\tau}\right)^2 - \frac{1}{1 - r_s/r}\left(\frac{dr}{d\tau}\right)^2 = 1[/tex] ...(3)
Eliminating [itex]dt/d\tau[/itex] between (1) and (3) gives
[tex] \frac {E^2 - (dr/d\tau)^2}{1 - r_s/r} = 1[/tex] ...(4)
[tex]\left(\frac{dr}{d\tau}\right)^2 = E^2 - 1 + \frac{r_s}{r}[/tex] ...(5)
[tex]\int_{r_1}^{r_2}\frac{dr}{\sqrt{E^2 - 1 + r_s/r}} = \tau[/tex] ...(6)
[tex]\left(\frac{dr}{d\tau}\right)^2 = E^2 - 1 + \frac{r_s}{r}[/tex] ...(5)
[tex]\int_{r_1}^{r_2}\frac{dr}{\sqrt{E^2 - 1 + r_s/r}} = \tau[/tex] ...(6)
Feed that into your favourite online integrator. (Good luck, I had problems with one of sites I tried! )
If the particle is released from "rest" ([itex]dr/d\tau = 0[/itex]) at radius r3, then (5) tells you the value of E, and so (6) becomes
[tex]\int_{r_1}^{r_2}\frac{dr}{\sqrt{r_s/r - r_s/r_3}} = \tau[/tex] ...(7)
In particular if [itex]r_3 = \infty[/itex], E = 1 (which makes the integral easy enough to do without computer assistance).
Coordinate time t can be calculated by combining (1) and (5):
[tex] t = \int \frac{dt}{d\tau} d\tau = \int \frac{dt/d\tau}{dr/d\tau} dr [/tex] ...(8)
As for the distance, the simple answer is zero! The particle passes through both events so the distance between events in the particle's frame is zero. But I guess what you really mean is, what is the distance between two hovering buoys at the start and end points? Off the top of my head, I don't know, I'll have to think about that.
References
1. Woodhouse, N.M.J. (2007), General Relativity, Springer, London, ISBN 978-1-84628-486-1, pages 100, 107, 124.
2. Woodhouse, N.M.J. (2003), http://people.maths.ox.ac.uk/~nwoodh/gr/index.html - online lecture notes on which the book (1) was closely based, pages 54-55 (Lecture 12), p59 (L13), p73 (L15).
Online Integration results
Integrate[1/Sqrt[a^2 - 1 + (2*m)/x], x] ==
(Sqrt[-1 + a^2 + (2*m)/x]*x)/(-1 + a^2) - (m*Log[m + (-1 + a^2 + Sqrt[-1 + a^2]* Sqrt[-1 + a^2 + (2*m)/x])*x])/(-1 + a^2)^(3/2)
Integrate[1/Sqrt[(-2*m)/r3 + (2*m)/x], x] ==
-(r3*(2*Sqrt[m*(-r3^(-1) + x^(-1))]*x + Sqrt[m]*Sqrt[r3]*ArcTan[ (Sqrt[r3]*Sqrt[m*(-r3^(-1) + x^(-1))]* (-r3 + 2*x))/(2*Sqrt[m]*(-r3 + x))]))/ (2*Sqrt[2]*m)
Last edited by a moderator: