skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: General relativistic MHD simulations of non-thermal flaring in Sagittarius A*
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
Award ID(s):
1743747 1911080
PAR ID:
10349390
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ;
Publisher / Repository:
Royal Astronomical Society
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
507
Issue:
4
ISSN:
0035-8711
Page Range / eLocation ID:
5281 to 5302
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT 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
  2. Sagittarius A* (Sgr A*), the supermassive black hole at the centre of the Milky Way, undergoes large-amplitude near-infrared (NIR) flares that can coincide with the continuous rotation of the NIR emission region. One promising explanation for this observed NIR behaviour is a magnetic flux eruption, which occurs in three-dimensional General Relativistic Magneto-Hydrodynamic (3D GRMHD) simulations of magnetically arrested accretion flows. After running two-temperature 3D GRMHD simulations, where the electron temperature is evolved self-consistently along with the gas temperature, it is possible to calculate ray-traced images of the synchotron emission from thermal electrons in the accretion flow. Changes in the gas-dominated (σ = b2/2ρ < 1) regions of the accretion flow during a magnetic flux eruption reproduce the NIR flaring and NIR emission region rotation of Sgr A* with durations consistent with observation. In this paper, we demonstrate that these models also predict that large (1.5x – 2x) size increases of the sub-millimeter (sub-mm) and millimeter (mm) emission region follow most NIR flares by 20–50 min. These size increases occur across a wide parameter space of black hole spin (a = 0.3, 0.5, −0.5, and 0.9375) and initial tilt angle between the accretion flow and black hole spin axes θ0 (θ0 = 0°, 16°, and 30°). We also calculate the sub-mm polarization angle rotation and the shift of the sub-mm spectral index from zero to –0.8 during a prominent NIR flare in our high spin (a = 0.9375) simulation. We show that, during a magnetic flux eruption, a large (∼10rg), magnetically dominated (σ > 1), low-density, and high-temperature ‘bubble’ forms in the accretion flow. The drop in density inside the bubble and additional electron heating in accretion flow between 15rg and 25rg leads to a sub-mm size increase in corresponding images. 
    more » « less
  3. 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. 
    more » « less
  4. null (Ed.)
    ABSTRACT Large-amplitude Sgr A* near-infrared (NIR) flares result from energy injection into electrons near the black hole event horizon. Astrometry data show continuous rotation of the emission region during bright flares, and corresponding rotation of the linear polarization angle. One broad class of physical flare models invokes magnetic reconnection. Here, we show that such a scenario can arise in a general relativistic magnetohydrodynamic simulation of a magnetically arrested disc. Saturation of magnetic flux triggers eruption events, where magnetically dominated plasma is expelled from near the horizon and forms a rotating, spiral structure. Dissipation occurs via reconnection at the interface of the magnetically dominated plasma and surrounding fluid. This dissipation is associated with large increases in NIR emission in models of Sgr A*, with durations and amplitudes consistent with the observed flares. Such events occur at roughly the time-scale to re-accumulate the magnetic flux from the inner accretion disc, ≃10 h for Sgr A*. We study NIR observables from one sample event to show that the emission morphology tracks the boundary of the magnetically dominated region. As the region rotates around the black hole, the NIR centroid and linear polarization angle both undergo continuous rotation, similar to the behaviour seen in Sgr A* flares. 
    more » « less
  5. ABSTRACT High-energy astrophysical systems frequently contain collision-less relativistic plasmas that are heated by turbulent cascades and cooled by emission of radiation. Understanding the nature of this radiative turbulence is a frontier of extreme plasma astrophysics. In this paper, we use particle-in-cell simulations to study the effects of external inverse Compton radiation on turbulence driven in an optically thin, relativistic pair plasma. We focus on the statistical steady state (where injected energy is balanced by radiated energy) and perform a parameter scan spanning from low magnetization to high magnetization (0.04 ≲ σ ≲ 11). We demonstrate that the global particle energy distributions are quasi-thermal in all simulations, with only a modest population of non-thermal energetic particles (extending the tail by a factor of ∼2). This indicates that non-thermal particle acceleration (observed in similar non-radiative simulations) is quenched by strong radiative cooling. The quasi-thermal energy distributions are well fit by analytic models in which stochastic particle acceleration (due to, e.g. second-order Fermi mechanism or gyroresonant interactions) is balanced by the radiation reaction force. Despite the efficient thermalization of the plasma, non-thermal energetic particles do make a conspicuous appearance in the anisotropy of the global momentum distribution as highly variable, intermittent beams (for high magnetization cases). The beamed high-energy particles are spatially coincident with intermittent current sheets, suggesting that localized magnetic reconnection may be a mechanism for kinetic beaming. This beaming phenomenon may explain rapid flares observed in various astrophysical systems (such as blazar jets, the Crab nebula, and Sagittarius A*). 
    more » « less