Search

The correlation between the broad line region radius and continuum luminosity (R-L relation) of active galactic nuclei (AGNs) is critical for single-epoch mass estimates of supermassive black holes (SMBHs). At z ∼ 1-2, where AGN activity peaks, the R-L relation is constrained by the reverberation mapping (RM) lags of the Mg II line. We present 25 Mg II lags from the Australian Dark Energy Survey RM project based on 6 yr of monitoring. We define quantitative criteria to select good lag measurements and verify their reliability with simulations based on both the damped random walk stochastic model and the rescaled, resampled versions of the observed light curves of local, well-measured AGN. Our sample significantly increases the number of Mg II lags and extends the R-L relation to higher redshifts and luminosities. The relative iron line strength RFe has little impact on the R-L relation. The best-fitting Mg II R-L relation has a slope α = 0.39 ± 0.08 with an intrinsic scatter σrl = 0.15+−000203. The slope is consistent with previous measurements and shallower than the H β R-L relation. The intrinsic scatter of the new R-L relation is substantially smaller than previous studies and comparable to the intrinsic scatter of the H β R-L relation. Our new R-L relation will enable more precise single-epoch mass estimates and SMBH demographic studies at cosmic noon.

We present the first Hubble diagram of superluminous supernovae (SLSNe) out to a redshift of two, together with constraints on the matter density, ωM, and the dark energy equation-of-state parameter, w(p/ρ). We build a sample of 20 cosmologically useful SLSNe I based on light curve and spectroscopy quality cuts. We confirm the robustness of the peak-decline SLSN I standardization relation with a larger data set and improved fitting techniques than previous works. We then solve the SLSN model based on the above standardization via minimization of the χ2 computed from a covariance matrix that includes statistical and systematic uncertainties. For a spatially flat Λ cold dark matter (ΛCDM) cosmological model, we find $\Omega{\rm M}=0.38^{+0.24}_{-0.19}$, with an rms of 0.27 mag for the residuals of the distance moduli. For a w0waCDM cosmological model, the addition of SLSNe I to a 'baseline' measurement consisting of Planck temperature together with Type Ia supernovae, results in a small improvement in the constraints of w0 and wa of 4 per cent. We present simulations of future surveys with 868 and 492 SLSNe I (depending on the configuration used) and show that such a sample can deliver cosmological constraints in a flat ΛCDM model with the same precision (considering only statistical uncertainties) as current surveys that use Type Ia supernovae, while providing a factor of 2-3 improvement in the precision of the constraints on the time variation of dark energy, w0 and wa. This paper represents the proof of concept for superluminous supernova cosmology, and demonstrates they can provide an independent test of cosmology in the high-redshift (z > 1) universe.

We present the first Hubble diagram of superluminous supernovae (SLSNe) out to a redshift of two, together with constraints on the matter density, ωM, and the dark energy equation-of-state parameter, w(p/ρ). We build a sample of 20 cosmologically useful SLSNe I based on light curve and spectroscopy quality cuts. We confirm the robustness of the peak-decline SLSN I standardization relation with a larger data set and improved fitting techniques than previous works. We then solve the SLSN model based on the above standardization via minimization of the χ2 computed from a covariance matrix that includes statistical and systematic uncertainties. For a spatially flat Λ cold dark matter (ΛCDM) cosmological model, we find $\Omega{\rm M}=0.38^{+0.24}_{-0.19}$, with an rms of 0.27 mag for the residuals of the distance moduli. For a w0waCDM cosmological model, the addition of SLSNe I to a 'baseline' measurement consisting of Planck temperature together with Type Ia supernovae, results in a small improvement in the constraints of w0 and wa of 4 per cent. We present simulations of future surveys with 868 and 492 SLSNe I (depending on the configuration used) and show that such a sample can deliver cosmological constraints in a flat ΛCDM model with the same precision (considering only statistical uncertainties) as current surveys that use Type Ia supernovae, while providing a factor of 2-3 improvement in the precision of the constraints on the time variation of dark energy, w0 and wa. This paper represents the proof of concept for superluminous supernova cosmology, and demonstrates they can provide an independent test of cosmology in the high-redshift (z > 1) universe.

Periodically variable quasars have been suggested as close binary supermassive black holes. We present a systematic search for periodic light curves in 625 spectroscopically confirmed quasars with a median redshift of 1.8 in a 4.6 deg2 overlapping region of the Dark Energy Survey Supernova (DES-SN) fields and the Sloan Digital Sky Survey Stripe 82 (SDSS-S82). Our sample has a unique 20-yr long multicolour (griz) light curve enabled by combining DES-SN Y6 observations with archival SDSS-S82 data. The deep imaging allows us to search for periodic light curves in less luminous quasars (down to r∼23.5 mag) powered by less massive black holes (with masses ≳ 108.5M⊙) at high redshift for the first time. We find five candidates with significant (at >99.74 per cent single-frequency significance in at least two bands with a global p-value of ∼7 × 10-4 -3 × 10-3 accounting for the look-elsewhere effect) periodicity with observed periods of ∼3-5 yr (i.e. 1-2 yr in rest frame) having ∼4-6 cycles spanned by the observations. If all five candidates are periodically variable quasars, this translates into a detection rate of ∼0.8+0.5-0.3 per cent or ∼1.1+0.7-0.5 quasar per deg2. Our detection rate is 4-80 times larger than those found by previous searches using shallower surveys over larger areas. This discrepancy is likely caused by differences in the quasar populations probed and the survey data qualities. We discuss implications on the future direct detection of low-frequency gravitational waves. Continued photometric monitoring will further assess the robustness and characteristics of these candidate periodic quasars to determine their physical origins.

The Dark Energy Spectroscopic Instrument Legacy Surveys, a combination of three ground-based imaging surveys, have mapped 16,000 deg2 in three optical bands (g, r, and z) to a depth 1-2 mag deeper than the Sloan Digital Sky Survey. Our work addresses one of the major challenges of wide-field imaging surveys conducted at ground-based observatories: the varying depth that results from varying observing conditions at Earth-bound sites. To mitigate these effects, the Legacy Surveys (the Dark Energy Camera Legacy Survey, or DECaLS; the Mayall z-band Legacy Survey, or MzLS; and the Beiijing-Arizona Sky Survey, or BASS) employed a unique strategy to dynamically adjust the exposure times as rapidly as possible in response to the changing observing conditions. We present the tiling and observing strategies used by the first two of these surveys. We demonstrate that the tiling and dynamic observing strategies jointly result in a more uniform-depth survey that has higher efficiency for a given total observing time compared with the traditional approach of using fixed exposure times.

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 of1500 km s-1, and the light curve shows an initial peak that fades after 50 d before slowly rebrightening over a further 100 d to reach an absolute brightness of Mr ∼-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 56Co, 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 56Ni production of 0.075 M· from an ∼15 M· red supergiant progenitor typical of other SNe II, or (ii) a relatively compact ∼40 M· star, explosion energy of 1 foe, and 0.08 M· of 56Ni. Both scenarios require additional energy input to explain the late-time light curve, which is consistent with fallback accretion at a rate of ∼0.5 × 10-8 M· s-1.

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 of1500 km s-1, and the light curve shows an initial peak that fades after 50 d before slowly rebrightening over a further 100 d to reach an absolute brightness of Mr ∼-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 56Co, 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 56Ni production of 0.075 M· from an ∼15 M· red supergiant progenitor typical of other SNe II, or (ii) a relatively compact ∼40 M· star, explosion energy of 1 foe, and 0.08 M· of 56Ni. Both scenarios require additional energy input to explain the late-time light curve, which is consistent with fallback accretion at a rate of ∼0.5 × 10-8 M· s-1.

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. For the DES-SN sample, when considering a 5D (z, x1, c, α, β) bias correction, we find evidence of a Hubble residual 'mass step', where SNe Ia in high-mass galaxies (>1010M☉) are intrinsically more luminous (after correction) than their low-mass counterparts by γ = 0.040 ± 0.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 ± 0.020 mag. The 1D−5D γ difference for DES-SN is 0.026 ± 0.009 mag. We show that this difference is due to a strong correlation between host galaxy stellar mass and the x1 component of the 5D distance-bias correction. Including an intrinsic correlation between the observed properties of SNe Ia, stretch and colour, 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 modelling correlations between galaxy properties and SN is necessary to ensure unbiased precision estimates of the dark energy equation of state as we enter the era of LSST.

We report the results of the STRong lensing Insights into the Dark Energy Survey (STRIDES) follow-up campaign of the late 2017/early 2018 season. We obtained spectra of 65 lensed quasar candidates with ESO Faint Object Spectrograph and Camera 2 on the NTT and Echellette Spectrograph and Imager onKeck, confirming 10 newlensed quasars and 10 quasar pairs. Eight lensed quasars are doubly imaged with source redshifts between 0.99 and 2.90, one is triply imaged (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 Southern Astrophysical Research Telescope or Near InfraRed Camera 2 on Keck, give total magnifications between 3.2 and 5.6, and Einstein radii between 0.49 and 1.97 arcsec. After spectroscopic follow-up, we extract multi-epoch grizY photometry of confirmed lensed quasars and contaminant quasar+star pairs from DES data using parametric multiband modelling, and compare variability in each system's components. By measuring the reduced χ2 associated with fitting all epochs to the samemagnitude, we find a simple cut on the less variable component that retains all confirmed lensed quasars, while removing 94 per cent of contaminant systems. Based on our spectroscopic follow-up, this variability information improves selection of lensed quasars and quasar pairs from 34-45 per cent to 51-70 per cent, with most remaining contaminants being star-forming galaxies. Using mock lensed quasar light curves 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, explained by the magnification bias and hence lower luminosity/more variable sources in quads.

We report the results of the STRong lensing Insights into the Dark Energy Survey (STRIDES) follow-up campaign of the late 2017/early 2018 season. We obtained spectra of 65 lensed quasar candidates with ESO Faint Object Spectrograph and Camera 2 on the NTT and Echellette Spectrograph and Imager onKeck, confirming 10 newlensed quasars and 10 quasar pairs. Eight lensed quasars are doubly imaged with source redshifts between 0.99 and 2.90, one is triply imaged (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 Southern Astrophysical Research Telescope or Near InfraRed Camera 2 on Keck, give total magnifications between 3.2 and 5.6, and Einstein radii between 0.49 and 1.97 arcsec. After spectroscopic follow-up, we extract multi-epoch grizY photometry of confirmed lensed quasars and contaminant quasar+star pairs from DES data using parametric multiband modelling, and compare variability in each system's components. By measuring the reduced χ2 associated with fitting all epochs to the samemagnitude, we find a simple cut on the less variable component that retains all confirmed lensed quasars, while removing 94 per cent of contaminant systems. Based on our spectroscopic follow-up, this variability information improves selection of lensed quasars and quasar pairs from 34-45 per cent to 51-70 per cent, with most remaining contaminants being star-forming galaxies. Using mock lensed quasar light curves 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, explained by the magnification bias and hence lower luminosity/more variable sources in quads.

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. For the DES-SN sample, when considering a 5D (z, x1, c, α, β) bias correction, we find evidence of a Hubble residual 'mass step', where SNe Ia in high-mass galaxies (>1010M☉) are intrinsically more luminous (after correction) than their low-mass counterparts by γ = 0.040 ± 0.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 ± 0.020 mag. The 1D−5D γ difference for DES-SN is 0.026 ± 0.009 mag. We show that this difference is due to a strong correlation between host galaxy stellar mass and the x1 component of the 5D distance-bias correction. Including an intrinsic correlation between the observed properties of SNe Ia, stretch and colour, 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 modelling correlations between galaxy properties and SN is necessary to ensure unbiased precision estimates of the dark energy equation of state as we enter the era of LSST.

Calibration uncertainties have been the leading systematic uncertainty in recent analyses using Type Ia supernovae (SNe Ia) to measure cosmological parameters. To improve the calibration, we present the application of spectral energy distribution-dependent 'chromatic corrections' to the SN light-curve photometry from the Dark Energy Survey (DES). These corrections depend on the combined atmospheric and instrumental transmission function for each exposure, and they affect photometry at the 0.01 mag (1 per cent) level, comparable to systematic uncertainties in calibration and photometry. Fitting our combined DES and low-z SN Ia sample with baryon acoustic oscillation (BAO) and cosmic microwave background (CMB) priors for the cosmological parameters m (the fraction of the critical density of the universe comprised of matter) and w (the dark energy equation of state parameter), we compare those parameters before and after applying the corrections. We find the change in w and m due to not including chromatic corrections is −0.002 and 0.000, respectively, for the DES-SN3YR sample with BAO and CMB priors, consistent with a larger DES-SN3YR-like simulation, which has a w-change of 0.0005 with an uncertainty of 0.008 and an m change of 0.000 with an uncertainty of 0.002. However, when considering samples on individual CCDs we find large redshift-dependent biases (∼0.02 in distance modulus) for SN distances.

We present a multi-messenger measurement of the Hubble constant H 0 using the binary-black-hole merger GW170814 as a standard siren, combined with a photometric redshift catalog from the Dark Energy Survey (DES). The luminosity distance is obtained from the gravitational wave signal detected by the Laser Interferometer Gravitational-Wave Observatory (LIGO)/Virgo Collaboration (LVC) on 2017 August 14, and the redshift information is provided by the DES Year 3 data. Black hole mergers such as GW170814 are expected to lack bright electromagnetic emission to uniquely identify their host galaxies and build an object-by-object Hubble diagram. However, they are suitable for a statistical measurement, provided that a galaxy catalog of adequate depth and redshift completion is available. Here we present the first Hubble parameter measurement using a black hole merger. Our analysis results in , which is consistent with both SN Ia and cosmic microwave background measurements of the Hubble constant. The quoted 68% credible region comprises 60% of the uniform prior range [20, 140] km s-1 Mpc-1, and it depends on the assumed prior range. If we take a broader prior of [10, 220] km s-1 Mpc-1, we find (57% of the prior range). Although a weak constraint on the Hubble constant from a single event is expected using the dark siren method, a multifold increase in the LVC event rate is anticipated in the coming years and combinations of many sirens will lead to improved constraints on H 0.

The DESI Legacy Imaging Surveys (http://legacysurvey.org/) are a combination of three public projects (the Dark Energy Camera Legacy Survey, the Beijing-Arizona Sky Survey, and the Mayall z-band Legacy Survey) that will jointly image ≈14,000 deg2 of the extragalactic sky visible from the northern hemisphere in three optical bands (g, r, and z) using telescopes at the Kitt Peak National Observatory and the Cerro Tololo Inter-American Observatory. The combined survey footprint is split into two contiguous areas by the Galactic plane. The optical imaging is conducted using a unique strategy of dynamically adjusting the exposure times and pointing selection during observing that results in a survey of nearly uniform depth. In addition to calibrated images, the project is delivering a catalog, constructed by using a probabilistic inference-based approach to estimate source shapes and brightnesses. The catalog includes photometry from the grz optical bands and from four mid-infrared bands (at 3.4, 4.6, 12, and 22 μm) observed by the Wide-field Infrared Survey Explorer satellite during its full operational lifetime. The project plans two public data releases each year. All the software used to generate the catalogs is also released with the data. This paper provides an overview of the Legacy Surveys project.

We present a multi-messenger measurement of the Hubble constant H 0 using the binary-black-hole merger GW170814 as a standard siren, combined with a photometric redshift catalog from the Dark Energy Survey (DES). The luminosity distance is obtained from the gravitational wave signal detected by the Laser Interferometer Gravitational-Wave Observatory (LIGO)/Virgo Collaboration (LVC) on 2017 August 14, and the redshift information is provided by the DES Year 3 data. Black hole mergers such as GW170814 are expected to lack bright electromagnetic emission to uniquely identify their host galaxies and build an object-by-object Hubble diagram. However, they are suitable for a statistical measurement, provided that a galaxy catalog of adequate depth and redshift completion is available. Here we present the first Hubble parameter measurement using a black hole merger. Our analysis results in , which is consistent with both SN Ia and cosmic microwave background measurements of the Hubble constant. The quoted 68% credible region comprises 60% of the uniform prior range [20, 140] km s-1 Mpc-1, and it depends on the assumed prior range. If we take a broader prior of [10, 220] km s-1 Mpc-1, we find (57% of the prior range). Although a weak constraint on the Hubble constant from a single event is expected using the dark siren method, a multifold increase in the LVC event rate is anticipated in the coming years and combinations of many sirens will lead to improved constraints on H 0.

We search Dark Energy Survey (DES) Year 3 imaging data for galaxy-galaxy strong gravitational lenses using convolutional neural networks. We generate 250 000 simulated lenses at redshifts > 0.8 from which we create a data set for training the neural networks with realistic seeing, sky and shot noise. Using the simulations as a guide, we build a catalogue of 1.1 million DES sources with 1.8 < g − i < 5, 0.6 < g − r < 3, r mag > 19, g mag > 20, and i mag > 18.2. We train two ensembles of neural networks on training sets consisting of simulated lenses, simulated non-lenses, and real sources. We use the neural networks to score images of each of the sources in our catalogue with a value from 0 to 1, and select those with scores greater than a chosen threshold for visual inspection, resulting in a candidate set of 7301 galaxies. During visual inspection, we rate 84 as 'probably' or 'definitely' lenses. Four of these are previously known lenses or lens candidates. We inspect a further 9428 candidates with a different score threshold, and identify four new candidates. We present 84 new strong lens candidates, selected after a few hours of visual inspection by astronomers. This catalogue contains a comparable number of high-redshift lenses to that predicted by simulations. Based on simulations, we estimate our sample to contain most discoverable lenses in this imaging and at this redshift range.

We search Dark Energy Survey (DES) Year 3 imaging data for galaxy-galaxy strong gravitational lenses using convolutional neural networks. We generate 250 000 simulated lenses at redshifts > 0.8 from which we create a data set for training the neural networks with realistic seeing, sky and shot noise. Using the simulations as a guide, we build a catalogue of 1.1 million DES sources with 1.8 < g − i < 5, 0.6 < g − r < 3, r mag > 19, g mag > 20, and i mag > 18.2. We train two ensembles of neural networks on training sets consisting of simulated lenses, simulated non-lenses, and real sources. We use the neural networks to score images of each of the sources in our catalogue with a value from 0 to 1, and select those with scores greater than a chosen threshold for visual inspection, resulting in a candidate set of 7301 galaxies. During visual inspection, we rate 84 as 'probably' or 'definitely' lenses. Four of these are previously known lenses or lens candidates. We inspect a further 9428 candidates with a different score threshold, and identify four new candidates. We present 84 new strong lens candidates, selected after a few hours of visual inspection by astronomers. This catalogue contains a comparable number of high-redshift lenses to that predicted by simulations. Based on simulations, we estimate our sample to contain most discoverable lenses in this imaging and at this redshift range.

We present griz light curves of 251 SNe Ia from the first 3 years of the Dark Energy Survey Supernova Program's (DES-SN) spectroscopically classified sample. The photometric pipeline described in this paper produces the calibrated fluxes and associated uncertainties used in the cosmological parameter analysis by employing a scene modeling approach that simultaneously models a variable transient flux and temporally constant host galaxy. We inject artificial point sources onto DECam images to test the accuracy of our photometric method. Upon comparison of input and measured artificial supernova fluxes, we find that flux biases peak at 3 mmag. We require corrections to our photometric uncertainties as a function of host galaxy surface brightness at the transient location, similar to that seen by the DES Difference Imaging Pipeline used to discover transients. The public release of the light curves can be found at https://des.ncsa.illinois.edu/releases/sn.

We present griz light curves of 251 SNe Ia from the first 3 years of the Dark Energy Survey Supernova Program's (DES-SN) spectroscopically classified sample. The photometric pipeline described in this paper produces the calibrated fluxes and associated uncertainties used in the cosmological parameter analysis by employing a scene modeling approach that simultaneously models a variable transient flux and temporally constant host galaxy. We inject artificial point sources onto DECam images to test the accuracy of our photometric method. Upon comparison of input and measured artificial supernova fluxes, we find that flux biases peak at 3 mmag. We require corrections to our photometric uncertainties as a function of host galaxy surface brightness at the transient location, similar to that seen by the DES Difference Imaging Pipeline used to discover transients. The public release of the light curves can be found at https://des.ncsa.illinois.edu/releases/sn.

Calibration uncertainties have been the leading systematic uncertainty in recent analyses using Type Ia supernovae (SNe Ia) to measure cosmological parameters. To improve the calibration, we present the application of spectral energy distribution-dependent 'chromatic corrections' to the SN light-curve photometry from the Dark Energy Survey (DES). These corrections depend on the combined atmospheric and instrumental transmission function for each exposure, and they affect photometry at the 0.01 mag (1 per cent) level, comparable to systematic uncertainties in calibration and photometry. Fitting our combined DES and low-z SN Ia sample with baryon acoustic oscillation (BAO) and cosmic microwave background (CMB) priors for the cosmological parameters m (the fraction of the critical density of the universe comprised of matter) and w (the dark energy equation of state parameter), we compare those parameters before and after applying the corrections. We find the change in w and m due to not including chromatic corrections is −0.002 and 0.000, respectively, for the DES-SN3YR sample with BAO and CMB priors, consistent with a larger DES-SN3YR-like simulation, which has a w-change of 0.0005 with an uncertainty of 0.008 and an m change of 0.000 with an uncertainty of 0.002. However, when considering samples on individual CCDs we find large redshift-dependent biases (∼0.02 in distance modulus) for SN distances.

The primary goals of the STRong lensing Insights into the Dark Energy Survey (STRIDES) collaboration are to measure the dark energy equation of state parameter and the free streaming length of dark matter. To this aim, STRIDES is discovering strongly lensed quasars in the imaging data of the Dark Energy Survey and following them up to measure time delays, high resolution imaging, and spectroscopy sufficient to construct accurate lens models. In this paper, we first present forecasts for STRIDES. Then, we describe the STRIDES classification scheme, and give an overview of the Fall 2016 follow-up campaign. We continue by detailing the results of two selection methods, the outlier selection technique and a morphological algorithm, and presenting lens models of a system that could possibly be a lensed quasar in an unusual configuration. We conclude with the summary statistics of the Fall 2016 campaign. Including searches presented in companion papers (Anguita et al.; Ostrovski et al.), STRIDES followed up 117 targets identifying 7 new strongly lensed systems, and 7 nearly identical quasars, which could be confirmed as lenses by the detection of the lens galaxy. 76 candidates were rejected and 27 remain otherwise inconclusive, for a success rate in the range of 6-35 per cent. This rate is comparable to that of previous searches like SDSS Quasar Lens Search even though the parent data set of STRIDES is purely photometric and our selection of candidates cannot rely on spectroscopic information.

The primary goals of the STRong lensing Insights into the Dark Energy Survey (STRIDES) collaboration are to measure the dark energy equation of state parameter and the free streaming length of dark matter. To this aim, STRIDES is discovering strongly lensed quasars in the imaging data of the Dark Energy Survey and following them up to measure time delays, high resolution imaging, and spectroscopy sufficient to construct accurate lens models. In this paper, we first present forecasts for STRIDES. Then, we describe the STRIDES classification scheme, and give an overview of the Fall 2016 follow-up campaign. We continue by detailing the results of two selection methods, the outlier selection technique and a morphological algorithm, and presenting lens models of a system that could possibly be a lensed quasar in an unusual configuration. We conclude with the summary statistics of the Fall 2016 campaign. Including searches presented in companion papers (Anguita et al.; Ostrovski et al.), STRIDES followed up 117 targets identifying 7 new strongly lensed systems, and 7 nearly identical quasars, which could be confirmed as lenses by the detection of the lens galaxy. 76 candidates were rejected and 27 remain otherwise inconclusive, for a success rate in the range of 6-35 per cent. This rate is comparable to that of previous searches like SDSS Quasar Lens Search even though the parent data set of STRIDES is purely photometric and our selection of candidates cannot rely on spectroscopic information.

We combine Dark Energy Survey Year 1 clustering and weak lensing data with baryon acoustic oscillations and Big Bang nucleosynthesis experiments to constrain the Hubble constant. Assuming a flat ΛCDM model with minimal neutrino mass (Σmυ = 0.06 eV), we find H0 = 67.4 -1.2+1.1 km s-1 Mpc-1 (68 per cent CL). This result is completely independent of Hubble constant measurements based on the distance ladder, cosmic microwave background anisotropies (both temperature and polarization), and strong lensing constraints. There are now five data sets that: (a) have no shared observational systematics; and (b) each constrains the Hubble constant with fractional uncertainty at the few-per cent level. We compare these five independent estimates, and find that, as a set, the differences between them are significant at the 2.5σ level (χ2/dof = 24/11, probability to exceed = 1.1 per cent). Having set the threshold for consistency at 3σ, we combine all five data sets to arrive at H0 = 69.3-0.6+0.4 km s-1 Mpc-1.

We present observations of DES16C2nm, the first spectroscopically confirmed hydrogen-free superluminous supernova (SLSN-I) at redshift z ≊ 2. DES16C2nm was discovered by the Dark Energy Survey (DES) Supernova Program, with follow-up photometric data from the Hubble Space Telescope, Gemini, and the European Southern Observatory Very Large Telescope supplementing the DES data. Spectroscopic observations confirm DES16C2nm to be at z = 1.998, and spectroscopically similar to Gaia16apd (a SLSN-I at z = 0.102), with a peak absolute magnitude of U = -22.26 ± 0.06. The high redshift of DES16C2nm provides a unique opportunity to study the ultraviolet (UV) properties of SLSNe-I. Combining DES16C2nm with 10 similar events from the literature, we show that there exists a homogeneous class of SLSNe-I in the UV (λrest ≊ 2500 Å), with peak luminosities in the (rest-frame) U band, and increasing absorption to shorter wavelengths. There is no evidence that the mean photometric and spectroscopic properties of SLSNe-I differ between low (z < 1) and high redshift (z > 1), but there is clear evidence of diversity in the spectrum at λrest < 2000 Å, possibly caused by the variations in temperature between events. No significant correlations are observed between spectral line velocities and photometric luminosity. Using these data, we estimate that SLSNe-I can be discovered to z = 3.8 by DES. While SLSNe-I are typically identified from their blue observed colors at low redshift (z < 1), we highlight that at z > 2 these events appear optically red, peaking in the observer-frame z-band. Such characteristics are critical to identify these objects with future facilities such as the Large Synoptic Survey Telescope, Euclid, and the Wide-field Infrared Survey Telescope, which should detect such SLSNe-I to z = 3.5, 3.7, and 6.6, respectively.

We present observations of DES16C2nm, the first spectroscopically confirmed hydrogen-free superluminous supernova (SLSN-I) at redshift z ≊ 2. DES16C2nm was discovered by the Dark Energy Survey (DES) Supernova Program, with follow-up photometric data from the Hubble Space Telescope, Gemini, and the European Southern Observatory Very Large Telescope supplementing the DES data. Spectroscopic observations confirm DES16C2nm to be at z = 1.998, and spectroscopically similar to Gaia16apd (a SLSN-I at z = 0.102), with a peak absolute magnitude of U = -22.26 ± 0.06. The high redshift of DES16C2nm provides a unique opportunity to study the ultraviolet (UV) properties of SLSNe-I. Combining DES16C2nm with 10 similar events from the literature, we show that there exists a homogeneous class of SLSNe-I in the UV (λrest ≊ 2500 Å), with peak luminosities in the (rest-frame) U band, and increasing absorption to shorter wavelengths. There is no evidence that the mean photometric and spectroscopic properties of SLSNe-I differ between low (z < 1) and high redshift (z > 1), but there is clear evidence of diversity in the spectrum at λrest < 2000 Å, possibly caused by the variations in temperature between events. No significant correlations are observed between spectral line velocities and photometric luminosity. Using these data, we estimate that SLSNe-I can be discovered to z = 3.8 by DES. While SLSNe-I are typically identified from their blue observed colors at low redshift (z < 1), we highlight that at z > 2 these events appear optically red, peaking in the observer-frame z-band. Such characteristics are critical to identify these objects with future facilities such as the Large Synoptic Survey Telescope, Euclid, and the Wide-field Infrared Survey Telescope, which should detect such SLSNe-I to z = 3.5, 3.7, and 6.6, respectively.

The discovery of a kilonova (KN) associated with the Advanced LIGO (aLIGO)/Virgo event GW170817 opens up new avenues of multi-messenger astrophysics. Here, using realistic simulations, we provide estimates of the number of KNe that could be found in data from past, present, and future surveys without a gravitational-wave trigger. For the simulation, we construct a spectral time-series model based on the DES-GW multi-band light curve from the single known KN event, and we use an average of BNS rates from past studies of 103Gpc-3 yr-1, consistent with the one event found so far. Examining past and current data sets from transient surveys, the number of KNe we expect to find for ASAS-SN, SDSS, PS1, SNLS, DES, and SMT is between 0 and 0.3. We predict the number of detections per future survey to be 8.3 from ATLAS, 10.6 from ZTF, 5.5/69 from LSST (the Deep Drilling/Wide Fast Deep), and 16.0 from WFIRST. The maximum redshift of KNe discovered for each survey is z=0.8 for WFIRST, z=0.8 for LSST, and z=0.8 for ZTF and ATLAS. This maximum redshift for WFIRST is well beyond the sensitivity of aLIGO and some future GW missions. For the LSST survey, we also provide contamination estimates from Type Ia and core-collapse supernovae: after light curve and template-matching requirements, we estimate a background of just two events. More broadly, we stress that future transient surveys should consider how to optimize their search strategies to improve their detection efficiency and to consider similar analyses for GW follow-up programs.

We present time-delay measurements for the new quadruple imaged quasar DES J04085354, the first quadruple imaged quasar found in the Dark Energy Survey (DES). Our result is made possible by implementing a new observational strategy using almost daily observations with the MPIA 2.2 m telescope at La Silla observatory and deep exposures reaching a signal-to-noise ratio of about 1000 per quasar image. This data qualityallows us to catch small photometric variations (a few mmag rms) of the quasar, acting on temporal scales much shorter than microlensing, and hence making the time delay measurement very robust against microlensing. In only seven months we very accurately measured one of the time delays in DES J04085354: - "t(AB) = 112.1 ± 2.1 days (1.8%) using only the MPIA 2.2 m data. In combination with data taken with the 1.2 m Euler Swiss telescope, we also measured two delays involving the D component of the system - "t(AD) = ± 12.8 days (8.2%) and - "t(BD) = 42.4 ± 17.6 days (41%), where all the error bars include systematics. Turning these time delays into cosmological constraints will require deep Hubble Space Telescope (HST) imaging or ground-based adaptive optics (AO), and information on the velocity field of the lensing galaxy.

We combine Dark Energy Survey Year 1 clustering and weak lensing data with baryon acoustic oscillations and Big Bang nucleosynthesis experiments to constrain the Hubble constant. Assuming a flat ΛCDM model with minimal neutrino mass (Σmυ = 0.06 eV), we find H0 = 67.4 -1.2+1.1 km s-1 Mpc-1 (68 per cent CL). This result is completely independent of Hubble constant measurements based on the distance ladder, cosmic microwave background anisotropies (both temperature and polarization), and strong lensing constraints. There are now five data sets that: (a) have no shared observational systematics; and (b) each constrains the Hubble constant with fractional uncertainty at the few-per cent level. We compare these five independent estimates, and find that, as a set, the differences between them are significant at the 2.5σ level (χ2/dof = 24/11, probability to exceed = 1.1 per cent). Having set the threshold for consistency at 3σ, we combine all five data sets to arrive at H0 = 69.3-0.6+0.4 km s-1 Mpc-1.

We present time-delay measurements for the new quadruple imaged quasar DES J04085354, the first quadruple imaged quasar found in the Dark Energy Survey (DES). Our result is made possible by implementing a new observational strategy using almost daily observations with the MPIA 2.2 m telescope at La Silla observatory and deep exposures reaching a signal-to-noise ratio of about 1000 per quasar image. This data qualityallows us to catch small photometric variations (a few mmag rms) of the quasar, acting on temporal scales much shorter than microlensing, and hence making the time delay measurement very robust against microlensing. In only seven months we very accurately measured one of the time delays in DES J04085354: - "t(AB) = 112.1 ± 2.1 days (1.8%) using only the MPIA 2.2 m data. In combination with data taken with the 1.2 m Euler Swiss telescope, we also measured two delays involving the D component of the system - "t(AD) = ± 12.8 days (8.2%) and - "t(BD) = 42.4 ± 17.6 days (41%), where all the error bars include systematics. Turning these time delays into cosmological constraints will require deep Hubble Space Telescope (HST) imaging or ground-based adaptive optics (AO), and information on the velocity field of the lensing galaxy.

The discovery of a kilonova (KN) associated with the Advanced LIGO (aLIGO)/Virgo event GW170817 opens up new avenues of multi-messenger astrophysics. Here, using realistic simulations, we provide estimates of the number of KNe that could be found in data from past, present, and future surveys without a gravitational-wave trigger. For the simulation, we construct a spectral time-series model based on the DES-GW multi-band light curve from the single known KN event, and we use an average of BNS rates from past studies of 103Gpc-3 yr-1, consistent with the one event found so far. Examining past and current data sets from transient surveys, the number of KNe we expect to find for ASAS-SN, SDSS, PS1, SNLS, DES, and SMT is between 0 and 0.3. We predict the number of detections per future survey to be 8.3 from ATLAS, 10.6 from ZTF, 5.5/69 from LSST (the Deep Drilling/Wide Fast Deep), and 16.0 from WFIRST. The maximum redshift of KNe discovered for each survey is z=0.8 for WFIRST, z=0.8 for LSST, and z=0.8 for ZTF and ATLAS. This maximum redshift for WFIRST is well beyond the sensitivity of aLIGO and some future GW missions. For the LSST survey, we also provide contamination estimates from Type Ia and core-collapse supernovae: after light curve and template-matching requirements, we estimate a background of just two events. More broadly, we stress that future transient surveys should consider how to optimize their search strategies to improve their detection efficiency and to consider similar analyses for GW follow-up programs.

We present results for the first three years of OzDES, a six year programme to obtain redshifts for objects in the Dark Energy Survey (DES) supernova fields using the 2dF fibre positioner and AAOmega spectrograph on the Anglo-Australian Telescope. OzDES is a multi-object spectroscopic survey targeting multiple types of targets at multiple epochs over a multiyear baseline and is one of the first multi-object spectroscopic surveys to dynamically include transients into the target list soon after their discovery. At the end of three years, OzDES has spectroscopically confirmed almost 100 supernovae, and has measured redshifts for 17 000 objects, including the redshifts of 2566 supernova hosts. We examine how our ability to measure redshifts for targets of various types depends on signal-to-noise ratio (S/N), magnitude and exposure time, finding that our redshift success rate increases significantly at a S/N of 2-3 per 1-Å bin. We also find that the change in S/N with exposure time closely matches the Poisson limit for stacked exposures as long as 10 h. We use these results to predict the redshift yield of the full OzDES survey, as well as the potential yields of future surveys on other facilities such as (i.e. the 4-m Multi-Object Spectroscopic Telescope, the Subaru Prime Focus Spectrograph and the Maunakea Spectroscopic Explorer). This work marks the first OzDES data release, comprising 14 693 redshifts. OzDES is on target to obtain over 30 000 redshifts over the 6-yr duration of the survey, including a yield of approximately 5700 supernova host-galaxy redshifts.

We present a study of NGC 4993, the host galaxy of the GW170817 gravitational-wave event, the GRB 170817A short gamma-ray burst (sGRB), and the AT 2017gfo kilonova. We use Dark Energy Camera imaging, AAT spectra, and publicly available data, relating our findings to binary neutron star (BNS) formation scenarios and merger delay timescales. NGC 4993 is a nearby early-type galaxy, with an i-band Sérsic index n = 4.0 and low asymmetry (A = 0.04 ±0.01). These properties are unusual for sGRB hosts. However, NGC 4993 presents shell-like structures and dust lanes indicative of a recent galaxy merger, with the optical transient located close to a shell. We constrain the star formation history (SFH) of the galaxy assuming that the galaxy merger produced a star formation burst, but find little to no ongoing star formation in either spatially resolved broadband SED or spectral fitting. We use the best-fit SFH to estimate the BNS merger rate in this type of galaxy, as . If star formation is the only considered BNS formation scenario, the expected number of BNS mergers from early-type galaxies detectable with LIGO during its first two observing seasons is , as opposed to ∼0.5 from all galaxy types. Hypothesizing that the binary formed due to dynamical interactions during the galaxy merger, the subsequent time elapsed can constrain the delay time of the BNS coalescence. By using velocity dispersion estimates and the position of the shells, we find that the galaxy merger occurred t mer ≲ 200 Myr prior to the BNS coalescence.

We present a study of NGC 4993, the host galaxy of the GW170817 gravitational-wave event, the GRB 170817A short gamma-ray burst (sGRB), and the AT 2017gfo kilonova. We use Dark Energy Camera imaging, AAT spectra, and publicly available data, relating our findings to binary neutron star (BNS) formation scenarios and merger delay timescales. NGC 4993 is a nearby early-type galaxy, with an i-band Sérsic index n = 4.0 and low asymmetry (A = 0.04 ±0.01). These properties are unusual for sGRB hosts. However, NGC 4993 presents shell-like structures and dust lanes indicative of a recent galaxy merger, with the optical transient located close to a shell. We constrain the star formation history (SFH) of the galaxy assuming that the galaxy merger produced a star formation burst, but find little to no ongoing star formation in either spatially resolved broadband SED or spectral fitting. We use the best-fit SFH to estimate the BNS merger rate in this type of galaxy, as . If star formation is the only considered BNS formation scenario, the expected number of BNS mergers from early-type galaxies detectable with LIGO during its first two observing seasons is , as opposed to ∼0.5 from all galaxy types. Hypothesizing that the binary formed due to dynamical interactions during the galaxy merger, the subsequent time elapsed can constrain the delay time of the BNS coalescence. By using velocity dispersion estimates and the position of the shells, we find that the galaxy merger occurred t mer ≲ 200 Myr prior to the BNS coalescence.

We present results for the first three years of OzDES, a six year programme to obtain redshifts for objects in the Dark Energy Survey (DES) supernova fields using the 2dF fibre positioner and AAOmega spectrograph on the Anglo-Australian Telescope. OzDES is a multi-object spectroscopic survey targeting multiple types of targets at multiple epochs over a multiyear baseline and is one of the first multi-object spectroscopic surveys to dynamically include transients into the target list soon after their discovery. At the end of three years, OzDES has spectroscopically confirmed almost 100 supernovae, and has measured redshifts for 17 000 objects, including the redshifts of 2566 supernova hosts. We examine how our ability to measure redshifts for targets of various types depends on signal-to-noise ratio (S/N), magnitude and exposure time, finding that our redshift success rate increases significantly at a S/N of 2-3 per 1-Å bin. We also find that the change in S/N with exposure time closely matches the Poisson limit for stacked exposures as long as 10 h. We use these results to predict the redshift yield of the full OzDES survey, as well as the potential yields of future surveys on other facilities such as (i.e. the 4-m Multi-Object Spectroscopic Telescope, the Subaru Prime Focus Spectrograph and the Maunakea Spectroscopic Explorer). This work marks the first OzDES data release, comprising 14 693 redshifts. OzDES is on target to obtain over 30 000 redshifts over the 6-yr duration of the survey, including a yield of approximately 5700 supernova host-galaxy redshifts.

We present the Dark Energy Camera (DECam) discovery of the optical counterpart of the first binary neutron star merger detected through gravitational-wave emission, GW170817. Our observations commenced 10.5 hr post-merger, as soon as the localization region became accessible from Chile. We imaged 70 deg2 in the i and z bands, covering 93% of the initial integrated localization probability, to a depth necessary to identify likely optical counterparts (e.g., a kilonova). At 11.4 hr post-merger we detected a bright optical transient located from the nucleus of NGC 4993 at redshift z = 0.0098, consistent (for km s-1 Mpc-1) with the distance of 40 ±8 Mpc reported by the LIGO Scientific Collaboration and the Virgo Collaboration (LVC). At detection the transient had magnitudes of and , and thus an absolute magnitude of , in the luminosity range expected for a kilonova. We identified 1500 potential transient candidates. Applying simple selection criteria aimed at rejecting background events such as supernovae, we find the transient associated with NGC 4993 as the only remaining plausible counterpart, and reject chance coincidence at the 99.5% confidence level. We therefore conclude that the optical counterpart we have identified near NGC 4993 is associated with GW170817. This discovery ushers in the era of multi-messenger astronomy with gravitational waves and demonstrates the power of DECam to identify the optical counterparts of gravitational-wave sources.

We present the Dark Energy Camera (DECam) discovery of the optical counterpart of the first binary neutron star merger detected through gravitational-wave emission, GW170817. Our observations commenced 10.5 hr post-merger, as soon as the localization region became accessible from Chile. We imaged 70 deg2 in the i and z bands, covering 93% of the initial integrated localization probability, to a depth necessary to identify likely optical counterparts (e.g., a kilonova). At 11.4 hr post-merger we detected a bright optical transient located from the nucleus of NGC 4993 at redshift z = 0.0098, consistent (for km s-1 Mpc-1) with the distance of 40 ±8 Mpc reported by the LIGO Scientific Collaboration and the Virgo Collaboration (LVC). At detection the transient had magnitudes of and , and thus an absolute magnitude of , in the luminosity range expected for a kilonova. We identified 1500 potential transient candidates. Applying simple selection criteria aimed at rejecting background events such as supernovae, we find the transient associated with NGC 4993 as the only remaining plausible counterpart, and reject chance coincidence at the 99.5% confidence level. We therefore conclude that the optical counterpart we have identified near NGC 4993 is associated with GW170817. This discovery ushers in the era of multi-messenger astronomy with gravitational waves and demonstrates the power of DECam to identify the optical counterparts of gravitational-wave sources.

We present the Dark Energy Survey (DES) discovery of DES15E2mlf, the most distant superluminous supernova (SLSN) spectroscopically confirmed to date. The light curves and Gemini spectroscopy of DES15E2mlf indicate that it is a Type I superluminous supernova (SLSN-I) at z = 1.861 (a lookback time of ~10 Gyr) and peaking at MAB = -22.3 ± 0.1 mag. Given the high redshift, our data probe the rest-frame ultraviolet (1400-3500 Å) properties of the SN, finding velocity of the C III feature changes by ~5600 km s-1 over 14 d around maximum light. We find the host galaxy of DES15E2mlf has a stellar mass of 3.5+3.6-2.4 × 109 M⊙, which is more massive than the typical SLSN-I host galaxy.

We report the results of searches for strong gravitational lens systems in the Dark Energy Survey (DES) Science Verification and Year 1 observations. The Science Verification data span approximately 250 sq. deg. with a median i-band limiting magnitude for extended objects (10σ) of 23.0. The Year 1 data span approximately 2000 sq. deg. and have an i-band limiting magnitude for extended objects (10σ) of 22.9. As these data sets are both wide and deep, they are particularly useful for identifying strong gravitational lens candidates. Potential strong gravitational lens candidate systems were initially identified based on a color and magnitude selection in the DES object catalogs or because the system is at the location of a previously identified galaxy cluster. Cutout images of potential candidates were then visually scanned using an object viewer and numerically ranked according to whether or not we judged them to be likely strong gravitational lens systems. Having scanned nearly 400,000 cutouts, we present 374 candidate strong lens systems, of which 348 are identified for the first time. We provide the R.A. and decl., the magnitudes and photometric properties of the lens and source objects, and the distance (radius) of the source(s) from the lens center for each system.

We report the results of searches for strong gravitational lens systems in the Dark Energy Survey (DES) Science Verification and Year 1 observations. The Science Verification data span approximately 250 sq. deg. with a median i-band limiting magnitude for extended objects (10σ) of 23.0. The Year 1 data span approximately 2000 sq. deg. and have an i-band limiting magnitude for extended objects (10σ) of 22.9. As these data sets are both wide and deep, they are particularly useful for identifying strong gravitational lens candidates. Potential strong gravitational lens candidate systems were initially identified based on a color and magnitude selection in the DES object catalogs or because the system is at the location of a previously identified galaxy cluster. Cutout images of potential candidates were then visually scanned using an object viewer and numerically ranked according to whether or not we judged them to be likely strong gravitational lens systems. Having scanned nearly 400,000 cutouts, we present 374 candidate strong lens systems, of which 348 are identified for the first time. We provide the R.A. and decl., the magnitudes and photometric properties of the lens and source objects, and the distance (radius) of the source(s) from the lens center for each system.

We report on SPT-CLJ2011-5228, a giant system of arcs created by a cluster at z = 1.06. The arc system is notable for the presence of a bright central image. The source is a Lyman break galaxy at z s = 2.39 and the mass enclosed within the Einstein ring of radius 14 arcsec is . We perform a full reconstruction of the light profile of the lensed images to precisely infer the parameters of the mass distribution. The brightness of the central image demands that the central total density profile of the lens be shallow. By fitting the dark matter as a generalized Navarro-Frenk-White profile - with a free parameter for the inner density slope - we find that the break radius is kpc, and that the inner density falls with radius to the power -0.38 ±0.04 at 68% confidence. Such a shallow profile is in strong tension with our understanding of relaxed cold dark matter halos; dark matter-only simulations predict that the inner density should fall as . The tension can be alleviated if this cluster is in fact a merger; a two-halo model can also reconstruct the data, with both clumps (density varying as and ) much more consistent with predictions from dark matter-only simulations. At the resolution of our Dark Energy Survey imaging, we are unable to choose between these two models, but we make predictions for forthcoming Hubble Space Telescope imaging that will decisively distinguish between them.

We report the discovery of two new candidate stellar systems in the constellation of Cetus using the data from the first two years of the Dark Energy Survey (DES). The objects, DES J0111-1341 and DES J0225+0304, are located at a heliocentric distance of ~25 kpc and appear to have old and metal-poor populations. Their distances to the Sagittarius orbital plane, ~1.73 kpc (DES J0111-1341) and ~0.50 kpc (DES J0225+0304), indicate that they are possibly associated with the Sagittarius dwarf stream. The half-light radius (rh ≃ 4.55 pc) and luminosity (MV ≃ +0.3) of DES J0111-1341 are consistent with it being an ultrafaint stellar cluster, while the half-light radius (rh ≃ 18.55 pc) and luminosity (MV ≃ -1.1) of DES J0225+0304 place it in an ambiguous region of size-luminosity space between stellar clusters and dwarf galaxies. Determinations of the characteristic parameters of the Sagittarius stream, metallicity spread (-2.18 ≲ [Fe/H] ≲ -0.95) and distance gradient (23 kpc ≲ D⊙ ≲ 29 kpc), within the DES footprint in the Southern hemisphere, using the same DES data, also indicate a possible association between these systems. If these objects are confirmed through spectroscopic follow-up to be gravitationally bound systems and to share a Galactic trajectory with the Sagittarius stream, DES J0111-1341 and DES J0225+0304 would be the first ultrafaint stellar systems associated with the Sagittarius stream. Furthermore, DES J0225+0304 would also be the first confirmed case of an ultrafaint satellite of a satellite.

We report the discovery of two new candidate stellar systems in the constellation of Cetus using the data from the first two years of the Dark Energy Survey (DES). The objects, DES J0111-1341 and DES J0225+0304, are located at a heliocentric distance of ~25 kpc and appear to have old and metal-poor populations. Their distances to the Sagittarius orbital plane, ~1.73 kpc (DES J0111-1341) and ~0.50 kpc (DES J0225+0304), indicate that they are possibly associated with the Sagittarius dwarf stream. The half-light radius (rh ≃ 4.55 pc) and luminosity (MV ≃ +0.3) of DES J0111-1341 are consistent with it being an ultrafaint stellar cluster, while the half-light radius (rh ≃ 18.55 pc) and luminosity (MV ≃ -1.1) of DES J0225+0304 place it in an ambiguous region of size-luminosity space between stellar clusters and dwarf galaxies. Determinations of the characteristic parameters of the Sagittarius stream, metallicity spread (-2.18 ≲ [Fe/H] ≲ -0.95) and distance gradient (23 kpc ≲ D⊙ ≲ 29 kpc), within the DES footprint in the Southern hemisphere, using the same DES data, also indicate a possible association between these systems. If these objects are confirmed through spectroscopic follow-up to be gravitationally bound systems and to share a Galactic trajectory with the Sagittarius stream, DES J0111-1341 and DES J0225+0304 would be the first ultrafaint stellar systems associated with the Sagittarius stream. Furthermore, DES J0225+0304 would also be the first confirmed case of an ultrafaint satellite of a satellite.

The coalescence of a binary neutron star pair is expected to produce gravitational waves (GW) and electromagnetic radiation, both of which may be detectable with currently available instruments. We describe a search for a predicted r-process optical transient from these mergers, dubbed the "kilonova" (KN), using griz broadband data from the Dark Energy Survey Supernova Program (DES-SN). Some models predict KNe to be redder, shorter-lived, and dimmer than supernovae (SNe), but the event rate of KNe is poorly constrained. We simulate KN and SN light curves with the Monte-Carlo simulation code SNANA to optimize selection requirements, determine search efficiency, and predict SN backgrounds. Our analysis of the first two seasons of DES-SN data results in 0 events, and is consistent with our prediction of 1.1 ± 0.2 background events based on simulations of SNe. From our prediction, there is a 33% chance of finding 0 events in the data. Assuming no underlying galaxy flux, our search sets 90% upper limits on the KN volumetric rate of 1.0 Gpc-3 yr-1 for the dimmest KN model we consider (peak i-band absolute magnitude mag) and 2.4 Gpc-3 yr-1 for the brightest ( mag). Accounting for anomalous subtraction artifacts on bright galaxies, these limits are ∼3 times higher. This analysis is the first untriggered optical KN search and informs selection requirements and strategies for future KN searches. Our upper limits on the KN rate are consistent with those measured by GW and gamma-ray burst searches.

The coalescence of a binary neutron star pair is expected to produce gravitational waves (GW) and electromagnetic radiation, both of which may be detectable with currently available instruments. We describe a search for a predicted r-process optical transient from these mergers, dubbed the "kilonova" (KN), using griz broadband data from the Dark Energy Survey Supernova Program (DES-SN). Some models predict KNe to be redder, shorter-lived, and dimmer than supernovae (SNe), but the event rate of KNe is poorly constrained. We simulate KN and SN light curves with the Monte-Carlo simulation code SNANA to optimize selection requirements, determine search efficiency, and predict SN backgrounds. Our analysis of the first two seasons of DES-SN data results in 0 events, and is consistent with our prediction of 1.1 ± 0.2 background events based on simulations of SNe. From our prediction, there is a 33% chance of finding 0 events in the data. Assuming no underlying galaxy flux, our search sets 90% upper limits on the KN volumetric rate of 1.0 Gpc-3 yr-1 for the dimmest KN model we consider (peak i-band absolute magnitude mag) and 2.4 Gpc-3 yr-1 for the brightest ( mag). Accounting for anomalous subtraction artifacts on bright galaxies, these limits are ∼3 times higher. This analysis is the first untriggered optical KN search and informs selection requirements and strategies for future KN searches. Our upper limits on the KN rate are consistent with those measured by GW and gamma-ray burst searches.

We present the Dark Energy Survey (DES) discovery of DES15E2mlf, the most distant superluminous supernova (SLSN) spectroscopically confirmed to date. The light curves and Gemini spectroscopy of DES15E2mlf indicate that it is a Type I superluminous supernova (SLSN-I) at z = 1.861 (a lookback time of ~10 Gyr) and peaking at MAB = -22.3 ± 0.1 mag. Given the high redshift, our data probe the rest-frame ultraviolet (1400-3500 Å) properties of the SN, finding velocity of the C III feature changes by ~5600 km s-1 over 14 d around maximum light. We find the host galaxy of DES15E2mlf has a stellar mass of 3.5+3.6-2.4 × 109 M⊙, which is more massive than the typical SLSN-I host galaxy.

Host galaxy identification is a crucial step for modern supernova (SN) surveys such as the Dark Energy Survey and the Large Synoptic Survey Telescope, which will discover SNe by the thousands. Spectroscopic resources are limited, and so in the absence of real-time SN spectra these surveys must rely on host galaxy spectra to obtain accurate redshifts for the Hubble diagram and to improve photometric classification of SNe. In addition, SN luminosities are known to correlate with host-galaxy properties. Therefore, reliable identification of host galaxies is essential for cosmology and SN science. We simulate SN events and their locations within their host galaxies to develop and test methods for matching SNe to their hosts. We use both real and simulated galaxy catalog data from the Advanced Camera for Surveys General Catalog and MICECATv2.0, respectively. We also incorporate "hostless" SNe residing in undetected faint hosts into our analysis, with an assumed hostless rate of 5%. Our fully automated algorithm is run on catalog data and matches SNe to their hosts with 91% accuracy. We find that including a machine learning component, run after the initial matching algorithm, improves the accuracy (purity) of the matching to 97% with a 2% cost in efficiency (true positive rate). Although the exact results are dependent on the details of the survey and the galaxy catalogs used, the method of identifying host galaxies we outline here can be applied to any transient survey.

The joint analysis of galaxy-galaxy lensing and galaxy clustering is a promising method for inferring the growth function of large-scale structure. Anticipating a near future application of this analysis to Dark Energy Survey (DES) measurements of galaxy positions and shapes, we develop a practical approach to modeling the assumptions and systematic effects affecting the joint analysis of small-scale galaxy-galaxy lensing 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 being subdominant. 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 cover 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.

The joint analysis of galaxy-galaxy lensing and galaxy clustering is a promising method for inferring the growth function of large-scale structure. Anticipating a near future application of this analysis to Dark Energy Survey (DES) measurements of galaxy positions and shapes, we develop a practical approach to modeling the assumptions and systematic effects affecting the joint analysis of small-scale galaxy-galaxy lensing 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 being subdominant. 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 cover 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.

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 colour cuts necessary to produce a luminosity-thresholded LRG sample of constant comoving density. We demonstrate that redMaGiC photo-zs 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 catalogue sampling the redshift range z ∈ [0.2, 0.8]. Our fiducial sample has a comoving space density of 10-3 (h-1 Mpc)-3, and a median photo-z bias (zspec - zphoto) and scatter (σz/(1 + z)) of 0.005 and 0.017, respectively. The corresponding 5σ outlier fraction is 1.4 per cent. We also test our algorithm with Sloan Digital Sky Survey Data Release 8 and Stripe 82 data, and discuss how spectroscopic training can be used to control photo-z biases at the 0.1 per cent level.

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 colour cuts necessary to produce a luminosity-thresholded LRG sample of constant comoving density. We demonstrate that redMaGiC photo-zs 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 catalogue sampling the redshift range z ∈ [0.2, 0.8]. Our fiducial sample has a comoving space density of 10-3 (h-1 Mpc)-3, and a median photo-z bias (zspec - zphoto) and scatter (σz/(1 + z)) of 0.005 and 0.017, respectively. The corresponding 5σ outlier fraction is 1.4 per cent. We also test our algorithm with Sloan Digital Sky Survey Data Release 8 and Stripe 82 data, and discuss how spectroscopic training can be used to control photo-z biases at the 0.1 per cent level.