skip to main content

Title: Imprints of temperature fluctuations on the z ∼ 5 Lyman-α forest: a view from radiation-hydrodynamic simulations of reionization

Reionization leads to large spatial fluctuations in the intergalactic temperature that can persist well after its completion. We study the imprints of such fluctuations on the $z$ ∼ 5 Ly α forest flux power spectrum using a set of radiation-hydrodynamic simulations that model different reionization scenarios. We find that large-scale coherent temperature fluctuations bring ${\sim}20\text{--}60{{\ \rm per\ cent}}$ extra power at k ∼ 0.002 s km−1, with the largest enhancements in the models where reionization is extended or ends the latest. On smaller scales (k ≳ 0.1 s km−1), we find that temperature fluctuations suppress power by ${\lesssim}10{{\ \rm per\ cent}}$. We find that the shape of the power spectrum is mostly sensitive to the reionization mid-point rather than temperature fluctuations from reionization’s patchiness. However, for all of our models with reionization mid-points of $z$ ≤ 8 ($z$ ≤ 12), the shape differences are ${\lesssim}20{{\ \rm per\ cent}}$ (${\lesssim}40{{\ \rm per\ cent}}$) because of a surprisingly well-matched cancellation between thermal broadening and pressure smoothing that occurs for realistic thermal histories. We also consider fluctuations in the ultraviolet background, finding their impact on the power spectrum to be much smaller than temperature fluctuations at k ≳ 0.01 s km−1. Furthermore, we compare our models to power spectrum more » measurements, finding that none of our models with reionization mid-points of $z$ < 8 is strongly preferred over another and that all of our models with mid-points of $z$ ≥ 8 are excluded at 2.5σ. Future measurements may be able to distinguish between viable reionization models if they can be performed at lower k or, alternatively, if the error bars on the high-k power can be reduced by a factor of 1.5.

« less
 ;  ;  ;  ;  ;  ;  ;  
Award ID(s):
Publication Date:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Page Range or eLocation-ID:
p. 3177-3195
Oxford University Press
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT Using a set of high resolution simulations, we quantify the effect of species-specific initial transfer functions on probes of the intergalactic medium (IGM) via the Lyman-α forest. We focus on redshifts 2–6, after H i reionization. We explore the effect of these initial conditions on measures of the thermal state of the low density IGM: the curvature, Doppler width cutoff, and Doppler width distribution. We also examine the matter and flux power spectrum, and potential consequences for constraints on warm dark matter models. We find that the curvature statistic is at most affected at the $\approx 2{{\ \rm per\ cent}}$ level at z = 6. The Doppler width cutoff parameters are affected by $\approx 5{{\ \rm per\ cent}}$ for the intercept, and $\approx 8{{\ \rm per\ cent}}$ for the fit slope, though this is subdominant to sample variation. The Doppler width distribution shows a $\approx 30{{\ \rm per\ cent}}$ effect at z = 3, however the distribution is not fully converged with simulation box size and resolution. The flux power spectrum is at most affected by $\approx 5{{\ \rm per\ cent}}$ at high redshift and small scales. We discuss numerical convergence with simulation parameters.

    Feedback from active galactic nuclei and stellar processes changes the matter distribution on small scales, leading to significant systematic uncertainty in weak lensing constraints on cosmology. We investigate how the observable properties of group-scale haloes can constrain feedback’s impact on the matter distribution using Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS). Extending the results of previous work to smaller halo masses and higher wavenumber, k, we find that the baryon fraction in haloes contains significant information about the impact of feedback on the matter power spectrum. We explore how the thermal Sunyaev Zel’dovich (tSZ) signal from group-scale haloes contains similar information. Using recent Dark Energy Survey weak lensing and Atacama Cosmology Telescope tSZ cross-correlation measurements and models trained on CAMELS, we obtain 10 per cent constraints on feedback effects on the power spectrum at $k \sim 5\, h\, {\rm Mpc}^{-1}$. We show that with future surveys, it will be possible to constrain baryonic effects on the power spectrum to $\mathcal {O}(\lt 1~{{\ \rm per\ cent}})$ at $k = 1\, h\, {\rm Mpc}^{-1}$ and $\mathcal {O}(3~{{\ \rm per\ cent}})$ at $k = 5\, h\, {\rm Mpc}^{-1}$ using the methods that we introduce here. Finally, we investigate the impact of feedbackmore »on the matter bispectrum, finding that tSZ observables are highly informative in this case.

    « less
  3. Abstract

    Recently, the Hydrogen Epoch of Reionization Array (HERA) has produced the experiment’s first upper limits on the power spectrum of 21 cm fluctuations atz∼ 8 and 10. Here, we use several independent theoretical models to infer constraints on the intergalactic medium (IGM) and galaxies during the epoch of reionization from these limits. We find that the IGM must have been heated above the adiabatic-cooling threshold byz∼ 8, independent of uncertainties about IGM ionization and the radio background. Combining HERA limits with complementary observations constrains the spin temperature of thez∼ 8 neutral IGM to 27 KT¯S630 K (2.3 KT¯S640 K) at 68% (95%) confidence. They therefore also place a lower bound on X-ray heating, a previously unconstrained aspects of early galaxies. For example, if the cosmic microwave background dominates thez∼ 8 radio background, the new HERA limits imply that the first galaxies produced X-rays more efficiently than local ones. Thez∼ 10 limits require even earlier heating if dark-matter interactions cool the hydrogen gas. If an extra radio background is produced by galaxies, we rule out (at 95% confidence) the combination of high radio and low X-raymore »luminosities ofLr,ν/SFR > 4 × 1024W Hz−1M1yr andLX/SFR < 7.6 × 1039erg s−1M1yr. The new HERA upper limits neither support nor disfavor a cosmological interpretation of the recent Experiment to Detect the Global EOR Signature (EDGES) measurement. The framework described here provides a foundation for the interpretation of future HERA results.

    « less

    We investigate the contribution of extended radio sources such as Centaurus A, and Galactic supernova remnants (SNRs) to our ability to detect the statistical 21-cm signal from the Epoch of Reionisation (EoR) with the Murchison Widefield Array (MWA). These sources are typically ignored because they are in highly attenuated parts of the MWA primary beam, however, in aggregate, these sources have apparent flux densities of $10\, \rm {Jy}$ on angular scales we expect to detect the 21-cm signal. We create bespoke multicomponent 2D Gaussian models for Galactic SNRs and for Centaurus A, and simulate the visibilities for two MWA snapshot observations. We grid those visibilities and then Fourier transform them with respect to frequency, averaging them both spherically and cylindrically to produce the 1D and 2D power spectra. We compare the simulated 1D power spectra to the expected 21-$\rm {cm}$ power spectrum. We find that although these extended sources are in highly attenuated parts of the MWA primary beam pattern, collectively they have enough power (∼104−105 $\rm {mK^2}\, {\it h^{-3}} \, \rm {Mpc^{3}}$) on EoR significant modes $(|{\boldsymbol k}| \lesssim 0.1\, h\, \rm {Mpc^{-1}})$ to prohibit detection of the 21-$\rm {cm}$ signal (∼104 $\rm {mK^2}\, {\it h^{-3}} \, \rm {Mpc^{3}}$). Wemore »find that $50{-}90{{\ \rm per\ cent}}$ of sources must be removed in order to reduce leakage to a level of $\sim 10{-}20{{\ \rm per\ cent}}$ of the 21-$\rm {cm}$ power spectrum on EoR significant modes. The effects of wide-field extended sources will have implications on the detectability of the 21-$\rm {cm}$ signal for the MWA and with the future Square Kilometre Array (SKA).

    « less

    The variety of star formation histories (SFHs) of z ≳ 6 galaxies provides important insights into early star formation, but has been difficult to systematically quantify. Some observations suggest that many z ∼ 6–9 galaxies are dominated by ≳200 Myr stellar populations, implying significant star formation at z ≳ 9, while others find that most reionization era galaxies are ≲10 Myr, consistent with little z ≳ 9 star formation. Here, we quantify the distribution of ages of UV-bright ($-22.5\lesssim M_{\rm \small UV}\lesssim -21$) galaxies colour-selected to lie at z ≃ 6.6–6.9, an ideal redshift range to systematically study the SFHs of reionization era galaxies with ground-based observatories and Spitzer. We infer galaxy properties with two SED modelling codes and compare results, finding that stellar masses are largely insensitive to the model, but the inferred ages can vary by an order of magnitude. We infer a distribution of ages assuming a simple, parametric SFH model, finding a median age of ∼30–70 Myr depending on SED model. We quantify the fractions of ≤10 and ≥250 Myr galaxies, finding that these systems comprise ∼15–30 per cent and ∼20–25 per cent of the population, respectively. With a flexible SFH model, the shapes of the SFHs are consistent with those implied bymore »the simple model (e.g. young galaxies have rapidly rising SFHs). However, stellar masses can differ significantly, with those of young systems sometimes being more than an order of magnitude larger with the flexible SFH. We quantify the implications of these results for z ≳ 9 stellar mass assembly and discuss improvements expected from JWST.

    « less