pyC${}^2$Ray: A flexible and GPU-accelerated Radiative Transfer Framework for Simulating the Cosmic Epoch of Reionization

Abstract

Detailed modeling of the evolution of neutral hydrogen in the intergalactic medium during the Epoch of Reionization, $5\le z\le 20$, is critical in interpreting the cosmological signals from current and upcoming 21-cm experiments such as the Low-Frequency Array (LOFAR) and the Square Kilometre Array (SKA). Numerical radiative transfer codes provide the most physically accurate models of the reionization process. However, they are computationally expensive as they must encompass enormous cosmological volumes while accurately capturing astrophysical processes occurring at small scales ($\lesssim \mathrm{Mpc}$). Here, we present pyC${}^2$Ray, an updated version of the massively parallel ray-tracing and chemistry code, C${}^2$Ray, which has been extensively employed in reionization simulations. The most time-consuming part of the code is calculating the hydrogen column density along the path of the ionizing photons. Here, we present the Accelerated Short-characteristics Octahedral ray-tracing (ASORA) method, a ray-tracing algorithm specifically designed to run on graphical processing units (GPUs). We include a modern Python interface, allowing easy and customized use of the code without compromising computational efficiency. We test pyC${}^2$Ray on a series of standard ray-tracing tests and a complete cosmological simulation with volume size $(349\mathrm{Mpc}{)}^3$, mesh size of $250^3$ and approximately $10^6$ sources. Compared to the original code, pyC${}^2$Ray achieves the same results with negligible fractional differences, $\sim 10^{-5}$, and a speedup factor of two orders of magnitude. Benchmark analysis shows that ASORA takes a few nanoseconds per source per voxel and scales linearly for an increasing number of sources and voxels within the ray-tracing radii.

20 pages, 13 figures, 1 table

References in corpus (15)

• Array Programming with NumPy
• The Astropy Project: Sustaining and Growing a Community-oriented Open-source Project and the Latest Major Release (v5.0) of the Core Package
• Cosmology at Low Frequencies: The 21 cm Transition and the High-Redshift Universe
• Grackle: a Chemistry and Cooling Library for Astrophysics
• Improved upper limits on the 21-cm signal power spectrum of neutral hydrogen at $\mathbf{z}\approx 9.1$ from LOFAR
• Deep multi-redshift limits on Epoch of Reionisation 21cm Power Spectra from Four Seasons of Murchison Widefield Array Observations
• Cosmological Radiative Transfer Codes Comparison Project I: The Static Density Field Tests
• A radiative transfer scheme for cosmological reionization based on a local Eddington tensor
• TRAPHIC - Radiative Transfer for Smoothed Particle Hydrodynamics Simulations
• Lyman-alpha radiative transfer during the Epoch of Reionization: contribution to 21-cm signal fluctuations
• Bubble size statistics during reionization from 21-cm tomography
• SPHRAY: A Smoothed Particle Hydrodynamics Ray Tracer for Radiative Transfer
• Turbulence in simulated HII regions
• The Diffuse Nature of Stromgren Spheres
• Physics of Cosmic Reionization

Cited by in corpus (5)

• Deep learning approach for identification of HII regions during reionization in 21-cm observations -- III. image recovery
• Impact of anisotropic photon emission from sources during the epoch of reionisation
• Mapping neutral islands during end stages of reionization with photometric intergalactic medium tomography
• CosmoUiT: A Vision Transformer-UNet Hybrid for Fast and Accurate Emulation of 21-cm Maps from the Epoch of Reionization
• Reionization in HESTIA: Studying reionization in the LG through zoom simulations