Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userRadiation GRMHD simulations of M87: funnel properties and prospects for gap acceleration
75
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We use the public code ebhlight to carry out 3D radiative general relativistic magnetohydrodynamics (GRMHD) simulations of accretion onto 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 sub-grid 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) gamma-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 2times 10^{40} , erg , s^{-1}$ for SANE models, which are comparable to measurements of M87s VHE flares. The time-dependence seen in our calculations is insufficient to explain the flaring behavior. Our results provide a step towards bridging theoretical models of near-horizon properties seen in black hole images with the VHE activity of M87.
rate research
Read More
The candidate supermassive black hole in the Galactic Centre, Sagittarius A* (Sgr A*), is known to be fed by a radiatively inefficient accretion flow (RIAF), inferred by its low accretion rate. Consequently, radiative cooling has in general been overlooked in the study of Sgr A*. However, the radiative properties of the plasma in RIAFs are poorly understood. In this work, using full 3D general-relativistic magneto-hydrodynamical simulations, we study the impact of radiative cooling on the dynamical evolution of the accreting plasma, presenting spectral energy distributions and synthetic sub-millimeter images generated from the accretion flow around Sgr A*. These simulations solve the approximated equations for radiative cooling processes self-consistently, including synchrotron, bremsstrahlung, and inverse Compton processes. We find that radiative cooling plays an increasingly important role in the dynamics of the accretion flow as the accretion rate increases: the mid-plane density grows and the infalling gas is less turbulent as cooling becomes stronger. The changes in the dynamical evolution become important when the accretion rate is larger than $10^{-8},M_{odot}~{rm yr}^{-1}$ ($gtrsim 10^{-7} dot{M}_{rm Edd}$, where $dot{M}_{rm Edd}$ is the Eddington accretion rate). The resulting spectra in the cooled models also differ from those in the non-cooled models: the overall flux, including the peak values at the sub-mm and the far-UV, is slightly lower as a consequence of a decrease in the electron temperature. Our results suggest that radiative cooling should be carefully taken into account in modelling Sgr A* and other low-luminosity active galactic nuclei that have a mass accretion rate of $dot{M} > 10^{-7},dot{M}_{rm Edd}$.
One of the main dissipation processes acting on all scales in relativistic jets is thought to be governed by magnetic reconnection. Such dissipation processes have been studied in idealized environments, such as reconnection layers, which evolve in merging islands and lead to the production of plasmoids, ultimately resulting in efficient particle acceleration. In accretion flows onto black holes, reconnection layers can be developed and destroyed rapidly during the turbulent evolution of the flow. We present a series of two-dimensional general-relativistic magnetohydrodynamic simulations of tori accreting onto rotating black holes focusing our attention on the formation and evolution of current sheets. Initially, the tori are endowed with a poloidal magnetic field having a multi-loop structure along the radial direction and with an alternating polarity. During reconnection processes, plasmoids and plasmoid chains are developed leading to a flaring activity and hence to a variable electromagnetic luminosity. We describe the methods developed to track automatically the plasmoids that are generated and ejected during the simulation, contrasting the behaviour of multi-loop initial data with that encountered in typical simulations of accreting black holes having initial dipolar field composed of one loop only. Finally, we discuss the implications that our results have on the variability to be expected in accreting supermassive black holes.
Using general relativistic magnetohydrodynamic simulations of accreting black holes, we show that a suitable subtraction of the linear polarization per pixel from total intensity images can enhance the photon ring features. We find that the photon ring is typically a factor of $simeq 2$ less polarized than the rest of the image. This is due to a combination of plasma and general relativistic effects, as well as magnetic turbulence. When there are no other persistently depolarized image features, adding the subtracted residuals over time results in a sharp image of the photon ring. We show that the method works well for sample, viable GRMHD models of Sgr A* and M87*, where measurements of the photon ring properties would provide new measurements of black hole mass and spin, and potentially allow for tests of the no-hair theorem of general relativity.
We present axisymmetric two-temperature general relativistic radiation magnetohydrodynamic (GRRMHD) simulations of the inner region of the accretion flow onto the supermassive black hole M87. We address uncertainties from previous modeling efforts through inclusion of models for (1) self-consistent dissipative and Coulomb electron heating (2) radiation transport (3) frequency-dependent synchrotron emission, self-absorption, and Compton scattering. We adopt a distance $D=16.7$ Mpc, an observer angle $theta = 20^{circ}$, and consider black hole masses $M/M_{odot} = (3.3times10^{9}, 6.2times10^{9})$ and spins $a_{star} = (0.5, 0.9375)$ in a four-simulation suite. For each $(M, a_{star})$, we identify the accretion rate that recovers the 230 GHz flux from VLBI measurements. We report on disk thermodynamics at these accretion rates ($dot{M}/dot{M}_{mathrm{Edd}} sim 10^{-5}$). The disk remains geometrically thick; cooling does not lead to a thin disk component. While electron heating is dominated by Coulomb rather than dissipation for $r gtrsim 10 GM/c^2$, the accretion disk remains two-temperature. Radiative cooling of electrons is not negligible, especially for $r lesssim 10 GM/c^2$. The Compton $y$ parameter is of order unity. We then compare derived and observed or inferred spectra, mm images, and jet powers. Simulations with $M/M_{odot} = 3.3times10^{9}$ are in conflict with observations. These simulations produce mm images that are too small, while the low-spin simulation also overproduces X-rays. For $M/M_{odot} = 6.2times10^{9}$, both simulations agree with constraints on radio/IR/X-ray fluxes and mm image sizes. Simulation jet power is a factor $10^2-10^3$ below inferred values, a possible consequence of the modest net magnetic flux in our models.
We present global radiation GRMHD simulations of strongly magnetized accretion onto a spinning, stellar mass black hole at sub-Eddington rates. Using a frequency-dependent Monte Carlo procedure for Compton scattering, we self-consistently evolve a two-temperature description of the ion-electron fluid and its radiation field. For an Eddington ratio $L/L_{rm Edd} gtrsim 10^{-3}$, the emergent spectrum forms an apparent power law shape from thermal Comptonization up to a cutoff at $simeq 100$ keV, characteristic of that seen in the hard spectral states of black hole X-ray binary systems. At these luminosities, the radiative efficiency is high ($approx 24%$) and results in a denser midplane region where magnetic fields are dynamically important. For $L/L_{rm Edd} sim 10^{-2}$, our hot accretion flow appears to undergo thermal runaway and collapse. Our simulations demonstrate that hot accretion flows can be radiatively efficient and provide an estimate of their maximum luminosity.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
