skip to main content

This content will become publicly available on June 23, 2023

Title: LyMAS reloaded: improving the predictions of the large-scale Lyman- α forest statistics from dark matter density and velocity fields
ABSTRACT We present LyMAS2, an improved version of the ‘Lyman-α Mass Association Scheme’ aiming at predicting the large-scale 3D clustering statistics of the Lyman-α forest (Ly α) from moderate-resolution simulations of the dark matter (DM) distribution, with prior calibrations from high-resolution hydrodynamical simulations of smaller volumes. In this study, calibrations are derived from the Horizon-AGN suite simulations, (100 Mpc h)−3 comoving volume, using Wiener filtering, combining information from DM density and velocity fields (i.e. velocity dispersion, vorticity, line-of-sight 1D-divergence and 3D-divergence). All new predictions have been done at z = 2.5 in redshift space, while considering the spectral resolution of the SDSS-III BOSS Survey and different DM smoothing (0.3, 0.5, and 1.0 Mpc h−1 comoving). We have tried different combinations of DM fields and found that LyMAS2, applied to the Horizon-noAGN DM fields, significantly improves the predictions of the Ly α 3D clustering statistics, especially when the DM overdensity is associated with the velocity dispersion or the vorticity fields. Compared to the hydrodynamical simulation trends, the two-point correlation functions of pseudo-spectra generated with LyMAS2 can be recovered with relative differences of ∼5 per cent even for high angles, the flux 1D power spectrum (along the light of sight) with ∼2 per cent and the flux 1D probability distribution function exactly. more » Finally, we have produced several large mock BOSS spectra (1.0 and 1.5 Gpc h−1) expected to lead to much more reliable and accurate theoretical predictions. « less
; ; ; ; ; ; ; ;
Award ID(s):
Publication Date:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Page Range or eLocation-ID:
3222 to 3245
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract We present a new investigation of the intergalactic medium (IGM) near the end of reionization using “dark gaps” in the Ly α forest. Using spectra of 55 QSOs at z em > 5.5, including new data from the XQR-30 VLT Large Programme, we identify gaps in the Ly α forest where the transmission averaged over 1 comoving h −1 Mpc bins falls below 5%. Nine ultralong ( L > 80 h −1 Mpc) dark gaps are identified at z < 6. In addition, we quantify the fraction of QSO spectra exhibiting gaps longer than 30 h −1 Mpc, F 30 , as a function of redshift. We measure F 30 ≃ 0.9, 0.6, and 0.15 at z = 6.0, 5.8, and 5.6, respectively, with the last of these long dark gaps persisting down to z ≃5.3. Comparing our results with predictions from hydrodynamical simulations, we find that the data are consistent with models wherein reionization extends significantly below redshift six. Models wherein the IGM is essentially fully reionized that retain large-scale fluctuations in the ionizing UV background at z ≲6 are also potentially consistent with the data. Overall, our results suggest that signatures of reionization in the form ofmore »islands of neutral hydrogen and/or large-scale fluctuations in the ionizing background remain present in the IGM until at least z ≃ 5.3.« less
  2. ABSTRACT Hydrodynamic simulations provide a powerful, but computationally expensive, approach to study the interplay of dark matter and baryons in cosmological structure formation. Here, we introduce the EMulating Baryonic EnRichment (EMBER) Deep Learning framework to predict baryon fields based on dark matter-only simulations thereby reducing computational cost. EMBER comprises two network architectures, U-Net and Wasserstein Generative Adversarial Networks (WGANs), to predict 2D gas and H i densities from dark matter fields. We design the conditional WGANs as stochastic emulators, such that multiple target fields can be sampled from the same dark matter input. For training we combine cosmological volume and zoom-in hydrodynamical simulations from the Feedback in Realistic Environments (FIRE) project to represent a large range of scales. Our fiducial WGAN model reproduces the gas and H i power spectra within 10 per cent accuracy down to ∼10 kpc scales. Furthermore, we investigate the capability of EMBER to predict high resolution baryon fields from low resolution dark matter inputs through upsampling techniques. As a practical application, we use this methodology to emulate high-resolution H i maps for a dark matter simulation of a $L=100\, \text{Mpc}\, h^{ -1}$ comoving cosmological box. The gas content of dark matter haloes and the H i column density distributions predicted by EMBERmore »agree well with results of large volume cosmological simulations and abundance matching models. Our method provides a computationally efficient, stochastic emulator for augmenting dark matter only simulations with physically consistent maps of baryon fields.« less

    Ly α forest decomposed into Voigt profile components allows us to study clustering properties of the intergalactic-medium and its dependence on various physical quantities. Here, we report the first detections of probability excess of low-z (i.e z < 0.48) Ly α absorber triplets over redshift-space scale of r∥ ≤ 8 pMpc (Mpc in physical units) with maximum amplitude of $8.76^{+1.96}_{-1.65}$ at a longitudinal separation of 1–2 pMpc. We measure non-zero three-point correlation ($\zeta = 4.76^{+1.98}_{-1.67}$) only at this scale with reduced three-point correlation Q = $0.95^{+0.39}_{-0.38}$. The measured ζ shows an increasing trend with increasing minimum H i column density (NH i) threshold while Q does not show any NH i dependence. About 88 per cent of the triplets contributing to ζ (at z ≤ 0.2) have nearby galaxies (whose distribution is known to be complete for ∼0.1L* at z < 0.1 and for ∼L* at z ∼ 0.25 within 20 arcsec to the quasar sightlines) within velocity separation of 500 km s−1 and median impact parameter of 405 pkpc. The measured impact parameters are consistent with majority of the identified triplets not originating from individual galaxies but tracing the underlying galaxy distribution. Frequency of occurrence of Broad-Ly α absorbers (b > 40 km s−1) in triplets (∼85 per cent) ismore »factor ∼3 higher than that found among the full sample (∼32 per cent). Using four different cosmological simulations, we quantify the effect of peculiar velocities and feedback and show that most of the observed trends are broadly reproduced. However, ζ at small scales (r∥ < 1 pMpc) and its b-dependence found in simulations are inconsistent with observations. This could either be related to the failure of these simulations to reproduce the observed b and NH i distributions for NH i > 1014 cm−2 self-consistently or to the wide spread of signal-to-noise ratio in the observed data.

    « less
  4. Abstract

    We probe what constraints for the global ionized hydrogen fraction the void probability function (VPF) clustering can give for the Lyman Alpha Galaxies in the Epoch of Reionization (LAGER) narrowband survey as a function of area. Neutral hydrogen acts like a fog for Lyαemission, and measuring the drop in the luminosity function of Lyαemitters (LAEs) has been used to constrain the ionization fraction in narrowband surveys. However, the clustering of LAEs is independent of the luminosity function’s inherent evolution, and can offer additional constraints for reionization under different models. The VPF measures how likely a given circle is to be empty. It is a volume-averaged clustering statistic that traces the behavior of higher-order correlations, and its simplicity offers helpful frameworks for planning surveys. Using the Jensen et al. simulations of LAEs within various amounts of ionized intergalactic medium, we predict the behavior of the VPF in one (301 × 150.5 × 30 Mpc3), four (5.44 × 106Mpc3), or eight (1.1 × 107Mpc3) fields of LAGER imaging. We examine the VPF at 5′ and 13′, corresponding to the minimum scale implied by the LAE density and the separation of the 2D VPF from random, and the maximum scale from themore »eight-field 15.5 deg2LAGER area. We find that even a single DECam field of LAGER (2–3 deg2) could discriminate between mostly neutral versus ionized. Additionally, we find four fields allow for the distinction between 30%, 50%, and 95% ionized, and eight fields could even distinguish between 30%, 50%, 73%, and 95% ionized.

    « less

    We describe our non-linear emulation (i.e. interpolation) framework that combines the halo occupation distribution (HOD) galaxy bias model with N-body simulations of non-linear structure formation, designed to accurately predict the projected clustering and galaxy–galaxy lensing signals from luminous red galaxies in the redshift range 0.16 < z < 0.36 on comoving scales 0.6 < rp < 30 $h^{-1} \, \text{Mpc}$. The interpolation accuracy is ≲ 1–2 per cent across the entire physically plausible range of parameters for all scales considered. We correctly recover the true value of the cosmological parameter S8 = (σ8/0.8228)(Ωm/0.3107)0.6 from mock measurements produced via subhalo abundance matching (SHAM)-based light-cones designed to approximately match the properties of the SDSS LOWZ galaxy sample. Applying our model to Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 14 (DR14) LOWZ galaxy clustering and galaxy-shear cross-correlation measurements made with Sloan Digital Sky Survey (SDSS) Data Release 8 (DR8) imaging, we perform a prototype cosmological analysis marginalizing over wCDM cosmological parameters and galaxy HOD parameters. We obtain a 4.4 per cent measurement of S8 = 0.847 ± 0.037, in 3.5σ tension with the Planck cosmological results of 1.00 ± 0.02. We discuss the possibility of underestimated systematic uncertainties or astrophysical effects that could explain this discrepancy.