No Arabic abstract
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}$.
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.
In quiescence, Sgr A* is surprisingly dim, shining 100,000 times less than expected for its environment. This problem has motivated a host of theoretical models to explain radiatively inefficient accretion flows (RIAFs). The Chandra Galactic Center (GC) X-ray Visionary Program obtained approximately 3 Ms (one month) of Chandra HETG data, offering the only opportunity to examine the quiescent X-ray emission of Sgr A* with high resolution spectroscopy. Utilizing custom background regions and filters for removing overlapping point sources, this work provides the first ever look at stacked HETG spectra of Sgr A*. We model the background datasets with a cubic spline and fit the unbinned Sgr A* spectra with a simple parametric model of a power law plus Gaussian lines under the effects of interstellar extinction. We detect a strong 6.7 keV iron emission line in the HEG spectra and a 3.1 keV emission line in the MEG spectra. In all cases, the line centroids and equivalent widths are consistent with those measured from low-resolution CCD spectra. An examination of the unbinned, stacked HEG+/-1 spectrum reveals fine structure in the iron line complex. In addition to resolving the resonant and forbidden lines from He-like iron, there are apparent emission features arising with higher statistical significance at lower energy, potentially associated with FeXX-XXIV ions in a ~1 keV plasma arising near the Bondi radius of Sgr A*. With this work, we release the cleaned and stacked Sgr A* and background HETG spectra to the public as a special legacy dataset.
In this study we present three-dimensional radiative cooling hydrodynamical simulations of galactic winds generated particularly in M82-like starburst galaxies. We have considered intermittent winds induced by SNe explosions within super star clusters randomly distributed in the central region of the galaxy and were able to reproduce the observed M82 wind conditions with its complex morphological outflow structure. We have found that the environmental conditions in the disk in nearly recent past are crucial to determine whether the wind will develop a large scale rich filamentary structure, as in M82 wind, or not. Also, the numerical evolution of the SN ejecta have allowed us to obtain the abundance distribution over the first 3 kpc extension of the wind and we have found that the SNe explosions change significantly the metallicity only of the hot, low-density wind component. Moreover, we have found that the SN-driven wind transports to outside the disk large amounts of energy, momentum and gas, but the more massive high-density component reaches only intermediate altitudes smaller than 1.5 kpc. Therefore, no significant amounts of gas mass are lost to the IGM and the mass evolution of the galaxy is not much affected by the starburst events occurring in the nuclear region.
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.