Source author record

M. Lima

M. Lima appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

81works
1topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

81 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: Calibration of Lens Sample Redshift Distributions using Clustering Redshifts with BOSS/eBOSS

We present clustering redshift measurements for Dark Energy Survey (DES) lens sample galaxies to be used in weak gravitational lensing and galaxy clustering studies. To perform this measurement, we cross-correlate with spectroscopic galaxies from the Baryon Acoustic Oscillation Survey (BOSS) and its extension, eBOSS. We validate our methodology in simulations, including a new technique to calibrate systematic errors due to the galaxy clustering bias, finding our method to be generally unbiased in calibrating the mean redshift. We apply our method to the data, and estimate the redshift distribution for eleven different photometrically-selected bins. We find general agreement between clustering redshift and photometric redshift estimates, with differences on the inferred mean redshift to be below $|Δz|=0.01$ in most of the bins. We also test a method to calibrate a width parameter for redshift distributions, which we found necessary to use for some of our samples. Our typical uncertainties on the mean redshift ranged from 0.003 to 0.008, while our uncertainties on the width ranged from 4 to 9\%. We discuss how these results calibrate the photometric redshift distributions used in companion DES Year 3 Results papers.

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: 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: 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: Galaxy-halo connection from galaxy-galaxy lensing

Galaxy-galaxy lensing is a powerful probe of the connection between galaxies and their host dark matter halos, which is important both for galaxy evolution and cosmology. We extend the measurement and modeling of the galaxy-galaxy lensing signal in the recent Dark Energy Survey Year 3 cosmology analysis to the highly nonlinear scales ($\sim 100$ kpc). This extension enables us to study the galaxy-halo connection via a Halo Occupation Distribution (HOD) framework for the two lens samples used in the cosmology analysis: a luminous red galaxy sample (redMaGiC) and a magnitude-limited galaxy sample (MagLim). We find that redMaGiC (MagLim) galaxies typically live in dark matter halos of mass $\log_{10}(M_{h}/M_{\odot}) \approx 13.7$ which is roughly constant over redshift ($13.3-13.5$ depending on redshift). We constrain these masses to $\sim 15\%$, approximately $1.5$ times improvement over previous work. We also constrain the linear galaxy bias more than 5 times better than what is inferred by the cosmological scales only. We find the satellite fraction for redMaGiC (MagLim) to be $\sim 0.1-0.2$ ($0.1-0.3$) with no clear trend in redshift. Our constraints on these halo properties are broadly consistent with other available estimates from previous work, large-scale constraints and simulations. The framework built in this paper will be used for future HOD studies with other galaxy samples and extensions for cosmological analyses.

preprint2022arXiv

Dark Energy Survey Year 3 results: imprints of cosmic voids and superclusters in the Planck CMB lensing map

The CMB lensing signal from cosmic voids and superclusters probes the growth of structure in the low-redshift cosmic web. In this analysis, we cross-correlated the Planck CMB lensing map with voids detected in the Dark Energy Survey Year 3 (Y3) data set ($\sim$5,000 deg$^{2}$), expanding on previous measurements that used Y1 catalogues ($\sim$1,300 deg$^{2}$). Given the increased statistical power compared to Y1 data, we report a $6.6σ$ detection of negative CMB convergence ($κ$) imprints using approximately 3,600 voids detected from a redMaGiC luminous red galaxy sample. However, the measured signal is lower than expected from the MICE N-body simulation that is based on the $Λ$CDM model (parameters $Ω_{\rm m} = 0.25$, $σ_8 = 0.8$), and the discrepancy is associated mostly with the void centre region. Considering the full void lensing profile, we fit an amplitude $A_κ=κ_{\rm DES}/κ_{\rm MICE}$ to a simulation-based template with fixed shape and found a moderate $2σ$ deviation in the signal with $A_κ\approx0.79\pm0.12$. We also examined the WebSky simulation that is based on a Planck 2018 $Λ$CDM cosmology, but the results were even less consistent given the slightly higher matter density fluctuations than in MICE. We then identified superclusters in the DES and the MICE catalogues, and detected their imprints at the $8.4σ$ level; again with a lower-than-expected $A_κ=0.84\pm0.10$ amplitude. The combination of voids and superclusters yields a $10.3σ$ detection with an $A_κ=0.82\pm0.08$ constraint on the CMB lensing amplitude, thus the overall signal is $2.3σ$ weaker than expected from MICE.

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'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 "source" galaxies and magnification biases of nearer "lens" 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: 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

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

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

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.

preprint2022arXiv

Velocity Dispersions of Clusters in the Dark Energy Survey Y3 redMaPPer Catalog

We measure the velocity dispersions of clusters of galaxies selected by the redMaPPer algorithm in the first three years of data from the Dark Energy Survey (DES), allowing us to probe cluster selection and richness estimation, $λ$, in light of cluster dynamics. Our sample consists of 126 clusters with sufficient spectroscopy for individual velocity dispersion estimates. We examine the correlations between cluster velocity dispersion, richness, X-ray temperature and luminosity as well as central galaxy velocity offsets. The velocity dispersion-richness relation exhibits a bimodal distribution. The majority of clusters follow scaling relations between velocity dispersion, richness, and X-ray properties similar to those found for previous samples; however, there is a significant population of clusters with velocity dispersions which are high for their richness. These clusters account for roughly 22\% of the $λ< 70$ systems in our sample, but more than half (55\%) of $λ< 70$ clusters at $z>0.5$. A couple of these systems are hot and X-ray bright as expected for massive clusters with richnesses that appear to have been underestimated, but most appear to have high velocity dispersions for their X-ray properties likely due to line-of-sight structure. These results suggest that projection effects contribute significantly to redMaPPer selection, particularly at higher redshifts and lower richnesses. The redMaPPer determined richnesses for the velocity dispersion outliers are consistent with their X-ray properties, but several are X-ray undetected and deeper data is needed to understand their nature.

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

Cross-correlation of DES Y3 lensing and ACT/${\it Planck}$ thermal Sunyaev Zel'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: cosmology from combined galaxy clustering and lensing -- validation on cosmological simulations

We present a validation of the Dark Energy Survey Year 3 (DES Y3) $3\times2$-point analysis choices by testing them on Buzzard v2.0, a new suite of cosmological simulations that is tailored for the testing and validation of combined galaxy clustering and weak lensing analyses. We show that the Buzzard v2.0 simulations accurately reproduce many important aspects of the DES Y3 data, including photometric redshift and magnitude distributions, and the relevant set of two-point clustering and weak lensing statistics. We then show that our model for the $3\times2$-point data vector is accurate enough to recover the true cosmology in simulated surveys assuming the true redshift distributions for our source and lens samples, demonstrating robustness to uncertainties in the modeling of the non-linear matter power spectrum, non-linear galaxy bias and higher-order lensing corrections. Additionally, we demonstrate for the first time that our photometric redshift calibration methodology, including information from photometry, spectroscopy, clustering cross-correlations, and galaxy-galaxy lensing ratios, is accurate enough to recover the true cosmology in simulated surveys in the presence of realistic photometric redshift uncertainties.

preprint2021arXiv

Dark Energy Survey Year 3 results: Cosmology with peaks using an emulator approach

We constrain the matter density $Ω_{\mathrm{m}}$ and the amplitude of density fluctuations $σ_8$ within the $Λ$CDM cosmological model with shear peak statistics and angular convergence power spectra using mass maps constructed from the first three years of data of the Dark Energy Survey (DES Y3). We use tomographic shear peak statistics, including cross-peaks: peak counts calculated on maps created by taking a harmonic space product of the convergence of two tomographic redshift bins. Our analysis follows a forward-modelling scheme to create a likelihood of these statistics using N-body simulations, using a Gaussian process emulator. We include the following lensing systematics: multiplicative shear bias, photometric redshift uncertainty, and galaxy intrinsic alignment. Stringent scale cuts are applied to avoid biases from unmodelled baryonic physics. We find that the additional non-Gaussian information leads to a tightening of the constraints on the structure growth parameter yielding $S_8~\equiv~σ_8\sqrt{Ω_{\mathrm{m}}/0.3}~=~0.797_{-0.013}^{+0.015}$ (68% confidence limits), with a precision of 1.8%, an improvement of ~38% compared to the angular power spectra only case. The results obtained with the angular power spectra and peak counts are found to be in agreement with each other and no significant difference in $S_8$ is recorded. We find a mild tension of $1.5 \thinspace σ$ between our study and the results from Planck 2018, with our analysis yielding a lower $S_8$. Furthermore, we observe that the combination of angular power spectra and tomographic peak counts breaks the degeneracy between galaxy intrinsic alignment $A_{\mathrm{IA}}$ and $S_8$, improving cosmological constraints. We run a suite of tests concluding that our results are robust and consistent with the results from other studies using DES Y3 data.

preprint2021arXiv

Dark Energy Survey Year 3 Results: Galaxy clustering and systematics treatment for lens galaxy samples

In this work we present the galaxy clustering measurements of the two DES lens galaxy samples: a magnitude-limited sample optimized for the measurement of cosmological parameters, MagLim, and a sample of luminous red galaxies selected with the redMaGiC algorithm. MagLim / redMaGiC sample contains over 10 million / 2.5 million galaxies and is divided into six / five photometric redshift bins spanning the range $z\in[0.20,1.05]$ / $z\in[0.15,0.90]$. Both samples cover 4143 deg$^2$ over which we perform our analysis blind, measuring the angular correlation function with a S/N $\sim 63$ for both samples. In a companion paper (DES Collaboration et al. 2021)), these measurements of galaxy clustering are combined with the correlation functions of cosmic shear and galaxy-galaxy lensing of each sample to place cosmological constraints with a 3$\times$2pt analysis. We conduct a thorough study of the mitigation of systematic effects caused by the spatially varying survey properties and we correct the measurements to remove artificial clustering signals. We employ several decontamination methods with different configurations to ensure the robustness of our corrections and to determine the systematic uncertainty that needs to be considered for the final cosmology analyses. We validate our fiducial methodology using log-normal mocks, showing that our decontamination procedure induces biases no greater than $0.5σ$ in the $(Ω_m, b)$ plane, where $b$ is galaxy bias. We demonstrate that failure to remove the artificial clustering would introduce strong biases up to $\sim 7 σ$ in $Ω_m$ and of more than $4 σ$ in galaxy bias.

preprint2021arXiv

DES Y1 results: Splitting growth and geometry to test $Λ$CDM

We analyze Dark Energy Survey (DES) data to constrain a cosmological model where a subset of parameters -- focusing on $Ω_m$ -- are split into versions associated with structure growth (e.g. $Ω_m^{\rm grow}$) and expansion history (e.g. $Ω_m^{\rm geo}$). Once the parameters have been specified for the $Λ$CDM cosmological model, which includes general relativity as a theory of gravity, it uniquely predicts the evolution of both geometry (distances) and the growth of structure over cosmic time. Any inconsistency between measurements of geometry and growth could therefore indicate a breakdown of that model. Our growth-geometry split approach therefore serves as both a (largely) model-independent test for beyond-$Λ$CDM physics, and as a means to characterize how DES observables provide cosmological information. We analyze the same multi-probe DES data as arXiv:1811.02375 : DES Year 1 (Y1) galaxy clustering and weak lensing, which are sensitive to both growth and geometry, as well as Y1 BAO and Y3 supernovae, which probe geometry. We additionally include external geometric information from BOSS DR12 BAO and a compressed Planck 2015 likelihood, and external growth information from BOSS DR12 RSD. We find no significant disagreement with $Ω_m^{\rm grow}=Ω_m^{\rm geo}$. When DES and external data are analyzed separately, degeneracies with neutrino mass and intrinsic alignments limit our ability to measure $Ω_m^{\rm grow}$, but combining DES with external data allows us to constrain both growth and geometric quantities. We also consider a parameterization where we split both $Ω_m$ and $w$, but find that even our most constraining data combination is unable to separately constrain $Ω_m^{\rm grow}$ and $w^{\rm grow}$. Relative to $Λ$CDM, splitting growth and geometry weakens bounds on $σ_8$ but does not alter constraints on $h$.

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 "needle in a haystack" 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

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 'contamination' 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 'non-Ia' 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 'BEAMS with Bias Correction' (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 DES view of the Eridanus supervoid and the CMB Cold Spot

The Cold Spot is a puzzling large-scale feature in the Cosmic Microwave Background temperature maps and its origin has been subject to active debate. As an important foreground structure at low redshift, the Eridanus supervoid was recently detected, but it was subsequently determined that, assuming the standard $Λ$CDM model, only about 10-20$\%$ of the observed temperature depression can be accounted for via its Integrated Sachs-Wolfe imprint. However, $R\gtrsim100~h^{-1}\mathrm{Mpc}$ supervoids elsewhere in the sky have shown ISW imprints $A_{\mathrm{ISW}}\approx5.2\pm1.6$ times stronger than expected from $Λ$CDM ($A_{\mathrm{ISW}}=1$), which warrants further inspection. Using the Year-3 redMaGiC catalogue of luminous red galaxies from the Dark Energy Survey, here we confirm the detection of the Eridanus supervoid as a significant under-density in the Cold Spot's direction at $z<0.2$. We also show, with $\mathrm{S/N}\gtrsim5$ significance, that the Eridanus supervoid appears as the most prominent large-scale under-density in the dark matter mass maps that we reconstructed from DES Year-3 gravitational lensing data. While we report no significant anomalies, an interesting aspect is that the amplitude of the lensing signal from the Eridanus supervoid at the Cold Spot centre is about $30\%$ lower than expected from similar peaks found in N-body simulations based on the standard $Λ$CDM model with parameters $Ω_{\rm m} = 0.279$ and $σ_8 = 0.82$. Overall, our results confirm the causal relation between these individually rare structures in the cosmic web and in the CMB, motivating more detailed future surveys in the Cold Spot region.

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 "local" 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'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'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's evolution. We confirm prior findings of negative evolution in the SMHM relation'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'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.

preprint2021arXiv

The WaZP galaxy cluster sample of the Dark Energy Survey Year 1

We present a new (2+1)D galaxy cluster finder based on photometric redshifts called Wavelet Z Photometric (WaZP) applied to DES first year (Y1A1) data. The results are compared to clusters detected by the South Pole Telescope (SPT) survey and the redMaPPer cluster finder, the latter based on the same photometric data. WaZP searches for clusters in wavelet-based density maps of galaxies selected in photometric redshift space without any assumption on the cluster galaxy populations. The comparison to other cluster samples was performed with a matching algorithm based on angular proximity and redshift difference of the clusters. It led to the development of a new approach to match two optical cluster samples, following an iterative approach to minimize incorrect associations. The WaZP cluster finder applied to DES Y1A1 galaxy survey (1,511 deg2 up to i=23 mag) led to the detection of 60,547 galaxy clusters with redshifts 0.05<z<0.9 and richness Ngals>5. Considering the overlapping regions and redshift ranges between the DES Y1A1 and SPT cluster surveys, all SZ based SPT clusters are recovered by the WaZP sample.The comparison between WaZP and redMaPPer cluster samples showed an excellent overall agreement for clusters with richness Ngals (λ for redMaPPer) greater than 25 (20), with 95% recovery on both directions. Based on the cluster cross-match we explorethe relative fragmentation of the two cluster samples and investigate the possible signatures of unmatched clusters.

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

Constraining Radio Mode Feedback in Galaxy Clusters with the Cluster Radio AGN Properties to z$\sim$1

We study the properties of the Sydney University Molonglo Sky Survey (SUMSS) 843~MHz radio AGN population in galaxy clusters from two large catalogs created using the Dark Energy Survey (DES): $\sim$11,800 optically selected RM-Y3 and $\sim$1,000 X-ray selected MARD-Y3 clusters. We show that cluster radio loud AGN are highly concentrated around cluster centers to $z\sim1$. We measure the halo occupation number for cluster radio AGN above a threshold luminosity, finding that the number of radio AGN per cluster increases with cluster halo mass as $N\propto M^{1.2\pm0.1}$ ($N\propto M^{0.68\pm0.34}$) for the RM-Y3 (MARD-Y3) sample. Together, these results indicate that radio mode feedback is favoured in more massive galaxy clusters. Using optical counterparts for these sources, we demonstrate weak redshift evolution in the host broad band colors and the radio luminosity at fixed host galaxy stellar mass. We use the redshift evolution in radio luminosity to break the degeneracy between density and luminosity evolution scenarios in the redshift trend of the radio AGN luminosity function (LF). The LF exhibits a redshift trend of the form $(1+z)^γ$ in density and luminosity, respectively, of $γ_{\rm D}=3.0\pm0.4$ and $γ_{\rm P}=0.21\pm0.15$ in the RM-Y3 sample, and $γ_{\rm D}=2.6\pm0.7$ and $γ_{\rm P}=0.31\pm0.15$ in MARD-Y3. We discuss the physical drivers of radio mode feedback in cluster AGN, and we use the cluster radio galaxy LF to estimate the average radio-mode feedback energy as a function of cluster mass and redshift and compare it to the core ($<0.1R_{500}$) X-ray radiative losses for clusters at $z<1$.

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

First Cosmology Results using Type Ia Supernovae from the Dark Energy Survey: The Effect of Host Galaxy Properties on Supernova Luminosity

We present improved photometric measurements for the host galaxies of 206 spectroscopically confirmed type Ia supernovae discovered by the Dark Energy Survey Supernova Program (DES-SN) and used in the first DES-SN cosmological analysis. Fitting spectral energy distributions to the $griz$ photometric measurements of the DES-SN host galaxies, we derive stellar masses and star-formation rates. For the DES-SN sample, when considering a 5D ($z$, $x_1$, $c$, $α$, $β$) bias correction, we find evidence of a Hubble residual `mass step', where SNe Ia in high mass galaxies ($>10^{10} \textrm{M}_{\odot}$) are intrinsically more luminous (after correction) than their low mass counterparts by $γ=0.040\pm0.019$mag. This value is larger by $0.031$mag than the value found in the first DES-SN cosmological analysis. This difference is due to a combination of updated photometric measurements and improved star formation histories and is not from host-galaxy misidentification. When using a 1D (redshift-only) bias correction the inferred mass step is larger, with $γ=0.066\pm0.020$mag. The 1D-5D $γ$ difference for DES-SN is $0.026\pm0.009$mag. We show that this difference is due to a strong correlation between host galaxy stellar mass and the $x_1$ component of the 5D distance-bias correction. To better understand this effect, we include an intrinsic correlation between light-curve width and stellar mass in simulated SN Ia samples. We show that a 5D fit recovers $γ$ with $-9$mmag bias compared to a $+2$mmag bias for a 1D fit. This difference can explain part of the discrepancy seen in the data. Improvements in modeling correlations between galaxy properties and SN is necessary to determine the implications for $γ$ and ensure unbiased precision estimates of the dark energy equation-of-state as we enter the era of LSST.

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

OzDES multi-object fibre spectroscopy for the Dark Energy Survey: Results and second data release

We present a description of the Australian Dark Energy Survey (OzDES) and summarise the results from its six years of operations. Using the 2dF fibre positioner and AAOmega spectrograph on the 3.9-metre Anglo-Australian Telescope, OzDES has monitored 771 AGN, classified hundreds of supernovae, and obtained redshifts for thousands of galaxies that hosted a transient within the 10 deep fields of the Dark Energy Survey. We also present the second OzDES data release, containing the redshifts of almost 30,000 sources, some as faint as $r_{\mathrm AB}=24$ mag, and 375,000 individual spectra. These data, in combination with the time-series photometry from the Dark Energy Survey, will be used to measure the expansion history of the Universe out to $z\sim1.2$ and the masses of hundreds of black holes out to $z\sim4$. OzDES is a template for future surveys that combine simultaneous monitoring of targets with wide-field imaging cameras and wide-field multi-object spectrographs.

preprint2020arXiv

Quasar Accretion Disk Sizes from Continuum Reverberation Mapping in the DES Standard Star Fields

Measurements of the physical properties of accretion disks in active galactic nuclei are important for better understanding the growth and evolution of supermassive black holes. We present the accretion disk sizes of 22 quasars from continuum reverberation mapping with data from the Dark Energy Survey (DES) standard star fields and the supernova C fields. We construct continuum lightcurves with the \textit{griz} photometry that span five seasons of DES observations. These data sample the time variability of the quasars with a cadence as short as one day, which corresponds to a rest frame cadence that is a factor of a few higher than most previous work. We derive time lags between bands with both JAVELIN and the interpolated cross-correlation function method, and fit for accretion disk sizes using the JAVELIN Thin Disk model. These new measurements include disks around black holes with masses as small as $\sim10^7$ $M_{\odot}$, which have equivalent sizes at 2500Å\, as small as $\sim 0.1$ light days in the rest frame. We find that most objects have accretion disk sizes consistent with the prediction of the standard thin disk model when we take disk variability into account. We have also simulated the expected yield of accretion disk measurements under various observational scenarios for the Large Synoptic Survey Telescope Deep Drilling Fields. We find that the number of disk measurements would increase significantly if the default cadence is changed from three days to two days or one day.

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 "effective" 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'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

Studying Type II supernovae as cosmological standard candles using the Dark Energy Survey

Despite vast improvements in the measurement of the cosmological parameters, the nature of dark energy and an accurate value of the Hubble constant (H$_0$) in the Hubble-Lemaître law remain unknown. To break the current impasse, it is necessary to develop as many independent techniques as possible, such as the use of Type II supernovae (SNe II). The goal of this paper is to demonstrate the utility of SNe II for deriving accurate extragalactic distances, which will be an asset for the next generation of telescopes where more-distant SNe II will be discovered. More specifically, we present a sample from the Dark Energy Survey Supernova Program (DES-SN) consisting of 15 SNe II with photometric and spectroscopic information spanning a redshift range up to 0.35. Combining our DES SNe with publicly available samples, and using the standard candle method (SCM), we construct the largest available Hubble diagram with SNe II in the Hubble flow (70 SNe II) and find an observed dispersion of 0.27 mag. We demonstrate that adding a colour term to the SN II standardisation does not reduce the scatter in the Hubble diagram. Although SNe II are viable as distance indicators, this work points out important issues for improving their utility as independent extragalactic beacons: find new correlations, define a more standard subclass of SNe II, construct new SN II templates, and dedicate more observing time to high-redshift SNe II. Finally, for the first time, we perform simulations to estimate the redshift-dependent distance-modulus bias due to selection effects.

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' galaxies that host multiple SNe Ia 'siblings'. 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 Diffuse Light Envelope of Luminous Red Galaxies

We use a stacking method to study the radial light profiles of luminous red galaxies (LRGs) at redshift $\sim 0.62$ and $\sim 0.25$, out to a radial range of 200 kpc. We do not find noticeable evolution of the profiles at the two redshifts. The LRG profiles appear to be well approximated by a single Sersic profile, although some excess light can be seen outside 60 kpc. We quantify the excess light by measuring the integrated flux and find that the excess is about 10\% -- a non-dominant but still nonnegligible component.

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

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

Monte Carlo Control Loops for cosmic shear cosmology with DES Year 1

Weak lensing by large-scale structure is a powerful probe of cosmology and of the dark universe. This cosmic shear technique relies on the accurate measurement of the shapes and redshifts of background galaxies and requires precise control of systematic errors. The Monte Carlo Control Loops (MCCL) is a forward modelling method designed to tackle this problem. It relies on the Ultra Fast Image Generator (UFig) to produce simulated images tuned to match the target data statistically, followed by calibrations and tolerance loops. We present the first end-to-end application of this method, on the Dark Energy Survey (DES) Year 1 wide field imaging data. We simultaneously measure the shear power spectrum $C_{\ell}$ and the redshift distribution $n(z)$ of the background galaxy sample. The method includes maps of the systematic sources, Point Spread Function (PSF), an Approximate Bayesian Computation (ABC) inference of the simulation model parameters, a shear calibration scheme, and the fast estimation of the covariance matrix. We find a close statistical agreement between the simulations and the DES Y1 data using an array of diagnostics. In a non-tomographic setting, we derive a set of $C_\ell$ and $n(z)$ curves that encode the cosmic shear measurement, as well as the systematic uncertainty. Following a blinding scheme, we measure the combination of $Ω_m$, $σ_8$, and intrinsic alignment amplitude $A_{\rm{IA}}$, defined as $S_8D_{\rm{IA}} = σ_8(Ω_m/0.3)^{0.5}D_{\rm{IA}}$, where $D_{\rm{IA}}=1-0.11(A_{\rm{IA}}-1)$. We find $S_8D_{\rm{IA}}=0.895^{+0.054}_{-0.039}$, where systematics are at the level of roughly 60\% of the statistical errors. We discuss these results in the context of earlier cosmic shear analyses of the DES Y1 data. Our findings indicate that this method and its fast runtime offer good prospects for cosmic shear measurements with future wide-field surveys.

preprint2019arXiv

Observation and Confirmation of Nine Strong Lensing Systems in Dark Energy Survey Year 1 Data

We describe the observation and confirmation of \nbconfirmtext\ new strong gravitational lenses discovered in Year 1 data from the Dark Energy Survey (DES). We created candidate lists based on a) galaxy group and cluster samples and b) photometrically selected galaxy samples. We selected 46 candidates through visual inspection and then used the Gemini Multi-Object Spectrograph (GMOS) at the Gemini South telescope to acquire spectroscopic follow-up of 21 of these candidates. Through analysis of this spectroscopic follow-up data, we confirmed nine new lensing systems and rejected 2 candidates, but the analysis was inconclusive on 10 candidates. For each of the confirmed systems, we report measured spectroscopic properties, estimated \einsteinradiussub, and estimated enclosed masses. The sources that we targeted have an i-band surface brightness range of iSB ~ 22 - 24 mag arcsec^2 and a spectroscopic redshift range of zspec ~0.8 - 2.6. The lens galaxies have a photometric redshift range of zlens ~ 0.3 - 0.7. The lensing systems range in image-lens separation 2 - 9 arcsec and in enclosed mass 10^12 - 10^13 Msol.

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

preprint2016arXiv

A Dark Energy Camera Search for an Optical Counterpart to the First Advanced LIGO Gravitational Wave Event GW150914

We report initial results of a deep search for an optical counterpart to the gravitational wave event GW150914, the first trigger from the Advanced LIGO gravitational wave detectors. We used the Dark Energy Camera (DECam) to image a 102 deg$^2$ area, corresponding to 38% of the initial trigger high-probability sky region and to 11% of the revised high-probability region. We observed in i and z bands at 4-5, 7, and 24 days after the trigger. The median $5σ$ point-source limiting magnitudes of our search images are i=22.5 and z=21.8 mag. We processed the images through a difference-imaging pipeline using templates from pre-existing Dark Energy Survey data and publicly available DECam data. Due to missing template observations and other losses, our effective search area subtends 40 deg$^{2}$, corresponding to 12% total probability in the initial map and 3% of the final map. In this area, we search for objects that decline significantly between days 4-5 and day 7, and are undetectable by day 24, finding none to typical magnitude limits of i= 21.5,21.1,20.1 for object colors (i-z)=1,0,-1, respectively. Our search demonstrates the feasibility of a dedicated search program with DECam and bodes well for future research in this emerging field.

preprint2016arXiv

A Dark Energy Camera Search for Missing Supergiants in the LMC After the Advanced LIGO Gravitational Wave Event GW150914

The collapse of the core of a star is expected to produce gravitational radiation. While this process will usually produce a luminous supernova, the optical signatue could be subluminous and a direct collapse to a black hole, with the star just disappearing, is possible. The gravitational wave event GW150914 reported by the LIGO Virgo Collaboration (LVC) on 2015 September 16, was detected by a burst analysis and whose high probability spatial localization included the Large Magellanic Cloud. Shortly after the announcement of the event, we used the Dark Energy Camera to observe 102 deg$^2$ of the localization area, including a 38 deg$^2$ area centered on the LMC. Using a catalog of 152 LMC luminous red supergiants, candidates to undergo a core collapse without a visible supernova, we find that the positions of 144 of these are inside our images, and that all are detected - none have disappeared. There are other classes of candidates: we searched existing catalogs of red supergiants, yellow supergiants, Wolf-Rayet stars, and luminous blue variable stars, recovering all that were inside the imaging area. Based on our observations, we conclude that it is unlikely that GW150914 was caused by the core collapse of a supergiant in the LMC, consistent with the LIGO Collaboration analyses of the gravitational wave form as best described by a binary black hole merger. We discuss how to generalize this search for future very nearby core collapse candidates.

preprint2016arXiv

A DECam Search for an Optical Counterpart to the LIGO Gravitational Wave Event GW151226

We report the results of a Dark Energy Camera (DECam) optical follow-up of the gravitational wave (GW) event GW151226, discovered by the Advanced LIGO detectors. Our observations cover 28.8 deg$^2$ of the localization region in the $i$ and $z$ bands (containing 3% of the BAYESTAR localization probability), starting 10 hours after the event was announced and spanning four epochs at $2-24$ days after the GW detection. We achieve $5σ$ point-source limiting magnitudes of $i\approx21.7$ and $z\approx21.5$, with a scatter of $0.4$ mag, in our difference images. Given the two day delay, we search this area for a rapidly declining optical counterpart with $\gtrsim 3σ$ significance steady decline between the first and final observations. We recover four sources that pass our selection criteria, of which three are cataloged AGN. The fourth source is offset by $5.8$ arcsec from the center of a galaxy at a distance of 187 Mpc, exhibits a rapid decline by $0.5$ mag over $4$ days, and has a red color of $i-z\approx 0.3$ mag. These properties roughly match the expectations for a kilonova. However, this source was detected several times, starting $94$ days prior to GW151226, in the Pan-STARRS Survey for Transients (dubbed as PS15cdi) and is therefore unrelated to the GW event. Given its long-term behavior, PS15cdi is likely a Type IIP supernova that transitioned out of its plateau phase during our observations, mimicking a kilonova-like behavior. We comment on the implications of this detection for contamination in future optical follow-up observations.

preprint2016arXiv

CMB lensing tomography with the DES Science Verification galaxies

We measure the cross-correlation between the galaxy density in the Dark Energy Survey (DES) Science Verification data and the lensing of the cosmic microwave background (CMB) as reconstructed with the Planck satellite and the South Pole Telescope (SPT). When using the DES main galaxy sample over the full redshift range $0.2 < z < 1.2$, a cross-correlation signal is detected at $6 σ$ and $4σ$ with SPT and Planck respectively. We then divide the DES galaxies into five photometric redshift bins, finding significant ($>$$2 σ$) detections in all bins. Comparing to the fiducial Planck cosmology, we find the redshift evolution of the signal matches expectations, although the amplitude is consistently lower than predicted across redshift bins. We test for possible systematics that could affect our result and find no evidence for significant contamination. Finally, we demonstrate how these measurements can be used to constrain the growth of structure across cosmic time. We find the data are fit by a model in which the amplitude of structure in the $z<1.2$ universe is $0.73 \pm 0.16$ times as large as predicted in the LCDM Planck cosmology, a $1.7σ$ deviation.

preprint2016arXiv

Comparing Dark Energy Survey and HST-CLASH observations of the galaxy cluster RXC J2248.7-4431: implications for stellar mass versus dark matter

We derive the stellar mass fraction in the galaxy cluster RXC J2248.7-4431 observed with the Dark Energy Survey (DES) during the Science Verification period. We compare the stellar mass results from DES (five filters) with those from the Hubble Space Telescope Cluster Lensing And Supernova Survey (CLASH; 17 filters). When the cluster spectroscopic redshift is assumed, we show that stellar masses from DES can be estimated within 25% of CLASH values. We compute the stellar mass contribution coming from red and blue galaxies, and study the relation between stellar mass and the underlying dark matter using weak lensing studies with DES and CLASH. An analysis of the radial profiles of the DES total and stellar mass yields a stellar-to-total fraction of f*=(6.8+-1.7)x10^-3 within a radius of r_200c~2 Mpc. Our analysis also includes a comparison of photometric redshifts and star/galaxy separation efficiency for both data sets. We conclude that space-based small field imaging can be used to calibrate the galaxy properties in DES for the much wider field of view. The technique developed to derive the stellar mass fraction in galaxy clusters can be applied to the ~100 000 clusters that will be observed within this survey and yield important information about galaxy evolution.

preprint2016arXiv

Cosmic Shear Measurements with DES Science Verification Data

We present measurements of weak gravitational lensing cosmic shear two-point statistics using Dark Energy Survey Science Verification data. We demonstrate that our results are robust to the choice of shear measurement pipeline, either ngmix or im3shape, and robust to the choice of two-point statistic, including both real and Fourier-space statistics. Our results pass a suite of null tests including tests for B-mode contamination and direct tests for any dependence of the two-point functions on a set of 16 observing conditions and galaxy properties, such as seeing, airmass, galaxy color, galaxy magnitude, etc. We furthermore use a large suite of simulations to compute the covariance matrix of the cosmic shear measurements and assign statistical significance to our null tests. We find that our covariance matrix is consistent with the halo model prediction, indicating that it has the appropriate level of halo sample variance. We compare the same jackknife procedure applied to the data and the simulations in order to search for additional sources of noise not captured by the simulations. We find no statistically significant extra sources of noise in the data. The overall detection significance with tomography for our highest source density catalog is 9.7sigma. Cosmological constraints from the measurements in this work are presented in a companion paper (DES et al. 2015).

preprint2016arXiv

Cosmic Voids and Void Lensing in the Dark Energy Survey Science Verification Data

Galaxies and their dark matter halos populate a complicated filamentary network around large, nearly empty regions known as cosmic voids. Cosmic voids are usually identified in spectroscopic galaxy surveys, where 3D information about the large-scale structure of the Universe is available. Although an increasing amount of photometric data is being produced, its potential for void studies is limited since photometric redshifts induce line-of-sight position errors of $\sim50$ Mpc/$h$ or more that can render many voids undetectable. In this paper we present a new void finder designed for photometric surveys, validate it using simulations, and apply it to the high-quality photo-$z$ redMaGiC galaxy sample of the Dark Energy Survey Science Verification (DES-SV) data. The algorithm works by projecting galaxies into 2D slices and finding voids in the smoothed 2D galaxy density field of the slice. Fixing the line-of-sight size of the slices to be at least twice the photo-$z$ scatter, the number of voids found in these projected slices of simulated spectroscopic and photometric galaxy catalogs is within 20% for all transverse void sizes, and indistinguishable for the largest voids of radius $\sim 70$ Mpc/$h$ and larger. The positions, radii, and projected galaxy profiles of photometric voids also accurately match the spectroscopic void sample. Applying the algorithm to the DES-SV data in the redshift range $0.2<z<0.8$, we identify 87 voids with comoving radii spanning the range 18-120 Mpc/$h$, and carry out a stacked weak lensing measurement. With a significance of $4.4σ$, the lensing measurement confirms the voids are truly underdense in the matter field and hence not a product of Poisson noise, tracer density effects or systematics in the data. It also demonstrates, for the first time in real data, the viability of void lensing studies in photometric surveys.

preprint2016arXiv

Cosmology constraints from shear peak statistics in Dark Energy Survey Science Verification data

Shear peak statistics has gained a lot of attention recently as a practical alternative to the two point statistics for constraining cosmological parameters. We perform a shear peak statistics analysis of the Dark Energy Survey (DES) Science Verification (SV) data, using weak gravitational lensing measurements from a 139 deg$^2$ field. We measure the abundance of peaks identified in aperture mass maps, as a function of their signal-to-noise ratio, in the signal-to-noise range $0<\mathcal S / \mathcal N<4$. To predict the peak counts as a function of cosmological parameters we use a suite of $N$-body simulations spanning 158 models with varying $Ω_{\rm m}$ and $σ_8$, fixing $w = -1$, $Ω_{\rm b} = 0.04$, $h = 0.7$ and $n_s=1$, to which we have applied the DES SV mask and redshift distribution. In our fiducial analysis we measure $σ_{8}(Ω_{\rm m}/0.3)^{0.6}=0.77 \pm 0.07$, after marginalising over the shear multiplicative bias and the error on the mean redshift of the galaxy sample. We introduce models of intrinsic alignments, blending, and source contamination by cluster members. These models indicate that peaks with $\mathcal S / \mathcal N>4$ would require significant corrections, which is why we do not include them in our analysis. We compare our results to the cosmological constraints from the two point analysis on the SV field and find them to be in good agreement in both the central value and its uncertainty. We discuss prospects for future peak statistics analysis with upcoming DES data.

preprint2016arXiv

Detection of the kinematic Sunyaev-Zel'dovich effect with DES Year 1 and SPT

We detect the kinematic Sunyaev-Zel'dovich (kSZ) effect with a statistical significance of $4.2 σ$ by combining a cluster catalogue derived from the first year data of the Dark Energy Survey (DES) with CMB temperature maps from the South Pole Telescope Sunyaev-Zel'dovich (SPT-SZ) Survey. This measurement is performed with a differential statistic that isolates the pairwise kSZ signal, providing the first detection of the large-scale, pairwise motion of clusters using redshifts derived from photometric data. By fitting the pairwise kSZ signal to a theoretical template we measure the average central optical depth of the cluster sample, $\barτ_e = (3.75 \pm 0.89)\cdot 10^{-3}$. We compare the extracted signal to realistic simulations and find good agreement with respect to the signal-to-noise, the constraint on $\barτ_e$, and the corresponding gas fraction. High-precision measurements of the pairwise kSZ signal with future data will be able to place constraints on the baryonic physics of galaxy clusters, and could be used to probe gravity on scales $ \gtrsim 100$ Mpc.

preprint2016arXiv

Galaxy bias from the Dark Energy Survey Science Verification data: combining galaxy density maps and weak lensing maps

We measure the redshift evolution of galaxy bias for a magnitude-limited galaxy sample by combining the galaxy density maps and weak lensing shear maps for a $\sim$116 deg$^{2}$ area of the Dark Energy Survey (DES) Science Verification data. This method was first developed in Amara et al. (2012) and later re-examined in a companion paper (Pujol et al. 2016) with rigorous simulation tests and analytical treatment of tomographic measurements. In this work we apply this method to the DES SV data and measure the galaxy bias for a i$<$22.5 galaxy sample. We find the galaxy bias and 1$σ$ error bars in 4 photometric redshift bins to be 1.12$\pm$0.19 (z=0.2-0.4), 0.97$\pm$0.15 (z=0.4-0.6), 1.38$\pm$0.39 (z=0.6-0.8)), and 1.45$\pm$0.56 (z=0.8-1.0). These measurements are consistent at the 2$σ$ level with measurements on the same dataset using galaxy clustering and cross-correlation of galaxies with CMB lensing, with most of the redshift bins consistent within the 1σ error bars. In addition, our method provides the only $σ_8$-independent constraint among the three. We forward-model the main observational effects using mock galaxy catalogs by including shape noise, photo-z errors and masking effects. We show that our bias measurement from the data is consistent with that expected from simulations. With the forthcoming full DES data set, we expect this method to provide additional constraints on the galaxy bias measurement from more traditional methods. Furthermore, in the process of our measurement, we build up a 3D mass map that allows further exploration of the dark matter distribution and its relation to galaxy evolution.

preprint2016arXiv

Joint Measurement of Lensing-Galaxy Correlations Using SPT and DES SV Data

We measure the correlation of galaxy lensing and cosmic microwave background lensing with a set of galaxies expected to trace the matter density field. The measurements are performed using pre-survey Dark Energy Survey (DES) Science Verification optical imaging data and millimeter-wave data from the 2500 square degree South Pole Telescope Sunyaev-Zel'dovich (SPT-SZ) survey. The two lensing-galaxy correlations are jointly fit to extract constraints on cosmological parameters, constraints on the redshift distribution of the lens galaxies, and constraints on the absolute shear calibration of DES galaxy lensing measurements. We show that an attractive feature of these fits is that they are fairly insensitive to the clustering bias of the galaxies used as matter tracers. The measurement presented in this work confirms that DES and SPT data are consistent with each other and with the currently favored $Λ$CDM cosmological model. It also demonstrates that joint lensing-galaxy correlation measurement considered here contains a wealth of information that can be extracted using current and future surveys.

preprint2016arXiv

Localization and broadband follow-up of the gravitational-wave transient GW150914

A gravitational-wave (GW) transient was identified in data recorded by the Advanced Laser Interferometer Gravitational-wave Observatory (LIGO) detectors on 2015 September 14. The event, initially designated G184098 and later given the name GW150914, is described in detail elsewhere. By prior arrangement, preliminary estimates of the time, significance, and sky location of the event were shared with 63 teams of observers covering radio, optical, near-infrared, X-ray, and gamma-ray wavelengths with ground- and space-based facilities. In this Letter we describe the low-latency analysis of the GW data and present the sky localization of the first observed compact binary merger. We summarize the follow-up observations reported by 25 teams via private Gamma-ray Coordinates Network circulars, giving an overview of the participating facilities, the GW sky localization coverage, the timeline and depth of the observations. As this event turned out to be a binary black hole merger, there is little expectation of a detectable electromagnetic (EM) signature. Nevertheless, this first broadband campaign to search for a counterpart of an Advanced LIGO source represents a milestone and highlights the broad capabilities of the transient astronomy community and the observing strategies that have been developed to pursue neutron star binary merger events. Detailed investigations of the EM data and results of the EM follow-up campaign are being disseminated in papers by the individual teams.

preprint2016arXiv

No galaxy left behind: accurate measurements with the faintest objects in the Dark Energy Survey

Accurate statistical measurement with large imaging surveys has traditionally required throwing away a sizable fraction of the data. This is because most measurements have have relied on selecting nearly complete samples, where variations in the composition of the galaxy population with seeing, depth, or other survey characteristics are small. We introduce a new measurement method that aims to minimize this wastage, allowing precision measurement for any class of stars or galaxies detectable in an imaging survey. We have implemented our proposal in Balrog, a software package which embeds fake objects in real imaging in order to accurately characterize measurement biases. We demonstrate this technique with an angular clustering measurement using Dark Energy Survey (DES) data. We first show that recovery of our injected galaxies depends on a wide variety of survey characteristics in the same way as the real data. We then construct a flux-limited sample of the faintest galaxies in DES, chosen specifically for their sensitivity to depth and seeing variations. Using the synthetic galaxies as randoms in the standard Landy-Szalay correlation function estimator suppresses the effects of variable survey selection by at least two orders of magnitude. With this correction, our measured angular clustering is found to be in excellent agreement with that of a matched sample drawn from much deeper, higher-resolution space-based Cosmological Evolution Survey (COSMOS) imaging; over angular scales of $0.004^{\circ} < θ< 0.2^{\circ}$, we find a best-fit scaling amplitude between the DES and COSMOS measurements of $1.00 \pm 0.09$. We expect this methodology to be broadly useful for extending the statistical reach of measurements in a wide variety of coming imaging surveys.

preprint2016arXiv

SDSS-IV eBOSS emission-line galaxy pilot survey

The Sloan Digital Sky Survey IV extended Baryonic Oscillation Spectroscopic Survey (SDSS-IV/eBOSS) will observe 195,000 emission-line galaxies (ELGs) to measure the Baryonic Acoustic Oscillation standard ruler (BAO) at redshift 0.9. To test different ELG selection algorithms, 9,000 spectra were observed with the SDSS spectrograph as a pilot survey based on data from several imaging surveys. First, using visual inspection and redshift quality flags, we show that the automated spectroscopic redshifts assigned by the pipeline meet the quality requirements for a reliable BAO measurement. We also show the correlations between sky emission, signal-to-noise ratio in the emission lines, and redshift error. Then we provide a detailed description of each target selection algorithm we tested and compare them with the requirements of the eBOSS experiment. As a result, we provide reliable redshift distributions for the different target selection schemes we tested. Finally, we determine an target selection algorithms that is best suited to be applied on DECam photometry because they fulfill the eBOSS survey efficiency requirements.

preprint2016arXiv

Supplement: Localization and broadband follow-up of the gravitational-wave transient GW150914

This Supplement provides supporting material for arXiv:1602.08492 . We briefly summarize past electromagnetic (EM) follow-up efforts as well as the organization and policy of the current EM follow-up program. We compare the four probability sky maps produced for the gravitational-wave transient GW150914, and provide additional details of the EM follow-up observations that were performed in the different bands.

preprint2016arXiv

Testing the lognormality of the galaxy and weak lensing convergence distributions from Dark Energy Survey maps

It is well known that the probability distribution function (PDF) of galaxy density contrast is approximately lognormal; whether the PDF of mass fluctuations derived from weak lensing convergence (kappa_WL) is lognormal is less well established. We derive PDFs of the galaxy and projected matter density distributions via the Counts in Cells (CiC) method. We use maps of galaxies and weak lensing convergence produced from the Dark Energy Survey (DES) Science Verification data over 139 deg^2. We test whether the underlying density contrast is well described by a lognormal distribution for the galaxies, the convergence and their joint PDF. We confirm that the galaxy density contrast distribution is well modeled by a lognormal PDF convolved with Poisson noise at angular scales from 10-40 arcmin (corresponding to physical scales of 3-10 Mpc). We note that as kappa_WL is a weighted sum of the mass fluctuations along the line of sight, its PDF is expected to be only approximately lognormal. We find that the kappa_WL distribution is well modeled by a lognormal PDF convolved with Gaussian shape noise at scales between 10 and 20 arcmin, with a best-fit chi^2/DOF of 1.11 compared to 1.84 for a Gaussian model, corresponding to p-values 0.35 and 0.07 respectively, at a scale of 10 arcmin. Above 20 arcmin a simple Gaussian model is sufficient. The joint PDF is also reasonably fitted by a bivariate lognormal. As a consistency check we compare the variances derived from the lognormal modelling with those directly measured via CiC. Our methods are validated against maps from the MICE Grand Challenge N-body simulation.

preprint2016arXiv

The Dark Energy Survey: more than dark energy - an overview

This overview article describes the legacy prospect and discovery potential of the Dark Energy Survey (DES) beyond cosmological studies, illustrating it with examples from the DES early data. DES is using a wide-field camera (DECam) on the 4m Blanco Telescope in Chile to image 5000 sq deg of the sky in five filters (grizY). By its completion the survey is expected to have generated a catalogue of 300 million galaxies with photometric redshifts and 100 million stars. In addition, a time-domain survey search over 27 sq deg is expected to yield a sample of thousands of Type Ia supernovae and other transients. The main goals of DES are to characterise dark energy and dark matter, and to test alternative models of gravity; these goals will be pursued by studying large scale structure, cluster counts, weak gravitational lensing and Type Ia supernovae. However, DES also provides a rich data set which allows us to study many other aspects of astrophysics. In this paper we focus on additional science with DES, emphasizing areas where the survey makes a difference with respect to other current surveys. The paper illustrates, using early data (from `Science Verification', and from the first, second and third seasons of observations), what DES can tell us about the solar system, the Milky Way, galaxy evolution, quasars, and other topics. In addition, we show that if the cosmological model is assumed to be Lambda+ Cold Dark Matter (LCDM) then important astrophysics can be deduced from the primary DES probes. Highlights from DES early data include the discovery of 34 Trans Neptunian Objects, 17 dwarf satellites of the Milky Way, one published z > 6 quasar (and more confirmed) and two published superluminous supernovae (and more confirmed).

preprint2016arXiv

The DES Science Verification Weak Lensing Shear Catalogues

We present weak lensing shear catalogues for 139 square degrees of data taken during the Science Verification (SV) time for the new Dark Energy Camera (DECam) being used for the Dark Energy Survey (DES). We describe our object selection, point spread function estimation and shear measurement procedures using two independent shear pipelines, IM3SHAPE and NGMIX, which produce catalogues of 2.12 million and 3.44 million galaxies respectively. We detail a set of null tests for the shear measurements and find that they pass the requirements for systematic errors at the level necessary for weak lensing science applications using the SV data. We also discuss some of the planned algorithmic improvements that will be necessary to produce sufficiently accurate shear catalogues for the full 5-year DES, which is expected to cover 5000 square degrees.

preprint2016arXiv

The redMaPPer Galaxy Cluster Catalog From DES Science Verification Data

We describe updates to the \redmapper{} algorithm, a photometric red-sequence cluster finder specifically designed for large photometric surveys. The updated algorithm is applied to $150\,\mathrm{deg}^2$ of Science Verification (SV) data from the Dark Energy Survey (DES), and to the Sloan Digital Sky Survey (SDSS) DR8 photometric data set. The DES SV catalog is locally volume limited, and contains 786 clusters with richness $λ>20$ (roughly equivalent to $M_{\rm{500c}}\gtrsim10^{14}\,h_{70}^{-1}\,M_{\odot}$) and $0.2<z<0.9$. The DR8 catalog consists of 26311 clusters with $0.08<z<0.6$, with a sharply increasing richness threshold as a function of redshift for $z\gtrsim 0.35$. The photometric redshift performance of both catalogs is shown to be excellent, with photometric redshift uncertainties controlled at the $σ_z/(1+z)\sim 0.01$ level for $z\lesssim0.7$, rising to $\sim0.02$ at $z\sim0.9$ in DES SV. We make use of \emph{Chandra} and \emph{XMM} X-ray and South Pole Telescope Sunyaev-Zeldovich data to show that the centering performance and mass--richness scatter are consistent with expectations based on prior runs of \redmapper{} on SDSS data. We also show how the \redmapper{} \photoz{} and richness estimates are relatively insensitive to imperfect star/galaxy separation and small-scale star masks.

preprint2015arXiv

Cross-correlation of gravitational lensing from DES Science Verification data with SPT and Planck lensing

We measure the cross-correlation between weak lensing of galaxy images and of the cosmic microwave background (CMB). The effects of gravitational lensing on different sources will be correlated if the lensing is caused by the same mass fluctuations. We use galaxy shape measurements from 139 deg$^{2}$ of the Dark Energy Survey (DES) Science Verification data and overlapping CMB lensing from the South Pole Telescope (SPT) and Planck. The DES source galaxies have a median redshift of $z_{\rm med} {\sim} 0.7$, while the CMB lensing kernel is broad and peaks at $z{\sim}2$. The resulting cross-correlation is maximally sensitive to mass fluctuations at $z{\sim}0.44$. Assuming the Planck 2015 best-fit cosmology, the amplitude of the DES$\times$SPT cross-power is found to be $A = 0.88 \pm 0.30$ and that from DES$\times$Planck to be $A = 0.86 \pm 0.39$, where $A=1$ corresponds to the theoretical prediction. These are consistent with the expected signal and correspond to significances of $2.9 σ$ and $2.2 σ$ respectively. We demonstrate that our results are robust to a number of important systematic effects including the shear measurement method, estimator choice, photometric redshift uncertainty and CMB lensing systematics. Significant intrinsic alignment of galaxy shapes would increase the cross-correlation signal inferred from the data; we calculate a value of $A = 1.08 \pm 0.36$ for DES$\times$SPT when we correct the observations with a simple IA model. With three measurements of this cross-correlation now existing in the literature, there is not yet reliable evidence for any deviation from the expected LCDM level of cross-correlation, given the size of the statistical uncertainties and the significant impact of systematic errors, particularly IAs. We provide forecasts for the expected signal-to-noise of the combination of the five-year DES survey and SPT-3G.

preprint2015arXiv

Galaxy clustering, photometric redshifts and diagnosis of systematics in the DES Science Verification data

We study the clustering of galaxies detected at $i<22.5$ in the Science Verification observations of the Dark Energy Survey (DES). Two-point correlation functions are measured using $2.3\times 10^6$ galaxies over a contiguous 116 deg$^2$ region in five bins of photometric redshift width $Δz = 0.2$ in the range $0.2 < z < 1.2.$ The impact of photometric redshift errors are assessed by comparing results using a template-based photo-$z$ algorithm (BPZ) to a machine-learning algorithm (TPZ). A companion paper (Leistedt et al 2015) presents maps of several observational variables (e.g. seeing, sky brightness) which could modulate the galaxy density. Here we characterize and mitigate systematic errors on the measured clustering which arise from these observational variables, in addition to others such as Galactic dust and stellar contamination. After correcting for systematic effects we measure galaxy bias over a broad range of linear scales relative to mass clustering predicted from the Planck $Λ$CDM model, finding agreement with CFHTLS measurements with $χ^2$ of 4.0 (8.7) with 5 degrees of freedom for the TPZ (BPZ) redshifts. We test a "linear bias" model, in which the galaxy clustering is a fixed multiple of the predicted non-linear dark-matter clustering. The precision of the data allow us to determine that the linear bias model describes the observed galaxy clustering to $2.5\%$ accuracy down to scales at least $4$ to $10$ times smaller than those on which linear theory is expected to be sufficient.

preprint2015arXiv

Joint Analysis of Galaxy-Galaxy Lensing and Galaxy Clustering: Methodology and Forecasts for DES

The joint analysis of galaxy-galaxy lensing and galaxy clustering is a promising method for inferring the growth function of large scale structure. This analysis will be carried out on data from the Dark Energy Survey (DES), with its measurements of both the distribution of galaxies and the tangential shears of background galaxies induced by these foreground lenses. We develop a practical approach to modeling the assumptions and systematic effects affecting small scale lensing, which provides halo masses, and large scale galaxy clustering. Introducing parameters that characterize the halo occupation distribution (HOD), photometric redshift uncertainties, and shear measurement errors, we study how external priors on different subsets of these parameters affect our growth constraints. Degeneracies within the HOD model, as well as between the HOD and the growth function, are identified as the dominant source of complication, with other systematic effects sub-dominant. The impact of HOD parameters and their degeneracies necessitate the detailed joint modeling of the galaxy sample that we employ. We conclude that DES data will provide powerful constraints on the evolution of structure growth in the universe, conservatively/optimistically constraining the growth function to 7.9\%/4.8\% with its first-year data that covered over 1000 square degrees, and to 3.9\%/2.3\% with its full five-year data that will survey 5000 square degrees, including both statistical and systematic uncertainties.

preprint2015arXiv

Mapping and simulating systematics due to spatially-varying observing conditions in DES Science Verification data

Spatially-varying depth and characteristics of observing conditions, such as seeing, airmass, or sky background, are major sources of systematic uncertainties in modern galaxy survey analyses, in particular in deep multi-epoch surveys. We present a framework to extract and project these sources of systematics onto the sky, and apply it to the Dark Energy Survey (DES) to map the observing conditions of the Science Verification (SV) data. The resulting distributions and maps of sources of systematics are used in several analyses of DES SV to perform detailed null tests with the data, and also to incorporate systematics in survey simulations. We illustrate the complementarity of these two approaches by comparing the SV data with the BCC-UFig, a synthetic sky catalogue generated by forward-modelling of the DES SV images. We analyse the BCC-UFig simulation to construct galaxy samples mimicking those used in SV galaxy clustering studies. We show that the spatially-varying survey depth imprinted in the observed galaxy densities and the redshift distributions of the SV data are successfully reproduced by the simulation and well-captured by the maps of observing conditions. The combined use of the maps, the SV data and the BCC-UFig simulation allows us to quantify the impact of spatial systematics on $N(z)$, the redshift distributions inferred using photometric redshifts. We conclude that spatial systematics in the SV data are mainly due to seeing fluctuations and are under control in current clustering and weak lensing analyses. The framework presented here is relevant to all multi-epoch surveys, and will be essential for exploiting future surveys such as the Large Synoptic Survey Telescope (LSST), which will require detailed null-tests and realistic end-to-end image simulations to correctly interpret the deep, high-cadence observations of the sky.

preprint2015arXiv

Observation and Confirmation of Six Strong Lensing Systems in The Dark Energy Survey Science Verification Data

We report the observation and confirmation of the first group- and cluster-scale strong gravitational lensing systems found in Dark Energy Survey (DES) data. Through visual inspection of data from the Science Verification (SV) season, we identified 53 candidate systems. We then obtained spectroscopic follow-up of 21 candidates using the Gemini Multi-Object Spectrograph (GMOS) at the Gemini South telescope and the Inamori-Magellan Areal Camera and Spectrograph (IMACS) at the Magellan/Baade telescope. With this follow-up, we confirmed six candidates as gravitational lenses: Three of the systems are newly discovered, and the remaining three were previously known. Of the 21 observed candidates, the remaining 15 were either not detected in spectroscopic observations, were observed and did not exhibit continuum emission (or spectral features), or were ruled out as lensing systems. The confirmed sample consists of one group-scale and five galaxy cluster-scale lenses. The lensed sources range in redshift z ~ 0.80-3.2, and in i-band surface brightness i_{SB} ~ 23-25 mag/sq.-arcsec. (2" aperture). For each of the six systems, we estimate the Einstein radius and the enclosed mass, which have ranges ~ 5.0 - 8.6" and ~ 7.5 x 10^{12} - 6.4 x 10^{13} solar masses, respectively.

preprint2015arXiv

Photometric redshifts and clustering of emission line galaxies selected jointly by DES and eBOSS

We present the results of the first test plates of the extended Baryon Oscillation Spectroscopic Survey. This paper focuses on the emission line galaxies (ELG) population targetted from the Dark Energy Survey (DES) photometry. We analyse the success rate, efficiency, redshift distribution, and clustering properties of the targets. From the 9000 spectroscopic redshifts targetted, 4600 have been selected from the DES photometry. The total success rate for redshifts between 0.6 and 1.2 is 71\% and 68\% respectively for a bright and faint, on average more distant, samples including redshifts measured from a single strong emission line. We find a mean redshift of 0.8 and 0.87, with 15 and 13\% of unknown redshifts respectively for the bright and faint samples. In the redshift range 0.6<z<1.2, for the most secure spectroscopic redshifts, the mean redshift for the bright and faint sample is 0.85 and 0.9 respectively. Star contamination is lower than 2\%. We measure a galaxy bias averaged on scales of 1 and 10~Mpc/h of 1.72 \pm 0.1 for the bright sample and of 1.78 \pm 0.12 for the faint sample. The error on the galaxy bias have been obtained propagating the errors in the correlation function to the fitted parameters. This redshift evolution for the galaxy bias is in agreement with theoretical expectations for a galaxy population with MB-5\log h < -21.0. We note that biasing is derived from the galaxy clustering relative to a model for the mass fluctuations. We investigate the quality of the DES photometric redshifts and find that the outlier fraction can be reduced using a comparison between template fitting and neural network, or using a random forest algorithm.

preprint2015arXiv

redMaGiC: Selecting Luminous Red Galaxies from the DES Science Verification Data

We introduce redMaGiC, an automated algorithm for selecting Luminous Red Galaxies (LRGs). The algorithm was specifically developed to minimize photometric redshift uncertainties in photometric large-scale structure studies. redMaGiC achieves this by self-training the color-cuts necessary to produce a luminosity-thresholded LRG sample of constant comoving density. We demonstrate that redMaGiC photozs are very nearly as accurate as the best machine-learning based methods, yet they require minimal spectroscopic training, do not suffer from extrapolation biases, and are very nearly Gaussian. We apply our algorithm to Dark Energy Survey (DES) Science Verification (SV) data to produce a redMaGiC catalog sampling the redshift range $z\in[0.2,0.8]$. Our fiducial sample has a comoving space density of $10^{-3}\ (h^{-1} Mpc)^{-3}$, and a median photoz bias ($z_{spec}-z_{photo}$) and scatter $(σ_z/(1+z))$ of 0.005 and 0.017 respectively. The corresponding $5σ$ outlier fraction is 1.4%. We also test our algorithm with Sloan Digital Sky Survey (SDSS) Data Release 8 (DR8) and Stripe 82 data, and discuss how spectroscopic training can be used to control photoz biases at the 0.1% level.

preprint2015arXiv

Redshift distributions of galaxies in the DES Science Verification shear catalogue and implications for weak lensing

We present photometric redshift estimates for galaxies used in the weak lensing analysis of the Dark Energy Survey Science Verification (DES SV) data. Four model- or machine learning-based photometric redshift methods -- ANNZ2, BPZ calibrated against BCC-Ufig simulations, SkyNet, and TPZ -- are analysed. For training, calibration, and testing of these methods, we construct a catalogue of spectroscopically confirmed galaxies matched against DES SV data. The performance of the methods is evaluated against the matched spectroscopic catalogue, focusing on metrics relevant for weak lensing analyses, with additional validation against COSMOS photo-zs. From the galaxies in the DES SV shear catalogue, which have mean redshift $0.72\pm0.01$ over the range $0.3<z<1.3$, we construct three tomographic bins with means of $z=\{0.45, 0.67, 1.00\}$. These bins each have systematic uncertainties $δz \lesssim 0.05$ in the mean of the fiducial SkyNet photo-z $n(z)$. We propagate the errors in the redshift distributions through to their impact on cosmological parameters estimated with cosmic shear, and find that they cause shifts in the value of $σ_8$ of approx. 3%. This shift is within the one sigma statistical errors on $σ_8$ for the DES SV shear catalog. We further study the potential impact of systematic differences on the critical surface density, $Σ_{\mathrm{crit}}$, finding levels of bias safely less than the statistical power of DES SV data. We recommend a final Gaussian prior for the photo-z bias in the mean of $n(z)$ of width $0.05$ for each of the three tomographic bins, and show that this is a sufficient bias model for the corresponding cosmology analysis.

preprint2015arXiv

Weak lensing by galaxy troughs in DES Science Verification data

We measure the weak lensing shear around galaxy troughs, i.e. the radial alignment of background galaxies relative to underdensities in projections of the foreground galaxy field over a wide range of redshift in Science Verification data from the Dark Energy Survey. Our detection of the shear signal is highly significant (10 to 15$σ$ for the smallest angular scales) for troughs with the redshift range z in [0.2,0.5] of the projected galaxy field and angular diameters of 10 arcmin...1°. These measurements probe the connection between the galaxy, matter density, and convergence fields. By assuming galaxies are biased tracers of the matter density with Poissonian noise, we find agreement of our measurements with predictions in a fiducial Lambda cold dark matter model. The prediction for the lensing signal on large trough scales is virtually independent of the details of the underlying model for the connection of galaxies and matter. Our comparison of the shear around troughs with that around cylinders with large galaxy counts is consistent with a symmetry between galaxy and matter over- and underdensities. In addition, we measure the two-point angular correlation of troughs with galaxies which, in contrast to the lensing signal, is sensitive to galaxy bias on all scales. The lensing signal of troughs and their clustering with galaxies is therefore a promising probe of the statistical properties of matter underdensities and their connection to the galaxy field.

preprint2014arXiv

Combining Dark Energy Survey Science Verification Data with Near Infrared Data from the ESO VISTA Hemisphere Survey

We present the combination of optical data from the Science Verification phase of the Dark Energy Survey (DES) with near infrared data from the ESO VISTA Hemisphere Survey (VHS). The deep optical detections from DES are used to extract fluxes and associated errors from the shallower VHS data. Joint 7-band ($grizYJK$) photometric catalogues are produced in a single 3 sq-deg DECam field centred at 02h26m$-$04d36m where the availability of ancillary multi-wavelength photometry and spectroscopy allows us to test the data quality. Dual photometry increases the number of DES galaxies with measured VHS fluxes by a factor of $\sim$4.5 relative to a simple catalogue level matching and results in a $\sim$1.5 mag increase in the 80\% completeness limit of the NIR data. Almost 70\% of DES sources have useful NIR flux measurements in this initial catalogue. Photometric redshifts are estimated for a subset of galaxies with spectroscopic redshifts and initial results, although currently limited by small number statistics, indicate that the VHS data can help reduce the photometric redshift scatter at both $z<0.5$ and $z>1$. We present example DES+VHS colour selection criteria for high redshift Luminous Red Galaxies (LRGs) at $z\sim0.7$ as well as luminous quasars. Using spectroscopic observations in this field we show that the additional VHS fluxes enable a cleaner selection of both populations with $<$10\% contamination from galactic stars in the case of spectroscopically confirmed quasars and $<0.5\%$ contamination from galactic stars in the case of spectroscopically confirmed LRGs. The combined DES+VHS dataset, which will eventually cover almost 5000 sq-deg, will therefore enable a range of new science and be ideally suited for target selection for future wide-field spectroscopic surveys.

preprint2014arXiv

J-PAS: The Javalambre-Physics of the Accelerated Universe Astrophysical Survey

The Javalambre-Physics of the Accelerated Universe Astrophysical Survey (J-PAS) is a narrow band, very wide field Cosmological Survey to be carried out from the Javalambre Observatory in Spain with a purpose-built, dedicated 2.5m telescope and a 4.7 sq.deg. camera with 1.2Gpix. Starting in late 2015, J-PAS will observe 8500sq.deg. of Northern Sky and measure $0.003(1+z)$ photo-z for $9\times10^7$ LRG and ELG galaxies plus several million QSOs, sampling an effective volume of $\sim 14$ Gpc$^3$ up to $z=1.3$ and becoming the first radial BAO experiment to reach Stage IV. J-PAS will detect $7\times 10^5$ galaxy clusters and groups, setting constrains on Dark Energy which rival those obtained from its BAO measurements. Thanks to the superb characteristics of the site (seeing ~0.7 arcsec), J-PAS is expected to obtain a deep, sub-arcsec image of the Northern sky, which combined with its unique photo-z precision will produce one of the most powerful cosmological lensing surveys before the arrival of Euclid. J-PAS unprecedented spectral time domain information will enable a self-contained SN survey that, without the need for external spectroscopic follow-up, will detect, classify and measure $σ_z\sim 0.5\%$ redshifts for $\sim 4000$ SNeIa and $\sim 900$ core-collapse SNe. The key to the J-PAS potential is its innovative approach: a contiguous system of 54 filters with $145Å$ width, placed $100Å$ apart over a multi-degree FoV is a powerful "redshift machine", with the survey speed of a 4000 multiplexing low resolution spectrograph, but many times cheaper and much faster to build. The J-PAS camera is equivalent to a 4.7 sq.deg. "IFU" and it will produce a time-resolved, 3D image of the Northern Sky with a very wide range of Astrophysical applications in Galaxy Evolution, the nearby Universe and the study of resolved stellar populations.

preprint2014arXiv

Photometric redshift analysis in the Dark Energy Survey Science Verification data

We present results from a study of the photometric redshift performance of the Dark Energy Survey (DES), using the early data from a Science Verification (SV) period of observations in late 2012 and early 2013 that provided science-quality images for almost 200 sq.~deg.~at the nominal depth of the survey. We assess the photometric redshift performance using about 15000 galaxies with spectroscopic redshifts available from other surveys. These galaxies are used, in different configurations, as a calibration sample, and photo-$z$'s are obtained and studied using most of the existing photo-$z$ codes. A weighting method in a multi-dimensional color-magnitude space is applied to the spectroscopic sample in order to evaluate the photo-$z$ performance with sets that mimic the full DES photometric sample, which is on average significantly deeper than the calibration sample due to the limited depth of spectroscopic surveys. Empirical photo-$z$ methods using, for instance, Artificial Neural Networks or Random Forests, yield the best performance in the tests, achieving core photo-$z$ resolutions $σ_{68} \sim 0.08$. Moreover, the results from most of the codes, including template fitting methods, comfortably meet the DES requirements on photo-$z$ performance, therefore, providing an excellent precedent for future DES data sets.