Researcher profile

D. W. Gerdes

D. W. Gerdes contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 13 - UnverifiedVerification L1Unclaimed author
53works
0followers
1topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

53 published item(s)

preprint2022arXiv

A DECam Search for Explosive Optical Transients Associated with IceCube Neutrinos

In this work, we investigate the likelihood of association between realtime, TeV-PeV energy neutrino alerts from IceCube and optical counterparts in the form of core-collapse supernovae (CC SNe). The optical follow-up of IceCube alerts requires two main instrumental capabilities: (1) deep imaging, since 73\% of neutrinos would come from CC SNe at redshifts $z > 0.3$, and (2) a large field of view (FoV), since typical IceCube muon neutrino pointing accuracy is on the order of $\sim1$~deg. With Blanco/DECam ($gri$ to 24th magnitude and $2.2$~deg diameter FoV), we performed a triggered optical follow-up observation of two IceCube alerts, IC170922A and IC171106A on $\sim6$~nights during the $\sim3$~weeks following each alert. For the IC170922A (IC171106A) follow-up observations, we expect that 12.1\% (9.5\%) of coincident CC SNe at $z \lesssim 0.3$ are detectable, and that on average, 0.23 (0.07) unassociated SNe in the neutrino 90\% containment regions also pass our selection criteria. We find two candidate CC SNe that are temporally coincident with the neutrino alerts in the FoV, but none in the 90\% containment regions, which is statistically consistent with expected rates of background CC SNe for these observations. If CC SNe are the dominant source of TeV-PeV neutrinos, we would expect an excess of coincident CC SNe to be detectable at the $3σ$ confidence level using DECam observations similar to those of this work for $\sim60$ ($\sim200$) neutrino alerts with (without) redshift information for all candidates.

preprint2022arXiv

Constraints on the Physical Properties of GW190814 through Simulations based on DECam Follow-up Observations by the Dark Energy Survey

On 14 August 2019, the LIGO and Virgo Collaborations detected gravitational waves from a black hole and a 2.6 solar mass compact object, possibly the first neutron star -- black hole (NSBH) merger. In search of an optical counterpart, the Dark Energy Survey (DES) obtained deep imaging of the entire 90 percent confidence level localization area with Blanco/DECam 0, 1, 2, 3, 6, and 16 nights after the merger. Objects with varying brightness were detected by the DES Pipeline and we systematically reduced the candidate counterparts through catalog matching, light curve properties, host-galaxy photometric redshifts, SOAR spectroscopic follow-up observations, and machine-learning-based photometric classification. All candidates were rejected as counterparts to the merger. To quantify the sensitivity of our search, we applied our selection criteria to full light curve simulations of supernovae and kilonovae as they would appear in the DECam observations. Since the source class of the merger was uncertain, we utilized an agnostic, three-component kilonova model based on tidally-disrupted NS ejecta properties to quantify our detection efficiency of a counterpart if the merger included a NS. We find that if a kilonova occurred during this merger, configurations where the ejected matter is greater than 0.07 solar masses, has lanthanide abundance less than $10^{-8.56}$, and has a velocity between $0.18c$ and $0.21c$ are disfavored at the $2σ$ level. Furthermore, we estimate that our background reduction methods are capable of associating gravitational wave signals with a detected electromagnetic counterpart at the $4σ$ level in $95\%$ of future follow-up observations.

preprint2022arXiv

Cosmological Constraints from DES Y1 Cluster Abundances and SPT Multi-wavelength data

We perform a joint analysis of the counts of redMaPPer clusters selected from the Dark Energy Survey (DES) Y1 data and multi-wavelength follow-up data collected within the 2500 deg$^2$ South Pole Telescope (SPT) SZ survey. The SPT follow-up data, calibrating the richness--mass relation of the optically selected redMaPPer catalog, enable the cosmological exploitation of the DES cluster abundance data. To explore possible systematics related to the modeling of projection effects, we consider two calibrations of the observational scatter on richness estimates: a simple Gaussian model which account only for the background contamination (BKG), and a model which further includes contamination and incompleteness due to projection effects (PRJ). Assuming either a $Λ$CDM+$\sum m_ν$ or $w$CDM+$\sum m_ν$ cosmology, and for both scatter models, we derive cosmological constraints consistent with multiple cosmological probes of the low and high redshift Universe, and in particular with the SPT cluster abundance data. This result demonstrates that the DES Y1 and SPT cluster counts provide consistent cosmological constraints, if the same mass calibration data set is adopted. It thus supports the conclusion of the DES Y1 cluster cosmology analysis which interprets the tension observed with other cosmological probes in terms of systematics affecting the stacked weak lensing analysis of optically--selected low--richness clusters. Finally, we analyse the first combined optically-SZ selected cluster catalogue obtained by including the SPT sample above the maximum redshift probed by the DES Y1 redMaPPer sample. Besides providing a mild improvement of the cosmological constraints, this data combination serves as a stricter test of our scatter models: the PRJ model, providing scaling relations consistent between the two abundance and multi-wavelength follow-up data, is favored over the BKG model.

preprint2022arXiv

Dark Energy Survey Year 3 Results: A 2.7% measurement of Baryon Acoustic Oscillation distance scale at redshift 0.835

We present angular diameter measurements obtained by measuring the position of Baryon Acoustic Oscillations (BAO) in an optimised sample of galaxies from the first three years of Dark Energy Survey data (DES Y3). The sample consists of 7 million galaxies distributed over a footprint of 4100 deg$^2$ with $0.6 < z_{\rm photo} < 1.1$ and a typical redshift uncertainty of $0.03(1+z)$. The sample selection is the same as in the BAO measurement with the first year of DES data, but the analysis presented here uses three times the area, extends to higher redshift and makes a number of improvements, including a fully analytical BAO template, the use of covariances from both theory and simulations, and an extensive pre-unblinding protocol. We used two different statistics: angular correlation function and power spectrum, and validate our pipeline with an ensemble of over 1500 realistic simulations. Both statistics yield compatible results. We combine the likelihoods derived from angular correlations and spherical harmonics to constrain the ratio of comoving angular diameter distance $D_M$ at the effective redshift of our sample to the sound horizon scale at the drag epoch. We obtain $D_M(z_{\rm eff}=0.835)/r_{\rm d} = 18.92 \pm 0.51$, which is consistent with, but smaller than, the Planck prediction assuming flat \lcdm, at the level of $2.3 σ$. The analysis was performed blind and is robust to changes in a number of analysis choices. It represents the most precise BAO distance measurement from imaging data to date, and is competitive with the latest transverse ones from spectroscopic samples at $z>0.75$. When combined with DES 3x2pt + SNIa, they lead to improvements in $H_0$ and $Ω_m$ constraints by $\sim 20\%$

preprint2022arXiv

Dark Energy Survey Year 3 Results: Cosmological Constraints from Galaxy Clustering and Weak Lensing

We present the first cosmology results from large-scale structure in the Dark Energy Survey (DES) spanning 5000 deg$^2$. We perform an analysis combining three two-point correlation functions (3$\times$2pt): (i) cosmic shear using 100 million source galaxies, (ii) galaxy clustering, and (iii) the cross-correlation of source galaxy shear with lens galaxy positions. The analysis was designed to mitigate confirmation or observer bias; we describe specific changes made to the lens galaxy sample following unblinding of the results. We model the data within the flat $Λ$CDM and $w$CDM cosmological models. We find consistent cosmological results between the three two-point correlation functions; their combination yields clustering amplitude $S_8=0.776^{+0.017}_{-0.017}$ and matter density $Ω_{\mathrm{m}} = 0.339^{+0.032}_{-0.031}$ in $Λ$CDM, mean with 68% confidence limits; $S_8=0.775^{+0.026}_{-0.024}$, $Ω_{\mathrm{m}} = 0.352^{+0.035}_{-0.041}$, and dark energy equation-of-state parameter $w=-0.98^{+0.32}_{-0.20}$ in $w$CDM. This combination of DES data is consistent with the prediction of the model favored by the Planck 2018 cosmic microwave background (CMB) primary anisotropy data, which is quantified with a probability-to-exceed $p=0.13$ to $0.48$. When combining DES 3$\times$2pt data with available baryon acoustic oscillation, redshift-space distortion, and type Ia supernovae data, we find $p=0.34$. Combining all of these data sets with Planck CMB lensing yields joint parameter constraints of $S_8 = 0.812^{+0.008}_{-0.008}$, $Ω_{\mathrm{m}} = 0.306^{+0.004}_{-0.005}$, $h=0.680^{+0.004}_{-0.003}$, and $\sum m_ν<0.13 \;\mathrm{eV\; (95\% \;CL)}$ in $Λ$CDM; $S_8 = 0.812^{+0.008}_{-0.008}$, $Ω_{\mathrm{m}} = 0.302^{+0.006}_{-0.006}$, $h=0.687^{+0.006}_{-0.007}$, and $w=-1.031^{+0.030}_{-0.027}$ in $w$CDM. (abridged)

preprint2022arXiv

Dark Energy Survey Year 3 results: cosmological constraints from the analysis of cosmic shear in harmonic space

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 $Λ$CDM, and using the same fiducial model as in the DES Y3 real space analysis, we find ${S_8 \equiv σ_8 \sqrt{Ω_{\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σ$ in the real space analysis to $1.5σ$ 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σ$ lower $S_8$ value.

preprint2022arXiv

Dark Energy Survey Year 3 Results: Cosmology from Cosmic Shear and Robustness to Modeling Uncertainty

This work and its companion paper, Amon et al. (2021), present cosmic shear measurements and cosmological constraints from over 100 million source galaxies in the Dark Energy Survey (DES) Year 3 data. We constrain the lensing amplitude parameter $S_8\equivσ_8\sqrt{Ω_\textrm{m}/0.3}$ at the 3% level in $Λ$CDM: $S_8=0.759^{+0.025}_{-0.023}$ (68% CL). Our constraint is at the 2% level when using angular scale cuts that are optimized for the $Λ$CDM analysis: $S_8=0.772^{+0.018}_{-0.017}$ (68% CL). With cosmic shear alone, we find no statistically significant constraint on the dark energy equation-of-state parameter at our present statistical power. We carry out our analysis blind, and compare our measurement with constraints from two other contemporary weak-lensing experiments: the Kilo-Degree Survey (KiDS) and Hyper-Suprime Camera Subaru Strategic Program (HSC). We additionally quantify the agreement between our data and external constraints from the Cosmic Microwave Background (CMB). Our DES Y3 result under the assumption of $Λ$CDM is found to be in statistical agreement with Planck 2018, although favors a lower $S_8$ than the CMB-inferred value by $2.3σ$ (a $p$-value of 0.02). This paper explores the robustness of these cosmic shear results to modeling of intrinsic alignments, the matter power spectrum and baryonic physics. We additionally explore the statistical preference of our data for intrinsic alignment models of different complexity. The fiducial cosmic shear model is tested using synthetic data, and we report no biases greater than 0.3$σ$ in the plane of $S_8\timesΩ_\textrm{m}$ caused by uncertainties in the theoretical models.

preprint2022arXiv

Dark Energy Survey Year 3 Results: Deep Field Optical + Near-Infrared Images and Catalogue

We describe the Dark Energy Survey (DES) Deep Fields, a set of images and associated multi-wavelength catalogue ($ugrizJHKs$) built from Dark Energy Camera (DECam) and Visible and Infrared Survey Telescope for Astronomy (VISTA) data. The DES Deep Fields comprise 11 fields (10 DES supernova fields plus COSMOS), with a total area of $\sim30~$ square degrees in $ugriz$ bands and reaching a maximum $i$-band depth of 26.75 (AB, $10σ$, 2 arcsec). We present a catalogue for the DES 3-year cosmology analysis of those four fields with full 8-band coverage, totalling $5.88~$ sq. deg. after masking. Numbering $2.8~$million objects ($1.6~$million post masking), our catalogue is drawn from images coadded to consistent depths of $r=25.7, i=25, z=24.3$ mag. We use a new model-fitting code, built upon established methods, to deblend sources and ensure consistent colours across the $u$-band to $Ks$-band wavelength range. We further detail the tight control we maintain over the point-spread function modelling required for the model fitting, astrometry and consistency of photometry between the four fields. The catalogue allows us to perform a careful star-galaxy separation and produces excellent photometric redshift performance (${\rm NMAD} = 0.023$ at $i<23$). The Deep-Fields catalogue will be made available as part of the cosmology data products release, following the completion of the DES 3-year weak lensing and galaxy clustering cosmology work.

preprint2022arXiv

Dark Energy Survey Year 3 Results: Exploiting small-scale information with lensing shear ratios

Using the first three years of data from the Dark Energy Survey, we use ratios of small-scale galaxy-galaxy lensing measurements around the same lens sample to constrain source redshift uncertainties, intrinsic alignments and other nuisance parameters of our model. Instead of using a simple geometric approach for the ratios, we use the full modeling of the galaxy-galaxy lensing measurements, including the corresponding integration over the power spectrum and the contributions from intrinsic alignments and lens magnification. We perform extensive testing of the small-scale shear ratio (SR) modeling by studying the impact of different effects such as the inclusion of baryonic physics, non-linear biasing, halo occupation distribution descriptions and lens magnification, among others, and using realistic $N$-body simulations. We validate the robustness of our constraints in the data by using two independent lens samples, and by deriving constraints using the corresponding large-scale ratios for which the modeling is simpler. The DES Y3 results demonstrate how the ratios provide significant improvements in constraining power for several nuisance parameters in our model, especially on source redshift calibration and intrinsic alignments (IA). For source redshifts, SR improves the constraints from the prior by up to 38\% in some redshift bins. Such improvements, and especially the constraints it provides on IA, translate to tighter cosmological constraints when SR is combined with cosmic shear and other 2pt functions. In particular, for the DES Y3 data, SR improves $S_8$ constraints from cosmic shear by up to 31\%, and for the full combination of probes (3$\times$2pt) by up to 10\%. The shear ratios presented in this work are used as an additional likelihood for cosmic shear, 2$\times$2pt and the full 3$\times$2pt in the fiducial DES Y3 cosmological analysis.

preprint2022arXiv

Dark Energy Survey Year 3 Results: High-precision measurement and modeling of galaxy-galaxy lensing

We present and characterize the galaxy-galaxy lensing signal measured using the first three years of data from the Dark Energy Survey (DES Y3) covering 4132 deg$^2$. These galaxy-galaxy measurements are used in the DES Y3 3$\times$2pt cosmological analysis, which combines weak lensing and galaxy clustering information. We use two lens samples: a magnitude-limited sample and the redMaGic sample, which span the redshift range $\sim 0.2-1$ with 10.7 M and 2.6 M galaxies respectively. For the source catalog, we use the Metacalibration shape sample, consisting of $\simeq$100 M galaxies separated into 4 tomographic bins. Our galaxy-galaxy lensing estimator is the mean tangential shear, for which we obtain a total S/N of $\sim$148 for MagLim ($\sim$120 for redMaGic), and $\sim$67 ($\sim$55) after applying the scale cuts of 6 Mpc/$h$. Thus we reach percent-level statistical precision, which requires that our modeling and systematic-error control be of comparable accuracy. The tangential shear model used in the 3$\times$2pt cosmological analysis includes lens magnification, a five-parameter intrinsic alignment model (TATT), marginalization over a point-mass to remove information from small scales and a linear galaxy bias model validated with higher-order terms. We explore the impact of these choices on the tangential shear observable and study the significance of effects not included in our model, such as reduced shear, source magnification and source clustering. We also test the robustness of our measurements to various observational and systematics effects, such as the impact of observing conditions, lens-source clustering, random-point subtraction, scale-dependent Metacalibration responses, PSF residuals, and B-modes.

preprint2022arXiv

Dark Energy Survey Year 3 Results: Measuring the Survey Transfer Function with Balrog

We describe an updated calibration and diagnostic framework, Balrog, used to directly sample the selection and photometric biases of the Dark Energy Survey&#39;s (DES) Year 3 (Y3) dataset. We systematically inject onto the single-epoch images of a random 20% subset of the DES footprint an ensemble of nearly 30 million realistic galaxy models derived from DES Deep Field observations. These augmented images are analyzed in parallel with the original data to automatically inherit measurement systematics that are often too difficult to capture with traditional generative models. The resulting object catalog is a Monte Carlo sampling of the DES transfer function and is used as a powerful diagnostic and calibration tool for a variety of DES Y3 science, particularly for the calibration of the photometric redshifts of distant &#34;source&#34; galaxies and magnification biases of nearer &#34;lens&#34; galaxies. The recovered Balrog injections are shown to closely match the photometric property distributions of the Y3 GOLD catalog, particularly in color, and capture the number density fluctuations from observing conditions of the real data within 1% for a typical galaxy sample. We find that Y3 colors are extremely well calibrated, typically within ~1-8 millimagnitudes, but for a small subset of objects we detect significant magnitude biases correlated with large overestimates of the injected object size due to proximity effects and blending. We discuss approaches to extend the current methodology to capture more aspects of the transfer function and reach full coverage of the survey footprint for future analyses.

preprint2022arXiv

Dark Energy Survey Year 3 Results: Point-Spread Function Modeling

We introduce a new software package for modeling the point-spread function (PSF) of astronomical images, called Piff (PSFs In the Full FOV), which we apply to the first three years (known as Y3) of the Dark Energy Survey (DES) data. We describe the relevant details about the algorithms used by Piff to model the PSF, including how the PSF model varies across the field of view (FOV). Diagnostic results show that the systematic errors from the PSF modeling are very small over the range of scales that are important for the DES Y3 weak lensing analysis. In particular, the systematic errors from the PSF modeling are significantly smaller than the corresponding results from the DES year one (Y1) analysis. We also briefly describe some planned improvements to Piff that we expect to further reduce the modeling errors in future analyses.

preprint2022arXiv

Dark Energy Survey Year 3 Results: Weak Lensing Shape Catalogue

We present and characterise the galaxy shape catalogue from the first 3 years of Dark Energy Survey (DES) observations, over an effective area of ~4143 deg$^2$ of the southern sky. We describe our data analysis process and our self-calibrating shear measurement pipeline METACALIBRATION, which builds and improves upon the pipeline used in the DES Year 1 analysis in several aspects. The DES Year 3 weak-lensing shape catalogue consists of 100,204,026 galaxies, measured in the $riz$ bands, resulting in a weighted source number density of $n_{\rm eff} = 5.59$ gal/arcmin$ ^{2}$ and corresponding shape noise $σ_e = 0.261$. We perform a battery of internal null tests on the catalogue, including tests on systematics related to the point-spread function (PSF) modelling, spurious catalogue B-mode signals, catalogue contamination, and galaxy properties.

preprint2022arXiv

DES Y3 results: Blending shear and redshift biases in image simulations

As the statistical power of galaxy weak lensing reaches percent level precision, large, realistic and robust simulations are required to calibrate observational systematics, especially given the increased importance of object blending as survey depths increase. To capture the coupled effects of blending in both shear and photometric redshift calibration, we define the effective redshift distribution for lensing, $n_γ(z)$, and describe how to estimate it using image simulations. We use an extensive suite of tailored image simulations to characterize the performance of the shear estimation pipeline applied to the Dark Energy Survey (DES) Year 3 dataset. We describe the multi-band, multi-epoch simulations, and demonstrate their high level of realism through comparisons to the real DES data. We isolate the effects that generate shear calibration biases by running variations on our fiducial simulation, and find that blending-related effects are the dominant contribution to the mean multiplicative bias of approximately $-2\%$. By generating simulations with input shear signals that vary with redshift, we calibrate biases in our estimation of the effective redshfit distribution, and demonstrate the importance of this approach when blending is present. We provide corrected effective redshift distributions that incorporate statistical and systematic uncertainties, ready for use in DES Year 3 weak lensing analyses.

preprint2022arXiv

From the Fire: A Deeper Look at the Phoenix Stream

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.

preprint2022arXiv

Superclustering with the Atacama Cosmology Telescope and Dark Energy Survey: I. Evidence for thermal energy anisotropy using oriented stacking

The cosmic web contains filamentary structure on a wide range of scales. On the largest scales, superclustering aligns multiple galaxy clusters along inter-cluster bridges, visible through their thermal Sunyaev-Zel&#39;dovich signal in the Cosmic Microwave Background. We demonstrate a new, flexible method to analyze the hot gas signal from multi-scale extended structures. We use a Compton-$y$ map from the Atacama Cosmology Telescope (ACT) stacked on redMaPPer cluster positions from the optical Dark Energy Survey (DES). Cutout images from the $y$ map are oriented with large-scale structure information from DES galaxy data such that the superclustering signal is aligned before being overlaid. We find evidence for an extended quadrupole moment of the stacked $y$ signal at the 3.5$σ$ level, demonstrating that the large-scale thermal energy surrounding galaxy clusters is anisotropically distributed. We compare our ACT$\times$DES results with the Buzzard simulations, finding broad agreement. Using simulations, we highlight the promise of this novel technique for constraining the evolution of anisotropic, non-Gaussian structure using future combinations of microwave and optical surveys.

preprint2022arXiv

Synthetic Galaxy Clusters and Observations Based on Dark Energy Survey Year 3 Data

We develop a novel data-driven method for generating synthetic optical observations of galaxy clusters. In cluster weak lensing, the interplay between analysis choices and systematic effects related to source galaxy selection, shape measurement and photometric redshift estimation can be best characterized in end-to-end tests going from mock observations to recovered cluster masses. To create such test scenarios, we measure and model the photometric properties of galaxy clusters and their sky environments from the Dark Energy Survey Year 3 (DES Y3) data in two bins of cluster richness $λ\in[30;\,45)$, $λ\in[45;\,60)$ and three bins in cluster redshift ($z\in[0.3;\,0.35)$, $z\in[0.45;\,0.5)$ and $z\in[0.6;\,0.65)$. Using deep-field imaging data we extrapolate galaxy populations beyond the limiting magnitude of DES Y3 and calculate the properties of cluster member galaxies via statistical background subtraction. We construct mock galaxy clusters as random draws from a distribution function, and render mock clusters and line-of-sight catalogs into synthetic images in the same format as actual survey observations. Synthetic galaxy clusters are generated from real observational data, and thus are independent from the assumptions inherent to cosmological simulations. The recipe can be straightforwardly modified to incorporate extra information, and correct for survey incompleteness. New realizations of synthetic clusters can be created at minimal cost, which will allow future analyses to generate the large number of images needed to characterize systematic uncertainties in cluster mass measurements.

preprint2022arXiv

The Dark Energy Survey 5-year photometrically identified Type Ia Supernovae

As part of the cosmology analysis using Type Ia Supernovae (SN Ia) in the Dark Energy Survey (DES), we present photometrically identified SN Ia samples using multi-band light-curves and host galaxy redshifts. For this analysis, we use the photometric classification framework SuperNNova (SNN; Möller et al. 2019) trained on realistic DES-like simulations. For reliable classification, we process the DES SN programme (DES-SN) data and introduce improvements to the classifier architecture, obtaining classification accuracies of more than 98 per cent on simulations. This is the first SN classification to make use of ensemble methods, resulting in more robust samples. Using photometry, host galaxy redshifts, and a classification probability requirement, we identify 1,863 SNe Ia from which we select 1,484 cosmology-grade SNe Ia spanning the redshift range of 0.07 < z < 1.14. We find good agreement between the light-curve properties of the photometrically-selected sample and simulations. Additionally, we create similar SN Ia samples using two types of Bayesian Neural Network classifiers that provide uncertainties on the classification probabilities. We test the feasibility of using these uncertainties as indicators for out-of-distribution candidates and model confidence. Finally, we discuss the implications of photometric samples and classification methods for future surveys such as Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST).

preprint2022arXiv

The DECam Local Volume Exploration Survey Data Release 2

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

preprint2022arXiv

The DES Bright Arcs Survey: Candidate Strongly Lensed Galaxy Systems from the Dark Energy Survey 5,000 Sq. Deg. Footprint

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σ$) 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 &#34;tiles&#34; 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.

preprint2022arXiv

The Evolution of AGN Activity in Brightest Cluster Galaxies

We present the results of an analysis of Wide-field Infrared Survey Explorer (WISE) observations on the full 2500 deg^2 South Pole Telescope (SPT)-SZ cluster sample. We describe a process for identifying active galactic nuclei (AGN) in brightest cluster galaxies (BCGs) based on WISE mid-infrared color and redshift. Applying this technique to the BCGs of the SPT-SZ sample, we calculate the AGN-hosting BCG fraction, which is defined as the fraction of BCGs hosting bright central AGNs over all possible BCGs. Assuming {\bf an evolving} single-burst stellar population model, we find statistically significant evidence (>99.9%) for a mid-IR excess at high redshift compared to low redshift, suggesting that the fraction of AGN-hosting BCGs increases with redshift over the range of 0 < z < 1.3. The best-fit redshift trend of the AGN-hosting BCG fraction has the form (1+z)^(4.1+/-1.0). These results are consistent with previous studies in galaxy clusters as well as field galaxies. One way to explain this result is that member galaxies at high redshift tend to have more cold gas. While BCGs in nearby galaxy clusters grow mostly by dry mergers with cluster members, leading to no increase in AGN activity, BCGs at high redshift could primarily merge with gas-rich satellites, providing fuel for feeding AGNs. If this observed increase in AGN activity is linked to gas-rich mergers, rather than ICM cooling, we would expect to see an increase in scatter in the P_cav vs L_cool relation at z > 1. Lastly, this work confirms that the runaway cooling phase, as predicted by the classical cooling flow model, in the Phoenix cluster is extremely rare and most BCGs have low (relative to Eddington) black hole accretion rates.

preprint2021arXiv

A machine learning approach to galaxy properties: joint redshift-stellar mass probability distributions with Random Forest

We demonstrate that highly accurate joint redshift-stellar mass probability distribution functions (PDFs) can be obtained using the Random Forest (RF) machine learning (ML) algorithm, even with few photometric bands available. As an example, we use the Dark Energy Survey (DES), combined with the COSMOS2015 catalogue for redshifts and stellar masses. We build two ML models: one containing deep photometry in the $griz$ bands, and the second reflecting the photometric scatter present in the main DES survey, with carefully constructed representative training data in each case. We validate our joint PDFs for $10,699$ test galaxies by utilizing the copula probability integral transform and the Kendall distribution function, and their univariate counterparts to validate the marginals. Benchmarked against a basic set-up of the template-fitting code BAGPIPES, our ML-based method outperforms template fitting on all of our predefined performance metrics. In addition to accuracy, the RF is extremely fast, able to compute joint PDFs for a million galaxies in just under $6$ min with consumer computer hardware. Such speed enables PDFs to be derived in real time within analysis codes, solving potential storage issues. As part of this work we have developed GALPRO, a highly intuitive and efficient Python package to rapidly generate multivariate PDFs on-the-fly. GALPRO is documented and available for researchers to use in their cosmology and galaxy evolution studies.

preprint2021arXiv

A search of the full six years of the Dark Energy Survey for outer Solar System objects

We present the results of a search for outer Solar System objects in the full six years of data (Y6) from the Dark Energy Survey (DES). The DES covered a contiguous $5000$ deg$^2$ of the southern sky with $\approx 80,000$ $3$ deg$^2$ exposures in the $grizY$ optical/IR filters between 2013 and 2019. This search yielded 815 trans-Neptunian objects (TNOs), one Centaur and one Oort cloud comet, with 461 objects reported for the first time in this paper. We present methodology that builds upon our previous search carried out on the first four years of data. Here, all DES images were reprocessed with an improved detection pipeline that leads to an average completeness gain of 0.47 mag per exposure, as well as an improved transient catalog production and optimized algorithms for linkage of detections into orbits. All objects were verified by visual inspection and by computing the sub-threshold significance, the total signal-to-noise ratio in the stack of images in which the object&#39;s presence is indicated by the orbit fit, but no detection was reported. This yields a highly pure catalog of TNOs complete to $r \approx 23.8$ mag and distances $29<d<2500$ au. The Y6 TNOs have minimum (median) of 7 (12) distinct nights&#39; detections and arcs of 1.1 (4.2) years, and will have $grizY$ magnitudes available in a further publication. We present software for simulating our observational biases that enable comparisons of population models to our detections. Initial inferences demonstrating the statistical power of the DES catalog are: the data are inconsistent with the CFEPS-L7 model for the classical Kuiper Belt; the 16 ``extreme&#39;&#39; TNOs ($a>150$ au, $q>30$ au) are consistent with the null hypothesis of azimuthal isotropy; and non-resonant TNOs with $q>38$ au, $a>50$ au show a highly significant tendency to be sunward of the major mean motion resonances, whereas this tendency is not present for $q<38$ au.

preprint2021arXiv

Cross-correlation of DES Y3 lensing and ACT/${\it Planck}$ thermal Sunyaev Zel&#39;dovich Effect I: Measurements, systematics tests, and feedback model constraints

We present a tomographic measurement of the cross-correlation between thermal Sunyaev-Zeldovich (tSZ) maps from ${\it Planck}$ and the Atacama Cosmology Telescope (ACT) and weak galaxy lensing shears measured during the first three years of observations of the Dark Energy Survey (DES Y3). This correlation is sensitive to the thermal energy in baryons over a wide redshift range, and is therefore a powerful probe of astrophysical feedback. We detect the correlation at a statistical significance of $21σ$, the highest significance to date. We examine the tSZ maps for potential contaminants, including cosmic infrared background (CIB) and radio sources, finding that CIB has a substantial impact on our measurements and must be taken into account in our analysis. We use the cross-correlation measurements to test different feedback models. In particular, we model the tSZ using several different pressure profile models calibrated against hydrodynamical simulations. Our analysis marginalises over redshift uncertainties, shear calibration biases, and intrinsic alignment effects. We also marginalise over $Ω_{\rm m}$ and $σ_8$ using ${\it Planck}$ or DES priors. We find that the data prefers the model with a low amplitude of the pressure profile at small scales, compatible with a scenario with strong AGN feedback and ejection of gas from the inner part of the halos. When using a more flexible model for the shear profile, constraints are weaker, and the data cannot discriminate between different baryonic prescriptions.

preprint2021arXiv

Dark Energy Survey Year 3 Results: Optimizing the Lens Sample in Combined Galaxy Clustering and Galaxy-Galaxy Lensing Analysis

We investigate potential gains in cosmological constraints from the combination of galaxy clustering and galaxy-galaxy lensing by optimizing the lens galaxy sample selection using information from Dark Energy Survey (DES) Year 3 data and assuming the DES Year 1 Metacalibration sample for the sources. We explore easily reproducible selections based on magnitude cuts in $i$-band as a function of (photometric) redshift, $z_{\rm phot}$, and benchmark the potential gains against those using the well established redMaGiC sample. We focus on the balance between density and photometric redshift accuracy, while marginalizing over a realistic set of cosmological and systematic parameters. Our optimal selection, the MagLim sample, satisfies $i < 4 \, z_{\rm phot} + 18$ and has $\sim 30\%$ wider redshift distributions but $\sim 3.5$ times more galaxies than redMaGiC. Assuming a wCDM model and equivalent scale cuts to mitigate nonlinear effects, this leads to $40\%$ increase in the figure of merit for the pair combinations of $Ω_m$, $w$, and $σ_8$, and gains of $16\%$ in $σ_8$, $10\%$ in $Ω_m$, and $12\%$ in $w$. Similarly, in LCDM we find an improvement of $19\%$ and $27\%$ on $σ_8$ and $Ω_m$, respectively. We also explore flux-limited samples with a flat magnitude cut finding that the optimal selection, $i < 22.2$, has $\sim 7$ times more galaxies and $\sim 20\%$ wider redshift distributions compared to MagLim, but slightly worse constraints. We show that our results are robust with respect to the assumed galaxy bias and photometric redshift uncertainties with only moderate further gains from increased number of tomographic bins or the inclusion of bin cross-correlations, except in the case of the flux-limited sample, for which these gains are more significant.

preprint2021arXiv

Finding quadruply imaged quasars with machine learning. I. Methods

Strongly lensed quadruply imaged quasars (quads) are extraordinary objects. They are very rare in the sky -- only a few tens are known to date -- and yet they provide unique information about a wide range of topics, including the expansion history and the composition of the Universe, the distribution of stars and dark matter in galaxies, the host galaxies of quasars, and the stellar initial mass function. Finding them in astronomical images is a classic &#34;needle in a haystack&#34; problem, as they are outnumbered by other (contaminant) sources by many orders of magnitude. To solve this problem, we develop state-of-the-art deep learning methods and train them on realistic simulated quads based on real images of galaxies taken from the Dark Energy Survey, with realistic source and deflector models, including the chromatic effects of microlensing. The performance of the best methods on a mixture of simulated and real objects is excellent, yielding area under the receiver operating curve in the range 0.86 to 0.89. Recall is close to 100% down to total magnitude i~21 indicating high completeness, while precision declines from 85% to 70% in the range i~17-21. The methods are extremely fast: training on 2 million samples takes 20 hours on a GPU machine, and 10^8 multi-band cutouts can be evaluated per GPU-hour. The speed and performance of the method pave the way to apply it to large samples of astronomical sources, bypassing the need for photometric pre-selection that is likely to be a major cause of incompleteness in current samples of known quads.

preprint2021arXiv

Is diffuse intracluster light a good tracer of the galaxy cluster matter distribution?

We explore the relation between diffuse intracluster light (central galaxy included) and the galaxy cluster (baryonic and dark) matter distribution using a sample of 528 clusters at $0.2\leq z \leq 0.35$ found in the Dark Energy Survey (DES) Year 1 data. The surface brightness of the diffuse light shows an increasing dependence on cluster total mass at larger radius, and appears to be self-similar with a universal radial dependence after scaling by cluster radius. We also compare the diffuse light radial profiles to the cluster (baryonic and dark) matter distribution measured through weak lensing and find them to be comparable. The IllustrisTNG galaxy formation simulation, TNG300, offers further insight into the connection between diffuse stellar mass and cluster matter distributions -- the simulation radial profile of the diffuse stellar component does not have a similar slope with the total cluster matter content, although that of the cluster satellite galaxies does. Regardless of the radial trends, the amount of diffuse stellar mass has a low-scatter scaling relation with cluster&#39;s total mass in the simulation, out-performing the total stellar mass of cluster satellite galaxies. We conclude that there is no consistent evidence yet on whether or not diffuse light is a faithful radial tracer of the cluster matter distribution. Nevertheless, both observational and simulation results reveal that diffuse light is an excellent indicator of the cluster&#39;s total mass.

preprint2021arXiv

Milky Way Satellite Census. III. Constraints on Dark Matter Properties from Observations of Milky Way Satellite Galaxies

We perform a comprehensive study of Milky Way (MW) satellite galaxies to constrain the fundamental properties of dark matter (DM). This analysis fully incorporates inhomogeneities in the spatial distribution and detectability of MW satellites and marginalizes over uncertainties in the mapping between galaxies and DM halos, the properties of the MW system, and the disruption of subhalos by the MW disk. Our results are consistent with the cold, collisionless DM paradigm and yield the strongest cosmological constraints to date on particle models of warm, interacting, and fuzzy dark matter. At $95\%$ confidence, we report limits on (i) the mass of thermal relic warm DM, $m_{\rm WDM} > 6.5\ \mathrm{keV}$ (free-streaming length, $λ_{\rm{fs}} \lesssim 10\,h^{-1}\ \mathrm{kpc}$), (ii) the velocity-independent DM-proton scattering cross section, $σ_{0} < 8.8\times 10^{-29}\ \mathrm{cm}^{2}$ for a $100\ \mathrm{MeV}$ DM particle mass (DM-proton coupling, $c_p \lesssim (0.3\ \mathrm{GeV})^{-2}$), and (iii) the mass of fuzzy DM, $m_ϕ> 2.9 \times 10^{-21}\ \mathrm{eV}$ (de Broglie wavelength, $λ_{\rm{dB}} \lesssim 0.5\ \mathrm{kpc}$). These constraints are complementary to other observational and laboratory constraints on DM properties.

preprint2021arXiv

No Evidence for Orbital Clustering in the Extreme Trans-Neptunian Objects

The apparent clustering in longitude of perihelion $\varpi$ and ascending node $Ω$ of extreme trans-Neptunian objects (ETNOs) has been attributed to the gravitational effects of an unseen 5-10 Earth-mass planet in the outer solar system. To investigate how selection bias may contribute to this clustering, we consider 14 ETNOs discovered by the Dark Energy Survey, the Outer Solar System Origins Survey, and the survey of Sheppard and Trujillo. Using each survey&#39;s published pointing history, depth, and TNO tracking selections, we calculate the joint probability that these objects are consistent with an underlying parent population with uniform distributions in $\varpi$ and $Ω$. We find that the mean scaled longitude of perihelion and orbital poles of the detected ETNOs are consistent with a uniform population at a level between $17\%$ and $94\%$, and thus conclude that this sample provides no evidence for angular clustering.

preprint2021arXiv

The Dark Energy Survey Supernova Program: Cosmological biases from supernova photometric classification

Cosmological analyses of samples of photometrically-identified Type Ia supernovae (SNe Ia) depend on understanding the effects of &#39;contamination&#39; from core-collapse and peculiar SN Ia events. We employ a rigorous analysis on state-of-the-art simulations of photometrically identified SN Ia samples and determine cosmological biases due to such &#39;non-Ia&#39; contamination in the Dark Energy Survey (DES) 5-year SN sample. As part of the analysis, we test on our DES simulations the performance of SuperNNova, a photometric SN classifier based on recurrent neural networks. Depending on the choice of non-Ia SN models in both the simulated data sample and training sample, contamination ranges from 0.8-3.5 %, with the efficiency of the classification from 97.7-99.5 %. Using the Bayesian Estimation Applied to Multiple Species (BEAMS) framework and its extension &#39;BEAMS with Bias Correction&#39; (BBC), we produce a redshift-binned Hubble diagram marginalised over contamination and corrected for selection effects and we use it to constrain the dark energy equation-of-state, $w$. Assuming a flat universe with Gaussian $Ω_M$ prior of $0.311\pm0.010$, we show that biases on $w$ are $<0.008$ when using SuperNNova and accounting for a wide range of non-Ia SN models in the simulations. Systematic uncertainties associated with contamination are estimated to be at most $σ_{w, \mathrm{syst}}=0.004$. This compares to an expected statistical uncertainty of $σ_{w,\mathrm{stat}}=0.039$ for the DES-SN sample, thus showing that contamination is not a limiting uncertainty in our analysis. We also measure biases due to contamination on $w_0$ and $w_a$ (assuming a flat universe), and find these to be $<$0.009 in $w_0$ and $<$0.108 in $w_a$, hence 5 to 10 times smaller than the statistical uncertainties expected from the DES-SN sample.

preprint2021arXiv

The Effect of Environment on Type Ia Supernovae in the Dark Energy Survey Three-Year Cosmological Sample

Analyses of type Ia supernovae (SNe Ia) have found puzzling correlations between their standardised luminosities and host galaxy properties: SNe Ia in high-mass, passive hosts appear brighter than those in lower-mass, star-forming hosts. We examine the host galaxies of SNe Ia in the Dark Energy Survey three-year spectroscopically-confirmed cosmological sample, obtaining photometry in a series of &#34;local&#34; apertures centred on the SN, and for the global host galaxy. We study the differences in these host galaxy properties, such as stellar mass and rest-frame $U-R$ colours, and their correlations with SN Ia parameters including Hubble residuals. We find all Hubble residual steps to be $>3σ$ in significance, both for splitting at the traditional environmental property sample median and for the step of maximum significance. For stellar mass, we find a maximal local step of $0.098\pm0.018$ mag; $\sim 0.03$ mag greater than the largest global stellar mass step in our sample ($0.070 \pm 0.017$ mag). When splitting at the sample median, differences between local and global $U-R$ steps are small, both $\sim 0.08$ mag, but are more significant than the global stellar mass step ($0.057\pm0.017$ mag). We split the data into sub-samples based on SN Ia light curve parameters: stretch ($x_1$) and colour ($c$), finding that redder objects ($c > 0$) have larger Hubble residual steps, for both stellar mass and $U-R$, for both local and global measurements, of $\sim0.14$ mag. Additionally, the bluer (star-forming) local environments host a more homogeneous SN Ia sample, with local $U-R$ r.m.s. scatter as low as $0.084 \pm 0.017$ mag for blue ($c < 0$) SNe Ia in locally blue $U-R$ environments.

preprint2021arXiv

The Observed Evolution of the Stellar Mass - Halo Mass Relation for Brightest Central Galaxies

We quantify evolution in the cluster scale stellar mass - halo mass (SMHM) relation&#39;s parameters using 2323 clusters and brightest central galaxies (BCGs) over the redshift range $0.03 \le z \le 0.60$. The precision on inferred SMHM parameters is improved by including the magnitude gap ($\rm m_{gap}$) between the BCG and fourth brightest cluster member (M14) as a third parameter in the SMHM relation. At fixed halo mass, accounting for $\rm m_{gap}$, through a stretch parameter, reduces the SMHM relation&#39;s intrinsic scatter. To explore this redshift range, we use clusters, BCGs, and cluster members identified using the Sloan Digital Sky Survey C4 and redMaPPer cluster catalogs and the Dark Energy Survey redMaPPer catalog. Through this joint analysis, we detect no systematic differences in BCG stellar mass, $\rm m_{gap}$, and cluster mass (inferred from richness) between the datsets. We utilize the Pareto function to quantify each parameter&#39;s evolution. We confirm prior findings of negative evolution in the SMHM relation&#39;s slope (3.5$σ$) and detect negative evolution in the stretch parameter (4.0$σ$) and positive evolution in the offset parameter (5.8$σ$). This observed evolution, combined with the absence of BCG growth, when stellar mass is measured within 50kpc, suggests that this evolution results from changes in the cluster&#39;s $\rm m_{gap}$. For this to occur, late-term growth must be in the intra-cluster light surrounding the BCG. We also compare the observed results to Illustris TNG 300-1 cosmological hydrodynamic simulations and find modest qualitative agreement. However, the simulations lack the evolutionary features detected in the real data.

preprint2020arXiv

A joint SZ-Xray-optical analysis of the dynamical state of 288 massive galaxy clusters

We use imaging from the first three years of the Dark Energy Survey to characterize the dynamical state of 288 galaxy clusters at $0.1 \lesssim z \lesssim 0.9$ detected in the South Pole Telescope (SPT) Sunyaev-Zeldovich (SZ) effect survey (SPT-SZ). We examine spatial offsets between the position of the brightest cluster galaxy (BCG) and the center of the gas distribution as traced by the SPT-SZ centroid and by the X-ray centroid/peak position from Chandra and XMM data. We show that the radial distribution of offsets provides no evidence that SPT SZ-selected cluster samples include a higher fraction of mergers than X-ray-selected cluster samples. We use the offsets to classify the dynamical state of the clusters, selecting the 43 most disturbed clusters, with half of those at $z \gtrsim 0.5$, a region seldom explored previously. We find that Schechter function fits to the galaxy population in disturbed clusters and relaxed clusters differ at $z>0.55$ but not at lower redshifts. Disturbed clusters at $z>0.55$ have steeper faint-end slopes and brighter characteristic magnitudes. Within the same redshift range, we find that the BCGs in relaxed clusters tend to be brighter than the BCGs in disturbed samples, while in agreement in the lower redshift bin. Possible explanations includes a higher merger rate, and a more efficient dynamical friction at high redshift. The red-sequence population is less affected by the cluster dynamical state than the general galaxy population.

preprint2020arXiv

Birds of a Feather? Magellan/IMACS Spectroscopy of the Ultra-Faint Satellites Grus II, Tucana IV, and Tucana V

We present Magellan/IMACS spectroscopy of three recently discovered ultra-faint Milky Way satellites, Grus II, Tucana IV, and Tucana V. We measure systemic velocities of V_hel = -110.0 +/- 0.5 km/s, V_hel = 15.9 +/- 1.8 km/s, and V_hel = -36.2 +/-2.5 km/s for the three objects, respectively. Their large relative velocity differences demonstrate that the satellites are unrelated despite their close physical proximity. We determine a velocity dispersion for Tuc IV of sigma = 4.3^+1.7_-1.0 km/s, but we cannot resolve the velocity dispersions of the other two systems. For Gru II we place an upper limit (90% confidence) on the dispersion of sigma < 1.9 km/s, and for Tuc V we do not obtain any useful limits. All three satellites have metallicities below [Fe/H] = -2.1, but none has a detectable metallicity spread. We determine proper motions for each satellite based on Gaia astrometry and compute their orbits around the Milky Way. Gru II is on a tightly bound orbit with a pericenter of 25 kpc and orbital eccentricity of 0.45. Tuc V likely has an apocenter beyond 100 kpc, and could be approaching the Milky Way for the first time. The current orbit of Tuc IV is similar to that of Gru II, with a pericenter of 25 kpc and an eccentricity of 0.36. However, a backward integration of the position of Tuc IV demonstrates that it collided with the Large Magellanic Cloud at an impact parameter of 4 kpc ~120 Myr ago, deflecting its trajectory and possibly altering its internal kinematics. Based on their sizes, masses, and metallicities, we classify Gru II and Tuc IV as likely dwarf galaxies, but the nature of Tuc V remains uncertain.

preprint2020arXiv

Blinding multi-probe cosmological experiments

The goal of blinding is to hide an experiment&#39;s critical results -- here the inferred cosmological parameters -- until all decisions affecting its analysis have been finalised. This is especially important in the current era of precision cosmology, when the results of any new experiment are closely scrutinised for consistency or tension with previous results. In analyses that combine multiple observational probes, like the combination of galaxy clustering and weak lensing in the Dark Energy Survey (DES), it is challenging to blind the results while retaining the ability to check for (in)consistency between different parts of the data. We propose a simple new blinding transformation that works by modifying the summary statistics that are input to parameter estimation, such as two-point correlation functions. The transformation shifts the measured statistics to new values that are consistent with (blindly) shifted cosmological parameters, while preserving internal (in)consistency. We apply the blinding transformation to simulated data for the projected DES Year 3 galaxy clustering and weak lensing analysis, demonstrating that practical blinding is achieved without significant perturbation of internal-consistency checks, as measured here by degradation of the $χ^2$ between data and best-fitting model. Our blinding method conserves $χ^2$ more precisely as experiments evolve to higher precision.

preprint2020arXiv

Chemical Analysis of the Ultra-Faint Dwarf Galaxy Grus~II. Signature of high-mass stellar nucleosynthesis

We present a detailed abundance analysis of the three brightest member stars at the top of the giant branch of the ultra-faint dwarf galaxy Grus~II. All stars exhibit a higher than expected $\mathrm{[Mg/Ca]}$ ratio compared to metal-poor stars in other ultra-faint dwarf galaxies and in the Milky Way halo. Nucleosynthesis in high mass ($\geqslant 20$M$_\odot$) core-collapse supernovae has been shown to create this signature. The abundances of this small sample (3) stars suggest the chemical enrichment of Grus~II could have occurred through substantial high-mass stellar evolution and is consistent with the framework of a top-heavy initial mass function. However, with only three stars it can not be ruled out that the abundance pattern is the result of a stochastic chemical enrichment at early times in the galaxy. The most metal-rich of the three stars also possesses a small enhancement in rapid neutron-capture ($r$-process) elements. The abundance pattern of the $r$-process elements in this star matches the scaled $r$-process pattern of the solar system and $r$-process enhanced stars in other dwarf galaxies and in the Milky Way halo, hinting at a common origin for these elements across a range of environments. All current proposed astrophysical sites of $r$-process element production are associated with high-mass stars, thus the possible top-heavy initial mass function of Grus~II would increase the likelihood of any of these events occurring. The time delay between the $α$ and $r$-process element enrichment of the galaxy favors a neutron star merger as the origin of the $r$-process elements in Grus~II.

preprint2020arXiv

DES16C3cje: A low-luminosity, long-lived supernova

We present DES16C3cje, a low-luminosity, long-lived type II supernova (SN II) at redshift 0.0618, detected by the Dark Energy Survey (DES). DES16C3cje is a unique SN. The spectra are characterized by extremely narrow photospheric lines corresponding to very low expansion velocities of $\lesssim1500$ km s$^{-1}$, and the light curve shows an initial peak that fades after 50 days before slowly rebrightening over a further 100 days to reach an absolute brightness of M$_r\sim -15.5$ mag. The decline rate of the late-time light curve is then slower than that expected from the powering by radioactive decay of $^{56}$Co but is comparable to that expected from accretion power. Comparing the bolometric light curve with hydrodynamical models, we find that DES16C3cje can be explained by either i) a low explosion energy (0.11 foe) and relatively large $^{56}$Ni production of 0.075 M$_{\odot}$ from a $\sim15$ M$_{\odot}$ red supergiant progenitor typical of other SNe II, or ii) a relatively compact $\sim40$ M$_{\odot}$ star, explosion energy of 1 foe, and 0.08 M$_{\odot}$ of $^{56}$Ni. Both scenarios require additional energy input to explain the late-time light curve, which is consistent with fallback accretion at a rate of $\sim0.5\times{10^{-8}}$ M$_{\odot}$ s$^{-1}$.

preprint2020arXiv

Detection of cross-correlation between gravitational lensing and gamma rays

In recent years, many gamma-ray sources have been identified, yet the unresolved component hosts valuable information on the faintest emission. In order to extract it, a cross-correlation with gravitational tracers of matter in the Universe has been shown to be a promising tool. We report here the first identification of a cross-correlation signal between gamma rays and the distribution of mass in the Universe probed by weak gravitational lensing. We use the Dark Energy Survey Y1 weak lensing catalogue and the Fermi Large Area Telescope 9-year gamma-ray data, obtaining a signal-to-noise ratio of 5.3. The signal is mostly localised at small angular scales and high gamma-ray energies, with a hint of correlation at extended separation. Blazar emission is likely the origin of the small-scale effect. We investigate implications of the large-scale component in terms of astrophysical sources and particle dark matter emission.

preprint2020arXiv

Dynamical Classification of Trans-Neptunian Objects Detected by the Dark Energy Survey

The outer Solar System contains a large number of small bodies (known as trans-Neptunian objects or TNOs) that exhibit diverse types of dynamical behavior. The classification of bodies in this distant region into dynamical classes -- sub-populations that experience similar orbital evolution -- aids in our understanding of the structure and formation of the Solar System. In this work, we propose an updated dynamical classification scheme for the outer Solar System. This approach includes the construction of a new (automated) method for identifying mean-motion resonances. We apply this algorithm to the current dataset of TNOs observed by the Dark Energy Survey (DES) and present a working classification for all of the DES TNOs detected to date. Our classification scheme yields 1 inner centaur, 19 outer centaurs, 21 scattering disk objects, 47 detached TNOs, 48 securely resonant objects, 7 resonant candidates, and 97 classical belt objects. Among the scattering and detached objects, we detect 8 TNOs with semi-major axes greater than 150 AU.

preprint2020arXiv

Is every strong lens model unhappy in its own way? Uniform modelling of a sample of 13 quadruply+ imaged quasars

Strong-gravitational lens systems with quadruply-imaged quasars (quads) are unique probes to address several fundamental problems in cosmology and astrophysics. Although they are intrinsically very rare, ongoing and planned wide-field deep-sky surveys are set to discover thousands of such systems in the next decade. It is thus paramount to devise a general framework to model strong-lens systems to cope with this large influx without being limited by expert investigator time. We propose such a general modelling framework (implemented with the publicly available software Lenstronomy) and apply it to uniformly model three-band Hubble Space Telescope Wide Field Camera 3 images of 13 quads. This is the largest uniformly modelled sample of quads to date and paves the way for a variety of studies. To illustrate the scientific content of the sample, we investigate the alignment between the mass and light distribution in the deflectors. The position angles of these distributions are well-aligned, except when there is strong external shear. However, we find no correlation between the ellipticity of the light and mass distributions. We also show that the observed flux-ratios between the images depart significantly from the predictions of simple smooth models. The departures are strongest in the bluest band, consistent with microlensing being the dominant cause in addition to millilensing. Future papers will exploit this rich dataset in combination with ground based spectroscopy and time delays to determine quantities such as the Hubble constant, the free streaming length of dark matter, and the normalization of the initial stellar mass function.

preprint2020arXiv

Milky Way Satellite Census. II. Galaxy--Halo Connection Constraints Including the Impact of the Large Magellanic Cloud

The population of Milky Way (MW) satellites contains the faintest known galaxies and thus provides essential insight into galaxy formation and dark matter microphysics. Here we combine a model of the galaxy--halo connection with newly derived observational selection functions based on searches for satellites in photometric surveys over nearly the entire high Galactic latitude sky. In particular, we use cosmological zoom-in simulations of MW-like halos that include realistic Large Magellanic Cloud (LMC) analogs to fit the position-dependent MW satellite luminosity function. We report decisive evidence for the statistical impact of the LMC on the MW satellite population due to an estimated $6\pm 2$ observed LMC-associated satellites, consistent with the number of LMC satellites inferred from Gaia proper-motion measurements, confirming the predictions of cold dark matter models for the existence of satellites within satellite halos. Moreover, we infer that the LMC fell into the MW within the last $2\ \rm{Gyr}$ at high confidence. Based on our detailed full-sky modeling, we find that the faintest observed satellites inhabit halos with peak virial masses below $3.2\times 10^{8}\ M_{\rm{\odot}}$ at $95\%$ confidence, and we place the first robust constraints on the fraction of halos that host galaxies in this regime. We predict that the faintest detectable satellites occupy halos with peak virial masses above $10^{6}\ M_{\rm{\odot}}$, highlighting the potential for powerful galaxy formation and dark matter constraints from future dwarf galaxy searches.

preprint2020arXiv

Modelling the Milky Way. I -- Method and first results fitting the thick disk and halo with DES-Y3 data

We present MWFitting, a method to fit the stellar components of the Galaxy by comparing Hess Diagrams (HDs) from TRILEGAL models to real data. We apply MWFitting to photometric data from the first three years of the Dark Energy Survey (DES). After removing regions containing known resolved stellar systems such as globular clusters, dwarf galaxies, nearby galaxies, the Large Magellanic Cloud and the Sagittarius Stream, our main sample spans a total area of $\sim$2,300 deg$^2$ distributed across the DES footprint. We further explore a smaller subset ($\sim$ 1,300 deg$^2$) that excludes all regions with known stellar streams and stellar overdensities. Validation tests on synthetic data possessing similar properties to the DES data show that the method is able to recover input parameters with a precision better than 3\%. Based on the best-fit models, we create simulated stellar catalogues covering the whole DES footprint down to $g = 24$ magnitude. Comparisons of data and simulations provide evidence for a break in the power law index describing the stellar density of the Milky Way (MW) halo. Several previously discovered stellar over-densities are recovered in the residual stellar density map, showing the reliability of MWFitting in determining the Galactic components. Simulations made with the best-fitting parameters are a promising way to predict MW star counts for surveys such as LSST and Euclid.

preprint2020arXiv

STRIDES: a 3.9 per cent measurement of the Hubble constant from the strong lens system DES J0408-5354

We present a blind time-delay cosmographic analysis for the lens system DES J0408$-$5354. This system is extraordinary for the presence of two sets of multiple images at different redshifts, which provide the opportunity to obtain more information at the cost of increased modelling complexity with respect to previously analyzed systems. We perform detailed modelling of the mass distribution for this lens system using three band Hubble Space Telescope imaging. We combine the measured time delays, line-of-sight central velocity dispersion of the deflector, and statistically constrained external convergence with our lens models to estimate two cosmological distances. We measure the &#34;effective&#34; time-delay distance corresponding to the redshifts of the deflector and the lensed quasar $D_{Δt}^{\rm eff}=3382^{+146}_{-115}$ Mpc and the angular diameter distance to the deflector $D_{\rm d}=1711^{+376}_{-280}$ Mpc, with covariance between the two distances. From these constraints on the cosmological distances, we infer the Hubble constant $H_0 = 74.2^{+2.7}_{-3.0}$ km s$^{-1}$ Mpc$^{-1}$ assuming a flat $Λ$CDM cosmology and a uniform prior for $Ω_{\rm m}$ as $Ω_{\rm m} \sim \mathcal{U}(0.05, 0.5)$. This measurement gives the most precise constraint on $H_0$ to date from a single lens. Our measurement is consistent with that obtained from the previous sample of six lenses analyzed by the $H_0$ Lenses in COSMOGRAIL&#39;s Wellspring (H0LiCOW) collaboration. It is also consistent with measurements of $H_0$ based on the local distance ladder, reinforcing the tension with the inference from early Universe probes, for example, with 2.2$σ$ discrepancy from the cosmic microwave background measurement.

preprint2020arXiv

Supernova Host Galaxies in the Dark Energy Survey: I. Deep Coadds, Photometry, and Stellar Masses

The five-year Dark Energy Survey supernova programme (DES-SN) is one of the largest and deepest transient surveys to date in terms of volume and number of supernovae. Identifying and characterising the host galaxies of transients plays a key role in their classification, the study of their formation mechanisms, and the cosmological analyses. To derive accurate host galaxy properties, we create depth-optimised coadds using single-epoch DES-SN images that are selected based on sky and atmospheric conditions. For each of the five DES-SN seasons, a separate coadd is made from the other 4 seasons such that each SN has a corresponding deep coadd with no contaminating SN emission. The coadds reach limiting magnitudes of order $\sim 27$ in $g$-band, and have a much smaller magnitude uncertainty than the previous DES-SN host templates, particularly for faint objects. We present the resulting multi-band photometry of host galaxies for samples of spectroscopically confirmed type Ia (SNe Ia), core-collapse (CCSNe), and superluminous (SLSNe) as well as rapidly evolving transients (RETs) discovered by DES-SN. We derive host galaxy stellar masses and probabilistically compare stellar-mass distributions to samples from other surveys. We find that the DES spectroscopically confirmed sample of SNe Ia selects preferentially fewer high mass hosts at high redshift compared to other surveys, while at low redshift the distributions are consistent. DES CCSNe and SLSNe hosts are similar to other samples, while RET hosts are unlike the hosts of any other transients, although these differences have not been disentangled from selection effects.

preprint2020arXiv

Supernova Siblings: Assessing the Consistency of Properties of Type Ia Supernovae that Share the Same Parent Galaxies

While many studies have shown a correlation between properties of the light curves of Type Ia SN (SNe Ia) and properties of their host galaxies, it remains unclear what is driving these correlations. We introduce a new direct method to study these correlations by analyzing `parent&#39; galaxies that host multiple SNe Ia &#39;siblings&#39;. Here, we search the Dark Energy Survey SN sample, one of the largest samples of discovered SNe, and find 8 galaxies that hosted two likely Type Ia SNe. Comparing the light-curve properties of these SNe and recovered distances from the light curves, we find no better agreement between properties of SNe in the same galaxy as any random pair of galaxies, with the exception of the SN light-curve stretch. We show at $2.8σ$ significance that at least 1/2 of the intrinsic scatter of SNe Ia distance modulus residuals is not from common host properties. We also discuss the robustness with which we could make this evaluation with LSST, which will find $100\times$ more pairs of galaxies, and pave a new line of study on the consistency of Type Ia supernovae in the same parent galaxies. Finally, we argue that it is unlikely some of these SNe are actually single, lensed SN with multiple images.

preprint2020arXiv

The impact of spectroscopic incompleteness in direct calibration of redshift distributions for weak lensing surveys

Obtaining accurate distributions of galaxy redshifts is a critical aspect of weak lensing cosmology experiments. One of the methods used to estimate and validate redshift distributions is apply weights to a spectroscopic sample so that their weighted photometry distribution matches the target sample. In this work we estimate the \textit{selection bias} in redshift that is introduced in this procedure. We do so by simulating the process of assembling a spectroscopic sample (including observer-assigned confidence flags) and highlight the impacts of spectroscopic target selection and redshift failures. We use the first year (Y1) weak lensing analysis in DES as an example data set but the implications generalise to all similar weak lensing surveys. We find that using colour cuts that are not available to the weak lensing galaxies can introduce biases of $Δ~z\sim0.015$ in the weighted mean redshift of different redshift intervals. To assess the impact of incompleteness in spectroscopic samples, we select only objects with high observer-defined confidence flags and compare the weighted mean redshift with the true mean. We find that the mean redshift of the DES Y1 weak lensing sample is typically biased at the $Δ~z=0.005-0.05$ level after the weighting is applied. The bias we uncover can have either sign, depending on the samples and redshift interval considered. For the highest redshift bin, the bias is larger than the uncertainties in the other DES Y1 redshift calibration methods, justifying the decision of not using this method for the redshift estimations. We discuss several methods to mitigate this bias.

preprint2020arXiv

The Mystery of Photometric Twins DES17X1boj and DES16E2bjy

We present an analysis of DES17X1boj and DES16E2bjy, two peculiar transients discovered by the Dark Energy Survey (DES). They exhibit nearly identical double-peaked light curves which reach very different maximum luminosities (M$_\mathrm{r}$ = -15.4 and M$_\mathrm{r}$ = -17.9, respectively). The light curve evolution of these events is highly atypical and has not been reported before. The transients are found in different host environments: DES17X1boj was found near the nucleus of a spiral galaxy, while DES16E2bjy is located in the outskirts of a passive red galaxy. Early photometric data is well fitted with a blackbody and the resulting moderate photospheric expansion velocities (1800 km/s for DES17X1boj and 4800 km/s for DES16E2bjy) suggest an explosive or eruptive origin. Additionally, a feature identified as high-velocity CaII absorption (v $\approx$ 9400km/s) in the near-peak spectrum of DES17X1boj may imply that it is a supernova. While similar light curve evolution suggests a similar physical origin for these two transients, we are not able to identify or characterise the progenitors.

preprint2020arXiv

Validation of Selection Function, Sample Contamination and Mass Calibration in Galaxy Cluster Samples

We construct and validate the selection function of the MARD-Y3 sample. This sample was selected through optical follow-up of the 2nd ROSAT faint source catalog (2RXS) with Dark Energy Survey year 3 (DES-Y3) data. The selection function is modeled by combining an empirically constructed X-ray selection function with an incompleteness model for the optical follow-up. We validate the joint selection function by testing the consistency of the constraints on the X-ray flux--mass and richness--mass scaling relation parameters derived from different sources of mass information: (1) cross-calibration using SPT-SZ clusters, (2) calibration using number counts in X-ray, in optical and in both X-ray and optical while marginalizing over cosmological parameters, and (3) other published analyses. We find that the constraints on the scaling relation from the number counts and SPT-SZ cross-calibration agree, indicating that our modeling of the selection function is adequate. Furthermore, we apply a largely cosmology independent method to validate selection functions via the computation of the probability of finding each cluster in the SPT-SZ sample in the MARD-Y3 sample and vice-versa. This test reveals no clear evidence for MARD-Y3 contamination, SPT-SZ incompleteness or outlier fraction. Finally, we discuss the prospects of the techniques presented here to limit systematic selection effects in future cluster cosmological studies.

preprint2020arXiv

Weak Lensing of Type Ia Supernovae from the Dark Energy Survey

We consider the effects of weak gravitational lensing on observations of 196 spectroscopically confirmed Type Ia Supernovae (SNe Ia) from years 1 to 3 of the Dark Energy Survey (DES). We simultaneously measure both the angular correlation function and the non-Gaussian skewness caused by weak lensing. This approach has the advantage of being insensitive to the intrinsic dispersion of SNe Ia magnitudes. We model the amplitude of both effects as a function of $σ_8$, and find $σ_8 = 1.2^{+0.9}_{-0.8}$. We also apply our method to a subsample of 488 SNe from the Joint Light-curve Analysis (JLA) (chosen to match the redshift range we use for this work), and find $σ_8 = 0.8^{+1.1}_{-0.7}$. The comparable uncertainty in $σ_8$ between DES-SN and the larger number of SNe from JLA highlights the benefits of homogeneity of the DES-SN sample, and improvements in the calibration and data analysis.

preprint2019arXiv

Dark Energy Survey Year 1 Results: Wide field mass maps via forward fitting in harmonic space

We present new wide-field weak lensing mass maps for the Year 1 Dark Energy Survey data, generated via a forward fitting approach. This method of producing maps does not impose any prior constraints on the mass distribution to be reconstructed. The technique is found to improve the map reconstruction on the edges of the field compared to the conventional Kaiser-Squires method, which applies a direct inversion on the data; our approach is in good agreement with the previous direct approach in the central regions of the footprint. The mapping technique is assessed and verified with tests on simulations; together with the Kaiser-Squires method, the technique is then applied to data from the Dark Energy Survey Year 1 data and the differences between the two methods are compared. We also produce the first DES measurements of the convergence Minkowski functionals and compare them to those measured in simulations.

preprint2019arXiv

Stellar mass as a galaxy cluster mass proxy: application to the Dark Energy Survey redMaPPer clusters

We introduce a galaxy cluster mass observable, $μ_\star$, based on the stellar masses of cluster members, and we present results for the Dark Energy Survey (DES) Year 1 observations. Stellar masses are computed using a Bayesian Model Averaging method, and are validated for DES data using simulations and COSMOS data. We show that $μ_\star$ works as a promising mass proxy by comparing our predictions to X-ray measurements. We measure the X-ray temperature-$μ_\star$ relation for a total of 150 clusters matched between the wide-field DES Year 1 redMaPPer catalogue, and Chandra and XMM archival observations, spanning the redshift range $0.1<z<0.7$. For a scaling relation which is linear in logarithmic space, we find a slope of $α= 0.488\pm0.043$ and a scatter in the X-ray temperature at fixed $μ_\star$ of $σ_{{\rm ln} T_X|μ_\star}=0.266^{+0.019}_{-0.020}$ for the joint sample. By using the halo mass scaling relations of the X-ray temperature from the Weighing the Giants program, we further derive the $μ_\star$-conditioned scatter in mass, finding $σ_{{\rm ln} M|μ_\star}=0.26^{+ 0.15}_{- 0.10}$. These results are competitive with well-established cluster mass proxies used for cosmological analyses, showing that $μ_\star$ can be used as a reliable and physically motivated mass proxy to derive cosmological constraints.

preprint2019arXiv

The SPTpol Extended Cluster Survey

We describe the observations and resultant galaxy cluster catalog from the 2770 deg$^2$ SPTpol Extended Cluster Survey (SPT-ECS). Clusters are identified via the Sunyaev-Zel&#39;dovich (SZ) effect, and confirmed with a combination of archival and targeted follow-up data, making particular use of data from the Dark Energy Survey (DES). With incomplete followup we have confirmed as clusters 244 of 266 candidates at a detection significance $ξ\ge 5$ and an additional 204 systems at $4<ξ<5$. The confirmed sample has a median mass of $M_{500c} \sim {4.4 \times 10^{14} M_\odot h_{70}^{-1}}$, a median redshift of $z=0.49$, and we have identified 44 strong gravitational lenses in the sample thus far. Radio data are used to characterize contamination to the SZ signal; the median contamination for confirmed clusters is predicted to be $\sim$1% of the SZ signal at the $ξ>4$ threshold, and $<4\%$ of clusters have a predicted contamination $>10\% $ of their measured SZ flux. We associate SZ-selected clusters, from both SPT-ECS and the SPT-SZ survey, with clusters from the DES redMaPPer sample, and find an offset distribution between the SZ center and central galaxy in general agreement with previous work, though with a larger fraction of clusters with significant offsets. Adopting a fixed Planck-like cosmology, we measure the optical richness-to-SZ-mass ($λ-M$) relation and find it to be 28% shallower than that from a weak-lensing analysis of the DES data---a difference significant at the 4 $σ$ level---with the relations intersecting at $λ=60$ . The SPT-ECS cluster sample will be particularly useful for studying the evolution of massive clusters and, in combination with DES lensing observations and the SPT-SZ cluster sample, will be an important component of future cosmological analyses.

preprint2019arXiv

The STRong lensing Insights into the Dark Energy Survey (STRIDES) 2017/2018 follow-up campaign: Discovery of 10 lensed quasars and 10 quasar pairs

We report the results of the STRong lensing Insights from the Dark Energy Survey (STRIDES) follow-up campaign of the late 2017/early 2018 season. We obtained spectra of 65 lensed quasar candidates either with EFOSC2 on the NTT or ESI on Keck, which confirm 10 new gravitationally lensed quasars and 10 quasar pairs with similar spectra, but which do not show a lensing galaxy in DES images. Eight lensed quasars are doubly imaged with source redshifts between 0.99 and 2.90, one is triply imaged by a group (DESJ0345-2545, $z=1.68$), and one is quadruply imaged (quad: DESJ0053-2012, $z=3.8$). Singular isothermal ellipsoid models for the doubles, based on high-resolution imaging from SAMI on SOAR or NIRC2 on Keck, give total magnifications between 3.2 and 5.6, and Einstein radii between 0.49 and 1.97 arcseconds. After spectroscopic follow-up, we extract multi-epoch $grizY$ photometry of confirmed lensed quasars and contaminant quasar+star pairs from the first 4 years of DES data using parametric multi-band modelling, and compare variability in each system&#39;s components. By measuring the reduced $χ^2$ associated with fitting all epochs to the same magnitude, we find a simple cut on the less variable component that retains all confirmed lensed quasars, while removing 94 per cent of contaminant systems with stellar components. Based on our spectroscopic follow-up, this variability information can improve selection of lensed quasars and quasar pairs from 34-45 per cent to 51-70 per cent, with the majority of remaining contaminants being compact star-forming galaxies. Using mock lensed quasar lightcurves we demonstrate that selection based only on variability will over-represent the quad fraction by 10 per cent over a complete DES magnitude-limited sample (excluding microlensing differences), explained by the magnification bias and hence lower luminosity (more variable) sources in quads.