grmonty: a Monte Carlo Code for Relativistic Radiative Transport

We describe a Monte Carlo radiative transport code intended for calculating spectra of hot, optically thin plasmas in full general relativity. The version we describe here is designed to model hot accretion flows in the Kerr metric and therefore incorporates synchrotron emission and absorption, and Compton scattering. The code can be readily generalized, however, to account for other radiative processes and an arbitrary spacetime. We describe a suite of test problems, and demonstrate the expected $N^{-1/2}$ convergence rate, where $N$ is the number of Monte Carlo samples. Finally we illustrate the capabilities of the code with a model calculation, a spectrum of the slowly accreting black hole Sgr A* based on data provided by a numerical general relativistic MHD model of the accreting plasma.

💡 Research Summary

The paper introduces grmonty, a Monte Carlo radiative‑transfer code designed to compute spectra from hot, optically thin plasmas in full general relativity. Unlike earlier 1‑D or non‑relativistic Monte Carlo tools, grmonty integrates photon geodesics directly in the Kerr metric, thereby capturing all relativistic effects (gravitational redshift, light bending, frame dragging) that are essential near rotating black holes. The code is built around three logical modules: photon generation, propagation, and interaction. Photon generation implements synchrotron emissivity and absorptivity assuming a locally thermal (Maxwell‑Jüttner) electron distribution; the emissivity depends on the magnetic field strength, electron temperature, and viewing angle. Absorption is treated as the inverse process, using the same coefficients to ensure detailed balance. Interaction is limited to Compton scattering, for which the full Klein‑Nishina cross‑section (the exact Cohn‑Simpson kernel) is sampled, allowing accurate energy and angle redistribution even in the ultra‑relativistic regime.

Propagation solves the null geodesic equations using a fourth‑order Runge‑Kutta integrator with adaptive step control. At each step the code queries the underlying fluid model (density, temperature, magnetic field, four‑velocity) and decides probabilistically whether an emission, absorption, or scattering event occurs, based on the local optical depth. This design makes it straightforward to replace the Kerr background with any arbitrary metric, or to import time‑dependent data from general‑relativistic magnetohydrodynamic (GRMHD) simulations.

The authors validate grmonty with a suite of four benchmark problems. (1) Vacuum geodesic propagation reproduces analytically known photon trajectories in Kerr spacetime. (2) Synchrotron spectra from a uniform magnetized plasma match the textbook power‑law dependence on frequency. (3) Single‑scattering tests confirm that the post‑scattering photon energy distribution follows the Klein‑Nishina prediction. (4) A combined synchrotron‑absorption‑Compton scenario demonstrates that the total emergent spectrum converges to the analytic solution as the number of Monte Carlo samples N increases. In all cases the statistical error scales as N⁻¹ᐟ², confirming the expected Monte Carlo convergence.

To showcase astrophysical relevance, the code is applied to a realistic model of the supermassive black hole Sgr A*. The authors import three‑dimensional GRMHD data (density, temperature, magnetic field, velocity) from a low‑accretion‑rate (ADAF‑like) simulation of the Galactic centre. Using these fields, grmonty computes a broadband spectrum spanning radio, sub‑millimetre, infrared, and X‑ray frequencies. The resulting spectral energy distribution reproduces the observed radio/sub‑mm peak (synchrotron‑dominated) and the high‑energy X‑ray tail (produced by inverse‑Compton up‑scattering of synchrotron photons). The agreement with observations demonstrates that the code can bridge the gap between sophisticated GRMHD dynamics and observable radiation signatures.

Key contributions of the work are: (i) a fully relativistic photon‑geodesic engine in the Kerr spacetime; (ii) an integrated treatment of synchrotron emission/absorption and Compton scattering within a Monte Carlo framework; (iii) a modular architecture that facilitates extension to other metrics, non‑thermal electron distributions, or additional radiative processes such as bremsstrahlung or pair production; (iv) thorough verification of numerical accuracy and statistical convergence; and (v) a proof‑of‑concept astrophysical application to Sgr A*. The authors suggest future extensions including non‑thermal electron power‑law tails, radiative transfer in non‑Kerr spacetimes (e.g., modified gravity or binary black‑hole backgrounds), and the inclusion of photon‑photon interactions. Such developments would enable high‑fidelity modeling of a wide range of relativistic accretion phenomena, from low‑luminosity galactic nuclei to luminous quasars and tidal‑disruption events.