We present a new Bayesian hierarchical model (BHM) named Steve for performing Type Ia supernova (SN Ia) cosmology fits. This advances previous works by including an improved treatment of Malmquist bias, accounting for additional sources of systematic uncertainty, and increasing numerical efficiency. Given light-curve fit parameters, redshifts, and host-galaxy masses, we fit Steve simultaneously for parameters describing cosmology, SN Ia populations, and systematic uncertainties. Selection effects are characterized using Monte Carlo simulations. We demonstrate its implementation by fitting realizations of SN Ia data sets where the SN Ia model closely follows that used in Steve. Next, we validate on more realistic SNANA simulations of SN Ia samples from the Dark Energy Survey and low-redshift surveys (DES Collaboration et al. 2018). These simulated data sets contain more than 60,000 SNe Ia, which we use to evaluate biases in the recovery of cosmological parameters, specifically the equation of state of dark energy, w. This is the most rigorous test of a BHM method applied to SN Ia cosmology fitting and reveals small w biases that depend on the simulated SN Ia properties, in particular the intrinsic SN Ia scatter model. This w bias is less than 0.03 on average, less than half the statistical uncertainty on w. These simulation test results are a concern for BHM cosmology fitting applications on large upcoming surveys; therefore, future development will focus on minimizing the sensitivity of Steve to the SN Ia intrinsic scatter model.

We describe an algorithm for identifying point-source transients and moving objects on reference-subtracted optical images containing artifacts of processing and instrumentation. The algorithm makes use of the supervised machine learning technique known as Random Forest. We present results from its use in the Dark Energy Survey Supernova program (DES-SN), where it was trained using a sample of 898,963 signal and background events generated by the transient detection pipeline. After reprocessing the data collected during the first DES-SN observing season (2013 September through 2014 February) using the algorithm, the number of transient candidates eligible for human scanning decreased by a factor of 13.4, while only 1.0% of the artificial Type Ia supernovae (SNe) injected into search images to monitor survey efficiency were lost, most of which were very faint events. Here we characterize the algorithm's performance in detail, and we discuss how it can inform pipeline design decisions for future time-domain imaging surveys, such as the Large Synoptic Survey Telescope and the Zwicky Transient Facility. An implementation of the algorithm and the training data used in this paper are available at http://portal.nserc.gov/project/dessn/autoscan.

We use numerical simulations to characterize the performance of a clustering-based method to calibrate photometric redshift biases. In particular, we cross-correlate the weak lensing source galaxies from the Dark Energy Survey Year 1 sample with redMaGiC galaxies (luminous red galaxies with secure photometric redshifts) to estimate the redshift distribution of the former sample. The recovered redshift distributions are used to calibrate the photometric redshift bias of standard photo-z methods applied to the same source galaxy sample.We apply the method to two photo-z codes run in our simulated data: Bayesian Photometric Redshift and Directional Neighbourhood Fitting. We characterize the systematic uncertainties of our calibration procedure, and find that these systematic uncertainties dominate our error budget. The dominant systematics are due to our assumption of unevolving bias and clustering across each redshift bin, and to differences between the shapes of the redshift distributions derived by clustering versus photo-zs. The systematic uncertainty in the mean redshift bias of the source galaxy sample is Δz ≲ 0.02, though the precise value depends on the redshift bin under consideration. We discuss possible ways to mitigate the impact of our dominant systematics in future analyses.

We present DES14X3taz, a new hydrogen-poor superluminous supernova (SLSN-I) discovered by the Dark Energy Survey (DES) supernova program, with additional photometric data provided by the Survey Using DECam for Superluminous Supernovae. Spectra obtained using Optical System for Imaging and low-Intermediate-Resolution Integrated Spectroscopy on the Gran Telescopio CANARIAS show DES14X3taz is an SLSN-I at z = 0.608. Multi-color photometry reveals a double-peaked light curve: a blue and relatively bright initial peak that fades rapidly prior to the slower rise of the main light curve. Our multi-color photometry allows us, for the first time, to show that the initial peak cools from 22,000 to 8000 K over 15 rest-frame days, and is faster and brighter than any published core-collapse supernova, reaching 30% of the bolometric luminosity of the main peak. No physical 56Ni-powered model can fit this initial peak. We show that a shock-cooling model followed by a magnetar driving the second phase of the light curve can adequately explain the entire light curve of DES14X3taz. Models involving the shock-cooling of extended circumstellar material at a distance of ≃400 R⊙ are preferred over the cooling of shock-heated surface layers of a stellar envelope. We compare DES14X3taz to the few double-peaked SLSN-I events in the literature. Although the rise times and characteristics of these initial peaks differ, there exists the tantalizing possibility that they can be explained by one physical interpretation.

We present details on the observing strategy, data-processing techniques, and spectroscopic targeting algorithms for the first three years of operation for the Dark Energy Survey Supernova Program (DES-SN). This five-year program using the Dark Energy Camera mounted on the 4 m Blanco telescope in Chile was designed to discover and follow supernovae (SNe) Ia over a wide redshift range (0.05 < z < 1.2) to measure the equation-of-state parameter of dark energy. We describe the SN program in full: Strategy, observations, data reduction, spectroscopic follow-up observations, and classification. From three seasons of data, we have discovered 12,015 likely SNe, 308 of which have been spectroscopically confirmed, including 251 SNe Ia over a redshift range of 0.017 < z < 0.85. We determine the effective spectroscopic selection function for our sample and use it to investigate the redshiftdependent bias on the distance moduli of SNe Ia we have classified. The data presented here are used for the first cosmology analysis by DES-SN ("DES-SN3YR"), the results of which are given in Dark Energy Survey Collaboration et al. The 489 spectra that are used to define the DES-SN3YR sample are publicly available 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.

We describe the derivation and validation of redshift distribution estimates and their uncertainties for the populations of galaxies used as weak-lensing sources in the Dark Energy Survey (DES) Year 1 cosmological analyses. The Bayesian Photometric Redshift (BPZ) code is used to assign galaxies to four redshift bins between z ≈ 0.2 and ≈1.3, and to produce initial estimates of the lensing-weighted redshift distributions nPZi(z) ∝ dni/dz for members of bin i. Accurate determination of cosmological parameters depends critically on knowledge of ni, but is insensitive to bin assignments or redshift errors for individual galaxies. The cosmological analyses allow for shifts ni (z) = nPZi(z - Δzi) to correct themean redshift of ni(z) for biases in nPZi. The Δzi are constrained by comparison of independently estimated 30-band photometric redshifts of galaxies in the Cosmic Evolution Survey (COSMOS) field to BPZ estimates made from the DES griz fluxes, for a sample matched in fluxes, pre-seeing size, and lensing weight to the DES weak-lensing sources. In companion papers, the Δzi of the three lowest redshift bins are further constrained by the angular clustering of the source galaxies around red galaxies with secure photometric redshifts at 0.15 < z < 0.9. This paper details the BPZ and COSMOS procedures, and demonstrates that the cosmological inference is insensitive to details of the ni(z) beyond the choice of Δzi. The clustering and COSMOS validation methods produce consistent estimates of Δzi in the bins where both can be applied, with combined uncertainties of σΔzi = 0.015, 0.013, 0.011, and 0.022 in the four bins. Repeating the photo-z procedure instead using the Directional Neighbourhood Fitting algorithm, or using the ni(z) estimated from the matched sample in COSMOS, yields no discernible difference in cosmological inferences.

We present the first cosmological parameter constraints using measurements of type Ia supernovae (SNe Ia) from the Dark Energy Survey Supernova Program (DES-SN). The analysis uses a subsample of 207 spectroscopically confirmed SNe Ia from the first three years of DES-SN, combined with a low-redshift sample of 122 SNe from the literature. Our "DES-SN3YR" result from these 329 SNe Ia is based on a series of companion analyses and improvements covering SN Ia discovery, spectroscopic selection, photometry, calibration, distance bias corrections, and evaluation of systematic uncertainties. For a flat ΛCDM model we find a matter density Ω m = 0.331 ± 0.038. For a flat wCDM model, and combining our SN Ia constraints with those from the cosmic microwave background (CMB), we find a dark energy equation of state w = -0.978 ± 0.059, and ω m = 0.321 ± 0.018. For a flat w 0 w a CDM model, and combining probes from SN Ia, CMB and baryon acoustic oscillations, we find w0 = -0.885 ± 0.114 and w a = -0.387 ± 0.430. These results are in agreement with a cosmological constant and with previous constraints using SNe Ia (Pantheon, JLA).