Find Geodesics in Dynamic Ellis Orbits Metric

In summary: The geodesic Lagrangian method will tell you, probably more easily for this metric than the brute force way of computing the geodesic equation. That's because your metric has only one function of the coordinates (the coefficient of ##d\phi^2## is a function of ##t## and ##p##), so most of the Christoffel symbols vanish. The geodesic Lagrangian method automatically...So it would be more efficient to use the geodesic Lagrangian method?Yes.
  • #1
Onyx
139
4
TL;DR Summary
Orbits in dynamic Ellis Wormhole
Does anyone see a way I can find geodesics in the metric ##ds^2=-dt^2+dp^2+(5p^2+4t^2)d\phi^2## (ones with nonzero angular momentum)? I'm hoping it can be done analytically, but that may be wishful thinking. FYI, this is the metric listed at the bottom of the Wikipedia article about Ellis Wormholes.
 
Physics news on Phys.org
  • #2
Is there a reason why you can't use the techniques already mentioned in the previous thread you had on geodesics?
 
  • #5
PeterDonis said:
Is there a reason why you can't use the techniques already mentioned in the previous thread you had on geodesics?
The techniques mentioned in the FLRW metric paper?
 
  • #6
Onyx said:
The techniques mentioned in the FLRW metric paper?
What techniques are you referring to?
 
  • #8
It may help to note that ##\partial_\phi## is manifestly a Killing field, so should yield a conserved quantity that you can probably use to eliminate derivatives of ##\phi##.
 
  • #9
Ibix said:
It may help to note that ##\partial_\phi## is manifestly a Killing field, so should yield a conserved quantity that you can probably use to eliminate derivatives of ##\phi##.
The geodesic Lagrangian method that I referred to in the previous thread automatically does this kind of elimination.
 
  • Like
Likes Ibix
  • #10
PeterDonis said:
The geodesic Lagrangian method that I referred to in the previous thread automatically does this kind of elimination.
I have not tried this yet. But there are no ##t## or ##p## killing fields, so I don't know what I would do with the ##t## and ##p## even if I did.
 
  • #11
Onyx said:
I have not tried this yet. But there are no ##t## or ##p## killing fields, so I don't know what I would do with the ##t## and ##p## even if I did.
Would doing the Lagrangian method produce equations that are more solvable than ##ds^2=-dt^2+dp^2+\frac{L^2}{5p^2+4t^2}##?
 
  • #12
Onyx said:
there are no ##t## or ##p## killing fields, so I don't know what I would do with the ##t## and ##p## even if I did.
The geodesic Lagrangian method will tell you, probably more easily for this metric than the brute force way of computing the geodesic equation. That's because your metric has only one function of the coordinates (the coefficient of ##d\phi^2## is a function of ##t## and ##p##), so most of the Christoffel symbols vanish. The geodesic Lagrangian method automatically filters out those cases and focuses attention on the small number of cases where the Christoffel symbols don't vanish and the geodesic equation has some meaningful content.
 
  • #13
Onyx said:
Would doing the Lagrangian method produce equations that are more solvable than ##ds^2=-dt^2+dp^2+\frac{L^2}{5p^2+4t^2}##?
I don't know what you mean by "more solvable". The metric isn't an equation you "solve" in the first place. You certainly can't just "solve" the metric if you want to find the geodesics. You have to solve the geodesic equation.

Basically it seems like you are looking for some magical shortcut that will tell you the answer without you actually having to do the math. There isn't one.
 
  • #14
Onyx said:
##ds^2=-dt^2+dp^2+\frac{L^2}{5p^2+4t^2}##?
This is a different metric from the one you gave in the OP. Which one are you using?
 
  • #15
PeterDonis said:
I don't know what you mean by "more solvable". The metric isn't an equation you "solve" in the first place. You certainly can't just "solve" the metric if you want to find the geodesics. You have to solve the geodesic equation.

Basically it seems like you are looking for some magical shortcut that will tell you the answer without you actually having to do the math. There isn't one
I should have said that I'm wondering if there is a way to integrate geodesics exactly instead of using Euler or Runge-Kutta.
PeterDonis said:
This is a different metric from the one you gave in the OP. Which one are you using?
This is the metric:
##ds^2=-dt^2+dp^2+(5p^2+4t^2)d\phi^2##
PeterDonis said:
.
 
  • #16
Onyx said:
I should have said that I'm wondering if there is a way to integrate geodesics exactly instead of using Euler or Runge-Kutta.
Some differential equations can be integrated exactly and some can't. The only way to find out which kind you're dealing with is to look at the equations.

Onyx said:
This is the metric:
##ds^2=-dt^2+dp^2+(5p^2+4t^2)d\phi^2##
Ok.
 
  • #17
PeterDonis said:
I don't know what you mean by "more solvable". The metric isn't an equation you "solve" in the first place. You certainly can't just "solve" the metric if you want to find the geodesics. You have to solve the geodesic equation.

Basically it seems like you are looking for some magical shortcut that will tell you the answer without you actually having to do the math. There isn't one
I should have said that I'm wondering if
PeterDonis said:
The geodesic Lagrangian method will tell you, probably more easily for this metric than the brute force way of computing the geodesic equation. That's because your metric has only one function of the coordinates (the coefficient of ##d\phi^2## is a function of ##t## and ##p##), so most of the Christoffel symbols vanish. The geodesic Lagrangian method automatically filters out those cases and focuses attention on the small number of cases where the Christoffel symbols don't vanish and the geodesic equation has some meaningful content.
Is the Lagrangian method only for timelike geodesics?
PeterDonis said:
.
 
  • #18
Onyx said:
Is the Lagrangian method only for timelike geodesics?
The Wikipedia article describes it that way in the section I referenced, but you can actually use it for any kind of geodesic (with some technical complications required to justify the method for null geodesics, not discussed in the article). The key point is that the quantity whose variation is set to zero is the same regardless of the kind of geodesic.
 
  • #19
PeterDonis said:
The Wikipedia article describes it that way in the section I referenced, but you can actually use it for any kind of geodesic (with some technical complications required to justify the method for null geodesics, not discussed in the article). The key point is that the quantity whose variation is set to zero is the same regardless of the kind of geodesic.
And you would recommend this method over just looking at ##dt^2=dp^2+(5p^2+4t^2)d\phi^2## in the null case?
 
  • #20
Onyx said:
And you would recommend this method over just looking at ##dt^2=dp^2+(5p^2+4t^2)d\phi^2## in the null case?
What would just looking at that equation do for you?
 
  • #21
@Onyx, it seems like you're spending an awful lot of time asking questions instead of actually trying anything. Why don't you just try some things? The metric you're talking about is simple enough that the difference in computational effort between different methods isn't that much, and none of them will take a huge amount of effort. You could just try them all and see what works best.
 
  • #22
PeterDonis said:
The geodesic Lagrangian method that I referred to in the previous thread automatically does this kind of elimination.
It is, indeed, quite neat. I just tried it for the Schwarzschild metric; you set ##\theta=\pi/2## and it immediately give you that ##g_{\phi\phi}\frac{d\phi}{d\tau}=\mathrm{const}## and ##g_{tt}\frac{dt}{d\tau}=\mathrm{const}##. The ##r## equation is messy but you can decline to solve it and just plug the above results into ##g_{ab}\dot x^a\dot x^b## to get ##dr/d\tau## in terms of the effective potential.

I guess it's messier for a spacetime that isn't so symmetric as Schwarzschild, or with a poorer choice of coordinates, but it's still nice.
 
Last edited:
  • #23
Ibix said:
It is, indeed, quite neat. I just tried it for the Schwarzschild metric; you set ##\theta=\pi/2## and it immediately give you that ##g_{\phi\phi}\frac{d\phi}{d\tau}=\mathrm{const}## and ##g_{tt}\frac{dt}{d\tau}=\mathrm{const}##. The ##r## equation is messy but you can decline to solve it and just plug the above results into ##g_{ab}\dot x^a\dot x^b## to get ##dr/d\tau## in terms of the effective potential.

I guess it's messier for a spacetime that isn't so symmetric as Schwarzschild, or with a poorer choice of coordinates, but it's still nice.
Excellent. Can you post a walkthrough of how you did it with Schwarzschild?
 
  • #24
Onyx said:
Excellent. Can you post a walkthrough of how you did it with Schwarzschild?
Do you know what the Lagrangian is? (Hint: see post #7.) Do you know the Euler-Lagrange equation? (Hint: google is your friend.)
 
Last edited:
  • #25
Ibix said:
It may help to note that ##\partial_\phi## is manifestly a Killing field, so should yield a conserved quantity that you can probably use to eliminate derivatives of ##\phi##.
For the non-dynamic Ellis wormhole (i.e., no ##t## dependence in the metric), this works easily. For the dynamic case (the one in the OP, with ##t## dependence in the metric), it's more complicated, because the "conserved" quantity you get is explicitly time-dependent.
 
  • #26
PeterDonis said:
the "conserved" quantity you get is explicitly time-dependent.
Presumably because the circumference of a circle at fixed ##p## (which I think is meant to be ##\rho##, by the way) is varying with time and that affects the conservation of angular momentum?
 
  • #27
Ibix said:
Presumably because the circumference of a circle at fixed ##p## (which I think is meant to be ##\rho##, by the way) is varying with time and that affects the conservation of angular momentum?
I think that would be a reasonable physical interpretation, yes.
 
  • Like
Likes Ibix
  • #28
Ibix said:
Do you know what the Lagrangian is? (Hint: see post #7.) Do you know the Euler-Lagrange equation? (Hint: google is your friend.)
Do you mean the method of defining the Lagrangian and then differentiate it w.r.t the one-forms (##dt, d\rho, d\phi##)?
 
  • #30
Onyx said:
Do you mean the method of defining the Lagrangian and then differentiate it w.r.t the one-forms (##dt, d\rho, d\phi##)?
You don't differentiate with respect to the one-forms. You differentiate with respect to the coordinates and their derivatives with respect to the affine parameter.

Ibix said:
Yes.
I don't think so. See above.
 
  • Like
Likes vanhees71 and Ibix
  • #31
Their derivatives with respect to the affine parameter.

This is what I meant by one-forms
 
  • #32
Onyx said:
This is what I meant by one-forms
Then you should be aware that that is not correct terminology.
 
  • Like
Likes vanhees71
  • #33
So, @Onyx, what is the relevant Lagrangian for GR, and for Schwarzschild spacetime specifically?
 
  • Like
Likes vanhees71
  • #34
Ibix said:
So, @Onyx, what is the relevant Lagrangian for GR, and for Schwarzschild spacetime specifically?
Ibix said:
So, @Onyx, what is the relevant Lagrangian for GR, and for Schwarzschild spacetime specifically?
##L=\sqrt{{g_{uv}}u'v'}##
 
  • #35
Onyx said:
##L=\sqrt{{g_{uv}}u'v'}##
Well, I wouldn't bother with the square root, but that's up to you. Set ##\theta=\pi/2## for simplicity's sake, which makes ##\partial\mathcal{L}/\partial\theta=\partial\mathcal{L}/\partial\dot\theta=0##. What do the three remaining Euler-Lagrange equations give you?
 
Back
Top