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

Painleve-Gullstrand coordinates

  1. Jul 16, 2006 #1

    pervect

    User Avatar
    Staff Emeritus
    Science Advisor

    This came up in another thread. I thought I'd make a few notes about them.

    It should be spelled Painleve, I can't edit the typo in the title :-(

    Google finds http://www.physics.umd.edu/grt/taj/776b/hw1soln.pdf, which, being a homework solution, will probably disappear soon.

    The line element is:

    [tex]
    -{d{{T}}}^{2}+ \left( d{{r}}+\sqrt {{\frac {2M}{r}}}d{{T}}
    \right) ^{2}+{r}^{2} \left( {d{{\theta}}}^{2}+ \sin ^2
    \theta d {{\phi}}}^{2} \right)
    [/tex]

    The following coordinate transformation will map PG coordinates into Schwarzschild coordinates:

    [tex]
    t = T -2\,\sqrt {2\,M\,r}+4\,M\,\mathrm{arctanh} \left(
    \sqrt {{\frac {r}{2 \, M}}} \right)
    [/tex]

    [add]
    arctanh(x) is defined only for x<1, while arctanh(x) = 1/2 (ln(1+x)-ln(1-x)), more work needs to be done to deal with the sign issues that arise when making x > 1.


    The above
    The metric is not a function of T, therfore [itex]u_0 = g_{0i} u^i =[/itex] = (-1+2M/r) [itex]dT/d\tau[/itex] + sqrt(2M/r) [itex]dr/d\tau[/itex] = constant, a conserved energy-like quantity of the orbit.

    Here [itex]u^i[/itex] is the 4-velocity [itex](dT/d\tau, dr/d\tau, d\theta/d\tau, d\phi / d\tau)[/itex]

    Similarly, since the metric is not a function of [itex]\phi[/itex]

    [itex]u_3 = g_{3i} u^i[/itex] = r^2 sin^2 [itex]\theta[/itex] [itex]d\phi / d\tau[/itex]= constant

    representing a conserved angular momentum-like quantity of the orbit.

    Generally, the orbit will be taken to be in the equatorial plane, [itex]\theta = \pi/2[/itex] and the above two conserved quantites plus the metric equation will be sufficient to calculate orbital motion.
     
    Last edited: Jul 17, 2006
  2. jcsd
  3. Jul 17, 2006 #2

    pervect

    User Avatar
    Staff Emeritus
    Science Advisor

    OK, it looks like

    [tex]
    t = T -2\,M\sqrt {2}\sqrt {{\frac {r}{M}}}-2\,M\ln \left( \left| \sqrt {{\frac {r}{2\, M}}}-1 \right| \right) +2\,M\ln \left( \sqrt {{\frac {r}{2\,M}}}+1 \right)
    [/tex]

    works for the transformation between Schwarzschild and Painleve-Gullstrand inside & outside the horizon, t being Schwarzschld and T being Painleve.

    Another source for the above info, including the conversion equation between Schw and PG:

    http://arxiv.org/abs/gr-qc/0509030

    less likely to go away and more complete than the homework solution (though it's only periphrial to the paper).
     
    Last edited: Aug 15, 2007
  4. Jul 17, 2006 #3

    pervect

    User Avatar
    Staff Emeritus
    Science Advisor

    Another worthwhile note: the solution for r(tau), T(tau) from a free-fall at infinty:

    let Tdot = [itex]dt/d\tau[/itex], rdot = [itex]dr/d\tau[/itex]

    The energy equation says

    (-1 + 2M/r) Tdot + sqrt(2M/r) rdot = constant

    If Tdot = 1 (zero velocity) at r=infinty, this constant must be -1

    The metric equation says

    [tex]
    \left( \sqrt{\frac{2M}{r}} Tdot + rdot \right)^2 - Tdot^2=-1
    [/tex]

    (just divide both sides of the Lorentz interval by ds^2, and note that [itex]d\tau^2 = -ds^2[/itex] to get the above eq).

    Solving, the energy equation and the metric equation gives us the formula for rdot and Tdot. Taking the minus sign for rdot (for radially infalling particles), we get:

    rdot = -sqrt(2M/r)
    Tdot = 1

    This can be integrated to find
    [tex]
    r = \left( \frac{9M}{2} \tau^2 \right) ^ \frac{1}{3}
    [/tex]

    and T = [itex]\tau[/itex]

    We can further note that [itex]d^2 r / d\tau^2 = -M/r^2[/itex]

    Because T=[itex]\tau[/itex], one can re-write r(tau) as R(T) in the above.
     
  5. Jul 24, 2006 #4

    pervect

    User Avatar
    Staff Emeritus
    Science Advisor

    Solving the previous two equations, we have a quadratic in tdot = dt/dtau
    [tex]
    \frac{r-2 m}{2 m} tdot^2 - \frac{E r}{m} tdot + \frac{E^2 r}{2 m} + 1 = 0
    [/tex]

    This has two solutions. Any horizon crossing solution must satisfy

    [tex]
    \frac{dt}{d\tau} = {\frac {2\,m+{E}^{2}r}{Er+\sqrt {2}\sqrt {m}\sqrt {-r+2\,m+{E}^{2}r}}}
    [/tex]

    The form of the required solution is dicated by the fact that the coefficient of r^2 vanishes at r=2m, therefore we must use the 2c/..... form of the solution of the quadratic, the ...... / 2a solution will have a zero denominator at r=2m

    For non-horizon crossing solutions

    [tex]
    \frac{dt}{dtau} = {\frac {Er+\sqrt {-2\,rm+4\,{m}^{2}+2\,m{E}^{2}r}}{r-2\,m}}
    [/tex]

    is a possible solution as well.

    Sticking with horizon crossing solutions for the moment, we have

    [tex]
    \frac{dt}{dtau} = {\frac {2\,m+{E}^{2}r}{Er+\sqrt {2}\sqrt {m}\sqrt {-r+2\,m+{E}^{2}r}}}
    [/tex]

    and

    [tex]
    \frac{dr}{d tau} = -{\frac { \left( -2\,rm+4\,{m}^{2}+2\,m{E}^{2}r+Er\sqrt {2}\sqrt
    {m}\sqrt {-r+2\,m+{E}^{2}r} \right) }{\sqrt {2 r m}
    \left( Er+\sqrt {2}\sqrt {m}\sqrt {-r+2\,m+{E}^{2}r} \right) }}

    [/tex]

    When dr/dtau = 0, one might need to know d^2 r / dtau^2. Solving the geodesic eq directly when dr/dtau = 0 gives

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

    note that the above DOES assume that dr/dtau = 0.
     
    Last edited: Jul 24, 2006
Know someone interested in this topic? Share this thread via Reddit, Google+, Twitter, or Facebook

Have something to add?