-
The Primordial Inflation Explorer (PIXIE): Mission Design and Science Goals
Authors:
Alan Kogut,
Eric Switzer,
Dale Fixsen,
Nabila Aghanim,
Jens Chluba,
Dave Chuss,
Jacques Delabrouille,
Cora Dvorkin,
Brandon Hensley,
Colin Hill,
Bruno Maffei,
Anthony Pullen,
Aditya Rotti,
Alina Sabyr,
Leander Thiele,
Ed Wollack,
Ioana Zelko
Abstract:
The Primordial Inflation Explorer (PIXIE) is an Explorer-class mission concept to measure the energy spectrum and linear polarization of the cosmic microwave background (CMB). A single cryogenic Fourier transform spectrometer compares the sky to an external blackbody calibration target, measuring the Stokes I, Q, U parameters to levels ~200 Jy/sr in each 2.65 degree diameter beam over the full sky…
▽ More
The Primordial Inflation Explorer (PIXIE) is an Explorer-class mission concept to measure the energy spectrum and linear polarization of the cosmic microwave background (CMB). A single cryogenic Fourier transform spectrometer compares the sky to an external blackbody calibration target, measuring the Stokes I, Q, U parameters to levels ~200 Jy/sr in each 2.65 degree diameter beam over the full sky, in each of 300 frequency channels from 28 GHz to 6 THz. With sensitivity over 1000 times greater than COBE/FIRAS, PIXIE opens a broad discovery space for the origin, contents, and evolution of the universe. Measurements of small distortions from a CMB blackbody spectrum provide a robust determination of the mean electron pressure and temperature in the universe while constraining processes including dissipation of primordial density perturbations, black holes, and the decay or annihilation of dark matter. Full-sky maps of linear polarization measure the optical depth to reionization at nearly the cosmic variance limit and constrain models of primordial inflation. Spectra with sub-percent absolute calibration spanning microwave to far-IR wavelengths provide a legacy data set for analyses including line intensity mapping of extragalactic emission and the cosmic infrared background amplitude and anisotropy. We describe the PIXIE instrument sensitivity, foreground subtraction, and anticipated science return from both the baseline 2-year mission and a potential extended mission.
△ Less
Submitted 30 May, 2024;
originally announced May 2024.
-
Cosmological constraints from weak lensing scattering transform using HSC Y1 data
Authors:
Sihao Cheng,
Gabriela A. Marques,
Daniela Grandón,
Leander Thiele,
Masato Shirasaki,
Brice Ménard,
Jia Liu
Abstract:
As weak lensing surveys go deeper, there is an increasing need for reliable characterization of non-Gaussian structures at small angular scales. Here we present the first cosmological constraints with weak lensing scattering transform, a statistical estimator that combines efficiency, robustness, and interpretability. With the Hyper Suprime-Cam survey (HSC) year 1 data, we obtain…
▽ More
As weak lensing surveys go deeper, there is an increasing need for reliable characterization of non-Gaussian structures at small angular scales. Here we present the first cosmological constraints with weak lensing scattering transform, a statistical estimator that combines efficiency, robustness, and interpretability. With the Hyper Suprime-Cam survey (HSC) year 1 data, we obtain $Ω_\text{m}=0.29_{-0.03}^{+0.04}$, $S_8\equiv σ_8(Ω_\text{m}/0.3)^{0.5}=0.83\pm0.02$, and intrinsic alignment strength $A_\text{IA}=1.0\pm0.4$ through simulation-based forward modeling. Our constraints are consistent with those derived from Planck. The error bar of $Ω_\text{m}$ is 2 times tighter than that obtained from the power spectrum when the same scale range is used. This constraining power is on par with that of convolutional neural networks, suggesting that further investment in spatial information extraction may not yield substantial benefits.
We also point out an internal tension of $S_8$ estimates linked to a redshift bin around z ~ 1 in the HSC data. We found that discarding that bin leads to a consistent decrease of $S_8$ from 0.83 to 0.79, for all statistical estimators. We argue that photometric redshift estimation is now the main limitation in the estimation of $S_8$ using HSC. This limitation is likely to affect other ground-based weak lensing surveys reaching redshifts greater than one. Alternative redshift estimation techniques, like clustering redshifts, may help alleviate this limitation.
△ Less
Submitted 24 April, 2024;
originally announced April 2024.
-
Impact of baryonic feedback on HSC Y1 weak lensing non-Gaussian statistics
Authors:
Daniela Grandón,
Gabriela A. Marques,
Leander Thiele,
Sihao Cheng,
Masato Shirasaki,
Jia Liu
Abstract:
Baryonic feedback is a major systematic in weak lensing cosmology. Its most studied effect is the suppression of the lensing power spectrum, a second-order statistic, on small scales. Motivated by the growing interest in statistics beyond the second order, we investigate the effect of baryons on lensing non-Gaussian statistics and the resulting biases in the matter clustering amplitude…
▽ More
Baryonic feedback is a major systematic in weak lensing cosmology. Its most studied effect is the suppression of the lensing power spectrum, a second-order statistic, on small scales. Motivated by the growing interest in statistics beyond the second order, we investigate the effect of baryons on lensing non-Gaussian statistics and the resulting biases in the matter clustering amplitude $S_8 = σ_8\sqrt{Ω_m/0.3}$. We focus on the Subaru Hyper Suprime-Cam Year 1 (HSC-Y1) data which, with its high source number density, closely resembles those expected from the upcoming Euclid and Rubin LSST. We study four non-Gaussian statistics -- peak counts, minimum counts, the probability distribution function, and the scattering transform -- in addition to the usual power spectrum. We first estimate the biases in $S_8$ using mock observations built from the IllustrisTNG and BAHAMAS hydrodynamical simulations and theoretical models built from dark matter-only simulations. We find up to $1σ$ bias in $S_8$ when the smallest scales (2 arcmin) and the highest feedback level are considered. We then analyze the HSC-Y1 data and compare the $S_8$ obtained for each statistic with different smoothing scales or scale cuts. As we expect that baryons mostly affect the small scales, comparing the results obtained from including and excluding small scales can indicate the level of impact from baryons. With HSC data, we find only minor ($\leq0.5σ$) differences in $S_8$ for all statistics, even when considering very small scales (2 arcmin). Our results suggest that the effect of baryons is insignificant at the level of HSC-Y1 down to 2~arcmin for all statistics examined here, or it is canceled by other scale-dependent systematics.
△ Less
Submitted 6 March, 2024;
originally announced March 2024.
-
Cosmology from weak lensing peaks and minima with Subaru Hyper Suprime-Cam survey first-year data
Authors:
Gabriela A. Marques,
Jia Liu,
Masato Shirasaki,
Leander Thiele,
Daniela Grandón,
Kevin M. Huffenberger,
Sihao Cheng,
Joachim Harnois-Déraps,
Ken Osato,
William R. Coulton
Abstract:
We present cosmological constraints derived from peak counts, minimum counts, and the angular power spectrum of the Subaru Hyper Suprime-Cam first-year (HSC Y1) weak lensing shear catalog. Weak lensing peak and minimum counts contain non-Gaussian information and hence are complementary to the conventional two-point statistics in constraining cosmology. In this work, we forward-model the three summ…
▽ More
We present cosmological constraints derived from peak counts, minimum counts, and the angular power spectrum of the Subaru Hyper Suprime-Cam first-year (HSC Y1) weak lensing shear catalog. Weak lensing peak and minimum counts contain non-Gaussian information and hence are complementary to the conventional two-point statistics in constraining cosmology. In this work, we forward-model the three summary statistics and their dependence on cosmology, using a suite of $N$-body simulations tailored to the HSC Y1 data. We investigate systematic and astrophysical effects including intrinsic alignments, baryon feedback, multiplicative bias, and photometric redshift uncertainties. We mitigate the impact of these systematics by applying cuts on angular scales, smoothing scales, statistic bins, and tomographic redshift bins. By combining peaks, minima, and the power spectrum, assuming a flat-$Λ$CDM model, we obtain $S_{8} \equiv σ_8\sqrt{Ω_m/0.3}= 0.810^{+0.022}_{-0.026}$, a 35\% tighter constraint than that obtained from the angular power spectrum alone. Our results are in agreement with other studies using HSC weak lensing shear data, as well as with Planck 2018 cosmology and recent CMB lensing constraints from the Atacama Cosmology Telescope and the South Pole Telescope.
△ Less
Submitted 21 August, 2023;
originally announced August 2023.
-
Neutrino mass constraint from an Implicit Likelihood Analysis of BOSS voids
Authors:
Leander Thiele,
Elena Massara,
Alice Pisani,
ChangHoon Hahn,
David N. Spergel,
Shirley Ho,
Benjamin Wandelt
Abstract:
Cosmic voids identified in the spatial distribution of galaxies provide complementary information to two-point statistics. In particular, constraints on the neutrino mass sum, $\sum m_ν$, promise to benefit from the inclusion of void statistics. We perform inference on the CMASS NGC sample of SDSS-III/BOSS with the aim of constraining $\sum m_ν$. We utilize the void size function, the void galaxy…
▽ More
Cosmic voids identified in the spatial distribution of galaxies provide complementary information to two-point statistics. In particular, constraints on the neutrino mass sum, $\sum m_ν$, promise to benefit from the inclusion of void statistics. We perform inference on the CMASS NGC sample of SDSS-III/BOSS with the aim of constraining $\sum m_ν$. We utilize the void size function, the void galaxy cross power spectrum, and the galaxy auto power spectrum. To extract constraints from these summary statistics we use a simulation-based approach, specifically implicit likelihood inference. We populate approximate gravity-only, particle neutrino cosmological simulations with an expressive halo occupation distribution model. With a conservative scale cut of kmax=0.15 h/Mpc and a Planck-inspired LCDM prior, we find upper bounds on $\sum m_ν$ of 0.43 and 0.35 eV from the galaxy auto power spectrum and the full data vector, respectively (95% credible interval). We observe hints that the void statistics may be most effective at constraining $\sum m_ν$ from below. We also substantiate the usual assumption that the void size function is Poisson distributed.
△ Less
Submitted 14 July, 2023;
originally announced July 2023.
-
Cosmological constraints from HSC Y1 lensing convergence PDF
Authors:
Leander Thiele,
Gabriela A. Marques,
Jia Liu,
Masato Shirasaki
Abstract:
We utilize the probability distribution function (PDF) of normalized convergence maps reconstructed from the Subaru Hyper Suprime-Cam (HSC) Y1 shear catalogue, in combination with the power spectrum, to measure the matter clustering amplitude $S_8=σ_8\sqrt{Ω_m/0.3}$. The large-scale structure's statistical properties are incompletely described by the traditional two-point statistics, motivating ou…
▽ More
We utilize the probability distribution function (PDF) of normalized convergence maps reconstructed from the Subaru Hyper Suprime-Cam (HSC) Y1 shear catalogue, in combination with the power spectrum, to measure the matter clustering amplitude $S_8=σ_8\sqrt{Ω_m/0.3}$. The large-scale structure's statistical properties are incompletely described by the traditional two-point statistics, motivating our investigation of the PDF -- a complementary higher-order statistic. By defining the PDF over the standard deviation-normalized convergence map we are able to isolate the non-Gaussian information. We use tailored simulations to compress the data vector and construct a likelihood approximation. We mitigate the impact of survey and astrophysical systematics with cuts on smoothing scales, redshift bins, and data vectors. We find $S_8=0.860^{+0.066}_{-0.109}$ from the PDF alone and $S_8=0.798^{+0.029}_{-0.042}$ from the combination of PDF and power spectrum (68% CL). The PDF improves the power spectrum-only constraint by about 10%.
△ Less
Submitted 21 November, 2023; v1 submitted 12 April, 2023;
originally announced April 2023.
-
Predicting the impact of feedback on matter clustering with machine learning in CAMELS
Authors:
Ana Maria Delgado,
Daniel Angles-Alcazar,
Leander Thiele,
Shivam Pandey,
Kai Lehman,
Rachel S. Somerville,
Michelle Ntampaka,
Shy Genel,
Francisco Villaescusa-Navarro,
Lars Hernquist
Abstract:
Extracting information from the total matter power spectrum with the precision needed for upcoming cosmological surveys requires unraveling the complex effects of galaxy formation processes on the distribution of matter. We investigate the impact of baryonic physics on matter clustering at $z=0$ using a library of power spectra from the Cosmology and Astrophysics with MachinE Learning Simulations…
▽ More
Extracting information from the total matter power spectrum with the precision needed for upcoming cosmological surveys requires unraveling the complex effects of galaxy formation processes on the distribution of matter. We investigate the impact of baryonic physics on matter clustering at $z=0$ using a library of power spectra from the Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS) project, containing thousands of $(25\,h^{-1}{\rm Mpc})^3$ volume realizations with varying cosmology, initial random field, stellar and AGN feedback strength and sub-grid model implementation methods. We show that baryonic physics affects matter clustering on scales $k \gtrsim 0.4\,h\,\mathrm{Mpc}^{-1}$ and the magnitude of this effect is dependent on the details of the galaxy formation implementation and variations of cosmological and astrophysical parameters. Increasing AGN feedback strength decreases halo baryon fractions and yields stronger suppression of power relative to N-body simulations, while stronger stellar feedback often results in weaker effects by suppressing black hole growth and therefore the impact of AGN feedback. We find a broad correlation between mean baryon fraction of massive halos ($M_{\rm 200c} > 10^{13.5}$\,\Msun) and suppression of matter clustering but with significant scatter compared to previous work owing to wider exploration of feedback parameters and cosmic variance effects. We show that a random forest regressor trained on the baryon content and abundance of halos across the full mass range $10^{10} \leq M_\mathrm{halo}/$\Msun$< 10^{15}$ can predict the effect of galaxy formation on the matter power spectrum on scales $k = 1.0$--20.0\,$h\,\mathrm{Mpc}^{-1}$.
△ Less
Submitted 5 October, 2023; v1 submitted 5 January, 2023;
originally announced January 2023.
-
The SZ flux-mass ($Y$-$M$) relation at low halo masses: improvements with symbolic regression and strong constraints on baryonic feedback
Authors:
Digvijay Wadekar,
Leander Thiele,
J. Colin Hill,
Shivam Pandey,
Francisco Villaescusa-Navarro,
David N. Spergel,
Miles Cranmer,
Daisuke Nagai,
Daniel Anglés-Alcázar,
Shirley Ho,
Lars Hernquist
Abstract:
Feedback from active galactic nuclei (AGN) and supernovae can affect measurements of integrated SZ flux of halos ($Y_\mathrm{SZ}$) from CMB surveys, and cause its relation with the halo mass ($Y_\mathrm{SZ}-M$) to deviate from the self-similar power-law prediction of the virial theorem. We perform a comprehensive study of such deviations using CAMELS, a suite of hydrodynamic simulations with exten…
▽ More
Feedback from active galactic nuclei (AGN) and supernovae can affect measurements of integrated SZ flux of halos ($Y_\mathrm{SZ}$) from CMB surveys, and cause its relation with the halo mass ($Y_\mathrm{SZ}-M$) to deviate from the self-similar power-law prediction of the virial theorem. We perform a comprehensive study of such deviations using CAMELS, a suite of hydrodynamic simulations with extensive variations in feedback prescriptions. We use a combination of two machine learning tools (random forest and symbolic regression) to search for analogues of the $Y-M$ relation which are more robust to feedback processes for low masses ($M\lesssim 10^{14}\, h^{-1} \, M_\odot$); we find that simply replacing $Y\rightarrow Y(1+M_*/M_\mathrm{gas})$ in the relation makes it remarkably self-similar. This could serve as a robust multiwavelength mass proxy for low-mass clusters and galaxy groups. Our methodology can also be generally useful to improve the domain of validity of other astrophysical scaling relations.
We also forecast that measurements of the $Y-M$ relation could provide percent-level constraints on certain combinations of feedback parameters and/or rule out a major part of the parameter space of supernova and AGN feedback models used in current state-of-the-art hydrodynamic simulations. Our results can be useful for using upcoming SZ surveys (e.g., SO, CMB-S4) and galaxy surveys (e.g., DESI and Rubin) to constrain the nature of baryonic feedback. Finally, we find that the an alternative relation, $Y-M_*$, provides complementary information on feedback than $Y-M$
△ Less
Submitted 28 April, 2023; v1 submitted 5 September, 2022;
originally announced September 2022.
-
An exploration of the properties of cluster profiles for the thermal and kinetic Sunyaev-Zel'dovich effects
Authors:
Billy K. K. Lee,
William R. Coulton,
Leander Thiele,
Shirley Ho
Abstract:
With the advent of high-resolution, low-noise CMB measurements, the ability to extract cosmological information from thermal Sunyaev-Zel'dovich effect and kinetic Sunyaev-Zel'dovich effect will be limited not by statistical uncertainties but rather by systematic and theoretical uncertainties. The theoretical uncertainty is driven by the lack of knowledge about the electron pressure and density. Th…
▽ More
With the advent of high-resolution, low-noise CMB measurements, the ability to extract cosmological information from thermal Sunyaev-Zel'dovich effect and kinetic Sunyaev-Zel'dovich effect will be limited not by statistical uncertainties but rather by systematic and theoretical uncertainties. The theoretical uncertainty is driven by the lack of knowledge about the electron pressure and density. Thus we explore the electron pressure and density distributions in the IllustrisTNG hydrodynamical simulations, and we demonstrate that the cluster properties exhibit a strong dependence on the halo concentration -- providing some of the first evidence of cluster assembly bias in the electron pressure and density. Further, our work shows evidence for a broken power-law mass dependence, with lower pressure in lower mass halos than previous work and a strong evolution with mass of the radial correlations in the electron density and pressure. Both of these effects highlight the differing impact of active galactic nuclei and supernova feedback on the gas in galaxy groups compared to massive clusters. We verified that we see qualitatively similar features in the SIMBA hydro-dynamical simulations, suggesting these effects could be generic features. Finally, we provide a parametric formula for the electron pressure and density profile as a function of dark matter halo mass, halo concentration, and redshift. These fitting formulae can reproduce the distribution of density and pressure of clusters and will be useful in extracting cosmological information from upcoming CMB surveys.
△ Less
Submitted 3 May, 2022;
originally announced May 2022.
-
Predicting the Thermal Sunyaev-Zel'dovich Field using Modular and Equivariant Set-Based Neural Networks
Authors:
Leander Thiele,
Miles Cranmer,
William Coulton,
Shirley Ho,
David N. Spergel
Abstract:
Theoretical uncertainty limits our ability to extract cosmological information from baryonic fields such as the thermal Sunyaev-Zel'dovich (tSZ) effect. Being sourced by the electron pressure field, the tSZ effect depends on baryonic physics that is usually modeled by expensive hydrodynamic simulations. We train neural networks on the IllustrisTNG-300 cosmological simulation to predict the continu…
▽ More
Theoretical uncertainty limits our ability to extract cosmological information from baryonic fields such as the thermal Sunyaev-Zel'dovich (tSZ) effect. Being sourced by the electron pressure field, the tSZ effect depends on baryonic physics that is usually modeled by expensive hydrodynamic simulations. We train neural networks on the IllustrisTNG-300 cosmological simulation to predict the continuous electron pressure field in galaxy clusters from gravity-only simulations. Modeling clusters is challenging for neural networks as most of the gas pressure is concentrated in a handful of voxels and even the largest hydrodynamical simulations contain only a few hundred clusters that can be used for training. Instead of conventional convolutional neural net (CNN) architectures, we choose to employ a rotationally equivariant DeepSets architecture to operate directly on the set of dark matter particles. We argue that set-based architectures provide distinct advantages over CNNs. For example, we can enforce exact rotational and permutation equivariance, incorporate existing knowledge on the tSZ field, and work with sparse fields as are standard in cosmology. We compose our architecture with separate, physically meaningful modules, making it amenable to interpretation. For example, we can separately study the influence of local and cluster-scale environment, determine that cluster triaxiality has negligible impact, and train a module that corrects for mis-centering. Our model improves by 70 % on analytic profiles fit to the same simulation data. We argue that the electron pressure field, viewed as a function of a gravity-only simulation, has inherent stochasticity, and model this property through a conditional-VAE extension to the network. This modification yields further improvement by 7 %, it is limited by our small training set however. (abridged)
△ Less
Submitted 10 June, 2022; v1 submitted 28 February, 2022;
originally announced March 2022.
-
Percent-level constraints on baryonic feedback with spectral distortion measurements
Authors:
Leander Thiele,
Digvijay Wadekar,
J. Colin Hill,
Nicholas Battaglia,
Jens Chluba,
Francisco Villaescusa-Navarro,
Lars Hernquist,
Mark Vogelsberger,
Daniel Anglés-Alcázar,
Federico Marinacci
Abstract:
High-significance measurements of the monopole thermal Sunyaev-Zel'dovich CMB spectral distortions have the potential to tightly constrain poorly understood baryonic feedback processes. The sky-averaged Compton-y distortion and its relativistic correction are measures of the total thermal energy in electrons in the observable universe and their mean temperature. We use the CAMELS suite of hydrodyn…
▽ More
High-significance measurements of the monopole thermal Sunyaev-Zel'dovich CMB spectral distortions have the potential to tightly constrain poorly understood baryonic feedback processes. The sky-averaged Compton-y distortion and its relativistic correction are measures of the total thermal energy in electrons in the observable universe and their mean temperature. We use the CAMELS suite of hydrodynamic simulations to explore possible constraints on parameters describing the subgrid implementation of feedback from active galactic nuclei and supernovae, assuming a PIXIE-like measurement. The small 25 Mpc/h CAMELS boxes present challenges due to the significant cosmic variance. We utilize machine learning to construct interpolators through the noisy simulation data. Using the halo model, we translate the simulation halo mass functions into correction factors to reduce cosmic variance where required. Our results depend on the subgrid model. In the case of IllustrisTNG, we find that the best-determined parameter combination can be measured to ~2% and corresponds to a product of AGN and SN feedback. In the case of SIMBA, the tightest constraint is ~0.2% on a ratio between AGN and SN feedback. A second orthogonal parameter combination can be measured to ~8%. Our results demonstrate the significant constraining power a measurement of the late-time spectral distortion monopoles would have for baryonic feedback models.
△ Less
Submitted 27 April, 2022; v1 submitted 5 January, 2022;
originally announced January 2022.
-
Augmenting astrophysical scaling relations with machine learning: application to reducing the Sunyaev-Zeldovich flux-mass scatter
Authors:
Digvijay Wadekar,
Leander Thiele,
Francisco Villaescusa-Navarro,
J. Colin Hill,
Miles Cranmer,
David N. Spergel,
Nicholas Battaglia,
Daniel Anglés-Alcázar,
Lars Hernquist,
Shirley Ho
Abstract:
Complex astrophysical systems often exhibit low-scatter relations between observable properties (e.g., luminosity, velocity dispersion, oscillation period). These scaling relations illuminate the underlying physics, and can provide observational tools for estimating masses and distances. Machine learning can provide a fast and systematic way to search for new scaling relations (or for simple exten…
▽ More
Complex astrophysical systems often exhibit low-scatter relations between observable properties (e.g., luminosity, velocity dispersion, oscillation period). These scaling relations illuminate the underlying physics, and can provide observational tools for estimating masses and distances. Machine learning can provide a fast and systematic way to search for new scaling relations (or for simple extensions to existing relations) in abstract high-dimensional parameter spaces. We use a machine learning tool called symbolic regression (SR), which models patterns in a dataset in the form of analytic equations. We focus on the Sunyaev-Zeldovich flux$-$cluster mass relation ($Y_\mathrm{SZ}-M$), the scatter in which affects inference of cosmological parameters from cluster abundance data. Using SR on the data from the IllustrisTNG hydrodynamical simulation, we find a new proxy for cluster mass which combines $Y_\mathrm{SZ}$ and concentration of ionized gas ($c_\mathrm{gas}$): $M \propto Y_\mathrm{conc}^{3/5} \equiv Y_\mathrm{SZ}^{3/5} (1-A\, c_\mathrm{gas})$. $Y_\mathrm{conc}$ reduces the scatter in the predicted $M$ by $\sim 20-30$\% for large clusters ($M\gtrsim 10^{14}\, h^{-1} \, M_\odot$), as compared to using just $Y_\mathrm{SZ}$. We show that the dependence on $c_\mathrm{gas}$ is linked to cores of clusters exhibiting larger scatter than their outskirts. Finally, we test $Y_\mathrm{conc}$ on clusters from CAMELS simulations and show that $Y_\mathrm{conc}$ is robust against variations in cosmology, subgrid physics, and cosmic variance. Our results and methodology can be useful for accurate multiwavelength cluster mass estimation from upcoming CMB and X-ray surveys like ACT, SO, eROSITA and CMB-S4.
△ Less
Submitted 17 March, 2023; v1 submitted 4 January, 2022;
originally announced January 2022.
-
The CAMELS project: public data release
Authors:
Francisco Villaescusa-Navarro,
Shy Genel,
Daniel Anglés-Alcázar,
Lucia A. Perez,
Pablo Villanueva-Domingo,
Digvijay Wadekar,
Helen Shao,
Faizan G. Mohammad,
Sultan Hassan,
Emily Moser,
Erwin T. Lau,
Luis Fernando Machado Poletti Valle,
Andrina Nicola,
Leander Thiele,
Yongseok Jo,
Oliver H. E. Philcox,
Benjamin D. Oppenheimer,
Megan Tillman,
ChangHoon Hahn,
Neerav Kaushal,
Alice Pisani,
Matthew Gebhardt,
Ana Maria Delgado,
Joyce Caliendo,
Christina Kreisch
, et al. (22 additional authors not shown)
Abstract:
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…
▽ More
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-$α$ 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}.
△ Less
Submitted 4 January, 2022;
originally announced January 2022.
-
BISOU: a balloon project to measure the CMB spectral distortions
Authors:
B. Maffei,
M. H. Abitbol,
N. Aghanim,
J. Aumont,
E. Battistelli,
J. Chluba,
X. Coulon,
P. De Bernardis,
M. Douspis,
J. Grain,
S. Gervasoni,
J. C. Hill,
A. Kogut,
S. Masi,
T. Matsumura,
C. O Sullivan,
L. Pagano,
G. Pisano,
M. Remazeilles,
A. Ritacco,
A. Rotti,
V. Sauvage,
G. Savini,
S. L. Stever,
A. Tartari
, et al. (2 additional authors not shown)
Abstract:
The BISOU (Balloon Interferometer for Spectral Observations of the Universe) project aims to study the viability and prospects of a balloon-borne spectrometer, pathfinder of a future space mission dedicated to the measurements of the CMB spectral distortions. We present here a preliminary concept based on previous space mission proposals, together with some sensitivity calculation results for the…
▽ More
The BISOU (Balloon Interferometer for Spectral Observations of the Universe) project aims to study the viability and prospects of a balloon-borne spectrometer, pathfinder of a future space mission dedicated to the measurements of the CMB spectral distortions. We present here a preliminary concept based on previous space mission proposals, together with some sensitivity calculation results for the observation goals, showing that a 5-sigma measurement of the y-distortions is achievable.
△ Less
Submitted 30 October, 2021;
originally announced November 2021.
-
The CAMELS Multifield Dataset: Learning the Universe's Fundamental Parameters with Artificial Intelligence
Authors:
Francisco Villaescusa-Navarro,
Shy Genel,
Daniel Angles-Alcazar,
Leander Thiele,
Romeel Dave,
Desika Narayanan,
Andrina Nicola,
Yin Li,
Pablo Villanueva-Domingo,
Benjamin Wandelt,
David N. Spergel,
Rachel S. Somerville,
Jose Manuel Zorrilla Matilla,
Faizan G. Mohammad,
Sultan Hassan,
Helen Shao,
Digvijay Wadekar,
Michael Eickenberg,
Kaze W. K. Wong,
Gabriella Contardo,
Yongseok Jo,
Emily Moser,
Erwin T. Lau,
Luis Fernando Machado Poletti Valle,
Lucia A. Perez
, et al. (3 additional authors not shown)
Abstract:
We present the Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS) Multifield Dataset, CMD, a collection of hundreds of thousands of 2D maps and 3D grids containing many different properties of cosmic gas, dark matter, and stars from 2,000 distinct simulated universes at several cosmic times. The 2D maps and 3D grids represent cosmic regions that span $\sim$100 million light year…
▽ More
We present the Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS) Multifield Dataset, CMD, a collection of hundreds of thousands of 2D maps and 3D grids containing many different properties of cosmic gas, dark matter, and stars from 2,000 distinct simulated universes at several cosmic times. The 2D maps and 3D grids represent cosmic regions that span $\sim$100 million light years and have been generated from thousands of state-of-the-art hydrodynamic and gravity-only N-body simulations from the CAMELS project. Designed to train machine learning models, CMD is the largest dataset of its kind containing more than 70 Terabytes of data. In this paper we describe CMD in detail and outline a few of its applications. We focus our attention on one such task, parameter inference, formulating the problems we face as a challenge to the community. We release all data and provide further technical details at https://camels-multifield-dataset.readthedocs.io.
△ Less
Submitted 22 September, 2021;
originally announced September 2021.
-
Robust marginalization of baryonic effects for cosmological inference at the field level
Authors:
Francisco Villaescusa-Navarro,
Shy Genel,
Daniel Angles-Alcazar,
David N. Spergel,
Yin Li,
Benjamin Wandelt,
Leander Thiele,
Andrina Nicola,
Jose Manuel Zorrilla Matilla,
Helen Shao,
Sultan Hassan,
Desika Narayanan,
Romeel Dave,
Mark Vogelsberger
Abstract:
We train neural networks to perform likelihood-free inference from $(25\,h^{-1}{\rm Mpc})^2$ 2D maps containing the total mass surface density from thousands of hydrodynamic simulations of the CAMELS project. We show that the networks can extract information beyond one-point functions and power spectra from all resolved scales ($\gtrsim 100\,h^{-1}{\rm kpc}$) while performing a robust marginalizat…
▽ More
We train neural networks to perform likelihood-free inference from $(25\,h^{-1}{\rm Mpc})^2$ 2D maps containing the total mass surface density from thousands of hydrodynamic simulations of the CAMELS project. We show that the networks can extract information beyond one-point functions and power spectra from all resolved scales ($\gtrsim 100\,h^{-1}{\rm kpc}$) while performing a robust marginalization over baryonic physics at the field level: the model can infer the value of $Ω_{\rm m} (\pm 4\%)$ and $σ_8 (\pm 2.5\%)$ from simulations completely different to the ones used to train it.
△ Less
Submitted 21 September, 2021;
originally announced September 2021.
-
Multifield Cosmology with Artificial Intelligence
Authors:
Francisco Villaescusa-Navarro,
Daniel Anglés-Alcázar,
Shy Genel,
David N. Spergel,
Yin Li,
Benjamin Wandelt,
Andrina Nicola,
Leander Thiele,
Sultan Hassan,
Jose Manuel Zorrilla Matilla,
Desika Narayanan,
Romeel Dave,
Mark Vogelsberger
Abstract:
Astrophysical processes such as feedback from supernovae and active galactic nuclei modify the properties and spatial distribution of dark matter, gas, and galaxies in a poorly understood way. This uncertainty is one of the main theoretical obstacles to extract information from cosmological surveys. We use 2,000 state-of-the-art hydrodynamic simulations from the CAMELS project spanning a wide vari…
▽ More
Astrophysical processes such as feedback from supernovae and active galactic nuclei modify the properties and spatial distribution of dark matter, gas, and galaxies in a poorly understood way. This uncertainty is one of the main theoretical obstacles to extract information from cosmological surveys. We use 2,000 state-of-the-art hydrodynamic simulations from the CAMELS project spanning a wide variety of cosmological and astrophysical models and generate hundreds of thousands of 2-dimensional maps for 13 different fields: from dark matter to gas and stellar properties. We use these maps to train convolutional neural networks to extract the maximum amount of cosmological information while marginalizing over astrophysical effects at the field level. Although our maps only cover a small area of $(25~h^{-1}{\rm Mpc})^2$, and the different fields are contaminated by astrophysical effects in very different ways, our networks can infer the values of $Ω_{\rm m}$ and $σ_8$ with a few percent level precision for most of the fields. We find that the marginalization performed by the network retains a wealth of cosmological information compared to a model trained on maps from gravity-only N-body simulations that are not contaminated by astrophysical effects. Finally, we train our networks on multifields -- 2D maps that contain several fields as different colors or channels -- and find that not only they can infer the value of all parameters with higher accuracy than networks trained on individual fields, but they can constrain the value of $Ω_{\rm m}$ with higher accuracy than the maps from the N-body simulations.
△ Less
Submitted 20 September, 2021;
originally announced September 2021.
-
Can small-scale baryon inhomogeneities resolve the Hubble tension? An investigation with ACT DR4
Authors:
Leander Thiele,
Yilun Guan,
J. Colin Hill,
Arthur Kosowsky,
David N. Spergel
Abstract:
Small-scale inhomogeneities in the baryon density around recombination have been proposed as a solution to the tension between local and global determinations of the Hubble constant. These baryon clumping models make distinct predictions for the cosmic microwave background anisotropy power spectra on small angular scales. We use recent data from the Atacama Cosmology Telescope to test these predic…
▽ More
Small-scale inhomogeneities in the baryon density around recombination have been proposed as a solution to the tension between local and global determinations of the Hubble constant. These baryon clumping models make distinct predictions for the cosmic microwave background anisotropy power spectra on small angular scales. We use recent data from the Atacama Cosmology Telescope to test these predictions. No evidence for baryon clumping is found, assuming a range of parameterizations for time-independent baryon density probability distribution functions. The inferred Hubble constant remains in significant tension with the SH0ES measurement.
△ Less
Submitted 10 August, 2021; v1 submitted 6 May, 2021;
originally announced May 2021.
-
Accurate Analytic Model for the Weak Lensing Convergence One-Point Probability Distribution Function and its Auto-Covariance
Authors:
Leander Thiele,
J. Colin Hill,
Kendrick M. Smith
Abstract:
The one-point probability distribution function (PDF) is a powerful summary statistic for non-Gaussian cosmological fields, such as the weak lensing (WL) convergence reconstructed from galaxy shapes or cosmic microwave background (CMB) maps. Thus far, no analytic model has been developed that successfully describes the high-convergence tail of the WL convergence PDF for small smoothing scales from…
▽ More
The one-point probability distribution function (PDF) is a powerful summary statistic for non-Gaussian cosmological fields, such as the weak lensing (WL) convergence reconstructed from galaxy shapes or cosmic microwave background (CMB) maps. Thus far, no analytic model has been developed that successfully describes the high-convergence tail of the WL convergence PDF for small smoothing scales from first principles. Here, we present a halo-model formalism to compute the WL convergence PDF, building upon our previous results for the thermal Sunyaev-Zel'dovich field. Furthermore, we extend our formalism to analytically compute the covariance matrix of the convergence PDF. Comparisons to numerical simulations generally confirm the validity of our formalism in the non-Gaussian, positive tail of the WL convergence PDF, but also reveal the convergence PDF's strong sensitivity to small-scale systematic effects in the simulations (e.g., due to finite resolution). Finally, we present a simple Fisher forecast for a Rubin-Observatory-like survey, based on our new analytic model. Considering the $\{A_s, Ω_m, Σm_ν\}$ parameter space and assuming a Planck CMB prior on $A_s$ only, we forecast a marginalized constraint $σ(Σm_ν) \approx 0.08$ eV from the WL convergence PDF alone, even after marginalizing over parameters describing the halo concentration-mass relation. This error bar on the neutrino mass sum is comparable to the minimum value allowed in the normal hierarchy, illustrating the strong constraining power of the WL convergence PDF. We make our code publicly available at https://github.com/leanderthiele/hmpdf.
△ Less
Submitted 23 December, 2020; v1 submitted 14 September, 2020;
originally announced September 2020.
-
Teaching neural networks to generate Fast Sunyaev Zel'dovich Maps
Authors:
Leander Thiele,
Francisco Villaescusa-Navarro,
David N. Spergel,
Dylan Nelson,
Annalisa Pillepich
Abstract:
The thermal Sunyaev-Zel'dovich (tSZ) and the kinematic Sunyaev-Zel'dovich (kSZ) effects trace the distribution of electron pressure and momentum in the hot Universe. These observables depend on rich multi-scale physics, thus, simulated maps should ideally be based on calculations that capture baryonic feedback effects such as cooling, star formation, and other complex processes. In this paper, we…
▽ More
The thermal Sunyaev-Zel'dovich (tSZ) and the kinematic Sunyaev-Zel'dovich (kSZ) effects trace the distribution of electron pressure and momentum in the hot Universe. These observables depend on rich multi-scale physics, thus, simulated maps should ideally be based on calculations that capture baryonic feedback effects such as cooling, star formation, and other complex processes. In this paper, we train deep convolutional neural networks with a U-Net architecture to map from the three-dimensional distribution of dark matter to electron density, momentum and pressure at ~ 100 kpc resolution. These networks are trained on a combination of the TNG300 volume and a set of cluster zoom-in simulations from the IllustrisTNG project. The neural nets are able to reproduce the power spectrum, one-point probability distribution function, bispectrum, and cross-correlation coefficients of the simulations more accurately than the state-of-the-art semi-analytical models. Our approach offers a route to capture the richness of a full cosmological hydrodynamical simulation of galaxy formation with the speed of an analytical calculation.
△ Less
Submitted 14 July, 2020;
originally announced July 2020.
-
Disentangling magnification in combined shear-clustering analyses
Authors:
Leander Thiele,
Christopher A. J. Duncan,
David Alonso
Abstract:
We investigate the sensitivity to the effects of lensing magnification on large-scale structure analyses combining photometric cosmic shear and galaxy clustering data (i.e. the now commonly called "3$\times$2-point" analysis). Using a Fisher matrix bias formalism, we disentangle the contribution to the bias on cosmological parameters caused by ignoring the effects of magnification in a theory fit…
▽ More
We investigate the sensitivity to the effects of lensing magnification on large-scale structure analyses combining photometric cosmic shear and galaxy clustering data (i.e. the now commonly called "3$\times$2-point" analysis). Using a Fisher matrix bias formalism, we disentangle the contribution to the bias on cosmological parameters caused by ignoring the effects of magnification in a theory fit from individual elements in the data vector, for Stage-III and Stage-IV surveys, assuming well-known redshift distributions, sample selection based on magnitude, and magnification strengths inferred from CFHTLenS and DES Y1 data. We show that the removal of elements of the data vectors that are dominated by magnification does not guarantee a reduction in the cosmological bias due to the magnification signal, but can instead increase the sensitivity to magnification. We find that the most sensitive elements of the data vector come from the shear-clustering cross-correlations, particularly between the highest-redshift shear bin and any lower-redshift lens sample, and that the parameters $Ω_{M}$, $S_8=σ_8\sqrt{Ω_M/0.3}$ and $w_0$ show the most significant biases for both survey models. Our forecasts predict that current analyses are not significantly biased by magnification, but this bias will become highly significant with the continued increase of statistical power in the near future. We therefore conclude that future surveys should measure and model the magnification as part of their flagship "3$\times$2-point" analysis.
△ Less
Submitted 7 December, 2019; v1 submitted 30 July, 2019;
originally announced July 2019.
-
An Accurate Analytic Model for the Thermal Sunyaev-Zel'dovich One-Point PDF
Authors:
Leander Thiele,
J. Colin Hill,
Kendrick M. Smith
Abstract:
Non-Gaussian statistics of late-time cosmological fields contain information beyond that captured in the power spectrum. Here we focus on one such example: the one-point probability distribution function (PDF) of the thermal Sunyaev-Zel'dovich (tSZ) signal in maps of the cosmic microwave background (CMB). It has been argued that the one-point PDF is a near-optimal statistic for cosmological constr…
▽ More
Non-Gaussian statistics of late-time cosmological fields contain information beyond that captured in the power spectrum. Here we focus on one such example: the one-point probability distribution function (PDF) of the thermal Sunyaev-Zel'dovich (tSZ) signal in maps of the cosmic microwave background (CMB). It has been argued that the one-point PDF is a near-optimal statistic for cosmological constraints from the tSZ signal, as most of the constraining power in tSZ $N$-point functions is contained in their amplitudes (rather than their shapes), which probe differently-weighted integrals over the halo mass function. In this paper, we develop a new analytic halo model for the tSZ PDF, discarding simplifying assumptions made in earlier versions of this approach. In particular, we account for effects due to overlaps of the tSZ profiles of different halos, as well as effects due to the clustering of halos. We verify the accuracy of our analytic model via comparison to numerical simulations. We demonstrate that this more accurate model is necessary for the analysis of the tSZ PDF in upcoming CMB experiments. The novel formalism developed here may be useful in modeling the one-point PDF of other cosmological observables, such as the weak lensing convergence field.
△ Less
Submitted 26 December, 2018; v1 submitted 13 December, 2018;
originally announced December 2018.