We present simulation-based cosmological $w$CDM inference using Dark Energy
Survey Year 3 weak-lensing maps, via neural data compression of weak-lensing
map summary statistics: power spectra, peak ...counts, and direct map-level
compression/inference with convolutional neural networks (CNN). Using
simulation-based inference, also known as likelihood-free or implicit
inference, we use forward-modelled mock data to estimate posterior probability
distributions of unknown parameters. This approach allows all statistical
assumptions and uncertainties to be propagated through the forward-modelled
mock data; these include sky masks, non-Gaussian shape noise, shape measurement
bias, source galaxy clustering, photometric redshift uncertainty, intrinsic
galaxy alignments, non-Gaussian density fields, neutrinos, and non-linear
summary statistics. We include a series of tests to validate our inference
results. This paper also describes the Gower Street simulation suite: 791
full-sky PKDGRAV dark matter simulations, with cosmological model parameters
sampled with a mixed active-learning strategy, from which we construct over
3000 mock DES lensing data sets. For $w$CDM inference, for which we allow
$-1<w<-\frac{1}{3}$, our most constraining result uses power spectra combined
with map-level (CNN) inference. Using gravitational lensing data only, this
map-level combination gives $\Omega_{\rm m} = 0.283^{+0.020}_{-0.027}$, ${S_8 =
0.804^{+0.025}_{-0.017}}$, and $w < -0.80$ (with a 68 per cent credible
interval); compared to the power spectrum inference, this is more than a factor
of two improvement in dark energy parameter ($\Omega_{\rm DE}, w$) precision.
Beyond-two-point statistics contain additional information on cosmological as well as astrophysical and observational (systematics) parameters. In this methodology paper we provide an end-to-end ...simulation-based analysis of a set of Gaussian and non-Gaussian weak lensing statistics using detailed mock catalogues of the Dark Energy Survey. We implement: 1) second and third moments; 2) wavelet phase harmonics (WPH); 3) the scattering transform (ST). Our analysis is fully based on simulations, it spans a space of seven \(\nu w\)CDM cosmological parameters, and it forward models the most relevant sources of systematics of the data (masks, noise variations, clustering of the sources, intrinsic alignments, and shear and redshift calibration). We implement a neural network compression of the summary statistics, and we estimate the parameter posteriors using a likelihood-free-inference approach. We validate the pipeline extensively, and we find that WPH exhibits the strongest performance when combined with second moments, followed by ST. and then by third moments. The combination of all the different statistics further enhances constraints with respect to second moments, up to 25 per cent, 15 per cent, and 90 per cent for \(S_8\), \(\Omega_{\rm m}\), and the Figure-Of-Merit \({\rm FoM_{S_8,\Omega_{\rm m}}}\), respectively. We further find that non-Gaussian statistics improve constraints on \(w\) and on the amplitude of intrinsic alignment with respect to second moments constraints. The methodological advances presented here are suitable for application to Stage IV surveys from Euclid, Rubin-LSST, and Roman with additional validation on mock catalogues for each survey. In a companion paper we present an application to DES Year 3 data.
We search for signatures of cosmological shocks in gas pressure profiles of galaxy clusters using the cluster catalogs from three surveys: the Dark Energy Survey (DES) Year 3, the South Pole ...Telescope (SPT) SZ survey, and the Atacama Cosmology Telescope (ACT) data releases 4, 5, and 6, and using thermal Sunyaev-Zeldovich (SZ) maps from SPT and ACT. The combined cluster sample contains around \(10^5\) clusters with mass and redshift ranges \(10^{13.7} < M_{\rm 200m}/M_\odot < 10^{15.5}\) and \(0.1 < z < 2\), and the total sky coverage of the maps is \(\approx 15,000 \,\,{\rm deg}^2\). We find a clear pressure deficit at \(R/R_{\rm 200m}\approx 1.1\) in SZ profiles around both ACT and SPT clusters, estimated at \(6\sigma\) significance, which is qualitatively consistent with a shock-induced thermal non-equilibrium between electrons and ions. The feature is not as clearly determined in profiles around DES clusters. We verify that measurements using SPT or ACT maps are consistent across all scales, including in the deficit feature. The SZ profiles of optically selected and SZ-selected clusters are also consistent for higher mass clusters. Those of less massive, optically selected clusters are suppressed on small scales by factors of 2-5 compared to predictions, and we discuss possible interpretations of this behavior. An oriented stacking of clusters -- where the orientation is inferred from the SZ image, the brightest cluster galaxy, or the surrounding large-scale structure measured using galaxy catalogs -- shows the normalization of the one-halo and two-halo terms vary with orientation. Finally, the location of the pressure deficit feature is statistically consistent with existing estimates of the splashback radius.
We present an alternative calibration of the MagLim lens sample redshift distributions from the Dark Energy Survey (DES) first three years of data (Y3). The new calibration is based on a combination ...of a Self-Organising Maps based scheme and clustering redshifts to estimate redshift distributions and inherent uncertainties, which is expected to be more accurate than the original DES Y3 redshift calibration of the lens sample. We describe in detail the methodology, we validate it on simulations and discuss the main effects dominating our error budget. The new calibration is in fair agreement with the fiducial DES Y3 redshift distributions calibration, with only mild differences (\(<3\sigma\)) in the means and widths of the distributions. We study the impact of this new calibration on cosmological constraints, analysing DES Y3 galaxy clustering and galaxy-galaxy lensing measurements, assuming a \(\Lambda\)CDM cosmology. We obtain \(\Omega_{\rm m} = 0.30\pm 0.04\), \(\sigma_8 = 0.81\pm 0.07 \) and \(S_8 = 0.81\pm 0.04\), which implies a \(\sim 0.4\sigma\) shift in the \(\Omega_{\rm}-S_8\) plane compared to the fiducial DES Y3 results, highlighting the importance of the redshift calibration of the lens sample in multi-probe cosmological analyses.
The kinematic and thermal Sunyaev-Zel'dovich (kSZ and tSZ) effects probe the abundance and thermodynamics of ionized gas in galaxies and clusters. We present a new hybrid estimator to measure the kSZ ...effect by combining cosmic microwave background temperature anisotropy maps with photometric and spectroscopic optical survey data. The method interpolates a velocity reconstruction from a spectroscopic catalog at the positions of objects in a photometric catalog, which makes it possible to leverage the high number density of the photometric catalog and the precision of the spectroscopic survey. Combining this hybrid kSZ estimator with a measurement of the tSZ effect simultaneously constrains the density and temperature of free electrons in the photometrically selected galaxies. Using the 1000 deg2 of overlap between the Atacama Cosmology Telescope (ACT) Data Release 5, the first three years of data from the Dark Energy Survey (DES), and the Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 12, we detect the kSZ signal at 4.8\({\sigma}\) and reject the null (no-kSZ) hypothesis at 5.1\({\sigma}\). This corresponds to 2.0\({\sigma}\) per 100,000 photometric objects with a velocity field based on a spectroscopic survey with 1/5th the density of the photometric catalog. For comparison, a recent ACT analysis using exclusively spectroscopic data from BOSS measured the kSZ signal at 2.1\({\sigma}\) per 100,000 objects. Our derived constraints on the thermodynamic properties of the galaxy halos are consistent with previous measurements. With future surveys, such as the Dark Energy Spectroscopic Instrument and the Rubin Observatory Legacy Survey of Space and Time, we expect that this hybrid estimator could result in measurements with significantly better signal-to-noise than those that rely on spectroscopic data alone.
Cross-correlation between weak lensing of the Cosmic Microwave Background (CMB) and weak lensing of galaxies offers a way to place robust constraints on cosmological and astrophysical parameters with ...reduced sensitivity to certain systematic effects affecting individual surveys. We measure the angular cross-power spectrum between the Atacama Cosmology Telescope (ACT) DR4 CMB lensing and the galaxy weak lensing measured by the Dark Energy Survey (DES) Y3 data. Our baseline analysis uses the CMB convergence map derived from ACT-DR4 and \(\textit{Planck}\) data, where most of the contamination due to the thermal Sunyaev Zel'dovich effect is removed, thus avoiding important systematics in the cross-correlation. In our modelling, we consider the nuisance parameters of the photometric uncertainty, multiplicative shear bias and intrinsic alignment of galaxies. The resulting cross-power spectrum has a signal-to-noise ratio \(= 7.1\) and passes a set of null tests. We use it to infer the amplitude of the fluctuations in the matter distribution (\(S_8 \equiv \sigma_8 (\Omega_{\rm m}/0.3)^{0.5} = 0.782\pm 0.059\)) with informative but well-motivated priors on the nuisance parameters. We also investigate the validity of these priors by significantly relaxing them and checking the consistency of the resulting posteriors, finding them consistent, albeit only with relatively weak constraints. This cross-correlation measurement will improve significantly with the new ACT-DR6 lensing map and form a key component of the joint 6x2pt analysis between DES and ACT.
We demonstrate and measure the impact of source galaxy clustering on higher-order summary statistics of weak gravitational lensing data. By comparing simulated data with galaxies that either trace or ...do not trace the underlying density field, we show this effect can exceed measurement uncertainties for common higher-order statistics for certain analysis choices. Source clustering effects are larger at small scales and for statistics applied to combinations of low and high redshift samples, and diminish at high redshift. We evaluate the impact on different weak lensing observables, finding that third moments and wavelet phase harmonics are more affected than peak count statistics. Using Dark Energy Survey Year 3 data we construct null tests for the source-clustering-free case, finding a \(p\)-value of \(p=4\times10^{-3}\) (2.6 \(\sigma\)) using third-order map moments and \(p=3\times10^{-11}\) (6.5 \(\sigma\)) using wavelet phase harmonics. The impact of source clustering on cosmological inference can be either be included in the model or minimized through \textit{ad-hoc} procedures (e.g. scale cuts). We verify that the procedures adopted in existing DES Y3 cosmological analyses (using map moments and peaks) were sufficient to render this effect negligible. Failing to account for source clustering can significantly impact cosmological inference from higher-order gravitational lensing statistics, e.g. higher-order N-point functions, wavelet-moment observables (including phase harmonics and scattering transforms), and deep learning or field level summary statistics of weak lensing maps. We provide recipes both to minimise the impact of source clustering and to incorporate source clustering effects into forward-modelled mock data.
Widefield surveys of the sky probe many clustered scalar fields -- such as galaxy counts, lensing potential, gas pressure, etc. -- that are sensitive to different cosmological and astrophysical ...processes. Our ability to constrain such processes from these fields depends crucially on the statistics chosen to summarize the field. In this work, we explore the cumulative distribution function (CDF) at multiple scales as a summary of the galaxy lensing convergence field. Using a suite of N-body lightcone simulations, we show the CDFs' constraining power is modestly better than that of the 2nd and 3rd moments of the field, as they approximately capture the information from all moments of the field in a concise data vector. We then study the practical aspects of applying the CDFs to observational data, using the first three years of the Dark Energy Survey (DES Y3) data as an example, and compute the impact of different systematics on the CDFs. The contributions from the point spread function are 2-3 orders of magnitude below the cosmological signal, while those from reduced shear approximation contribute \(\lesssim 1\%\) to the signal. Source clustering effects and baryon imprints contribute \(1-10\%\). Enforcing scale cuts to limit systematics-driven biases in parameter constraints degrades these constraints a noticeable amount, and this degradation is similar for the CDFs and the moments. We also detect correlations between the observed convergence field and the shape noise field at \(13\sigma\). We find that the non-Gaussian correlations in the noise field must be modeled accurately to use the CDFs, or other statistics sensitive to all moments, as a rigorous cosmology tool.
We present a measurement of the cross-correlation between the MagLim galaxies selected from the Dark Energy Survey (DES) first three years of observations (Y3) and cosmic microwave background (CMB) ...lensing from the Atacama Cosmology Telescope (ACT) Data Release 4 (DR4), reconstructed over \(\sim 436\) sq.deg. of the sky. Our galaxy sample, which covers \(\sim 4143\) sq.deg., is divided into six redshift bins spanning the redshift range of \(0.20<z<1.05\). We adopt a blinding procedure until passing all consistency and systematics tests. After imposing scale cuts for the cross-power spectrum measurement, we reject the null hypothesis of no correlation at 9.1\sigma. We constrain cosmological parameters from a joint analysis of galaxy and CMB lensing-galaxy power spectra considering a flat \LCDM model, marginalized over 23 astrophysical and systematic nuisance parameters. We find the clustering amplitude \(S_8\equiv \sigma_8 (\Omega_m/0.3)^{0.5} = 0.75^{+0.04}_{-0.05}\). In addition, we constrain the linear growth of cosmic structure as a function of redshift. Our results are consistent with recent DES Y3 analyses and suggest a preference for a lower \(S_8\) compared to results from measurements of CMB anisotropies by the Planck satellite, although at a mild level (\(< 2 \sigma\)) of statistical significance.
We study the effect of magnification in the Dark Energy Survey Year 3 analysis of galaxy clustering and galaxy-galaxy lensing, using two different lens samples: a sample of Luminous red galaxies, ...redMaGiC, and a sample with a redshift-dependent magnitude limit, MagLim. We account for the effect of magnification on both the flux and size selection of galaxies, accounting for systematic effects using the Balrog image simulations. We estimate the impact of magnification on the galaxy clustering and galaxy-galaxy lensing cosmology analysis, finding it to be a significant systematic for the MagLim sample. We show cosmological constraints from the galaxy clustering auto-correlation and galaxy-galaxy lensing signal with different magnifications priors, finding broad consistency in cosmological parameters in \(\Lambda\)CDM and \(w\)CDM. However, when magnification bias amplitude is allowed to be free, we find the two-point correlations functions prefer a different amplitude to the fiducial input derived from the image simulations. We validate the magnification analysis by comparing the cross-clustering between lens bins with the prediction from the baseline analysis, which uses only the auto-correlation of the lens bins, indicating systematics other than magnification may be the cause of the discrepancy. We show adding the cross-clustering between lens redshift bins to the fit significantly improves the constraints on lens magnification parameters and allows uninformative priors to be used on magnification coefficients, without any loss of constraining power or prior volume concerns.