A WENO Algorithm of the Temperature and Ionization Profiles around a Point Source

Abstract

We develop a numerical solver for radiative transfer problems based on the weighted essentially nonoscillatory (WENO) scheme modified with anti-diffusive flux corrections, in order to solve the temperature and ionization profiles around a point source of photons in the reionization epoch. Algorithms for such simulation must be able to handle the following two features: 1. the sharp profiles of ionization and temperature at the ionizing front (I-front) and the heating front (T-front), and 2. the fraction of neutral hydrogen within the ionized sphere is extremely small due to the stiffness of the rate equations of atom processes. The WENO scheme can properly handle these two features, as it has been shown to have high order of accuracy and good convergence in capturing discontinuities and complicated structures in fluid as well as to be significantly superior over piecewise smooth solutions containing discontinuities. With this algorithm, we show the time-dependence of the preheated shell around a UV photon source. In the first stage the I-front and T-front are coincident, and propagate with almost the speed of light. In later stage, when the frequency spectrum of UV photons is hardened, the speeds of propagation of the ionizing and heating fronts are both significantly less than the speed of light, and the heating front is always beyond the ionizing front. In the spherical shell between the I- and T-fronts, the IGM is heated, while atoms keep almost neutral. The time scale of the preheated shell evolution is dependent on the intensity of the photon source. We also find that the details of the pre-heated shell and the distribution of neutral hydrogen remained in the ionized sphere are actually sensitive to the parameters used. The WENO algorithm can provide stable and robust solutions to study these details.

24 pages, 7 figures, accepted in New Astronomy

References in corpus (13)

• P3M-SPH simulations of the Lyman-alpha Forest
• Photoevaporation of Cosmological Minihalos during Reionization
• Radiation Hydrodynamical Evolution of Primordial H II Regions
• The Structure and Evolution of Early Cosmological HII Regions
• Radio Signatures of HI at High Redshift: Mapping the End of the ``Dark Ages''
• The H II Region of the First Star
• Cosmological Radiative Transfer Codes Comparison Project I: The Static Density Field Tests
• Cosmological Reionization Around the First Stars: Monte Carlo Radiative Transfer
• Cosmic Hydrogen Was Significantly Neutral a Billion Years After the Big Bang
• Cosmological Hydrogen Reionization with Three Dimensional Radiative Transfer
• Prospects for Redshifted 21-cm observations of quasar HII regions
• 3D numerical cosmological radiative transfer in an inhomogeneous medium
• The apparent shape of the "Strömgren sphere'' around the highest-redshift QSOs with Gunn-Peterson troughs

Cited by in corpus (9)

• The Molecular Hydrogen Deficit in Gamma-Ray Burst Afterglows
• Regulated star formation in forming disk galaxies under ultraviolet radiation background
• Time evolution of Wouthuysen-Field coupling
• Resonant Scattering and Ly-alpha Radiation Emergent from Neutral Hydrogen Halos
• A WENO algorithm for radiative transfer with resonant scattering: the time scale of the Wouthuysen-Field Coupling
• Wouthuysen-Field Coupling in the 21 cm Region Around High Redshift Sources
• Cosmological radiative transfer for the line-of-sight proximity effect
• 21 cm Signals from Early Ionizing Sources
• Effect of Dust on Lyman-alpha Photon Transfer in Optically Thick Halo