-
Symbolic Implementation of Extensions of the $\texttt{PyCosmo}$ Boltzmann Solver
Authors:
Beatrice Moser,
Christiane S. Lorenz,
Uwe Schmitt,
Alexandre Refregier,
Janis Fluri,
Raphael Sgier,
Federica Tarsitano,
Lavinia Heisenberg
Abstract:
$\texttt{PyCosmo}$ is a Python-based framework for the fast computation of cosmological model predictions. One of its core features is the symbolic representation of the Einstein-Boltzmann system of equations. Efficient $\texttt{C/C++}$ code is generated from the $\texttt{SymPy}$ symbolic expressions making use of the $\texttt{sympy2c}…
▽ More
$\texttt{PyCosmo}$ is a Python-based framework for the fast computation of cosmological model predictions. One of its core features is the symbolic representation of the Einstein-Boltzmann system of equations. Efficient $\texttt{C/C++}$ code is generated from the $\texttt{SymPy}$ symbolic expressions making use of the $\texttt{sympy2c}$ package. This enables easy extensions of the equation system for the implementation of new cosmological models. We illustrate this with three extensions of the $\texttt{PyCosmo}$ Boltzmann solver to include a dark energy component with a constant equation of state, massive neutrinos and a radiation streaming approximation. We describe the $\texttt{PyCosmo}$ framework, highlighting new features, and the symbolic implementation of the new models. We compare the $\texttt{PyCosmo}$ predictions for the $Λ$CDM model extensions with $\texttt{CLASS}$, both in terms of accuracy and computational speed. We find a good agreement, to better than 0.1% when using high-precision settings and a comparable computational speed. Links to the Python Package Index (PyPI) page of the code release and to the PyCosmo Hub, an online platform where the package is installed, are available at: https://cosmology.ethz.ch/research/software-lab/PyCosmo.html.
△ Less
Submitted 17 June, 2022; v1 submitted 15 December, 2021;
originally announced December 2021.
-
Rapid Simulations of Halo and Subhalo Clustering
Authors:
Pascale Berner,
Alexandre Refregier,
Raphael Sgier,
Tomasz Kacprzak,
Luca Tortorelli,
Pierluigi Monaco
Abstract:
The analysis of cosmological galaxy surveys requires realistic simulations for their interpretation. Forward modelling is a powerful method to simulate galaxy clustering without the need for an underlying complex model. This approach requires fast cosmological simulations with a high resolution and large volume, to resolve small dark matter halos associated to single galaxies. In this work, we pre…
▽ More
The analysis of cosmological galaxy surveys requires realistic simulations for their interpretation. Forward modelling is a powerful method to simulate galaxy clustering without the need for an underlying complex model. This approach requires fast cosmological simulations with a high resolution and large volume, to resolve small dark matter halos associated to single galaxies. In this work, we present fast halo and subhalo clustering simulations based on the Lagrangian perturbation theory code PINOCCHIO, which generates halos and merger trees. The subhalo progenitors are extracted from the merger history and the survival of subhalos is modelled. We introduce a new fitting function for the subhalo merger time, which includes a redshift dependence of the fitting parameters. The spatial distribution of subhalos within their hosts is modelled using a number density profile. We compare our simulations with the halo finder ROCKSTAR applied to the full N-body code GADGET-2. The subhalo velocity function and the correlation function of halos and subhalos are in good agreement. We investigate the effect of the chosen number density profile on the resulting subhalo clustering. Our simulation is approximate yet realistic and significantly faster compared to a full N-body simulation combined with a halo finder. The fast halo and subhalo clustering simulations offer good prospects for galaxy forward models using subhalo abundance matching.
△ Less
Submitted 24 October, 2022; v1 submitted 15 December, 2021;
originally announced December 2021.
-
Dark Energy Survey Year 3 results: Cosmology with peaks using an emulator approach
Authors:
D. Zürcher,
J. Fluri,
R. Sgier,
T. Kacprzak,
M. Gatti,
C. Doux,
L. Whiteway,
A. Refregier,
C. Chang,
N. Jeffrey,
B. Jain,
P. Lemos,
D. Bacon,
A. Alarcon,
A. Amon,
K. Bechtol,
M. Becker,
G. Bernstein,
A. Campos,
R. Chen,
A. Choi,
C. Davis,
J. Derose,
S. Dodelson,
F. Elsner
, et al. (97 additional authors not shown)
Abstract:
We constrain the matter density $Ω_{\mathrm{m}}$ and the amplitude of density fluctuations $σ_8$ within the $Λ$CDM cosmological model with shear peak statistics and angular convergence power spectra using mass maps constructed from the first three years of data of the Dark Energy Survey (DES Y3). We use tomographic shear peak statistics, including cross-peaks: peak counts calculated on maps create…
▽ More
We constrain the matter density $Ω_{\mathrm{m}}$ and the amplitude of density fluctuations $σ_8$ within the $Λ$CDM cosmological model with shear peak statistics and angular convergence power spectra using mass maps constructed from the first three years of data of the Dark Energy Survey (DES Y3). We use tomographic shear peak statistics, including cross-peaks: peak counts calculated on maps created by taking a harmonic space product of the convergence of two tomographic redshift bins. Our analysis follows a forward-modelling scheme to create a likelihood of these statistics using N-body simulations, using a Gaussian process emulator. We include the following lensing systematics: multiplicative shear bias, photometric redshift uncertainty, and galaxy intrinsic alignment. Stringent scale cuts are applied to avoid biases from unmodelled baryonic physics. We find that the additional non-Gaussian information leads to a tightening of the constraints on the structure growth parameter yielding $S_8~\equiv~σ_8\sqrt{Ω_{\mathrm{m}}/0.3}~=~0.797_{-0.013}^{+0.015}$ (68% confidence limits), with a precision of 1.8%, an improvement of ~38% compared to the angular power spectra only case. The results obtained with the angular power spectra and peak counts are found to be in agreement with each other and no significant difference in $S_8$ is recorded. We find a mild tension of $1.5 \thinspace σ$ between our study and the results from Planck 2018, with our analysis yielding a lower $S_8$. Furthermore, we observe that the combination of angular power spectra and tomographic peak counts breaks the degeneracy between galaxy intrinsic alignment $A_{\mathrm{IA}}$ and $S_8$, improving cosmological constraints. We run a suite of tests concluding that our results are robust and consistent with the results from other studies using DES Y3 data.
△ Less
Submitted 21 October, 2021; v1 submitted 19 October, 2021;
originally announced October 2021.
-
Combined $13\times2$-point analysis of the Cosmic Microwave Background and Large-Scale Structure: implications for the $S_8$-tension and neutrino mass constraints
Authors:
Raphael Sgier,
Christiane Lorenz,
Alexandre Refregier,
Janis Fluri,
Dominik Zürcher,
Federica Tarsitano
Abstract:
We present cosmological constraints for the flat $Λ$CDM model, including the sum of neutrino masses, by performing a multi-probe analysis of a total of 13 tomographic auto- and cross-angular power spectra. This is achieved by combining, at map level, the latest primary CMB and CMB-lensing measurements from the Planck 2018 data release, as well as spectroscopic galaxy samples from BOSS DR12, and th…
▽ More
We present cosmological constraints for the flat $Λ$CDM model, including the sum of neutrino masses, by performing a multi-probe analysis of a total of 13 tomographic auto- and cross-angular power spectra. This is achieved by combining, at map level, the latest primary CMB and CMB-lensing measurements from the Planck 2018 data release, as well as spectroscopic galaxy samples from BOSS DR12, and the latest Kilo-Degree Survey (KiDS-1000) tomographic weak lensing shear data release. Our analysis includes auto- and cross-correlations as well as calibration parameters for all cosmological probes, thus providing a self-calibration of the combined data sets. We find a good fit (reduced $χ^2$=1.7) for the combined probes with calibration parameters only moderately different from their nominal value, thus giving a possible interpretation of the tension between the early- and late-Universe probes. The resulting value for the structure growth parameter is $S_8 = 0.754 \pm 0.016$ (68\% CL). We also obtain a $\sim$2.3$σ$ constraint on the neutrino mass sum of $\sum m_ν= 0.51^{+0.21}_{-0.24}$ eV (68\% CL), which is compatible with current particle physics limits. We perform several tests by fixing the neutrino mass sum to a low value, considering narrower priors on the multiplicative bias parameters for cosmic shear, and by fixing all calibration parameters to their expected values. These tests result in worse fits compared to our fiducial run, especially for the case when all calibration parameters are fixed. This latter test also yields a lower upper limit of the neutrino mass sum. We discuss how the interplay between the cosmological and calibration parameters impact the $S_8$-tension and the constraints on the neutrino mass sum. [abridged]
△ Less
Submitted 7 October, 2021;
originally announced October 2021.
-
Baryonic effects in the Effective Field Theory of Large-Scale Structure and an analytic recipe for lensing in CMB-S4
Authors:
Diogo P. L. Bragança,
Matthew Lewandowski,
David Sekera,
Leonardo Senatore,
Raphael Sgier
Abstract:
Upcoming Large-Scale Structure surveys will likely become the next leading sources of cosmological information, making it crucial to have a precise understanding of the influence of baryons on cosmological observables. The Effective Field Theory of Large-Scale Structure (EFTofLSS) provides a consistent way to predict the clustering of dark matter and baryons on large scales, where their leading co…
▽ More
Upcoming Large-Scale Structure surveys will likely become the next leading sources of cosmological information, making it crucial to have a precise understanding of the influence of baryons on cosmological observables. The Effective Field Theory of Large-Scale Structure (EFTofLSS) provides a consistent way to predict the clustering of dark matter and baryons on large scales, where their leading corrections in perturbation theory are given by a simple and calculable functional form even after the onset of baryonic processes. In this paper, we extend the two-fluid-like system up to two-loop order in perturbation theory. Along the way, we show that a new linear counterterm proportional to the relative velocity of the fluids could generically be present, but we show that its effects are expected to be small in our universe. Regardless, we show how to consistently perform perturbation theory in the presence of this new term. We find that the EFTofLSS at two-loop order can accurately account for the details of baryonic processes on large scales. We compare our results to a hydrodynamical $N$-body simulation at many redshifts and find that the counterterms associated with the leading corrections to dark matter and baryons start to differ between redshifts $z \approx 3$ and $z \approx 2$, signaling the onset of star-formation physics. We then use these fits to compute the lensing power spectrum, show that the understanding of baryonic processes will be important for analyzing CMB-S4 data, and show that the two-loop EFTofLSS accurately captures these effects for $\ell \lesssim 2000$. Our results are also potentially of interest for current and future weak lensing surveys.
△ Less
Submitted 6 October, 2020;
originally announced October 2020.
-
Fast Lightcones for Combined Cosmological Probes
Authors:
Raphael Sgier,
Janis Fluri,
Jörg Herbel,
Alexandre Réfrégier,
Adam Amara,
Tomasz Kacprzak,
Andrina Nicola
Abstract:
The combination of different cosmological probes offers stringent tests of the $Λ$CDM model and enhanced control of systematics. For this purpose, we present an extension of the lightcone generator UFalcon first introduced in Sgier et al. 2019 (arXiv:1801.05745), enabling the simulation of a self-consistent set of maps for different cosmological probes. Each realization is generated from the same…
▽ More
The combination of different cosmological probes offers stringent tests of the $Λ$CDM model and enhanced control of systematics. For this purpose, we present an extension of the lightcone generator UFalcon first introduced in Sgier et al. 2019 (arXiv:1801.05745), enabling the simulation of a self-consistent set of maps for different cosmological probes. Each realization is generated from the same underlying simulated density field, and contains full-sky maps of different probes, namely weak lensing shear, galaxy overdensity including RSD, CMB lensing, and CMB temperature anisotropies from the ISW effect. The lightcone generation performed by UFalcon is parallelized and based on the replication of a large periodic volume simulated with the GPU-accelerated $N$-Body code PkdGrav3. The post-processing to construct the lightcones requires only a runtime of about 1 walltime-hour corresponding to about 100 CPU-hours. We use a randomization procedure to increase the number of quasi-independent full-sky UFalcon map-realizations, which enables us to compute an accurate multi-probe covariance matrix. Using this framework, we forecast cosmological parameter constraints by performing a multi-probe likelihood analysis for a combination of simulated future stage-IV-like surveys. We find that the inclusion of the cross-correlations between the probes significantly increases the information gain in the parameter constraints. We also find that the use of a non-Gaussian covariance matrix is increasingly important, as more probes and cross-correlation power spectra are included. A version of the UFalcon package currently including weak gravitational lensing is publicly available.
△ Less
Submitted 14 May, 2021; v1 submitted 11 July, 2020;
originally announced July 2020.
-
Cosmological Forecast for non-Gaussian Statistics in large-scale weak Lensing Surveys
Authors:
Dominik Zürcher,
Janis Fluri,
Raphael Sgier,
Tomasz Kacprzak,
Alexandre Refregier
Abstract:
Cosmic shear data contains a large amount of cosmological information encapsulated in the non-Gaussian features of the weak lensing mass maps. This information can be extracted using non-Gaussian statistics. We compare the constraining power in the $Ω_{\mathrm{m}} - σ_8$ plane of three map-based non-Gaussian statistics with the angular power spectrum, namely; peak/minimum counts and Minkowski func…
▽ More
Cosmic shear data contains a large amount of cosmological information encapsulated in the non-Gaussian features of the weak lensing mass maps. This information can be extracted using non-Gaussian statistics. We compare the constraining power in the $Ω_{\mathrm{m}} - σ_8$ plane of three map-based non-Gaussian statistics with the angular power spectrum, namely; peak/minimum counts and Minkowski functionals. We further analyze the impact of tomography and systematic effects originating from galaxy intrinsic alignments, multiplicative shear bias and photometric redshift systematics. We forecast the performance of the statistics for a stage-3-like weak lensing survey and restrict ourselves to scales $\geq$ 10 arcmin. We find, that in our setup, the considered non-Gaussian statistics provide tighter constraints than the angular power spectrum. The peak counts show the greatest potential, increasing the Figure-of-Merit (FoM) in the $Ω_{\mathrm{m}} - σ_8$ plane by a factor of about 4. A combined analysis using all non-Gaussian statistics in addition to the power spectrum increases the FoM by a factor of 5 and reduces the error on $S_8$ by $\approx$ 25\%. We find that the importance of tomography is diminished when combining non-Gaussian statistics with the angular power spectrum. The non-Gaussian statistics indeed profit less from tomography and the minimum counts and Minkowski functionals add some robustness against galaxy intrinsic alignment in a non-tomographic setting. We further find that a combination of the angular power spectrum and the non-Gaussian statistics allows us to apply conservative scale cuts in the analysis, thus helping to minimize the impact of baryonic and relativistic effects, while conserving the cosmological constraining power. We make the code that was used to conduct this analysis publicly available.
△ Less
Submitted 22 June, 2020;
originally announced June 2020.
-
Predicting Cosmological Observables with PyCosmo
Authors:
F. Tarsitano,
U. Schmitt,
A. Refregier,
J. Fluri,
R. Sgier,
A. Nicola,
J. Herbel,
A. Amara,
T. Kacprzak,
L. Heisenberg
Abstract:
Current and upcoming cosmological experiments open a new era of precision cosmology, thus demanding accurate theoretical predictions for cosmological observables. Because of the complexity of the codes delivering such predictions, reaching a high level of numerical accuracy is challenging. Among the codes already fulfilling this task, $\textsf{PyCosmo}$ is a Python based framework providing soluti…
▽ More
Current and upcoming cosmological experiments open a new era of precision cosmology, thus demanding accurate theoretical predictions for cosmological observables. Because of the complexity of the codes delivering such predictions, reaching a high level of numerical accuracy is challenging. Among the codes already fulfilling this task, $\textsf{PyCosmo}$ is a Python based framework providing solutions to the Einstein-Boltzmann equations and accurate predictions for cosmological observables. In this work, we first describe how the observables are implemented. Then, we check the accuracy of the theoretical predictions for background quantities, power spectra and Limber and beyond-Limber angular power spectra by comparison with other codes: the Core Cosmology Library ($\texttt{CCL}$), $\texttt{CLASS}$, $\texttt{HMCode}$ and $\texttt{iCosmo}$. In our analysis we quantify the agreement of $\textsf{PyCosmo}$ with the other codes, for a range of cosmological models, monitored through a series of $\textit{unit tests}$. $\textsf{PyCosmo}$, conceived as a multi purpose cosmology calculation tool in $\texttt{Python}$, is designed to be interactive and user friendly. A current version of the code (without the Boltzmann Solver) is publicly available and can be used interactively on the platform $\textsf{PyCosmo Hub}$, all accessible from this link: https://cosmology.ethz.ch/research/software-lab/PyCosmo.html . On the hub the users can perform their own computations using $\texttt{Jupyter Notebooks}$ without the need of installing any software, access to the results presented in this work and benefit from tutorial notebooks illustrating the usage of the code. The link above also redirects to the code release and documentation.
△ Less
Submitted 1 May, 2020;
originally announced May 2020.
-
Monte Carlo Control Loops for cosmic shear cosmology with DES Year 1
Authors:
T. Kacprzak,
J. Herbel,
A. Nicola,
R. Sgier,
F. Tarsitano,
C. Bruderer,
A. Amara,
A. Refregier,
S. L. Bridle,
A. Drlica-Wagner,
D. Gruen,
W. G. Hartley,
B. Hoyle,
L. F. Secco,
J. Zuntz,
J. Annis,
S. Avila,
E. Bertin,
D. Brooks,
E. Buckley-Geer,
A. Carnero Rosell,
M. Carrasco Kind,
J. Carretero,
L. N. da Costa,
J. De Vicente
, et al. (30 additional authors not shown)
Abstract:
Weak lensing by large-scale structure is a powerful probe of cosmology and of the dark universe. This cosmic shear technique relies on the accurate measurement of the shapes and redshifts of background galaxies and requires precise control of systematic errors. The Monte Carlo Control Loops (MCCL) is a forward modelling method designed to tackle this problem. It relies on the Ultra Fast Image Gene…
▽ More
Weak lensing by large-scale structure is a powerful probe of cosmology and of the dark universe. This cosmic shear technique relies on the accurate measurement of the shapes and redshifts of background galaxies and requires precise control of systematic errors. The Monte Carlo Control Loops (MCCL) is a forward modelling method designed to tackle this problem. It relies on the Ultra Fast Image Generator (UFig) to produce simulated images tuned to match the target data statistically, followed by calibrations and tolerance loops. We present the first end-to-end application of this method, on the Dark Energy Survey (DES) Year 1 wide field imaging data. We simultaneously measure the shear power spectrum $C_{\ell}$ and the redshift distribution $n(z)$ of the background galaxy sample. The method includes maps of the systematic sources, Point Spread Function (PSF), an Approximate Bayesian Computation (ABC) inference of the simulation model parameters, a shear calibration scheme, and the fast estimation of the covariance matrix. We find a close statistical agreement between the simulations and the DES Y1 data using an array of diagnostics. In a non-tomographic setting, we derive a set of $C_\ell$ and $n(z)$ curves that encode the cosmic shear measurement, as well as the systematic uncertainty. Following a blinding scheme, we measure the combination of $Ω_m$, $σ_8$, and intrinsic alignment amplitude $A_{\rm{IA}}$, defined as $S_8D_{\rm{IA}} = σ_8(Ω_m/0.3)^{0.5}D_{\rm{IA}}$, where $D_{\rm{IA}}=1-0.11(A_{\rm{IA}}-1)$. We find $S_8D_{\rm{IA}}=0.895^{+0.054}_{-0.039}$, where systematics are at the level of roughly 60\% of the statistical errors. We discuss these results in the context of earlier cosmic shear analyses of the DES Y1 data. Our findings indicate that this method and its fast runtime offer good prospects for cosmic shear measurements with future wide-field surveys.
△ Less
Submitted 3 June, 2019;
originally announced June 2019.
-
Effects of baryons on weak lensing peak statistics
Authors:
Andreas J. Weiss,
Aurel Schneider,
Raphael Sgier,
Tomasz Kacprzak,
Adam Amara,
Alexandre Refregier
Abstract:
Upcoming weak-lensing surveys have the potential to become leading cosmological probes provided all systematic effects are under control. Recently, the ejection of gas due to feedback energy from active galactic nuclei (AGN) has been identified as major source of uncertainty, challenging the success of future weak-lensing probes in terms of cosmology. In this paper we investigate the effects of ba…
▽ More
Upcoming weak-lensing surveys have the potential to become leading cosmological probes provided all systematic effects are under control. Recently, the ejection of gas due to feedback energy from active galactic nuclei (AGN) has been identified as major source of uncertainty, challenging the success of future weak-lensing probes in terms of cosmology. In this paper we investigate the effects of baryons on the number of weak-lensing peaks in the convergence field. Our analysis is based on full-sky convergence maps constructed via light-cones from $N$-body simulations, and we rely on the baryonic correction model of Schneider et al. (2019) to model the baryonic effects on the density field. As a result we find that the baryonic effects strongly depend on the Gaussian smoothing applied to the convergence map. For a DES-like survey setup, a smoothing of $θ_k\gtrsim8$ arcmin is sufficient to keep the baryon signal below the expected statistical error. Smaller smoothing scales lead to a significant suppression of high peaks (with signal-to-noise above 2), while lower peaks are not affected. The situation is more severe for a Euclid-like setup, where a smoothing of $θ_k\gtrsim16$ arcmin is required to keep the baryonic suppression signal below the statistical error. Smaller smoothing scales require a full modelling of baryonic effects since both low and high peaks are strongly affected by baryonic feedback.
△ Less
Submitted 1 October, 2019; v1 submitted 28 May, 2019;
originally announced May 2019.
-
DeepSphere: Efficient spherical Convolutional Neural Network with HEALPix sampling for cosmological applications
Authors:
Nathanaël Perraudin,
Michaël Defferrard,
Tomasz Kacprzak,
Raphael Sgier
Abstract:
Convolutional Neural Networks (CNNs) are a cornerstone of the Deep Learning toolbox and have led to many breakthroughs in Artificial Intelligence. These networks have mostly been developed for regular Euclidean domains such as those supporting images, audio, or video. Because of their success, CNN-based methods are becoming increasingly popular in Cosmology. Cosmological data often comes as spheri…
▽ More
Convolutional Neural Networks (CNNs) are a cornerstone of the Deep Learning toolbox and have led to many breakthroughs in Artificial Intelligence. These networks have mostly been developed for regular Euclidean domains such as those supporting images, audio, or video. Because of their success, CNN-based methods are becoming increasingly popular in Cosmology. Cosmological data often comes as spherical maps, which make the use of the traditional CNNs more complicated. The commonly used pixelization scheme for spherical maps is the Hierarchical Equal Area isoLatitude Pixelisation (HEALPix). We present a spherical CNN for analysis of full and partial HEALPix maps, which we call DeepSphere. The spherical CNN is constructed by representing the sphere as a graph. Graphs are versatile data structures that can act as a discrete representation of a continuous manifold. Using the graph-based representation, we define many of the standard CNN operations, such as convolution and pooling. With filters restricted to being radial, our convolutions are equivariant to rotation on the sphere, and DeepSphere can be made invariant or equivariant to rotation. This way, DeepSphere is a special case of a graph CNN, tailored to the HEALPix sampling of the sphere. This approach is computationally more efficient than using spherical harmonics to perform convolutions. We demonstrate the method on a classification problem of weak lensing mass maps from two cosmological models and compare the performance of the CNN with that of two baseline classifiers. The results show that the performance of DeepSphere is always superior or equal to both of these baselines. For high noise levels and for data covering only a smaller fraction of the sphere, DeepSphere achieves typically 10% better classification accuracy than those baselines. Finally, we show how learned filters can be visualized to introspect the neural network.
△ Less
Submitted 26 March, 2019; v1 submitted 29 October, 2018;
originally announced October 2018.
-
Weak lensing peak statistics in the era of large scale cosmological surveys
Authors:
Janis Fluri,
Tomasz Kacprzak,
Raphael Sgier,
Alexandre Réfrégier,
Adam Amara
Abstract:
Weak lensing peak counts are a powerful statistical tool for constraining cosmological parameters. So far, this method has been applied only to surveys with relatively small areas, up to several hundred square degrees. As future surveys will provide weak lensing datasets with size of thousands of square degrees, the demand on the theoretical prediction of the peak statistics will become heightened…
▽ More
Weak lensing peak counts are a powerful statistical tool for constraining cosmological parameters. So far, this method has been applied only to surveys with relatively small areas, up to several hundred square degrees. As future surveys will provide weak lensing datasets with size of thousands of square degrees, the demand on the theoretical prediction of the peak statistics will become heightened. In particular, large simulations of increased cosmological volume are required. In this work, we investigate the possibility of using simulations generated with the fast Comoving-Lagrangian acceleration (COLA) method, coupled to the convergence map generator Ufalcon, for predicting the peak counts. We examine the systematics introduced by the COLA method by comparing it with a full TreePM code. We find that for a 2000 deg$^2$ survey, the systematic error is much smaller than the statistical error. This suggests that the COLA method is able to generate promising theoretical predictions for weak lensing peaks. We also examine the constraining power of various configurations of data vectors, exploring the influence of splitting the sample into tomographic bins and combining different smoothing scales. We find the combination of smoothing scales to have the most constraining power, improving the constraints on the $S_8$ amplitude parameter by at least 40% compared to a single smoothing scale, with tomography brining only limited increase in measurement precision.
△ Less
Submitted 31 October, 2018; v1 submitted 22 March, 2018;
originally announced March 2018.
-
Fast cosmic web simulations with generative adversarial networks
Authors:
Andres C. Rodriguez,
Tomasz Kacprzak,
Aurelien Lucchi,
Adam Amara,
Raphael Sgier,
Janis Fluri,
Thomas Hofmann,
Alexandre Réfrégier
Abstract:
Dark matter in the universe evolves through gravity to form a complex network of halos, filaments, sheets and voids, that is known as the cosmic web. Computational models of the underlying physical processes, such as classical N-body simulations, are extremely resource intensive, as they track the action of gravity in an expanding universe using billions of particles as tracers of the cosmic matte…
▽ More
Dark matter in the universe evolves through gravity to form a complex network of halos, filaments, sheets and voids, that is known as the cosmic web. Computational models of the underlying physical processes, such as classical N-body simulations, are extremely resource intensive, as they track the action of gravity in an expanding universe using billions of particles as tracers of the cosmic matter distribution. Therefore, upcoming cosmology experiments will face a computational bottleneck that may limit the exploitation of their full scientific potential. To address this challenge, we demonstrate the application of a machine learning technique called Generative Adversarial Networks (GAN) to learn models that can efficiently generate new, physically realistic realizations of the cosmic web. Our training set is a small, representative sample of 2D image snapshots from N-body simulations of size 500 and 100 Mpc. We show that the GAN-generated samples are qualitatively and quantitatively very similar to the originals. For the larger boxes of size 500 Mpc, it is very difficult to distinguish them visually. The agreement of the power spectrum $P_k$ is 1-2\% for most of the range, between $k=0.06$ and $k=0.4$. An important advantage of generating cosmic web realizations with a GAN is the considerable gains in terms of computation time. Each new sample generated by a GAN takes a fraction of a second, compared to the many hours needed by traditional N-body techniques. We anticipate that the use of generative models such as GANs will therefore play an important role in providing extremely fast and precise simulations of cosmic web in the era of large cosmological surveys, such as Euclid and Large Synoptic Survey Telescope (LSST).
△ Less
Submitted 29 November, 2018; v1 submitted 27 January, 2018;
originally announced January 2018.
-
Fast Generation of Covariance Matrices for Weak Lensing
Authors:
Raphael Sgier,
Alexandre Réfrégier,
Adam Amara,
Andrina Nicola
Abstract:
Upcoming weak lensing surveys will probe large fractions of the sky with unprecedented accuracy. To infer cosmological constraints, a large ensemble of survey simulations are required to accurately model cosmological observables and their covariances. We develop a parallelized multi-lens-plane pipeline called UFalcon, designed to generate full-sky weak lensing maps from lightcones within a minimal…
▽ More
Upcoming weak lensing surveys will probe large fractions of the sky with unprecedented accuracy. To infer cosmological constraints, a large ensemble of survey simulations are required to accurately model cosmological observables and their covariances. We develop a parallelized multi-lens-plane pipeline called UFalcon, designed to generate full-sky weak lensing maps from lightcones within a minimal runtime. It makes use of L-PICOLA, an approximate numerical code, which provides a fast and accurate alternative to cosmological $N$-Body simulations. The UFalcon maps are constructed by nesting 2 simulations covering a redshift-range from $z=0.1$ to $1.5$ without replicating the simulation volume. We compute the convergence and projected overdensity maps for L-PICOLA in the lightcone or snapshot mode. The generation of such a map, including the L-PICOLA simulation, takes about 3 hours walltime on 220 cores. We use the maps to calculate the spherical harmonic power spectra, which we compare to theoretical predictions and to UFalcon results generated using the full $N$-Body code GADGET-2. We then compute the covariance matrix of the full-sky spherical harmonic power spectra using 150 UFalcon maps based on L-PICOLA in lightcone mode. We consider the PDF, the higher-order moments and the variance of the smoothed field variance to quantify the accuracy of the covariance matrix, which we find to be a few percent for scales $\ell \sim 10^2$ to $10^3$. We test the impact of this level of accuracy on cosmological constraints using an optimistic survey configuration, and find that the final results are robust to this level of uncertainty. The speed and accuracy of our developed pipeline provides a basis to also include further important features such as masking, varying noise and will allow us to compute covariance matrices for models beyond $Λ$CDM. [abridged]
△ Less
Submitted 17 January, 2018;
originally announced January 2018.
-
Cosmological model discrimination with Deep Learning
Authors:
Jorit Schmelzle,
Aurelien Lucchi,
Tomasz Kacprzak,
Adam Amara,
Raphael Sgier,
Alexandre Réfrégier,
Thomas Hofmann
Abstract:
We demonstrate the potential of Deep Learning methods for measurements of cosmological parameters from density fields, focusing on the extraction of non-Gaussian information. We consider weak lensing mass maps as our dataset. We aim for our method to be able to distinguish between five models, which were chosen to lie along the $σ_8$ - $Ω_m$ degeneracy, and have nearly the same two-point statistic…
▽ More
We demonstrate the potential of Deep Learning methods for measurements of cosmological parameters from density fields, focusing on the extraction of non-Gaussian information. We consider weak lensing mass maps as our dataset. We aim for our method to be able to distinguish between five models, which were chosen to lie along the $σ_8$ - $Ω_m$ degeneracy, and have nearly the same two-point statistics. We design and implement a Deep Convolutional Neural Network (DCNN) which learns the relation between five cosmological models and the mass maps they generate. We develop a new training strategy which ensures the good performance of the network for high levels of noise. We compare the performance of this approach to commonly used non-Gaussian statistics, namely the skewness and kurtosis of the convergence maps. We find that our implementation of DCNN outperforms the skewness and kurtosis statistics, especially for high noise levels. The network maintains the mean discrimination efficiency greater than $85\%$ even for noise levels corresponding to ground based lensing observations, while the other statistics perform worse in this setting, achieving efficiency less than $70\%$. This demonstrates the ability of CNN-based methods to efficiently break the $σ_8$ - $Ω_m$ degeneracy with weak lensing mass maps alone. We discuss the potential of this method to be applied to the analysis of real weak lensing data and other datasets.
△ Less
Submitted 18 July, 2017; v1 submitted 17 July, 2017;
originally announced July 2017.