Abstract:Data derived from general relativistic MHD simulations of accretion onto black holes can be used as input to a post-processing scheme that predicts the radiated spectrum. Combining a relativistic Compton-scattering radiation transfer solution in the corona with detailed local atmosphere solutions incorporating local ionization and thermal balance within the disk photosphere, it is possible to study both spectral formation and intrinsic spectral variability in the radiation from relativistic accretion disks. With this method, we find that radiatively-efficient systems with black holes of $10M_\odot$ accreting at $\approx 0.01$ in Eddington units produce spectra very similar to those observed in the hard states of X-ray binaries. The spectral shape above 10keV is well-described by a power-law with an exponential cutoff. Intrinsic turbulent variations lead to order-unity fluctuations in bolometric luminosity, logarithmic spectral slope fluctuations over a range $\approx 0.15$, and factor of 2 changes in the cut-off energy on timescales $\sim 50\, (M_{\rm BH}/10 M_\odot)$ms. Within the corona, the range of gas temperature spans more than an order of magnitude. The wide distribution of temperatures is central to defining the spectrum: the logarithmic spectral slope is harder by $\sim 0.3$ and the cut-off energy larger by a factor $\sim 10 - 30$ than if the coronal temperature everywhere were its mass-weighted mean.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to predict the X - ray spectra of black holes and their spectral variations through simulation - based methods. Specifically, the authors use the data from general - relativistic magnetohydrodynamics (GRMHD) simulations as input, combined with the relativistic Compton scattering radiative transfer solution and the local atmosphere solution of the disk, to study the spectral formation process in the relativistic accretion disk and the inherent spectral variation characteristics. This method can predict the spectra produced when a black hole with a mass of about 10 solar masses accretes at about 0.01 times the Eddington unit in a highly efficient radiation system, and these spectra are very similar to the hard - state observations of X - ray binaries.
### Main problems:
1. **Spectral prediction**: How to directly predict the spectral shape of this radiation and its variation characteristics from dynamic calculations?
2. **Spectral variation**: How to explain the variation of these spectra over time and space, especially the spectral variation of hard X - ray radiation?
### Method overview:
- **GRMHD simulation**: Use the HARM3D code to generate a three - dimensional map of the accretion material around the black hole, including the rest - mass density, internal energy density, 4 - velocity, and cooling rate.
- **Post - processing**: Use Pandurata to post - process the 3D simulation snapshots, calculate the electron temperature and photon spectral intensity, and solve the local thermal equilibrium through the Monte Carlo radiative transfer code.
- **Radiative transfer**: Combine the one - dimensional radiative transfer solution to calculate the opacity and emissivity of the entire disk, and finally obtain the radiative spectrum of the system.
### Key findings:
- **Spectral characteristics**: The predicted spectra above 10 keV can be described by a power law plus an exponential cutoff, and the power - law index and cutoff energy vary with time.
- **Spectral variation**: The intrinsic turbulence variation leads to the fluctuation of the luminosity across the entire band. The fluctuation range of the logarithmic spectral slope is about 0.15, the cutoff energy varies by a factor of about 2, and the time scale is about 50 (MBH/10M⊙) milliseconds.
- **Temperature distribution**: The gas temperature range in the corona region exceeds one order of magnitude, which determines the shape of the spectrum: compared with the case where the corona region temperature is assumed to be uniform, the logarithmic spectral slope is about 0.3 steeper, and the cutoff energy is about 10 - 30 times larger.
### Conclusions:
- **Physical mechanism**: By considering the inhomogeneity of the internal structure of the accretion flow and the variability driven by dynamics, this method can more accurately predict the spectra of hard X - ray radiation and their variations.
- **Comparison with observations**: The predicted spectral slope and cutoff energy are consistent with the observations of X - ray binaries and active galactic nuclei, supporting the view that a highly efficient radiation accretion disk is accompanied by a corona region.
Through these studies, the authors not only provide a new understanding of the X - ray radiation mechanism of black holes but also provide new tools and methods for future research.