Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userPlasmoid formation in global GRMHD simulations and AGN flares
59
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
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.
rate research
Read More
The Event Horizon Telescope (EHT) collaboration, an Earth-size sub-millimetre radio interferometer, recently captured the first images of the central supermassive black hole in M87. These images were interpreted as gravitationally-lensed synchrotron emission from hot plasma orbiting around the black hole. In the accretion flows around low-luminosity active galactic nuclei such as M87, electrons and ions are not in thermal equilibrium. Therefore, the electron temperature, which is important for the thermal synchrotron radiation at EHT frequencies of 230 GHz, is not independently determined. In this work, we investigate the commonly used parameterised ion-to-electron temperature ratio prescription, the so-called R-$beta$ model, considering images at 230 GHz by comparing with electron-heating prescriptions obtained from general-relativistic magnetohydrodynamical (GRMHD) simulations of magnetised accretion flows in a Magnetically Arrested Disc (MAD) regime with different recipes for the electron thermodynamics. When comparing images at 230 GHz, we find a very good match between images produced with the R-$beta$ prescription and those produced with the turbulent- and magnetic reconnection- heating prescriptions. Indeed, this match is on average even better than that obtained when comparing the set of images built with the R-$beta$ prescription with either a randomly chosen image or with a time-averaged one. From this comparative study of different physical aspects, which include the image, visibilities, broadband spectra, and light curves, we conclude that, within the context of images at 230 GHz relative to MAD accretion flows around supermassive black holes, the commonly-used and simple R-$beta$ model is able to reproduce well the various and more complex electron-heating prescriptions considered here.
We present results of a set of three-dimensional, general relativistic radiation magnetohydrodynamics simulations of thin accretion discs around a non-rotating black hole to test their thermal stability. We consider two cases, one that is initially radiation pressure dominated and expected to be thermally unstable and another that is initially gas-pressure dominated and expected to remain stable. Indeed, we find that cooling dominates over heating in the radiation pressure dominated model, causing the disc to collapse vertically on roughly the local cooling timescale. We also find that heating and cooling within the disc have a different dependence on the mid-plane pressure, a prerequisite of thermal instability. Comparison of our data with the relevant thin-disc thermal equilibrium curve suggests that our disc may be headed for the thermally stable, gas-pressure-dominated branch. However, because the disc collapses to the point that we are no longer able to resolve it, we had to terminate the simulation. On the other hand, the gas pressure dominated model, which was run for twice as long as the radiation pressure dominated one, remains stable, with heating and cooling roughly in balance. Finally, the radiation pressure dominated simulation shows some evidence of viscous instability. The strongest evidence is in plots of surface density, which show the disc breaking up into rings.
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.
We present the results of 2D particle-in-cell (PIC) simulations of relativistic magnetic reconnection (RMR) in electron-positron plasma, including the dynamical influence of the synchrotron radiation process, and integrating the observable emission signatures. The simulations are initiated with a single Harris current layer with a central gap that triggers the RMR process. We achieve a steady-state reconnection with unrestricted outflows by means of open boundary conditions. The radiative cooling efficiency is regulated by the choice of initial plasma temperature Theta. We explore different values of Theta and of the background magnetisation sigma_0. Throughout the simulations, plasmoids are generated in the central region of the layer, and they evolve at different rates, achieving a wide range of sizes. The gaps between plasmoids are filled by smooth relativistic outflows called minijets, whose contribution to the observed radiation is very limited due to their low particle densities. Small-sized plasmoids are rapidly accelerated, however, they have lower contributions to the observed emission, despite stronger relativistic beaming. Large-sized plasmoids are slow, but produce most of the observed synchrotron emission, with major part of their radiation produced within the central cores, the density of which is enhanced by radiative cooling. Synchrotron lightcurves show rapid bright flares that can be identified as originating from mergers between small/fast plasmoids and large/slow targets moving in the same direction. In the high-magnetisation case, the accelerated particles form a broken power-law energy distribution with a soft tail produced by particles accelerated in the minijets.
Luminous active galactic nuclei (AGN) and X-Ray binaries (XRBs) tend to be surrounded by geometrically thin, radiatively cooled accretion discs. According to both theory and observations, these are -- in many cases -- highly misaligned with the black hole spin axis. In this work we present the first general relativistic magnetohydrodynamic simulations of very thin ($h/r sim 0.015-0.05$) accretion discs around rapidly spinning ($a sim 0.9$) black holes and tilted by 45-65 degrees. We show that the inner regions of the discs with $h/r lesssim 0.03$ align with the black hole equator, though at smaller radii than predicted by theoretical work. The inner aligned and outer misaligned disc regions are separated by a sharp break in tilt angle accompanied by a sharp drop in density. We find that frame-dragging by the spinning black hole overpowers the disc viscosity, which is self-consistently produced by magnetized turbulence, tearing the disc apart and forming a rapidly precessing inner sub-disc surrounded by a slowly precessing outer sub-disc. We find that at all tilt values the system produces a pair of relativistic jets. At small distances the jets precess rapidly together with the inner sub-disc, whereas at large distances they partially align with the outer sub-disc and precess more slowly. If the tearing radius can be modeled accurately in future work, emission model independent measurements of black hole spin based on precession-driven quasi-periodic oscillations may become possible.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
