Likelihood-free Cosmological Constraints with Artificial Neural Networks: An Application on Hubble Parameters and SNe Ia

Yu-Chen Wang,Yuan-Bo Xie,Tong-Jie Zhang,Hui-Chao Huang,Tingting Zhang,Kun Liu
DOI: https://doi.org/10.3847/1538-4365/abf8aa
2021-04-15
Abstract:The errors of cosmological data generated from complex processes, such as the observational Hubble parameter data (OHD) and the Type Ia supernova (SN Ia) data, cannot be accurately modeled by simple analytical probability distributions, e.g. Gaussian distribution. To constrain cosmological parameters from these data, likelihood-free inference is usually used to bypass the direct calculation of the likelihood. In this paper, we propose a new procedure to perform likelihood-free cosmological inference using two artificial neural networks (ANN), the Masked Autoregressive Flow (MAF) and the denoising autoencoder (DAE). Our procedure is the first to use DAE to extract features from data, in order to simplify the structure of MAF needed to estimate the posterior. Tested on simulated Hubble parameter data with a simple Gaussian likelihood, the procedure shows the capability of extracting features from data and estimating posterior distributions without the need of tractable likelihood. We demonstrate that it can accurately approximate the real posterior, achieve performance comparable to the traditional MCMC method, and the MAF gets better training results for small number of simulation when the DAE is added. We also discuss the application of the proposed procedure to OHD and Pantheon SN Ia data, and use them to constrain cosmological parameters from the non-flat $\Lambda$CDM model. For SNe Ia, we use fitted light curve parameters to find constraints on $H_0,\Omega_m,\Omega_\Lambda$ similar to relevant work, using less empirical distributions. In addition, this work is also the first to use Gaussian process in the procedure of OHD simulation.
Cosmology and Nongalactic Astrophysics
What problem does this paper attempt to address?
The problem that this paper attempts to solve is how to constrain cosmological parameters without using traditional likelihood functions when dealing with complex astronomical data. Specifically, for the observed Hubble parameter data (OHD) and Type Ia supernovae (SN Ia) data, the errors of these data cannot be accurately modeled by simple analytical probability distributions (such as Gaussian distribution). Therefore, traditional likelihood - based methods are difficult to be directly applied to these data. To solve this problem, the author proposes a new method, using two artificial neural networks (ANN) - Masked Autoregressive Flow (MAF) and Denoising Autoencoder (DAE) to perform likelihood - free inference. ### Main contributions 1. **Propose a new likelihood - free inference process**: For the first time, the author uses Denoising Autoencoder (DAE) to extract features from data to simplify the MAF structure used for estimating posterior distributions. 2. **Verify the effectiveness of the method**: By testing this method on simulated Hubble parameter data, it shows its ability in feature extraction and posterior distribution estimation without relying on an analytical likelihood function. 3. **Comparison with traditional methods**: This method can accurately approximate the true posterior distribution, and its performance is comparable to that of the traditional MCMC method. Moreover, when the number of samples is small, adding DAE can improve the training effect of MAF. 4. **Practical application**: Discusses the application of this method in real - observed data (OHD and Pantheon SN Ia data) and uses it to constrain cosmological parameters in the non - flat ΛCDM model. 5. **Innovation points**: For the first time, Gaussian process is used in the OHD simulation process to better simulate the uncertainty of redshift variation, and inference is directly carried out from the data of SN Ia light - curve fitting, reducing the number of prior assumptions. ### Specific problems solved - **Error modeling of complex data**: Traditional likelihood functions are difficult to accurately describe the error distribution of complex data. - **Processing of high - dimensional data**: Through DAE, the data dimension is reduced, enabling MAF to be trained and inferred more efficiently. - **Optimization of computing resources**: Compared with the traditional MCMC method, this method is more computationally efficient, especially when the number of samples is small. In conclusion, by combining DAE and MAF, this paper provides a new likelihood - free inference method, which can effectively constrain cosmological parameters when dealing with complex astronomical data, and at the same time performs well in terms of computational efficiency and accuracy.