期刊
JOURNAL OF COSMOLOGY AND ASTROPARTICLE PHYSICS
卷 -, 期 2, 页码 -出版社
IOP Publishing Ltd
DOI: 10.1088/1475-7516/2019/02/031
关键词
cosmological parameters from LSS; cosmological simulations; Lyman alpha forest
资金
- National Science Foundation (NSF) [PHY-1607611]
- Science Research Council (VR) of Sweden
- European Research Council (ERC) under the European Community's Seventh Framework Programme (FP7/2007-2013)/ERC grant [306478-CosmicDawn]
- research project grant Fundamental Physics from Cosmological Surveys - Swedish Research Council (VR) [2017-04212]
- Royal Society
- Simons Foundation
- NSF [AST-1817256]
- European Union's Horizon 2020 research and innovation programme ERC (BePreSySe) [725327]
- Spanish MINECO AEI/FEDER, UE [AYA2014-58747-P]
- ICCUB (Unidad de Excelencia Maria de Maeztu) [MDM-2014-0369]
- Science and Technology Facilities Council (STFC) Ernest Rutherford Fellowship [ST/N003853/1]
- STFC Consolidated Grant [ST/R000476/1]
- University College London (UCL) Cosmoparticle Initiative
- STFC [ST/R000476/1, ST/N003853/1] Funding Source: UKRI
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.
作者
我是这篇论文的作者
点击您的名字以认领此论文并将其添加到您的个人资料中。
推荐
暂无数据