A typical Bayesian inference on the values of some parameters of interest q from some data D involves running a Markov Chain (MC) to sample from the posterior p(q,n divided by D)proportional to L(D divided by q,n)p(q)p(n), where n are some nuisance parameters with a separable prior. In some cases, the nuisance parameters are high-dimensional, and their prior p(n) is itself defined only by a set of samples that have been drawn from some other MC. The MC for the posterior will typically require evaluation of p(n) at arbitrary values of n, i.e., one needs to provide a density estimator over the full n space from the provided samples. But the high dimensionality of n hinders both the density estimation and the efficiency of the MC for the posterior. We describe a solution to this problem: a linear compression of the n space into a much lower-dimensional space u, which projects away directions in n space that cannot appreciably alter L. The algorithm for doing so is a slight modification to principal components analysis, and is less restrictive on p(n) than other proposed solutions to this issue. We demonstrate this "mode projection" technique using the analysis of 2-point correlation functions of weak lensing fields and galaxy density in the Dark Energy Survey, where n is a binned representation of the redshift distribution n(z) of the galaxies.

Characterization of the redshift distribution of ensembles of galaxies is pivotal for large scale structure cosmological studies. In this work, we focus on improving the Self-Organizing Map (SOM) methodology for photometric redshift estimation (SOMPZ), specifically in anticipation of the Dark Energy Survey Year 6 (DES Y6) data. This data set, featuring deeper and fainter galaxies than DES Year 3 (DES Y3), demands adapted techniques to ensure accurate recovery of the underlying redshift distribution. We investigate three strategies for enhancing the existing SOM-based approach used in DES Y3: 1) Replacing the Y3 SOM algorithm with one tailored for redshift estimation challenges; 2) Incorporating $\textit{g}$-band flux information to refine redshift estimates (i.e. using $\textit{griz}$ fluxes as opposed to only $\textit{riz}$); 3) Augmenting redshift data for galaxies where available. These methods are applied to DES Y3 data, and results are compared to the Y3 fiducial ones. Our analysis indicates significant improvements with the first two strategies, notably reducing the overlap between redshift bins. By combining strategies 1 and 2, we have successfully managed to reduce redshift bin overlap in DES Y3 by up to 66$\%$. Conversely, the third strategy, involving the addition of redshift data for selected galaxies as an additional feature in the method, yields inferior results and is abandoned. Our findings contribute to the advancement of weak lensing redshift characterization and lay the groundwork for better redshift characterization in DES Year 6 and future stage IV surveys, like the Rubin Observatory.

The orientation of triaxial galaxy clusters with respect to the line-of-sight is expected to be one of the prime sources of scatter and potential bias in optical observables (e.g., richness and weak-lensing signal) of galaxy clusters. In this work, we use the observed shape of the central Brightest Cluster Galaxy (BCG) as proxy for the orientation along the line-of-sight for clusters selected via the Sunyaev-Zel'dovich (SZ) effect from the South Pole Telescope (SPT) and Atacama Cosmology Telescope (ACT) surveys, matched to optically selected clusters from the Dark Energy Survey Year 3 (DES). We construct two samples of clusters that are designed to be identical in SZ mass estimate and redshift but with the roundest vs. the most elliptical BCGs, which we expect to correspond to BCGs (and clusters) with major axes aligned along the line-of-sight vs. in the plane of the sky, respectively. We find that the optical richness of round-BCG clusters is ∼ 10% larger than that of elliptical-BCG clusters, in agreement with the expectation from projection effects and presenting the first such detection in data. The density profiles, however, are not in agreement with the expectation from projection effects: the 1-halo term (below 6 h^-1Mpc) of both the weak-lensing and galaxy density profiles are the same for the subsamples, contrary to previous studies based on X-ray selected clusters. In the 2-halo regime (above 6 h^-1Mpc), we find a significant excess of the elliptical-BCG cluster profiles compared to the round-BCG cluster profiles, which is the opposite of the expectation from numerical simulations. We hypothesize that the intrinsic shape of the BCG reflects not just the orientation angle, but also intrinsic properties of the cluster which can affect both the SZ signal and the amplitude of the 2-halo term.

NASA’s Nancy Grace Roman Space Telescope (Roman) will provide an opportunity to study dark energy with unprecedented precision and accuracy using several techniques, including measurements of high- z Type Ia Supernovae (SNe Ia; z ≲ 3.0) via the High-Latitude Time Domain Survey (HLTDS). In this work, we do an initial “benchmark” characterization of the photometric repeatability of stellar fluxes, which must be below 1% when sky noise is subdominant in order to enable a number of calibration requirements. Achieving this level of flux precision requires attention to Roman’s highly structured, spatially varying, undersampled point-spread function (PSF). In this work, we build a library of effective PSFs (ePSFs) compatible with the OpenUniverse HLTDS simulations. Using our library of ePSFs, we recover stellar flux to between 0.6% and 1.2% photometric precision, finding that redder bands perform better by this metric. We also find that flux recovery is improved by up to 20% when a chip (sensor chip assembly, SCA) is divided into eight sub-SCAs in order to account for the spatial variation of the PSF. With our optimized algorithm, we measure nonlinearity due to photometry (magnitude dependence) of ∣ s _NL ∣ < 1.93 × 10 ^−3 dex ^−1 , which is still larger than stated requirements for detector effects and indicates that further work is necessary. We also measure the dependence of photometric residuals on stellar color and find the largest possible dependence in R062, implying a color-dependent PSF model may be needed. Finally, we characterize the detection efficiency function of each OpenUniverse Roman filter, which will inform future studies.

Cosmological constraints using weak gravitational lensing measurements from the Roman Space Telescope will require a powerful method for modelling uncertainties in the galaxy redshift distribution. In this work, we use an optimized version of the principal component analysis (PCA) to model uncertainties in the full shape of the redshift distributions, a method proposed by and recently used in the Dark Energy Survey Y6 analysis. Here, we implement this new approach within the Roman High Latitude Imaging Survey (HLIS) Cosmology Project Infrastructure Team (PIT) pipeline, namely Cobaya-Cosmolike Joint Architecture (). To validate the PCA in mitigating biases on cosmological parameters, S_8 and Ω_m, we use a set of redshift distributions from generated for a variety of Roman configurations. Overall, when the simulated cosmic shear data vector is not strongly miscalibrated relative to the fiducial one, both the mean-shift and the PCA-based approaches produce consistent cosmological constraints when marginalizing over nuisance parameters. For mild to strong miscalibration, including additional PCs progressively mitigates biases in S_8 and Ω_m, and can achieve comparable performance with fewer parameters than the nine tomographic-bin mean-shift model.

Modeling the intrinsic alignment (IA) of galaxies poses a challenge to weak lensing analyses. Using the Dark Energy Survey Year 3 shape catalog, we expect less impact from IA when we limit the sample to blue, star-forming galaxies. The cosmological parameter constraints from this BLUE cosmic shear sample are stable to IA model choice, unlike passive galaxies in the full DES Y3 sample, the goodness-of-fit is improved and the Omega(m) and S-8 better agree with the observations from Planck on the cosmic microwave background. Mitigating IA with sample selection in DES, rather than flexible model choices, can reduce uncertainty in S-8 by a factor of 1.5.

Cosmic shear, galaxy clustering, and the abundance of massive halos each probe the large-scale structure of the Universe in complementary ways. We present cosmological constraints from the joint analysis of the three probes, building on the latest analyses of the lensing-informed abundance of clusters identified by the South Pole Telescope (SPT) and of the auto- and cross-correlation of galaxy position and weak lensing measurements (3 x 2pt) in the Dark Energy Survey (DES). We consider the cosmological correlation between the different tracers and we account for the systematic uncertainties that are shared between the large-scale lensing correlation functions and the small-scale lensing-based cluster mass calibration. Marginalized over the remaining Lambda cold dark matter (Lambda CDM) parameters (including the sum of neutrino masses) and 52 astrophysical modeling parameters, we measure Omega(m) = 0.300 +/- 0.017 and sigma(8) = 0.797 +/- 0.026. Compared to constraints from Planck primary cosmic microwave background (CMB) anisotropies, our constraints are only 15% wider with a probability to exceed of 0.22 (1.2 sigma) for the two-parameter difference. We further obtain S-8 equivalent to sigma(8)(Omega(m)/0.3)(0.5) = 0.796 +/- 0.013 which is lower than the Planck measurement at the 1.6 sigma level. The combined SPT cluster, DES 3 x 2pt, and Planck datasets mildly prefer a nonzero positive neutrino mass, with a 95% upper limit Sigma m(nu) < 0.25 eV on the sum of neutrino masses. Assuming a wCDM model, we constrain the dark energy equation of state parameter w = -1.15(-0.17)(+0.23) and when combining with Planck primary CMB anisotropies, we recover w = -1.20(-0.09)(+0.15), a 1.7 sigma difference with a cosmological constant. The precision of our results highlights the benefits of multiwavelength multiprobe cosmology and our analysis paves the way for upcoming joint analyses of next-generation datasets.

NASA's Nancy Grace Roman Space Telescope (Roman) will provide an opportunity to study dark energy with unprecedented precision using several techniques, including measurements of Type Ia Supernovae (SNe Ia). Here, we present `phrosty` (PHotometry for ROman with SFFT for tYpe Ia supernovae): a difference imaging pipeline for measuring the brightness of transient point sources in the sky, primarily SNe Ia, using Roman data. `phrosty` is written in Python. We implement a GPU-accelerated version of the Saccadic Fast Fourier Transform (SFFT) method for difference imaging.

The Nancy Grace Roman Space Telescope, set to launch in 2026, will bring unprecedented precision to measurements of weak gravitational lensing. Because weak lensing is an inherently small signal, it is imperative to minimize systematic errors in measurements as completely as possible; this will ensure that the lensing measurements can be used to their full potential when extracting cosmological information. In this paper, we use laboratory tests of the Roman detectors, simulations of the Roman High Latitude Survey observations, and the proposed Roman image combination pipeline to investigate the magnitude of detector read noise biasing on weak lensing measurements from Roman. First, we combine lab-measured detector noise fields with simulated observations and propagate these images through the Roman image combination pipeline, Imcom. We characterize the specific signatures of the noise fields in the resultant images and find that noise contributes to the combined images most strongly at scales relevant to physical characteristics of the detector including PSF shape, chip boundaries, and roll angles. We then measure shapes of simulated stars and galaxies and determine the magnitude of noise-induced shear bias on these measurements. We find that star shape correlations satisfy the system noise requirements as defined by the Roman Science Requirements Document. However, for galaxies fainter than mAB similar or equal to 24, correction for noise correlations will be needed in order to ensure confidence in shape measurements in any observation band.

We present simulation-based cosmological wcold dark matter (wCDM) 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 pkdgrav3 dark matter simulations, with cosmological model parameters sampled with a mixed active-learning strategy, from which we construct over 3000 mock dark energy survey lensing data sets. For wCDM inference, for which we allow -1<w<-(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(m)=0.283(-0.027)(+0.020), S-8=0.804(-0.017)(+0.025), 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(DE),w) precision.

ABSTRACT We 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 that this effect can exceed measurement uncertainties for common higher order statistics for certain analysis choices. 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 (DES) Year 3 (Y3) data, we construct null tests for the source-clustering-free case, finding a p-value of p = 4 × 10−3 (2.6σ) using third-order map moments and p = 3 × 10−11 (6.5σ) using wavelet phase harmonics. The impact of source clustering on cosmological inference can be either included in the model or minimized through ad hoc procedures (e.g. scale cuts). We verify that the procedures adopted in existing DES Y3 cosmological analyses 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, and deep learning or field-level summary statistics of weak lensing maps.

We present cosmological constraints from the sample of Type Ia supernovae (SN Ia) discovered during the full five years of the Dark Energy Survey (DES) Supernova Program. In contrast to most previous cosmological samples, in which SN are classified based on their spectra, we classify the DES SNe using a machine learning algorithm applied to their light curves in four photometric bands. Spectroscopic redshifts are acquired from a dedicated follow-up survey of the host galaxies. After accounting for the likelihood of each SN being a SN Ia, we find 1635 DES SNe in the redshift range $0.100.5$ SNe compared to the previous leading compilation of Pantheon+, and results in the tightest cosmological constraints achieved by any SN data set to date. To derive cosmological constraints we combine the DES supernova data with a high-quality external low-redshift sample consisting of 194 SNe Ia spanning $0.025

We present galaxy-galaxy lensing measurements using a sample of low surface brightness galaxies (LSBGs) drawn from the Dark Energy Survey Year 3 (Y3) data as lenses. LSBGs are diffuse galaxies with a surface brightness dimmer than the ambient night sky. These dark-matter-dominated objects are intriguing due to potentially unusual formation channels that lead to their diffuse stellar component. Given the faintness of LSBGs, using standard observational techniques to characterize their total masses proves challenging. Weak gravitational lensing, which is less sensitive to the stellar component of galaxies, could be a promising avenue to estimate the masses of LSBGs. Our LSBG sample consists of 23,790 galaxies separated into red and blue color types at $g-i\ge 0.60$ and $g-i< 0.60$, respectively. Combined with the DES Y3 shear catalog, we measure the tangential shear around these LSBGs and find signal-to-noise ratios of 6.67 for the red sample, 2.17 for the blue sample, and 5.30 for the full sample. We use the clustering redshifts method to obtain redshift distributions for the red and blue LSBG samples. Assuming all red LSBGs are satellites, we fit a simple model to the measurements and estimate the host halo mass of these LSBGs to be $\log(M_{\rm host}/M_{\odot}) = 12.98 ^{+0.10}_{-0.11}$. We place a 95% upper bound on the subhalo mass at $\log(M_{\rm sub}/M_{\odot})<11.51$. By contrast, we assume the blue LSBGs are centrals, and place a 95% upper bound on the halo mass at $\log(M_\mathrm{host}/M_\odot) < 11.84$. We find that the stellar-to-halo mass ratio of the LSBG samples is consistent with that of the general galaxy population. This work illustrates the viability of using weak gravitational lensing to constrain the halo masses of LSBGs.

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.

One challenge for applying current weak lensing analysis tools to the Nancy Grace Roman Space Telescope is that individual images will be undersampled. Our companion paper presented an initial application of Imcom - an algorithm that builds an optimal mapping from input to output pixels to reconstruct a fully sampled combined image - on the Roman image simulations. In this paper, we measure the output noise power spectra, identify the sources of the major features in the power spectra, and show that simple analytic models that ignore sampling effects underestimate the power spectra of the coadded noise images. We compute the moments of both idealized injected stars and fully simulated stars in the coadded images, and their 1- and 2-point statistics. We show that the idealized injected stars have root-mean-square ellipticity errors (1 - 6) x 10-4 per component depending on the band; the correlation functions are >= 2 orders of magnitude below requirements, indicating that the image combination step itself is using a small fraction of the overall Roman 2nd moment error budget, although the 4th moments are larger and warrant further investigation. The stars in the simulated sky images, which include blending and chromaticity effects, have correlation functions near the requirement level (and below the requirement level in a wide-band image constructed by stacking all 4 filters). We evaluate the noise-induced biases in the ellipticities of injected stars, and explain the resulting trends with an analytical model. We conclude by enumerating the next steps in developing an image coaddition pipeline for Roman.

We present the angular diameter distance measurement obtained with the Baryonic Acoustic Oscillation feature from galaxy clustering in the completed Dark Energy Survey, consisting of six years (Y6) of observations. We use the Y6 BAO galaxy sample, optimized for BAO science in the redshift range 0.6<z<1.2, with an effective redshift at z_ eff=0.85 and split into six tomographic bins. The sample has nearly 16 million galaxies over 4,273 square degrees. Our consensus measurement constrains the ratio of the angular distance to sound horizon scale to D_M(z_ eff)/r_d = 19.51±0.41 (at 68.3

We present cosmological constraints from the abundance of galaxy clusters selected via the thermal Sunyaev-Zel'dovich (SZ) effect in South Pole Telescope (SPT) data with a simultaneous mass calibration using weak gravitational lensing data from the Dark Energy Survey (DES) and the Hubble Space Telescope (HST). The cluster sample is constructed from the combined SPT-SZ, SPTpol ECS, and SPTpol 500d surveys, and comprises 1,005 confirmed clusters in the redshift range $0.25-1.78$ over a total sky area of 5,200 deg$^2$. We use DES Year 3 weak-lensing data for 688 clusters with redshifts $z<0.95$ and HST weak-lensing data for 39 clusters with $0.6<z<1.7$. The weak-lensing measurements enable robust mass measurements of sample clusters and allow us to empirically constrain the SZ observable--mass relation. For a flat $\Lambda$CDM cosmology, and marginalizing over the sum of massive neutrinos, we measure $\Omega_\mathrm{m}=0.286\pm0.032$, $\sigma_8=0.817\pm0.026$, and the parameter combination $\sigma_8\,(\Omega_\mathrm{m}/0.3)^{0.25}=0.805\pm0.016$. Our measurement of $S_8\equiv\sigma_8\,\sqrt{\Omega_\mathrm{m}/0.3}=0.795\pm0.029$ and the constraint from Planck CMB anisotropies (2018 TT,TE,EE+lowE) differ by $1.1\sigma$. In combination with that Planck dataset, we place a 95% upper limit on the sum of neutrino masses $\sum m_\nu<0.18$ eV. When additionally allowing the dark energy equation of state parameter $w$ to vary, we obtain $w=-1.45\pm0.31$ from our cluster-based analysis. In combination with Planck data, we measure $w=-1.34^{+0.22}_{-0.15}$, or a $2.2\sigma$ difference with a cosmological constant. We use the cluster abundance to measure $\sigma_8$ in five redshift bins between 0.25 and 1.8, and we find the results to be consistent with structure growth as predicted by the $\Lambda$CDM model fit to Planck primary CMB data.

We present a simulation-based cosmological analysis using a combination of Gaussian and non-Gaussian statistics of the weak lensing mass (convergence) maps from the first three years (Y3) of the Dark Energy Survey (DES). We implement: 1) second and third moments; 2) wavelet phase harmonics; 3) the scattering transform. Our analysis is fully based on simulations, spans a space of seven $\nu w$CDM cosmological parameters, and forward models the most relevant sources of systematics inherent in 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 simulation-based inference approach. Including and combining different non-Gaussian statistics is a powerful tool that strongly improves constraints over Gaussian statistics (in our case, the second moments); in particular, the Figure of Merit $\textrm{FoM}(S_8, \Omega_{\textrm{m}})$ is improved by 70 percent ($\Lambda$CDM) and 90 percent ($w$CDM). When all the summary statistics are combined, we achieve a 2 percent constraint on the amplitude of fluctuations parameter $S_8 \equiv \sigma_8 (\Omega_{\textrm{m}}/0.3)^{0.5}$, obtaining $S_8 = 0.794 \pm 0.017$ ($\Lambda$CDM) and $S_8 = 0.817 \pm 0.021$ ($w$CDM). The constraints from different statistics are shown to be internally consistent (with a $p$-value>0.1 for all combinations of statistics examined). We compare our results to other weak lensing results from the DES Y3 data, finding good consistency; we also compare with results from external datasets, such as \planck{} constraints from the Cosmic Microwave Background, finding statistical agreement, with discrepancies no greater than $<2.2\sigma$.

We present a Bayesian population modeling method to analyze the abundance of galaxy clusters identified by the South Pole Telescope (SPT) with a simultaneous mass calibration using weak gravitational lensing data from the Dark Energy Survey (DES) and the Hubble Space Telescope (HST). We discuss and validate the modeling choices with a particular focus on a robust, weak-lensing-based mass calibration using DES data. For the DES Year 3 data, we report a systematic uncertainty in weak-lensing mass calibration that increases from 1% at z = 0.25 to 10% at z = 0.95, to which we add 2% in quadrature to account for uncertainties in the impact of baryonic effects. We implement an analysis pipeline that joins the cluster abundance likelihood with a multiobservable likelihood for the Sunyaev-Zel'dovich effect, optical richness, and weak-lensing measurements for each individual cluster. We validate that our analysis pipeline can recover unbiased cosmological constraints by analyzing mocks that closely resemble the cluster sample extracted from the SPT-SZ, SPTpol ECS, and SPTpol 500d surveys and the DES Year 3 and HST-39 weak-lensing datasets. This work represents a crucial prerequisite for the subsequent cosmological analysis of the real dataset.

The upcoming Nancy Grace Roman Space Telescope will carry out a wide-area survey in the near infrared. A key science objective is the measurement of cosmic structure via weak gravitational lensing. Roman data will be undersampled, which introduces new challenges in the measurement of source galaxy shapes; a potential solution is to use linear algebra-based coaddition techniques such as Imcom that combine multiple undersampled images to produce a single oversampled output mosaic with a desired "target" point spread function (PSF). We present here an initial application of Imcom to 0.64 square degrees of simulated Roman data, based on the Roman branch of the Legacy Survey of Space and Time (LSST) Dark Energy Science Collaboration (DESC) Data Challenge 2 (DC2) simulation. We show that Imcom runs successfully on simulated data that includes features such as plate scale distortions, chip gaps, detector defects, and cosmic ray masks. We simultaneously propagate grids of injected sources and simulated noise fields as well as the full simulation. We quantify the residual deviations of the PSF from the target (the "leakage"), as well as noise properties of the output images; we discuss how the overall tiling pattern as well as Moiré patterns appear in the final leakage and noise maps. We include appendices on interpolation algorithms and the interaction of undersampling with image processing operations that may be of broader applicability. The companion paper ("Paper II") explores the implications for weak lensing analyses.