We present the photometric and spectroscopic evolution of the Type II supernova (SN II) SN 2017ivv (also known as ASASSN-17qp). Located in an extremely faint galaxy (M$_r=-10.3$ mag), SN 2017ivv shows an unprecedented evolution during the two years o
f observations. At early times, the light curve shows a fast rise ($sim6-8$ days) to a peak of ${rm M}^{rm max}_{g}= -17.84$ mag, followed by a very rapid decline of $7.94pm0.48$ mag per 100 days in the $V-$band. The extensive photometric coverage at late phases shows that the radioactive tail has two slopes, one steeper than that expected from the decay of $^{56}$Co (between 100 and 350 days), and another slower (after 450 days), probably produced by an additional energy source. From the bolometric light curve, we estimated that the amount of ejected $^{56}$Ni is $sim0.059pm0.003$ M$odot$. The nebular spectra of SN 2017ivv show a remarkable transformation that allows the evolution to be split into three phases: (1) H$alpha$ strong phase ($<200$ days); (2) H$alpha$ weak phase (between 200 and 350 days); and (3) H$alpha$ broad phase ($>500$ days). We find that the nebular analysis favours a binary progenitor and an asymmetric explosion. Finally, comparing the nebular spectra of SN 2017ivv to models suggests a progenitor with a zero-age main-sequence mass of 15 -- 17 Msun.
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 observat
ions 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.
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 o
therwise 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.
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 el
ectromagnetic 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 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 f
ollowing 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.
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 nucl
eosynthesis 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.4pm1.2$ times solar, we find that burning of a fuel with neutron excess $eta approx 6times 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.
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.
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 l
ightcurve 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] lambdalambda 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 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 late
r 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.
