Search

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

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

We report constraints on a variety of non-standard cosmological models using the full 5-year photometrically-classified type Ia supernova sample from the Dark Energy Survey (DES-SN5YR). Both Akaike Information Criterion (AIC) and Suspiciousness calculations find no strong evidence for or against any of the non-standard models we explore. When combined with external probes, the AIC and Suspiciousness agree that 11 of the 15 models are moderately preferred over Flat-$\Lambda$CDM suggesting additional flexibility in our cosmological models may be required beyond the cosmological constant. We also provide a detailed discussion of all cosmological assumptions that appear in the DES supernova cosmology analyses, evaluate their impact, and provide guidance on using the DES Hubble diagram to test non-standard models. An approximate cosmological model, used to perform bias corrections to the data holds the biggest potential for harbouring cosmological assumptions. We show that even if the approximate cosmological model is constructed with a matter density shifted by $\Delta\Omega_m\sim0.2$ from the true matter density of a simulated data set the bias that arises is sub-dominant to statistical uncertainties. Nevertheless, we present and validate a methodology to reduce this bias., Comment: Published to MNRAS on 20 August 2024; v2 updates to the accepted version

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

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

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

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

We cross-match and compare characteristics of galaxy clusters identified in observations from two sky surveys using two completely different techniques. One sample is optically selected from the analysis of three years of Dark Energy Survey observations using the redMaPPer cluster detection algorithm. The second is X-ray selected from XMM observations analysed by the XMM Cluster Survey. The samples comprise a total area of 57.4 deg$^2$, bounded by the area of 4 contiguous XMM survey regions that overlap the DES footprint. We find that the X-ray selected sample is fully matched with entries in the redMaPPer catalogue, above $\lambda>$20 and within 0.1$< z <$0.9. Conversely, only 38\% of the redMaPPer catalogue is matched to an X-ray extended source. Next, using 120 optically clusters and 184 X-ray selected clusters, we investigate the form of the X-ray luminosity-temperature ($L_{X}-T_{X}$), luminosity-richness ($L_{X}-\lambda$) and temperature-richness ($T_{X}-\lambda$) scaling relations. We find that the fitted forms of the $L_{X}-T_{X}$ relations are consistent between the two selection methods and also with other studies in the literature. However, we find tentative evidence for a steepening of the slope of the relation for low richness systems in the X-ray selected sample. When considering the scaling of richness with X-ray properties, we again find consistency in the relations (i.e., $L_{X}-\lambda$ and $T_{X}-\lambda$) between the optical and X-ray selected samples. This is contrary to previous similar works that find a significant increase in the scatter of the luminosity scaling relation for X-ray selected samples compared to optically selected samples., Comment: Accepted for publication to MNRAS

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

We present a sample of 19,583 ultracool dwarf candidates brighter than z $\leq 23$ selected from the Dark Energy Survey DR2 coadd data matched to VHS DR6, VIKING DR5 and AllWISE covering $\sim$ 4,800 $deg^2$. The ultracool candidates were first pre-selected based on their (i-z), (z-Y), and (Y-J) colours. They were further classified using a method that compares their optical, near-infrared and mid-infrared colours against templates of M, L and T dwarfs. 14,099 objects are presented as new L and T candidates and the remaining objects are from the literature, including 5,342 candidates from our previous work. Using this new and deeper sample of ultracool dwarf candidates we also present: 20 new candidate members to nearby young moving groups (YMG) and associations, variable candidate sources and four new wide binary systems composed of two ultracool dwarfs. Finally, we also show the spectra of twelve new ultracool dwarfs discovered by our group and presented here for the first time. These spectroscopically confirmed objects are a sanity check of our selection of ultracool dwarfs and photometric classification method., Comment: 18 pages, 10 figures, 7 tables. Accepted for publication in MNRAS

We address the problem of optimally identifying all kilonovae detected via gravitational wave emission in the upcoming LIGO/Virgo/KAGRA Collaboration observing run, O4, which is expected to be sensitive to a factor of $\sim 7$ more Binary Neutron Stars alerts than previously. Electromagnetic follow-up of all but the brightest of these new events will require $>1$ meter telescopes, for which limited time is available. We present an optimized observing strategy for the Dark Energy Camera during O4. We base our study on simulations of gravitational wave events expected for O4 and wide-prior kilonova simulations. We derive the detectabilities of events for realistic observing conditions. We optimize our strategy for confirming a kilonova while minimizing telescope time. For a wide range of kilonova parameters, corresponding to a fainter kilonova compared to GW170817/AT2017gfo we find that, with this optimal strategy, the discovery probability for electromagnetic counterparts with the Dark Energy Camera is $\sim 80\%$ at the nominal binary neutron star gravitational wave detection limit for the next LVK observing run (190 Mpc), which corresponds to a $\sim 30\%$ improvement compared to the strategy adopted during the previous observing run. For more distant events ($\sim 330$ Mpc), we reach a $\sim 60\%$ probability of detection, a factor of $\sim 2$ increase. For a brighter kilonova model dominated by the blue component that reproduces the observations of GW170817/AT2017gfo, we find that we can reach $\sim 90\%$ probability of detection out to 330 Mpc, representing an increase of $\sim 20 \%$, while also reducing the total telescope time required to follow-up events by $\sim 20\%$., Comment: 26 pages, 11 figures, accepted by ApJ

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 ultra-faint dwarf galaxy Reticulum II (Ret II) exhibits a unique chemical evolution history, with 72 +10/-12% of its stars strongly enhanced in r-process elements. We present deep Hubble Space Telescope photometry of Ret II and analyze its star formation history. As in other ultra-faint dwarfs, the color-magnitude diagram is best fit by a model consisting of two bursts of star formation. If we assume that the bursts were instantaneous, then the older burst occurred around the epoch of reionization and formed ~80% of the stars in the galaxy, while the remainder of the stars formed ~3 Gyr later. When the bursts are allowed to have nonzero durations we obtain slightly better fits. The best-fitting model in this case consists of two bursts beginning before reionization, with approximately half the stars formed in a short (100 Myr) burst and the other half in a more extended period lasting 2.6 Gyr. Considering the full set of viable star formation history models, we find that 28% of the stars formed within 500 +/- 200 Myr of the onset of star formation. The combination of the star formation history and the prevalence of r-process-enhanced stars demonstrates that the r-process elements in Ret II must have been synthesized early in its initial star-forming phase. We therefore constrain the delay time between the formation of the first stars in Ret II and the r-process nucleosynthesis to be less than 500 Myr. This measurement rules out an r-process source with a delay time of several Gyr or more such as GW170817., Comment: 14 pages, 5 figures, 1 table. Accepted for publication in ApJ

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

Type Ia supernovae (SNe Ia) are useful distance indicators in cosmology, provided their luminosity is standardized by applying empirical corrections based on light-curve properties. One factor behind these corrections is dust extinction, accounted for in the color-luminosity relation of the standardization. This relation is usually assumed to be universal, which could potentially introduce systematics into the standardization. The ``mass-step'' observed for SNe Ia Hubble residuals has been suggested as one such systematic. We seek to obtain a completer view of dust attenuation properties for a sample of 162 SN Ia host galaxies and to probe their link to the ``mass-step''. We infer attenuation laws towards hosts from both global and local (4 kpc) Dark Energy Survey photometry and Composite Stellar Population model fits. We recover a optical depth/attenuation slope relation, best explained by differing star/dust geometry for different galaxy orientations, which is significantly different from the optical depth/extinction slope relation observed directly for SNe. We obtain a large variation of attenuation slopes and confirm these change with host properties, like stellar mass and age, meaning a universal SN Ia correction should ideally not be assumed. Analyzing the cosmological standardization, we find evidence for a ``mass-step'' and a two dimensional ``dust-step'', both more pronounced for red SNe. Although comparable, the two steps are found no to be completely analogous. We conclude that host galaxy dust data cannot fully account for the ``mass-step'', using either an alternative SN standardization with extinction proxied by host attenuation or a ``dust-step'' approach., Comment: 20 pages, 10 figues, 9 tables. Supplementary material included (10 pages). Accepted for publication on A&A

The Jupiter Trojans are a large group of asteroids that are co-orbiting with Jupiter near its L4 and L5 Lagrange points. The study of Jupiter Trojans is crucial for testing different models of planet formation that are directly related to our understanding of solar system evolution. In this work, we select known Jupiter Trojans listed by the Minor Planet Center (MPC) from the full six years dataset (Y6) of the Dark Energy Survey (DES) to analyze their photometric properties. The DES data allow us to study Jupiter Trojans with a fainter magnitude limit than previous studies in a homogeneous survey with $griz$ band measurements. We extract a final catalog of 573 unique Jupiter Trojans. Our sample include 547 asteroids belonging to L5. This is one of the largest analyzed samples for this group. By comparing with the data reported by other surveys we found that the color distribution of L5 Trojans is similar to that of L4 Trojans. We find that L5 Trojans' $g - i$ and $g - r$ colors become less red with fainter absolute magnitudes, a trend also seen in L4 Trojans. Both the L4 and L5 clouds consistently show such a color-size correlation over an absolute magnitude range $11 < H < 18$. We also use DES colors to perform taxonomic classifications. C and P-type asteroids outnumber D-type asteroids in the L5 Trojans DES sample, which have diameters in the 5 - 20 km range. This is consistent with the color-size correlation., Comment: This manuscript is accepted for publication in PSJ. There is a full version of table 1

We perform a search for galaxy-galaxy strong lens systems using a convolutional neural network (CNN) applied to imaging data from the first public data release of the DECam Local Volume Exploration Survey (DELVE), which contains $\sim 520$ million astronomical sources covering $\sim 4,000$ $\mathrm{deg}^2$ of the southern sky to a $5\sigma$ point-source depth of $g=24.3$, $r=23.9$, $i=23.3$, and $z=22.8$ mag. Following the methodology of similar searches using DECam data, we apply color and magnitude cuts to select a catalog of $\sim 11$ million extended astronomical sources. After scoring with our CNN, the highest scoring 50,000 images were visually inspected and assigned a score on a scale from 0 (definitely not a lens) to 3 (very probable lens). We present a list of 581 strong lens candidates, 562 of which are previously unreported. We categorize our candidates using their human-assigned scores, resulting in 55 Grade A candidates, 149 Grade B candidates, and 377 Grade C candidates. We additionally highlight eight potential quadruply lensed quasars from this sample. Due to the location of our search footprint in the northern Galactic cap ($b > 10$ deg) and southern celestial hemisphere (${\rm Dec.}<0$ deg), our candidate list has little overlap with other existing ground-based searches. Where our search footprint does overlap with other searches, we find a significant number of high-quality candidates which were previously unidentified, indicating a degree of orthogonality in our methodology. We report properties of our candidates including apparent magnitude and Einstein radius estimated from the image separation., Comment: 24 pages; published version (ApJ)

We cross-correlate positions of galaxies measured in data from the first three years of the Dark Energy Survey with Compton-$y$-maps generated using data from the South Pole Telescope (SPT) and the {\it Planck} mission. We model this cross-correlation measurement together with the galaxy auto-correlation to constrain the distribution of gas in the Universe. We measure the hydrostatic mass bias or, equivalently, the mean halo bias-weighted electron pressure $\langle b_{h}P_{e}\rangle$, using large-scale information. We find $\langle b_{h}P_{e}\rangle$ to be $[0.16^{+0.03}_{-0.04},0.28^{+0.04}_{-0.05},0.45^{+0.06}_{-0.10},0.54^{+0.08}_{-0.07},0.61^{+0.08}_{-0.06},0.63^{+0.07}_{-0.08}]$ meV cm$^{-3}$ at redshifts $z \sim [0.30, 0.46, 0.62,0.77, 0.89, 0.97]$. These values are consistent with previous work where measurements exist in the redshift range. We also constrain the mean gas profile using small-scale information, enabled by the high-resolution of the SPT data. We compare our measurements to different parametrized profiles based on the cosmo-OWLS hydrodynamical simulations. We find that our data are consistent with the simulation that assumes an AGN heating temperature of $10^{8.5}$K but are incompatible with the model that assumes an AGN heating temperature of $10^{8.0}$K. These comparisons indicate that the data prefer a higher value of electron pressure than the simulations within $r_{500c}$ of the galaxies' halos., Comment: 20 pages, 9 figures. Submitted to MNRAS

We present the first systematic follow-up of Planck Sunyaev-Zeldovich effect (SZE) selected candidates down to signal-to-noise (S/N) of 3 over the 5000 deg$^2$ covered by the Dark Energy Survey. Using the MCMF cluster confirmation algorithm, we identify optical counterparts, determine photometric redshifts and richnesses and assign a parameter, $f_{\rm cont}$, that reflects the probability that each SZE-optical pairing represents a random superposition of physically unassociated systems rather than a real cluster. The new PSZ-MCMF cluster catalogue consists of 853 MCMF confirmed clusters and has a purity of 90%. We present the properties of subsamples of the PSZ-MCMF catalogue that have purities ranging from 90% to 97.5%, depending on the adopted $f_{\rm cont}$ threshold. Halo mass estimates $M_{500}$, redshifts, richnesses, and optical centers are presented for all PSZ-MCMF clusters. The PSZ-MCMF catalogue adds 589 previously unknown Planck identified clusters over the DES footprint and provides redshifts for an additional 50 previously published Planck selected clusters with S/N>4.5. Using the subsample with spectroscopic redshifts, we demonstrate excellent cluster photo-$z$ performance with an RMS scatter in $\Delta z/(1+z)$ of 0.47%. Our MCMF based analysis allows us to infer the contamination fraction of the initial S/N>3 Planck selected candidate list, which is ~50%. We present a method of estimating the completeness of the PSZ-MCMF cluster sample. In comparison to the previously published Planck cluster catalogues. this new S/N>3 MCMF confirmed cluster catalogue populates the lower mass regime at all redshifts and includes clusters up to z$\sim$1.3., Comment: 20 pages, 5 Appendices, 17 figures. The newest submission matches the accepted MNRAS version

The three-dimensional correlation function offers an effective way to summarize the correlation of the large-scale structure even for imaging galaxy surveys. We have applied the projected three-dimensional correlation function, $\xi_{\rm p}$ to measure the Baryonic Acoustic Oscillations (BAO) scale on the first-three years Dark Energy Survey data. The sample consists of about 7 million galaxies in the redshift range $ 0.6 < z_{\rm p } < 1.1 $ over a footprint of $4108 \, \mathrm{deg}^2 $. Our theory modeling includes the impact of realistic true redshift distributions beyond Gaussian photo-$z$ approximation. To increase the signal-to-noise of the measurements, a Gaussian stacking window function is adopted in place of the commonly used top-hat. Using the full sample, $ D_{\rm M}(z_{\rm eff} ) / r_{\rm s} $, the ratio between the comoving angular diameter distance and the sound horizon, is constrained to be $ 19.00 \pm 0.67 $ (top-hat) and $ 19.15 \pm 0.58 $ (Gaussian) at $z_{\rm eff} = 0.835$. The constraint is weaker than the angular correlation $w$ constraint ($18.84 \pm 0.50$) because the BAO signals are heterogeneous across redshift. When a homogeneous BAO-signal sub-sample in the range $ 0.7 < z_{\rm p } < 1.0 $ ($z_{\rm eff} = 0.845$) is considered, $\xi_{\rm p} $ yields $ 19.80 \pm 0.67 $ (top-hat) and $ 19.84 \pm 0.53 $ (Gaussian). The latter is mildly stronger than the $w$ constraint ($19.86 \pm 0.55 $). We find that the $\xi_{\rm p} $ results are more sensitive to photo-$z$ errors than $w$ because $\xi_{\rm p}$ keeps the three-dimensional clustering information causing it to be more prone to photo-$z$ noise. The Gaussian window gives more robust results than the top-hat as the former is designed to suppress the low signal modes. $\xi_{\rm p}$ and the angular statistics such as $w$ have their own pros and cons, and they serve an important crosscheck with each other., Comment: 20 pages, 12 figures, minor changes to match published version

Reverberation mapping measurements have been used to constrain the relationship between the size of the broad-line region and luminosity of active galactic nuclei (AGN). This $R-L$ relation is used to estimate single-epoch virial black hole masses, and has been proposed for use to standardise AGN to determine cosmological distances. We present reverberation measurements made with H$\beta$ from the six-year Australian Dark Energy Survey (OzDES) Reverberation Mapping Program. We successfully recover reverberation lags for eight AGN at $0.12

We present a method for mapping variations between probability distribution functions and apply this method within the context of measuring galaxy redshift distributions from imaging survey data. This method, which we name PITPZ for the probability integral transformations it relies on, uses a difference in curves between distribution functions in an ensemble as a transformation to apply to another distribution function, thus transferring the variation in the ensemble to the latter distribution function. This procedure is broadly applicable to the problem of uncertainty propagation. In the context of redshift distributions, for example, the uncertainty contribution due to certain effects can be studied effectively only in simulations, thus necessitating a transfer of variation measured in simulations to the redshift distributions measured from data. We illustrate the use of PITPZ by using the method to propagate photometric calibration uncertainty to redshift distributions of the Dark Energy Survey Year 3 weak lensing source galaxies. For this test case, we find that PITPZ yields a lensing amplitude uncertainty estimate due to photometric calibration error within 1 per cent of the truth, compared to as much as a 30 per cent underestimate when using traditional methods.

We study the effect of magnification in the Dark Energy Survey Year 3 analysis of galaxy clustering and galaxy-galaxy lensing, using two different lens samples: a sample of Luminous red galaxies, redMaGiC, and a sample with a redshift-dependent magnitude limit, MagLim. We account for the effect of magnification on both the flux and size selection of galaxies, accounting for systematic effects using the Balrog image simulations. We estimate the impact of magnification on the galaxy clustering and galaxy-galaxy lensing cosmology analysis, finding it to be a significant systematic for the MagLim sample. We show cosmological constraints from the galaxy clustering auto-correlation and galaxy-galaxy lensing signal with different magnifications priors, finding broad consistency in cosmological parameters in $\Lambda$CDM and $w$CDM. However, when magnification bias amplitude is allowed to be free, we find the two-point correlations functions prefer a different amplitude to the fiducial input derived from the image simulations. We validate the magnification analysis by comparing the cross-clustering between lens bins with the prediction from the baseline analysis, which uses only the auto-correlation of the lens bins, indicating systematics other than magnification may be the cause of the discrepancy. We show adding the cross-clustering between lens redshift bins to the fit significantly improves the constraints on lens magnification parameters and allows uninformative priors to be used on magnification coefficients, without any loss of constraining power or prior volume concerns., Comment: Version accepted for publication in MNRAS. 21 pages, 13 figures, See this https://www.darkenergysurvey.org/des-year-3-cosmology-results-papers/ URL for the full DES Y3 cosmology release

We compare the two largest galaxy morphology catalogues, which separate early and late type galaxies at intermediate redshift. The two catalogues were built by applying supervised deep learning (convolutional neural networks, CNNs) to the Dark Energy Survey data down to a magnitude limit of $\sim$21 mag. The methodologies used for the construction of the catalogues include differences such as the cutout sizes, the labels used for training, and the input to the CNN - monochromatic images versus $gri$-band normalized images. In addition, one catalogue is trained using bright galaxies observed with DES ($i<18$), while the other is trained with bright galaxies ($r<17.5$) and `emulated' galaxies up to $r$-band magnitude $22.5$. Despite the different approaches, the agreement between the two catalogues is excellent up to $i<19$, demonstrating that CNN predictions are reliable for samples at least one magnitude fainter than the training sample limit. It also shows that morphological classifications based on monochromatic images are comparable to those based on $gri$-band images, at least in the bright regime. At fainter magnitudes, $i>19$, the overall agreement is good ($\sim$95\%), but is mostly driven by the large spiral fraction in the two catalogues. In contrast, the agreement within the elliptical population is not as good, especially at faint magnitudes. By studying the mismatched cases we are able to identify lenticular galaxies (at least up to $i<19$), which are difficult to distinguish using standard classification approaches. The synergy of both catalogues provides an unique opportunity to select a population of unusual galaxies., Comment: 17 pages, 14 figures (1 appendix for galaxy examples including 3 figures)

We present an alternative calibration of the MagLim lens sample redshift distributions from the Dark Energy Survey (DES) first three years of data (Y3). The new calibration is based on a combination of a Self-Organising Maps based scheme and clustering redshifts to estimate redshift distributions and inherent uncertainties, which is expected to be more accurate than the original DES Y3 redshift calibration of the lens sample. We describe in detail the methodology, we validate it on simulations and discuss the main effects dominating our error budget. The new calibration is in fair agreement with the fiducial DES Y3 redshift distributions calibration, with only mild differences ($<3\sigma$) in the means and widths of the distributions. We study the impact of this new calibration on cosmological constraints, analysing DES Y3 galaxy clustering and galaxy-galaxy lensing measurements, assuming a $\Lambda$CDM cosmology. We obtain $\Omega_{\rm m} = 0.30\pm 0.04$, $\sigma_8 = 0.81\pm 0.07 $ and $S_8 = 0.81\pm 0.04$, which implies a $\sim 0.4\sigma$ shift in the $\Omega_{\rm}-S_8$ plane compared to the fiducial DES Y3 results, highlighting the importance of the redshift calibration of the lens sample in multi-probe cosmological analyses.

We characterise the properties and evolution of Bright Central Galaxies (BCGs) and the surrounding intracluster light (ICL) in galaxy clusters identified in overlapping regions of the Dark Energy Survey and Atacama Cosmology Telescope Survey (DES-ACT), covering the redshift range $0.20$14.4. We also measure the stellar mass - halo mass (SMHM) relation for the BCG+ICL system and find that the slope, $\beta$, which characterises the dependence of $M_{\rm 200m,SZ}$ on the BCG+ICL stellar mass, increases with radius. The outskirts are more strongly correlated with the halo than the core, which supports that the BCG+ICL system follows a two-phase growth, where recent growth ($z<2$) occurs beyond the BCG's core. Additionally, we compare our observed SMHM relation results to the IllustrisTNG 300-1 cosmological hydrodynamic simulations and find moderate qualitative agreement in the amount of diffuse light. However, the SMHM relation's slope is steeper in TNG300-1 and the intrinsic scatter is lower, likely from the absence of projection effects in TNG300-1. Additionally, we find that the ICL exhibits a colour gradient such that the outskirts are bluer than the core. Moreover, for the lower halo mass clusters (log$_{10}(M_{\rm 200m,SZ}$/M$_{\odot})<$14.59 ), we detect a modest change in the colour gradient's slope with lookback time, which combined with the absence of stellar mass growth may suggest that lower mass clusters have been involved in growth via tidal stripping more recently than their higher mass counterparts., Comment: 19 pages, 12 Figures, 3 Tables. Submitted to MNRAS on 9/7/2022 -- Metadata typo corrected

Recent analyses have found intriguing correlations between the colour ($c$) of type Ia supernovae (SNe Ia) and the size of their 'mass-step', the relationship between SN Ia host galaxy stellar mass ($M_\mathrm{stellar}$) and SN Ia Hubble residual, and suggest that the cause of this relationship is dust. Using 675 photometrically-classified SNe Ia from the Dark Energy Survey 5-year sample, we study the differences in Hubble residual for a variety of global host galaxy and local environmental properties for SN Ia subsamples split by their colour. We find a $3\sigma$ difference in the mass-step when comparing blue ($c<0$) and red ($c>0$) SNe. We observe the lowest r.m.s. scatter ($\sim0.14$ mag) in the Hubble residual for blue SNe in low mass/blue environments, suggesting that this is the most homogeneous sample for cosmological analyses. By fitting for $c$-dependent relationships between Hubble residuals and $M_\mathrm{stellar}$, approximating existing dust models, we remove the mass-step from the data and find tentative $\sim 2\sigma$ residual steps in rest-frame galaxy $U-R$ colour. This indicates that dust modelling based on $M_\mathrm{stellar}$ may not fully explain the remaining dispersion in SN Ia luminosity. Instead, accounting for a $c$-dependent relationship between Hubble residuals and global $U-R$, results in $\leq1\sigma$ residual steps in $M_\mathrm{stellar}$ and local $U-R$, suggesting that $U-R$ provides different information about the environment of SNe Ia compared to $M_\mathrm{stellar}$, and motivating the inclusion of galaxy $U-R$ colour in SN Ia distance bias correction., Comment: 19 pages, 8 figures. Published in MNRAS

We infer the mean optical depth of a sample of optically-selected galaxy clusters from the Dark Energy Survey (DES) via the pairwise kinematic Sunyaev-Zel'dovich (kSZ) effect. The pairwise kSZ signal between pairs of clusters drawn from the DES Year-3 cluster catalog is detected at $4.1 \sigma$ in cosmic microwave background (CMB) temperature maps from two years of observations with the SPT-3G camera on the South Pole Telescope. After cuts, there are 24,580 clusters in the $\sim 1,400$ deg$^2$ of the southern sky observed by both experiments. We infer the mean optical depth of the cluster sample with two techniques. The optical depth inferred from the pairwise kSZ signal is $\bar{\tau}_e = (2.97 \pm 0.73) \times 10^{-3}$, while that inferred from the thermal SZ signal is $\bar{\tau}_e = (2.51 \pm 0.55^{\text{stat}} \pm 0.15^{\rm syst}) \times 10^{-3}$. The two measures agree at $0.6 \sigma$. We perform a suite of systematic checks to test the robustness of the analysis.

We present the luminosity functions and host galaxy properties of the Dark Energy Survey (DES) core-collapse supernova (CCSN) sample, consisting of 69 Type II and 50 Type Ibc spectroscopically and photometrically-confirmed supernovae over a redshift range $0.045

We constrain extensions to the $\Lambda$CDM model using measurements from the Dark Energy Survey's first three years of observations and external data. The DES data are the two-point correlation functions of weak gravitational lensing, galaxy clustering, and their cross-correlation. We use simulated data and blind analyses of real data to validate the robustness of our results. In many cases, constraining power is limited by the absence of nonlinear predictions that are reliable at our required precision. The models are: dark energy with a time-dependent equation of state, non-zero spatial curvature, sterile neutrinos, modifications of gravitational physics, and a binned $\sigma_8(z)$ model which serves as a probe of structure growth. For the time-varying dark energy equation of state evaluated at the pivot redshift we find $(w_{\rm p}, w_a)= (-0.99^{+0.28}_{-0.17},-0.9\pm 1.2)$ at 68% confidence with $z_{\rm p}=0.24$ from the DES measurements alone, and $(w_{\rm p}, w_a)= (-1.03^{+0.04}_{-0.03},-0.4^{+0.4}_{-0.3})$ with $z_{\rm p}=0.21$ for the combination of all data considered. Curvature constraints of $\Omega_k=0.0009\pm 0.0017$ and effective relativistic species $N_{\rm eff}=3.10^{+0.15}_{-0.16}$ are dominated by external data. For massive sterile neutrinos, we improve the upper bound on the mass $m_{\rm eff}$ by a factor of three compared to previous analyses, giving 95% limits of $(\Delta N_{\rm eff},m_{\rm eff})\leq (0.28, 0.20\, {\rm eV})$. We also constrain changes to the lensing and Poisson equations controlled by functions $\Sigma(k,z) = \Sigma_0 \Omega_{\Lambda}(z)/\Omega_{\Lambda,0}$ and $\mu(k,z)=\mu_0 \Omega_{\Lambda}(z)/\Omega_{\Lambda,0}$ respectively to $\Sigma_0=0.6^{+0.4}_{-0.5}$ from DES alone and $(\Sigma_0,\mu_0)=(0.04\pm 0.05,0.08^{+0.21}_{-0.19})$ for the combination of all data. Overall, we find no significant evidence for physics beyond $\Lambda$CDM., Comment: Updated to match published version and fix a citation reference. 46 pages, 25 figures, data available at https://dev.des.ncsa.illinois.edu/releases/y3a2/Y3key-extensions

We use stacked spectra of the host galaxies of photometrically identified type Ia supernovae (SNe Ia) from the Dark Energy Survey (DES) to search for correlations between Hubble diagram residuals and the spectral properties of the host galaxies. Utilising full spectrum fitting techniques on stacked spectra binned by Hubble residual, we find no evidence for trends between Hubble residuals and properties of the host galaxies that rely on spectral absorption features ($< 1.3\sigma$), such as stellar population age, metallicity, and mass-to-light ratio. However, we find significant trends between the Hubble residuals and the strengths of [OII] ($4.4\sigma$) and the Balmer emission lines ($3\sigma$). These trends are weaker than the well known trend between Hubble residuals and host galaxy stellar mass ($7.2\sigma$) that is derived from broad band photometry. After light curve corrections, we see fainter SNe Ia residing in galaxies with larger line strengths. We also find a trend (3$\sigma$) between Hubble residual and the Balmer decrement (a measure of reddening by dust) using H${\beta}$ and H${\gamma}$. The trend, quantified by correlation coefficients, is slightly more significant in the redder SNe Ia, suggesting that bluer SNe Ia are relatively unaffected by dust in the interstellar medium of the host and that dust contributes to current Hubble diagram scatter impacting the measurement of cosmological parameters., Comment: 15 pages, 10 figures. Accepted for publication in MNRAS

Gravitational time delays provide a powerful one step measurement of $H_0$, independent of all other probes. One key ingredient in time delay cosmography are high accuracy lens models. Those are currently expensive to obtain, both, in terms of computing and investigator time (10$^{5-6}$ CPU hours and $\sim$ 0.5-1 year, respectively). Major improvements in modeling speed are therefore necessary to exploit the large number of lenses that are forecast to be discovered over the current decade. In order to bypass this roadblock, building on the work by Shajib et al. (2019), we develop an automated modeling pipeline and apply it to a sample of 30 quadruply imaged quasars and one lensed compact galaxy, observed by the Hubble Space Telescope in multiple bands. Our automated pipeline can derive models for 30/31 lenses with few hours of human time and <100 CPU hours of computing time for a typical system. For each lens, we provide measurements of key parameters and predictions of magnification as well as time delays for the multiple images. We characterize the cosmography-readiness of our models using the stability of differences in Fermat potential (proportional to time delay) w.r.t. modeling choices. We find that for 10/30 lenses our models are cosmography or nearly cosmography grade (<3% and 3-5% variations). For 6/30 lenses the models are close to cosmography grade (5-10%). These results are based on informative priors and will need to be confirmed by further analysis. However, they are also likely to improve by extending the pipeline modeling sequence and options. In conclusion, we show that uniform cosmography grade modeling of large strong lens samples is within reach., Comment: 40 pages, 24 figures, 11 tables. Submitted to MNRAS

Gravitationally lensed supernovae (LSNe) are important probes of cosmic expansion, but they remain rare and difficult to find. Current cosmic surveys likely contain and 5-10 LSNe in total while next-generation experiments are expected to contain several hundreds to a few thousands of these systems. We search for these systems in observed Dark Energy Survey (DES) 5-year SN fields -- 10 3-sq. deg. regions of sky imaged in the $griz$ bands approximately every six nights over five years. To perform the search, we utilize the DeepZipper approach: a multi-branch deep learning architecture trained on image-level simulations of LSNe that simultaneously learns spatial and temporal relationships from time series of images. We find that our method obtains a LSN recall of 61.13% and a false positive rate of 0.02% on the DES SN field data. DeepZipper selected 2,245 candidates from a magnitude-limited ($m_i$ $<$ 22.5) catalog of 3,459,186 systems. We employ human visual inspection to review systems selected by the network and find three candidate LSNe in the DES SN fields., Comment: Accepted by ApJ

We present the second public data release (DR2) from the DECam Local Volume Exploration survey (DELVE). DELVE DR2 combines new DECam observations with archival DECam data from the Dark Energy Survey, the DECam Legacy Survey, and other DECam community programs. DELVE DR2 consists of ~160,000 exposures that cover >21,000 deg^2 of the high Galactic latitude (|b| > 10 deg) sky in four broadband optical/near-infrared filters (g, r, i, z). DELVE DR2 provides point-source and automatic aperture photometry for ~2.5 billion astronomical sources with a median 5{\sigma} point-source depth of g=24.3, r=23.9, i=23.5, and z=22.8 mag. A region of ~17,000 deg^2 has been imaged in all four filters, providing four-band photometric measurements for ~618 million astronomical sources. DELVE DR2 covers more than four times the area of the previous DELVE data release and contains roughly five times as many astronomical objects. DELVE DR2 is publicly available via the NOIRLab Astro Data Lab science platform., Comment: 23 pages, 10 figures, 5 tables; to be submitted to AAS Journals; public data release at https://datalab.noirlab.edu/delve/. arXiv admin note: substantial text overlap with arXiv:2103.07476

Cross-correlations of galaxy positions and galaxy shears with maps of gravitational lensing of the cosmic microwave background (CMB) are sensitive to the distribution of large-scale structure in the Universe. Such cross-correlations are also expected to be immune to some of the systematic effects that complicate correlation measurements internal to galaxy surveys. We present measurements and modeling of the cross-correlations between galaxy positions and galaxy lensing measured in the first three years of data from the Dark Energy Survey with CMB lensing maps derived from a combination of data from the 2500 deg$^2$ SPT-SZ survey conducted with the South Pole Telescope and full-sky data from the Planck satellite. The CMB lensing maps used in this analysis have been constructed in a way that minimizes biases from the thermal Sunyaev Zel'dovich effect, making them well suited for cross-correlation studies. The total signal-to-noise of the cross-correlation measurements is 23.9 (25.7) when using a choice of angular scales optimized for a linear (nonlinear) galaxy bias model. We use the cross-correlation measurements to obtain constraints on cosmological parameters. For our fiducial galaxy sample, which consist of four bins of magnitude-selected galaxies, we find constraints of $\Omega_{m} = 0.272^{+0.032}_{-0.052}$ and $S_{8} \equiv \sigma_8 \sqrt{\Omega_{m}/0.3}= 0.736^{+0.032}_{-0.028}$ ($\Omega_{m} = 0.245^{+0.026}_{-0.044}$ and $S_{8} = 0.734^{+0.035}_{-0.028}$) when assuming linear (nonlinear) galaxy bias in our modeling. Considering only the cross-correlation of galaxy shear with CMB lensing, we find $\Omega_{m} = 0.270^{+0.043}_{-0.061}$ and $S_{8} = 0.740^{+0.034}_{-0.029}$. Our constraints on $S_8$ are consistent with recent cosmic shear measurements, but lower than the values preferred by primary CMB measurements from Planck., Comment: 25 pages, 19 figures, submitted to PRD

We present cosmological constraints from the analysis of angular power spectra of cosmic shear maps based on data from the first three years of observations by the Dark Energy Survey (DES Y3). Our measurements are based on the pseudo-$C_\ell$ method and offer a view complementary to that of the two-point correlation functions in real space, as the two estimators are known to compress and select Gaussian information in different ways, due to scale cuts. They may also be differently affected by systematic effects and theoretical uncertainties, such as baryons and intrinsic alignments (IA), making this analysis an important cross-check. In the context of $\Lambda$CDM, and using the same fiducial model as in the DES Y3 real space analysis, we find ${S_8 \equiv \sigma_8 \sqrt{\Omega_{\rm m}/0.3} = 0.793^{+0.038}_{-0.025}}$, which further improves to ${S_8 = 0.784\pm 0.026 }$ when including shear ratios. This constraint is within expected statistical fluctuations from the real space analysis, and in agreement with DES~Y3 analyses of non-Gaussian statistics, but favors a slightly higher value of $S_8$, which reduces the tension with the Planck cosmic microwave background 2018 results from $2.3\sigma$ in the real space analysis to $1.5\sigma$ in this work. We explore less conservative IA models than the one adopted in our fiducial analysis, finding no clear preference for a more complex model. We also include small scales, using an increased Fourier mode cut-off up to $k_{\rm max}={5}{h{\rm Mpc}^{-1}}$, which allows to constrain baryonic feedback while leaving cosmological constraints essentially unchanged. Finally, we present an approximate reconstruction of the linear matter power spectrum at present time, which is found to be about 20\% lower than predicted by Planck 2018, as reflected by the $1.5\sigma$ lower $S_8$ value.

Current and future cosmological analyses with Type Ia Supernovae (SNe Ia) face three critical challenges: i) measuring redshifts from the supernova or its host galaxy; ii) classifying SNe without spectra; and iii) accounting for correlations between the properties of SNe Ia and their host galaxies. We present here a novel approach that addresses each challenge. In the context of the Dark Energy Survey (DES), we analyze a SNIa sample with host galaxies in the redMaGiC galaxy catalog, a selection of Luminous Red Galaxies. Photo-$z$ estimates for these galaxies are expected to be accurate to $\sigma_{\Delta z/(1+z)}\sim0.02$. The DES-5YR photometrically classified SNIa sample contains approximately 1600 SNe and 125 of these SNe are in redMaGiC galaxies. We demonstrate that redMaGiC galaxies almost exclusively host SNe Ia, reducing concerns with classification uncertainties. With this subsample, we find similar Hubble scatter (to within $\sim0.01$ mag) using photometric redshifts in place of spectroscopic redshifts. With detailed simulations, we show the bias due to using photo-$z$s from redMaGiC host galaxies on the measurement of the dark energy equation-of-state $w$ is up to $\Delta w \sim 0.01-0.02$. With real data, we measure a difference in $w$ when using redMaGiC photometric redshifts versus spectroscopic redshifts of $\Delta w = 0.005$. Finally, we discuss how SNe in redMaGiC galaxies appear to be a more standardizable population due to a weaker relation between color and luminosity ($\beta$) compared to the DES-3YR population by $\sim5\sigma$; this finding is consistent with predictions that redMaGiC galaxies exhibit lower reddening ratios ($\textrm{R}_\textrm{V}$) than the general population of SN host galaxies. These results establish the feasibility of performing redMaGiC SN cosmology with photometric survey data in the absence of spectroscopic data., Comment: 15 pages, 8 figures. Submitting to ApJ, comments welcome

We evaluate the consistency between lensing and clustering probes of large-scale structure based on measurements of projected galaxy clustering from BOSS combined with overlapping galaxy-galaxy lensing from three surveys: DES Y3, HSC Y1, and KiDS-1000. An intra-lensing-survey study finds good agreement between these lensing data. We model the observations using the Dark Emulator and fit the data at two fixed cosmologies: Planck, with $S_8=0.83$, and a Lensing cosmology with $S_8=0.76$. For a joint analysis limited to scales with $R>5.25h^{-1}$Mpc, we find that both cosmologies provide an acceptable fit to the data. Full utilisation of the small-scale clustering and lensing measurements is hindered by uncertainty in the impact of baryon feedback and assembly bias, which we account for with a reasoned theoretical error budget. We incorporate a systematic scaling parameter for each redshift bin, $A$, that decouples the lensing and clustering to capture any inconsistency. When a wide range of scales ($0.15

We use a recent census of the Milky Way (MW) satellite galaxy population to constrain the lifetime of particle dark matter (DM). We consider two-body decaying dark matter (DDM) in which a heavy DM particle decays with lifetime $\tau$ comparable to the age of the Universe to a lighter DM particle (with mass splitting $\epsilon$) and to a dark radiation species. These decays impart a characteristic "kick velocity," $V_{\mathrm{kick}}=\epsilon c$, on the DM daughter particles, significantly depleting the DM content of low-mass subhalos and making them more susceptible to tidal disruption. We fit the suppression of the present-day DDM subhalo mass function (SHMF) as a function of $\tau$ and $V_{\mathrm{kick}}$ using a suite of high-resolution zoom-in simulations of MW-mass halos, and we validate this model on new DDM simulations of systems specifically chosen to resemble the MW. We implement our DDM SHMF predictions in a forward model that incorporates inhomogeneities in the spatial distribution and detectability of MW satellites and uncertainties in the mapping between galaxies and DM halos, the properties of the MW system, and the disruption of subhalos by the MW disk using an empirical model for the galaxy--halo connection. By comparing to the observed MW satellite population, we conservatively exclude DDM models with $\tau < 18\ \mathrm{Gyr}$ ($29\ \mathrm{Gyr}$) for $V_{\mathrm{kick}}=20\ \mathrm{km}\, \mathrm{s}^{-1}$ ($40\ \mathrm{km}\, \mathrm{s}^{-1}$) at $95\%$ confidence. These constraints are among the most stringent and robust small-scale structure limits on the DM particle lifetime and strongly disfavor DDM models that have been proposed to alleviate the Hubble and $S_8$ tensions., Comment: 26 pages, 11 figures, 1 table. Updated to published version

We present high signal-to-noise measurements of three-point shear correlations and the third moment of the mass aperture statistic using the first 3 years of data from the Dark Energy Survey. We additionally obtain the first measurements of the configuration and scale dependence of the four three-point shear correlations which carry cosmological information. With the third-order mass aperture statistic, we present tomographic measurements over angular scales of 4 to 60 arcminutes with a combined statistical significance of 15.0$\sigma$. Using the tomographic information and measuring also the second-order mass aperture, we additionally obtain a skewness parameter and its redshift evolution. We find that the amplitudes and scale-dependence of these shear 3pt functions are in qualitative agreement with measurements in a mock galaxy catalog based on N-body simulations, indicating promise for including them in future cosmological analyses. We validate our measurements by showing that B-modes, parity-violating contributions and PSF modeling uncertainties are negligible, and determine that the measured signals are likely to be of astrophysical and gravitational origin., Comment: Minor changes. Matches version published in PRD

We study the optical $gri$ photometric variability of a sample of 190 quasars within the SDSS Stripe 82 region that have long-term photometric coverage during $\sim 1998-2020$ with SDSS, PanSTARRS-1, the Dark Energy Survey, and dedicated follow-up monitoring with Blanco 4m/DECam. With on average $\sim 200$ nightly epochs per quasar per filter band, we improve the parameter constraints from a Damped Random Walk (DRW) model fit to the light curves over previous studies with 10-15 yr baselines and $\lesssim 100$ epochs. We find that the average damping timescale $\tau_{\rm DRW}$ continues to rise with increased baseline, reaching a median value of $\sim 750$ days ($g$ band) in the rest-frame of these quasars using the 20-yr light curves. Some quasars may have gradual, long-term trends in their light curves, suggesting that either the DRW fit requires very long baselines to converge, or that the underlying variability is more complex than a single DRW process for these quasars. Using a subset of quasars with better-constrained $\tau_{\rm DRW}$ (less than 20\% of the baseline), we confirm a weak wavelength dependence of $\tau_{\rm DRW}\propto \lambda^{0.51\pm0.20}$. We further quantify optical variability of these quasars over days to decades timescales using structure function (SF) and power spectrum density (PSD) analyses. The SF and PSD measurements qualitatively confirm the measured (hundreds of days) damping timescales from the DRW fits. However, the ensemble PSD is steeper than that of a DRW on timescales less than $\sim$ a month for these luminous quasars, and this second break point correlates with the longer DRW damping timescale., Comment: 20 pages, 24 figures

Lensing Without Borders is a cross-survey collaboration created to assess the consistency of galaxy-galaxy lensing signals ($\Delta\Sigma$) across different data-sets and to carry out end-to-end tests of systematic errors. We perform a blind comparison of the amplitude of $\Delta\Sigma$ using lens samples from BOSS and six independent lensing surveys. We find good agreement between empirically estimated and reported systematic errors which agree to better than 2.3$\sigma$ in four lens bins and three radial ranges. For lenses with $z_{\rm L}>0.43$ and considering statistical errors, we detect a 3-4$\sigma$ correlation between lensing amplitude and survey depth. This correlation could arise from the increasing impact at higher redshift of unrecognised galaxy blends on shear calibration and imperfections in photometric redshift calibration. At $z_{\rm L}>0.54$ amplitudes may additionally correlate with foreground stellar density. The amplitude of these trends is within survey-defined systematic error budgets which are designed to include known shear and redshift calibration uncertainty. Using a fully empirical and conservative method, we do not find evidence for large unknown systematics. Systematic errors greater than 15% (25%) ruled out in three lens bins at 68% (95%) confidence at $z<0.54$. Differences with respect to predictions based on clustering are observed to be at the 20-30% level. Our results therefore suggest that lensing systematics alone are unlikely to fully explain the "lensing is low" effect at $z<0.54$. This analysis demonstrates the power of cross-survey comparisons and provides a promising path for identifying and reducing systematics in future lensing analyses., Comment: 41 page, 20 figures

We perform a cosmic shear analysis in harmonic space using the first year of data collected by the Dark Energy Survey (DES-Y1). We measure the cosmic weak lensing shear power spectra using the Metacalibration catalogue and perform a likelihood analysis within the framework of CosmoSIS. We set scale cuts based on baryonic effects contamination and model redshift and shear calibration uncertainties as well as intrinsic alignments. We adopt as fiducial covariance matrix an analytical computation accounting for the mask geometry in the Gaussian term, including non-Gaussian contributions. A suite of 1200 lognormal simulations is used to validate the harmonic space pipeline and the covariance matrix. We perform a series of stress tests to gauge the robustness of the harmonic space analysis. Finally, we use the DES-Y1 pipeline in configuration space to perform a similar likelihood analysis and compare both results, demonstrating their compatibility in estimating the cosmological parameters $S_8$, $\sigma_8$ and $\Omega_m$. The methods implemented and validated in this paper will allow us to perform a consistent harmonic space analysis in the upcoming DES data., Comment: 16 pages, 11 figures. This version matches the published one

We use six years of data from the Dark Energy Survey to perform a detailed photometric characterization of the Phoenix stellar stream, a 15-degree long, thin, dynamically cold, low-metallicity stellar system in the southern hemisphere. We use natural splines, a non-parametric modeling technique, to simultaneously fit the stream track, width, and linear density. This updated stream model allows us to improve measurements of the heliocentric distance ($17.4 \pm 0.1\,{\rm (stat.)} \pm 0.8\,{\rm (sys.)}$ kpc) and distance gradient ($-0.009 \pm 0.006$ kpc deg$^{-1}$) of Phoenix, which corresponds to a small change of $0.13 \pm 0.09$ kpc in heliocentric distance along the length of the stream. We measure linear intensity variations on degree scales, as well as deviations in the stream track on $\sim 2$-degree scales, suggesting that the stream may have been disturbed during its formation and/or evolution. We recover three peaks and one gap in linear intensity along with fluctuations in the stream track. Compared to other thin streams, the Phoenix stream shows more fluctuations and, consequently, the study of Phoenix offers a unique perspective on gravitational perturbations of stellar streams. We discuss possible sources of perturbations to Phoenix including baryonic structures in the Galaxy and dark matter subhalos., Comment: 18 Pages, 7 Figures, Matches published version

We report the combined results of eight searches for strong gravitational lens systems in the full 5,000 sq. deg. of Dark Energy Survey (DES) observations. The observations accumulated by the end of the third observing season fully covered the DES footprint in 5 filters (grizY), with an $i-$band limiting magnitude (at $10\sigma$) of 23.44. In four searches, a list of potential candidates was identified using a color and magnitude selection from the object catalogs created from the first three observing seasons. Three other searches were conducted at the locations of previously identified galaxy clusters. Cutout images of potential candidates were then visually scanned using an object viewer. An additional set of candidates came from a data-quality check of a subset of the color-coadd "tiles" created from the full DES six-season data set. A short list of the most promising strong lens candidates was then numerically ranked according to whether or not we judged them to be bona fide strong gravitational lens systems. These searches discovered a diverse set of 247 strong lens candidate systems, of which 81 are identified for the first time. We provide the coordinates, magnitudes, and photometric properties of the lens and source objects, and an estimate of the Einstein radius for 81 new systems and 166 previously reported. This catalog will be of use for selecting interesting systems for detailed follow-up, studies of galaxy cluster and group mass profiles, as well as a training/validation set for automated strong lens searches., Comment: 38 pages, 17 figures, 4 tables, accepted by ApJS