thor: a massively-parallel GPU-accelerated MCRT code

Volume rendering of the temperature field for a massive node in a galaxy formation simulation using thor.

Emission and absorption line features are important diagnostics for the physics underlying extragalactic astronomy. The interpretation of observed signatures involves comparing against forward modeled spectra from galaxy formation simulations as well as more simplified geometries, while including the complex scattering radiative transfer (RT) of resonant emission lines. We introduce thor, a modern C++ radiative transfer code focused initially on resonant emission lines. thor is a high-performance, distributed memory MPI-parallel, multi-target code, running on CPUs, GPUs and other accelerators, yielding large $\sim 10-50\rm{x}$ speed-ups compared to previous CPU-only codes.

Benchmark on a common Lyman-alpha radiative transfer ‘Neufeld solution’ (T=10⁎K, 𝜏₀=10⁶) toy problem for different resonant emission MCRT codes. Compiled with Clang 18. CPU: Intel 8360Y with 36 cores at 2.40GHz, GPU: Nvidia A100 40GB. Tested for latest code versions in July 2025.
Lyman-alpha radiative transfer weak scaling test for number of volume elements for a cosmological galaxy formation simulation for up to ~200 billion volume elements, memory distributed with up to 15k CPU cores and 200 GPUs.

To be published. Pre-print available.

Chris Byrohl
Chris Byrohl
Postdoc

My research interests include cosmoslogical galaxy formation simulations, Lyman-$\alpha$ radiation to study galaxies and the large-scale structure, supernovae type Ia and high-performance computing.