The Last Journey is a large-volume, gravity-only, cosmological N-body simulation evolving more than 1.24 trillion particles in a periodic box with a side-length of 5.025Gpc. It was implemented using the HACC simulation and analysis framework on the BG/Q system, Mira. The cosmological parameters are chosen to be consistent with the results from the Planck satellite. A range of analysis tools have been run in situ to enable a diverse set of science projects, and at the same time, to keep the resulting data amount manageable. Analysis outputs have been generated starting at redshift z~10 to allow for construction of synthetic galaxy catalogs using a semi-analytic modeling approach in post-processing. As part of our in situ analysis pipeline we employ a new method for tracking halo sub-structures, introducing the concept of subhalo cores. The production of multi-wavelength synthetic sky maps is facilitated by generating particle lightcones in situ, also beginning at z~10. We provide an overview of the simulation set-up and the generated data products; a first set of analysis results is presented. A subset of the data is publicly available.

Many scientific investigations of photometric galaxy surveys require redshift estimates, whose uncertainty properties are best encapsulated by photometric redshift (photo-z) posterior probability density functions (PDFs). A plethora of photo-z PDF estimation methodologies abound, producing discrepant results with no consensus on a preferred approach. We present the results of a comprehensive experiment comparing twelve photo-z algorithms applied to mock data produced for the Large Synoptic Survey Telescope (LSST) Dark Energy Science Collaboration (DESC). By supplying perfect prior information, in the form of the complete template library and a representative training set as inputs to each code, we demonstrate the impact of the assumptions underlying each technique on the output photo-z PDFs. In the absence of a notion of true, unbiased photo-z PDFs, we evaluate and interpret multiple metrics of the ensemble properties of the derived photo-z PDFs as well as traditional reductions to photo-z point estimates. We report systematic biases and overall over/under-breadth of the photo-z PDFs of many popular codes, which may indicate avenues for improvement in the algorithms or implementations. Furthermore, we raise attention to the limitations of established metrics for assessing photo-z PDF accuracy; though we identify the conditional density estimate (CDE) loss as a promising metric of photo-z PDF performance in the case where true redshifts are available but true photo-z PDFs are not, we emphasize the need for science-specific performance metrics.

As part of the effort to meet the needs of the Large Synoptic Survey Telescope Dark Energy Science Collaboration (LSST DESC) for accurate, realistically complex mock galaxy catalogs, we have developed GalSampler, an open-source python package that assists in generating large volumes of synthetic cosmological data. The key idea behind GalSampler is to recast hydrodynamical simulations and semi-analytic models as physically-motivated galaxy libraries. GalSampler populates a new, larger-volume halo catalog with galaxies drawn from the baseline library; by using weighted sampling guided by empirical modeling techniques, GalSampler inherits statistical accuracy from the empirical model and physically-motivated complexity from the baseline library. We have recently used GalSampler to produce the cosmoDC2 extragalactic catalog made for the LSST DESC Data Challenge 2. Using cosmoDC2 as a guiding example, we outline how GalSampler can continue to support ongoing and near-future galaxy surveys such as the Dark Energy Survey (DES), the Dark Energy Spectroscopic Instrument (DESI), WFIRST, and Euclid.

This paper introduces cosmoDC2, a large synthetic galaxy catalog designed to support precision dark energy science with the Large Synoptic Survey Telescope (LSST). CosmoDC2 is the starting point for the second data challenge (DC2) carried out by the LSST Dark Energy Science Collaboration (LSST DESC). The catalog is based on a trillion-particle, 4.225 Gpc^3 box cosmological N-body simulation, the `Outer Rim' run. It covers 440 deg^2 of sky area to a redshift of z=3 and is complete to a magnitude depth of 28 in the r-band. Each galaxy is characterized by a multitude of properties including stellar mass, morphology, spectral energy distributions, broadband filter magnitudes, host halo information and weak lensing shear. The size and complexity of cosmoDC2 requires an efficient catalog generation methodology; our approach is based on a new hybrid technique that combines data-driven empirical approaches with semi-analytic galaxy modeling. A wide range of observation-based validation tests has been implemented to ensure that cosmoDC2 enables the science goals of the planned LSST DESC DC2 analyses. This paper also represents the official release of the cosmoDC2 data set, including an efficient reader that facilitates interaction with the data.

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 LCDM model we find a matter density Omega_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 Omega_m = 0.321 +_ 0.018. For a flat w0waCDM model, and combining probes from SN Ia, CMB and baryon acoustic oscillations, we find w0 = -0.885 +_ 0.114 and wa = -0.387 +_ 0.430. These results are in agreement with a cosmological constant and with previous constraints using SNe Ia (Pantheon, JLA).

Supernova (SN) classification and redshift estimation using photometric data only have become very important for the Large Synoptic Survey Telescope (LSST), given the large number of SNe that LSST will observe and the impossibility of spectroscopically following up all the SNe. We investigate the performance of a SN classifier that uses SN colors to classify LSST SNe with the Random Forest classification algorithm. Our classifier results in an AUC of 0.98 which represents excellent classification. We are able to obtain a photometric SN sample containing 99$\%$ SNe Ia by choosing a probability threshold. We estimate the photometric redshifts (photo-z) of SNe in our sample by fitting the SN light curves using the SALT2 model with nested sampling. We obtain a mean bias ($\left<z_\mathrm{phot}-z_\mathrm{spec}\right>$) of 0.012 with $\sigma\left( \frac{z_\mathrm{phot}-z_\mathrm{spec}}{1+z_\mathrm{spec}}\right) = 0.0294$ without using a host-galaxy photo-z prior, and a mean bias ($\left<z_\mathrm{phot}-z_\mathrm{spec}\right>$) of 0.0017 with $\sigma\left( \frac{z_\mathrm{phot}-z_\mathrm{spec}}{1+z_\mathrm{spec}}\right) = 0.0116$ using a host-galaxy photo-z prior. Assuming a flat $\Lambda CDM$ model with $\Omega_m=0.3$, we obtain $\Omega_m$ of $0.305\pm0.008$ (statistical errors only), using the simulated LSST sample of photometric SNe Ia (with intrinsic scatter $\sigma_\mathrm{int}=0.11$) derived using our methodology without using host-galaxy photo-z prior. Our method will help boost the power of SNe from the LSST as cosmological probes.

The use of high-quality simulated sky catalogs is essential for the success of cosmological surveys. The catalogs have diverse applications, such as investigating signatures of fundamental physics in cosmological observables, understanding the effect of systematic uncertainties on measured signals and testing mitigation strategies for reducing these uncertainties, aiding analysis pipeline development and testing, and survey strategy optimization. The list of applications is growing with improvements in the quality of the catalogs and the details that they can provide. Given the importance of simulated catalogs, it is critical to provide rigorous validation protocols that enable both catalog providers and users to assess the quality of the catalogs in a straightforward and comprehensive way. For this purpose, we have developed the DESCQA framework for the Large Synoptic Survey Telescope Dark Energy Science Collaboration as well as for the broader community. The goal of DESCQA is to enable the inspection, validation, and comparison of an inhomogeneous set of synthetic catalogs via the provision of a common interface within an automated framework. In this paper, we present the design concept and first implementation of DESCQA. In order to establish and demonstrate its full functionality we use a set of interim catalogs and validation tests. We highlight several important aspects, both technical and scientific, that require thoughtful consideration when designing a validation framework, including validation metrics and how these metrics impose requirements on the synthetic sky catalogs.

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 multi-year 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 2,566 supernova hosts. We examine how our ability to measure redshifts for targets of various types depends on signal-to-noise, magnitude, and exposure time, finding that our redshift success rate increases significantly at a signal-to-noise of 2 to 3 per 1-Angstrom bin. We also find that the change in signal-to-noise with exposure time closely matches the Poisson limit for stacked exposures as long as 10 hours. 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 the 4m Multi-Object Spectroscopic Telescope (4MOST), the Subaru Prime Focus Spectrograph (PFS), and the Maunakea Spectroscopic Explorer (MSE). This work marks the first OzDES data release, comprising 14,693 redshifts. OzDES is on target to obtain over a yield of approximately 5,700 supernova host-galaxy redshifts.

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 M_AB = -22.3 +/- 0.1 mag. Given the high redshift, our data probe the rest-frame ultraviolet (1400-3500 A) properties of the SN, finding velocity of the C III feature changes by ~5600 km/s over 14 days around maximum light. We find the host galaxy of DES15E2mlf has a stellar mass of 3.5^+3.6_-2.4 x 10^9 M_sun, which is more massive than the typical SLSN-I host galaxy.

DESI (Dark Energy Spectropic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. The DESI instrument is a robotically-actuated, fiber-fed spectrograph capable of taking up to 5,000 simultaneous spectra over a wavelength range from 360 nm to 980 nm. The fibers feed ten three-arm spectrographs with resolution $R= \lambda/\Delta\lambda$ between 2000 and 5500, depending on wavelength. The DESI instrument will be used to conduct a five-year survey designed to cover 14,000 deg$^2$. This powerful instrument will be installed at prime focus on the 4-m Mayall telescope in Kitt Peak, Arizona, along with a new optical corrector, which will provide a three-degree diameter field of view. The DESI collaboration will also deliver a spectroscopic pipeline and data management system to reduce and archive all data for eventual public use.

DESI (Dark Energy Spectroscopic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations (BAO) and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. To trace the underlying dark matter distribution, spectroscopic targets will be selected in four classes from imaging data. We will measure luminous red galaxies up to $z=1.0$. To probe the Universe out to even higher redshift, DESI will target bright [O II] emission line galaxies up to $z=1.7$. Quasars will be targeted both as direct tracers of the underlying dark matter distribution and, at higher redshifts ($ 2.1 < z < 3.5$), for the Ly-$\alpha$ forest absorption features in their spectra, which will be used to trace the distribution of neutral hydrogen. When moonlight prevents efficient observations of the faint targets of the baseline survey, DESI will conduct a magnitude-limited Bright Galaxy Survey comprising approximately 10 million galaxies with a median $z\approx 0.2$. In total, more than 30 million galaxy and quasar redshifts will be obtained to measure the BAO feature and determine the matter power spectrum, including redshift space distortions.

Host galaxy identification is a crucial step for modern supernova (SN) surveys such as the Dark Energy Survey (DES) and the Large Synoptic Survey Telescope (LSST), which will discover SNe by the thousands. Spectroscopic resources are limited, 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.

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_U_peak = -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 (sub-solar), low stellar-mass host galaxy (log(M/M_sun) = 9.3 +/- 0.3); consistent with what is seen for other SLSNe-I. We compare the bolometric light curve of DES13S2cmm to fourteen similarly well-observed SLSNe-I in the literature and find it possesses one of the slowest declining tails (beyond +30 days 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 magnitudes between +25 and +30 days after peak (rest frame) depending on redshift range studied; this could be important for 'standardising' 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 days in the rest frame of the supernova.

We present an analysis of supernova light curves simulated for the upcoming Dark Energy Survey (DES) supernova search. The simulations employ a code suite that generates and fits realistic light curves in order to obtain distance modulus/redshift pairs that are passed to a cosmology fitter. We investigated several different survey strategies including field selection, supernova selection biases, and photometric redshift measurements. Using the results of this study, we chose a 30 square degree search area in the griz filter set. We forecast 1) that this survey will provide a homogeneous sample of up to 4000 Type Ia supernovae in the redshift range 0.05<z<1.2, and 2) that the increased red efficiency of the DES camera will significantly improve high-redshift color measurements. The redshift of each supernova with an identified host galaxy will be obtained from spectroscopic observations of the host. A supernova spectrum will be obtained for a subset of the sample, which will be utilized for control studies. In addition, we have investigated the use of combined photometric redshifts taking into account data from both the host and supernova. We have investigated and estimated the likely contamination from core-collapse supernovae based on photometric identification, and have found that a Type Ia supernova sample purity of up to 98% is obtainable given specific assumptions. Furthermore, we present systematic uncertainties due to sample purity, photometric calibration, dust extinction priors, filter-centroid shifts, and inter-calibration. We conclude by estimating the uncertainty on the cosmological parameters that will be measured from the DES supernova data.

We present the results of a study of selection criteria to identify Type Ia supernovae photometrically in a simulated mixed sample of Type Ia supernovae and core collapse supernovae. The simulated sample is a mockup of the expected results of the Dark Energy Survey. Fits to the MLCS2k2 and SALT2 Type Ia supernova models are compared and used to help separate the Type Ia supernovae from the core collapse sample. The Dark Energy Task Force Figure of Merit (modified to include core collapse supernovae systematics) is used to discriminate among the various selection criteria. This study of varying selection cuts for Type Ia supernova candidates is the first to evaluate core collapse contamination using the Figure of Merit. Different factors that contribute to the Figure of Merit are detailed. With our analysis methods, both SALT2 and MLCS2k2 Figures of Merit improve with tighter selection cuts and higher purities, peaking at 98% purity.

We describe a dedicated cosmic-ray telescope that explores a new method for detecting Cerenkov radiation from high-energy primary cosmic rays and the large particle air shower they induce upon entering the atmosphere. Using a camera comprising 16 multi-anode photomultiplier tubes for a total of 256 pixels, the Track Imaging Cerenkov Experiment (TrICE) resolves substructures in particle air showers with 0.086 degree resolution. Cerenkov radiation is imaged using a novel two-part optical system in which a Fresnel lens provides a wide-field optical trigger and a mirror system collects delayed light with four times the magnification. TrICE records well-resolved cosmic-ray air showers at rates ranging between 0.01-0.1 Hz.

The Track Imaging Cerenov Experiment (TrICE) is an air Cerenkov prototype telescope designed to use multi-anode photomultiplier to acheive a high angular resolution for measuring cosmic-ray composition at TeV-PeV energies. The TrICE camera, composed of 16 Hamamatsu R8900 16-channel multi-anode photomultiplier tubes, achieves 0.086 degree angular width per pixel over 1.5 degree wide field of view. We present a description of the TrICE camera design, calibration and performance.