2012ApJ...752...18K

Radiation spectra of supercritical black hole accretion flows are computed using a Monte Carlo method by post-processing the results of axisymmetric radiation hydrodynamic simulations. We take into account thermal/bulk Comptonization, free-free absorption, and photon trapping. We found that a shock-heated region (∼10⁸ K) appears at the funnel wall near the black hole where the supersonic inflow is reflected by the centrifugal barrier of the potential. Both thermal and bulk Comptonization significantly harden photon spectra although most of the photons upscattered above 40 keV are swallowed by the black hole due to the photon trapping. When the accretion rate onto the black hole is {dot}M ~ 200L_E/c ², where L_E is the Eddington luminosity, the spectrum has a power-law component which extends up to ∼10 keV by upscattering of photons in the shock-heated region. In higher mass accretion rates, the spectra roll over around 5 keV due to downscattering of the photons by cool electrons in the dense outflow surrounding the jet. Our results are consistent with the spectral features of ultraluminous X-ray sources, which typically show either a hard power-law component extending up to 10 keV or a rollover around 5 keV. We found that the spectrum of NGC 1313 X-2 is quite similar to the spectrum numerically obtained for high accretion rate ({dot}M ~ 200L_E/c ²) source observed with low viewing angle (i = 10°-20°). Our numerical results also demonstrate that the face-on luminosity of supercritically accreting stellar mass black holes (10 M_☉) can significantly exceed 10⁴⁰ erg/s.