Pair-Instability Supernovae (PISNe) may signal the deaths of extremely massive stars in the local Universe or massive primordial stars after the end of the Cosmic Dark Ages. Hydrodynamic simulations of these explosions, performed in 1D, 2D, and 3D geometry, have revealed the strong dependence of mixing in the PISN ejecta on dimensionality. This chemical rearrangement is mainly driven by Rayleigh-Taylor instabilities that start to grow shortly after the collapse of the carbon-oxygen core. We investigate the effects of such mixing on the spectroscopic evolution of PISNe by post-processing explosion profiles with the radiation diffusion-equilibrium code SNEC and the implicit Monte Carlo-discrete diffusion Monte Carlo (IMC-DDMC) radiation transport code SuperNu. The first 3D radiation transport calculation of a PISN explosion is presented yielding viewing angle-dependent synthetic spectra and lightcurves. We find that while 2D and 3D mixing does not significantly affect the lightcurves of PISNe, their spectroscopic and color evolution is impacted. Strong features of intermediate mass elements dominated by silicon, magnesium and oxygen appear at different phases and reach different intensities depending on the extent of mixing in the silicon/oxygen interface of the PISN ejecta. On the other hand, we do not find a significant dependence of PISN lightcurves and spectra on viewing angle. Our results showcase the capabilities of SuperNu to handle 3D radiation transport and highlight the importance of modeling time-series of spectra in identifying PISNe with future missions.

A very massive star with a carbon-oxygen core in the range of $64$ M$_{\odot}<M_{\mathrm{CO}}<133$ M$_{\odot}$ is expected to undergo a very different kind of explosion known as a pair instability supernova. Pair instability supernovae are candidates for superluminous supernovae due to the prodigious amounts of radioactive elements they create. While the basic mechanism for the explosion is understood, how a star reaches a state is not, thus observations of a nearby pair-instability supernova would allow us to test current models of stellar evolution at the extreme of stellar masses. Much will be sought within the electromagnetic radiation we detect from such a supernova but we should not forget that the neutrinos from a pair-instability supernova contain unique signatures of the event that unambiguously identify this type of explosion. We calculate the expected neutrino flux at Earth from two, one-dimensional pair-instability supernova simulations which bracket the mass range of stars which explode by this mechanism taking into account the full time and energy dependence of the neutrino emission and the flavor evolution through the outer layers of the star. We calculate the neutrino signals in five different detectors chosen to represent present or near future designs. We find the more massive progenitors explode as pair-instability supernova which can easily be detected in multiple different neutrino detectors at the `standard' supernova distance of $10\;{\rm kpc}$ producing several events in DUNE, JUNE and SuperKamiokande, while the lightest progenitors only produce a handful of events (if any) in the same detectors. The proposed HyperKamiokande detector would detect neutrinos from a large pair-instability supernova as far as $\sim 50\;{\rm kpc}$ allowing it to reach the Megallanic Clouds and the several very high mass stars known to exist there.

In recent years, the viability of the pair-instability supernova (PISN) scenario for explaining superluminous supernovae has all but disappeared except for a few slowly-evolving examples. However, PISN are not predicted to be superluminous throughout the bulk of their mass range. In fact, it is more likely that the first PISN we see (if we have not seen one already) will not be superluminous. Here, we present hydrodynamic simulations of PISNe for four stellar models with unique envelope properties spanning the PISN mass range. In addition, we compute synthetic light curves for comparison with current and future observations. We also investigate, in the context of our most massive model, the prospect of mixing in the supernova ejecta alleviating discrepancies between current PISN models and the remaining superluminous candidate events. To this end, we present the first published 3D hydrodynamic simulations of PISNe. After achieving convergence between 1D, 2D, and 3D simulations we examine mixing in the supernova ejecta and its affect on the bolometric light curve. We observe slight deviations from spherical symmetry which increase with the number of dimensions. We find no significant effects on the bolometric light curve, however we conclude that mixing between the silicon and oxygen rich layers caused by the Rayleigh-Taylor instability may affect spectra.

With an increasing number of superluminous supernovae (SLSNe) discovered the question of their origin remains open and causes heated debates in the supernova community. Currently, there are three proposed mechanisms for SLSNe: (1) pair-instability supernovae (PISN), (2) magnetar-driven supernovae, and (3) models in which the supernova ejecta interacts with a circumstellar material ejected before the explosion. Based on current observations of SLSNe, the PISN origin has been disfavoured for a number of reasons. Many PISN models provide overly broad light curves and too reddened spectra, because of massive ejecta and a high amount of nickel. In the current study we re-examine PISN properties using progenitor models computed with the GENEC code. We calculate supernova explosions with FLASH and light curve evolution with the radiation hydrodynamics code STELLA. We find that high-mass models (200 and 250 solar masses) at relatively high metallicity (Z=0.001) do not retain hydrogen in the outer layers and produce relatively fast evolving PISNe Type I and might be suitable to explain some SLSNe. We also investigate uncertainties in light curve modelling due to codes, opacities, the nickel-bubble effect and progenitor structure and composition.

Since the emergence of the new class of extremely bright transients, super-luminous supernovae (SLSNe), three main mechanisms to power their light curves (LCs) have been discussed. They are the spin-down of a magnetar, interaction with circumstellar material, and the decay of large amounts of radioactive nickel in pair-instability supernovae (PISNe). Given the high degree of diversity seen within the class, it is possible that all three mechanisms are at work. PISN models can be self- consistently simulated from the main sequence phase of very massive stars (VMS) through to their explosion. These models robustly predict large amounts of radioactive nickel and thus very luminous SN events. However, PISN model LCs evolve more slowly than even the slowest evolving SLSNe. Multidimensional effects on the ejecta structure, specifically the mixing of radioactive nickel out to large radii, could alleviate this discrepancy with observation. Here we explore the multidimensional effects on the LC evolution by simulating the explosion phase in 1D, 2D, and 3D. We find that the ejecta from the multidimensional simulations have slightly shallower abundance gradients due to mixing at shell boundaries. We compute synthetic LCs whose shapes show no discernible differences due to the multidimensional effects.