-1
$\begingroup$

I am currently working on some work with neutron stars so this led me to calculations of photon trajectories in curved spacetime.

enter image description here

These are the steps I followed to write code that would produce cartesian coordinates for me to graph.

enter image description here

I utilized this equation to form the steps listed above.

enter image description here

This is the graph of various trajectories when the photon is initialized at phi=0. In these scenarios the only variable adjusted is initial direction of photon motion(alpha). When alpha=pi/2 I get this circular trajectory. When alpha=pi/3, pi/4, and so on, I get the rotating spiral path which appears to be correct. Why is it that the graph is circular? Is this incorrect physically? It is important to note that I used polar coordinates in my calculations then converted back to cartesian.

$\endgroup$
4
  • 4
    $\begingroup$ There's a circular orbit at $r=1.5r_s$ (the photon sphere), but it's unstable, and when using numerical integration you need to use extremely tiny $\Delta\varphi$ steps to get accurate trajectories near the photon sphere. See nature.com/articles/s41598-021-93595-w $\endgroup$
    – PM 2Ring
    Commented Oct 6, 2023 at 6:04
  • $\begingroup$ I have some answers about photon trajectories in the Schwarzschild metric that you might find helpful: physics.stackexchange.com/a/766180/123208 & physics.stackexchange.com/a/680961/123208 $\endgroup$
    – PM 2Ring
    Commented Oct 6, 2023 at 6:05
  • 1
    $\begingroup$ BTW, you should explain what your $\alpha_0$ parameter is. And for goodness sake, set your aspect ratio to 1 so that your circles look like circles. ;) $\endgroup$
    – PM 2Ring
    Commented Oct 6, 2023 at 6:12
  • $\begingroup$ possible follow-up question: what spacetime has a stable circular orbit for photons, and is it a spacetime that could come about? $\endgroup$ Commented Jan 18 at 22:07

1 Answer 1

1
$\begingroup$

I think your result is correct: for $ \phi_0=0 $, if $ \alpha_0={\pi\over 2} $ the path is circular.

This is due to the fact that the equation $ u''=3u^2-u $ can be integrated if you multiply it by $ 2u' $, that is

$ 2u'u''=6u'u^2-2u'u $ which gives the integration

$ u'^2=2u^3-u^2+cte\ \ \ \ [1] $

If $ \phi_0=0 $ and $ \alpha_0={\pi\over 2} $, you have $ u'_0=0 $ (because $ \tan\phi_0=0 $ and $ \tan\alpha_0=\infty $).

This means $ u_1=u_0 $ and $ [1] $ leads to $ u'_1=u'_0 $ which is $ 0 $.

So, you will see that $ \forall n\ u'_n=0 $ and $ u_n=u_0 $ which results in a circular trajectory.

Hoping to have answered your question,

Best regards.

$\endgroup$

Not the answer you're looking for? Browse other questions tagged or ask your own question.