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.

The Einstein equations relate perturbations in the metric to perturbations in the matter and radiation. Both kinds of perturbations can be decomposed in to scalar, vector, and tensor types which evolve independently of each other. We derive equations for the evolution of the two scalar metric potentials. While vector modes decay rapidly if they are not sourced, tensor modes are important as they describe gravitational waves. There are two functions that characterize the two polarization states of tensor perturbations, which propagate as waves and whose amplitudes fall off due to expansion once they enter the horizon.

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.

The evidence for dark matter (DM) remains compelling, although attempts to understand its particle nature remain inconclusive. One promising method to study DM is detecting DM subhalos through their gravitational interactions with stellar streams. In this study, we apply Neural Posterior Estimation (NPE) to constrain subhalo interaction parameters, including mass, scale radius, velocity, and encounter geometry, from stellar stream kinematics. We generate particle spray simulations based on the Lagrange Cloud stripping technique, focusing on the ATLAS-Aliqa Uma stream as a test case. We train multiple NPE models across multiple observational scenarios, quantifying how kinematic completeness affects inference and forecasting constraints from upcoming surveys including LSST, 4MOST, and 10-year Gaia data. Our results demonstrate that NPE can produce accurate and well-calibrated posteriors. In the idealized case with full 6D coordinates, we achieve subhalo mass uncertainties of 15-20

Weak gravitational lensing (WL) is a key cosmological probe that requires precise measurement of galaxy images to infer shape distortions, or shear, and constrain cosmology. Accurate estimation of the point spread function (PSF) is crucial for shear measurement, but the wavelength dependence of the PSF introduces chromatic biases that can systematically impact shear inference. We focus on biases arising from spectral energy distribution (SED) differences between stars, used for PSF modelling, and galaxies, used for shear measurement. We investigate these effects in Roman's four design reference mission WL bands (Y106, J129, H158, F184) and wide filter (W146). Using Roman-like image simulations, we quantify the induced shear biases and compare them to requirements on those biases. Multiplicative biases over all galaxies hover around similar to 0.2 per cent in the WL bands and 2 per cent in the wide filter, exceeding the mission requirement of |m| < 0.032 and relaxed requirement of |m| < 0.1. In individual redshift bins, biases can reach 0.4-0.9 per cent for the WL bands and 3-6 per cent for the wide filter. Additive biases remain acceptable in the WL bands but exceed systematic limits in the wide filter. We develop and test PSF-level corrections, showing that a first-order correction reduces biases within survey requirements for the WL bands; however, higher-order terms are necessary for the wide filter. Our results highlight the necessity of chromatic corrections for precision WL with Roman and provide a framework for mitigating these biases. Finally, we compare analytical colour-based corrections to self-organizing maps and find that both methods effectively reduce biases.

Precision cosmology benefits from extracting maximal information from cosmic structures, motivating the use of higher-order statistics (HOS) at small spatial scales. However, predicting how baryonic processes modify matter statistics at these scales has been challenging. The baryonic correction model (BCM) addresses this by modifying dark-matteronly simulations to mimic baryonic effects, providing a flexible, simulation-based framework for predicting both two-point and HOS. We show that a 3-parameter version of the BCM can jointly fit weak lensing maps' two-point statistics, wavelet phase harmonics coefficients, scattering coefficients, and the third and fourth moments to within 2% accuracy across all scales l < 2000 and tomographic bins for a DES-Y3-like redshift distribution (z less than or similar to 2), using the FLAMINGO simulations. These results demonstrate the viability of BCM-assisted, simulation-based weak lensing inference of two-point and HOS, paving the way for robust cosmological constraints that fully exploit non-Gaussian information on small spatial scales.

We test the isotropy of Hubble expansion by combining several probes for the first time, constructing full-sky maps of expansion rate variation using Type Ia supernovae, fundamental plane galaxies, and cosmic microwave background (CMB) temperature fluctuations. We find no hint of anisotropy or correlation between early- and late-Universe expansion across all systematic models. The 99% confidence upper limits on expansion rate anisotropy are 0.39% for low-redshift supernovae, 0.95% for high-redshift CMB, and 0.37% when combined at a 60-degree smoothing scale. A significant anomaly in the fundamental plane residual map may reflect systematics in the current dark energy spectroscopic instrument dataset, as evidenced by the absence of cross-correlation with other tracers and its correlation with spatial density variations.

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.

The analyses of the next generation cosmological surveys demand an accurate, efficient, and differentiable method for simulating the universe and its observables across cosmological volumes. We present Hamiltonian ray tracing (HRT) - the first post-Born (accounting for lens-lens coupling and without relying on the Born approximation), threedimensional (without assuming the thin-lens approximation), and on-the-fly (applicable to any structure formation simulations) ray tracing algorithm based on the Hamiltonian formalism. HRT performs symplectic integration of the photon geodesics in a weak gravitational field, and can integrate tightly with any gravity solver, enabling co-evolution of matter particles and light rays with minimal additional computations. We implement HRT in the particle-mesh library pmwd, leveraging hardware accelerators such as GPUs and automatic differentiation capabilities based on JAX. When tested on a point-mass lens, HRT achieves sub-percent accuracy in deflection angles above the resolution limit across both weak and moderately strong lensing regimes. We also test HRT in cosmological simulations on the convergence maps and their power spectra.

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.

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-circle dot over cosmological distances. Using 1532 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 per cent confidence level after marginalization 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 similar to 0.04.

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

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