-
The Lyman-limit photon mean free path at the end of late reionisation in the Sherwood-Relics simulations
Authors:
Jennifer Feron,
Luke Conaboy,
James S. Bolton,
Emma Chapman,
Martin G. Haehnelt,
Laura C. Keating,
Girish Kulkarni,
Ewald Puchwein
Abstract:
Recent evidence supporting reionisation ending at redshift z<6 includes the rapid redshift evolution of the mean free path for Lyman-limit photons through the intergalactic medium at 5<z<6. Here we investigate the mean free path predicted by the Sherwood-Relics suite of hybrid radiation hydrodynamical simulations. Simulations with comoving volumes of 40^3 h^-3 cMpc^3 (160^3 h^-3 cMpc^3), calibrate…
▽ More
Recent evidence supporting reionisation ending at redshift z<6 includes the rapid redshift evolution of the mean free path for Lyman-limit photons through the intergalactic medium at 5<z<6. Here we investigate the mean free path predicted by the Sherwood-Relics suite of hybrid radiation hydrodynamical simulations. Simulations with comoving volumes of 40^3 h^-3 cMpc^3 (160^3 h^-3 cMpc^3), calibrated to match the observed Lyman-alpha forest transmission with a late end to reionisation at z<6, are consistent with recent mean free path measurements at z<5.9, and are 1.2σ(1.8σ) above the highest redshift mean free path measurement at z=5.93. The majority of the Lyman-limit opacity at the end of reionisation is attributable to highly ionised Lyman-alpha forest absorbers with neutral hydrogen column densities N_HI=10^16--10^17 cm^-2. Resolving these systems is critical for capturing the redshift evolution of the mean free path during the final stages of reionisation. After reionisation completes, overdense gas will reduce the mean free path by up to 20 per cent around haloes with masses M_h ~ 10^9--10^11 h^-1 M_sol, but during reionisation ionised bubbles will instead boost the mean free path around haloes by up to an order of magnitude when the IGM is as much as 90 per cent neutral by volume. This effect will play an important role in the visibility of Lyman-alpha emitting galaxies at z>10 discovered with JWST.
△ Less
Submitted 27 June, 2024;
originally announced June 2024.
-
Relative baryon-dark matter velocities in cosmological zoom simulations
Authors:
Luke Conaboy,
Ilian T. Iliev,
Anastasia Fialkov,
Keri L. Dixon,
David Sullivan
Abstract:
Supersonic relative motion between baryons and dark matter due to the decoupling of baryons from the primordial plasma after recombination affects the growth of the first small-scale structures. Large box sizes (greater than a few hundred Mpc) are required to sample the full range of scales pertinent to the relative velocity, while the effect of the relative velocity is strongest on small scales (…
▽ More
Supersonic relative motion between baryons and dark matter due to the decoupling of baryons from the primordial plasma after recombination affects the growth of the first small-scale structures. Large box sizes (greater than a few hundred Mpc) are required to sample the full range of scales pertinent to the relative velocity, while the effect of the relative velocity is strongest on small scales (less than a few hundred kpc). This separation of scales naturally lends itself to the use of `zoom' simulations, and here we present our methodology to self-consistently incorporate the relative velocity in zoom simulations, including its cumulative effect from recombination through to the start time of the simulation. We apply our methodology to a large-scale cosmological zoom simulation, finding that the inclusion of relative velocities suppresses the halo baryon fraction by $46$--$23$ per cent between $z=13.6$ and $11.2$, in qualitative agreement with previous works. In addition, we find that including the relative velocity delays the formation of star particles by $\sim 20 {~\rm Myr}$ Myr on average (of the order of the lifetime of a $\sim 9~{\rm M}_\odot$ Population III star) and suppresses the final stellar mass by as much as $79$ per cent at $z=11.2$.
△ Less
Submitted 14 September, 2023; v1 submitted 23 July, 2022;
originally announced July 2022.
-
The short ionizing photon mean free path at z=6 in Cosmic Dawn III, a new fully-coupled radiation-hydrodynamical simulation of the Epoch of Reionization
Authors:
Joseph S. W. Lewis,
Pierre Ocvirk,
Jenny G. Sorce,
Yohan Dubois,
Dominique Aubert,
Luke Conaboy,
Paul R. Shapiro,
Taha Dawoodbhoy,
Romain Teyssier,
Gustavo Yepes,
Stefan Gottlöber,
Yann Rasera,
Kyungjin Ahn,
Ilian T. Iliev,
Hyunbae Park,
Émilie Thélie
Abstract:
Recent determinations of the mean free path of ionising photons (mfp) in the intergalactic medium (IGM) at $\rm z=6$ are lower than many theoretical predictions. To gain insight into this issue, we investigate the evolution of the mfp in our new massive fully coupled radiation hydrodynamics cosmological simulation of reionization: Cosmic Dawn III (CoDaIII). CoDaIII's scale ($\rm 94^3 \, cMpc^3$) a…
▽ More
Recent determinations of the mean free path of ionising photons (mfp) in the intergalactic medium (IGM) at $\rm z=6$ are lower than many theoretical predictions. To gain insight into this issue, we investigate the evolution of the mfp in our new massive fully coupled radiation hydrodynamics cosmological simulation of reionization: Cosmic Dawn III (CoDaIII). CoDaIII's scale ($\rm 94^3 \, cMpc^3$) and resolution ($\rm 8192^3$ grid) make it particularly suitable to study the evolution of the IGM during Reionization. The simulation was performed with RAMSES-CUDATON on Summit, and used 131072 processors coupled to 24576 GPUs, making it the largest Reionization simulation, and largest RAMSES simulation ever performed. A superior agreement with global constraints on Reionization is obtained in CoDaIII over CoDaII especially for the evolution of the neutral hydrogen fraction and the cosmic photo-ionization rate, thanks to an improved calibration, later end of reionization ($\rm z=5.6$), and higher spatial resolution. Analyzing the mfp, we find that CoDaIII reproduces the most recent observations very well, from $\rm z=6$ to $\rm z=4.6$. We show that the distribution of the mfp in CoDaIII is bimodal, with short (neutral) and long (ionized) mfp modes, respectively, due to the patchiness of reionization and the co-existence of neutral versus ionized regions during Reionization. The neutral mode peaks at sub-kpc to kpc scales of mfp, while the ionized mode peak evolves from $\rm 0.1 Mpc/h$ at $\rm z=7$ to $\sim 10$ Mpc/h at $\rm z=5.2$. Computing the mfp as the average of the ionized mode provides the best match to the recent observational determinations. The distribution reduces to a single neutral (ionized) mode at $\rm z>13$ ($\rm z<5$).
△ Less
Submitted 22 August, 2022; v1 submitted 11 February, 2022;
originally announced February 2022.