Horizon-scale observations of the jetted active galactic nucleus M87 are compared with simulations spanning a broad range of dissipation mechanisms and plasma content in three-dimensional general relativistic flows around spinning black holes. Observations of synchrotron radiation from radio to X-ray frequencies can be compared with simulations by adding prescriptions specifying the relativistic electron-plus-positron distribution function and associated radiative transfer coefficients. A suite of time-varying simulations with various spins, plasma magnetizations and turbulent heating and equipartition-based emission prescriptions (and piecewise combinations thereof) is chosen to represent distinct possibilities for the M87 jet/accretion flow/black hole system. Simulation jet morphology, polarization, and variation are then ‘observed’ and compared with real observations to infer the rules that govern the polarized emissivity. Our models support several possible spin/emission model/plasma composition combinations supplying the jet in M87, whose black hole shadow has been observed down to the photon ring at 230 GHz by the Event Horizon Telescope (EHT). Net linear polarization and circular polarization constraints favour magnetically arrested disc (MAD) models whereas resolved linear polarization favours standard and normal evolution (SANE) in our parameter space. We also show that some MAD cases dominated by intrinsic circular polarization have near-linear V/I dependence on un-paired electron or positron content while SANE polarization exhibits markedly greater positron-dependent Faraday effects – future probes of the SANE/MAD dichotomy and plasma content with the EHT. This is the second work in a series also applying the ‘observing’ simulations methodology to near-horizon regions of supermassive black holes in Sgr A* and 3C 279.
First Sagittarius A* Event Horizon Telescope Results. V. Testing Astrophysical Models of the Galactic Center Black Hole
Abstract In this paper we provide a first physical interpretation for the Event Horizon Telescope's (EHT) 2017 observations of Sgr A*. Our main approach is to compare resolved EHT data at 230 GHz and unresolved non-EHT observations from radio to X-ray wavelengths to predictions from a library of models based on time-dependent general relativistic magnetohydrodynamics simulations, including aligned, tilted, and stellar-wind-fed simulations; radiative transfer is performed assuming both thermal and nonthermal electron distribution functions. We test the models against 11 constraints drawn from EHT 230 GHz data and observations at 86 GHz, 2.2 μ m, and in the X-ray. All models fail at least one constraint. Light-curve variability provides a particularly severe constraint, failing nearly all strongly magnetized (magnetically arrested disk (MAD)) models and a large fraction of weakly magnetized models. A number of models fail only the variability constraints. We identify a promising cluster of these models, which are MAD and have inclination i ≤ 30°. They have accretion rate (5.2–9.5) × 10 −9 M ⊙ yr −1 , bolometric luminosity (6.8–9.2) × 10 35 erg s −1 , and outflow power (1.3–4.8) × 10 38 erg s −1 . We also find that all models with i ≥ 70° fail at least two constraints, as do all models with equal ion and electron temperature; exploratory, nonthermal model sets tend to have higher 2.2 μ m flux density; and the population of cold electrons is limited by X-ray constraints due to the risk of bremsstrahlung overproduction. Finally, we discuss physical and numerical limitations of the models, highlighting the possible importance of kinetic effects and duration of the simulations.
more »
« less
- NSF-PAR ID:
- 10330547
- Author(s) / Creator(s):
- ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; more »
- Date Published:
- Journal Name:
- The Astrophysical Journal Letters
- Volume:
- 930
- Issue:
- 2
- ISSN:
- 2041-8205
- Page Range / eLocation ID:
- L16
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
ABSTRACT -
Abstract We report a timing analysis of near-infrared (NIR), X-ray, and submillimeter data during a 3 day coordinated campaign observing Sagittarius A*. Data were collected at 4.5 μ m with the Spitzer Space Telescope, 2–8 keV with the Chandra X-ray Observatory, 3–70 keV with NuSTAR, 340 GHz with ALMA, and 2.2 μ m with the GRAVITY instrument on the Very Large Telescope Interferometer. Two dates show moderate variability with no significant lags between the submillimeter and the infrared at 99% confidence. A moderately bright NIR flare ( F K ∼ 15 mJy) was captured on July 18 simultaneous with an X-ray flare ( F 2−10 keV ∼ 0.1 counts s −1 ) that most likely preceded bright submillimeter flux ( F 340 GHz ∼ 5.5 Jy) by about + 34 − 33 + 14 minutes at 99% confidence. The uncertainty in this lag is dominated by the fact that we did not observe the peak of the submillimeter emission. A synchrotron source cooled through adiabatic expansion can describe a rise in the submillimeter once the synchrotron self-Compton NIR and X-ray peaks have faded. This model predicts high GHz and THz fluxes at the time of the NIR/X-ray peak and electron densities well above those implied from average accretion rates for Sgr A*. However, the higher electron density postulated in this scenario would be in agreement with the idea that 2019 was an extraordinary epoch with a heightened accretion rate. Since the NIR and X-ray peaks can also be fit by a nonthermal synchrotron source with lower electron densities, we cannot rule out an unrelated chance coincidence of this bright submillimeter flare with the NIR/X-ray emission.more » « less
-
ABSTRACT In large-scale hydrodynamical cosmological simulations, the fate of massive galaxies is mainly dictated by the modelling of feedback from active galactic nuclei (AGNs). The amount of energy released by AGN feedback is proportional to the mass that has been accreted on to the black holes (BHs), but the exact subgrid modelling of AGN feedback differs in all simulations. While modern simulations reliably produce populations of quiescent massive galaxies at z ≤ 2, it is also crucial to assess the similarities and differences of the responsible AGN populations. Here, we compare the AGN populations of the Illustris, TNG100, TNG300, Horizon-AGN, EAGLE, and SIMBA simulations. The AGN luminosity function (LF) varies significantly between simulations. Although in agreement with current observational constraints at z = 0, at higher redshift the agreement of the LFs deteriorates with most simulations producing too many AGNs of $L_{\rm x, 2\!-\!10 \, keV}\sim 10^{43\!-\!44}\, \rm erg\, s^{-1}$. AGN feedback in some simulations prevents the existence of any bright AGN with $L_{\rm x, 2\!-\!10 \, keV}\geqslant 10^{45}\rm \,erg\, s^{-1}$ (although this is sensitive to AGN variability), and leads to smaller fractions of AGN in massive galaxies than in the observations at z ≤ 2. We find that all the simulations fail at producing a number density of AGN in good agreement with observational constraints for both luminous ($L_{\rm x, 2\!-\!10 \, keV}\sim 10^\text{43-45}\, \rm erg\, s^{-1}$) and fainter ($L_{\rm x, 2\!-\!10 \, keV}\sim 10^\text{42-43}\, \rm erg\, s^{-1}$) AGNs and at both low and high redshifts. These differences can aid us in improving future BH and galaxy subgrid modelling in simulations. Upcoming X-ray missions (e.g. Athena, AXIS, and LynX) will bring faint AGNs to light and new powerful constraints. After accounting for AGN obscuration, we find that the predicted number density of detectable AGNs in future surveys spans at least one order of magnitude across the simulations, at any redshift.more » « less
-
ABSTRACT Sgr A* exhibits regular variability in its multiwavelength emission, including daily X-ray flares and roughly continuous near-infrared (NIR) flickering. The origin of this variability is still ambiguous since both inverse Compton and synchrotron emission are possible radiative mechanisms. The underlying particle distributions are also not well constrained, particularly the non-thermal contribution. In this work, we employ the GPU-accelerated general relativistic magnetohydrodynamics code H-AMR to perform a study of flare flux distributions, including the effect of particle acceleration for the first time in high-resolution 3D simulations of Sgr A*. For the particle acceleration, we use the general relativistic ray-tracing code bhoss to perform the radiative transfer, assuming a hybrid thermal+non-thermal electron energy distribution. We extract ∼60 h light curves in the sub-millimetre, NIR and X-ray wavebands, and compare the power spectra and the cumulative flux distributions of the light curves to statistical descriptions for Sgr A* flares. Our results indicate that non-thermal populations of electrons arising from turbulence-driven reconnection in weakly magnetized accretion flows lead to moderate NIR and X-ray flares and reasonably describe the X-ray flux distribution while fulfilling multiwavelength flux constraints. These models exhibit high rms per cent amplitudes, $\gtrsim 150{{\ \rm per\ cent}}$ both in the NIR and the X-rays, with changes in the accretion rate driving the 230 GHz flux variability, in agreement with Sgr A* observations.more » « less
-
more » « less
ABSTRACT Wind-fed models offer a unique way to form predictive models of the accretion flow surrounding Sagittarius A*. We present 3D wind-fed magnetohydrodynamic (MHD) and general relativistic magnetohydrodynamic (GRMHD) simulations spanning the entire dynamic range of accretion from parsec scales to the event horizon. We expand on previous work by including non-zero black hole spin and dynamically evolved electron thermodynamics. Initial conditions for these simulations are generated from simulations of the observed Wolf–Rayet stellar winds in the Galactic Centre. The resulting flow tends to be highly magnetized (β ≈ 2) with an ∼r−1 density profile independent of the strength of magnetic fields in the winds. Our simulations reach the magnetically arrested disc (MAD) state for some, but not all cases. In tilted flows, standard and normal evolution (SANE) jets tend to align with the angular momentum of the gas at large scales, even if that direction is perpendicular to the black hole spin axis. Conversely, MAD jets tend to align with the black hole spin axis. The gas angular momentum shows similar behaviour: SANE flows tend to only partially align while MAD flows tend to fully align. With a limited number of dynamical free parameters, our models can produce accretion rates, 230 GHz flux, and unresolved linear polarization fractions roughly consistent with observations for several choices of electron heating fraction. Absent another source of large-scale magnetic field, winds with a higher degree of magnetization (e.g. where the magnetic pressure is 1/100 of the ram pressure in the winds) may be required to get a sufficiently large rotation measure with consistent sign.
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.3847/2041-8213/ac6672
