The CAMELS Project: Public Data Release

The Cosmology and Astrophysics with Machine Learning Simulations (CAMELS) project was developed to combine cosmology with astrophysics through thousands of cosmological hydrodynamic simulations and machine learning. CAMELS contains 4233 cosmological simulations, 2049 N-body simulations, and 2184 state-of-the-art hydrodynamic simulations that sample a vast volume in parameter space. In this paper, we present the CAMELS public data release, describing the characteristics of the CAMELS simulations and a variety of data products generated from them, including halo, subhalo, galaxy, and void catalogs, power spectra, bispectra, Lyα spectra, probability distribution functions, halo radial profiles, and X-rays photon lists. We also release over 1000 catalogs that contain billions of galaxies from CAMELS-SAM: a large collection of N-body simulations that have been combined with the Santa Cruz semianalytic model. We release all the data, comprising more than 350 terabytes and containing 143,922 snapshots, millions of halos, galaxies, and summary statistics. We provide further technical details on how to access, download, read, and process the data at https://camels.readthedocs.io

Non-linear perturbation theory extension of the Boltzmann code CLASS

We present a new open-source code that calculates one-loop power spectra and cross spectra for matter fields and biased tracers in real and redshift space. These spectra incorporate all ingredients required for a direct application to data: nonlinear bias and redshift-space distortions, infrared resummation, counterterms, and the Alcock-Paczynski effect. Our code is based on the Boltzmann solver class and inherits its advantageous properties: user friendliness, ease of modification, high speed, and simple interface with other software. We present detailed descriptions of the theoretical model, the code structure, approximations, and accuracy tests. A typical end-to-end run for one cosmology takes 0.3 seconds, which is sufficient for Markov chain Monte Carlo parameter extraction. As an example, we apply the code to the Baryon Oscillation Spectroscopic Survey (BOSS) data and infer cosmological parameters from the shape of the galaxy power spectrum.We present a new open-source code that calculates one-loop power auto- and cross-power spectra for matter fields and biased tracers in real and redshift space. These spectra incorporate all ingredients required for a direct application to data: non-linear bias, redshift-space distortions, infra-red resummation, counterterms, and the Alcock-Paczynski effect. Our code is based on the Boltzmann solver CLASS and inherits its advantage: user friendliness, ease of modification, high speed, and simple interface with other software. We present detailed descriptions of the theoretical model, the code structure, approximations, and accuracy tests. A typical end-to-end run for one cosmology takes $\sim 0.3$ seconds, which is sufficient for Markov Chain Monte Carlo parameter extraction. As an example, we apply the code to data from the Baryon Oscillation Spectroscopic Survey (BOSS) and infer cosmological parameters from the shape of the galaxy power spectrum

Constraints on Single-Field Inflation from the BOSS Galaxy Survey

Nonlocal primordial non-Gaussianity (NLPNG) is a smoking gun of interactions in single-field inflationary models and can be written as a combination of the equilateral and orthogonal templates. We present the first constraints on these from the redshift-space galaxy power spectra and bispectra of the BOSS data. These are the first such measurements independent of the cosmic microwave background fluctuations. We perform a consistent analysis that includes all necessary nonlinear corrections generated by NLPNG and vary all relevant cosmological and nuisance parameters in a global fit to the data. Our conservative analysis yields joint limits on the amplitudes of the equilateral and orthogonal shapes, fNLequil=940±600 and fNLortho=-170±170 (both at 68% CL). These can be used to derive constraints on coefficients of the effective single-field inflationary Lagrangian; in particular, we find that the sound speed of inflaton fluctuations has the bound cs≥0.013 at 95% CL. Fixing the quadratic galaxy bias and cosmological parameters, the constraints can be tightened to fNLequil=260±300 and fNLortho=-23±120 (68% CL).Non-local primordial non-Gaussianity (NLPNG) is a smoking gun of interactions in single-field inflationary models, and can be written as a combination of the equilateral and orthogonal templates. We present the first constraints on these from the redshift-space galaxy power spectra and bispectra of the Baryon Oscillation Spectroscopic Survey (BOSS) data. These are the first such measurements independent of the cosmic microwave background fluctuations. We perform a consistent analysis that includes all necessary nonlinear corrections generated by NLPNG, and vary all relevant cosmological and nuisance parameters in a global fit to the data. Our conservative analysis yields joint limits on the amplitudes of the equilateral and orthogonal shapes, $f_{\rm NL}^{\rm equil}=940\pm 600$, $f_{\rm NL}^{\rm ortho}= -170\pm 170$ (both at 68% CL). These can be used to derive constraints on coefficients of the effective single-field inflationary Lagrangian; in particular, we find that the sound speed of inflaton fluctuations has the bound $c_s\geq 0.013$ at 95% CL. Fixing the quadratic galaxy bias and cosmological parameters, the constraints can be tightened to $f_{\rm NL}^{\rm equil}=260\pm 300$, $f_{\rm NL}^{\rm ortho}= -23\pm 120$ (68% CL)

Cosmology with the Galaxy Bispectrum Multipoles: Optimal Estimation and Application to BOSS Data

We present a framework for self-consistent cosmological analyses of the full-shape anisotropic bispectrum, including the quadrupole (ℓ=2) and hexadecapole (ℓ=4) moments. This features a novel window-free algorithm for extracting the latter quantities from data, derived using a maximum-likelihood prescription. Furthermore, we introduce a theoretical model for the bispectrum multipoles (which does not introduce new free parameters), and test both aspects of the pipeline on several high-fidelity mocks, including the PT Challenge suite of gigantic cumulative volume. This establishes that the systematic error is significantly below the statistical threshold, both for the measurement and modeling. As a realistic example, we extract the large-scale bispectrum multipoles from BOSS DR12 and analyze them in combination with the power spectrum data. Assuming a minimal ΛCDM model, with a BBN prior on the baryon density and a Planck prior on ns, we can extract the remaining cosmological parameters directly from the clustering data. The inclusion of the unwindowed higher-order (ℓ>0) large-scale bispectrum multipoles is found to moderately improve one-dimensional cosmological parameter posteriors (at the 5%–10% level), though these multipoles are detected only in three out of four BOSS data segments at ≈5σ. Combining information from the power spectrum and bispectrum multipoles, the real space power spectrum, and the postreconstructed BAO data, we find H0=68.2±0.8  km s-1 Mpc-1, Ωm=0.33±0.01 and σ8=0.736±0.033 (the tightest yet found in perturbative full-shape analyses). Our estimate of the growth parameter S8=0.77±0.04 agrees with both weak lensing and CMB results. The estimators and data used in this work have been made publicly available.We present a framework for self-consistent cosmological analyses of the full-shape anisotropic bispectrum, including the quadrupole $(\ell=2)$ and hexadecapole $(\ell=4)$ moments. This features a novel window-free algorithm for extracting the latter quantities from data, derived using a maximum-likelihood prescription. Furthermore, we introduce a theoretical model for the bispectrum multipoles (which does not introduce new free parameters), and test both aspects of the pipeline on several high-fidelity mocks, including the PT Challenge suite of gigantic cumulative volume. This establishes that the systematic error is significantly below the statistical threshold, both for the measurement and modeling. As a realistic example, we extract the large-scale bispectrum multipoles from BOSS DR12 and analyze them in combination with the power spectrum data. Assuming a minimal $\Lambda$CDM model, with a BBN prior on the baryon density and a \textit{Planck} prior on $n_s$, we can extract the remaining cosmological parameters directly from the clustering data. The inclusion of the unwindowed higher-order $(\ell>0)$ large-scale bispectrum multipoles is found to moderately improve one-dimensional cosmological parameter posteriors (at the $5\%-10\%$ level), though these multipoles are detected only in three out of four BOSS data segments at $\approx 5\sigma$. Combining information from the power spectrum and bispectrum multipoles, the real space power spectrum, and the post-reconstructed BAO data, we find $H_0 = 68.2\pm 0.8~\mathrm{km}\,\mathrm{s}^{-1}\mathrm{Mpc}^{-1}$, $\Omega_m =0.33\pm 0.01$ and $\sigma_8 = 0.736\pm 0.033$ (the tightest yet found in perturbative full-shape analyses). Our estimate of the growth parameter $S_8=0.77\pm 0.04$ agrees with both weak lensing and CMB results

Cosmological constraints without nonlinear redshift-space distortions

Nonlinear redshift-space distortions (“fingers of God”) are challenging to model analytically, a fact that limits the applicability of perturbation theory (PT) in redshift space as compared to real space. We show how this problem can be mitigated using a new observable, Q0, which can be easily estimated from the redshift-space clustering data and is approximately equal to the real-space power spectrum. The new statistic does not suffer from fingers of God and can be accurately described with PT down to kmax≃0.4  h Mpc-1. It can be straightforwardly included in the likelihood at negligible additional computational cost and yields noticeable improvements on cosmological parameters compared to standard power spectrum multipole analyses. Using both simulations and observational data from the Baryon Oscillation Spectroscopic Survey, we show that improvements vary from 10% to 100% depending on the cosmological parameter considered, the galaxy sample, and the survey volume.Non-linear redshift-space distortions ("fingers of God") are challenging to model analytically, a fact that limits the applicability of perturbation theory in redshift space as compared to real space. We show how this problem can be mitigated using a new observable, $Q_0$, which can be easily estimated from the redshift space clustering data and is approximately equal to the real space power spectrum. The new statistic does not suffer from fingers of God and can be accurately described with perturbation theory down to $k_{\rm max}\simeq 0.4~h~\text{Mpc}^{-1}$. It can be straightforwardly included in the likelihood at negligible additional computational cost, and yields noticeable improvements on cosmological parameters compared to standard power spectrum multipole analyses. Using both simulations and observational data from the Baryon Oscillation Spectroscopic Survey, we show that improvements vary from $10\%$ to $100\%$ depending on the cosmological parameter considered, the galaxy sample and the survey volume

Precision analysis of the redshift-space galaxy bispectrum

We study the information content of the angle-averaged redshift space galaxy bispectrum. The main novelty of our approach is the use of a systematic tree-level perturbation theory model that includes galaxy bias, IR resummation, and also accounts for nonlinear redshift space distortions, binning, and projection effects. We analyze data from the perturbation theory challenge simulations, whose cumulative volume of 566  h-3 Gpc3 allows for a precise comparison to theoretical predictions. Fitting the power spectrum and bispectrum of our simulated data, and varying all necessary cosmological and nuisance parameters in a consistent Markov chain Monte Carlo analysis, we find that our tree-level bispectrum model is valid up to kmax=0.08  h Mpc-1 (at z=0.61). We also find that inclusion of the bispectrum monopole improves constraints on cosmological parameters by (5–15)% relative to the power spectrum. The improvement is more significant for the quadratic bias parameters of our simulated galaxies, which we also show to deviate from biases of the host dark matter halos at the ∼3σ level. Finally, we adjust the covariance and scale cuts to match the volume of the BOSS survey, and estimate that within the minimal ΛCDM model the bispectrum data can tighten the constraint on the mass fluctuation amplitude σ8 by roughly 10%.We study the information content of the angle-averaged (monopole) redshift space galaxy bispectrum. The main novelty of our approach is the use of a systematic tree-level perturbation theory model that includes galaxy bias, IR resummation, and also accounts for nonlinear redshift space distortions, binning, and projection effects. We analyze data from the PT challenge simulations, whose cumulative volume of 566 $h^{-3}$Gpc$^3$ allows for a precise comparison to theoretical predictions. Fitting the power spectrum and bispectrum of our simulated data, and varying all necessary cosmological and nuisance parameters in a consistent Markov chain Monte Carlo analysis, we find that our tree-level bispectrum model is valid up to $k_{\max}=0.08~h{\rm Mpc}^{-1}$ (at $z=0.61$). We also find that inclusion of the bispectrum monopole improves constraints on cosmological parameters by $(5-15)\%$ relative to the power spectrum. The improvement is more significant for the quadratic bias parameters of our simulated galaxies, which we also show to deviate from biases of the host dark matter halos at the $\sim 3\sigma$ level. Finally, we adjust the covariance and scale cuts to match the volume of the BOSS survey, and estimate that within the minimal $\Lambda$CDM model the bispectrum data can tighten the constraint on the mass fluctuation amplitude $\sigma_8$ by roughly $10\%$

A Parameter-Masked Mock Data Challenge for Beyond-Two-Point Galaxy Clustering Statistics

International audienceThe last few years have seen the emergence of a wide array of novel techniques for analyzing high-precision data from upcoming galaxy surveys, which aim to extend the statistical analysis of galaxy clustering data beyond the linear regime and the canonical two-point (2pt) statistics. We test and benchmark some of these new techniques in a community data challenge "Beyond-2pt", initiated during the Aspen 2022 Summer Program "Large-Scale Structure Cosmology beyond 2-Point Statistics," whose first round of results we present here. The challenge dataset consists of high-precision mock galaxy catalogs for clustering in real space, redshift space, and on a light cone. Participants in the challenge have developed end-to-end pipelines to analyze mock catalogs and extract unknown ("masked") cosmological parameters of the underlying $\Lambda$CDM models with their methods. The methods represented are density-split clustering, nearest neighbor statistics, BACCO power spectrum emulator, void statistics, LEFTfield field-level inference using effective field theory (EFT), and joint power spectrum and bispectrum analyses using both EFT and simulation-based inference. In this work, we review the results of the challenge, focusing on problems solved, lessons learned, and future research needed to perfect the emerging beyond-2pt approaches. The unbiased parameter recovery demonstrated in this challenge by multiple statistics and the associated modeling and inference frameworks supports the credibility of cosmology constraints from these methods. The challenge data set is publicly available and we welcome future submissions from methods that are not yet represented

The CAMELS project: public data release

The Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS) project was developed to combine cosmology with astrophysics through thousands of cosmological hydrodynamic simulations and machine learning. CAMELS contains 4,233 cosmological simulations, 2,049 N-body and 2,184 state-of-the-art hydrodynamic simulations that sample a vast volume in parameter space. In this paper we present the CAMELS public data release, describing the characteristics of the CAMELS simulations and a variety of data products generated from them, including halo, subhalo, galaxy, and void catalogues, power spectra, bispectra, Lyman-$\alpha$ spectra, probability distribution functions, halo radial profiles, and X-rays photon lists. We also release over one thousand catalogues that contain billions of galaxies from CAMELS-SAM: a large collection of N-body simulations that have been combined with the Santa Cruz Semi-Analytic Model. We release all the data, comprising more than 350 terabytes and containing 143,922 snapshots, millions of halos, galaxies and summary statistics. We provide further technical details on how to access, download, read, and process the data at \url{https://camels.readthedocs.io}

Cosmology intertwined: A review of the particle physics, astrophysics, and cosmology associated with the cosmological tensions and anomalies

The standard Λ Cold Dark Matter (ΛCDM) cosmological model provides a good description of a wide range of astrophysical and cosmological data. However, there are a few big open questions that make the standard model look like an approximation to a more realistic scenario yet to be found. In this paper, we list a few important goals that need to be addressed in the next decade, taking into account the current discordances between the different cosmological probes, such as the disagreement in the value of the Hubble constant H0, the σ8–S8 tension, and other less statistically significant anomalies. While these discordances can still be in part the result of systematic errors, their persistence after several years of accurate analysis strongly hints at cracks in the standard cosmological scenario and the necessity for new physics or generalisations beyond the standard model. In this paper, we focus on the 5.0σ tension between the Planck CMB estimate of the Hubble constant H0 and the SH0ES collaboration measurements. After showing the H0 evaluations made from different teams using different methods and geometric calibrations, we list a few interesting new physics models that could alleviate this tension and discuss how the next decade's experiments will be crucial. Moreover, we focus on the tension of the Planck CMB data with weak lensing measurements and redshift surveys, about the value of the matter energy density Ωm, and the amplitude or rate of the growth of structure (σ8,fσ8). We list a few interesting models proposed for alleviating this tension, and we discuss the importance of trying to fit a full array of data with a single model and not just one parameter at a time. Additionally, we present a wide range of other less discussed anomalies at a statistical significance level lower than the H0–S8 tensions which may also constitute hints towards new physics, and we discuss possible generic theoretical approaches that can collectively explain the non-standard nature of these signals. Finally, we give an overview of upgraded experiments and next-generation space missions and facilities on Earth that will be of crucial importance to address all these open questions