Search for: All records

Abstract The black holes in the Event Horizon Telescope sources Messier 87* and Sagittarius A* (Sgr A*) are embedded in a hot, collisionless plasma that is fully described in kinetic theory yet is usually modeled as an ideal, magnetized fluid. In this Letter, we present results from a new set of weakly collisional fluid simulations in which leading-order kinetic effects are modeled as viscosity and heat conduction. Consistent with earlier, lower-resolution studies, we find that overall flow dynamics remain very similar between ideal and nonideal models. For the first time, we synthesize images and spectra of Sgr A* from weakly collisional models—assuming an isotropic, thermal population of electrons—and find that these remain largely indistinguishable from ideal fluid predictions. However, most weakly collisional models exhibit lower light-curve variability, with all magnetically dominated models showing a small but systematic decrease in variability. 

Abstract Models of highly sub-Eddington accretion onto black holes commonly use a single-fluid model for the collisionless, near-horizon plasma. These models must specify an equation of state. It is common to use an ideal gas withp = (γ − 1)uandγ = 4/3, 13/9, or 5/3, but these produce significantly different outcomes. We discuss the origins of this discrepancy and the assumptions underlying the single-fluid model. The main result of this investigation is that under conditions relevant to low-luminosity black hole accretion the best choice of single-fluid adiabatic index is close to but slightly less than 5/3. Along the way we provide a simple equilibrium model for the relation between the ion-to-electron dissipation ratio and the ion-to-electron temperature ratio, and explore the implications for electron temperature fluctuations in Event Horizon Telescope sources. 

Abstract General relativistic magnetohydrodynamics (GRMHD) simulations are an indispensable tool in studying accretion onto compact objects. The Event Horizon Telescope (EHT) frequently uses libraries of ideal GRMHD simulations to interpret polarimetric, event-horizon-scale observations of supermassive black holes at the centers of galaxies. In this work, we present a library of 10 nonradiative, ideal GRMHD simulations that were utilized by the EHT Collaboration in their analysis of Sagittarius A*. The parameter survey explores both low (SANE) and high (MAD) magnetization states across five black hole spinsa_* = −15/16, −1/2, 0, +1/2, +15/16 where each simulation was run out to 30,000GM/c⁻³. We find the angular momentum and energy flux in SANE simulations closely matches the thin-disk value, with minor deviations in prograde models due to fluid forces. This leads to spin equilibrium arounda_* ∼ 0.94, consistent with previous studies. We study the flow of conserved quantities in our simulations and find mass, angular momentum, and energy transport in SANE accretion flows to be primarily inward and fluid dominated. MAD models produce powerful jets with outflow efficiency >1 fora_* = + 0.94, leading to black hole spin-down in prograde cases. We observe outward directed energy and angular momentum fluxes on the horizon, as expected for the Blandford–Znajek mechanism. MAD accretion flows are sub-Keplerian and exhibit greater variability than their SANE counterpart. They are also hotter than SANE disks withinr ≲ 10GM/c⁻². This study is accompanied by a public release of simulation data athttp://thz.astro.illinois.edu/. 

Abstract Synchrotron emissivities, absorptivities, and Faraday rotation and conversion coefficients are needed in modeling a variety of astrophysical sources, including Event Horizon Telescope (EHT) sources. We develop a method for estimating transfer coefficients that exploits their linear dependence on the electron distribution function, decomposing the distribution function into a sum of parts each of whose emissivity can be calculated easily. We refer to this procedure as stochastic averaging and apply it in two contexts. First, we use it to estimate the emissivity of an isotropicκdistribution function with a high-energy cutoff. The resulting coefficients can be evaluated efficiently enough to be used directly in ray-tracing calculations, and we provide an example calculation. Second, we use stochastic averaging to assess the effect of subgrid turbulence on the volume-averaged emissivity and along the way provide a prescription for a turbulent emissivity. We find that for parameters appropriate to EHT sources turbulence reduces the emissivity slightly. In the infrared, turbulence can dramatically increase the emissivity. 

Abstract 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. 

Abstract The variance and fractional variance on a fixed time window (variously known as “rms percent” or “modulation index”) are commonly used to characterize the variability of astronomical sources. We summarize properties of this statistic for a Gaussian process. 

Abstract We develop a “dual-cone” model for millimeter wavelength emission near a spinning black hole. The model consists of optically thin, luminous cones of emission, centered on the spin axis, which are meant to represent jet walls. The resulting image is dominated by a thin ring. We first consider the effect of the black hole’s spin on the image and show that the dominant effect is to displace the ring perpendicular to the projection of the spin axis on the sky by $2a_{*}\sin i+\mathit{}\left(a_{*}^3\right)$. This effect is lower order ina_*than changes in the shape and size of the photon ring itself but is undetectable without a positional reference. We then show that the centerline of the jet can provide a suitable reference: its location is exactly independent of spin if the observer is outside the cone and nearly independent of spin if the observer is inside the cone. If astrophysical uncertainties can be controlled, then spin displacement is large enough to be detectable by future space very long baseline interferometry missions. Finally, we consider ring substructure in the dual-cone model and show that features in total intensity are not universal and depend on the cone-opening angle. 

Abstract The centers of our Galaxy and the nearby Messier 87 are known to contain supermassive black holes, which support accretion flows that radiate across the electromagnetic spectrum. Although the composition of the accreting gas is unknown, it is likely a mix of ionized hydrogen and helium. We use a simple analytic model and a suite of numerical general relativistic magnetohydrodynamic accretion simulations to study how polarimetric images and spectral energy distributions of the source are influenced by the hydrogen/helium content of the accreting matter. We aim to identify general trends rather than make quantitatively precise predictions, since it is not possible to fully explore the parameter space of accretion models. If the ion-to-electron temperature ratio is fixed, then increasing the helium fraction increases the gas temperature; to match the observational flux density constraints, the number density of electrons and magnetic field strengths must therefore decrease. In our numerical simulations, emission shifts from regions of low to high plasma β —both altering the morphology of the image and decreasing the variability of the light curve—especially in strongly magnetized models with emission close to the midplane. In polarized images, we find that the model gas composition influences the degree to which linear polarization is (de)scrambled and therefore affects estimates for the resolved linear polarization fraction. We also find that the spectra of helium-composition flows peak at higher frequencies and exhibit higher luminosities. We conclude that gas composition may play an important role in predictive models for black hole accretion.