In general relativistic magnetohydrodynamic (GRMHD) simulations, accreted magnetic flux on the black hole horizon episodically decays, during which magnetic reconnection heats up the plasma near the horizon, potentially powering high-energy flares like those observed in M87* and Sgr A*. We study the mm observational counterparts of such flaring episodes in very high resolution GRMHD simulations. The change in 230 GHz flux during the expected high energy flares depends primarily on the efficiency of accelerating γ ≳ 100 (Te ≳ 1011 K) electrons. For models in which the electrons are heated to Te ∼ 1011 K during flares, the hot plasma produced by reconnection significantly enhances 230 GHz emission and increases the size of the 230 GHz image. By contrast, for models in which the electrons are heated to higher temperatures (which we argue are better motivated), the reconnection-heated plasma is too hot to produce significant 230 GHz synchrotron emission, and the 230 GHz flux decreases during high energy flares. We do not find a significant change in the mm polarization during flares as long as the emission is Faraday thin. We also present expectations for the ring-shaped image as observed by the Event Horizon Telescope during flares, as well as multiwavelength synchrotron spectra. Our results highlight several limitations of standard post-processing prescriptions for the electron temperature in GRMHD simulations. We also discuss the implications of our results for current and future observations of flares in Sgr A*, M87*, and related systems. Appendices contain detailed convergence studies with respect to resolution and plasma magnetization.
Sgr A* exhibits flares in the near-infrared and X-ray bands, with the luminosity in these bands increasing by factors of 10–100 for ≈60 min. One of the models proposed to explain these flares is synchrotron emission of non-thermal particles accelerated by magnetic reconnection events in the accretion flow. We use the results from particle-in-cell simulations of magnetic reconnection to post-process 3D two-temperature GRMHD simulations of a magnetically arrested disc (MAD). We identify current sheets, retrieve their properties, estimate their potential to accelerate non-thermal particles, and compute the expected non-thermal synchrotron emission. We find that the flux eruptions of MADs can provide suitable conditions for accelerating non-thermal particles to energies γe ≲ 106 and producing simultaneous X-ray and near-infrared flares. For a suitable choice of current-sheet parameters and a simplified synchrotron cooling prescription, the model can simultaneously reproduce the quiescent and flaring X-ray luminosities as well as the X-ray spectral shape. While the near-infrared flares are mainly due to an increase in the temperature near the black hole during the MAD flux eruptions, the X-ray emission comes from narrow current sheets bordering highly magnetized, low-density regions near the black hole, and equatorial current sheets where the flux on the black hole reconnects. As a result, not all infrared flares are accompanied by X-ray ones. The non-thermal flaring emission can extend to very hard (≲ 100 keV) X-ray energies.
more » « less- Award ID(s):
- 1903335
- NSF-PAR ID:
- 10363238
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 511
- Issue:
- 3
- ISSN:
- 0035-8711
- Page Range / eLocation ID:
- p. 3536-3547
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
ABSTRACT -
more » « less
Abstract Active galactic nuclei in general, and the supermassive black hole in M87 in particular, show bright and rapid gamma-ray flares up to energies of 100 GeV and above. For M87, the flares show multiwavelength components, and the variability timescale is comparable to the dynamical time of the event horizon, suggesting that the emission may come from a compact region near the nucleus. However, the emission mechanism for these flares is not well understood. Recent high-resolution general-relativistic magnetohydrodynamic simulations show the occurrence of episodic magnetic reconnection events that can power flares near the black hole event horizon. In this work, we analyze the radiative properties of the reconnecting current layer under the extreme plasma conditions applicable to the black hole in M87 from first principles. We show that abundant pair production is expected in the vicinity of the reconnection layer, to the extent that the produced secondary pair plasma dominates the reconnection dynamics. Using analytic estimates backed by two-dimensional particle-in-cell simulations, we demonstrate that in the presence of strong synchrotron cooling, reconnection can produce a hard power-law distribution of pair plasma imprinted in the outgoing synchrotron (up to a few tens of MeV) and the inverse-Compton signal (up to TeV). We produce synthetic radiation spectra from our simulations, which can be directly compared with the results of future multiwavelength observations of M87* flares.
-
Abstract We use the public code ebhlight to carry out 3D radiative general relativistic magnetohydrodynamics (GRMHD) simulations of accretion on to the supermassive black hole in M87. The simulations self-consistently evolve a frequency-dependent Monte Carlo description of the radiation field produced by the accretion flow. We explore two limits of accumulated magnetic flux at the black hole (SANE and MAD), each coupled to several subgrid prescriptions for electron heating that are motivated by models of turbulence and magnetic reconnection. We present convergence studies for the radiation field and study its properties. We find that the near-horizon photon energy density is an order of magnitude higher than is predicted by simple isotropic estimates from the observed luminosity. The radially dependent photon momentum distribution is anisotropic and can be modeled by a set of point-sources near the equatorial plane. We draw properties of the radiation and magnetic field from the simulation and feed them into an analytic model of gap acceleration to estimate the very high energy (VHE) γ-ray luminosity from the magnetized jet funnel, assuming that a gap is able to form. We find luminosities of $\rm \sim 10^{41} \, erg \, s^{-1}$ for MAD models and $\rm \sim 2\times 10^{40} \, erg \, s^{-1}$ for SANE models, which are comparable to measurements of M87’s VHE flares. The time-dependence seen in our calculations is insufficient to explain the flaring behaviour. Our results provide a step towards bridging theoretical models of near-horizon properties seen in black hole images with the VHE activity of M87.more » « less
-
more » « less
Abstract It is commonly believed that blazar jets are relativistic magnetized plasma outflows from supermassive black holes. One key question is how the jets dissipate magnetic energy to accelerate particles and drive powerful multiwavelength flares. Relativistic magnetic reconnection has been proposed as the primary plasma physical process in the blazar emission region. Recent numerical simulations have shown strong acceleration of nonthermal particles that may lead to multiwavelength flares. Nevertheless, previous works have not directly evaluated
γ -ray signatures from first-principles simulations. In this paper, we employ combined particle-in-cell and polarized radiation transfer simulations to study multiwavelength radiation and optical polarization signatures under the leptonic scenario from relativistic magnetic reconnection. We find harder-when-brighter trends in optical and Fermi-LATγ -ray bands as well as closely correlated optical andγ -ray flares. The swings in optical polarization angle are also accompanied byγ -ray flares with trivial time delays. Intriguingly, we find highly variable synchrotron self-Compton signatures due to inhomogeneous particle distributions during plasmoid mergers. This feature may result in fastγ -ray flares or orphanγ -ray flares under the leptonic scenario, complementary to the frequently considered minijet scenario. It may also imply neutrino emission with low secondary synchrotron flux under the hadronic scenario, if plasmoid mergers can accelerate protons to very high energy. -
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
