Avery Meiksin

The strength of the 21 cm Cosmic Dawn radio absorption signature against the Cosmic Microwave Background or against a bright background radio source depends on the temperature of the neutral hydrogen in the intergalactic medium. While X-rays from forming galaxies will likely dominate the heating rate, recent models suggest scenarios in which heating by the scattering of Ly α photons sourced by the galaxies contributes non-negligibly as well, and may even dominate. The efficiency of Ly α photon heating for a point source at high redshift using a numerical solution to the exact radiative transfer equation is provided here. It is found to be a factor of several smaller than the rate computed in the commonly used diffusion approximation to the radiative transfer equation.

We assess the impact of QSOs on the high redshift (z > 4) Intergalactic Medium using Monte Carlo realisations of QSO populations and the HeIII regions they generate, applied to the Sherwood-Relics simulations, allowing for uncertainties in the QSO luminosity function, its evolution, and QSO spectra and ages. While QSO luminosity functions based on optical-infrared selection are unable to reproduce the broadening HI Lyman-alpha optical depth distributions at z > 5, much broader distributions are found for the higher numbers of QSOs based on x-ray selection, suggesting a large QSO contribution to the ultra-violet background at z > 5 may offer an alternative to late reionization models to account for the broad HI Lyman-alpha optical depth distributions. Realisations using QSOs based on the higher QSO counts also much better recover the measured pixel flux auto-correlation function at z > 5. The HeIII regions from QSO sources according to both types of luminosity function suppress the pixel flux power spectrum on small scales, k > 0.02 s/km, while enhancing it on larger, both by amounts of up to tens of percent at z > 4, with the magnitude increasing with the intergalactic HeIII filling factor and the boost in temperature within the HeIII regions.

The scattering of Lyα photons from the first radiating sources in the universe plays a pivotal role in 21 cm radio detections of Cosmic Dawn and the Epoch of Reionization through the Wouthuysen–Field effect. New data from JWST show the Lyα photon scattering rate exceeds that required to decouple the intergalactic hydrogen spin temperature from that of the Cosmic Microwave Background up to z ∼ 14 and render the neutral hydrogen visible.

We present the Sherwood-Relics simulations, a new suite of large cosmological hydrodynamical simulations aimed at modelling the intergalactic medium (IGM) during and after the cosmic reionization of hydrogen. The suite consists of over 200 simulations that cover a wide range of astrophysical and cosmological parameters. It also includes simulations that use a new lightweight hybrid scheme for treating radiative transfer effects. This scheme follows the spatial variations in the ionizing radiation field, as well as the associated fluctuations in IGM temperature and pressure smoothing. It is computationally much cheaper than full radiation hydrodynamics simulations and circumvents the difficult task of calibrating a galaxy formation model to observational constraints on cosmic reionization. Using this hybrid technique, we study the spatial fluctuations in IGM properties that are seeded by patchy cosmic reionization. We investigate the relevant physical processes and assess their impact on the z > 4 Lyman-alpha forest. Our main findings are: (i) Consistent with previous studies patchy reionization causes large scale temperature fluctuations that persist well after the end of reionization, (ii) these increase the Lyman-alpha forest flux power spectrum on large scales, and (iii) result in a spatially varying pressure smoothing that correlates well with the local reionization redshift. (iv) Structures evaporated or puffed up by photoheating cause notable features in the Lyman-alpha forest, such as flat-bottom or double-dip absorption profiles.

We compare numerical methods for solving the radiative transfer equation in the context of the photoionization of intergalactic gaseous hydrogen and helium by a central radiating source. Direct integration of the radiative transfer equation and solutions using photon packets are examined, both for solutions to the time-dependent radiative transfer equation and in the infinite-speed-of-light approximation. The photon packet schemes are found to be more generally computationally efficient than a direct integration scheme. Whilst all codes accurately describe the growth rate of hydrogen and helium ionization zones, it is shown that a fully time-dependent method is required to capture the gas temperature and ionization structure in the near zone of a source when an ionization front expands at a speed close to the speed of light. Applied to Quasi-Stellar Objects in the Epoch of Reionization (EoR), temperature differences as high as 5 × 10⁴ K result in the near-zone for solutions of the time-dependent radiative transfer equation compared with solutions in the infinite-speed-of-light approximation. Smaller temperature differences are found following the nearly full photoionization of helium in gas in which the hydrogen was already ionized and the helium was singly ionized. Variations found in the temperature and ionization structure far from the source, where the gas is predominantly neutral, may affect some predictions for 21-cm EoR experiments.

We present the Sherwood-Relics simulations, a new suite of large cosmological hydrodynamical simulations aimed at modelling the intergalactic medium (IGM) during and after the cosmic reionization of hydrogen. The suite consists of over 200 simulations that cover a wide range of astrophysical and cosmological parameters. It also includes simulations that use a new lightweight hybrid scheme for treating radiative transfer effects. This scheme follows the spatial variations in the ionizing radiation field, as well as the associated fluctuations in IGM temperature and pressure smoothing. It is computationally much cheaper than full radiation hydrodynamics simulations and circumvents the difficult task of calibrating a galaxy formation model to observational constraints on cosmic reionization. Using this hybrid technique, we study the spatial fluctuations in IGM properties that are seeded by patchy cosmic reionization. We investigate the relevant physical processes and assess their impact on the z > 4 Lyman-α forest. Our main findings are: (i) Consistent with previous studies patchy reionization causes large scale temperature fluctuations that persist well after the end of reionization, (ii) these increase the Lyman-α forest flux power spectrum on large scales, and (iii) result in a spatially varying pressure smoothing that correlates well with the local reionization redshift. (iv) Structures evaporated or puffed up by photoheating cause notable features in the Lyman-α forest, such as flat-bottom or double-dip absorption profiles.