The brightest point in accretion disk and black hole spin: implication to the image of black hole M87*

TL;DR

This paper introduces a new method to estimate black hole spin from high-resolution images by analyzing the brightest point in the accretion disk relative to the black hole silhouette, applied to M87*.

Contribution

It presents a novel approach to determine black hole spin using the position of the brightest accretion disk point in high-resolution images.

Findings

Estimated M87* black hole spin as a=0.75±0.15

Established a numerical relation between silhouette position and spin

Demonstrated the method's applicability to Event Horizon Telescope data

Abstract

We propose the simple new method for extracting the value of the black hole spin from the direct high-resolution image of black hole by using a thin accretion disk model. In this model the observed dark region on the first image of the supermassive black hole in the galaxy M87, obtained by the Event Horizon Telescope, is a silhouette of the black hole event horizon. The outline of this silhouette is the equator of the event horizon sphere. The dark silhouette of the black hole event horizon is placed within the expected position of the black hole shadow, which is not revealed on the first image. We calculated numerically the relation between the observed position of the black hole silhouette and the brightest point in the thin accretion disk, depending on the black hole spin. From this relation we derive the spin of the supermassive black hole M87*, $a=0.75\pm0.15$.

Equations

$$R(r)=[E(r^2+a^2)-La{]}^2-(r^2-2r+a^2)[{\mu }^2{r}^2+(L-aE{)}^2+Q]$$

$$\Theta (\theta )=Q-{\cos }^2\theta [a^2({\mu }^2-E^2)+L^2/{\sin }^2\theta ].$$

$$\alpha =-\frac{\lambda }{\sin {\theta }_0},\beta =\pm \sqrt{\Theta ({\theta }_0)},$$

$${\text{fint}}^r\frac{dr}{\sqrt{R(r)}}={\text{fint}}^{\theta }\frac{d\theta }{\sqrt{\Theta (\theta )}},$$

$$\phi ={\text{fint}}^r\frac{a(r^2+a^2-\lambda a)}{(r^2-2r+a^2)\sqrt{R(r)}}dr+{\text{fint}}^{\theta }\frac{\lambda -a{\sin }^2\theta }{{\sin }^2\theta \sqrt{\Theta (\theta )}}d\theta ,$$

$$t={\text{fint}}^r\frac{(r^2+a^2)P}{(r^2-2r+a^2)\sqrt{R(r)}}dr+{\text{fint}}^{\theta }\frac{(L-aE{\sin }^2\theta )a}{\sqrt{\Theta (\theta )}}d\theta ,$$

$${\int }_{r_s}^{r_0}\frac{dr}{\sqrt{R(r)}}={\int }_{{\theta }_s}^{{\theta }_0}\frac{d\theta }{\sqrt{\Theta (\theta )}}.$$

$${\int }_{r_s}^{r_0}\frac{dr}{\sqrt{R(r)}}={\int }_{{\theta }_{\min }}^{{\theta }_s}\frac{d\theta }{\sqrt{\Theta (\theta )}}+{\int }_{{\theta }_{\min }}^{{\theta }_0}\frac{d\theta }{\sqrt{\Theta (\theta )}}.$$

$$\lambda =\frac{-r^3+3r^2-a^2(r+1)}{a(r-1)},q^2=\frac{r^3[4a^2-r(r-3{)}^2]}{a^2(r-1{)}^2}$$

$$r_{\mathrm{ISCO}}=3+Z_2-\sqrt{(3-Z_1)(3+Z_1+2Z_2)},$$

$$Z_1=1+(1-a^2{)}^{1/3}[(1+a{)}^{1/3}+(1-a{)}^{1/3}],Z_2=\sqrt{3a^2+Z_1^2}.$$

$$\frac{E}{\mu }$$

$$\frac{L}{\mu }$$

$$\omega =\frac{2ar}{(r^2+a^2{)}^2-a^2(r^2-2r+a^2){\sin }^2\theta }.$$

$$\mathcal{E}(\lambda ,q)=\frac{p^{(t)}-V^{(\phi )}{p}^{(\phi )}-V^{(r)}{p}^{(r)}}{\sqrt{1-[V^{(r)}{]}^2-[V^{(\phi )}{]}^2}}.$$

$$V^{(\phi )}=\frac{r\sqrt{r^2-2r+a^2}L}{[r^3+a^2(r+2)]E-2aL}.$$

$$V^{(r)}=-\sqrt{\frac{r^3+a^2(r+2)}{r}}\frac{\sqrt{R(r)}}{[r^3+a^2(r+2)]E-2aL},$$

$$p^{(\phi )}=\sqrt{\frac{r}{r^3+a^2(r+2)}}\lambda ,p^{(t)}=\sqrt{\frac{r^3+a^2(r+2)}{r(r^2-2r+a^2)}}(1-\omega \lambda ),$$

$$p^{(r)}=-\frac{1}{r}\sqrt{\frac{(r^2+a^2-a\lambda {)}^2}{r^2-2r+a^2}-[(a-\lambda {)}^2+q^2]}.$$

Full text

I Introduction

The Event Horizon Telescope (EHT) recently presented for the first time the impressive direct image of the supermassive black hole in the galaxy M87 EHT1 ; EHT2 ; EHT3 ; EHT4 ; EHT5 ; EHT6 . There are two specific features on this image: the dark silhouette and the bright region.

It is obvious to suppose that the bright region is related to the emission of accreting matter around the supermassive black hole M87*. However, what is the physical origin of the dark silhouette, viewed on the presented image? It depends on the black hole astrophysical environment.

The black hole shadow is viewed by the distant observer in the case of a stationary luminous background, placed beyond the unstable circular photon orbit $r_{\rm ph}$ Bardeen73 ; Luminet79 ; Chandra ; Falcke00 ; ZakhPaoIngrNuc05 ; JohPsaltis10 ; Grenzebach14 ; Grenzebach15 ; Gralla15 ; Strom16 ; Gralla16 ; Strom17 ; Strom18 ; Cunha18a ; Cunha18b ; Huang18 . In the pioneering work of J. M. Bardeen Bardeen73 , the black hole shadow was called the “apparent boundary” of the black hole. The outline of the black hole shadow is shown by the dark red closed curves in our figures.

The quite different case is the black hole image produced by the luminous non-stationary matter plunging into black hole in a region between the unstable circular photon orbit (photon ring) $r_{\rm ph}$ and the black hole event horizon $r_{\rm ph}$. The evident example of the luminous non-stationary matter is the inner region of thin accretion disk between the Innermost Stable Circular Orbit $r_{\rm ISCO}$ (see definitions and details, e.g., in BPT ) and the event horizon $r_{\rm h}$. A non-stationary luminous matter in this region is spiralling down toward the event horizon and produces a dark image of the event horizon silhouette instead of the black hole shadow Luminet79 ; Dexter09 ; Bromley97 ; Fanton97 ; Fukue03 ; Fukue03b ; Tamburini11 ; Ru-SenLu16 ; Luminet19 . In addition, the outline (contour) of this silhouette is the image of the event horizon equator Dokuch19 ; doknaz19 ; doknazsm19 .

We calculate below the form of the event horizon silhouette of the rotating Kerr black hole by modelling the emission of the non-stationary luminous matter spiralling down into black hole in the inner region of thin accretion disk adjoining the event horizon. The trajectories of photons, reaching the distant observer from the inner region of the thin accretion disk, are calculated in the geometric optic approximation. It is also supposed that the accretion disk is optically thin. The event horizon silhouette is smaller than the black hole shadow and is placed on the celestial sphere within the awaited position of black hole shadow.

We identify the dark region on the EHT image with the event horizon silhouette and, respectively, the bright region on the EHT image with the brightest point in the accretion disk. For the case of the supermassive black hole M87*, we find the angular distance of the brightest point in the thin accretion disk from the center of the observed event horizon silhouette depending on the black hole spin parameter $a$. As a result, from this dependence, we find the spin of the supermassive black hole M87*, $a=0.75\pm 0.15$.

II Photon Trajectories in the Kerr Metric

We describe the trajectories of photons (null geodesics) in the Kerr metric by using the standard Boyer–Lindquist coordinate system BoyerLindquist with coordinates $(t,r,\theta,\varphi)$ and with units $G=c=1$. Additionally, we put $mG/c^{2}=1$. With these units, the dimensionless radius of the black hole event horizon is $r_{\rm h}=1+\sqrt{1-a^{2}}$, where the dimensionless spin parameter of the Kerr black hole is $0\leq a\leq 1$.

Trajectories of particles with a rest mass $\mu$ in the Kerr space-time are determined by three constants of motion: a total energy $E$, a component of angular momentum parallel to symmetry axis (azimuth angular momentum) $L$ and the Carter constant $Q$, which is related to the non-equatorial motion of particles Carter68 ; Chandra . Particle trajectories (geodesics) in the radial and latitudinal directions are defined, correspondingly, by a radial effective potential

[TABLE]

and a latitudinal effective potential

[TABLE]

Trajectories of photons (null geodesics with $\mu=0$) in the Kerr space-time are determined only by two dimensionless parameters, $\lambda=L/E$ and $q^{2}=Q/E^{2}$. These parameters are related to the impact parameters on the celestial sphere $\alpha$ and $\beta$ seen by the distant observer placed at a given radius $r_{0}>>r_{\rm h}$ (i.e., practically at infinity), at a given latitude $\theta_{0}$ and at a given azimuth $\varphi_{0}$ (see, e.g., Bardeen73 ; CunnBardeen73 for more details):

[TABLE]

where $\Theta(\theta)$ is from Equation (2).

We use integral equations of motion for photons Carter68 ; Chandra for numerical calculations of the gravitational lensing by the Kerr black hole

[TABLE]

[TABLE]

[TABLE]

where the effective potentials $R(r)$ and $\Theta(\theta)$ are from Equations (1) and (2). The integrals in Equations (4), and (5) are understood to be path integrals along the trajectory.

The path integrals in (4), (5) and (6) are the ordinary ones for photon trajectories without the turning points. For example, the path integral integral Equation (4) can be written in this case through the ordinary integrals as

[TABLE]

In the case of photon trajectories with one turning point $\theta_{\rm min}(\lambda,q)$ (an extremum of latitudinal potential $\Theta(\theta)$), integral Equation (4) can be written through the ordinary integrals as

[TABLE]

In general, a lensed black hole produces an infinite number of images CunnBardeen73 ; Viergutz93 ; RauchBlandf94 ; GralHolzWald19 . With the exception of the very special orientation cases, the most luminous image is a so-called direct or prime image, produced by photons that do not intersect the black hole equatorial plane on the way to a distant observer. In the meantime, the secondary images (named also like the higher order images or light echoes) are produced by photons that intersect the black hole equatorial plane several times. The energy flux from secondary images as a rule is very small in comparison to one from the direct image.

III Black Hole Shadow

A black hole shadow is the gravitational capture cross-section of photons from the stationary luminous background placed at a radial distance from a black hole exceeding the radius of unstable photon circular orbit $r_{\rm ph}$ (see the definition in BPT ).

A black hole shadow in the Kerr metric, projected on the celestial sphere and seen by a distant observer in the equatorial plane of the black hole, is determined from the simultaneous solution of equations $R(r)=0$ and $[rR(r)]^{\prime}=0$, where the effective radial potential $R(r)$ is from Eq. (17). The corresponding solution for the black hole shadow (for a distant observer in the black hole equatorial plane) in the parametric form $(\lambda,q)=(\lambda(r),q(r))$ is

[TABLE]

(see, e.g., Bardeen73 ; Chandra for more details).

Quite a different black hole image is produced in the case of a black hole highlighted by the non-stationary luminous matter plunging into a black hole inside the radius of unstable photon circular orbit $r_{\rm ph}$.

IV Event Horizon Silhouette

The observed dark event horizon silhouette is recovered by gravitational lensing of photons emitted in the innermost part of the accretion disk adjoining to the event horizon. In a geometrically thin accretion disk, placed in the equatorial plane of the black hole, there is an inner boundary for stable circular motion, named the marginally stable radius or the Inner Stable Circular Orbit (ISCO), $r=r_{\rm ISCO}$ (see, e.g., BPT for more details):

[TABLE]

where

[TABLE]

The corresponding values of parameters $E$ and $L$ for particles in the accretion disk co-rotating with the black hole at a circular orbit with a radius $r$ are obtained from the simultaneous solution of equations $R=0$ and $dR/dr=0$, where the effective radial potential $R$ is from Equation (1):

[TABLE]

We use a simple model for describing the non-stationary motion of the small gas elements of accreting matter in the region $r_{\rm h}\leq r\leq r_{\rm ISCO}$ and suppose the pure geodesic motion of the separate small gas element (or compact gas clump) in the accretion flow with the conserved orbital parameters $E$ and $L$ from (12) and (13), corresponding to the radius $r=r_{\rm ISCO}$.

See in Figures 1 and 2 the examples of numerically calculated $2D$ trajectories of small disc elements spiralling down into the black holes in the region $r_{\rm h}\leq r\leq r_{\rm ISCO}$.

See also in Figure 4 two $3D$ photon trajectories starting from the thin accretion disk (light green oval) at $r=1.01r_{\rm h}$ in the equatorial plane of the black hole with $a=0.9982$ and reaching the distant observer very near the outline (contour) of the event horizon silhouette (light black region). In the case of M87*, the viewed dark silhouette is the image of the southern hemisphere of the black hole horizon.

We suppose also that the thin accretion disk is transparent and the energy flux in the rest frame of small gas elements is isotropic and conserved during their spiralling down into the black hole.

See in Figure 3 the parameters of photon trajectories $\lambda$ and $q$, reaching the distant observer at $r_{0}\gg r_{\rm h}$ from the rings in thin accretion disk in the black hole equatorial plane at $\theta_{s}=\pi/2$ and $r_{s}=0.01r_{\rm h}$ and $r_{s}=r_{\rm ISCO}$. The pairs $(\lambda,q)$ are the numerically calculated solutions of integral Eqs. (4) and (8). Namely, the pairs of solutions $(\lambda,q)$ of Equation (4), corresponding to the photon trajectories without the turning points, are shown by the blue colors. Respectively, the pairs of solutions $(\lambda,q)$ of Equation (8), corresponding to the photon trajectories with one turning point in latitudinal direction at $\theta=\theta_{\rm min}(\lambda,q)$, are shown by the red colors. The case of the rotating black hole M87* with spin $a=0.75$ is shown. It is physically reasonable to suggest that the black hole spin axis is aligned with the large-scale jet. In the case of the supermassive black hole M87*, the line of sight makes a $\theta_{0}=163^{\circ}$ angle with the spin axis of the black hole. This value of the inclination angle is derived from the detailed VLBI (Very Long Baseline Interferometry) observations of the famous relativistic jet from the M87* Walker18 . Additionally, these VLBI observations indicate the clockwise flow rotation about the jet axis.

To calculate the energy shift of photons and energy flux from the lensed image of accretion disk, it is necessary to take into account both a red-shift in the black hole gravitational field and Doppler effect. It is convenient to use in these calculations the orthonormal Locally Non-Rotating Frames (LNRF) BPT ; Bardeen70 , for which the observers’ world lines are $r=const$, $\theta=const$, $\varphi=\omega t+const$, where a frame dragging angular velocity

[TABLE]

The requested photon energy in the frame, comoving with the small disk element (or compact gas cloud) at $\theta=\pi/2$, is doknazsm19 ; Viergutz93

[TABLE]

In this equation, the azimuth velocity $V^{(\varphi)}$ of the small disk element with orbital parameters $E$, $L$ and $Q=0$ relative the LNRF, falling in the equatorial plane onto a black hole, is

[TABLE]

A corresponding radial velocity in the equatorial plane relative the LNRF is

[TABLE]

where $R(r)$ is defined in (1) with $Q=0$. The components of photon 4-momentum in the LNRF are

[TABLE]

[TABLE]

Respectively, the energy shift (the ratio of photon frequency at infinity to a corresponding one in the rest frame of the small disk element is $g(\lambda,q)=1/{\cal{E}}(\lambda,q)$. This energy shift is used in numerical calculations of the energy flux from the accretion disk elements measured by the distant observer, following formalism by C. T. Cunnungham and J. M. Bardeen CunnBardeen73 . The corresponding values of the observed local energy flux from the non-stationary region $r_{\rm ISCO}\leq r\leq r_{\rm h}$ in the accretion disk are marked in Figures 5–7 by the different levels of color intensity: a more intensive color corresponds to a more intensive energy flux. In the meantime, the different colors from the dark blue to red in these Figures correspond to the red-shifts of observed photons. Near the brightest points in these Figures, the blue color of observed photons is chosen.

The form of the black hole silhouette does not depend on the observed frequency and is defined only by the properties of the black hole gravitational field. For this reason, we are free to choose any frequency (or color) for emission from the brightest point in the accretion disk. We choose the dark blue color for the brightest point, viewed by the distant observer. The observed frequencies (colors) from the other points in the non-stationary region $r_{\rm ISCO}\leq r\leq r_{\rm h}$ are scaled (normalized) relative to the chosen frequency at the brightest point in accordance with the calculated energy shift of the observed photons from Eq. (15).

V Spin of the Black Hole M87*

Direct images of the inner non-stationary part of thin accretion disk in the region $r_{\rm h}\leq r\leq r_{\rm ISCO}$, adjoining to the black hole event horizon, are presented in Figures 2, 3, 6 and 4 for the case of M87*. In these Figures the black region is the event horizon silhouette. The viewed black silhouette in the case of M87* is the image of the southern hemisphere of the black hole horizon. The outline (contour) of this black silhouette is the equator of the event horizon. The closed dark red curve is a border of the black hole shadow. The dashed red circle is the observed position of the black hole event horizon in the imaginary Euclidian space. The magenta arrow is the black hole rotation axes.

The gravitational red-shift and Doppler effect are taken into account in these images. Local artificial colors of the thin accretion disk images are related with an effective local black-body temperature of the accreting gas elements. We made numerical calculations for 10 different values of the black hole spin in the range $0\leq a\leq 1$ with the step 0.1. All these numerical calculations demonstrate that the observed brightest point in the thin accretion disk is always placed at radius $r=r_{\rm ISCO}$ at the point corresponding to photon trajectory without turning point and with the maximum permissible azimuth angular momentum $\lambda>0$. In Figure 3, the position of the brightest point is marked by a star “{\color[rgb]{.75,0,.25}\star}” at the right graph, corresponding to the largest (positive) azimuth angular momentum $\lambda$ of photon with the direct orbit, starting from $r=r_{\rm ISCO}$ and reaching the distant observer without the turning points.

The derived relation between the position of the brightest point in the thin accretion disk and the observed black hole silhouette depending on the black hole spin is shown in Figure 8. From the comparison of this relation with the image obtained by the EHT, it follows that the best fit for the spin of the supermassive black hole of M87* is $a=0.75\pm 0.15$. In this fit, the value $6\;10^{9}M_{\odot}$ for the mass of the supermassive black hole of M87* Gebhardt11 is used. The $1\sigma$ error in this fit corresponds to the observed width of the asymmetric bright ring on the M87* image EHT1 ; EHT2 ; EHT3 ; EHT4 ; EHT5 ; EHT6 . The derived value of the M87* spin is in a general agreement with the other similar estimations Broderick08 ; Li09 ; Feng17 ; Sobyanin18 ; Nokhrina19 ; Tamburini19 ; Bambi19 ; Nemmen19 ; Davoudiasl19 . See in Figure 6 the superposition of the M87* image and the modelled image of the thin accretion disk in the case of $a=0.75$.

VI Conclusions

We calculated numerically the form of the event horizon silhouette of the rotating Kerr black hole by modelling the emission of non-stationary luminous matter spiralling down into the black hole in the inner region of thin accretion disk at $r_{\rm ISCO}\leq r\leq r_{\rm h}$. It was supposed (i) that the motion of the separate small gas element (or compact gas clump) in the accretion flow is purely geodesic; (ii) that the thin accretion disk is transparent and (iii) that the energy flux in the rest frame of small gas elements is isotropic and conserved during their spiralling down into the black hole. The resulting form of the event horizon silhouette does not depend on the local emission of accretion disk and governed completely by the gravitational field of the black hole. In the case of the supermassive black hole M87*, viewed at the inclination angle $163^{\circ}$, a dark silhouette on the EHT image is the southern hemisphere of the black hole event horizon. The contour of this silhouette is the equator of the event horizon.

The brightest point in accretion disk corresponds to the largest (positive) azimuth angular momentum $\lambda$ of photon with the direct orbit, starting from $r=r_{\rm ISCO}$ and reaching the distant observer without the turning points. By using the first image of the black hole M87*, we find the value of the black hole spin $a=0.75\pm 0.15$.

It is clearly seen in the Figures 5, 6 and 7 that the effective horizontal and vertical diameters of the black hole shadow for the case of M87*, with the black hole rotation axis oriented nearly opposite to the direction toward the distant observer, very weakly depend on the black hole spin. These diameters are approximately 2 times greater than the corresponding ones for the dark silhouette of the black hole event horizon viewed on the first image of the supermassive black hole M87*. This means that the large awaited size of the black hole shadow is not reconciled with the size of the dark black region on the first black image presented by the EHT.

Acknowledgments

We are grateful to E. O. Babichev, V. A. Berezin, Yu. N. Eroshenko and A. L. Smirnov for stimulating discussions. The authors acknowledge support from Russian Foundation for Basic Research grant 18-52-15001 NCNIa.