Search

Models of the resolved Event Horizon Telescope (EHT) sources Sgr A* and M87* are constrained by observations at multiple wavelengths, resolutions, polarizations, and time cadences. In this paper, we compare unresolved circular polarization (CP) measurements to a library of models, where each model is characterized by a distribution of CP over time. In the library, we vary the spin of the black hole, the magnetic field strength at the horizon (i.e., both SANE and magnetically arrested disk or MAD models), the observer inclination, a parameter for the maximum ion–electron temperature ratio assuming a thermal plasma, and the direction of the magnetic field dipole moment. We find that Atacama Large Millimeter/submillimeter Array (ALMA) observations of Sgr A* are inconsistent with all edge-on ( i = 90°) models. Restricting attention to the MAD models favored by earlier EHT studies of Sgr A*, we find that only models with magnetic dipole moment pointing away from the observer are consistent with ALMA data. We also note that in 26 of the 27 passing MAD models, the accretion flow rotates clockwise on the sky. We provide a table of the means and standard deviations of the CP distributions for all model parameters, along with their trends.

The 230 GHz lightcurves of Sagittarius A* (Sgr A*) predicted by general relativistic magnetohydrodynamics and general relativistic ray-tracing (GRRT) models by the Event Horizon Telescope Collaboration have higher variability M _Δ _T compared to observations. In this series of papers, we explore the origin of such large brightness variability. In this first paper, we performed large GRRT parameter surveys that span from the optically thin to the optically thick regimes, covering the ion-to-electron temperature ratio under strongly magnetized conditions, R _Low , from 1 to 60. We find that increasing R _Low can lead to either an increase or a reduction in M _Δ _T depending on the other model parameters, making it consistent with the observed variability of Sgr A* in some cases. Our analysis of GRRT image snapshots finds that the major contribution to the large M _Δ _T for the R _Low = 1 models comes from the photon ring. However, secondary contributions from the accretion flow are also visible depending on the spin parameter. Our work demonstrates the importance of the electron temperature used for modeling radiatively inefficient accretion flows and places new constraints on the ion-to-electron temperature ratio. A more in-depth analysis for understanding the dependencies of M _Δ _T on R _Low will be performed in subsequent papers.

Fueling and feedback couple supermassive black holes (SMBHs) to their host galaxies across many orders of magnitude in spatial and temporal scales, making this problem notoriously challenging to simulate. We use a multi-zone computational method based on the general relativistic magnetohydrodynamic (GRMHD) code KHARMA that allows us to span 7 orders of magnitude in spatial scale, to simulate accretion onto a non-spinning SMBH from an external medium with a Bondi radius of R _B ≈ 2 × 10 ^5 GM _• / c ^2 , where M _• is the SMBH mass. For the classic idealized Bondi problem, spherical gas accretion without magnetic fields, our simulation results agree very well with the general relativistic analytic solution. Meanwhile, when the accreting gas is magnetized, the SMBH magnetosphere becomes saturated with a strong magnetic field. The density profile varies as ∼ r ^−1 rather than r ^−3/2 and the accretion rate $\dot{M}$ is consequently suppressed by over 2 orders of magnitude below the Bondi rate ${\dot{M}}_{{\rm{B}}}$ . We find continuous energy feedback from the accretion flow to the external medium at a level of $\sim {10}^{-2}\dot{M}{c}^{2}\sim 5\,\times \,{10}^{-5}{\dot{M}}_{{\rm{B}}}{c}^{2}$ . Energy transport across these widely disparate scales occurs via turbulent convection triggered by magnetic field reconnection near the SMBH. Thus, strong magnetic fields that accumulate on horizon scales transform the flow dynamics far from the SMBH and naturally explain observed extremely low accretion rates compared to the Bondi rate, as well as at least part of the energy feedback.

We present an in-depth analysis of the newly proposed correlation function in visibility space, between the E and B modes of linear polarization, hereafter the EB correlation, for a set of time-averaged general relativistic magnetohydrodynamical simulations compared with the phase map from different semianalytic models and the Event Horizon Telescope (EHT) 2017 data for M87*. We demonstrate that the phase map of time-averaged EB correlation contains novel information that might be linked to black hole (BH) spin, accretion state, and electron temperature. A detailed comparison with a semianalytic approach with different azimuthal expansion modes shows that to recover the morphology of real/imaginary part of the correlation function and its phase, we require higher orders of azimuthal modes. To extract the phase features, we use Zernike polynomial reconstruction developing an empirical metric to break degeneracies between models with different BH spins that are qualitatively similar. We use a set of geometrical ring models with various magnetic and velocity field morphologies, showing that both the image space and visibility-based EB -correlation morphologies in magnetically arrested disk simulations can be explained with simple fluid and magnetic field geometries as used in ring models. Standard and normal evolutions by contrast are harder to model, demonstrating that the simple fluid and magnetic field geometries of ring models are not sufficient to describe them owing to higher Faraday rotation depths. A qualitative comparison with the EHT data demonstrates that some of the features in the phase of EB correlation might be well explained by the current models for BH spins and electron temperatures, while others require larger theoretical surveys.

Interpretation of resolved polarized images of black holes by the Event Horizon Telescope (EHT) requires predictions of the polarized emission observable by an Earth-based instrument for a particular model of the black hole accretion system. Such predictions are generated by general relativistic radiative transfer (GRRT) codes, which integrate the equations of polarized radiative transfer in curved spacetime. A selection of ray-tracing GRRT codes used within the EHT Collaboration is evaluated for accuracy and consistency in producing a selection of test images, demonstrating that the various methods and implementations of radiative transfer calculations are highly consistent. When imaging an analytic accretion model, we find that all codes produce images similar within a pixel-wise normalized mean squared error (NMSE) of 0.012 in the worst case. When imaging a snapshot from a cell-based magnetohydrodynamic simulation, we find all test images to be similar within NMSEs of 0.02, 0.04, 0.04, and 0.12 in Stokes I , Q , U , and V , respectively. We additionally find the values of several image metrics relevant to published EHT results to be in agreement to much better precision than measurement uncertainties.

We investigate general relativistic magnetohydrodynamic simulations to determine the physical origin of the twisty patterns of linear polarization seen in spatially resolved black hole images and explain their morphological dependence on black hole spin. By characterizing the observed emission with a simple analytic ring model, we find that the twisty morphology is determined by the magnetic field structure in the emitting region. Moreover, the dependence of this twisty pattern on spin can be attributed to changes in the magnetic field geometry that occur due to the frame dragging. By studying an analytic ring model, we find that the roles of Doppler boosting and lensing are subdominant. Faraday rotation may cause a systematic shift in the linear polarization pattern, but we find that its impact is subdominant for models with strong magnetic fields and modest ion-to-electron temperature ratios. Models with weaker magnetic fields are much more strongly affected by Faraday rotation and have more complicated emission geometries than can be captured by a ring model. However, these models are currently disfavoured by the recent EHT observations of M87*. Our results suggest that linear polarization maps can provide a probe of the underlying magnetic field structure around a black hole, which may then be usable to indirectly infer black hole spins. The generality of these results should be tested with alternative codes, initial conditions, and plasma physics prescriptions.