We present DES16C3cje, a low-luminosity, long-lived type II supernova (SN II) at redshift 0.0618, detected by the Dark Energy Survey (DES). DES16C3cje is a unique SN. The spectra are characterized by extremely narrow photospheric lines corresponding to very low expansion velocities of $\lesssim1500$ km s$^{-1}$, and the light curve shows an initial peak that fades after 50 days before slowly rebrightening over a further 100 days to reach an absolute brightness of M$_r\sim -15.5$ mag. The decline rate of the late-time light curve is then slower than that expected from the powering by radioactive decay of $^{56}$Co but is comparable to that expected from accretion power. Comparing the bolometric light curve with hydrodynamical models, we find that DES16C3cje can be explained by either i) a low explosion energy (0.11 foe) and relatively large $^{56}$Ni production of 0.075 M$_{\odot}$ from a $\sim15$ M$_{\odot}$ red supergiant progenitor typical of other SNe II, or ii) a relatively compact $\sim40$ M$_{\odot}$ star, explosion energy of 1 foe, and 0.08 M$_{\odot}$ of $^{56}$Ni. Both scenarios require additional energy input to explain the late-time light curve, which is consistent with fallback accretion at a rate of $\sim0.5\times{10^{-8}}$ M$_{\odot}$ s$^{-1}$.

We present X-ray and radio observations of the Fast Blue Optical Transient (FBOT) CRTS-CSS161010 J045834-081803 (CSS161010 hereafter) at t=69-531 days. CSS161010 shows luminous X-ray ($L_x\sim5\times 10^{39}\,\rm{erg\,s^{-1}}$) and radio ($L_{\nu}\sim10^{29}\,\rm{erg\,s^{-1}Hz^{-1}}$) emission. The radio emission peaked at ~100 days post transient explosion and rapidly decayed. We interpret these observations in the context of synchrotron emission from an expanding blastwave. CSS161010 launched a mildly relativistic outflow with velocity $\Gamma\beta c\ge0.55c$ at ~100 days. This is faster than the non-relativistic AT2018cow ($\Gamma\beta c\sim0.1c$) and closer to ZTF18abvkwla ($\Gamma\beta c\ge0.3c$ at 63 days). The inferred initial kinetic energy of CSS161010 ($E_k\gtrsim10^{51}$ erg) is comparable to that of long Gamma Ray Bursts (GRBs), but the ejecta mass that is coupled to the mildly relativistic outflow is significantly larger ($\sim0.01-0.1\,\rm{M_{\odot}}$). This is consistent with the lack of observed gamma-rays. The luminous X-rays were produced by a different emission component to the synchrotron radio emission. CSS161010 is located at ~150 Mpc in a dwarf galaxy with stellar mass $M_{*}\sim10^{7}\,\rm{M_{\odot}}$ and specific star formation rate $sSFR\sim 0.3\,\rm{Gyr^{-1}}$. This mass is among the lowest inferred for host-galaxies of explosive transients from massive stars. Our observations of CSS161010 are consistent with an engine-driven aspherical explosion from a rare evolutionary path of a H-rich stellar progenitor, but we cannot rule out a stellar tidal disruption event on a centrally-located intermediate mass black hole. Regardless of the physical mechanism, CSS161010 establishes the existence of a new class of rare (rate $<0.4\%$ of the local core-collapse supernova rate) H-rich transients that can launch mildly relativistic outflows.

Rapidly evolving transients (RETs), also termed fast blue optical transients (FBOTs), are a distinct class of astrophysical event. They are characterised by lightcurves that decline much faster than common classes supernovae (SNe), span vast ranges in peak luminosity and can be seen to redshifts greater than 1. Their evolution on fast timescales has hindered high quality follow-up observations, and thus their origin and explosion/emission mechanism remains unexplained. In this paper we define the largest sample of RETs to date, comprising 106 objects from the Dark Energy Survey, and perform the most comprehensive analysis of RET host galaxies. Using deep-stacked photometry and emission-lines from OzDES spectroscopy, we derive stellar masses and star-formation rates (SFRs) for 49 host galaxies, and metallicities for 37. We find that RETs explode exclusively in star-forming galaxies and are thus likely associated with massive stars. Comparing RET hosts to samples of host galaxies of other explosive transients as well as field galaxies, we find that RETs prefer galaxies with high specific SFRs, indicating a link to young stellar populations, similar to stripped-envelope SNe. RET hosts appear to show a lack of chemical enrichment, their metallicities akin to long duration gamma-ray bursts and superluminous SN host galaxies. There are no clear relationships between properties of the host galaxies and the peak magnitudes or decline rates of the transients themselves.

We determine the viability of exploiting lensing time delays to observe strongly gravitationally lensed supernovae (gLSNe) from first light. Assuming a plausible discovery strategy, the Legacy Survey of Space and Time (LSST) and the Zwicky Transient Facility (ZTF) will discover $\sim$ 110 and $\sim$ 1 systems per year before the supernova (SN) explosion in the final image respectively. Systems will be identified $11.7^{+29.8}_{-9.3}$ days before the final explosion. We then explore the possibility of performing early-time observations for Type IIP and Type Ia SNe in LSST-discovered systems. Using a simulated Type IIP explosion, we predict that the shock breakout in one trailing image per year will peak at $\lesssim$ 24.1 mag ($\lesssim$ 23.3) in the $B$-band ($F218W$), however evolving over a timescale of $\sim$ 30 minutes. Using an analytic model of Type Ia companion interaction, we find that in the $B$-band we should observe at least one shock cooling emission event per year that peaks at $\lesssim$ 26.3 mag ($\lesssim$ 29.6) assuming all Type Ia gLSNe have a 1 M$_\odot$ red giant (main sequence) companion. We perform Bayesian analysis to investigate how well deep observations with 1 hour exposures on the European Extremely Large Telescope would discriminate between Type Ia progenitor populations. We find that if all Type Ia SNe evolved from the double-degenerate channel, then observations of the lack of early blue flux in 10 (50) trailing images would rule out more than 27% (19%) of the population having 1 M$_\odot$ main sequence companions at 95% confidence.

We present the first Hubble diagram of superluminous supernovae (SLSNe) out to a redshift of two, together with constraints on the matter density, $\Omega_{\rm M}$, and the dark energy equation-of-state parameter, $w(\equiv p/\rho)$. We build a sample of 20 cosmologically useful SLSNe~I based on light curve and spectroscopy quality cuts. We confirm the robustness of the peak decline SLSN~I standardization relation with a larger dataset and improved fitting techniques than previous works. We then solve the SLSN model based on the above standardisation via minimisation of the $\chi^2$ computed from a covariance matrix which includes statistical and systematic uncertainties. For a spatially flat $\Lambda$CDM cosmological model, we find $\Omega_{\rm M}=0.44^{+0.21}_{-0.21}$, with a rms of 0.28 mag for the residuals of the distance moduli. For an $w_0w_a$CDM cosmological model, the addition of SLSNe~I to a `baseline' measurement consisting of Planck temperature and WMAP polarization of the cosmic microwave background (CMB) fluctuation together with type Ia supernovae, results in a small improvement in the constraints of $w_0$ and $w_a$ of 4\%. We present simulations of future surveys with 847 SLSNe I and show that such a sample can deliver cosmological constraints in a flat $\Lambda$CDM model with the same precision (considering only statistical uncertainties) as current surveys that use type Ia supernovae, while providing an improvement of 15\% in the constraints on the time variation of dark energy, $w_0$ and $w_a$. This paper represents the proof-of-concept for superluminous supernova cosmology, and demonstrates they can provide an independent test of cosmology in the high-redshift ($z>1$) universe.

We present an analysis of DES17X1boj and DES16E2bjy, two peculiar transients discovered by the Dark Energy Survey (DES). They exhibit nearly identical double-peaked light curves which reach very different maximum luminosities (M$_\mathrm{r}$ = -15.4 and M$_\mathrm{r}$ = -17.9, respectively). The light curve evolution of these events is highly atypical and has not been reported before. The transients are found in different host environments: DES17X1boj was found near the nucleus of a spiral galaxy, while DES16E2bjy is located in the outskirts of a passive red galaxy. Early photometric data is well fitted with a blackbody and the resulting moderate photospheric expansion velocities (1800 km/s for DES17X1boj and 4800 km/s for DES16E2bjy) suggest an explosive or eruptive origin. Additionally, a feature identified as high-velocity CaII absorption (v $\approx$ 9400km/s) in the near-peak spectrum of DES17X1boj may imply that it is a supernova. While similar light curve evolution suggests a similar physical origin for these two transients, we are not able to identify or characterise the progenitors.

This work presents the observations and analysis of ATLAS19dqr/SN 2019bkc, an extraordinary rapidly evolving transient event located in an isolated environment, tens of kiloparsecs from any likely host. Its light curves rise to maximum light in $5-6$ d and then display a decline of $\Delta m_{15} \sim5$ mag. With such a pronounced decay, it has one of the most rapidly evolving light curves known for a stellar explosion. The early spectra show similarities to normal and `ultra-stripped' type Ic SNe, but the early nebular phase spectra, which were reached just over two weeks after explosion, display prominent calcium lines, marking SN 2019bkc as a Ca-rich transient. The Ca emission lines at this phase show an unprecedented and unexplained blueshift of 10,000 -- 12,000 km/s. Modelling of the light curve and the early spectra suggests that the transient had a low ejecta mass of $0.2 - 0.4$ M$_\odot$ and a low kinetic energy of $ (2-4)\times 10^{50}$ erg, giving a specific kinetic energy $\sim1$ [$10^{51}$ erg]/M$_\odot$. The origin of this event cannot be unambiguously defined. While the abundance distribution used to model the spectra marginally favours a progenitor of white dwarf origin through the tentative identification of \ArII, the specific kinetic energy, which is defined by the explosion mechanism, is found to be more similar to an ultra-stripped core-collapse events. SN 2019bkc adds to the diverse range of physical properties shown by Ca-rich events.

We present improved photometric measurements for the host galaxies of 206 spectroscopically confirmed type Ia supernovae discovered by the Dark Energy Survey Supernova Program (DES-SN) and used in the first DES-SN cosmological analysis. Fitting spectral energy distributions to the $griz$ photometric measurements of the DES-SN host galaxies, we derive stellar masses and star-formation rates. For the DES-SN sample, when considering a 5D ($z$, $x_1$, $c$, $\alpha$, $\beta$) bias correction, we find evidence of a Hubble residual `mass step', where SNe Ia in high mass galaxies ($>10^{10} \textrm{M}_{\odot}$) are intrinsically more luminous (after correction) than their low mass counterparts by $\gamma=0.040\pm0.019$mag. This value is larger by $0.031$mag than the value found in the first DES-SN cosmological analysis. This difference is due to a combination of updated photometric measurements and improved star formation histories and is not from host-galaxy misidentification. When using a 1D (redshift-only) bias correction the inferred mass step is larger, with $\gamma=0.066\pm0.020$mag. The 1D-5D $\gamma$ difference for DES-SN is $0.026\pm0.009$mag. We show that this difference is due to a strong correlation between host galaxy stellar mass and the $x_1$ component of the 5D distance-bias correction. To better understand this effect, we include an intrinsic correlation between light-curve width and stellar mass in simulated SN Ia samples. We show that a 5D fit recovers $\gamma$ with $-9$mmag bias compared to a $+2$mmag bias for a 1D fit. This difference can explain part of the discrepancy seen in the data. Improvements in modeling correlations between galaxy properties and SN is necessary to determine the implications for $\gamma$ and ensure unbiased precision estimates of the dark energy equation-of-state as we enter the era of LSST.

The five-year Dark Energy Survey supernova programme (DES-SN) is one of the largest and deepest transient surveys to date in terms of volume and number of supernovae. Identifying and characterising the host galaxies of transients plays a key role in their classification, the study of their formation mechanisms, and the cosmological analyses. To derive accurate host galaxy properties, we create depth-optimised coadds using single-epoch DES-SN images that are selected based on sky and atmospheric conditions. For each of the five DES-SN seasons, a separate coadd is made from the other 4 seasons such that each SN has a corresponding deep coadd with no contaminating SN emission. The coadds reach limiting magnitudes of order $\sim 27$ in $g$-band, and have a much smaller magnitude uncertainty than the previous DES-SN host templates, particularly for faint objects. We present the resulting multi-band photometry of host galaxies for samples of spectroscopically confirmed type Ia (SNe Ia), core-collapse (CCSNe), and superluminous (SLSNe) as well as rapidly evolving transients (RETs) discovered by DES-SN. We derive host galaxy stellar masses and probabilistically compare stellar-mass distributions to samples from other surveys. We find that the DES spectroscopically confirmed sample of SNe Ia selects preferentially fewer high mass hosts at high redshift compared to other surveys, while at low redshift the distributions are consistent. DES CCSNe and SLSNe hosts are similar to other samples, while RET hosts are unlike the hosts of any other transients, although these differences have not been disentangled from selection effects.

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.

We present a sample of 21 hydrogen-free superluminous supernovae (SLSNe-I), and one hydrogen-rich SLSN (SLSN-II) detected during the five-year Dark Energy Survey (DES). These SNe, located in the redshift range 0.220<z<1.998, represent the largest homogeneously-selected sample of SLSN events at high redshift. We present the observed g,r, i, z light curves for these SNe, which we interpolate using Gaussian Processes. The resulting light curves are analysed to determine the luminosity function of SLSN-I, and their evolutionary timescales. The DES SLSN-I sample significantly broadens the distribution of SLSN-I light curve properties when combined with existing samples from the literature. We fit a magnetar model to our SLSNe, and find that this model alone is unable to replicate the behaviour of many of the bolometric light curves. We search the DES SLSN-I light curves for the presence of initial peaks prior to the main light-curve peak. Using a shock breakout model, our Monte Carlo search finds that 3 of our 14 events with pre-max data display such initial peaks. However, 10 events show no evidence for such peaks, in some cases down to an absolute magnitude of <-16, suggesting that such features are not ubiquitous to all SLSN-I events. We also identify a red pre-peak feature within the light curve of one SLSN, which is comparable to that observed within SN2018bsz.

We present the volumetric rate of normal type Ia supernovae (SNe Ia) discovered by the Palomar Transient Factory (PTF). Using strict data-quality cuts, and considering only periods when the PTF maintained a regular cadence, PTF discovered 90 SNe Ia at $z\le0.09$ in a well-controlled sample over three years of operation (2010-2012). We use this to calculate the volumetric rate of SN Ia events by comparing this sample to simulations of hundreds of millions of SN Ia light curves produced in statistically representative realisations of the PTF survey. This quantifies the recovery efficiency of each PTF SN Ia event, and thus the relative weighting of each event. From this, the volumetric SN Ia rate was found to be $r_v=2.43\pm0.29\,\text{(stat)}_{-0.19}^{+0.33}\text{(sys)}\times10^{-5}\,\text{SNe yr}^{-1}\,\text{Mpc}^{-3}\, h_{70}^{3}$. This represents the most precise local measurement of the SN Ia rate. We fit a simple SN Ia delay-time distribution model, $\propto\mathrm{t}^{-\beta}$, to our PTF rate measurement combined with a literature sample of rate measurements from surveys at higher-redshifts. We find $\beta{\sim}1$, consistent with a progenitor channel governed by the gravitational in-spiral of binary white dwarfs.

We present an exquisite, 30-min cadence Kepler (K2) light curve of the Type Ia supernova (SN Ia) 2018oh (ASASSN-18bt), starting weeks before explosion, covering the moment of explosion and the subsequent rise, and continuing past peak brightness. These data are supplemented by multi-color Pan-STARRS1 and CTIO 4-m DECam observations obtained within hours of explosion. The K2 light curve has an unusual two-component shape, where the flux rises with a steep linear gradient for the first few days, followed by a quadratic rise as seen for typical SNe Ia. This "flux excess" relative to canonical SN Ia behavior is confirmed in our $i$-band light curve, and furthermore, SN 2018oh is especially blue during the early epochs. The flux excess peaks 2.14$\pm0.04$ days after explosion, has a FWHM of 3.12$\pm0.04$ days, a blackbody temperature of $T=17,500^{+11,500}_{-9,000}$ K, a peak luminosity of $4.3\pm0.2\times10^{37}\,{\rm erg\,s^{-1}}$, and a total integrated energy of $1.27\pm0.01\times10^{43}\,{\rm erg}$. We compare SN 2018oh to several models that may provide additional heating at early times, including collision with a companion and a shallow concentration of radioactive nickel. While all of these models generally reproduce the early K2 light curve shape, we slightly favor a companion interaction, at a distance of $\sim$$2\times10^{12}\,{\rm cm}$ based on our early color measurements, although the exact distance depends on the uncertain viewing angle. Additional confirmation of a companion interaction in future modeling and observations of SN 2018oh would provide strong support for a single-degenerate progenitor system.

We present a measurement of the volumetric rate of `calcium-rich' optical transients in the local universe, using a sample of three events from the Palomar Transient Factory (PTF). This measurement builds on a detailed study of the PTF transient detection efficiencies, and uses a Monte Carlo simulation of the PTF survey. We measure the volumetric rate of calcium-rich transients to be higher than previous estimates: $1.21^{+1.13}_{-0.39}\times10^{-5}$ events yr$^{-1}$ Mpc$^{-3}$. This is equivalent to 33-94% of the local volumetric type Ia supernova rate. This calcium-rich transient rate is sufficient to reproduce the observed calcium abundances in galaxy clusters, assuming an asymptotic calcium yield per calcium-rich event of ~0.05$\mathrm{M}_{\odot}$. We also study the PTF detection efficiency of these transients as a function of position within their candidate host galaxies. We confirm as a real physical effect previous results that suggest calcium-rich transients prefer large physical offsets from their host galaxies.

The late-time spectra of Type Ia supernovae (SNe Ia) are powerful probes of the underlying physics of their explosions. We investigate the late-time optical and near-infrared spectra of seven SNe Ia obtained at the VLT with XShooter at $>$200 d after explosion. At these epochs, the inner Fe-rich ejecta can be studied. We use a line-fitting analysis to determine the relative line fluxes, velocity shifts, and line widths of prominent features contributing to the spectra ([Fe II], [Ni II], and [Co III]). By focussing on [Fe II] and [Ni II] emission lines in the ~7000-7500 \AA\ region of the spectrum, we find that the ratio of stable [Ni II] to mainly radioactively-produced [Fe II] for most SNe Ia in the sample is consistent with Chandrasekhar-mass delayed-detonation explosion models, as well as sub-Chandrasekhar mass explosions that have metallicity values above solar. The mean measured Ni/Fe abundance of our sample is consistent with the solar value. The more highly ionised [Co III] emission lines are found to be more centrally located in the ejecta and have broader lines than the [Fe II] and [Ni II] features. Our analysis also strengthens previous results that SNe Ia with higher Si II velocities at maximum light preferentially display blueshifted [Fe II] 7155 \AA\ lines at late times. Our combined results lead us to speculate that the majority of normal SN Ia explosions produce ejecta distributions that deviate significantly from spherical symmetry.

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.

Merging neutron stars offer an exquisite laboratory for simultaneously studying strong-field gravity and matter in extreme environments. We establish the physical association of an electromagnetic counterpart EM170817 to gravitational waves (GW170817) detected from merging neutron stars. By synthesizing a panchromatic dataset, we demonstrate that merging neutron stars are a long-sought production site forging heavy elements by r-process nucleosynthesis. The weak gamma-rays seen in EM170817 are dissimilar to classical short gamma-ray bursts with ultra-relativistic jets. Instead, we suggest that breakout of a wide-angle, mildly-relativistic cocoon engulfing the jet elegantly explains the low-luminosity gamma-rays, the high-luminosity ultraviolet-optical-infrared and the delayed radio/X-ray emission. We posit that all merging neutron stars may lead to a wide-angle cocoon breakout; sometimes accompanied by a successful jet and sometimes a choked jet.

We present the transient source detection efficiencies of the Palomar Transient Factory (PTF), parameterizing the number of transients that PTF found, versus the number of similar transients that occurred over the same period in the survey search area but that were missed. PTF was an optical sky survey carried out with the Palomar 48-inch telescope over 2009-2012, observing more than 8000 square degrees of sky with cadences of between 1 and 5 days, locating around 50,000 non-moving transient sources, and spectroscopically confirming around 1900 supernovae. We assess the effectiveness with which PTF detected transient sources, by inserting ~7 million artificial point sources into real PTF data. We then study the efficiency with which the PTF real-time pipeline recovered these sources as a function of the source magnitude, host galaxy surface brightness, and various observing conditions (using proxies for seeing, sky brightness, and transparency). The product of this study is a multi-dimensional recovery efficiency grid appropriate for the range of observing conditions that PTF experienced, and that can then be used for studies of the rates, environments, and luminosity functions of different transient types using detailed Monte Carlo simulations. We illustrate the technique using the observationally well-understood class of type Ia supernovae.

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.