| 1. |
- Bird, Simeon, et al.
(author)
-
An emulator for the Lyman-alpha forest
- 2019
-
In: Journal of Cosmology and Astroparticle Physics. - : IOP Publishing. - 1475-7516. ; :2
-
Journal article (peer-reviewed)abstract
- We present methods for interpolating between the 1-D flux power spectrum of the Lyman-alpha forest, as output by cosmological hydrodynamic simulations. Interpolation is necessary for cosmological parameter estimation due to the limited number of simulations possible. We construct an emulator for the Lyman-alpha forest flux power spectrum from 21 small simulations using Latin hypercube sampling and Gaussian process interpolation. We show that this emulator has a typical accuracy of 1 : 5% and a worst-case accuracy of 4%, which compares well to the current statistical error of 3-5% at z < 3 from BOSS DR9. We compare to the previous state of the art, quadratic polynomial interpolation. The Latin hypercube samples the entire volume of parameter space, while quadratic polynomial emulation samples only lower-dimensional subspaces. The Gaussian process provides an estimate of the emulation error and we show using test simulations that this estimate is reasonable. We construct a likelihood function and use it to show that the posterior constraints generated using the emulator are unbiased. We show that our Gaussian process emulator has lower emulation error than quadratic polynomial interpolation and thus produces tighter posterior confidence intervals, which will be essential for future Lyman-alpha surveys such as DESI.
|
|
| 2. |
- Pedersen, Christian, et al.
(author)
-
Massive neutrinos and degeneracies in Lyman-alpha forest simulations
- 2020
-
In: Journal of Cosmology and Astroparticle Physics. - : IOP Publishing. - 1475-7516. ; :4
-
Journal article (peer-reviewed)abstract
- Using a suite of hydrodynamical simulations with cold dark matter, baryons, and neutrinos, we present a detailed study of the effect of massive neutrinos on the 1-D and 3-D flux power spectra of the Lyman-alpha (Ly alpha) forest. The presence of massive neutrinos in cosmology induces a scale- and time-dependent suppression of structure formation that is strongest on small scales. Measuring this suppression is a key method for inferring neutrino masses from cosmological data, and is one of the main goals of ongoing and future surveys like eBOSS, DES, LSST, Euclid or DESI. The clustering in the Ly alpha forest traces the quasilinear power at late times and on small scales. In combination with observations of the cosmic microwave background, the forest therefore provides some of the tightest constraints on the sum of the neutrino masses. However there is a well-known degeneracy between Sigma m(v) and the amplitude of perturbations in the linear matter power spectrum. We study the corresponding degeneracy in the 1-D flux power spectrum of the Ly alpha forest, and for the first time also study this degeneracy in the 3-D flux power spectrum. We show that the non-linear effects of massive neutrinos on the Lycy forest, beyond the effect of linear power amplitude suppression, are negligible, and this degeneracy persists in the Ly alpha forest observables to a high precision. We discuss the implications of this degeneracy for choosing parametrisations of the Ly alpha forest for cosmological analysis.
|
|
| 3. |
- Rogers, Keir K., et al.
(author)
-
Bayesian emulator optimisation for cosmology : application to the Lyman-alpha forest
- 2019
-
In: Journal of Cosmology and Astroparticle Physics. - : IOP Publishing. - 1475-7516. ; :2
-
Journal article (peer-reviewed)abstract
- The Lyman-alpha forest provides strong constraints on both cosmological parameters and intergalactic medium astrophysics, which are forecast to improve further with the next generation of surveys including eBOSS and DESI. As is generic in cosmological inference, extracting this information requires a likelihood to be computed throughout a high-dimensional parameter space. Evaluating the likelihood requires a robust and accurate mapping between the parameters and observables, in this case the 1D flux power spectrum. Cosmological simulations enable such a mapping, but due to computational time constraints can only be evaluated at a handful of sample points; emulators are designed to interpolate between these. The problem then reduces to placing the sample points such that an accurate mapping is obtained while minimising the number of expensive simulations required. To address this, we introduce an emulation procedure that employs Bayesian optimisation of the training set for a Gaussian process interpolation scheme. Starting with a Latin hypercube sampling (other schemes with good space-filling properties can be used), we iteratively augment the training set with extra simulations at new parameter positions which balance the need to reduce interpolation error while focussing on regions of high likelihood. We show that smaller emulator error from the Bayesian optimisation propagates to smaller widths on the posterior distribution. Even with fewer simulations than a Latin hypercube, Bayesian optimisation shrinks the 95% credible volume by 90% and, e.g., the 1 sigma error on the amplitude of small-scale primordial fluctuations by 38%. This is the first demonstration of Bayesian optimisation applied to large-scale structure emulation, and we anticipate the technique will generalise to many other probes such as galaxy clustering, weak lensing and 21cm.
|
|
| 4. |
- Rogers, Keir K., et al.
(author)
-
Correlations in the three-dimensional Lyman-alpha forest contaminated by high column density absorbers
- 2018
-
In: Monthly notices of the Royal Astronomical Society. - : Oxford University Press (OUP). - 0035-8711 .- 1365-2966. ; 476:3, s. 3716-3728
-
Journal article (peer-reviewed)abstract
- Correlations measured in three dimensions in the Lyman-alpha forest are contaminated by the presence of the damping wings of high column density (HCD) absorbing systems of neutral hydrogen (H I; having column densities N(H I) > 1.6 x 10(17) atoms cm(-2)), which extend significantly beyond the redshift-space location of the absorber. We measure this effect as a function of the column density of the HCD absorbers and redshift by measuring threedimensional (3D) flux power spectra in cosmological hydrodynamical simulations from the Illustris project. Survey pipelines exclude regions containing the largest damping wings. We find that, even after this procedure, there is a scale-dependent correction to the 3D Lyman alpha forest flux power spectrum from residual contamination. We model this residual using a simple physical model of the HCD absorbers as linearly biased tracers of the matter density distribution, convolved with their Voigt profiles and integrated over the column density distribution function. We recommend the use of this model over existing models used in data analysis, which approximate the damping wings as top-hats and so miss shape information in the extended wings. The simple 'linear Voigt model' is statistically consistent with our simulation results for a mock residual contamination up to small scales (vertical bar k vertical bar < 1 h Mpc(-1)). It does not account for the effect of the highest column density absorbers on the smallest scales (e.g. vertical bar k vertical bar > 0.4 h Mpc(-1) for small damped Lyman-alpha absorbers; HCD absorbers with N(H-I) similar to 10(21) atoms cm(-2)). However, these systems are in any case preferentially removed from survey data. Our model is appropriate for an accurate analysis of the baryon acoustic oscillations feature. It is additionally essential for reconstructing the full shape of the 3D flux power spectrum.
|
|
| 5. |
- Rogers, Keir K., et al.
(author)
-
Simulating the effect of high column density absorbers on the one-dimensional Lyman alpha forest flux power spectrum
- 2018
-
In: Monthly notices of the Royal Astronomical Society. - : Oxford University Press (OUP). - 0035-8711 .- 1365-2966. ; 474:3, s. 3032-3042
-
Journal article (peer-reviewed)abstract
- We measure the effect of high column density absorbing systems of neutral hydrogen (H I) on the one-dimensional (1D) Lyman alpha forest flux power spectrum using cosmological hydrodynamical simulations from the Illustris project. High column density absorbers (which we define to be those with HI column densities N(H I) > 1.6 x 10(17) atoms cm(-2)) cause broadened absorption lines with characteristic damping wings. These damping wings bias the 1D Lyman alpha forest flux power spectrum by causing absorption in quasar spectra away from the location of the absorber itself. We investigate the effect of high column density absorbers on the Lyman alpha forest using hydrodynamical simulations for the first time. We provide templates as a function of column density and redshift, allowing the flexibility to accurately model residual contamination, i.e. if an analysis selectively clips out the largest damping wings. This flexibility will improve cosmological parameter estimation, for example, allowing more accurate measurement of the shape of the power spectrum, with implications for cosmological models containing massive neutrinos or a running of the spectral index. We provide fitting functions to reproduce these results so that they can be incorporated straightforwardly into a data analysis pipeline.
|
|
