Discrete diffusion Lyman-alpha radiative transfer

Abstract

Due to its accuracy and generality, Monte Carlo radiative transfer (MCRT) has emerged as the prevalent method for Ly$\alpha$ radiative transfer in arbitrary geometries. The standard MCRT encounters a significant efficiency barrier in the high optical depth, diffusion regime. Multiple acceleration schemes have been developed to improve the efficiency of MCRT but the noise from photon packet discretization remains a challenge. The discrete diffusion Monte Carlo (DDMC) scheme has been successfully applied in state-of-the-art radiation hydrodynamics (RHD) simulations. Still, the established framework is not optimal for resonant line transfer. Inspired by the DDMC paradigm, we present a novel extension to resonant DDMC (rDDMC) in which diffusion in space and frequency are treated on equal footing. We explore the robustness of our new method and demonstrate a level of performance that justifies incorporating the method into existing Ly$\alpha$ codes. We present computational speedups of $\sim 10^2$-$10^6$ relative to contemporary MCRT implementations with schemes that skip scattering in the core of the line profile. This is because the rDDMC runtime scales with the spatial and frequency resolution rather than the number of scatterings - the latter is typically $\propto {\tau }_0$ for static media, or $\propto (a{\tau }_0{)}^{2/3}$ with core-skipping. We anticipate new frontiers in which on-the-fly Ly$\alpha$ radiative transfer calculations are feasible in 3D RHD. More generally, rDDMC is transferable to any computationally demanding problem amenable to a Fokker-Planck approximation of frequency redistribution.

15 pages, 13 figures, MNRAS, in press

References in corpus (10)

• Lyman Alpha Emitting Galaxies as a Probe of Reionization
• The Formation of the First Stars II. Radiative Feedback Processes and Implications for the Initial Mass Function
• Lyman alpha Radiative Transfer in Cosmological Simulations using Adaptive Mesh Refinement
• The Lyman-alpha signature of the first galaxies
• Radiation Transport for Explosive Outflows: Opacity Regrouping
• Lyman Alpha Driven Outflows Around Star Forming Galaxies
• Lyman-alpha radiation hydrodynamics of galactic winds before cosmic reionization
• Radiation pressure in super star cluster formation
• Assembly of supermassive black hole seeds
• Radiative effects during the assembly of direct collapse black holes

Cited by in corpus (16)

• Radiative Stellar Feedback in Galaxy Formation: Methods and Physics
• The physics of Lyman-alpha escape from high-redshift galaxies
• Axisymmetric Radiative Transfer Models of Kilonovae
• Impact of Lyman alpha pressure on metal-poor dwarf galaxies
• Intensity mapping with SDSS/BOSS Lyman-alpha emission, quasars and their Lyman-alpha forest
• Supermassive Black Holes in the Early Universe
• Radiation pressure in super star cluster formation
• Arepo-MCRT: Monte Carlo Radiation Hydrodynamics on a Moving Mesh
• Lyman-$\alpha$ feedback prevails at Cosmic Dawn: Implications for the first galaxies, stars, and star clusters
• SILVERRUSH. XI. Intensity Mapping for Lya Emission Extending over $100-1000$ comoving kpc around $z\sim 2-7$ LAEs with Subaru HSC-SSP and CHORUS Data
• Resonant-line radiative transfer within power-law density profiles
• Lyman-alpha radiation pressure: an analytical exploration
• Lyman-alpha resonant-line radiative transfer in expanding media
• A Novel Solution for Resonant Scattering Using Self-Consistent Boundary Conditions
• Influence of turbulence on Lyman-alpha scattering
• Lyman-$\alpha$ radiation pressure regulates star formation efficiency