Spectral Calculations of 3D Radiation Magnetohydrodynamic Simulations of Super-Eddington Accretion onto a Stellar-mass Black Hole

Abstract

We use the Athena++ Monte Carlo (MC) radiation transfer module to postprocess simulation snapshots from nonrelativistic Athena++ radiation magnetohydrodynamic (RMHD) simulations. These simulations were run using a gray (frequency-integrated) approach but were also restarted and ran with a multigroup approach that accounts for Compton scattering with a Kompaneets operator. These simulations produced moderately super-Eddington accretion rates onto a 6.62 M ⊙ black hole. Since we only achieve inflow equilibrium out to 20–25 gravitational radii, we focus on the hard X-ray emission. We provide a comparison between the MC and RMHD simulations, showing that the treatment of Compton scattering in the gray RMHD simulations underestimates the gas temperature in the funnel regions above and below the accretion disk. In contrast, the restarted multigroup snapshots provide a treatment for the radiation field that is more consistent with the MC calculations, and result in postprocessed spectra with harder X-ray emission compared to their gray snapshot counterparts. We characterize these MC postprocessed spectra using commonly employed phenomenological spectral fitting models. We also attempt to fit our MC spectra directly to observations of the ultraluminous X-ray source (ULX) NGC 1313 X-1, finding best-fit values that are competitive to phenomenological model fits, indicating that first principle models of super-Eddington accretion may adequately explain the observed hard X-ray spectra in some ULX sources.