When surrounded by a transparent emission region, black holes are expected to reveal a dark shadow caused by gravitational light bending and photon capture at the event horizon. To image and study this phenomenon, we have assembled the Event Horizon Telescope, a global very long baseline interferometry array observing at a wavelength of 1.3 mm. This allows us to reconstruct event-horizon-scale images of the supermassive black hole candidate in the center of the giant elliptical galaxy M87. We have resolved the central compact radio source as an asymmetric bright emission ring with a diameter of 42 +/- 3 mu as, which is circular and encompasses a central depression in brightness with a flux ratio greater than or similar to 10: 1. The emission ring is recovered using different calibration and imaging schemes, with its diameter and width remaining stable over four different observations carried out in different days. Overall, the observed image is consistent with expectations for the shadow of a Kerr black hole as predicted by general relativity. The asymmetry in brightness in the ring can be explained in terms of relativistic beaming of the emission from a plasma rotating close to the speed of light around a black hole. We compare our images to an extensive library of ray-traced general-relativistic magnetohydrodynamic simulations of black holes and derive a central mass of M = (6.5 +/- 0.7) x 10(9) M-circle dot. Our radio-wave observations thus provide powerful evidence for the presence of supermassive black holes in centers of galaxies and as the central engines of active galactic nuclei. They also present a new tool to explore gravity in its most extreme limit and on a mass scale that was so far not accessible.

/pdf/first-m87-event-horizon-telescope-results-i-the-shadow-of-1ezidlb1ou.pdf

First M87 Event Horizon Telescope Results. I. The Shadow of the Supermassive Black Hole

Black hole accretion flows can be divided into two broad classes: cold and hot. Whereas cold accretion flows consist of cool optically thick gas and are found at relatively high mass accretion rates, hot accretion flows, the topic of this review, are virially hot and optically thin, and occur at lower mass accretion rates. They are described by accretion solutions such as the advection-dominated accretion flow and luminous hot accretion flow. Because of energy advection, the radiative efficiency of these flows is in general lower than that of a standard thin accretion disk. Moreover, the efficiency decreases with decreasing mass accretion rate. Observations show that hot accretion flows are associated with jets. In addition, theoretical arguments suggest that hot flows should produce strong winds. Hot accretion flows are believed to be present in low-luminosity active galactic nuclei and in black hole X-ray binaries in the hard and quiescent states. The prototype is Sgr A*, the ultralow-luminosity supermassive black hole at our Galactic center. The jet, wind, and radiation from a supermassive black hole with a hot accretion flow can interact with the external interstellar medium and modify the evolution of the host galaxy.

Hot Accretion Flows Around Black Holes

We describe new optically thin solutions for rotating accretion flows around black holes and neutron stars. These solutions are advection dominated, so that most of the viscously dissipated energy is advected radially with the flow. We model the accreting gas as a two temperature plasma and include cooling by bremsstrahlung, synchrotron, and Comptonization. We obtain electron temperatures $T_e\sim 10^{8.5}-10^{10}$K. The new solutions are present only for mass accretion rates $\dot M$ less than a critical rate $\dot M_{crit}$ which we calculate as a function of radius $R$ and viscosity parameter $\alpha$. For $\dot M<\dot M_{crit}$ we show that there are three equilibrium branches of solutions. One of the branches corresponds to a cool optically thick flow which is the well-known thin disk solution of Shakura \& Sunyaev (1973). Another branch corresponds to a hot optically thin flow, discovered originally by Shapiro, Lightman \& Eardley (1976, SLE). This solution is thermally unstable. The third branch corresponds to our new advection-dominated solution. This solution is hotter and more optically thin than the SLE solution, but is viscously and thermally stable. It is related to the ion torus model of Rees et al. (1982) and may potentially explain the hard X-ray and $\gamma$-ray emission from X-ray binaries and active galactic nuclei.

/pdf/advection-dominated-accretion-underfed-black-holes-and-cc41v8tdhl.pdf

Advection-Dominated Accretion: Underfed Black Holes and Neutron Stars

Black hole (BH) accretion flows and jets are qualitatively affected by the presence of ordered magnetic fields. We study fully three-dimensional global general relativistic magnetohydrodynamic (MHD) simulations of radially extended and thick (height H to cylindrical radius R ratio of |H/R| {approx} 0.2-1) accretion flows around BHs with various dimensionless spins (a/M, with BH mass M) and with initially toroidally-dominated ({phi}-directed) and poloidally-dominated (R-z directed) magnetic fields. Firstly, for toroidal field models and BHs with high enough |a/M|, coherent large-scale (i.e. >> H) dipolar poloidal magnetic flux patches emerge, thread the BH, and generate transient relativistic jets. Secondly, for poloidal field models, poloidal magnetic flux readily accretes through the disk from large radii and builds-up to a natural saturation point near the BH. While models with |H/R| {approx} 1 and |a/M| {le} 0.5 do not launch jets due to quenching by mass infall, for sufficiently high |a/M| or low |H/R| the polar magnetic field compresses the inflow into a geometrically thin highly non-axisymmetric 'magnetically choked accretion flow' (MCAF) within which the standard linear magneto-rotational instability is suppressed. The condition of a highly-magnetized state over most of the horizon is optimal for the Blandford-Znajek mechanism that generates persistent relativistic jets with and 100% efficiency for |a/M| {approx}> 0.9. A magnetic Rayleigh-Taylor and Kelvin-Helmholtz unstable magnetospheric interface forms between the compressed inflow and bulging jet magnetosphere, which drives a new jet-disk oscillation (JDO) type of quasi-periodic oscillation (QPO) mechanism. The high-frequency QPO has spherical harmonic |m| = 1 mode period of {tau} {approx} 70GM/c{sup 3} for a/M {approx} 0.9 with coherence quality factors Q {approx}> 10. Overall, our models are qualitatively distinct from most prior MHD simulations (typically, |H/R| << 1 and poloidal flux is limited by initial conditions), so they should prove useful for testing accretion-jet theories and measuring a/M in systems such as SgrA* and M87.

General relativistic magnetohydrodynamic simulations of magnetically choked accretion flows around black holes

The Event Horizon Telescope (EHT) has mapped the central compact radio source of the elliptical galaxy M87 at 1.3 mm with unprecedented angular resolution. Here we consider the physical implications of the asymmetric ring seen in the 2017 EHT data. To this end, we construct a large library of models based on general relativistic magnetohydrodynamic (GRMHD) simulations and synthetic images produced by general relativistic ray tracing. We compare the observed visibilities with this library and confirm that the asymmetric ring is consistent with earlier predictions of strong gravitational lensing of synchrotron emission from a hot plasma orbiting near the black hole event horizon. The ring radius and ring asymmetry depend on black hole mass and spin, respectively, and both are therefore expected to be stable when observed in future EHT campaigns. Overall, the observed image is consistent with expectations for the shadow of a spinning Kerr black hole as predicted by general relativity. If the black hole spin and M87's large scale jet are aligned, then the black hole spin vector is pointed away from Earth. Models in our library of non-spinning black holes are inconsistent with the observations as they do not produce sufficiently powerful jets. At the same time, in those models that produce a sufficiently powerful jet, the latter is powered by extraction of black hole spin energy through mechanisms akin to the Blandford-Znajek process. We briefly consider alternatives to a black hole for the central compact object. Analysis of existing EHT polarization data and data taken simultaneously at other wavelengths will soon enable new tests of the GRMHD models, as will future EHT campaigns at 230 and 345 GHz.

/pdf/first-m87-event-horizon-telescope-results-v-physical-origin-2jx2usrxew.pdf

First M87 Event Horizon Telescope Results. V. Physical Origin of the Asymmetric Ring

We present results from two long-duration general relativistic magneto-hydrodynamic (GRMHD) simulations of advection-dominated accretion around a non-spinning black hole. The first simulation was designed to avoid significant accumulation of magnetic flux around the black hole. This simulation was run for a time of 200 000 GM/c3 and achieved inflow equilibrium out to a radius ∼90 GM/c2. Even at this relatively large radius, the mass outflow rate is found to be only 60 per cent of the net mass inflow rate into the black hole. The second simulation was designed to achieve substantial magnetic flux accumulation around the black hole in a magnetically arrested disc. This simulation was run for a shorter time of 100 000 GM/c3. Nevertheless, because the mean radial velocity was several times larger than in the first simulation, it reached inflow equilibrium out to a radius ∼170 GM/c2. Here, becomes equal to at r ∼ 160 GM/c2. Since the mass outflow rates in the two simulations do not show robust convergence with time, it is likely that the true outflow rates are lower than our estimates. The effect of black hole spin on mass outflow remains to be explored. Neither simulation shows strong evidence for convection, though a complete analysis including the effect of magnetic fields is left for the future.

/pdf/grmhd-simulations-of-magnetized-advection-dominated-89ew17x0ww.pdf

GRMHD simulations of magnetized advection‐dominated accretion on a non‐spinning black hole: role of outflows

We present results from two long-duration GRMHD simulations of advection-dominated accretion around a non-spinning black hole. The first simulation was designed to avoid significant accumulation of magnetic flux around the black hole. This simulation was run for a time of 200,000GM/c^3 and achieved inflow equilibrium out to a radius \sim90GM/c^2. Even at this relatively large radius, the mass outflow rate \dot{M}_{out} is found to be only 60% of the net mass inflow rate \dot{M}_{BH} into the black hole. The second simulation was designed to achieve substantial magnetic flux accumulation around the black hole in a magnetically arrested disc. This simulation was run for a shorter time of 100,000GM/c^3. Nevertheless, because the mean radial velocity was several times larger than in the first simulation, it reached inflow equilibrium out to a radius \sim170GM/c^2. Here, \dot{M}_{out} becomes equal to \dot{M}_{BH} at r\sim 160GM/c^2. Since the mass outflow rates in the two simulations do not show robust convergence with time, it is likely that the true outflow rates are lower than our estimates. The effect of black hole spin on mass outflow remains to be explored. Neither simulation shows strong evidence for convection, though a complete analysis including the effect of magnetic fields is left for the future.

GRMHD Simulations of Magnetized Advection Dominated Accretion on a Non-Spinning Black Hole: Role of Outflows

We report on the first global three-dimensional (3D) MHD simulations of disk accretion onto a rotating magnetized star through the Rayleigh-Taylor instability. The star has a dipole field misaligned relative to the rotation axis by a small angle Θ. Simulations show that, depending on the accretion rate, a star may be in the stable or unstable regime of accretion. In the unstable regime, matter penetrates deep into the magnetosphere through several elongated tongues which deposit matter at random places on the surface of the star, leading to stochastic light curves. In the stable regime, matter accretes in ordered funnel streams and the light curves are almost periodic. A star may switch between these two regimes depending on the accretion rate and may thus show alternate episodes of ordered pulsations and stochastic light curves. In the intermediate regime, both stochastic and ordered pulsations are observed. For -->Θ > 30°, the instability is suppressed and stable accretion through funnel streams dominates.

https://iopscience.iop.org/article/10.1086/527298/pdf

Unstable disk accretion onto magnetized stars : first global three-dimensional magnetohydrodynamic simulations

We report on the first global three-dimensional (3D) MHD simulations of disk accretion onto a rotating magnetized star through the Rayleigh-Taylor instability. The star has a dipole field misaligned relative to the rotation axis by a small angle Theta. Simulations show that, depending on the accretion rate, a star may be in the stable or unstable regime of accretion. In the unstable regime, matter penetrates deep into the magnetosphere through several elongated "tongues" which deposit matter at random places on the surface of the star, leading to stochastic light-curves. In the stable regime, matter accretes in ordered funnel streams and the light-curves are almost periodic. A star may switch between these two regimes depending on the accretion rate and may thus show alternate episodes of ordered pulsations and stochastic light-curves. In the intermediate regime, both stochastic and ordered pulsations are observed. For Theta > 30 degrees, the instability is suppressed and stable accretion through funnel streams dominates.

Unstable Disk Accretion to Magnetized Stars: First Global 3D MHD Simulations

We present results of 3D simulations of magnetohydrodynamics (MHD) instabilities at the accretion disc‐magnetosphere boundary. The instability is Rayleigh‐Taylor, and develops for a fairly broad range of accretion rates and stellar rotation rates and magnetic fields. It manifests itself in the form of tall, thin tongues of plasma that penetrate the magnetosphere in the equatorial plane. The shape and number of the tongues changes with time on the inner disc dynamical time-scale. In contrast with funnel flows, which deposit matter mainly in the polar region, the tongues deposit matter much closer to the stellar equator. The instability appears for relatively small misalignment angles, � � 30 ◦ , between the star’s rotation and magnetic axes, and is associated with higher accretion rates. The hotspots and light curves during accretion through instability are generally much more chaotic than during stable accretion. The unstable state of accretion has possible implications for quasi-periodic oscillations and intermittent pulsations from accreting systems, as well as planet migration.

Akshay K. Kulkarni

Papers

GRMHD simulations of magnetized advection‐dominated accretion on a non‐spinning black hole: role of outflows

GRMHD Simulations of Magnetized Advection Dominated Accretion on a Non-Spinning Black Hole: Role of Outflows

Unstable disk accretion onto magnetized stars : first global three-dimensional magnetohydrodynamic simulations

Unstable Disk Accretion to Magnetized Stars: First Global 3D MHD Simulations

Accretion to magnetized stars through the Rayleigh–Taylor instability: global 3D simulations