Search

Modeling the intrinsic alignment (IA) of galaxies poses a challenge to weak lensing analyses. The Dark Energy Survey is expected to be less impacted by IA when limited 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 cosmic microwave background. Mitigating IA with sample selection, instead of flexible model choices, can reduce uncertainty in $S_8$ by a factor of 1.5., Comment: Data access available at https://jamiemccullough.github.io/data/blueshear/

The number density and redshift evolution of optically selected galaxy clusters offer an independent measurement of the amplitude of matter fluctuations, $S_8$. However, recent results have shown that clusters chosen by the redMaPPer algorithm show richness-dependent biases that affect the weak lensing signals and number densities of clusters, increasing uncertainty in the cluster mass calibration and reducing their constraining power. In this work, we evaluate an alternative cluster proxy, outskirt stellar mass, $M_{\textrm{out}}$, defined as the total stellar mass within a $[50,100]$ kpc envelope centered on a massive galaxy. This proxy exhibits scatter comparable to redMaPPer richness, $\lambda$, but is less likely to be subject to projection effects. We compare the Dark Energy Survey Year 3 redMaPPer cluster catalog with a $M_{\textrm{out}}$ selected cluster sample from the Hyper-Suprime Camera survey. We use weak lensing measurements to quantify and compare the scatter of $M_{\textrm{out}}$ and $\lambda$ with halo mass. Our results show $M_{\textrm{out}}$ has a scatter consistent with $\lambda$, with a similar halo mass dependence, and that both proxies contain unique information about the underlying halo mass. We find $\lambda$-selected samples introduce features into the measured $\Delta \Sigma$ signal that are not well fit by a log-normal scatter only model, absent in $M_{\textrm{out}}$ selected samples. Our findings suggest that $M_{\textrm{out}}$ offers an alternative for cluster selection with more easily calibrated selection biases, at least at the generally lower richnesses probed here. Combining both proxies may yield a mass proxy with a lower scatter and more tractable selection biases, enabling the use of lower mass clusters in cosmology. Finally, we find the scatter and slope in the $\lambda-M_{\textrm{out}}$ scaling relation to be $0.49 \pm 0.02$ and $0.38 \pm 0.09$., Comment: 22 pages, 8 figures, 4 tables, submitted to PRD

Gravitational lensing magnification of Type Ia supernovae (SNe Ia) allows information to be obtained about the distribution of matter on small scales. In this paper, we derive limits on the fraction $\alpha$ of the total matter density in compact objects (which comprise stars, stellar remnants, small stellar groupings and primordial black holes) of mass $M > 0.03 M_{\odot}$ over cosmological distances. Using 1,532 SNe Ia from the Dark Energy Survey Year 5 sample (DES-SN5YR) combined with a Bayesian prior for the absolute magnitude $M$, we obtain $\alpha < 0.12$ at the 95\% confidence level after marginalisation over cosmological parameters, lensing due to large-scale structure, and intrinsic non-Gaussianity. Similar results are obtained using priors from the cosmic microwave background, baryon acoustic oscillations and galaxy weak lensing, indicating our results do not depend on the background cosmology. We argue our constraints are likely to be conservative (in the sense of the values we quote being higher than the truth), but discuss scenarios in which they could be weakened by systematics of the order of $\Delta \alpha \sim 0.04$, Comment: Submitted to MNRAS

We present constraints on the $f(R)$ gravity model using a sample of 1,005 galaxy clusters in the redshift range $0.25 - 1.78$ that have been selected through the thermal Sunyaev-Zel'dovich effect (tSZE) from South Pole Telescope (SPT) data and subjected to optical and near-infrared confirmation with the Multi-component Matched Filter (MCMF) algorithm. We employ weak gravitational lensing mass calibration from the Dark Energy Survey (DES) Year 3 data for 688 clusters at $z < 0.95$ and from the Hubble Space Telescope (HST) for 39 clusters with $0.6 < z < 1.7$. Our cluster sample is a powerful probe of $f(R)$ gravity, because this model predicts a scale-dependent enhancement in the growth of structure, which impacts the halo mass function (HMF) at cluster mass scales. To account for these modified gravity effects on the HMF, our analysis employs a semi-analytical approach calibrated with numerical simulations. Combining calibrated cluster counts with primary cosmic microwave background (CMB) temperature and polarization anisotropy measurements from the Planck2018 release, we derive robust constraints on the $f(R)$ parameter $f_{R0}$. Our results, $\log_{10} |f_{R0}| < -5.32$ at the 95 % credible level, are the tightest current constraints on $f(R)$ gravity from cosmological scales. This upper limit rules out $f(R)$-like deviations from general relativity that result in more than a $\sim$20 % enhancement of the cluster population on mass scales $M_\mathrm{200c}>3\times10^{14}M_\odot$., Comment: 21 pages, 6 figures, submitted to Phys. Rev. D

Statistics that capture the directional dependence of the baryon distribution in the cosmic web enable unique tests of cosmology and astrophysical feedback. We use constrained oriented stacking of thermal Sunyaev-Zel'dovich (tSZ) maps to measure the anisotropic distribution of hot gas $2.5-40$ Mpc away from galaxy clusters embedded in massive filaments and superclusters. The cluster selection and orientation (at a scale of $\sim15$ Mpc) use Dark Energy Survey (DES) Year 3 data, while expanded tSZ maps from the Atacama Cosmology Telescope Data Release 6 enable a $\sim3\times$ more significant measurement of the extended gas compared to the technique's proof-of-concept. Decomposing stacks into cosine multipoles of order $m$, we detect a dipole ($m=1$) and quadrupole ($m=2$) at $8-10\sigma$, as well as evidence for $m=4$ signal at up to $6\sigma$, indicating sensitivity to late-time non-Gaussianity. We compare to the Cardinal simulations with spherical gas models pasted onto dark matter halos. The fiducial tSZ data can discriminate between two models that deplete pressure differently in low-mass halos (mimicking astrophysical feedback), preferring higher average pressure in extended structures. However, uncertainty in the amount of cosmic infrared background contamination reduces the constraining power. Additionally, we apply the technique to DES galaxy density and weak lensing to study for the first time their oriented relationships with tSZ. In the tSZ-to-lensing relation, averaged on 7.5 Mpc (transverse) scales, we observe dependence on redshift but not shape or radial distance. Thus, on large scales, the superclustering of gas pressure, galaxies, and total matter is coherent in shape and extent., Comment: 45 pages, 18 figures, submitted to ApJ

Using 1533 type Ia supernovae (SNe Ia) from the five-year sample of the Dark Energy Survey (DES), we investigate the effects of projected galactocentric separation between the SNe and their host galaxies on their light curves and standardization. We show, for the first time, that the difference in SN Ia post-standardization brightnesses between high and low-mass hosts reduces from $0.078\pm0.011$ mag in the full sample to $0.036 \pm 0.018$ mag for SNe Ia located in the outer regions of their host galaxies, while increasing to $0.100 \pm 0.014$ mag for SNe in the inner regions. In these inner regions, the step can be reduced (but not removed) using a model where the $R_V$ of dust along the line-of-sight to the SN changes as a function of galaxy properties. To explain the remaining difference, we use the distributions of the SN Ia stretch parameter to test whether the inferred age of SN progenitors are more varied in the inner regions of galaxies. We find that the proportion of high-stretch SNe Ia in red (older) environments is more prominent in outer regions and that the outer regions stretch distributions are overall more homogeneous compared to inner regions, but conclude that this effect cannot explain the reduction in significance of any Hubble residual step in outer regions. We conclude that the standardized distances of SNe Ia located in the outer regions of galaxies are less affected by their global host galaxy properties than those in the inner regions., Comment: 17 pages, 13 figures

The present state of cosmology is facing a crisis where there is a fundamental disagreement in measurements of the Hubble constant ($H_{0}$), with significant tension between the early and late universe methods. Type Ia supernovae (SNe Ia) are important to measuring $H_{0}$ through the astronomical distance ladder. However, there remains potential to better standardise SN Ia light curves by using known dependencies on host galaxy properties after the standard light curve width and colour corrections have been applied to the peak SN Ia luminosities. To explore this, we use the 5-year photometrically identified SNe Ia sample obtained by the Dark Energy Survey, along with host galaxy spectra obtained by the Australian Dark Energy Survey. Using host galaxy spectroscopy, we find a significant trend with the equivalent width (EW) of the [OII] $\lambda\lambda$ 3727, 29 doublet, a proxy for specific star formation rate, and Hubble residuals. We find that the correlation with [OII] EW is a powerful alternative to the commonly used mass step after initial light curve corrections. We applied our [OII] EW correction to a sample of 20 SN Ia hosted by calibrator galaxies observed using WiFeS, and examined the impact on both the SN Ia absolute magnitude and $H_{0}$. We then explored different [OII] EW corrections and found $H_{0}$ values ranging between $72.80$ to $73.28~\mathrm{km} \mathrm{s}^{-1} \mathrm{Mpc}^{-1}$. Notably, even after using an additional [OII] EW correction, the impact of host galaxy properties in standardising SNe Ia appears limited in reducing the current tension ($\sim$5$\sigma$) with the Cosmic Microwave Background result for $H_{0}$., Comment: 16 pages, 13 figures. Submitting to MNRAS

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.

Cosmological analyses with Type Ia Supernovae (SNe Ia) have traditionally been reliant on spectroscopy for both classifying the type of supernova and obtaining reliable redshifts to measure the distance-redshift relation. While obtaining a host-galaxy spectroscopic redshift for most SNe is feasible for small-area transient surveys, it will be too resource intensive for upcoming large-area surveys such as the Vera Rubin Observatory Legacy Survey of Space and Time, which will observe on the order of millions of SNe. Here we use data from the Dark Energy Survey (DES) to address this problem with photometric redshifts (photo-z) inferred directly from the SN light-curve in combination with Gaussian and full p(z) priors from host-galaxy photo-z estimates. Using the DES 5-year photometrically-classified SN sample, we consider several photo-z algorithms as host-galaxy photo-z priors, including the Self-Organizing Map redshifts (SOMPZ), Bayesian Photometric Redshifts (BPZ), and Directional-Neighbourhood Fitting (DNF) redshift estimates employed in the DES 3x2 point analyses. With detailed catalog-level simulations of the DES 5-year sample, we find that the simulated w can be recovered within $\pm$0.02 when using SN+SOMPZ or DNF prior photo-z, smaller than the average statistical uncertainty for these samples of 0.03. With data, we obtain biases in w consistent with simulations within ~1$\sigma$ for three of the five photo-z variants. We further evaluate how photo-z systematics interplay with photometric classification and find classification introduces a subdominant systematic component. This work lays the foundation for next-generation fully photometric SNe Ia cosmological analyses., Comment: 19 pages, 9 figures. Submitting to MNRAS, comments welcome

We present a study of the weak lensing matter profiles of 698 South Pole Telescope (SPT) thermal Sunyaev-Zel'dovich effect (tSZE) selected galaxy clusters in the redshift range $0.25

We present the joint tomographic analysis of galaxy-galaxy lensing and galaxy clustering in harmonic space, using galaxy catalogues from the first three years of observations by the Dark Energy Survey (DES Y3). We utilise the redMaGiC and MagLim catalogues as lens galaxies and the METACALIBRATION catalogue as source galaxies. The measurements of angular power spectra are performed using the pseudo-$C_\ell$ method, and our theoretical modelling follows the fiducial analyses performed by DES Y3 in configuration space, accounting for galaxy bias, intrinsic alignments, magnification bias, shear magnification bias and photometric redshift uncertainties. We explore different approaches for scale cuts based on non-linear galaxy bias and baryonic effects contamination. Our fiducial covariance matrix is computed analytically, accounting for mask geometry in the Gaussian term, and including non-Gaussian contributions and super-sample covariance terms. To validate our harmonic space pipelines and covariance matrix, we used a suite of 1800 log-normal simulations. We also perform a series of stress tests to gauge the robustness of our harmonic space analysis. In the $\Lambda$CDM model, the clustering amplitude $S_8 =\sigma_8(\Omega_m/0.3)^{0.5}$ is constrained to $S_8 = 0.704\pm 0.029$ and $S_8 = 0.753\pm 0.024$ ($68\%$ C.L.) for the redMaGiC and MagLim catalogues, respectively. For the $w$CDM, the dark energy equation of state is constrained to $w = -1.28 \pm 0.29$ and $w = -1.26^{+0.34}_{-0.27}$, for redMaGiC and MagLim catalogues, respectively. These results are compatible with the corresponding DES Y3 results in configuration space and pave the way for harmonic space analyses using the DES Y6 data., Comment: To be submitted to MNRAS

The residuals of the distance moduli of Type Ia supernovae (SN Ia) relative to a Hubble diagram fit contain information about the inhomogeneity of the universe, due to weak lensing magnification by foreground matter. By correlating the residuals of the Dark Energy Survey Year 5 SN Ia sample (DES-SN5YR) with extra-galactic foregrounds from the DES Y3 Gold catalog, we detect the presence of lensing at $6.0 \sigma$ significance. This is the first detection with a significance level above $5\sigma$. Constraints on the effective mass-to-light ratios and radial profiles of dark-matter haloes surrounding individual galaxies are also obtained. We show that the scatter of SNe Ia around the Hubble diagram is reduced by modifying the standardisation of the distance moduli to include an easily calculable de-lensing (i.e., environmental) term. We use the de-lensed distance moduli to recompute cosmological parameters derived from SN Ia, finding in Flat $w$CDM a difference of $\Delta \Omega_{\rm M} = +0.036$ and $\Delta w = -0.056$ compared to the unmodified distance moduli, a change of $\sim 0.3\sigma$. We argue that our modelling of SN Ia lensing will lower systematics on future surveys with higher statistical power. We use the observed dispersion of lensing in DES-SN5YR to constrain $\sigma_8$, but caution that the fit is sensitive to uncertainties at small scales. Nevertheless, our detection of SN Ia lensing opens a new pathway to study matter inhomogeneity that complements galaxy-galaxy lensing surveys and has unrelated systematics., Comment: Submitted to MNRAS

We present $griz$ photometric light curves for the full 5 years of the Dark Energy Survey Supernova program (DES-SN), obtained with both forced Point Spread Function (PSF) photometry on Difference Images (DIFFIMG) performed during survey operations, and Scene Modelling Photometry (SMP) on search images processed after the survey. This release contains $31,636$ DIFFIMG and $19,706$ high-quality SMP light curves, the latter of which contains $1635$ photometrically-classified supernovae that pass cosmology quality cuts. This sample spans the largest redshift ($z$) range ever covered by a single SN survey ($0.1

Type Ia Supernovae (SNe Ia) are a critical tool in measuring the accelerating expansion of the universe. Recent efforts to improve these standard candles have focused on incorporating the effects of dust on distance measurements with SNe Ia. In this paper, we use the state-of-the-art Dark Energy Survey 5 year sample to evaluate two different families of dust models: empirical extinction models derived from SNe Ia data, and physical attenuation models from the spectra of galaxies. Among the SNe Ia-derived models, we find that a logistic function of the total-to-selective extinction RV best recreates the correlations between supernova distance measurements and host galaxy properties, though an additional 0.02 magnitudes of grey scatter are needed to fully explain the scatter in SNIa brightness in all cases. These empirically-derived extinction distributions are highly incompatible with the physical attenuation models from galactic spectral measurements. From these results, we conclude that SNe Ia must either preferentially select extreme ends of galactic dust distributions, or that the characterisation of dust along the SNe Ia line-of-sight is incompatible with that of galactic dust distributions.

We present a precise measurement of cosmological time dilation using the light curves of 1504 type Ia supernovae from the Dark Energy Survey spanning a redshift range $0.1\lesssim z\lesssim 1.2$. We find that the width of supernova light curves is proportional to $(1+z)$, as expected for time dilation due to the expansion of the Universe. Assuming type Ia supernovae light curves are emitted with a consistent duration $\Delta t_{\rm em}$, and parameterising the observed duration as $\Delta t_{\rm obs}=\Delta t_{\rm em}(1+z)^b$, we fit for the form of time dilation using two methods. Firstly, we find that a power of $b \approx 1$ minimises the flux scatter in stacked subsamples of light curves across different redshifts. Secondly, we fit each target supernova to a stacked light curve (stacking all supernovae with observed bandpasses matching that of the target light curve) and find $b=1.003\pm0.005$ (stat) $\pm\,0.010$ (sys). Thanks to the large number of supernovae and large redshift-range of the sample, this analysis gives the most precise measurement of cosmological time dilation to date, ruling out any non-time-dilating cosmological models at very high significance., Comment: 14 pages, 13 figures. Updated in response to reviewer feedback. Accepted for publication in the Monthly Notices of the Royal Astronomical Society (MNRAS)

We measure the current expansion rate of the Universe, Hubble's constant $H_0$, by calibrating the absolute magnitudes of supernovae to distances measured by Baryon Acoustic Oscillations. This `inverse distance ladder' technique provides an alternative to calibrating supernovae using nearby absolute distance measurements, replacing the calibration with a high-redshift anchor. We use the recent release of 1829 supernovae from the Dark Energy Survey spanning $0.01\lt z \lt1.13$ anchored to the recent Baryon Acoustic Oscillation measurements from DESI spanning $0.30 \lt z_{\mathrm{eff}} \lt 2.33$. To trace cosmology to $z=0$, we use the third-, fourth- and fifth-order cosmographic models, which, by design, are agnostic about the energy content and expansion history of the universe. With the inclusion of the higher-redshift DESI-BAO data, the third-order model is a poor fit to both data sets, with the fourth-order model being preferred by the Akaike Information Criterion. Using the fourth-order cosmographic model, we find $H_0=67.19^{+0.66}_{-0.64}\mathrm{~km} \mathrm{~s}^{-1} \mathrm{~Mpc}^{-1}$, in agreement with the value found by Planck without the need to assume Flat-$\Lambda$CDM. However the best-fitting expansion history differs from that of Planck, providing continued motivation to investigate these tensions.

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$., Comment: 24 pages, 13 figures, to be submitted to PRD. Comments welcome!

Extracting precise cosmology from weak lensing surveys requires modelling the non-linear matter power spectrum, which is suppressed at small scales due to baryonic feedback processes. However, hydrodynamical galaxy formation simulations make widely varying predictions for the amplitude and extent of this effect. We use measurements of Dark Energy Survey Year 3 weak lensing (WL) and Atacama Cosmology Telescope DR5 kinematic Sunyaev-Zel'dovich (kSZ) to jointly constrain cosmological and astrophysical baryonic feedback parameters using a flexible analytical model, `baryonification'. First, using WL only, we compare the $S_8$ constraints using baryonification to a simulation-calibrated halo model, a simulation-based emulator model and the approach of discarding WL measurements on small angular scales. We find that model flexibility can shift the value of $S_8$ and degrade the uncertainty. The kSZ provides additional constraints on the astrophysical parameters and shifts $S_8$ to $S_8=0.823^{+0.019}_{-0.020}$, a higher value than attained using the WL-only analysis. We measure the suppression of the non-linear matter power spectrum using WL + kSZ and constrain a mean feedback scenario that is more extreme than the predictions from most hydrodynamical simulations. We constrain the baryon fractions and the gas mass fractions and find them to be generally lower than inferred from X-ray observations and simulation predictions. We conclude that the WL + kSZ measurements provide a new and complementary benchmark for building a coherent picture of the impact of gas around galaxies across observations.

We measure the stacked lensing signal in the direction of galaxy clusters in the Dark Energy Survey Year 3 (DES Y3) redMaPPer sample, using cosmic microwave background (CMB) temperature data from SPT-3G, the third-generation CMB camera on the South Pole Telescope (SPT). We estimate the lensing signal using temperature maps constructed from the initial 2 years of data from the SPT-3G 'Main' survey, covering 1500 deg$^2$ of the Southern sky. We then use this signal as a proxy for the mean cluster mass of the DES sample. In this work, we employ three versions of the redMaPPer catalogue: a Flux-Limited sample containing 8865 clusters, a Volume-Limited sample with 5391 clusters, and a Volume&Redshift-Limited sample with 4450 clusters. For the three samples, we find the mean cluster masses to be ${M}_{200{\rm{m}}}=1.66\pm0.13$ [stat.]$\pm0.03$ [sys.], $1.97\pm0.18$ [stat.]$\pm0.05$ [sys.], and $2.11\pm0.20$ [stat.]$\pm0.05$ [sys.]$\times{10}^{14}\ {\rm{M}}_{\odot }$, respectively. This is a factor of $\sim2$ improvement relative to the precision of measurements with previous generations of SPT surveys and the most constraining cluster mass measurements using CMB cluster lensing to date. Overall, we find no significant tensions between our results and masses given by redMaPPer mass-richness scaling relations of previous works, which were calibrated using CMB cluster lensing, optical weak lensing, and velocity dispersion measurements from various combinations of DES, SDSS and Planck data. We then divide our sample into 3 redshift and 3 richness bins, finding no significant tensions with optical weak-lensing calibrated masses in these bins. We forecast a $5.7\%$ constraint on the mean cluster mass of the DES Y3 sample with the complete SPT-3G surveys when using both temperature and polarization data and including an additional $\sim1400$ deg$^2$ of observations from the 'Extended' SPT-3G survey., Comment: 23 pages, 9 figures, accepted for publication in JCAP. Minor changes and corrections have been made relative to v1

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

Current and future Type Ia Supernova (SN Ia) surveys will need to adopt new approaches to classifying SNe and obtaining their redshifts without spectra if they wish to reach their full potential. We present here a novel approach that uses only photometry to identify SNe Ia in the 5-year Dark Energy Survey (DES) dataset using the SuperNNova classifier. Our approach, which does not rely on any information from the SN host-galaxy, recovers SNe Ia that might otherwise be lost due to a lack of an identifiable host. We select 2,298 high-quality SNe Ia from the DES 5-year dataset an almost complete sample of detected SNe Ia. More than 700 of these have no spectroscopic host redshift and are potentially new SNIa compared to the DES-SN5YR cosmology analysis. To analyse these SNe Ia, we derive their redshifts and properties using only their light-curves with a modified version of the SALT2 light-curve fitter. Compared to other DES SN Ia samples with spectroscopic redshifts, our new sample has in average higher redshift, bluer and broader light-curves, and fainter host-galaxies. Future surveys such as LSST will also face an additional challenge, the scarcity of spectroscopic resources for follow-up. When applying our novel method to DES data, we reduce the need for follow-up by a factor of four and three for host-galaxy and live SN respectively compared to earlier approaches. Our novel method thus leads to better optimisation of spectroscopic resources for follow-up., Comment: Accepted MNRAS

In this paper we present and validate the galaxy sample used for the analysis of the baryon acoustic oscillation (BAO) signal in the Dark Energy Survey (DES) Y6 data. The definition is based on a color and redshift-dependent magnitude cut optimized to select galaxies at redshifts higher than 0.6, while ensuring a high-quality photo-$z$ determination. The optimization is performed using a Fisher forecast algorithm, finding the optimal $i$-magnitude cut to be given by $i$<19.64+2.894$z_{\rm ph}$. For the optimal sample, we forecast an increase in precision in the BAO measurement of $\sim$25% with respect to the Y3 analysis. Our BAO sample has a total of 15,937,556 galaxies in the redshift range 0.6<$z_{\rm ph}$<1.2, and its angular mask covers 4,273.42 deg${}^2$ to a depth of $i$=22.5. We validate its redshift distributions with three different methods: directional neighborhood fitting algorithm (DNF), which is our primary photo-$z$ estimation; direct calibration with spectroscopic redshifts from VIPERS; and clustering redshift using SDSS galaxies. The fiducial redshift distribution is a combination of these three techniques performed by modifying the mean and width of the DNF distributions to match those of VIPERS and clustering redshift. In this paper we also describe the methodology used to mitigate the effect of observational systematics, which is analogous to the one used in the Y3 analysis. This paper is one of the two dedicated to the analysis of the BAO signal in DES Y6. In its companion paper, we present the angular diameter distance constraints obtained through the fitting to the BAO scale., Comment: 23 pages, 10 figures. Submitted to PRD

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_{\rm 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_{\rm eff})/r_d$ = 19.51$\pm$0.41 (at 68.3% confidence interval), resulting from comparing the BAO position in our data to that predicted by Planck $\Lambda$CDM via the BAO shift parameter $\alpha=(D_M/r_d)/(D_M/r_d)_{\rm Planck}$. To achieve this, the BAO shift is measured with three different methods, Angular Correlation Function (ACF), Angular Power Spectrum (APS), and Projected Correlation Function (PCF) obtaining $\alpha=$ 0.952$\pm$0.023, 0.962$\pm$0.022, and 0.955$\pm$0.020, respectively, which we combine to $\alpha=$ 0.957$\pm$0.020, including systematic errors. When compared with the $\Lambda$CDM model that best fits Planck data, this measurement is found to be 4.3% and 2.1$\sigma$ below the angular BAO scale predicted. To date, it represents the most precise angular BAO measurement at $z$>0.75 from any survey and the most precise measurement at any redshift from photometric surveys. The analysis was performed blinded to the BAO position and it is shown to be robust against analysis choices, data removal, redshift calibrations and observational systematics., Comment: Submitted to PRD, 39 pages, 12 figures

Number counts of galaxy clusters across redshift are a powerful cosmological probe, if a precise and accurate reconstruction of the underlying mass distribution is performed -- a challenge called mass calibration. With the advent of wide and deep photometric surveys, weak gravitational lensing by clusters has become the method of choice to perform this measurement. We measure and validate the weak gravitational lensing (WL) signature in the shape of galaxies observed in the first 3 years of the DES Y3 caused by galaxy clusters selected in the first all-sky survey performed by SRG/eROSITA. These data are then used to determine the scaling between X-ray photon count rate of the clusters and their halo mass and redshift. We empirically determine the degree of cluster member contamination in our background source sample. The individual cluster shear profiles are then analysed with a Bayesian population model that self-consistently accounts for the lens sample selection and contamination, and includes marginalization over a host of instrumental and astrophysical systematics. To quantify the accuracy of the mass extraction of that model, we perform mass measurements on mock cluster catalogs with realistic synthetic shear profiles. This allows us to establish that hydro-dynamical modelling uncertainties at low lens redshifts ($z<0.6$) are the dominant systematic limitation. At high lens redshift the uncertainties of the sources' photometric redshift calibration dominate. With regard to the X-ray count rate to halo mass relation, we constrain all its parameters. This work sets the stage for a joint analysis with the number counts of eRASS1 clusters to constrain a host of cosmological parameters. We demonstrate that WL mass calibration of galaxy clusters can be performed successfully with source galaxies whose calibration was performed primarily for cosmic shear experiments., Comment: 27 pages, 18 figures, 2 appendices, submitted to A\&A

Cosmological analyses using galaxy clusters in optical/NIR photometric surveys require robust characterization of their galaxy content. Precisely determining which galaxies belong to a cluster is crucial. In this paper, we present the COlor Probabilistic Assignment of Clusters And BAyesiaN Analysis (Copacabana) algorithm. Copacabana computes membership probabilities for {\it all} galaxies within an aperture centred on the cluster using photometric redshifts, colours, and projected radial probability density functions. We use simulations to validate Copacabana and we show that it achieves up to 89\% membership accuracy with a mild dependency on photometric redshift uncertainties and choice of aperture size. We find that the precision of the photometric redshifts has the largest impact on the determination of the membership probabilities followed by the choice of the cluster aperture size. We also quantify how much these uncertainties in the membership probabilities affect the stellar mass--cluster mass scaling relation, a relation that directly impacts cosmology. Using the sum of the stellar masses weighted by membership probabilities ($\mu_{\star}$) as the observable, we find that Copacabana can reach an accuracy of 0.06 dex in the measurement of the scaling relation. These results indicate the potential of Copacabana and $\mu_{\star}$ to be used in cosmological analyses of optically selected clusters in the future., Comment: 14 pages, 10 figures, submitted to MNRAS

We present the full Hubble diagram of photometrically-classified Type Ia supernovae (SNe Ia) from the Dark Energy Survey supernova program (DES-SN). DES-SN discovered more than 20,000 SN candidates and obtained spectroscopic redshifts of 7,000 host galaxies. Based on the light-curve quality, we select 1635 photometrically-identified SNe Ia with spectroscopic redshift 0.10$< z <$1.13, which is the largest sample of supernovae from any single survey and increases the number of known $z>0.5$ supernovae by a factor of five. In a companion paper, we present cosmological results of the DES-SN sample combined with 194 spectroscopically-classified SNe Ia at low redshift as an anchor for cosmological fits. Here we present extensive modeling of this combined sample and validate the entire analysis pipeline used to derive distances. We show that the statistical and systematic uncertainties on cosmological parameters are $\sigma_{\Omega_M,{\rm stat+sys}}^{\Lambda{\rm CDM}}=$0.017 in a flat $\Lambda$CDM model, and $(\sigma_{\Omega_M},\sigma_w)_{\rm stat+sys}^{w{\rm CDM}}=$(0.082, 0.152) in a flat $w$CDM model. Combining the DES SN data with the highly complementary CMB measurements by Planck Collaboration (2020) reduces uncertainties on cosmological parameters by a factor of 4. In all cases, statistical uncertainties dominate over systematics. We show that uncertainties due to photometric classification make up less than 10% of the total systematic uncertainty budget. This result sets the stage for the next generation of SN cosmology surveys such as the Vera C. Rubin Observatory's Legacy Survey of Space and Time., Comment: 39 pages, 19 figures; Submitted to ApJ; companion paper Dark Energy Collaboration et al. on consecutive arxiv number 2401.02929

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 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

Context. The determination of accurate photometric redshifts (photo-zs) in large imaging galaxy surveys is key for cosmological studies. One of the most common approaches are machine learning techniques. These methods require a spectroscopic or reference sample to train the algorithms. Attention has to be paid to the quality and properties of these samples since they are key factors in the estimation of reliable photo-zs. Aims. The goal of this work is to calculate the photo-zs for the Y3 DES Deep Fields catalogue using the DNF machine learning algorithm. Moreover, we want to develop techniques to assess the incompleteness of the training sample and metrics to study how incompleteness affects the quality of photometric redshifts. Finally, we are interested in comparing the performance obtained with respect to the EAzY template fitting approach on Y3 DES Deep Fields catalogue. Methods. We have emulated -- at brighter magnitude -- the training incompleteness with a spectroscopic sample whose redshifts are known to have a measurable view of the problem. We have used a principal component analysis to graphically assess incompleteness and to relate it with the performance parameters provided by DNF. Finally, we have applied the results about the incompleteness to the photo-z computation on Y3 DES Deep Fields with DNF and estimated its performance. Results. The photo-zs for the galaxies on DES Deep Fields have been computed with the DNF algorithm and added to the Y3 DES Deep Fields catalogue. They are available at https://des.ncsa.illinois.edu/releases/y3a2/Y3deepfields. Some techniques have been developed to evaluate the performance in the absence of "true" redshift and to assess completeness. We have studied... (Partial abstract), Comment: 14 pages, 17 figures

We present a catalog of 689 galaxy cluster candidates detected at significance $\xi>4$ via their thermal Sunyaev-Zel'dovich (SZ) effect signature in 95 and 150 GHz data from the 500-square-degree SPTpol survey. We use optical and infrared data from the Dark Energy Camera and the Wide-field Infrared Survey Explorer (WISE) and \spitzer \ satellites, to confirm 544 of these candidates as clusters with $\sim94\%$ purity. The sample has an approximately redshift-independent mass threshold at redshift $z>0.25$ and spans $1.5 \times 10^{14} < M_{500c} < 9.1 \times 10^{14}$ $M_\odot/h_{70}$ \ and $0.031$. We use external radio data from the Sydney University Molonglo Sky Survey (SUMSS) to estimate contamination to the SZ signal from synchrotron sources. The contamination reduces the recovered $\xi$ by a median value of 0.032, or $\sim0.8\%$ of the $\xi=4$ threshold value, and $\sim7\%$ of candidates have a predicted contamination greater than $\Delta \xi = 1$. With the exception of a small number of systems $(<1\%)$, an analysis of clusters detected in single-frequency 95 and 150 GHz data shows no significant contamination of the SZ signal by emission from dusty or synchrotron sources. This cluster sample will be a key component in upcoming astrophysical and cosmological analyses of clusters. The SPTpol millimeter-wave maps and associated data products used to produce this sample are available at https://pole.uchicago.edu/public/data/sptpol_500d_clusters/index.html, and the NASA LAMBDA website. An interactive sky server with the SPTpol maps and Dark Energy Survey data release 2 images is also available at NCSA https://skyviewer.ncsa.illinois.edu., Comment: Matches version accepted by OJA. 19 pages + references, 14 figures, cluster candidate table provided in Appendix. Data products available at https://pole.uchicago.edu/public/data/sptpol_500d_clusters/index.html and an interactive sky server at https://skyviewer.ncsa.illinois.edu

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., Comment: 25 pages, 18 figures. Comments welcome!

We use Dark Energy Survey Year 3 (DES Y3) clusters with archival X-ray data from XMM-Newton and Chandra to assess the centering performance of the redMaPPer cluster finder and to measure key richness observable scaling relations. In terms of centering, we find that 10-20% of redMaPPer clusters are miscentered with no significant difference in bins of low versus high richness ($20<\lambda<40$ and $\lambda>40$) or redshift ($0.2

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 multi-observable 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., Comment: Accepted for publication in Phys. Rev. D. arXiv v2 corresponds to published article

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., Comment: [v2]: Version accepted to MNRAS

We present an extension to a Sunyaev-Zel'dovich Effect (SZE) selected cluster catalog based on observations from the South Pole Telescope (SPT); this catalog extends to lower signal-to-noise than the previous SPT-SZ catalog and therefore includes lower mass clusters. Optically derived redshifts, centers, richnesses and morphological parameters together with catalog contamination and completeness statistics are extracted using the multi-component matched filter algorithm (MCMF) applied to the S/N>4 SPT-SZ candidate list and the Dark Energy Survey (DES) photometric galaxy catalog. The main catalog contains 811 sources above S/N=4, has 91% purity and is 95% complete with respect to the original SZE selection. It contains 50% more total clusters and twice as many clusters above z=0.8 in comparison to the original SPT-SZ sample. The MCMF algorithm allows us to define subsamples of the desired purity with traceable impact on catalog completeness. As an example, we provide two subsamples with S/N>4.25 and S/N>4.5 for which the sample contamination and cleaning-induced incompleteness are both as low as the expected Poisson noise for samples of their size. The subsample with S/N>4.5 has 98% purity and 96% completeness, and will be included in a combined SPT cluster and DES weak-lensing cosmological analysis. We measure the number of false detections in the SPT-SZ candidate list as function of S/N, finding that it follows that expected from assuming Gaussian noise, but with a lower amplitude compared to previous estimates from simulations., Comment: 16 pages, 17 figures, submitted to MNRAS

We introduce a software suite developed for galaxy cluster cosmological analysis with the Dark Energy Survey Data. Cosmological analyses based on galaxy cluster number counts and weak-lensing measurements need efficient software infrastructure to explore an increasingly large parameter space, and account for various cosmological and astrophysical effects. Our software package is designed to model the cluster observables in a wide-field optical survey, including galaxy cluster counts, their averaged weak-lensing masses, or the cluster's averaged weak-lensing radial signals. To ensure maximum efficiency, this software package is developed in C++ in the CosmoSIS software framework, making use of the CUBA integration library. We also implement a testing and validation scheme to ensure the quality of the package. We demonstrate the effectiveness of this development by applying the software to the Dark Energy Survey Year 1 galaxy cluster cosmological data sets, and acquired cosmological constraints that are consistent with the fiducial Dark Energy Survey analysis.

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., Comment: 26 pages, 30 figures (including appendices). Data associated with this article is available at https://github.com/itrharrison/actdr4kappa-x-desy3gamma-data

Using the full six years of imaging data from the Dark Energy Survey, we study the surface brightness profiles of galaxy cluster central galaxies and intra-cluster light. We apply a ``stacking'' method to over four thousand galaxy clusters identified by the redMaPPer cluster finding algorithm in the redshift range of 0.2 to 0.5. This yields high signal-to-noise radial profile measurements of the central galaxy and intra-cluster light out to 1 Mpc from the cluster center. Using redMaPPer richness as a cluster mass indicator, we find that the intra-cluster light brightness has a strong mass dependence throughout the 0.2 to 0.5 redshift range, and the dependence grows stronger at a larger radius. In terms of redshift evolution, we find some evidence that the central galaxy, as well as the diffuse light within the transition region between the cluster central galaxy and intra-cluster light within 80 kpc from the center, may be growing over time. At larger radii, more than 80 kpc away from the cluster center, we do not find evidence of additional redshift evolution beyond the cluster mass dependence, which is consistent with the findings from the IllustrisTNG hydrodynamic simulation. We speculate that the major driver of intra-cluster light growth, especially at large radii, is associated with cluster mass growth. Finally, we find that the color of the cluster central galaxy and intra-cluster light displays a radial gradient that becomes bluer at a larger radius, which is consistent with a stellar stripping and disruption origin of intra-cluster light as suggested by simulation studies., Comment: Submitted to MNRAS

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., Comment: 21 pages, 12 figures

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., Comment: 5 pages, 2 figures, submitted to MNRAS Letters

Redshift measurements, primarily obtained from host galaxies, are essential for inferring cosmological parameters from type Ia supernovae (SNe Ia). Matching SNe to host galaxies using images is non-trivial, resulting in a subset of SNe with mismatched hosts and thus incorrect redshifts. We evaluate the host galaxy mismatch rate and resulting biases on cosmological parameters from simulations modeled after the Dark Energy Survey 5-Year (DES-SN5YR) photometric sample. For both DES-SN5YR data and simulations, we employ the directional light radius method for host galaxy matching. In our SN Ia simulations, we find that 1.7% of SNe are matched to the wrong host galaxy, with redshift difference between the true and matched host of up to 0.6. Using our analysis pipeline, we determine the shift in the dark energy equation of state parameter (Dw) due to including SNe with incorrect host galaxy matches. For SN Ia-only simulations, we find Dw = 0.0013 +/- 0.0026 with constraints from the cosmic microwave background (CMB). Including core-collapse SNe and peculiar SNe Ia in the simulation, we find that Dw ranges from 0.0009 to 0.0032 depending on the photometric classifier used. This bias is an order of magnitude smaller than the expected total uncertainty on w from the DES-SN5YR sample of around 0.03. We conclude that the bias on w from host galaxy mismatch is much smaller than the uncertainties expected from the DES-SN5YR sample, but we encourage further studies to reduce this bias through better host-matching algorithms or selection cuts., Comment: Accepted by ApJ

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

We present a joint cosmic shear analysis of the Dark Energy Survey (DES Y3) and the Kilo-Degree Survey (KiDS-1000) in a collaborative effort between the two survey teams. We find consistent cosmological parameter constraints between DES Y3 and KiDS-1000 which, when combined in a joint-survey analysis, constrain the parameter $S_8 = \sigma_8 \sqrt{\Omega_{\rm m}/0.3}$ with a mean value of $0.790^{+0.018}_{-0.014}$. The mean marginal is lower than the maximum a posteriori estimate, $S_8=0.801$, owing to skewness in the marginal distribution and projection effects in the multi-dimensional parameter space. Our results are consistent with $S_8$ constraints from observations of the cosmic microwave background by Planck, with agreement at the $1.7\sigma$ level. We use a Hybrid analysis pipeline, defined from a mock survey study quantifying the impact of the different analysis choices originally adopted by each survey team. We review intrinsic alignment models, baryon feedback mitigation strategies, priors, samplers and models of the non-linear matter power spectrum., Comment: 40 pages, 21 figures, 15 tables, accepted Open Journal of Astrophysics. Download the chains from https://des.ncsa.illinois.edu/releases/y3a2/Y3key-joint-des-kids or create your own chains with CosmoSIS using https://github.com/joezuntz/cosmosis-standard-library/blob/main/examples/des-y3_and_kids-1000.ini Watch the core team discuss this analysis at https://cosmologytalks.com/2023/05/26/des-kids

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., Comment: 19 pages, 15 figures - matches published version

We search for the presence of cool cores in optically-selected galaxy clusters from the Dark Energy Survey (DES) and investigate their prevalence as a function of redshift and cluster richness. Clusters were selected from the redMaPPer analysis of three years of DES observations that have archival Chandra X-ray observations, giving a sample of 99 clusters with a redshift range of $0.11 < z < 0.87$ and a richness range of $25 < \lambda < 207$. Using the X-ray data, the core temperature was compared to the outer temperature to identify clusters where the core temperature is a factor of 0.7 or less than the outer temperature. We found a cool core fraction of approximately 20% with no significant trend in the cool core fraction with either redshift or richness., Comment: shortened version accepted to RNAAS

Wavelength-dependent atmospheric effects impact photometric supernova flux measurements for ground-based observations. We present corrections on supernova flux measurements from the Dark Energy Survey Supernova Program's 5YR sample (DES-SN5YR) for differential chromatic refraction (DCR) and wavelength-dependent seeing, and we show their impact on the cosmological parameters $w$ and $\Omega_m$. We use $g-i$ colors of Type Ia supernovae (SNe Ia) to quantify astrometric offsets caused by DCR and simulate point spread functions (PSFs) using the GalSIM package to predict the shapes of the PSFs with DCR and wavelength-dependent seeing. We calculate the magnitude corrections and apply them to the magnitudes computed by the DES-SN5YR photometric pipeline. We find that for the DES-SN5YR analysis, not accounting for the astrometric offsets and changes in the PSF shape cause an average bias of $+0.2$ mmag and $-0.3$ mmag respectively, with standard deviations of $0.7$ mmag and $2.7$ mmag across all DES observing bands (\textit{griz}) throughout all redshifts. When the DCR and seeing effects are not accounted for, we find that $w$ and $\Omega_m$ are lower by less than $0.004\pm0.02$ and $0.001\pm0.01$ respectively, with $0.02$ and $0.01$ being the $1\sigma$ statistical uncertainties. Although we find that these biases do not limit the constraints of the DES-SN5YR sample, future surveys with much higher statistics, lower systematics, and especially those that observe in the $u$ band will require these corrections as wavelength-dependent atmospheric effects are larger at shorter wavelengths. We also discuss limitations of our method and how they can be better accounted for in future surveys., Comment: 15 pages, 13 figures, accepted by the Astronomical Journal

Clusters of galaxies are sensitive to the most nonlinear peaks in the cosmic density field. The weak gravitational lensing of background galaxies by clusters can allow us to infer their masses. However, galaxies associated with the local environment of the cluster can also be intrinsically aligned due to the local tidal gradient, contaminating any cosmology derived from the lensing signal. We measure this intrinsic alignment in Dark Energy Survey (DES) Year 1 redMaPPer clusters. We find evidence of a non-zero mean radial alignment of galaxies within clusters between redshift 0.1-0.7. We find a significant systematic in the measured ellipticities of cluster satellite galaxies that we attribute to the central galaxy flux and other intracluster light. We attempt to correct this signal, and fit a simple model for intrinsic alignment amplitude ($A_{\textrm{IA}}$) to the measurement, finding $A_{\textrm{IA}}=0.15\pm 0.04$, when excluding data near the edge of the cluster. We find a significantly stronger alignment of the central galaxy with the cluster dark matter halo at low redshift and with higher richness and central galaxy absolute magnitude (proxies for cluster mass). This is an important demonstration of the ability of large photometric data sets like DES to provide direct constraints on the intrinsic alignment of galaxies within clusters. These measurements can inform improvements to small-scale modeling and simulation of the intrinsic alignment of galaxies to help improve the separation of the intrinsic alignment signal in weak lensing studies., Comment: 14 pages, 13 figures. Accepted to MNRAS

We present direct constraints on galaxy intrinsic alignments using the Dark Energy Survey Year 3 (DES Y3), the Extended Baryon Oscillation Spectroscopic Survey (eBOSS) and its precursor, the Baryon Oscillation Spectroscopic Survey (BOSS). Our measurements incorporate photometric red sequence (redMaGiC) galaxies from DES with median redshift $z\sim0.2-1.0$, luminous red galaxies (LRGs) from eBOSS at $z\sim0.8$, and also a SDSS-III BOSS CMASS sample at $z\sim0.5$. We measure two point intrinsic alignment correlations, which we fit using a model that includes lensing, magnification and photometric redshift error. Fitting on scales $6

Recent cosmological analyses with large-scale structure and weak lensing measurements, usually referred to as 3$\times$2pt, had to discard a lot of signal-to-noise from small scales due to our inability to accurately model non-linearities and baryonic effects. Galaxy-galaxy lensing, or the position-shear correlation between lens and source galaxies, is one of the three two-point correlation functions that are included in such analyses, usually estimated with the mean tangential shear. However, tangential shear measurements at a given angular scale $\theta$ or physical scale $R$ carry information from all scales below that, forcing the scale cuts applied in real data to be significantly larger than the scale at which theoretical uncertainties become problematic. Recently there have been a few independent efforts that aim to mitigate the non-locality of the galaxy-galaxy lensing signal. Here we perform a comparison of the different methods, including the Y-transformation, the Point-Mass marginalization methodology and the Annular Differential Surface Density statistic. We do the comparison at the cosmological constraints level in a combined galaxy clustering and galaxy-galaxy lensing analysis. We find that all the estimators yield equivalent cosmological results assuming a simulated Rubin Observatory Legacy Survey of Space and Time (LSST) Year 1 like setup and also when applied to DES Y3 data. With the LSST Y1 setup, we find that the mitigation schemes yield $\sim$1.3 times more constraining $S_8$ results than applying larger scale cuts without using any mitigation scheme., Comment: 11+4 pages, 4+4 figures. Matches the accepted version in MNRAS

The fiducial cosmological analyses of imaging galaxy surveys like the Dark Energy Survey (DES) typically probe the Universe at redshifts $z < 1$. This is mainly because of the limited depth of these surveys, and also because such analyses rely heavily on galaxy lensing, which is more efficient at low redshifts. In this work we present the selection and characterization of high-redshift galaxy samples using DES Year 3 data, and the analysis of their galaxy clustering measurements. In particular, we use galaxies that are fainter than those used in the previous DES Year 3 analyses and a Bayesian redshift scheme to define three tomographic bins with mean redshifts around $z \sim 0.9$, $1.2$ and $1.5$, which significantly extend the redshift coverage of the fiducial DES Year 3 analysis. These samples contain a total of about 9 million galaxies, and their galaxy density is more than 2 times higher than those in the DES Year 3 fiducial case. We characterize the redshift uncertainties of the samples, including the usage of various spectroscopic and high-quality redshift samples, and we develop a machine-learning method to correct for correlations between galaxy density and survey observing conditions. The analysis of galaxy clustering measurements, with a total signal-to-noise $S/N \sim 70$ after scale cuts, yields robust cosmological constraints on a combination of the fraction of matter in the Universe $\Omega_m$ and the Hubble parameter $h$, $\Omega_m h = 0.195^{+0.023}_{-0.018}$, and 2-3% measurements of the amplitude of the galaxy clustering signals, probing galaxy bias and the amplitude of matter fluctuations, $b \sigma_8$. A companion paper $\textit{(in preparation)}$ will present the cross-correlations of these high-$z$ samples with CMB lensing from Planck and SPT, and the cosmological analysis of those measurements in combination with the galaxy clustering presented in this work., Comment: 28 pages, 25 figures. To be submitted to MNRAS. Comments welcome