As the scale of cosmological surveys increases, so does the complexity in the analyses. This complexity can often make it difficult to derive the underlying principles, necessitating statistically rigorous testing to ensure the results of an analysis are consistent and reasonable. This is particularly important in multi-probe cosmological analyses like those used in the Dark Energy Survey and the upcoming Legacy Survey of Space and Time, where accurate uncertainties are vital. In this paper, we present a statistically rigorous method to test the consistency of contours produced in these analyses, and apply this method to the Pippin cosmological pipeline used for Type Ia supernova cosmology with the Dark Energy Survey. We make use of the Neyman construction, a frequentist methodology that leverages extensive simulations to calculate confidence intervals, to perform this consistency check. A true Neyman construction is too computationally expensive for supernova cosmology, so we develop a method for approximating a Neyman construction with far fewer simulations. We find that for a simulated data-set, the 68% contour reported by the Pippin pipeline and the 68% confidence region produced by our approximate Neyman construction differ by less than a percent near the input cosmology, however show more significant differences far from the input cosmology, with a maximal difference of 0.05 in ωM, and 0.07 in w. This divergence is most impactful for analyses of cosmological tensions, but its impact is mitigated when combining supernovae with other cross-cutting cosmological probes, such as the Cosmic Microwave Background.

We present a near-infrared spectral sequence of the electromagnetic counterpart to the binary neutron star merger GW170817 detected by Advanced Laser Interferometer Gravitational-wave Observatory (LIGO)/Virgo. Our data set comprises seven epochs of J+H spectra taken with FLAMINGOS-2 on Gemini-South between 1.5 and 10.5 days after the merger. In the initial epoch, the spectrum is dominated by a smooth blue continuum due to a high-velocity, lanthanide-poor blue kilonova component. Starting the following night, all of the subsequent spectra instead show features that are similar to those predicted in model spectra of material with a high concentration of lanthanides, including spectral peaks near 1.07 and 1.55 μm. Our fiducial model with 0.04 M o of ejecta, an ejection velocity of v = 0.1c, and a lanthanide concentration of X lan = 10-2 provides a good match to the spectra taken in the first five days, although it over-predicts the late-time fluxes. We also explore models with multiple fitting components, in each case finding that a significant abundance of lanthanide elements is necessary to match the broad spectral peaks that we observe starting at 2.5 days after the merger. These data provide direct evidence that binary neutron star mergers are significant production sites of even the heaviest r-process elements.

This document presents the results from the Distances subgroup of the Cosmic Frontier Community Planning Study (Snowmass 2013). We summarize the current state of the field as well as future prospects and challenges. In addition to the established probes using Type Ia supernovae and baryon acoustic oscillations, we also consider prospective methods based on clusters, active galactic nuclei, gravitational wave sirens and strong lensing time delays.

In the next decade Type Ia supernovae (SNe Ia) will be used to test theories predicting changes in the Dark Energy equation of state with time. Ultimately this requires a dedicated space mission like JDEM. SNe Ia are mature cosmological probes --- their limitations are well characterized, and a path to improvement is clear. Dominant systematic errors include photometric calibration, selection effects, reddening, and population-dependent differences. Building on past lessons, well-controlled new surveys are poised to make strides in these areas: the Palomar Transient Factory, Skymapper, La Silla QUEST, Pan-STARRS, the Dark Energy Survey, LSST, and JDEM. They will obviate historical calibrations and selection biases, and allow comparisons via large subsamples. Some systematics follow from our ignorance of SN Ia progenitors, which there is hope of determining with SN Ia rate studies from 0

The Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) will observe several Deep Drilling Fields (DDFs) to a greater depth and with a more rapid cadence than the main survey. In this paper, we describe the 'DeepDrill' survey, which used the Spitzer Space Telescope Infrared Array Camera (IRAC) to observe three of the four currently defined DDFs in two bands, centred on 3.6 and 4.5 μm. These observations expand the area that was covered by an earlier set of observations in these three fields by the Spitzer Extragalactic Representative Volume Survey (SERVS). The combined DeepDrill and SERVS data cover the footprints of the LSST DDFs in the Extended Chandra Deep Field-South (ECDFS) field, the ELAIS-S1 field (ES1), and the XMM-Large-Scale Structure Survey field (XMM-LSS). The observations reach an approximate 5σ point-source depth of 2 μJy (corresponding to an AB magnitude of 23.1; sufficient to detect a 1011 M⊙ galaxy out to z ≈ 5) in each of the two bands over a total area of ≈ 29 deg2. The dual-band catalogues contain a total of 2.35 million sources. In this paper, we describe the observations and data products from the survey, and an overview of the properties of galaxies in the survey. We compare the source counts to predictions from the Shark semi-analytic model of galaxy formation. We also identify a population of sources with extremely red ([3.6]-[4.5] >1.2) colours which we show mostly consists of highly obscured active galactic nuclei.

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 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 present spectroscopic confirmation of two new lensed quasars via data obtained at the 6.5m Magellan/Baade Telescope. The lens candidates have been selected from the Dark Energy Survey (DES) and WISE based on their multi-band photometry and extended morphology in DES images. Images of DES J0115-5244 show two blue point sources at either side of a red galaxy. Our long-slit data confirm that both point sources are images of the same quasar at $z_{s}=1.64.$ The Einstein Radius estimated from the DES images is $0.51$". DES J2200+0110 is in the area of overlap between DES and the Sloan Digital Sky Survey (SDSS). Two blue components are visible in the DES and SDSS images. The SDSS fiber spectrum shows a quasar component at $z_{s}=2.38$ and absorption compatible with Mg II and Fe II at $z_{l}=0.799$, which we tentatively associate with the foreground lens galaxy. The long-slit Magellan spectra show that the blue components are resolved images of the same quasar. The Einstein Radius is $0.68$" corresponding to an enclosed mass of $1.6\times10^{11}\,M_{\odot}.$ Three other candidates were observed and rejected, two being low-redshift pairs of starburst galaxies, and one being a quasar behind a blue star. These first confirmation results provide an important empirical validation of the data-mining and model-based selection that is being applied to the entire DES dataset.

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 present DES13S2cmm, the first spectroscopically-confirmed superluminous supernova (SLSN) from the Dark Energy Survey (DES).We briefly discuss the data and search algorithm used to find this event in the first year of DES operations, and outline the spectroscopic data obtained from the European Southern Observatory (ESO) Very Large Telescope to confirm its redshift (z = 0.663 ± 0.001 based on the host-galaxy emission lines) and likely spectral type (Type I). Using this redshift, we find M peakU = -21.05+0.10-0.09 for the peak, rest-frame U-band absolute magnitude, and find DES13S2cmm to be located in a faint, low-metallicity (subsolar), low stellar-mass host galaxy (log (M/M⊙) = 9.3 ± 0.3), consistent with what is seen for other SLSNe-I. We compare the bolometric light curve of DES13S2cmm to 14 similarly well-observed SLSNe-I in the literature and find that it possesses one of the slowest declining tails (beyond +30 d rest-frame past peak), and is the faintest at peak. Moreover, we find the bolometric light curves of all SLSNe-I studied herein possess a dispersion of only 0.2-0.3 mag between +25 and +30 d after peak (rest frame) depending on redshift range studied; this could be important for 'standardizing' such supernovae, as is done with the more common Type Ia. We fit the bolometric light curve of DES13S2cmm with two competing models for SLSNe-I - the radioactive decay of 56Ni, and a magnetar - and find that while the magnetar is formally a better fit, neither model provides a compelling match to the data. Although we are unable to conclusively differentiate between these two physical models for this particular SLSN-I, further DES observations of more SLSNe-I should break this degeneracy, especially if the light curves of SLSNe-I can be observed beyond 100 d in the rest frame of the supernova.