The design and analysis of time-domain sky surveys requires the ability to simulate accurately realistic populations of core collapse supernova (SN) events. We present a set of spectral time-series templates designed for this purpose, for both hydrogen-rich (type II, IIn, IIb) and stripped envelope (types Ib, Ic, Ic-BL) core collapse supernovae. We use photometric and spectroscopic data for 67 core collapse supernovae from the literature, and for each generate a time-series spectral template. The techniques used to build the templates are fully data-driven with no assumption of any parametric form or model for the light curves. The template-building code is open-source, and can be applied to any transient for which well-sampled multi-band photometry and multiple spectroscopic observations are available. We extend these spectral templates into the near-ultraviolet to $\lambda \lambda \sim 1600 \AA$ using observer-frame ultraviolet photometry. We also provide a set of templates corrected for host galaxy dust extinction, and provide a set of luminosity functions that can be used with our spectral templates in simulations. We give an example of how these templates can be used by integrating them within the popular SN simulation package $\textsc{snana}$, and simulating core collapse supernovae in photometrically-selected cosmological type Ia supernova samples, prone to contamination from core collapse events.

Gravitational waves were discovered with the detection of binary black hole mergers and they should also be detectable from lower mass neutron star mergers. These are predicted to eject material rich in heavy radioactive isotopes that can power an electromagnetic signal called a kilonova. The gravitational wave source GW170817 arose from a binary neutron star merger in the nearby Universe with a relatively well confined sky position and distance estimate. Here we report observations and physical modelling of a rapidly fading electromagnetic transient in the galaxy NGC4993, which is spatially coincident with GW170817 and a weak short gamma-ray burst. The transient has physical parameters broadly matching the theoretical predictions of blue kilonovae from neutron star mergers. The emitted electromagnetic radiation can be explained with an ejected mass of 0.04 +/- 0.01 Msol, with an opacity of kappa <= 0.5 cm2/gm at a velocity of 0.2 +/- 0.1c. The power source is constrained to have a power law slope of beta = -1.2 +/- 0.3, consistent with radioactive powering from r-process nuclides. We identify line features in the spectra that are consistent with light r-process elements (90 < A < 140). As it fades, the transient rapidly becomes red, and emission may have contribution by a higher opacity, lanthanide-rich ejecta component. This indicates that neutron star mergers produce gravitational waves, radioactively powered kilonovae, and are a nucleosynthetic source of the r-process elements.

The Large Synoptic Survey Telescope is designed to provide an unprecedented optical imaging dataset that will support investigations of our Solar System, Galaxy and Universe, across half the sky and over ten years of repeated observation. However, exactly how the LSST observations will be taken (the observing strategy or "cadence") is not yet finalized. In this dynamically-evolving community white paper, we explore how the detailed performance of the anticipated science investigations is expected to depend on small changes to the LSST observing strategy. Using realistic simulations of the LSST schedule and observation properties, we design and compute diagnostic metrics and Figures of Merit that provide quantitative evaluations of different observing strategies, analyzing their impact on a wide range of proposed science projects. This is work in progress: we are using this white paper to communicate to each other the relative merits of the observing strategy choices that could be made, in an effort to maximize the scientific value of the survey. The investigation of some science cases leads to suggestions for new strategies that could be simulated and potentially adopted. Notably, we find motivation for exploring departures from a spatially uniform annual tiling of the sky: focusing instead on different parts of the survey area in different years in a "rolling cadence" is likely to have significant benefits for a number of time domain and moving object astronomy projects. The communal assembly of a suite of quantified and homogeneously coded metrics is the vital first step towards an automated, systematic, science-based assessment of any given cadence simulation, that will enable the scheduling of the LSST to be as well-informed as possible.

We present photometry and time-series spectroscopy of the nearby type Ia supernova (SN Ia) SN 2015F over $-16$ days to $+80$ days relative to maximum light, obtained as part of the Public ESO Spectroscopic Survey of Transient Objects (PESSTO). SN 2015F is a slightly sub-luminous SN Ia with a decline rate of $\Delta m15(B)=1.35 \pm 0.03$ mag, placing it in the region between normal and SN 1991bg-like events. Our densely-sampled photometric data place tight constraints on the epoch of first light and form of the early-time light curve. The spectra exhibit photospheric C II $\lambda 6580$ absorption until $-4$ days, and high-velocity Ca II is particularly strong at $<-10$ days at expansion velocities of $\simeq$23000\kms. At early times, our spectral modelling with syn++ shows strong evidence for iron-peak elements (Fe II, Cr II, Ti II, and V II) expanding at velocities $>14000$ km s$^{-1}$, suggesting mixing in the outermost layers of the SN ejecta. Although unusual in SN Ia spectra, including V II in the modelling significantly improves the spectral fits. Intriguingly, we detect an absorption feature at $\sim$6800 \AA\ that persists until maximum light. Our favoured explanation for this line is photospheric Al II, which has never been claimed before in SNe Ia, although detached high-velocity C II material could also be responsible. In both cases the absorbing material seems to be confined to a relatively narrow region in velocity space. The nucleosynthesis of detectable amounts of Al II would argue against a low-metallicity white dwarf progenitor. We also show that this 6800 \AA\ feature is weakly present in other normal SN Ia events, and common in the SN 1991bg-like sub-class.

We present a search for an electromagnetic counterpart of the gravitational wave source GW151226. Using the Pan-STARRS1 telescope we mapped out 290 square degrees in the optical i_ps filter starting 11.5hr after the LIGO information release and lasting for a further 28 days. The first observations started 49.5hr after the time of the GW151226 detection. We typically reached sensitivity limits of i_ps = 20.3-20.8 and covered 26.5% of the LIGO probability skymap. We supplemented this with ATLAS survey data, reaching 31% of the probability region to shallower depths of m~19. We found 49 extragalactic transients (that are not obviously AGN), including a faint transient in a galaxy at 7Mpc (a luminous blue variable outburst) plus a rapidly decaying M-dwarf flare. Spectral classification of 20 other transient events showed them all to be supernovae. We found an unusual transient, PS15dpn, with an explosion date temporally coincident with GW151226 which evolved into a type Ibn supernova. The redshift of the transient is secure at z=0.1747 +/- 0.0001 and we find it unlikely to be linked, since the luminosity distance has a negligible probability of being consistent with that of GW151226. In the 290 square degrees surveyed we therefore do not find a likely counterpart. However we show that our survey strategy would be sensitive to NS-NS mergers producing kilonovae at D < 100 Mpc which is promising for future LIGO/Virgo searches.

We present an analysis of the early, rising light curves of 18 Type Ia supernovae (SNe Ia) discovered by the Palomar Transient Factory (PTF) and the La Silla-QUEST variability survey (LSQ). We fit these early data flux using a simple power-law $(f(t) = {\alpha\times t^n})$ to determine the time of first light $({t_0})$, and hence the rise-time $({t_{rise}})$ from first light to peak luminosity, and the exponent of the power-law rise ($n$). We find a mean uncorrected rise time of $18.98 {\pm} 0.54$ days, with individual SN rise-times ranging from $15.98$ to $24.7$ days. The exponent n shows significant departures from the simple 'fireball model' of $n = 2$ (or ${f(t) \propto t^2}$) usually assumed in the literature. With a mean value of $n = 2.44 {\pm} 0.13$, our data also show significant diversity from event to event. This deviation has implications for the distribution of 56Ni throughout the SN ejecta, with a higher index suggesting a lesser degree of 56Ni mixing. The range of n found also confirms that the 56Ni distribution is not standard throughout the population of SNe Ia, in agreement with earlier work measuring such abundances through spectral modelling. We also show that the duration of the very early light curve, before the luminosity has reached half of its maximal value, does not correlate with the light curve shape or stretch used to standardise SNe Ia in cosmological applications. This has implications for the cosmological fitting of SN Ia light curves.