Low-luminosity type II supernovae (LL SNe~II) make up the low explosion energy end of core-collapse SNe, but their study and physical understanding remain limited. We present SN\,2016aqf, a LL SN~II with extensive spectral and photometric coverage. We measure a $V$-band peak magnitude of $-14.58$\,mag, a plateau duration of $\sim$100\,days, and an inferred $^{56}$Ni mass of $0.008 \pm 0.002$\,\msun. The peak bolometric luminosity, L$_{\rm bol} \approx 10^{41.4}$\,erg\,s$^{-1}$, and its spectral evolution is typical of other SNe in the class. Using our late-time spectra, we measure the [\ion{O}{i}] $\lambda\lambda6300, 6364$ lines, which we compare against SN II spectral synthesis models to constrain the progenitor zero-age main-sequence mass. We find this to be 12 $\pm$ 3\,\msun. Our extensive late-time spectral coverage of the [\ion{Fe}{ii}] $\lambda7155$ and [\ion{Ni}{ii}] $\lambda7378$ lines permits a measurement of the Ni/Fe abundance ratio, a parameter sensitive to the inner progenitor structure and explosion mechanism dynamics. We measure a constant abundance ratio evolution of $0.081^{+0.009}_{-0.010}$, and argue that the best epochs to measure the ratio are at $\sim$200 -- 300\,days after explosion. We place this measurement in the context of a large sample of SNe II and compare against various physical, light-curve and spectral parameters, in search of trends which might allow indirect ways of constraining this ratio. We do not find correlations predicted by theoretical models; however, this may be the result of the exact choice of parameters and explosion mechanism in the models, the simplicity of them and/or primordial contamination in the measured abundance ratio.

On 2019 August 14, the LIGO and Virgo interferometers detected a high-significance event labelled S190814bv. Preliminary analysis of the GW data suggests that the event was likely due to the merger of a compact binary system formed by a BH and a NS. ElectromagNetic counterparts of GRAvitational wave sources at the VEry Large Telescope (ENGRAVE) collaboration members carried out an intensive multi-epoch, multi-instrument observational campaign to identify the possible optical/near infrared counterpart of the event. In addition, the ATLAS, GOTO, GRAWITA-VST, Pan-STARRS and VINROUGE projects also carried out a search on this event. Our observations allow us to place limits on the presence of any counterpart and discuss the implications for the kilonova (KN) possibly generated by this NS-BH merger, and for the strategy of future searches. Altogether, our observations allow us to exclude a KN with large ejecta mass $M\gtrsim 0.1\,\mathrm{M_\odot}$ to a high ($>90\%$) confidence, and we can exclude much smaller masses in a subsample of our observations. This disfavours the tidal disruption of the neutron star during the merger. Despite the sensitive instruments involved in the campaign, given the distance of S190814bv we could not reach sufficiently deep limits to constrain a KN comparable in luminosity to AT 2017gfo on a large fraction of the localisation probability. This suggests that future (likely common) events at a few hundreds Mpc will be detected only by large facilities with both high sensitivity and large field of view. Galaxy-targeted observations can reach the needed depth over a relevant portion of the localisation probability with a smaller investment of resources, but the number of galaxies to be targeted in order to get a fairly complete coverage is large, even in the case of a localisation as good as that of this event.

At 66 Mpc, AT2019qiz is the closest optical tidal disruption event (TDE) to date, with a luminosity intermediate between the bulk of the population and the faint-and-fast event iPTF16fnl. Its proximity allowed a very early detection and triggering of multiwavelength and spectroscopic follow-up well before maximum light. The velocity dispersion of the host galaxy and fits to the TDE light curve indicate a black hole mass $\approx 10^6$ M $_\odot$, disrupting a star of $\approx$ 1 M$_\odot$. By analysing our comprehensive UV, optical and X-ray data, we show that the early optical emission is dominated by an outflow, with a luminosity evolution $L\propto t^2$, consistent with a photosphere expanding at constant velocity ($\gtrsim 2000$ km/s), and a line-forming region producing initially blueshifted H and He II profiles with $v=3000-10000$ km/s. The fastest optical ejecta approach the velocity inferred from radio detections, thus the same outflow is likely responsible for both the fast optical rise and the radio emission -- the first time this connection has been determined in a TDE. The light curve rise begins $35\pm1.5$ days before maximum light, peaking when AT2019qiz reaches the radius where optical photons can escape. The photosphere then undergoes a sudden transition, first cooling at constant radius then contracting at constant temperature. At the same time, the blueshifts disappear from the spectrum and Bowen fluorescence lines (N III) become prominent, implying a source of far-UV photons, while the X-ray light curve peaks at $\approx 10^{41}$ erg/s. Thus accretion began promptly in this event, favouring accretion-powered over collision-powered outflow models. The size and mass of the outflow are consistent with the reprocessing layer needed to explain the large optical to X-ray ratio in this and other optical TDEs.

Comparison of theoretical line profiles to observations provides important tests for supernova explosion models. We study the shapes of radioactive decay lines predicted by current 3D core-collapse explosion simulations, and compare these to observations of SN 1987A and Cas A. Both the widths and shifts of decay lines vary by several thousand kilometers per second depending on viewing angle. The line profiles can be complex with multiple peaks. By combining observational constraints from 56Co decay lines, 44Ti decay lines, and Fe IR lines, we delineate a picture of the morphology of the explosive burning ashes in SN 1987A. For M_ZAMS=15-20 Msun progenitors exploding with ~1.5 *10^51 erg, ejecta structures suitable to reproduce the observations involve a bulk asymmetry of the 56Ni of at least ~400 km/s and a bulk velocity of at least ~1500 km/s. By adding constraints to reproduce the UVOIR bolometric light curve of SN 1987A up to 600d, an ejecta mass around 14 Msun is favoured. We also investigate whether observed decay lines can constrain the neutron star (NS) kick velocity. The model grid provides a constraint V_NS > V_redshift, and applying this to SN 1987A gives a NS kick of at least 500 km/s. For Cas A, our single model provides a satisfactory fit to the NuSTAR observations and reinforces the result that current neutrino-driven core-collapse SN models can achieve enough bulk asymmetry in the explosive burning material. Finally, we investigate the internal gamma-ray field and energy deposition, and compare the 3D models to 1D approximations.

Superluminous supernovae radiate up to 100 times more energy than normal supernovae. The origin of this energy and the nature of their stellar progenitors are poorly understood. We identify neutral iron lines in the spectrum of one such transient, SN 2006gy, and show that they require a large mass of iron (>~0.3 Msun) expanding at 1500 km/s. We demonstrate that a model of a standard Type Ia supernova hitting a shell of circumstellar material produces a light curve and late-time iron-dominated spectrum that match SN 2006gy. In such a scenario, common envelope evolution of the progenitor system can synchronize envelope ejection and supernova explosion and may explain these bright transients.

We present SN2018kzr, the fastest declining supernova-like transient, second only to the kilonova, AT2017gfo. SN2018kzr is characterized by a peak magnitude of $M_r = -17.98$, peak bolometric luminosity of ${\sim} 1.4 \times 10^{43}$erg s$^{\mathrm{-1}}$ and a rapid decline rate of $0.48 \pm 0.03$ mag day$^{\textrm{-1}}$ in the $r$ band. The bolometric luminosity evolves too quickly to be explained by pure $^{\mathrm{56}}$Ni heating, necessitating the inclusion of an alternative powering source. Incorporating the spin-down of a magnetized neutron star adequately describes the lightcurve and we estimate a small ejecta mass of $M_\mathrm{ej} = 0.10 \pm 0.05$ $\textrm{M}_{\odot}$. Our spectral modelling suggests the ejecta is composed of intermediate mass elements including O, Si and Mg and trace amounts of Fe-peak elements, which disfavours a binary neutron star merger. We discuss three explosion scenarios for SN2018kzr, given the low ejecta mass, intermediate mass element composition and the high likelihood of additional powering - core collapse of an ultra-stripped progenitor, the accretion induced collapse of a white dwarf and the merger of a white dwarf and neutron star. The requirement for an alternative input energy source favours either the accretion induced collapse with magnetar powering or a white dwarf - neutron star merger with energy from disk wind shocks.

During the first few hundred days after the explosion, core-collapse supernovae (SNe) emit down-scattered X-rays and gamma-rays originating from radioactive line emissions, primarily from the $^{56}$Ni $\rightarrow$ $^{56}$Co $\rightarrow$ $^{56}$Fe chain. We use SN models based on three-dimensional neutrino-driven explosion simulations of single stars and mergers to compute this emission and compare the predictions with observations of SN 1987A. A number of models are clearly excluded, showing that high-energy emission is a powerful way of discriminating between models. The best models are almost consistent with the observations, but differences that cannot be matched by a suitable choice of viewing angle are evident. Therefore, our self-consistent models suggest that neutrino-driven explosions are able to produce, in principle, sufficient mixing, although remaining discrepancies may require small changes to the progenitor structures. The soft X-ray cutoff is primarily determined by the metallicity of the progenitor envelope. The main effect of asymmetries is to vary the flux level by a factor of ${\sim}$3. For the more asymmetric models, the shapes of the light curves also change. In addition to the models of SN 1987A, we investigate two models of Type II-P SNe and one model of a stripped-envelope Type IIb SN. The Type II-P models have similar observables as the models of SN 1987A, but the stripped-envelope SN model is significantly more luminous and evolves faster. Finally, we make simple predictions for future observations of nearby SNe.

Gamma-ray astronomy explores the most energetic photons in nature to address some of the most pressing puzzles in contemporary astrophysics. It encompasses a wide range of objects and phenomena: stars, supernovae, novae, neutron stars, stellar-mass black holes, nucleosynthesis, the interstellar medium, cosmic rays and relativistic-particle acceleration, and the evolution of galaxies. MeV gamma-rays provide a unique probe of nuclear processes in astronomy, directly measuring radioactive decay, nuclear de-excitation, and positron annihilation. The substantial information carried by gamma-ray photons allows us to see deeper into these objects, the bulk of the power is often emitted at gamma-ray energies, and radioactivity provides a natural physical clock that adds unique information. New science will be driven by time-domain population studies at gamma-ray energies. This science is enabled by next-generation gamma-ray instruments with one to two orders of magnitude better sensitivity, larger sky coverage, and faster cadence than all previous gamma-ray instruments. This transformative capability permits: (a) the accurate identification of the gamma-ray emitting objects and correlations with observations taken at other wavelengths and with other messengers; (b) construction of new gamma-ray maps of the Milky Way and other nearby galaxies where extended regions are distinguished from point sources; and (c) considerable serendipitous science of scarce events -- nearby neutron star mergers, for example. Advances in technology push the performance of new gamma-ray instruments to address a wide set of astrophysical questions.

Hydrogen-rich, core-collapse supernovae are typically divided into four classes: IIP, IIL, IIn, and IIb. In general, interaction with circumstellar material is only considered for Type IIn supernovae. However, recent hydrodynamic modeling of IIP and IIL supernovae requires circumstellar material to reproduce their early light curves. In this scenario, IIL supernovae experience large amounts of mass loss before exploding. We test this hypothesis on ASASSN-15oz, a Type IIL supernova. With extensive follow-up in the X- ray, UV, optical, IR, and radio we present our search for signs of interaction, and the mass-loss history indicated by their detection. We find evidence of short-lived intense mass-loss just prior to explosion from light curve modeling, amounting in 1.5 M$_{\odot}$ of material within 1800 R$_{\odot}$ of the progenitor. We also detect the supernova in the radio, indicating mass-loss rates of $10^{-6}-10^{-7}$ M$_{\odot}$ yr$^{-1}$ prior to the extreme mass-loss period. Our failure to detect the supernova in the X-ray and the lack of narrow emission lines in the UV, optical, and NIR do not contradict this picture and place an upper limit on the mass-loss rate outside the extreme period of $<10^{-4}$ M$_{\odot}$ yr$^{-1}$. This paper highlights the importance gathering comprehensive data on more Type II supernovae to enable detailed modeling of the progenitor and supernova which can elucidate their mass-loss histories and envelope structures and thus inform stellar evolution models.

We present observations of supernova (SN) 2017ens, discovered by the ATLAS survey and identified as a hot blue object through the GREAT program. The redshift z=0.1086 implies a peak brightness of M_g=-21.1 mag, placing the object within the regime of superluminous supernovae. We observe a dramatic spectral evolution, from initially being blue and featureless, to later developing features similar to those of the broadlined Type Ic SN 1998bw, and finally showing ~2000 km s^-1 wide H-alpha and H-beta emission. Relatively narrow Balmer emission (reminiscent of a SN IIn) is present at all times. We also detect coronal lines, indicative of a dense circumstellar medium. We constrain the progenitor wind velocity to ~50-60 km s^-1 based on P-Cygni profiles, which is far slower than those present in Wolf-Rayet stars. This may suggest that the progenitor passed through a luminous blue variable phase, or that the wind is instead from a binary companion red supergiant star. At late times we see the ~2000 km s^-1 wide H-alpha emission persisting at high luminosity (~3x10^40 erg s^-1) for at least 100 day, perhaps indicative of additional mass loss at high velocities that could have been ejected by a pulsational pair instability.

The detection of the binary neutron star GW170817 together with the observation of electromagnetic counterparts across the entire spectrum inaugurated a new era of multi-messenger astronomy. In this study we incorporate wavelength-dependent opacities and emissivities calculated from atomic-structure data enabling us to model both the measured lightcurves and spectra of the electromagnetic transient AT2017gfo. Best-fits of the observational data are obtained by Gaussian Process Regression, which allows us to present posterior samples for the kilonova and source properties connected to GW170817. Incorporating constraints obtained from the gravitational wave signal obtained by the LIGO-Virgo Scientific Collaboration, we present a $90\%$ upper bound on the mass ratio $q \lesssim 1.38$ and a lower bound on the tidal deformability of $\tilde{\Lambda} \gtrsim 197$, which rules out sufficiently soft equations of state. Our analysis is a path-finder for more realistic kilonova models and shows how the combination of gravitational wave and electromagnetic measurements allow for stringent constraints on the source parameters and the supranuclear equation of state.

We present JEKYLL, a new code for modelling of supernova (SN) spectra and lightcurves based on Monte-Carlo (MC) techniques for the radiative transfer. The code assumes spherical symmetry, homologous expansion and steady state for the matter, but is otherwise capable of solving the time-dependent radiative transfer problem in non-local-thermodynamic-equilibrium (NLTE). The method used was introduced in a series of papers by Lucy, but the full time-dependent NLTE capabilities of it have never been tested. Here, we have extended the method to include non-thermal excitation and ionization as well as charge-transfer and two-photon processes. Based on earlier work, the non-thermal rates are calculated by solving the Spencer-Fano equation. Using a method previously developed for the SUMO code, macroscopic mixing of the material is taken into account in a statistical sense. In addition, a statistical Markov-chain model is used to sample the emission frequency, and we introduce a method to control the sampling of the radiation field. Except for a description of JEKYLL, we provide comparisons with the ARTIS, SUMO and CMFGEN codes, which show good agreement in the calculated spectra as well as the state of the gas. In particular, the comparison with CMFGEN, which is similar in terms of physics but uses a different technique, shows that the Lucy method does indeed converge in the time-dependent NLTE case. Finally, as an example of the time-dependent NLTE capabilities of JEKYLL, we present a model of a Type IIb SN, taken from a set of models presented and discussed in detail in an accompanying paper. Based on this model we investigate the effects of NLTE, in particular those arising from non-thermal excitation and ionization, and find strong effects even on the bolometric lightcurve. This highlights the need for full NLTE calculations when simulating the spectra and lightcurves of SNe.

Despite more than 30 years of searches, the compact object in Supernova (SN) 1987A has not yet been detected. We present new limits on the compact object in SN 1987A using millimeter, near-infrared, optical, ultraviolet, and X-ray observations from ALMA, VLT, HST, and Chandra. The limits are approximately 0.1 mJy ($0.1\times 10^{-26}$ erg s$^{-1}$ cm$^{-2}$ Hz$^{-1}$) at 213 GHz, 1 Lsun ($6\times 10^{-29}$ erg s$^{-1}$ cm$^{-2}$ Hz$^{-1}$) in optical if our line-of-sight is free of ejecta dust, and $10^{36}$ erg s$^{-1}$ ($2\times 10^{-30}$ erg s$^{-1}$ cm$^{-2}$ Hz$^{-1}$) in 2-10 keV X-rays. Our X-ray limits are an order of magnitude less constraining than previous limits because we use a more realistic ejecta absorption model based on three-dimensional neutrino-driven SN explosion models (presented in an accompanying article). The allowed bolometric luminosity of the compact object is 22 Lsun if our line-of-sight is free of ejecta dust, or 138 Lsun if dust-obscured. Depending on assumptions, these values limit the effective temperature of a neutron star to <4-8 MK and do not exclude models, which typically are in the range 3-4 MK. For the simplest accretion model, the accretion rate for an efficiency $\eta$ is limited to $< 10^{-11} \eta^{-1}$ Msun yr$^{-1}$, which excludes most predictions. For pulsar activity modeled by a rotating magnetic dipole in vacuum, the limit on the magnetic field strength ($B$) for a given spin period ($P$) is $B < 10^{14} P^2$ G s$^{-2}$. By combining information about radiation reprocessing and geometry, it is likely that the compact object is a dust-obscured thermally-emitting neutron star, which may appear as a region of higher-temperature ejecta dust emission.

While interaction with circumstellar material is known to play an important role in Type IIn supernovae (SNe), analyses of the more common SNe IIP and IIL have not traditionally included interaction as a significant power source. However, recent campaigns to observe SNe within days of explosion have revealed narrow emission lines of high-ionization species in the earliest spectra of luminous SNe II of all subclasses. These "flash ionization" features indicate the presence of a confined shell of material around the progenitor star. Here we present the first low-luminosity (LL) SN to show flash ionization features, SN 2016bkv. This SN peaked at $M_V = -16$ mag and has H{\alpha} expansion velocities under 1350 km/s around maximum light, placing it at the faint/slow end of the distribution of SNe IIP (similar to SN 2005cs). The light-curve shape of SN 2016bkv is also extreme among SNe IIP. A very strong initial peak could indicate additional luminosity from circumstellar interaction. A very small fall from the plateau to the nickel tail indicates unusually large production of radioactive nickel compared to other LL SNe IIP. A comparison between nebular spectra of SN 2016bkv and models raises the possibility that SN 2016bkv is an electron-capture supernova.

Red supergiants have been confirmed as the progenitor stars of the majority of hydrogen-rich type II supernovae. However, while such stars are observed with masses >25M$_\odot$, detections of >18M$_\odot$ progenitors remain elusive. Red supergiants are also expected to form at all metallicities, but discoveries of explosions from low-metallicity progenitors are scarce. Here, we report observations of the type II supernova, SN 2015bs, for which we infer a progenitor metallicity of $\leq$0.1Z$_\odot$ from comparison to photospheric-phase spectral models, and a Zero Age Main-Sequence mass of 17-25M$_\odot$ through comparison to nebular-phase spectral models. SN 2015bs displays a normal 'plateau' light-curve morphology, and typical spectral properties, implying a red supergiant progenitor. This is the first example of such a high mass progenitor for a 'normal' type II supernova, suggesting a link between high mass red supergiant explosions and low-metallicity progenitors.

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.

We present one year of optical and near-infrared photometry and spectroscopy of the Type IIP SNe 2013K and 2013am. Both objects are affected by significant extinction, due to their location in dusty regions of their respective host galaxies, ESO 009-10 and NGC 3623 (M65). From the photospheric to nebular phases, these objects display spectra congruent with those of underluminous Type IIP SNe (i.e. the archetypal SNe 1997D or 2005cs), showing low photospheric velocities (~2 X 10**3 km/s at 50d) together with features arising from Ba II which are particularly prominent in faint SNe IIP. The peak V-band magnitudes of SN 2013K (-15.6 mag) and SN 2013am (-16.2 mag) are fainter than standard-luminosity Type IIP SNe. The ejected Nickel masses are 0.012+-0.010 and 0.015+-0.006 Msol for SN 2013K and SN 2013am, respectively. The physical properties of the progenitors at the time of explosion are derived through hydrodynamical modelling. Fitting the bolometric curves, the expansion velocity and the temperature evolution, we infer total ejected masses of 12 and 11.5 Msol, pre-SN radii of ~460 and ~360 Rsol, and explosion energies of 0.34 foe and 0.40 foe for SN 2013K and SN 2013am. Late time spectra are used to estimate the progenitor masses from the strength of nebular emission lines, which turn out to be consistent with red supergiant progenitors of ~15 Msol. For both SNe, a low-energy explosion of a moderate-mass red supergiant star is therefore the favoured scenario.

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.

A large fraction of core-collapse supernovae (CCSNe), 30-50%, are expected to originate from the low-mass end of progenitors with $M_{\rm ZAMS}~= 8-12~M_\odot$. However, degeneracy effects make stellar evolution modelling of such stars challenging, and few predictions for their supernova light curves and spectra have been presented. Here we calculate synthetic nebular spectra of a 9 $M_\odot$ Fe CCSN model exploded with the neutrino mechanism. The model predicts emission lines with FWHM$\sim$1000 km/s, including signatures from each deep layer in the metal core. We compare this model to observations of the three subluminous IIP SNe with published nebular spectra; SN 1997D, SN 2005cs, and SN 2008bk. The prediction of both line profiles and luminosities are in good agreement with SN 1997D and SN 2008bk. The close fit of a model with no tuning parameters provides strong evidence for an association of these objects with low-mass Fe CCSNe. For SN 2005cs, the interpretation is less clear, as the observational coverage ended before key diagnostic lines from the core had emerged. We perform a parameterised study of the amount of explosively made stable nickel, and find that none of these three SNe show the high $^{58}$Ni/$^{56}$Ni ratio predicted by current models of electron capture SNe (ECSNe) and ECSN-like explosions. Combined with clear detection of lines from O and He shell material, these SNe rather originate from Fe core progenitors. We argue that the outcome of self-consistent explosion simulations of low-mass stars, which gives fits to many key observables, strongly suggests that the class of subluminous Type IIP SNe is the observational counterpart of the lowest mass CCSNe.

We present our study of OGLE-2014-SN-073, one of the brightest Type II SN ever discovered, with an unusually broad lightcurve combined with high ejecta velocities. From our hydrodynamical modelling we infer a remarkable ejecta mass of $60^{+42}_{-16}$~M$_\odot$, and a relatively high explosion energy of $12.4^{+13.0}_{-5.9} \times10^{51}$~erg. We show that this object belongs, with a very small number of other hydrogen-rich SNe, to an energy regime that is not explained by standard core-collapse (CC) neutrino-driven explosions. We compare the quantities inferred by the hydrodynamical modelling with the expectations of various exploding scenarios, trying to explain the high energy and luminosity released. We find some qualitative similarities with pair-instabilities SNe, although a prompt injection of energy by a magnetar seems also a viable alternative to explain such extreme event.

Type Ia supernovae are widely accepted to be the outcomes of thermonuclear explosions in white dwarf stars. However, many details of these explosions remain uncertain (e.g. the mass, ignition mechanism, and flame speed). Theory predicts that at very late times (beyond 1000 d) it might be possible to distinguish between explosion models. Few very nearby supernovae can be observed that long after the explosion. The Type Ia supernova SN 2011fe located in M101 and along a line of sight with negligible extinction, provides us with the once-in-a-lifetime chance to obtain measurements that may distinguish between theoretical models. In this work, we present the analysis of photometric data of SN 2011fe taken between 900 and 1600 days after explosion with Gemini and HST. At these extremely late epochs theory suggests that the light curve shape might be used to measure isotopic abundances which is a useful model discriminant. However, we show in this work that there are several currently not well constrained physical processes introducing large systematic uncertainties to the isotopic abundance measurement. We conclude that without further detailed knowledge of the physical processes at this late stage one cannot reliably exclude any models on the basis of this dataset.

Results of calculations for inelastic e+Mg effective collision strengths for the lowest 25 physical states of Mg I (up to 3s6p 1P), and thus 300 transitions, from the convergent close-coupling (CCC) and the B-spline R-matrix (BSR) methods are presented. At temperatures of interest, ~5000 K, the results of the two calculations differ on average by only 4%, with a scatter of 27%. As the methods are independent, this suggests that the calculations provide datasets for e+Mg collisions accurate to this level. Comparison with the commonly used dataset compiled by Mauas et al. (1988), covering 25 transitions among 12 states, suggests the Mauas et al. data are on average ~57% too low, and with a very large scatter of a factor of ~6.5. In particular the collision strength for the transition corresponding to the Mg I intercombination line at 457 nm is significantly underestimated by Mauas et al., which has consequences for models that employ this dataset. In giant stars the new data leads to a stronger line compared to previous non-LTE calculations, and thus a reduction in the non-LTE abundance correction by ~0.1 dex (~25%). A non-LTE calculation in a supernova ejecta model shows this line becomes significantly stronger, by a factor of around two, alleviating the discrepancy where the 457 nm line in typical models with Mg/O ratios close to solar tended to be too weak compared to observations.

A small group of the newly discovered superluminous supernovae show broad and slowly evolving light curves. Here we present extensive observational data for the slow-evolving superluminous supernova LSQ14an, which brings this group of transients to four in total in the low redshift Universe (z$<$0.2; SN 2007bi, PTF12dam, SN 2015bn). We particularly focus on the optical and near-infrared evolution during the period from 50 days up to 400 days from peak, showing that they are all fairly similar in their light curve and spectral evolution. LSQ14an shows broad, blue-shifted [O III] $\lambda\lambda$4959, 5007 lines, as well as a blue-shifted [O II] $\lambda\lambda$7320, 7330 and [Ca II] $\lambda\lambda$7291, 7323. Furthermore, the sample of these four objects shows common features. Semi-forbidden and forbidden emission lines appear surprisingly early at 50-70 days and remain visible with almost no variation up to 400 days. The spectra remain blue out to 400 days. There are small, but discernible light curve fluctuations in all of them. The light curve of each shows a faster decline than $^{56}$Co after 150 days and it further steepens after 300 days. We also expand our analysis presenting X-ray limits for LSQ14an and SN2015bn and discuss their diagnostic power in interpreting their common features. These features are quite distinct from the faster evolving superluminous supernovae and are not easily explained in terms of only a variation in ejecta mass. While a central engine is still the most likely luminosity source, it appears that the ejecta structure is complex, with multiple emitting zones and at least some interaction between the expanding ejecta and surrounding material.

We present and analyse an extensive dataset of the superluminous supernova (SLSN) LSQ14mo (z = 0.256), consisting of a multi-colour lightcurve from -30 d to +70 d in the rest-frame and a series of 6 spectra from PESSTO covering -7 d to +50 d. This is among the densest spectroscopic coverage, and best-constrained rising lightcurve, for a fast-declining hydrogen-poor SLSN. The bolometric lightcurve can be reproduced with a millisecond magnetar model with ~ 4 M_sol ejecta mass, and the temperature and velocity evolution is also suggestive of a magnetar as the power source. Spectral modelling indicates that the SN ejected ~ 6 M_sol of CO-rich material with a kinetic energy of ~ 7 x 10^51 erg, and suggests a partially thermalised additional source of luminosity between -2 d and +22 d. This may be due to interaction with a shell of material originating from pre-explosion mass loss. We further present a detailed analysis of the host galaxy system of LSQ14mo. PESSTO and GROND imaging show three spatially resolved bright regions, and we used the VLT and FORS2 to obtain a deep (five-hour exposure) spectra of the SN position and the three star-forming regions, which are at a similar redshift. The FORS spectrum at +300 days shows no trace of SN emission lines and we place limits on the strength of [O I] from comparisons with other Ic SNe. The deep spectra provides a unique chance to investigate spatial variations in the host star-formation activity and metallicity. The specific star-formation rate is similar in all three components, as is the presence of a young stellar population. However, the position of LSQ14mo exhibits a lower metallicity, with 12 + log(O/H) = 8.2 in both the R23 and N2 scales (corresponding to ~ 0.3 Z_sol). We propose that the three bright regions in the host system are interacting, which thus triggers star-formation and forms young stellar populations.

When supernovae enter the nebular phase after a few months, they reveal spectral fingerprints of their deep interiors, glowing by radioactivity produced in the explosion. We are given a unique opportunity to see what an exploded star looks like inside. The line profiles and luminosities encode information about physical conditions, explosive and hydrostatic nucleosynthesis, and ejecta morphology, which link to the progenitor properties and the explosion mechanism. Here, the fundamental properties of spectral formation of supernovae in the nebular phase are reviewed. The formalism between ejecta morphology and line profile shapes is derived, including effects of scattering and absorption. Line luminosity expressions are derived in various physical limits, with examples of applications from the literature. The physical processes at work in the supernova ejecta, including gamma-ray deposition, non-thermal electron degradation, ionization and excitation, and radiative transfer are described and linked to the computation and application of advanced spectral models. Some of the results derived so far from nebular-phase supernova analysis are discussed.

SN 2014J in M82 is the closest Type Ia supernova (SN Ia) in decades. The proximity allows for detailed studies of supernova physics and provides insights into the circumstellar and interstellar environment. In this work we analyze Spitzer mid-IR data of SN 2014J in the 3.6 and 4.5 {\mu}m wavelength range, together with several other nearby and well-studied SNe Ia. We compile the first composite mid-IR light-curve templates from our sample of SNe~Ia, spanning the range from before peak brightness well into the nebular phase. Our observations indicate that SNe Ia form a very homogeneous class of objects at these wavelengths. Using the low-reddening supernovae for comparison, we constrain possible thermal emission from circumstellar dust around the highly reddened SN 2014J. We also study SNe 2006X and 2007le, where the presence of matter in the circumstellar environment has been suggested. No significant mid-IR excess is detected, allowing us to place upper limits on the amount of pre-existing dust in the circumstellar environment. For SN 2014J, $M_{dust} < 10^{-5}$ M$_{\odot}$ within $r_{dust} \sim 10^{17}$ cm, which is insufficient to account for the observed extinction. Similar limits are obtained for SNe 2006X and 2007le.

Nebular-phase observations and spectral models of Type Ic superluminous supernovae are presented. LSQ14an and SN 2015bn both display late-time spectra similar to galaxy-subtracted spectra of SN 2007bi, and the class shows strong similarity with broad-lined Type Ic SNe such as SN 1998bw. Near-infrared observations of SN 2015bn show a strong Ca II triplet, O I 9263, O I 1.13 um and Mg I 1.50 um, but no distinct He, Si, or S emission. The high Ca II NIR/[Ca II] 7291,7323 ratio of ~2 indicates a high electron density of n_e >~ 10^8 cm^{-3}. Spectral models of oxygen-zone emission are investigated to put constraints on the emitting region. Models require M(O-zone) >~ 10 Msun to produce enough [O I] 6300,6364 luminosity, irrespective of the powering situation and the density. The high oxygen-zone mass, supported by high estimated magnesium masses, points to explosions of massive CO cores, requiring M_ZAMS >~ 40 Msun. Collisions of pair-instability pulsations do not provide enough mass to account for the emission. [O II] and [O III] lines emerge naturally in many models, which strengthens the identification of broad [O II] 7320,7330, [O III] 4363, and [O III] 4959,5007 in some spectra. A small filling factor f <~ 0.01 for the O/Mg zone is needed to produce enough luminosity in Mg I] 4571, Mg I 1.504 um, and O I recombination lines, which shows that the ejecta is clumped. We review the constraints from the nebular spectral modelling in the context of the various scenarios proposed for superluminous supernovae.

Super-luminous supernovae that radiate more than 10^44 ergs per second at their peak luminosity have recently been discovered in faint galaxies at redshifts of 0.1-4. Some evolve slowly, resembling models of 'pair-instability' supernovae. Such models involve stars with original masses 140-260 times that of the Sun that now have carbon-oxygen cores of 65-30 solar masses. In these stars, the photons that prevent gravitational collapse are converted to electron-positron pairs, causing rapid contraction and thermonuclear explosions. Many solar masses of 56Ni are synthesized; this isotope decays to 56Fe via 56Co, powering bright light curves. Such massive progenitors are expected to have formed from metal-poor gas in the early Universe. Recently, supernova 2007bi in a galaxy at redshift 0.127 (about 12 billion years after the Big Bang) with a metallicity one-third that of the Sun was observed to look like a fading pair-instability supernova. Here we report observations of two slow-to-fade super-luminous supernovae that show relatively fast rise times and blue colours, which are incompatible with pair-instability models. Their late-time light-curve and spectral similarities to supernova 2007bi call the nature of that event into question. Our early spectra closely resemble typical fast-declining super-luminous supernovae, which are not powered by radioactivity. Modelling our observations with 10-16 solar masses of magnetar-energized ejecta demonstrates the possibility of a common explosion mechanism. The lack of unambiguous nearby pair-instability events suggests that their local rate of occurrence is less than 6x10^-6 times that of the core-collapse rate.

Due to its proximity, SN 1987A offers a unique opportunity to directly observe the geometry of a stellar explosion as it unfolds. Here we present spectral and imaging observations of SN 1987A obtained ~10,000 days after the explosion with HST/STIS and VLT/SINFONI at optical and near-infrared wavelengths. These observations allow us to produce the most detailed 3D map of H-alpha to date, the first 3D maps for [Ca II] \lambda \lambda 7292, 7324, [O I] \lambda \lambda 6300, 6364 and Mg II \lambda \lambda 9218, 9244, as well as new maps for [Si I]+[Fe II] 1.644 \mu m and He I 2.058 \mu m. A comparison with previous observations shows that the [Si I]+[Fe II] flux and morphology have not changed significantly during the past ten years, providing evidence that it is powered by 44Ti. The time-evolution of H-alpha shows that it is predominantly powered by X-rays from the ring, in agreement with previous findings. All lines that have sufficient signal show a similar large-scale 3D structure, with a north-south asymmetry that resembles a broken dipole. This structure correlates with early observations of asymmetries, showing that there is a global asymmetry that extends from the inner core to the outer envelope. On smaller scales, the two brightest lines, H-alpha and [Si I]+[Fe II] 1.644 \mu m, show substructures at the level of ~ 200 - 1000 km/s and clear differences in their 3D geometries. We discuss these results in the context of explosion models and the properties of dust in the ejecta.

We present nebular-phase imaging and spectroscopy for the hydrogen-poor superluminous supernova SN 2015bn, at redshift z=0.1136, spanning +250-400 d after maximum light. The light curve exhibits a steepening in the decline rate from 1.4 mag/(100 d) to 1.7 mag/(100 d), suggestive of a significant decrease in the opacity. This change is accompanied by a transition from a blue continuum superposed with photospheric absorption lines to a nebular spectrum dominated by emission lines of oxygen, calcium and magnesium. There are no obvious signatures of circumstellar interaction or large nickel mass. We show that the spectrum at +400 d is virtually identical to a number of energetic Type Ic supernovae such as SN 1997dq, SN 2012au, and SN 1998bw, indicating similar core conditions and strengthening the link between `hypernovae'/long gamma-ray bursts and superluminous supernovae. A single explosion mechanism may unify these events that span absolute magnitudes of -22 < M_B < -17. Both the light curve and spectrum of SN 2015bn are consistent with an engine-driven explosion ejecting 7-30 M$_\odot$ of oxygen-dominated ejecta (for reasonable choices in temperature and opacity). A strong and relatively narrow O II $\lambda$7774 line, seen in a number of these energetic events but not in normal supernovae, may point to an inner shell that is the signature of a central engine.

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 observations of SN 2015bn (= PS15ae = CSS141223-113342+004332 = MLS150211-113342+004333), a Type I superluminous supernova (SLSN) at redshift $z=0.1136$. As well as being one of the closest SLSNe I yet discovered, it is intrinsically brighter ($M_U\approx-23.1$) and in a fainter galaxy ($M_B\approx-16.0$) than other SLSNe at $z\sim0.1$. We used this opportunity to collect the most extensive dataset for any SLSN I to date, including densely-sampled spectroscopy and photometry, from the UV to the NIR, spanning $-$50 to +250 days from optical maximum. SN 2015bn fades slowly, but exhibits surprising undulations in the light curve on a timescale of 30-50 days, especially in the UV. The spectrum shows extraordinarily slow evolution except for a rapid transformation between +7 and +20-30 days. No narrow emission lines from slow-moving material are observed at any phase. We derive physical properties including the bolometric luminosity, and find slow velocity evolution and non-monotonic temperature and radial evolution. A deep radio limit rules out a healthy off-axis gamma-ray burst, and places constraints on the pre-explosion mass loss. The data can be consistently explained by a $\gtrsim10\,{\rm M}_\odot$ stripped progenitor exploding with $\sim 10^{51}\,$erg kinetic energy, forming a magnetar with a spin-down timescale of $\sim20$ days (thus avoiding a gamma-ray burst) that reheats the ejecta and drives ionization fronts. The most likely alternative scenario -- interaction with $\sim20\,{\rm M}_\odot$ of dense, inhomogeneous circumstellar material -- can be tested with continuing radio follow-up.

We investigate two stripped-envelope supernovae (SNe) discovered in the nearby galaxy NGC 5806 by the (i)PTF. These SNe, designated PTF12os/SN 2012P and iPTF13bvn, exploded at a similar distance from the host-galaxy center. We classify PTF12os as a Type IIb SN based on our spectral sequence; iPTF13bvn has previously been classified as Type Ib having a likely progenitor with zero age main sequence (ZAMS) mass below ~17 solar masses. Our main objective is to constrain the explosion parameters of iPTF12os and iPTF13bvn, and to put constraints on the SN progenitors. We present comprehensive datasets on the SNe, and introduce a new reference-subtraction pipeline (FPipe) currently in use by the iPTF. We perform a detailed study of the light curves (LCs) and spectral evolution of the SNe. The bolometric LCs are modeled using the hydrodynamical code HYDE. We use nebular models and late-time spectra to constrain the ZAMS mass of the progenitors. We perform image registration of ground-based images of PTF12os to archival HST images of NGC 5806 to identify a potential progenitor candidate. Our nebular spectra of iPTF13bvn indicate a low ZAMS mass of ~12 solar masses for the progenitor. The late-time spectra of PTF12os are consistent with a ZAMS mass of ~15 solar masses. We successfully identify a progenitor candidate to PTF12os using archival HST images. This source is consistent with being a cluster of massive stars. Our hydrodynamical modeling suggests that the progenitor of PTF12os had a compact He core with a mass of 3.25 solar masses, and that 0.063 solar masses of strongly mixed 56Ni was synthesized. Spectral comparisons to the Type IIb SN 2011dh indicate that the progenitor of PTF12os was surrounded by a hydrogen envelope with a mass lower than 0.02 solar masses. We also find tentative evidence that the progenitor of iPTF13bvn could have been surrounded by a small amount of hydrogen.

We present multi-band ultraviolet, optical, and near-infrared photometry, along with visual-wavelength spectroscopy, of supernova (SN) 2014G in the nearby galaxy NGC 3448 (25 Mpc). The early-phase spectra show strong emission lines of the high ionisation species He II/N IV/C IV during the first 2-3 d after explosion, traces of a metal-rich CSM probably due to pre-explosion mass loss events. These disappear by day 9 and the spectral evolution then continues matching that of normal Type II SNe. The post-maximum light curve declines at a rate typical of Type II-L class. The extensive photometric coverage tracks the drop from the photospheric stage and constrains the radioactive tail, with a steeper decline rate than that expected from the $^{56}$Co decay if $\gamma$-rays are fully trapped by the ejecta. We report the appearance of an unusual feature on the blue-side of H$\alpha$ after 100 d, which evolves to appear as a flat spectral feature linking H$\alpha$ and the O I doublet. This may be due to interaction of the ejecta with a strongly asymmetric, and possibly bipolar CSM. Finally, we report two deep spectra at ~190 and 340 d after explosion, the latter being arguably one of the latest spectra for a Type II-L SN. By modelling the spectral region around the Ca II, we find a supersolar Ni/Fe production. The strength of the O I $\lambda\lambda$6300,6363 doublet, compared with synthetic nebular spectra, suggests a progenitor with a zero-age main-sequence mass between 15 and 19 M$_\odot$.

We present optical and near-infrared photometric and spectroscopic observations of SN 2013ej, in galaxy M74, from 1 to 450 days after the explosion. SN 2013ej is a hydrogen-rich supernova, classified as a Type IIL due to its relatively fast decline following the initial peak. It has a relatively high peak luminosity (absolute magnitude M$_\rm{V}$ = -17.6) but a small $^{56}$Ni production of ~0.023 M$_\odot$. Its photospheric evolution is similar to other Type II SNe, with shallow absorption in the H{\alpha} profile typical for a Type IIL. During transition to the radioactive decay tail at ~100 days, we find the SN to grow bluer in B - V colour, in contrast to some other Type II supernovae. At late times, the bolometric light curve declined faster than expected from $^{56}$Co decay and we observed unusually broad and asymmetric nebular emission lines. Based on comparison of nebular emission lines most sensitive to the progenitor core mass, we find our observations are best matched to synthesized spectral models with a M$_\rm{ZAMS}$ = 12 - 15 M$_\odot$ progenitor. The derived mass range is similar to but not higher than the mass estimated for Type IIP progenitors. This is against the idea that Type IIL are from more massive stars. Observations are consistent with the SN having a progenitor with a relatively low-mass envelope.

High-quality collections of Type II supernova (SN) light curves are scarce because they evolve for hundreds of days, making follow-up observations time consuming and often extending over multiple observing seasons. In light of these difficulties, the diversity of SNe II is not fully understood. Here we present ultraviolet and optical photometry of 12 SNe II monitored by the Las Cumbres Observatory Global Telescope Network (LCOGT) during 2013-2014, and compare them with previously studied SNe having well-sampled light curves. We explore SN II diversity by searching for correlations between the slope of the linear light-curve decay after maximum light (historically used to divide SNe II into IIL and IIP) and other measured physical properties. While SNe IIL are found to be on average more luminous than SNe IIP, SNe IIL do not appear to synthesize more 56Ni than SNe IIP. Finally, optical nebular spectra obtained for several SNe in our sample are found to be consistent with models of red supergiant progenitors in the 12-16 Msun range. Consequently, SNe IIL appear not to account for the deficit of massive red supergiants as SN II progenitors.

Both CO and SiO have been observed at early and late phases in SN 1987A. H_2 was predicted to form at roughly the same time as these molecules, but was not detected at early epochs. Here we report the detection of NIR lines from H_2 at 2.12 mu and 2.40 mu in VLT/SINFONI spectra obtained between days 6489 and 10,120. The emission is concentrated to the core of the supernova in contrast to H-alpha and approximately coincides with the [Si I]/[Fe II] emission detected previously in the ejecta. Different excitation mechanisms and power sources of the emission are discussed. From the nearly constant H_2 luminosities we favour excitation resulting from the 44Ti decay.

We assemble a sample of 24 hydrogen-poor super-luminous supernovae (SLSNe). Parameterizing the light curve shape through rise and decline timescales shows that the two are highly correlated. Magnetar-powered models can reproduce the correlation, with the diversity in rise and decline rates driven by the diffusion timescale. Circumstellar interaction models can exhibit a similar rise-decline relation, but only for a narrow range of densities, which may be problematic for these models. We find that SLSNe are approximately 3.5 magnitudes brighter and have light curves 3 times broader than SNe Ibc, but that the intrinsic shapes are similar. There are a number of SLSNe with particularly broad light curves, possibly indicating two progenitor channels, but statistical tests do not cleanly separate two populations. The general spectral evolution is also presented. Velocities measured from Fe II are similar for SLSNe and SNe Ibc, suggesting that diffusion time differences are dominated by mass or opacity. Flat velocity evolution in most SLSNe suggests a dense shell of ejecta. If opacities in SLSNe are similar to other SNe Ibc, the average ejected mass is higher by a factor 2-3. Assuming $\kappa=0.1\,$cm$^2\,$g$^{-1}$, we estimate a mean (median) SLSN ejecta mass of 10$\,$M$_\odot$ (6$\,$M$_\odot$), with a range of 3-30$\,$M$_\odot$. Doubling the assumed opacity brings the masses closer to normal SNe Ibc, but with a high-mass tail. The most probable mechanism for generating SLSNe seems to be the core-collapse of a very massive hydrogen-poor star, forming a millisecond magnetar.

The observational effects of the 'Infrared Catastrophe' are discussed in view of the very late observations of the Type Ia SN 2011fe. Our model spectra at 1000d take non-local radiative transfer into account, and find that this has a crucial impact on the spectral formation. Although rapid cooling of the ejecta to a few 100 K occurs also in these models, the late-time optical/NIR flux is brighter by 1-2 magnitudes due to redistribution of UV emissivity, resulting from non-thermal excitation and ionization. This effect brings models into better agreement with late-time observations of SN 2011fe and other Type Ia supernovae, and offers a solution to the long standing discrepancy between models and observations. The models show that spectral formation shifts from Fe II and Fe III at 300d to Fe I at 1000d, which explains the apparent wavelength shifts seen in SN2011fe. We discuss effects of time dependence and energy input from 57Co, finding both to be important at 1000d.

If very massive stars (M >~ 100 Msun) can form and avoid too strong mass loss during their evolution, they are predicted to explode as pair-instability supernovae (PISNe). One critical test for candidate events is whether their nucleosynthesis yields and internal ejecta structure, being revealed through nebular-phase spectra at t >~ 1 yr, match those of model predictions. Here we compute theoretical spectra based on model PISN ejecta at 1-3 years post-explosion to allow quantitative comparison with observations. The high column densities of PISNe lead to complete gamma-ray trapping for t >~ 2 years which, combined with fulfilled conditions of steady state, leads to bolometric supernova luminosities matching the 56Co decay. Most of the gamma-rays are absorbed by the deep-lying iron and silicon/sulphur layers. The ionization balance shows a predominantly neutral gas state, which leads to emission lines of Fe I, Si I, and S I. For low-mass PISNe the metal core expands slowly enough to produce a forest of distinct lines, whereas high-mass PISNe expand faster and produce more featureless spectra. Line blocking is complete below ~5000 A for several years, and the model spectra are red. The strongest line is typically [Ca II] 7291,7323, one of few lines from ionized species. We compare our models with proposed PISN candidates SN 2007bi and PTF12dam, finding discrepancies for several key observables and thus no support for a PISN interpretation. We discuss distinct spectral features predicted by the models, and the possibility of detecting pair-instability explosions among non-superluminous supernovae.

We present optical and near infrared observations of the type IIb supernova (SN) 2011fu from a few days to $\sim300$ d after explosion. The SN presents a double-peaked light curve (LC) similar to that of SN 1993J, although more luminous and with a longer cooling phase after the primary peak. The spectral evolution is also similar to SN 1993J's, with hydrogen dominating the spectra to $\sim40$ d, then helium gaining strength, and nebular emission lines appearing from $\sim60$ d post-explosion. The velocities derived from the P-Cygni absorptions are overall similar to those of other type IIb SNe. We have found a strong similarity between the oxygen and magnesium line profiles at late times, which suggests that these lines are forming at the same location within the ejecta. The hydrodynamical modelling of the pseudo-bolometric LC and the observed photospheric velocities suggest that SN 2011fu was the explosion of an extended star ($\rm R\sim450$ R$_\odot$), in which 1.3 $\times 10^{51}$ erg of kinetic energy were released and 0.15 M$_{\rm \odot}$ of $^{56}$Ni were synthesised. In addition, a better reproduction of the observed early pseudo-bolometric LC is achieved if a more massive H-rich-envelope than for other type IIb SNe is considered (0.3 M$_{\rm \odot}$). The hydrodynamical modelling of the LC and the comparison of our late-time spectra with nebular spectral models for type IIb SNe, point to a progenitor for SN 2011fu with a ZAMS mass of 13-18 M$_{\rm \odot}$.

We present our findings on a supernova (SN) impostor, SNHunt248, based on optical and near-IR data spanning $\sim$15 yrs before discovery, to $\sim$1 yr post-discovery. The light curve displays three distinct peaks, the brightest of which is at $M_{R} \sim -15.0$ mag. The post-discovery evolution is consistent with the ejecta from the outburst interacting with two distinct regions of circumstellar material. The 0.5 - 2.2 $\mu$m spectral energy distribution at -740 d is well-matched by a single 6700 K blackbody with $\log(L/L_\odot) \sim 6.1$. This temperature and luminosity support previous suggestions of a yellow hypergiant progenitor; however, we find it to be brighter than the brightest and most massive Galactic late-F to early-G spectral type hypergiants. Overall the historical light curve displays variability of up to $\sim \pm1$ mag. At current epochs ($\sim$1 yr post-outburst), the absolute magnitude ($M_{R} \sim -9$ mag) is just below the faintest observed historical absolute magnitude $\sim$10 yrs before discovery.

We report on our findings based on the analysis of observations of the Type II-L supernova LSQ13cuw within the framework of currently accepted physical predictions of core-collapse supernova explosions. LSQ13cuw was discovered within a day of explosion, hitherto unprecedented for Type II-L supernovae. This motivated a comparative study of Type II-P and II-L supernovae with relatively well-constrained explosion epochs and rise times to maximum (optical) light. From our sample of twenty such events, we find evidence of a positive correlation between the duration of the rise and the peak brightness. On average, SNe II-L tend to have brighter peak magnitudes and longer rise times than SNe II-P. However, this difference is clearest only at the extreme ends of the rise time versus peak brightness relation. Using two different analytical models, we performed a parameter study to investigate the physical parameters that control the rise time behaviour. In general, the models qualitatively reproduce aspects of the observed trends. We find that the brightness of the optical peak increases for larger progenitor radii and explosion energies, and decreases for larger masses. The dependence of the rise time on mass and explosion energy is smaller than the dependence on the progenitor radius. We find no evidence that the progenitors of SNe II-L have significantly smaller radii than those of SNe II-P.

Super-luminous supernovae of type Ic have a tendency to occur in faint host galaxies which are likely to have low mass and low metallicity. PTF12dam is one of the closest and best studied super-luminous explosions that has a broad and slowly fading lightcurve similar to SN 2007bi. Here we present new photometry and spectroscopy for PTF12dam from 200-500 days (rest-frame) after peak and a detailed analysis of the host galaxy (SDSS J142446.21+461348.6 at z = 0.107). Using deep templates and image subtraction we show that the full lightcurve can be fit with a magnetar model if escape of high-energy gamma rays is taken into account. The full bolometric lightcurve from -53 to +399 days (with respect to peak) cannot be fit satisfactorily with the pair-instability models. An alternative model of interaction with a dense CSM produces a good fit to the data although this requires a very large mass (~ 13 M_sun) of hydrogen free CSM. The host galaxy is a compact dwarf (physical size ~ 1.9 kpc) and with M_g = -19.33 +/- 0.10, it is the brightest nearby SLSN Ic host discovered so far. The host is a low mass system (2.8 x 10^8 M_sun) with a star-formation rate (5.0 M_sun/year), which implies a very high specific star-formation rate (17.9 Gyr^-1). The remarkably strong nebular lines provide detections of the [O III] \lambda 4363 and [O II] \lambda\lambda 7320,7330 auroral lines and an accurate oxygen abundance of 12 + log(O/H) = 8.05 +/- 0.09. We show here that they are at the extreme end of the metallicity distribution of dwarf galaxies and propose that low metallicity is a requirement to produce these rare and peculiar supernovae.

We present data for LSQ14bdq, a hydrogen-poor super-luminous supernova (SLSN) discovered by the La Silla QUEST survey and classified by the Public ESO Spectroscopic Survey of Transient Objects. The spectrum and light curve are very similar to slow-declining SLSNe such as PTF12dam. However, detections within ~1 day after explosion show a bright and relatively fast initial peak, lasting for ~15 days, prior to the usual slow rise to maximum light. The broader, main peak can be fit with either central engine or circumstellar interaction models. We discuss the implications of the precursor peak in the context of these models. It is too bright and narrow to be explained as a normal 56Ni-powered SN, and we suggest that interaction models may struggle to fit the two peaks simultaneously. We propose that the initial peak may arise from the post-shock cooling of extended stellar material, and reheating by a central engine drives the second peak. In this picture, we show that an explosion energy of ~2x10^{52} erg and a progenitor radius of a few hundred solar radii would be required to power the early emission. The competing engine models involve rapidly spinning magnetars (neutron stars) or fall-back accretion onto a central black hole. The prompt energy required may favour the black hole scenario. The bright initial peak effectively rules out a compact Wolf-Rayet star as a progenitor, since the inferred energies and ejected masses become unphysical.

We present imaging and spectroscopic observations with HST and VLT of the ring of SN 1987A from 1994 to 2014. After an almost exponential increase of the shocked emission from the hotspots up to day ~8,000 (~2009), both this and the unshocked emission are now fading. From the radial positions of the hotspots we see an acceleration of these up to 500-1000 km/s, consistent with the highest spectroscopic shock velocities from the radiative shocks. In the most recent observations (2013 and 2014), we find several new hotspots outside the inner ring, excited by either X-rays from the shocks or by direct shock interaction. All of these observations indicate that the interaction with the supernova ejecta is now gradually dissolving the hotspots. We predict, based on the observed decay, that the inner ring will be destroyed by ~2025.

We present optical and near-infrared (NIR) photometry and spectroscopy as well as modelling of the lightcurves of the Type IIb supernova (SN) 2011dh. Our extensive dataset spans two years, and complemented with Spitzer mid-infrared (MIR) data, we use it to build an optical-to-MIR bolometric lightcurve between days 3 and 732. To model the bolometric lightcurve before day 400 we use a grid of hydrodynamical SN models and a bolometric correction determined with steady-state NLTE modelling. Using this method we find a helium core mass of 3.1 (+0.7-0.4) solar masses for SN 2011dh, consistent within error bars with previous results. We present bolometric and broad-band lightcurves between days 100 and 500 for the Jerkstrand et al. (2014) steady-state NLTE models, and the preferred 12 solar masses (initial mass) model shows a good overall agreement with the observed lightcurves. We find an excess in the K and the MIR bands developing between days 100 and 250, during which an increase in the optical decline rate is also observed. A local origin of the excess is suggested by the depth of the HeI 2.058 micron absorption. Steady-state NLTE models with a modest dust opacity in the core, turned on during this period, reproduce the observed behaviour, but an additional excess in the Spitzer 4.5 micron band remains. CO first-overtone band emission is detected at day 206, and possibly at day 89, and assuming the additional excess to be dominated by CO fundamental band emission, we find fundamental to first-overtone band ratios considerably higher than observed in SN 1987A. This paper concludes our extensive observational and modelling work on SN 2011dh. The results from hydrodynamical modelling, steady-state NLTE modelling, and stellar evolutionary progenitor analysis are all consistent, and suggest an initial mass of 12 solar masses for the progenitor.

Measurements of explosive nucleosynthesis yields in core-collapse supernovae provide tests for explosion models. We investigate constraints on explosive conditions derivable from measured amounts of nickel and iron after radioactive decays using nucleosynthesis networks with parameterized thermodynamic trajectories. The Ni/Fe ratio is for most regimes dominated by the production ratio of 58Ni/(54Fe + 56Ni), which tends to grow with higher neutron excess and with higher entropy. For SN 2012ec, a supernova that produced a Ni/Fe ratio of $3.4\pm1.2$ times solar, we find that burning of a fuel with neutron excess $\eta \approx 6\times 10^{-3}$ is required. Unless the progenitor metallicity is over 5 times solar, the only layer in the progenitor with such a neutron excess is the silicon shell. Supernovae producing large amounts of stable nickel thus suggest that this deep-lying layer can be, at least partially, ejected in the explosion. We find that common spherically symmetric models of $M_{\rm ZAMS} \lesssim 13$ Msun stars exploding with a delay time of less than one second ($M_{\rm cut} < 1.5$ Msun) are able to achieve such silicon-shell ejection. Supernovae that produce solar or sub-solar Ni/Fe ratios, such as SN 1987A, must instead have burnt and ejected only oxygen-shell material, which allows a lower limit to the mass cut to be set. Finally, we find that the extreme Ni/Fe value of 60-75 times solar derived for the Crab cannot be reproduced by any realistic-entropy burning outside the iron core, and neutrino-neutronization obtained in electron-capture models remains the only viable explanation.

The Public European Southern Observatory Spectroscopic Survey of Transient Objects (PESSTO) began as a public spectroscopic survey in April 2012. We describe the data reduction strategy and data products which are publicly available through the ESO archive as the Spectroscopic Survey Data Release 1 (SSDR1). PESSTO uses the New Technology Telescope with EFOSC2 and SOFI to provide optical and NIR spectroscopy and imaging. We target supernovae and optical transients brighter than 20.5mag for classification. Science targets are then selected for follow-up based on the PESSTO science goal of extending knowledge of the extremes of the supernova population. The EFOSC2 spectra cover 3345-9995A (at resolutions of 13-18 Angs) and SOFI spectra cover 0.935-2.53 micron (resolutions 23-33 Angs) along with JHK imaging. This data release contains spectra from the first year (April 2012 - 2013), consisting of all 814 EFOSC2 spectra and 95 SOFI spectra (covering 298 distinct objects), in standard ESO Phase 3 format. We estimate the accuracy of the absolute flux calibrations for EFOSC2 to be typically 15%, and the relative flux calibration accuracy to be about 5%. The PESSTO standard NIR reduction process does not yet produce high accuracy absolute spectrophotometry but the SOFI JHK imaging will improve this. Future data releases will focus on improving the automated flux calibration of the data products.

SN 2012ec is a Type IIP supernova (SN) with a progenitor detection and comprehensive photospheric-phase observational coverage. Here, we present Very Large Telescope and PESSTO observations of this SN in the nebular phase. We model the nebular [O I] 6300, 6364 lines and find their strength to suggest a progenitor main-sequence mass of 13-15 Msun. SN 2012ec is unique among hydrogen-rich SNe in showing a distinct and unblended line of stable nickel [Ni II] 7378. This line is produced by 58Ni, a nuclear burning ash whose abundance is a sensitive tracer of explosive burning conditions. Using spectral synthesis modelling, we use the relative strengths of [Ni II] 7378 and [Fe II] 7155 (the progenitor of which is 56Ni) to derive a Ni/Fe production ratio of 0.20pm0.07 (by mass), which is a factor 3.4pm1.2 times the solar value. High production of stable nickel is confirmed by a strong [Ni II] 1.939 micron line. This is the third reported case of a core-collapse supernova producing a Ni/Fe ratio far above the solar value, which has implications for core-collapse explosion theory and galactic chemical evolution models.

We present photometric and spectroscopic observations of the interacting transient SN 2009ip taken during the 2013 and 2014 observing seasons. We characterise the photometric evolution as a steady and smooth decline in all bands, with a decline rate that is slower than expected for a solely $^{56}$Co-powered supernova at late phases. No further outbursts or eruptions were seen over a two year period from 2012 December until 2014 December. SN 2009ip remains brighter than its historic minimum from pre-discovery images. Spectroscopically, SN 2009ip continues to be dominated by strong, narrow ($\lesssim$2000 km~s$^{-1}$) emission lines of H, He, Ca, and Fe. While we make tenuous detections of [Fe~{\sc ii}] $\lambda$7155 and [O~{\sc i}] $\lambda\lambda$6300,6364 lines at the end of 2013 June and the start of 2013 October respectively, we see no strong broad nebular emission lines that could point to a core-collapse origin. In general, the lines appear relatively symmetric, with the exception of our final spectrum in 2014 May, when we observe the appearance of a redshifted shoulder of emission at +550 km~s$^{-1}$. The lines are not blue-shifted, and we see no significant near- or mid-infrared excess. From the spectroscopic and photometric evolution of SN 2009ip until 820 days after the start of the 2012a event, we still see no conclusive evidence for core-collapse, although whether any such signs could be masked by ongoing interaction is unclear.

We present the results of a photometric and spectroscopic monitoring campaign of SN 2012ec, which exploded in the spiral galaxy NGC 1084, during the photospheric phase. The photometric light curve exhibits a plateau with luminosity L= 0.9 x 10^{42} erg/s and duration ~90 days, which is somewhat shorter than standard Type II-P supernovae. We estimate the nickel mass as 0.040 +/- 0.015 Msun from the luminosity at the beginning of the radioactive tail of the light curve. The explosion parameters of SN 2012ec were estimated from the comparison of the bolometric light curve and the observed temperature and velocity evolution of the ejecta with predictions from hydrodynamical models. We derived an envelope mass of 12.6 Msun, an initial progenitor radius of 1.6 x 10^{13} cm and an explosion energy of 1.2 foe. These estimates agree with an independent study of the progenitor star identified in pre-explosion images, for which an initial mass of M=14-22 Msun was determined. We have applied the same analysis to two other type II-P supernovae (SNe 2012aw and 2012A), and carried out a comparison with the properties of SN 2012ec derived in this paper. We find a reasonable agreement between the masses of the progenitors obtained from pre-explosion images and masses derived from hydrodynamical models. We estimate the distance to SN 2012ec with the Standardized Candle Method (SCM) and compare it with other estimates based on other primary and secondary indicators. SNe 2012A, 2012aw and 2012ec all follow the standard relations for the SCM for the use of Type II-P SNe as distance indicators.

We present observational data for a peculiar supernova discovered by the OGLE-IV survey and followed by the Public ESO Spectroscopic Survey for Transient Objects. The inferred redshift of $z=0.07$ implies an absolute magnitude in the rest-frame $I$-band of M$_{I}\sim-17.6$ mag. This places it in the luminosity range between normal Type Ia SNe and novae. Optical and near infrared spectroscopy reveal mostly Ti and Ca lines, and an unusually red color arising from strong depression of flux at rest wavelengths $<5000$ \AA. To date, this is the only reported SN showing Ti-dominated spectra. The data are broadly consistent with existing models for the pure detonation of a helium shell around a low-mass CO white dwarf and "double-detonation" models that include a secondary detonation of a CO core following a primary detonation in an overlying helium shell.

We investigate line formation processes in Type IIb supernovae (SNe) from 100 to 500 days post-explosion using spectral synthesis calculations. The modeling identifies the nuclear burning layers and physical mechanisms that produce the major emission lines, and the diagnostic potential of these. We compare the model calculations with data on the three best observed Type IIb SNe to-date - SN 1993J, SN 2008ax, and SN 2011dh. Oxygen nucleosynthesis depends sensitively on the main-sequence mass of the star and modeling of the [O I] 6300, 6364 lines constrains the progenitors of these three SNe to the M_ZAMS=12-16 M_sun range (ejected oxygen masses 0.3-0.9 M_sun), with SN 2011dh towards the lower end and SN 1993J towards the upper end of the range. The high ejecta masses from M_ZAMS >= 17 M_sun progenitors give rise to brighter nebular phase emission lines than observed. Nucleosynthesis analysis thus supports a scenario of low/moderate mass progenitors for Type IIb SNe, and by implication an origin in binary systems. We demonstrate how oxygen and magnesium recombination lines may be combined to diagnose the magnesium mass in the SN ejecta. For SN 2011dh, a magnesium mass of of 0.02-0.14 M_sun is derived, which gives a Mg/O production ratio consistent with the solar value. Nitrogen left in the He envelope from CNO-burning gives strong [N II] 6548, 6583 emission lines that dominate over H-alpha emission in our models. The hydrogen envelopes of Type IIb SNe are too small and dilute to produce any noticeable H-alpha emission or absorption after ~150 days, and nebular phase emission seen around 6550 A is in many cases likely caused by [N II] 6548, 6583. Finally, the influence of radiative transport on the emergent line profiles is investigated...(abridged)

We present new Herschel photometric and spectroscopic observations of Supernova 1987A, carried out in 2012. Our dedicated photometric measurements provide new 70 micron data and improved imaging quality at 100 and 160 micron compared to previous observations in 2010. Our Herschel spectra show only weak CO line emission, and provide an upper limit for the 63 micron [O I] line flux, eliminating the possibility that line contaminations distort the previously estimated dust mass. The far-infrared spectral energy distribution (SED) is well fitted by thermal emission from cold dust. The newly measured 70 micron flux constrains the dust temperature, limiting it to nearly a single temperature. The far-infrared emission can be fitted by 0.5+-0.1 Msun of amorphous carbon, about a factor of two larger than the current nucleosynthetic mass prediction for carbon. The observation of SiO molecules at early and late phases suggests that silicates may also have formed and we could fit the SED with a combination of 0.3 Msun of amorphous carbon and 0.5 Msun of silicates, totalling 0.8 Msun of dust. Our analysis thus supports the presence of a large dust reservoir in the ejecta of SN 1987A. The inferred dust mass suggests that supernovae can be an important source of dust in the interstellar medium, from local to high-redshift galaxies.

We present optical spectra and light curves for three hydrogen-poor super-luminous supernovae followed by the Public ESO Spectroscopic Survey of Transient Objects (PESSTO). Time series spectroscopy from a few days after maximum light to 100 days later shows them to be fairly typical of this class, with spectra dominated by Ca II, Mg II, Fe II and Si II, which evolve slowly over most of the post-peak photospheric phase. We determine bolometric light curves and apply simple fitting tools, based on the diffusion of energy input by magnetar spin-down, 56Ni decay, and collision of the ejecta with an opaque circumstellar shell. We investigate how the heterogeneous light curves of our sample (combined with others from the literature) can help to constrain the possible mechanisms behind these events. We have followed these events to beyond 100-200 days after peak, to disentangle host galaxy light from fading supernova flux and to differentiate between the models, which predict diverse behaviour at this phase. Models powered by radioactivity require unrealistic parameters to reproduce the observed light curves, as found by previous studies. Both magnetar heating and circumstellar interaction still appear to be viable candidates. A large diversity is emerging in observed tail-phase luminosities, with magnetar models failing in some cases to predict the rapid drop in flux. This would suggest either that magnetars are not responsible, or that the X-ray flux from the magnetar wind is not fully trapped. The light curve of one object shows a distinct re-brightening at around 100d after maximum light. We argue that this could result either from multiple shells of circumstellar material, or from a magnetar ionisation front breaking out of the ejecta.

We present nebular phase optical and near-infrared spectroscopy of the Type IIP supernova SN 2012aw combined with NLTE radiative transfer calculations applied to ejecta from stellar evolution/explosion models. Our spectral synthesis models generally show good agreement with the ejecta from a MZAMS = 15 Msun progenitor star. The emission lines of oxygen, sodium, and magnesium are all consistent with the nucleosynthesis in a progenitor in the 14 - 18 Msun range. We also demonstrate how the evolution of the oxygen cooling lines of [O I] 5577 A, [O I] 6300 A, and [O I] 6364 A can be used to constrain the mass of oxygen in the non-molecularly cooled ashes to < 1 Msun, independent of the mixing in the ejecta. This constraint implies that any progenitor model of initial mass greater than 20 Msun would be difficult to reconcile with the observed line strengths. A stellar progenitor of around MZAMS = 15 Msun can consistently explain the directly measured luminosity of the progenitor star, the observed nebular spectra, and the inferred pre-supernova mass-loss rate. We conclude that there is still no convincing example of a Type IIP explosion showing the nucleosynthesis expected from a MZAMS > 20 Msun progenitor.

We present optical and near-infrared (NIR) photometry and spectroscopy of the Type IIb supernova (SN) 2011dh for the first 100 days. We complement our extensive dataset with SWIFT ultra-violet (UV) and Spitzer mid-infrared (MIR) data to build a UV to MIR bolometric lightcurve using both photometric and spectroscopic data. Hydrodynamical modelling of the SN based on this bolometric lightcurve have been presented in Bersten (2012). We find that the absorption minimum for the hydrogen lines is never seen below ~11000 km/s but approaches this value as the lines get weaker. This suggests that the interface between the helium core and hydrogen rich envelope is located near this velocity in agreement with the Bersten et al. (2012) He4R270 ejecta model. Spectral modelling of the hydrogen lines using this ejecta model supports the conclusion and we find a hydrogen mass of 0.01-0.04 solar masses to be consistent with the observed spectral evolution. We estimate that the photosphere reaches the helium core at 5-7 days whereas the helium lines appear between ~10 and ~15 days, close to the photosphere and then move outward in velocity until ~40 days. This suggests that increasing non-thermal excitation due to decreasing optical depth for the gamma-rays is driving the early evolution of these lines. We also provide and discuss pre- and post-explosion observations of the SN site which shows a reduction by 75 percent in flux at the position of the yellow supergiant coincident with SN 2011dh. The B, V and r band decline rates of 0.0073, 0.0090 and 0.0053 mag/day respectively are consistent with the remaining flux being emitted by the SN. Hence we find that the star was indeed the progenitor of SN 2011dh as previously suggested by Maund et al. (2011) and which is also consistent with the results from the hydrodynamical modelling.

We present optical photometric and spectroscopic coverage of the superluminous supernova (SLSN) PS1-11ap, discovered with the Pan-STARRS1 Medium Deep Survey at z = 0.524. This intrinsically blue transient rose slowly to reach a peak magnitude of M_u = -21.4 mag and bolometric luminosity of 8 x 10^43 ergs^-1 before settling onto a relatively shallow gradient of decline. The observed decline is significantly slower than those of the superluminous type Ic SNe which have been the focus of much recent attention. Spectroscopic similarities with the lower redshift SN2007bi and a decline rate similar to 56Co decay timescale initially indicated that this transient could be a candidate for a pair instability supernova (PISN) explosion. Overall the transient appears quite similar to SN2007bi and the lower redshift object PTF12dam. The extensive data set, from 30 days before peak to 230 days after, allows a detailed and quantitative comparison with published models of PISN explosions. We find that the PS1-11ap data do not match these model explosion parameters well, supporting the recent claim that these SNe are not pair instability explosions. We show that PS1-11ap has many features in common with the faster declining superluminous Ic supernovae and the lightcurve evolution can also be quantitatively explained by the magnetar spin down model. At a redshift of z = 0.524 the observer frame optical coverage provides comprehensive restframe UV data and allows us to compare it with the superluminous SNe recently found at high redshifts between z = 2-4. While these high-z explosions are still plausible PISN candidates, they match the photometric evolution of PS1-11ap and hence could be counterparts to this lower redshift transient.

We report optical and near-infrared observations of SN 2012ca with the Public ESO Spectroscopy Survey of Transient Objects (PESSTO), spread over one year since discovery. The supernova (SN) bears many similarities to SN 1997cy and to other events classified as Type IIn but which have been suggested to have a thermonuclear origin with narrow hydrogen lines produced when the ejecta impact a hydrogen-rich circumstellar medium (CSM). Our analysis, especially in the nebular phase, reveals the presence of oxygen, magnesium and carbon features. This suggests a core collapse explanation for SN2012ca, in contrast to the thermonuclear interpretation proposed for some members of this group. We suggest that the data can be explained with a hydrogen and helium deficient SN ejecta (Type I) interacting with a hydrogen-rich CSM, but that the explosion was more likely a Type Ic core-collapse explosion than a Type Ia thermonuclear one. This suggests two channels (both thermonuclear and stripped envelope core-collapse) may be responsible for these SN 1997cy-like events.

We present the one-year long observing campaign of SN 2012A which exploded in the nearby (9.8 Mpc) irregular galaxy NGC 3239. The photometric evolution is that of a normal type IIP supernova. The absolute maximum magnitude, with MB = -16.23 +- 0.16 mag. SN2012A reached a peak luminosity of about 2X10**42 erg/s, which is brighter than those of other SNe with a similar 56Ni mass. The latter was estimated from the luminosity in the exponential tail of the light curve and found to be M(56Ni) = 0.011 +-0.004 Msun. The spectral evolution of SN 2012A is also typical of SN IIP, from the early spectra dominated by a blue continuum and very broad (~10**4 km/s) Balmer lines, to the late-photospheric spectra characterized by prominent P-Cygni features of metal lines (Fe II, Sc II, Ba II, Ti II, Ca II, Na ID). The photospheric velocity is moderately low, ~3X10**3 km/s at 50 days, for the low optical depth metal lines. The nebular spectrum obtained 394 days after the shock breakout shows the typical features of SNe IIP and the strength of the [O I] doublet suggests a progenitor of intermediate mass, similar to SN 2004et (~15 Msun). A candidate progenitor for SN 2012A has been identified in deep, pre-explosion K'-band Gemini North (NIRI) images, and found to be consistent with a star with a bolometric magnitude -7.08+-0.36 (log L/Lsun = 4.73 +- 0.14$ dex). The magnitude of the recovered progenitor in archival images points toward a moderate-mass 10.5 (-2/+4.5) Msun star as the precursor of SN 2012A. The explosion parameters and progenitor mass were also estimated by means of a hydrodynamical model, fitting the bolometric light curve, the velocity and the temperature evolution. We found a best fit for a kinetic energy of 0.48 foe, an initial radius of 1.8X10**13 cm and ejecta mass of 12.5 Msun.

We report extensive observational data for five of the lowest redshift Super-Luminous Type Ic Supernovae (SL-SNe Ic) discovered to date, namely PTF10hgi, SN2011ke, PTF11rks, SN2011kf and SN2012il. Photometric imaging of the transients at +50 to +230 days after peak combined with host galaxy subtraction reveals a luminous tail phase for four of these SL-SNe. A high resolution, optical and near infrared spectrum from xshooter provides detection of a broad He I $\lambda$10830 emission line in the spectrum (+50d) of SN2012il, revealing that at least some SL-SNe Ic are not completely helium free. At first sight, the tail luminosity decline rates that we measure are consistent with the radioactive decay of \co, and would require 1-4M of \ni to produce the luminosity. These \ni masses cannot be made consistent with the short diffusion times at peak, and indeed are insufficient to power the peak luminosity. We instead favour energy deposition by newborn magnetars as the power source for these objects. A semi-analytical diffusion model with energy input from the spin-down of a magnetar reproduces the extensive lightcurve data well. The model predictions of ejecta velocities and temperatures which are required are in reasonable agreement with those determined from our observations. We derive magnetar energies of $0.4\lesssim E$($10^{51}$erg) $\lesssim6.9$ and ejecta masses of $2.3\lesssim M_{ej}$(\M) $\lesssim 8.6$. The sample of five SL-SNe Ic presented here, combined with SN 2010gx - the best sampled SL-SNe Ic so far - point toward an explosion driven by a magnetar as a viable explanation for all SL-SNe Ic.

We present observations of the interacting transient SN 2009ip, from the start of the outburst in October 2012 until the end of the 2012 observing season. The transient reached a peak of $M_V$=-17.7 mag before fading rapidly, with a total integrated luminosity of 1.9$\times10^{49}$ erg over the period of August-December 2012. The optical and near infrared spectra are dominated by narrow emission lines, signaling a dense circumstellar environment, together with multiple components of broad emission and absorption in H and He at velocities between 0.5-1.2$\times10^4$ km s$^{-1}$\. We see no evidence for nucleosynthesized material in SN 2009ip, even in late-time pseudo-nebular spectra. We set a limit of $<$0.02 M$_{\odot}$\ on the mass of any synthesized $^{56}$Ni from the late time lightcurve. A simple model for the narrow Balmer lines is presented, and used to derive number densities for the circumstellar medium of between $\sim 10^{9}-10^{10}$ cm$^{-3}$. Our near-infrared data does not show any excess at longer wavelengths. Our last data, taken in December 2012, shows that SN 2009ip has spectroscopically evolved to something quite similar to its appearance in late 2009, albeit with higher velocities. It is possible that neither of the eruptive and high luminosity events of SN 2009ip were induced by a core-collapse. We show that the peak and total integrated luminosity can be due to the efficient conversion of kinetic energy from colliding ejecta, and that around 0.05-0.1 M$_{\odot}$\ of material moving at 0.5-1$\times10^4$ km s$^{-1}$\ could comfortably produce the observed luminosity. The ejection of multiple shells, lack of evidence for nucleosynthesied elements and broad nebular lines, are all consistent with the pulsational-pair instability scenario. In this case the progenitor star may still exist, and will be observed after the current outburst fades.

We present a study of the morphology of the ejecta in Supernova 1987A based on images and spectra from the HST as well as integral field spectroscopy from VLT/SINFONI. The HST observations were obtained between 1994 - 2011 and primarily probe the outer hydrogen-rich zones of the ejecta. The SINFONI observations were obtained in 2005 and 2011 and instead probe the [Si I]/[Fe II] emission from the inner regions. We find a strong temporal evolution of the morphology in the HST images, from a roughly elliptical shape before ~5,000 days, to a more irregular, edge-brightened morphology thereafter. This transition is a natural consequence of the change in the dominant energy source powering the ejecta, from radioactive decay before ~5,000 days to X-ray input from the circumstellar interaction thereafter. The [Si I]/[Fe II] images display a more uniform morphology, which may be due to a remaining significant contribution from radioactivity in the inner ejecta and the higher abundance of these elements in the core. Both the H-alpha and the [Si I]/[Fe II] line profiles show that the ejecta are distributed fairly close to the plane of the inner circumstellar ring, which is assumed to define the rotational axis of the progenitor. The H-alpha emission extends to higher velocities than [Si I]/[Fe II] as expected. There is no clear symmetry axis for all the emission and we are unable to model the ejecta distribution with a simple ellipsoid model with a uniform distribution of dust. Instead, we find that the emission is concentrated to clumps and that the emission is distributed somewhat closer to the ring in the north than in the south. This north-south asymmetry may be partially explained by dust absorption. We compare our results with explosion models and find some qualitative agreement, but note that the observations show a higher degree of large-scale asymmetry.

We present observations with VLT and HST of the broad emission lines from the inner ejecta and reverse shock of SN 1987A from 1999 until 2012 (days 4381 -- 9100 after explosion). We detect broad lines from H-alpha, H-beta, Mg I], Na I, [O I], [Ca II] and a feature at 9220 A. We identify the latter line with Mg II 9218, 9244,most likely pumped by Ly-alpha fluorescence. H-alpha, and H-beta both have a centrally peaked component, extending to 4500 km/s and a very broad component extending to 11,000 km/s, while the other lines have only the central component. The low velocity component comes from unshocked ejecta, heated mainly by X-rays from the circumstellar ring collision, whereas the broad component comes from faster ejecta passing through the reverse shock. The reverse shock flux in H-alpha has increased by a factor of 4-6 from 2000 to 2007. After that there is a tendency of flattening of the light curve, similar to what may be seen in soft X-rays and in the optical lines from the shocked ring. The core component seen in H-alpha, [Ca II] and Mg II has experienced a similar increase, consistent with that found from HST photometry. The ring-like morphology of the ejecta is explained as a result of the X-ray illumination, depositing energy outside of the core of the ejecta. The energy deposition in the ejecta of the external X-rays illumination is calculated using explosion models for SN 1987A and we predict that the outer parts of the unshocked ejecta will continue to brighten because of this. We finally discuss evidence for dust in the ejecta from line asymmetries.

SN 2004et is one of the nearest and best-observed Type IIP supernovae, with a progenitor detection as well as good photometric and spectroscopic observational coverage well into the nebular phase. Based on nucleosynthesis from stellar evolution/explosion models we apply spectral modeling to analyze its 140-700 day evolution from ultraviolet to mid-infrared. We find a M_ZAMS= 15 Msun progenitor star (with an oxygen mass of 0.8 Msun) to satisfactorily reproduce [O I] 6300, 6364 {\AA} and other emission lines of carbon, sodium, magnesium, and silicon, while 12 Msun and 19 Msun models under- and overproduce most of these lines, respectively. This result is in fair agreement with the mass derived from the progenitor detection, but in disagreement with hydrodynamical modeling of the early-time light curve. From modeling of the mid-infrared iron-group emission lines, we determine the density of the "Ni-bubble" to rho(t) = 7E-14*(t/100d)^-3 g cm^-3, corresponding to a filling factor of f = 0.15 in the metal core region (V = 1800 km/s). We also confirm that silicate dust, CO, and SiO emission are all present in the spectra.

Massive stars live fast and die young. They shine furiously for a few million years, during which time they synthesize most of the heavy elements in the universe in their cores. They end by blowing themselves up in a powerful explosion known as a supernova. During this process, the core collapses to a neutron star or a black hole, while the outer layers are expelled with velocities of thousands of kilometers per second. The resulting fireworks often outshine the entire host galaxy for many weeks. The explosion energy is eventually radiated away, but powering of the newborn nebula continues by radioactive isotopes synthesized in the explosion. The ejecta are now quite transparent, and we can see the material produced in the deep interiors of the star. To interpret the observations, detailed spectral modeling is needed. This thesis aims to develop and apply state-of-the-art computational tools for interpreting and modeling supernova observations in the nebular phase. This requires calculation of the physical conditions throughout the nebula, including non-thermal processes from the radioactivity, thermal and statistical equilibrium, as well as radiative transport. The inclusion of multi-line radiative transfer, which we compute with a Monte Carlo technique, represents one of the major advancements presented in this thesis.

We present a study of the nebular phase spectra of a sample of Type II-Plateau supernovae with identified progenitors or restrictive limits. The evolution of line fluxes, shapes, and velocities are compared within the sample, and interpreted by the use of a spectral synthesis code. The small diversity within the dataset can be explained by strong mixing occurring during the explosion, and by recognising that most lines have significant contributions from primordial metals in the H envelope, which dominates the total ejecta mass in these type of objects. In particular, when using the [O I] 6300, 6364 Angstrom doublet for estimating the core mass of the star, care has to be taken to account for emission from primordial O in the envelope. Finally, a correlation between the H-alpha line width and the mass of 56Ni is presented, suggesting that higher energy explosions are associated with higher 56Ni production.

We present the late-time optical light curve of the ejecta of SN 1987A measured from HST imaging observations spanning the past 17 years. We find that the flux from the ejecta declined up to around year 2001, powered by the radioactive decay of 44Ti. Then the flux started to increase, more than doubling by the end of 2009. We show that the increase is the result of energy deposited by X-rays produced in the interaction with the circumstellar medium. We suggest that the change of the dominant energy input to the ejecta, from internal to external, marks the transition from supernova to supernova remnant. The details of the observations and the modelling are described in the accompanying supplementary information.

SN 1987A provides a unique opportunity to study the evolution of a supernova from explosion into very late phases. Due to the rich chemical structure, the multitude of physical process involved, and extensive radiative transfer effects, detailed modeling is needed to interpret the emission from this and other supernovae. In this paper, we analyze the late-time (~8 years) HST spectrum of the SN 1987A ejecta, where 44Ti is the dominant power source. Based on an explosion model for a 19 Msun progenitor, we compute a model spectrum by calculating the degradation of positrons and gamma-rays from the radioactive decays, solving the equations governing temperature, ionization balance and NLTE level populations, and treating the radiative transfer with a Monte Carlo technique. We obtain a UV/optical/NIR model spectrum which is found to reproduce most of the lines in the observed spectrum to good accuracy. We find non-local radiative transfer in atomic lines to be an important process also at this late stage of the supernova, with ~30% of the emergent flux in the optical and NIR coming from scattering/fluorescence. We investigate the question of where the positrons deposit their energy, and favor the scenario where they are locally trapped in the Fe/He clumps by a magnetic field. Energy deposition into these largely neutral Fe/He clumps makes Fe I lines prominent in the emergent spectrum. Using the best available estimates for the dust extinction, we determine the amount of 44Ti produced in the explosion to 1.5\pm0.5 * 10^-4 Msun.

The interpretation of supernova (SN) spectra is essential for deriving SN ejecta properties such as density and composition, which in turn can tell us about their progenitors and the explosion mechanism. A very large number of atomic processes are important for spectrum formation. Several tools for calculating SN spectra exist, but they mainly focus on the very early or late epochs. The intermediate phase, which requires a NLTE treatment of radiation transport has rarely been studied. In this paper we present a new SN radiation transport code, NERO, which can look at those epochs. All the atomic processes are treated in full NLTE, under a steady-state assumption. This is a valid approach between roughly 50 and 500 days after the explosion depending on SN type. This covers the post-maximum photospheric and the early and the intermediate nebular phase. As a test, we compare NERO to the radiation transport code of Jerkstrand et al. (2011) and to the nebular code of Mazzali et al. (2001). All three codes have been developed independently and a comparison provides a valuable opportunity to investigate their reliability. Currently, NERO is one-dimensional and can be used for predicting spectra of synthetic explosion models or for deriving SN properties by spectral modelling. To demonstrate this, we study the spectra of the 'normal' SN Ia 2005cf between 50 and 350 days after the explosion and identify most of the common SN Ia line features at post maximum epochs.

Twenty years after the explosion of SN 1987A, we are now able to observe the three-dimensional spatially resolved inner ejecta. Detailed mapping of newly synthesised material and its radioactive decay daughter products sheds light on the explosion mechanism. This may reveal the geometry of the explosion and its connection to the equatorial ring and the outer rings around SN 1987A. We have used integral field spectroscopy to image the supernova ejecta and the equatorial ring in the emission lines of [Si I]+[Fe II] and He I. The spectral information can be mapped into a radial velocity image revealing the expansion of the ejecta both as projected onto the sky and perpendicular to the sky plane. The inner ejecta are spatially resolved in a North-South direction and are clearly asymmetric. We argue that the bulk of the ejecta is situated in the same plane as defined by the equatorial ring and does not form a bipolar structure as has been suggested. The exact shape of the ejecta is modelled and we find that an elongated triaxial ellipsoid fits the observations best. From our spectral analyses of the ejecta spectrum we find that most of the He I, [Si I] and [Fe I-II] emission originates in the core material which has undergone explosive nucleosynthesis. The He I emission may be the result of alpha-rich freeze-out if the positron energy is deposited locally. Our observations clearly indicate a non-symmetric explosion mechanism for SN 1987A. The elongation and velocity asymmetries point towards a large-scale spatial non-spherical distribution as predicted in recent explosion models. The orientation of the ejecta in the plane of the equatorial ring argues against a jet-induced explosion through the poles due to stellar rotation.