-
The impact of varying inhomogeneous reionization histories on metrics of Ly$α$ opacity
Authors:
Caitlin C. Doughty,
Joseph F. Hennawi,
Jose Oñorbe,
Frederick B. Davies,
Zarija Lukić
Abstract:
The epoch of hydrogen reionization is complete by $z=5$, but its progression at higher redshifts is uncertain. Measurements of Ly$α$ forest opacity show large scatter at $z<6$, suggestive of spatial fluctuations in neutral fraction ($x_\mathrm{HI}$), temperature, or ionizing background, either individually or in combination. However, these effects are degenerate, necessitating modeling these physi…
▽ More
The epoch of hydrogen reionization is complete by $z=5$, but its progression at higher redshifts is uncertain. Measurements of Ly$α$ forest opacity show large scatter at $z<6$, suggestive of spatial fluctuations in neutral fraction ($x_\mathrm{HI}$), temperature, or ionizing background, either individually or in combination. However, these effects are degenerate, necessitating modeling these physics in tandem in order to properly interpret the observations. We begin this process by developing a framework for modeling the reionization history and associated temperature fluctuations, with the intention of incorporating ionizing background fluctuations at a later time. To do this, we generate several reionization histories using semi-numerical code AMBER, selecting histories with volume-weighted neutral fractions that adhere to the observed CMB optical depth and dark pixel fractions. Implementing these histories in the \texttt{Nyx} cosmological hydrodynamics code, we examine the evolution of gas within the simulation, and the associated metrics of the Ly$α$ forest opacity. We find that the pressure smoothing scale within the IGM is strongly correlated with the adiabatic index of the temperature-density relation. We find that while models with 20,000 K photoheating at reionization are better able to reproduce the shape of the observed $z=5$ 1D flux power spectrum than those with 10,000 K, they fail to match the highest wavenumbers. The simulated autocorrelation function and optical depth distributions are systematically low and narrow, respectively, compared to the observed values, but are in better agreement when the reionization history is longer in duration, more symmetric in its distribution of reionization redshifts, or if there are remaining neutral regions at $z<6$. The systematically low variance likely requires the addition of a fluctuating UVB.
△ Less
Submitted 3 December, 2023;
originally announced December 2023.
-
Convergence of small scale Ly$α$ structure at high-$z$ under different reionization scenarios
Authors:
Caitlin C. Doughty,
Joseph F. Hennawi,
Frederick B. Davies,
Zarija Lukić,
Jose Oñorbe
Abstract:
The Ly$α$ forest (LAF) at $z>5$ probes the thermal and reionization history of the intergalactic medium (IGM) and the nature of dark matter, but its interpretation requires comparison to cosmological hydrodynamical simulations. At high-$z$, convergence of these simulations is more exacting since transmission is dominated by underdense voids that are challenging to resolve. With evidence mounting f…
▽ More
The Ly$α$ forest (LAF) at $z>5$ probes the thermal and reionization history of the intergalactic medium (IGM) and the nature of dark matter, but its interpretation requires comparison to cosmological hydrodynamical simulations. At high-$z$, convergence of these simulations is more exacting since transmission is dominated by underdense voids that are challenging to resolve. With evidence mounting for a late end to reionization, small structures down to the sub-kpc level may survive to later times than conventionally thought due to the reduced time for pressure smoothing to impact the gas, further tightening simulation resolution requirements. We perform a suite of simulations using the Eulerian cosmological hydrodynamics code Nyx, spanning domain sizes of 1.25-10 $h^{-1}$ Mpc and 5-80 $h^{-1}$ kpc cells, and explore the interaction of these variables with the timing of reionization on the properties of the matter distribution and the simulated LAF at $z=5.5$. In observable Ly$α$ power, convergence within 10% is achieved for $k< 0.1$ s/km, but larger $k$ shows deviation of up to 20 percent. While a later reionization retains more small structure in the density field, because of the greater thermal broadening there is little difference in the convergence of LAF power between early ($z=9$) and later ($z=6$) reionizations. We conclude that at $z\sim5.5$, resolutions of 10 kpc are necessary for convergence of LAF power at $k<0.1$ s/km, while higher-$k$ modes require higher resolution, and that the timing of reionization does not significantly impact convergence given realistic photoheating.
△ Less
Submitted 25 May, 2023;
originally announced May 2023.
-
The environments and hosts of metal absorption at z>5
Authors:
Caitlin C. Doughty,
Kristian M. Finlator
Abstract:
A growing population of metal absorbers are observed at z>5, many showing strong evolution in incidence approaching the epoch of hydrogen reionization. Follow-up surveys examining fields around these metals have resulted in galaxy detections but the direct physical relationship between the detected galaxies and absorbers is unclear. Upcoming observations will illuminate this galaxy-absorber relati…
▽ More
A growing population of metal absorbers are observed at z>5, many showing strong evolution in incidence approaching the epoch of hydrogen reionization. Follow-up surveys examining fields around these metals have resulted in galaxy detections but the direct physical relationship between the detected galaxies and absorbers is unclear. Upcoming observations will illuminate this galaxy-absorber relationship, but the theoretical framework for interpreting these observations is lacking. To inform future z>5 studies, we define the expected relationship between metals and galaxies using the Technicolor Dawn simulation to model metal absorption from z=5-7, encompassing the end of reionization. We find that metal absorber types and strengths are slightly better associated with their environment than with the traits of their host galaxies, as absorption system strengths are more strongly correlated with the local galaxy overdensity than the stellar mass of their host galaxy. For redshifts prior to the end of the epoch of reionization, strong high ionization transitions like C IV are more spatially correlated with brighter galaxies on scales of a few hundred proper kpc than are low ionization systems, due to the former's preference for environments with higher UVB amplitudes and those ions' relative rarity at z>6. Post-reionization, the galaxy counts near these high-ionization ions are reduced, and increase surrounding certain low-ionization ions due to a combination of their relative abundances and preferred denser gas phase. We conclude that galaxy-absorber relationships are expected to evolve rapidly such that high-ionization absorbers are better tracers of galaxies pre-reionization while low-ionization absorbers are better post-reionization.
△ Less
Submitted 14 November, 2022;
originally announced November 2022.
-
The effects of binary stars on galaxies and metal-enriched gas during reionization
Authors:
Caitlin C. Doughty,
Kristian M. Finlator
Abstract:
Binary stars are abundant in nearby galaxies, but are typically unaccounted for in simulations of the high redshift Universe. Stellar population synthesis models that include the effects of binary evolution result in greater relative abundances of ionizing photons that could significantly affect the ambient ionizing background during the epoch of hydrogen reionization, additionally leading to diff…
▽ More
Binary stars are abundant in nearby galaxies, but are typically unaccounted for in simulations of the high redshift Universe. Stellar population synthesis models that include the effects of binary evolution result in greater relative abundances of ionizing photons that could significantly affect the ambient ionizing background during the epoch of hydrogen reionization, additionally leading to differences in galaxy gas content and star formation. We use hydrodynamic cosmological simulations including in situ multifrequency radiative transfer to evaluate the effects of a high binary fraction in reionization-era galaxies on traits of the early intergalactic medium and the abundance of H I and He II ionizing photons. We further extend this to analyze the traits of enriched gas. In comparing metrics generated using a fiducial simulation assuming single stars with one incorporating a high binary fraction, we find that binary stars cause H I reionization to complete earlier and at an accelerated pace, while also increasing the abundances of high-ionization metals (C IV and Si IV) in simulated absorption spectra while reducing the abundance of low-ionization states (O I, Si II, and C II). However, through increased photoheating of galactic and circumgalactic gas, they simultaneously reduce the rate of star formation in low-mass galaxies, slowing the ongoing process of enrichment and suppressing their own ionizing background. This potentially contributes to a slower He II reionization process at $z\geq5$, and further indicates that self-regulation of galaxies could be underestimated when neglecting binary stellar evolution.
△ Less
Submitted 20 May, 2021;
originally announced May 2021.