Tuesday, February 12, 2019

general relativity - Expression for distance of closest approach in Schwarzschild Geodesics



The Wikipedia article Two-body problem in General Relativity uses two length-scale variables, $a$ and $b$, to simplify the math. For some information about these, consider these statements from the article:



They are constants of the motion and depend on the initial conditions (position and velocity) of the test particle.



and



Here, b can be interpreted as the distance of closest approach.



They first occur in this equation:


$$ \left( \frac{dr}{d\varphi} \right)^{2} = \frac{r^{4}}{b^{2}} - \left( 1 - \frac{r_{s}}{r} \right) \left( \frac{r^{4}}{a^{2}} + r^{2} \right) $$



Some of the prerequisite math can be found in other articles, and particularly Schwarzschild geodesics, where in the orbits of test particles section they give two equations that ostensibly define $a$ and $b$.


$$ \left( 1 - \frac{r_{s}}{r} \right) \left( \frac{dt}{d\tau} \right) = \frac{a}{b} $$


$$ r^{2} \left( \frac{d\varphi}{d\tau} \right) = a c $$


My first problem is that I don't understand how these define $a$ and $b$. As a previous quote implied, they should be calculable given initial conditions of a test particle, not depend on dynamic derivatives like the above equations imply. Let's return to the Two-body article. It mentions the following in-text, implying the straightforward expression of $a = h/c$.



The effective potential V can be re-written in terms of the length a = h/c



While not defined in that article, other articles on the subject give a consistent definition of $h$. Here, $L$ is angular momentum and $\mu$ is the reduced mass. Since we're only dealing with a test particle with low mass compared to the black hole, the reduced mass is its mass.


$$ h = L/ \mu = L / m$$


If you don't want it, we don't even have to include the particle's mass. Let's denote angular momentum as $L = \vec{r} \times m \vec{v}$, so we can write the following. This also implies an interpretation that $a$ is the specific relative angular momentum in the geometrized unit system.



$$ h = \vec{r} \times \vec{v} $$


$$ a = (\vec{r} \times \vec{v}) / c $$


This is even more complicated than necessary, since most of the work on the Two-body page takes motion to be in a single plane, which would eliminate the vector notation. To complete this question, let me explain that I would just solve for $b$ from the known expression of $a$ using one of the equations I've posted, but I don't have enough information. I can see that it would be necessary to have an expression relating $d \phi$ and $d \tau$, so I've searched through the information of Schwarzschild geodesics for this, I can't find it, and I suspect that's because of the simplifications introduced to do math for a test particle.


In summary, we have an interpretation for $b$ and an explicit equation for $a$. What is an explicit equation for $b$?



Answer



Wikipedia states:



Here, $b$ can be interpreted as the distance of closest approach.



Yes, but that interpretation on Wikipedia is made for massless particles with closest approach much bigger than the Schwarzschild radius $r_s$. In general that interpretation is not true.



On the other hand, OP's above equations (v1) are for massive particle, so let us focus on the massive case in what follows.


As OP writes $$ ac~=~h~=~\frac{L}{m_0} $$ is the specific angular momentum $h$. In general,


$$b~=~c\frac{L}{E}~=~\frac{m_0c}{E}h~=~a\frac{m_0c^2}{E}$$


is the ratio between the (specific) angular momentum and the (specific) energy times the speed of light. The quantities $L$ and $E$ are constants of motion, which in turn reflect (some of) the Killing symmetries of the Schwarzschild metric.


The equation for the closest distance to the black hole can be deduced from the radial geodesic equation


$$\frac{1}{c^2}\left( \frac{dr}{d\tau} \right)^{2} = \left(\frac{E}{m_0c^2}\right)^{2} - \left( 1 - \frac{r_{s}}{r} \right) \left( \left(\frac{a}{r}\right)^2 + 1 \right) $$


for a massive particle in the equatorial plane by putting


$$\frac{dr}{d\tau}~=~0$$


and solve this third-order equation for $r$.


References:




  1. S. Carroll, Lecture Notes on General Relativity, Chapter 7, p.172-179. The pdf file is available from his website.


No comments:

Post a Comment

classical mechanics - Moment of a force about a given axis (Torque) - Scalar or vectorial?

I am studying Statics and saw that: The moment of a force about a given axis (or Torque) is defined by the equation: $M_X = (\vec r \times \...