Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

Some general formulae for circular orbits in symmetric spacetimes

  1. Apr 17, 2013 #1

    pervect

    User Avatar
    Staff Emeritus
    Science Advisor

    Consider the equatorial plane of a spherically symmetric space-time. Then we can write the metric in the equatorial plane (theta=pi/2) in terms of three coordinates - [t, r, phi]

    ds^2 = -f(r) dt^2 + g(r) dr^2 + h(r) dphi^2

    For the Schwarzschild metric we can write:
    [tex]f = c^2(1 -\frac{2 G M}{c^2 r}) \quad g = 1 / (1 -\frac{2 G M}{c^2 r}) \quad h = r^2[/tex]

    For the isotropic Schwarzschild metric we can write:
    [tex]f = c^2 \left( \frac{1-\frac{GM}{2 c^2 r}}{1+\frac{GM}{2 c^2 r}} \right)^2 \quad g = (1 + \frac{GM}{2 c^2 r})^4 \quad h = r^2 (1 + \frac{GM}{2 c^2 r})^4 [/tex]

    For the non-cartesian PPN apprxomiation to isotropic Schwarzschild we can write:
    [tex]f = c^2 \left( 1 - \frac{2GM}{c^2 r} + \frac{2 G^2 M^2}{c^4 r^2} \right) \quad g = \left( 1 + \frac{2GM}{c^2 r} \right) \quad h = r^2 \left( 1 + \frac{2GM}{c^2 r} \right)
    [/tex]

    The geodesic equation for r is
    [tex]
    \frac{d^2r}{d \tau^2} + \frac{(\frac{df}{dr})}{ 2g} (\frac{d^2t}{d \tau^2}) + \frac{(\frac{dg}{dr})}{ 2g} (\frac{d^2 r}{d \tau^2}) + \frac{(\frac{dh}{dr})}{ 2g} (\frac{d^2\phi }{d \tau^2}) = 0
    [/tex]

    Various useful quantities are (the answers all appear to be sensible to me in Schwarzschild coordinates):

    The magnitude of the proper acceleration of a static observer (constant r) as measured with an onboard accelerometer

    [tex]c^2 \frac{(\frac{d f}{d r})}{2 f \sqrt{g}}[/tex]

    The coordinate acceleration (d^2 r / dt^2) of a geodesic observer who starts out at rest, (dr/dt=0, dphi/dt=0):

    [tex]-\frac{(\frac{d f}{d r})}{2g}[/tex]

    The coordinate angular velocity [itex]d \phi / dt[/itex]
    [tex]\sqrt{\frac{(\frac{d f}{d r})}{(\frac{d h}{d r})}}[/tex]

    The proper angular velocity [itex]d \phi / d \tau[/itex]

    [tex]c \frac{(\frac{d f}{d r}) } { f (\frac{d h}{d r}) - h (\frac{d f}{d r})} [/tex]

    coordinate and proper orbital periods - 2 pi over the respective angular velociteis

    The orbital period as measured by a co-located static observer
    *correction*
    [tex] \left( \frac{2 \pi}{c} \right){ \sqrt{f (\frac{d h}{d r}) / (\frac{d f}{d r}) } }[/tex]

    circumference of the circle with coordinate values r
    [tex]2 \pi \sqrt{h} [/tex]

    orbital velocity measured by a colocated static observer
    [tex]c \sqrt{ \frac{h (\frac{d f}{d r}) }{ f (\frac{d h}{d r})}}[/tex]

    The concept of velocity from "remote points of view" is not well-defined in GR, due to the problems of parallel tranpsorting velocities.

    http://math.ucr.edu/home/baez/einstein/node2.html

    though various ad-hoc approaches are possible.

    Note: The factors of c are annoyig to keep tract of - I took the simple approach of making sure the units come out correctly.
     
    Last edited: Apr 17, 2013
  2. jcsd
  3. Apr 17, 2013 #2

    PAllen

    User Avatar
    Science Advisor
    Gold Member

    Nice. But .. what do you mean by coordinate acceleration of a static observer? Even for circular orbit shouldn't dr/dt and the second derivative both be zero?
     
  4. Apr 17, 2013 #3

    Mentz114

    User Avatar
    Gold Member

    Useful. There is a small typo - the equatorial plane is ##\theta=\pi/2##, I think.
     
  5. Apr 17, 2013 #4

    pervect

    User Avatar
    Staff Emeritus
    Science Advisor

    That was wrongly worded - I meant the coordinate acceleration of a geodesic observer who was initially at rest. I'll fix it up, if I can.

    And add the missing minus sign...
     
    Last edited: Apr 17, 2013
  6. Apr 17, 2013 #5

    pervect

    User Avatar
    Staff Emeritus
    Science Advisor

    Fixed - thanks
     
  7. Apr 17, 2013 #6

    pervect

    User Avatar
    Staff Emeritus
    Science Advisor

    Some sanity check values, for the Schwarzschild metric:

    proper acceleration of a static observer:
    [tex]\frac{GM}{r^2 \sqrt{1-2GM/c^2r}}[/tex]

    notes: matches wiki

    coordinate acceleration of a geodesic observer who was initially at rest
    [tex]
    \frac{GM}{r^2} \left(1 - \frac{2GM}{c^2 r} \right)[/tex]

    coordinate angular velocity [itex]d \phi / dt[/itex]
    [tex]\sqrt{\frac{GM}{r^3}}[/tex]

    (notes: matches Newtonian)

    proper angular velocity [itex]d \phi / d \tau[/itex]
    [tex] \sqrt{\frac{GM}{r^2 \left(r - 3GM/c^2\right) }}[/tex]

    notes: goes to infinity at photon sphere

    orbital period measured by co-located static observer
    (correction!)

    [tex]2 \pi \sqrt{1-2GM/c^2r)}\sqrt{\frac{r^3}{GM}}[/tex]

    notes: time dilated from coordinate period, goes to infinity at the event horizon

    orbital velocity measured by co-located static observer
    [tex]\sqrt{\frac{GM}{r-2GM/c^2}}[/tex]

    matches Wiki
     
  8. Apr 17, 2013 #7

    pervect

    User Avatar
    Staff Emeritus
    Science Advisor

    One more correction, for the orbital period of a static observer, the time dilation was done wrongly. (It's fixed)

    The results seem to match wiki's http://en.wikipedia.org/w/index.php?title=Circular_orbit&oldid=540533714

    for the Schwarzschild metric. It's non-intuitive that [itex]d \phi / dt[/itex] remains finite at the event horizon (this scared me for a bit!) but this seems to match wiki's results.

    I should add that circular orbits are impossible below the photon sphere at r=3GM/c^2, so applying it at this altitude is wrong in any event, even if the results don't appear to "blow up".
     
    Last edited: Apr 17, 2013
  9. Apr 17, 2013 #8

    Mentz114

    User Avatar
    Gold Member

    I'm not sure that ##d \phi / dt## represents anything physically meaningful. It is an approximation that would break down anywhere near the horizon, like r < 3m. There's no doubt though that ##(d\phi/d\tau)/(dt/d\tau)## gives the Newtonian value.
     
  10. Apr 17, 2013 #9

    pervect

    User Avatar
    Staff Emeritus
    Science Advisor

    Another correction, but this one is too late for me to fix. Perhaps someone else can do it.

    The geodesic equation currently reads

    [tex]
    \frac{d^2r}{d \tau^2} + \frac{(\frac{df}{dr})}{ 2g} (\frac{d^2t}{d \tau^2}) + \frac{(\frac{dg}{dr})}{ 2g} (\frac{d^2 r}{d \tau^2}) + \frac{(\frac{dh}{dr})}{ 2g} (\frac{d^2\phi }{d \tau^2}) = 0
    [/tex]

    it should really read (fixes sign error and some typos)

    [tex]
    \frac{d^2r}{d \tau^2} + \frac{(\frac{df}{dr})}{ 2g} \left(\frac{dt}{d \tau} \right) ^2 + \frac{(\frac{dg}{dr})}{ 2g} \left(\frac{d r}{d \tau} \right)^2 - \frac{(\frac{dh}{dr})}{ 2g} \left(\frac{d \phi }{d \tau} \right)^2 = 0
    [/tex]

    Symbolically, this is just the geodesic equation, with the non-zero Christoffel symbols for the given metric written out explicitly. The symbolic version is:

    [tex]
    \frac{d^2r}{d \tau^2} + \Gamma^r{}_{tt} \left( \frac{dt}{d \tau} \right) ^2 + \Gamma^r{}_{rr} \left( \frac{d r}{d \tau} \right)^2 + \Gamma^r{}_{\phi \phi} \left( \frac{d \phi }{d \tau} \right)^2 = 0
    [/tex]

    thus
    [tex]\Gamma^r{}_{tt} = \frac{(\frac{df}{dr})}{ 2g} \quad \Gamma^r{}_{rr} = \frac{(\frac{dg}{dr})}{ 2g} \quad \Gamma^r{}_{\phi \phi} = - \frac{(\frac{dh}{dr})}{ 2g} [/tex]
     
    Last edited: Apr 18, 2013
  11. Apr 19, 2013 #10
    Very nice. I could probably rearrange the terms of the second expression of the last post above and get the expression myself, but it would be useful to know the expressions for ##dt/d\tau## for isotropic Schwarzschild and for the low order expansion of isotropic Schwarzschild (the PPN- approximation). The expression for the velocity of light in coordinate time in the PPN-approximation would be interesting to have to. For some reason people seem to "pick and choose", take the velocity of light and the proper time from anisotropic coordinates, but get the orbits from isotropic coordinates, which seems a bit odd to me.
     
  12. Apr 19, 2013 #11

    pervect

    User Avatar
    Staff Emeritus
    Science Advisor

    For a circular orbit

    *note missing square root in previously posted expression*

    [tex]\frac{d \phi}{d \tau} = c \sqrt{\frac{(\frac{d f}{d r}) } { f (\frac{d h}{d r}) - h (\frac{d f}{d r})}}[/tex]
    [tex]\frac{d t}{d \tau} = c \sqrt{\frac{(\frac{d h}{d r}) } { f (\frac{d h}{d r}) - h (\frac{d f}{d r})}}[/tex]

    Utilizing the facts that [itex]\frac{dr}{d \tau} = \frac{d^2r}{d \tau^2} = 0 [/itex] for a circular orbit, this comes from solving the following two equations.

    [tex]
    (1) \quad -f \left( \frac{d t}{d \tau} \right)^2 + h \left( \frac{d \phi}{d \tau} \right) ^2 = -c^2[/tex]



    [tex] (2) \quad \frac{(\frac{df}{dr})}{ 2g} \left(\frac{dt}{d \tau} \right) ^2 - \frac{(\frac{dh}{dr})}{ 2g} \left(\frac{d \phi }{d \tau} \right)^2 = 0[/tex]

    If you multiply both sides of eq (1) by [itex]d \tau[/itex]^2, you can see this eq (1) is just the metric equation.

    Eq (2) is the geodesic equaton for r, symbolically:

    [tex]\Gamma^r{}_{tt} \left( \frac{dt}{d \tau} \right) ^2 +\Gamma^r{}_{\phi \phi} \left( \frac{d \phi }{d \tau} \right)^2 = 0[/tex]

    The equation of light can be solved for by simply setting ds^2 = 0 in the metric - rearranging this will give the coordinate velocity .

    In the radial direction (dr/dt)^2 = f/g. In the [itex]\phi[/itex] direction, (d[itex]\phi[/itex]/dt)^2 = f/h.
     
  13. Apr 19, 2013 #12

    pervect

    User Avatar
    Staff Emeritus
    Science Advisor

    *corrected version*
    Consider the equatorial plane of a spherically symmetric space-time. Then we can write the metric in the equatorial plane (theta=pi/2) in terms of three coordinates - [t, r, [itex]\phi[/itex]]

    [tex]
    ds^2 = -f(r)\, dt^2 + g(r)\, dr^2 + h(r)\, d\phi^2
    [/tex]

    For the Schwarzschild metric we can write:
    [tex]f = c^2(1 -\frac{2 G M}{c^2 r}) \quad g = 1 / (1 -\frac{2 G M}{c^2 r}) \quad h = r^2[/tex]

    For the isotropic Schwarzschild metric we can write:
    [tex]f = c^2 \left( \frac{1-\frac{GM}{2 c^2 r}}{1+\frac{GM}{2 c^2 r}} \right)^2 \quad g = (1 + \frac{GM}{2 c^2 r})^4 \quad h = r^2 (1 + \frac{GM}{2 c^2 r})^4 [/tex]

    For the non-cartesian PPN apprxomiation to isotropic Schwarzschild we can write:
    [tex]f = c^2 \left( 1 - \frac{2GM}{c^2 r} + \frac{2 G^2 M^2}{c^4 r^2} \right) \quad g = \left( 1 + \frac{2GM}{c^2 r} \right) \quad h = r^2 \left( 1 + \frac{2GM}{c^2 r} \right)
    [/tex]

    We will need only one additional equation (other than the metric) to solve for circular orbits - the r component of the geodesic equation:

    [tex]
    \frac{d^2r}{d \tau^2} + \Gamma^r{}_{tt} \left( \frac{dt}{d \tau} \right) ^2 + \Gamma^r{}_{rr} \left( \frac{d r}{d \tau} \right)^2 + \Gamma^r{}_{\phi \phi} \left( \frac{d \phi }{d \tau} \right)^2 = 0
    [/tex]

    Computing the values of the Christoffel symbols:
    [tex]\Gamma^r{}_{tt} = \frac{(\frac{df}{dr})}{ 2g} \quad \Gamma^r{}_{rr} = \frac{(\frac{dg}{dr})}{ 2g} \quad \Gamma^r{}_{\phi \phi} = - \frac{(\frac{dh}{dr})}{ 2g} [/tex]

    we write:

    [tex]
    \frac{d^2r}{d \tau^2} + \frac{(\frac{df}{dr})}{ 2g} \left(\frac{dt}{d \tau} \right) ^2 + \frac{(\frac{dg}{dr})}{ 2g} \left(\frac{d r}{d \tau} \right)^2 - \frac{(\frac{dh}{dr})}{ 2g} \left(\frac{d \phi }{d \tau} \right)^2 = 0
    [/tex]

    Note that the first and second derivatives of r with respect to time are zero for a circular orbit. Additionally, [itex]d^2 t / d \tau^2 = d^2 \phi / d \tau^2 = 0[/itex], so [itex]dt/d\tau[/itex] and [itex]d\phi / d\tau[/itex] are constants of motion for a circular orbit.

    The metric equation implies that

    [tex]-f(r) \, \left( \frac{dt}{d\tau} \right) ^2
    + g(r) \, \left( \frac{dr}{d\tau} \right)^2 + h(r) \, \left( \frac{d \phi}{d\tau} \right)^2 = -c^2
    [/tex]

    This, combined with the geodesic equation, enables us to solve for all quantities of interest.


    Various useful quantities:

    The magnitude of the proper acceleration of a static observer (constant coordinates) as measured with an onboard accelerometer

    [tex]c^2 \frac{(\frac{d f}{d r})}{2 f \sqrt{g}}[/tex]

    The coordinate acceleration (d^2 r / dt^2) of a geodesic observer who starts out at rest, (dr/dt=0,[itex] d\phi/dt=0[/itex]) is

    [tex]-\frac{(\frac{d f}{d r})}{2g}[/tex]

    The coordinate angular velocity [itex]d \phi / dt[/itex]
    [tex]\sqrt{\frac{(\frac{d f}{d r})}{(\frac{d h}{d r})}}[/tex]

    The proper angular velocity [itex]d \phi / d \tau[/itex]
    [tex]\frac{d \phi}{d \tau} = c \sqrt{\frac{(\frac{d f}{d r}) } { f (\frac{d h}{d r}) - h (\frac{d f}{d r})}}[/tex]

    Time dilation for the orbiting observer, [itex]dt / d\tau[/itex]
    [tex]\frac{d t}{d \tau} = c \sqrt{\frac{(\frac{d h}{d r}) } { f (\frac{d h}{d r}) - h (\frac{d f}{d r})}}[/tex]

    Coordinate and proper orbital periods are just2 pi over the respective angular velocities

    The orbital period as measured by a co-located static observer is:
    [tex] \left( \frac{2 \pi}{c} \right){ \sqrt{f (\frac{d h}{d r}) / (\frac{d f}{d r}) } }[/tex]

    The set of points with coordinate values r=constant is a circle. It's circumference is:

    [tex]2 \pi \sqrt{h} [/tex]

    The orbital velocity measured by a co-located static observer is:
    [tex]c \sqrt{ \frac{h (\frac{d f}{d r}) }{ f (\frac{d h}{d r})}}[/tex]

    The concept of velocity from "remote points of view" is not well-defined in GR, due to the problems of parallel tranpsorting velocities.

    http://math.ucr.edu/home/baez/einstein/node2.html

    though various ad-hoc approaches (mostly based on specific coordinate systems) are possible.

    Note: The factors of c are annoyig to keep tract of - I took the simple approach of making sure the units come out correctly.
     
Know someone interested in this topic? Share this thread via Reddit, Google+, Twitter, or Facebook




Similar Discussions: Some general formulae for circular orbits in symmetric spacetimes
Loading...