On the accuracy of common moment-based radiative transfer methods for simulating reionization

Abstract

Cosmological simulations of reionization often treat radiative transfer by solving for the monopole and dipoles of the intensity field and by making ansatz for the quadrupole moments to close the system of equations. We investigate the accuracy of the most common closure methods, i.e. Eddington tensor choices. We argue that these algorithms are most likely to err after reionization and study quasi-analytic test problems mimicking these situations: large-scale post-reionization ionizing background fluctuations and radiative transfer in a predominantly ionized medium with discrete absorbers. We show that OTVET and M1 over-ionize self-shielding absorbers when fixing the background photoionization rate, leading to 30-40% higher emissivity to balance the increased recombination. This over-ionization results in a simulation run with these algorithms having a factor of ~2 lower average metagalactic photoionization rate relative to truth given an ionizing emissivity. Furthermore, these algorithms are unlikely to reproduce ionizing background fluctuations on scales below the photon mean path: OTVET tends to overpredict the fluctuations there when the simulation box is smaller than twice the mean free path and underpredict otherwise, while M1 drastically underpredicts these fluctuations. As a result, these numerical methods are likely not sufficiently accurate to interpret the Ly$\alpha$ forest opacity fluctuations observed after reionzation. We show that a high number of angular directions need to be followed to capture the post-reionization ionizing background fluctuations accurately with ray-tracing codes. Lastly, we argue that the strong dependence of the post-reionization ionizing background on the value of the reduced speed of light found in many simulations signals that the ionizing photon mean free path is several times larger in such simulations than the observationally measured value.

revision 1, submitted to JCAP

References in corpus (25)

• Model-independent evidence in favor of an end to reionization by z~6
• Evidence of patchy hydrogen reionization from an extreme Ly$\alpha$ trough below redshift six
• The morphology of HII regions during reionization
• On the (Non)Evolution of HI Disks over Cosmic Time
• Radiative transfer simulations of cosmic reionization I: methodology and initial results
• Cosmological Radiative Transfer Codes Comparison Project I: The Static Density Field Tests
• Simulating Cosmic Reionization at Large Scales II: the 21-cm Emission Features and Statistical Signals
• A scheme for radiation pressure and photon diffusion with the M1 closure in RAMSES-RT
• The Neutral Hydrogen Cosmological Mass Density at z=5
• A radiative transfer scheme for cosmological reionization based on a local Eddington tensor
• Imprint of Inhomogeneous Hydrogen Reionization on the Temperature Distribution of the Intergalactic Medium
• TRAPHIC - Radiative Transfer for Smoothed Particle Hydrodynamics Simulations
• An implementation of radiative transfer in the cosmological simulation code GADGET
• Scale-dependent bias in the BAO-scale intergalactic neutral hydrogen
• Hydrodynamic Response of the Intergalactic Medium to Reionization
• Multi-frequency, thermally coupled radiative transfer with TRAPHIC: Method and tests
• A New Moment Method for Continuum Radiative Transfer in Cosmological Reionization
• Statistical properties of damped Lyman-alpha systems from Sloan Digital Sky Survey DR12
• The ionizing background at the end of overlap
• The Acoustic Peak in the Lyman Alpha Forest
• Simulations of galaxy formation with radiative transfer: Hydrogen reionization and radiative feedback
• Simulating the effect of photoheating feedback during reionization
• On the Proper Use of the Reduced Speed of Light Approximation
• The influence of metagalactic ultra-violet background fluctuations on the high-redshift Ly$\alpha$ forest
• Revised estimates of CMB $B$-mode polarization induced by patchy reionization

Cited by in corpus (8)

• Unveiling Dark Matter free-streaming at the smallest scales with high redshift Lyman-alpha forest
• A short mean free path at $z=6$ favors late and rapid reionization by faint galaxies
• Lyman-alpha opacities at z=4-6 require low mass, radiatively-suppressed galaxies to drive cosmic reionization
• Reionization with Simba: How much does astrophysics matter in modeling cosmic reionization?
• Late-end reionization with ATON-HE: towards constraints from Lyman-$\alpha$ emitters observed with JWST
• The impact and response of minihalos and the inter-halo medium on cosmic reionization
• Measuring photo-ionization rate and mean free path of HeII ionizing photons at $2.5\le z\le 3.6$: Evidence for late and rapid HeII reionization Part-II
• The 21-cm signature of X-ray heated halos around galaxies during cosmic dawn