What is the capture cross-section of a black hole region for ultra-relativistic particles?

Question

What is the capture cross-section of a black hole region for ultra-relativistic particles? I have read that it is

$$\sigma ~=~ \frac{27}{4}\pi R^{2}_{s}$$

for a Schwarzschild BH in the geometric optics limit. Where does the coefficient come from?

Edit - Sources:

1. Absorption and emission spectra of a Schwarzschild black hole.
2. Fermion absorption cross section of a Schwarzschild black hole.

Answer 1

What follows is a very rough adaption of chapter 25 in Gravitation by Misner, Thorne, and Wheeler.

Begin with the Schwarzschild metric with polar angle $\theta$ fixed at $\pi/2$: $$ ds^2 = -\left(1 - \frac{R_\mathrm{S}}{r}\right) \mathrm{d}t^2 + \frac{1}{1-R_\mathrm{S}/r} \mathrm{d}r^2 + r^2 \mathrm{d}\phi^2. $$ For a test particle of rest mass $m$,¹ we know by definition $$ g_{\mu\nu} p^\mu p^\nu + m^2 = 0, $$ where $\vec{p}$ is the 4-momentum of the particle. For an affine parameter $\lambda$ parametrizing the worldline of the particle, these two equations can be combined to give $$ -\left(1 - \frac{R_\mathrm{S}}{r}\right) \left(\frac{\mathrm{d}t}{\mathrm{d}\lambda}\right)^2 + \frac{1}{1-R_\mathrm{S}/r} \left(\frac{\mathrm{d}r}{\mathrm{d}\lambda}\right)^2 + r^2 \left(\frac{\mathrm{d}\phi}{\mathrm{d}\lambda}\right)^2 + m^2 = 0 $$

Now the derivative in the first term is simply the energy $E$, which is related to the conserved energy at infinity $E_\infty$ by $E = E_\infty/(1-R_\mathrm{S}/r)$. Furthermore, the definition of angular momentum is $L = r^2 (\mathrm{d}\phi/\mathrm{d}\lambda)$, which is also conserved. Inserting these definitions gives $$ -\frac{E_\infty^2}{1-R_\mathrm{S}/r} + \frac{1}{1-R_\mathrm{S}/r} \left(\frac{\mathrm{d}r}{\mathrm{d}\lambda}\right)^2 + \frac{L^2}{r^2} + m^2 = 0, $$ or $$ \left(\frac{\mathrm{d}r}{\mathrm{d}\lambda}\right)^2 = E_\infty^2 - \left(1 - \frac{R_\mathrm{S}}{r}\right) \left(\frac{L^2}{r^2} + m^2\right). $$ Using again the definition of $L$, we can write $$ \left(\frac{\mathrm{d}r}{\mathrm{d}\phi}\right)^2 = \frac{r^4}{L^2} \left(E_\infty^2 - \left(1 - \frac{R_\mathrm{S}}{r}\right) \left(\frac{L^2}{r^2} + m^2\right)\right), $$ which we can rewrite to be \begin{align} \left(\frac{\mathrm{d}r}{\mathrm{d}\phi}\right)^2 & = r^4 \frac{E_\infty^2-m^2}{L^2} \left(\frac{E_\infty^2}{E_\infty^2-m^2} - \left(1 - \frac{R_\mathrm{S}}{r}\right) \left(\frac{1}{r^2} \left(\frac{L^2}{E_\infty^2-m^2}\right) + \frac{m^2}{E_\infty^2-m^2}\right)\right) \\ & = \frac{r^4}{b^2} \left(\frac{E_\infty^2}{E_\infty^2-m^2} - \left(1 - \frac{R_\mathrm{S}}{r}\right) \left(\frac{b^2}{r^2} + \frac{m^2}{E_\infty^2-m^2}\right)\right), \end{align} where $$ b = \frac{L}{\sqrt{E_\infty^2-m^2}} $$ is the impact parameter, defined to be the ratio of angular to linear momentum.

At this point, we have a nice general result, but to apply it to photons we take the limit $m \to 0$, which gives $$ \left(\frac{\mathrm{d}r}{\mathrm{d}\phi}\right)^2 = \frac{r^4}{b^2} \left(1 - \frac{b^2}{r^2} \left(1 - \frac{R_\mathrm{S}}{r}\right)\right). $$ The radius of closest approach will be the value $r = r_\text{min}$ for which $\mathrm{d}r/\mathrm{d}\phi$ vanishes: \begin{gather} \frac{r_\text{min}^4}{b^2} \left(1 - \frac{b^2}{r_\text{min}^2} \left(1 - \frac{R_\mathrm{S}}{r_\text{min}}\right)\right) = 0; \\ \frac{b^2}{r_\text{min}^2} \left(1 - \frac{R_\mathrm{S}}{r_\text{min}}\right) = 1; \\ b^2 (r_\text{min} - R_\mathrm{S}) = r_\text{min}^3. \end{gather}

The only thing left is to decide what $r_\text{min}$ is allowed to be if the light is to escape. Since at this point in the trajectory the photon's velocity is entirely in the tangential direction by construction, the question becomes how close in can a photon have a circular orbit? The answer is $r_\text{min} = 3R_\mathrm{S}/2$. Thus $b_\text{max}$, the maximum impact parameter for a photon to be captured, obeys $$ b_\text{max}^2 \left(\frac{3}{2} R_\mathrm{S} - R_\mathrm{S}\right) = \left(\frac{3}{2} R_\mathrm{S}\right)^3, $$ or $$ b_\text{max}^2 = \frac{27}{4} R_\mathrm{S}^2. $$ This is exactly what we sought, since it immediately tells us $$ \sigma = \pi b_\text{max}^2 = \frac{27}{4} \pi R_\mathrm{S}^2. $$

---

¹ We could set $m = 0$ here, but by holding off we get slightly more general intermediate results.

Answer 2

The cross section referred to in the papers is the overall cross section for light shining on the black hole. You might think that the cross section of the black hole is just the area of the event horizon i.e. $\pi r_s^2$. However any light approaching more closely than $r = 3r_s/2$ (the last stable orbit) will spiral into the black hole, so this increases the cross section. And finally, originally parallel light rays travelling from infinity get focused in towards the black hole so even a light ray further away than $r = 3r_s/2$ will be focused in and hit the last stable orbit. The end result is that for parallel light from infinity the black hole appears to have a cross section of $27\pi r_s^2/4$.

I say all this with great confidence, but I have absolutely no idea how to calculate this figure. I'm just quoting from this book. The book provides references is you want to dig in to how the calculation is done, but I would guess the calculations are messy and complicated.

Any particle moving fast enough will see a similar cross section to light.

Answer 3

1. Let us work in units where the speed-of-light $c=1$ is one.

   In Ref. 1 is derived the radial geodesic equation for a particle in the equatorial plane $$ (\frac{dr}{d\lambda})^2+2V(r)~=~E^2, \tag{7.47}$$ with potential $$ 2V(r)~:=~(1-\frac{r_s}{r})((\frac{L}{r})^2+\epsilon). \tag{7.48}$$ Here $\epsilon=0$ for a massless particle and $\epsilon=1$ for a massive particle. The energy $E$ and angular momentum $L$ are constants of motion (which reflect Killing-symmetries of the Schwarzschild metric); $\lambda$ is the affine parameter of the geodesic; and $r_s\equiv\frac{2GM}{c^2}$ is the Schwarzschild-radius. (More precisely, in the massive case $\epsilon=1$, the quantities $E$ and $L$ are specific quantities, i.e. quantities per unit rest mass; and $\lambda$ is proper time.)

2. By differentiating eq. (7.47) wrt. $\lambda$, we find that the condition for a circular orbit $$r(\lambda)~\equiv~ r_{*} \qquad\Rightarrow\qquad \frac{dr}{d\lambda}~\equiv~0$$
   is $$V'(r_{*})~=~0\qquad\Leftrightarrow\qquad \frac{2r_{*}}{r_s}~=~3+\epsilon(\frac{r_{*}}{L})^2.\tag{1}$$

3. Let us next investigate an incoming particle, which has non-constant radial coordinate $\lambda\mapsto r(\lambda)$, and that is precisely on the critical border between being captured and not being captured by the black hole. It would have a radial turning point $\frac{dr}{d\lambda}=0$ precisely at the radius $r=r_{*}$, so that $$2V(r_{*})~=~E^2\qquad\Leftrightarrow\qquad (1-\frac{r_s}{r_{*}})((\frac{L}{r_{*}})^2+\epsilon)~=~E^2.\tag{2} $$

4. The massless case $\epsilon=0$. Eq. (1) yields $$r_{*}~=~\frac{3}{2}r_s.\tag{3}$$ Plugging eq. (3) into eq. (2) then yields the ratio $$ \frac{L}{E}~=~\frac{3}{2}\sqrt{3}r_s.\tag{4}$$ We next use that $L$ and $E$ are constants of motion, so that we can easily identify them at spatial infinity $r=\infty$, where special relativistic formulas apply. The critical impact parameter $b$ is precisely this ratio $$ b~=~\frac{L}{p}~=~\frac{L}{E}~\stackrel{(4)}{=}~\underline{\underline{\frac{3}{2}\sqrt{3}r_s}}, \tag{5}$$ which is OP's sought-for formula.

5. The non-relativistic case $v_{\infty}\ll 1$. The specific energy $E\approx 1$ consists mostly of rest energy. Solving eqs. (1) and (2) then leads to a unique solution $$r_{*}~\approx~ 2r_s~\approx~ L.\tag{6}$$ The critical impact parameter $b$ becomes $$ b~=~\frac{L}{v_{\infty}}~\approx~\underline{\underline{2r_s\frac{c}{v_{\infty}}}}, \tag{7}$$ cf. Ref. 2. The cross section is $\sigma=\pi b^2$.

References:

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

2. V.P. Frolov and I.D. Novikov, Black Hole Physics: Basic Concepts and New Developments, p.48.