We describe the simulated data sample for the "Photometric LSST Astronomical Time Series Classification Challenge" (PLAsTiCC), a publicly available challenge to classify transient and variable events that will be observed by the Large Synoptic Survey Telescope (LSST), a new facility expected to start in the early 2020s. The challenge was hosted by Kaggle, ran from 2018 September 28 to 2018 December 17, and included 1,094 teams competing for prizes. Here we provide details of the 18 transient and variable source models, which were not revealed until after the challenge, and release the model libraries at https://doi.org/10.5281/zenodo.2612896. We describe the LSST Operations Simulator used to predict realistic observing conditions, and we describe the publicly available SNANA simulation code used to transform the models into observed fluxes and uncertainties in the LSST passbands (ugrizy). Although PLAsTiCC has finished, the publicly available models and simulation tools are being used within the astronomy community to further improve classification, and to study contamination in photometrically identified samples of type Ia supernova used to measure properties of dark energy. Our simulation framework will continue serving as a platform to improve the PLAsTiCC models, and to develop new models.

Recently, there have been two landmark discoveries of gravitationally lensed supernovae: the first multiply-imaged SN, "Refsdal", and the first Type Ia SN resolved into multiple images, SN iPTF16geu. Fitting the multiple light curves of such objects can deliver measurements of the lensing time delays, which are the difference in arrival times for the separate images. These measurements provide precise tests of lens models or constraints on the Hubble constant and other cosmological parameters that are independent of the local distance ladder. Over the next decade, accurate time delay measurements will be needed for the tens to hundreds of lensed SNe to be found by wide-field time-domain surveys such as LSST and WFIRST. We have developed an open source software package for simulations and time delay measurements of multiply-imaged SNe, including an improved characterization of the uncertainty caused by microlensing. We describe simulations using the package that suggest a before-peak detection of the leading image enables a more accurate and precise time delay measurement (by ~1 and ~2 days, respectively), when compared to an after-peak detection. We also conclude that fitting the effects of microlensing without an accurate prior often leads to biases in the time delay measurement and over-fitting to the data, but that employing a Gaussian Process Regression (GPR) technique is sufficient for determining the uncertainty due to microlensing.

Empirical models of supernova (SN) spectral energy distributions (SEDs) are widely used for SN survey simulations and photometric classifications. The existing library of SED models has excellent optical templates but limited, poorly constrained coverage of ultraviolet (UV) and infrared (IR) wavelengths. However, both regimes are critical for the design and operation of future SN surveys, particularly at IR wavelengths that will be accessible with the James Webb Space Telescope (JWST) and the Wide-Field Infrared Survey Telescope (WFIRST). We create a public repository of improved empirical SED templates using a sampling of Type Ia and core-collapse (CC) photometric light curves to extend the Type Ia parameterized SALT2 model and a set of SN Ib, SN Ic, and SN II SED templates into the UV and near-IR. We apply this new repository of extrapolated SN SED models to examine how future surveys can discriminate between CC and Type Ia SNe at UV and IR wavelengths, and present an open-source software package written in Python, SNSEDextend, that enables a user to generate their own extrapolated SEDs.