Massive stars have a strong impact on their surroundings, in particular when they produce a core-collapse supernova at the end of their evolution. In these proceedings, we review the general evolution of massive stars and their properties at collapse as well as the transition between massive and intermediate-mass stars. We also summarise the effects of metallicity and rotation. We then discuss some of the major uncertainties in the modelling of massive stars, with a particular emphasis on the treatment of convection in 1D stellar evolution codes. Finally, we present new 3D hydrodynamic simulations of convection in carbon burning and list key points to take from 3D hydrodynamic studies for the development of new prescriptions for convective boundary mixing in 1D stellar evolution codes.

We show how the dense shells of circumstellar gas immediately outside the red supergiants(RSGs) can affect the early optical light curves of Type II-P SNe taking the example of 2013ej. The peak in V, R and I bands, decline rate after peak and plateau length are found to be strongly influenced by the dense CSM formed due to enhanced mass loss during the oxygen and silicon burning stage of the progenitor. We find that the required explosion energy for the progenitors with CSM is reduced by almost a factor of 2.

Observational evidence from archival, pre-explosion images, suggests that progenitors of type-IIP SNe (SNe-IIP) have 8 ⩽ MP ⩽ 17 M⊙. However, the post-explosion temporal evolution of the event suggests that even in this mass range, the stellar evolutionary paths, the ensuing mass loss, and the eventual interaction of the supernova shock with the resulting CSM can show considerable diversity. Here we present the results from our program on multi-waveband (mainly optical) observations of SNe-IIP. Mass loss in their progenitors, with a massive and extended H-envelopes, is seen to occur via both strong stellar winds, or episodic mass ejections. Moreover, some type-IIP SNe also show unusually steep decline, characteristic of type-IIL (e.g. SN-IIP 2013ej). Our early and late-time spectrophotometry of these events shows CSM- shock interaction to varying degree among progenitors of comparable mass. Combined with X-ray data, our findings suggest that SNe-IIP progenitors can lose mass via strong stellar winds (e.g. SN2013ej, and SN2014cx), have episodic mass loss (SN2011ja), or have negligible mass loss (SN2012aw, SN2013ab).

Here we discuss the observational properties of supernovae exploding in extremely dense environments, namely Type IIn supernovae (SNe IIn). In SNe IIn, the surrounding environments play significant role in the supernovae energetics and evolution. Thus they are different than other classes of core collapse supernovae, whose energetics are not significantly altered by their environments. Though high density of medium is a prerequisite for radio and X-ray emission, less than 10% on SNe IIn are bright in these bands. This has important implications for their progenitor models. I will discuss the radio and X-ray observations of SNe IIn, which are crucial to unravel their complex environments. We also discuss some individual supernovae belonging to this class and discuss as to how they have refined our understanding of SNe IIn. Finally the importance of well sampled long term light curves in radio and X-ray bands cannot be stressed enough.

In their final stages, massive stars can show large eruptions which can resemble core-collapse IIn SNe. Here we present SN 2015bh in NGC 2770, a IIn/impostor, where archival data show variabilities for at least 21 years before the main event in 2015. Serendipitous spectra during an outburst are the only SN progenitor spectra available since SN 1987A and show an LBV with a fast, dense outflow. Analogues to SN 2015bh are SN 2009ip and SNhunt 248 while the SN 2000ch impostor could be equivalent to the outburst phase of SN 2015bh. It is currently unclear whether SN 2015bh (and SN 2009ip) were final core-collapse events. Alternatively, they might be large outbursts shedding the outer envelope and creating a Wolf-Rayet star in only a matter of decades. Future large-scale high-cadence surveys such as LSST will detect many more of these events and allow us a unique insight into the largely unknown late stages of massive stellar evolution.

Observing the supernovae (SNe) associated to the different types of gamma-ray bursts (GRBs) is one of the few means to study their progenitors. In the past years, it has become clear that GRB-like events are more heterogeneous than previously thought. There is a marked difference between long GRBs, which are produced by the collapse of very massive stars and are normally associated with broad-lined type Ic SNe, and short bursts, which occur when two compact objects merge and that, at least in some cases, can produce an associated kilonova. Moreover, the SNe associated with different sub-types of long GRBs are also seen to differ, especially those associated with ultra-long duration GRBs. To address this issue in a systematic way we started an observing programme in 2010 at the 10.4m GTC telescope. Here we present some results of our programme, including the detection of 12 new GRB-SNe. Highlights of our sample are the discovery of the first spectroscopic SN associated with a highly energetic (Eγ, iso ~ 1054 erg) “cosmological” burst (GRB 130427A), the study of the SN associated with a shock-breakout GRB (GRB 140606B) and the SN associated with the peculiar ultra-long GRB 101225A at z = 0.85. The sample includes also the follow-up of several short GRBs in search for kilonovae emission (GRB 130603B and GRB 160821B are important examples). Amongst our latest results we present the photometric and spectroscopic observations of the SNe associated with GRB 150818A and GRB 161219B.

We present the late-time Hubble Space Telescope observations of two Gamma Ray Burst (GRB) associated supernovae (SNe), GRB 030329/SN 2003dh and XRF 060218/SN 2006aj. Using the multi-color data up to ~320 days after the burst, we constrain the late-time decay nature of these SNe. The decay rates of SN 2003dh are steeper than SN 2006aj. A comparison with two other GRB SNe, GRB 980425/SN 1998bw and the SN associated with XRF 020903, shows that the decay rates of SN 2003dh are similar to XRF 020903 and those of SN 2006aj are similar to SN 1998bw. The late-time decay rates are steeper than the 56Co→56Fe radioactive decay rate indicating that there is some leakage of gamma-rays. We also compare the late-time decay rates of nine type Ic SNe, including the SNe of long GRBs, Ic broad lined and normal Ics. The decay rates of the SNe sample show a remarkable similarity in I band at late-times with a scatter of ~10%.

We investigate the relation between the emission properties of supernova shock breakout in the circumstellar matter (CSM) and the behavior of the shock. Using a Monte-Carlo method, we examine how the light curve and spectrum depends on the asphericity of the shock and bulk-Compton scattering, and compare the results with the observed properties of X-ray outburst (XRO) 080109/SN 2008D. We found that the rise and decay time of the X-ray light curve do not significantly depend on the degree of shock asphericity and the viewing angle in a steady and spherically symmetric CSM. The observed X-light curve and spectrum of XRO 080109 can be reproduced by considering the shock with a radial velocity of 60% of the speed of light, and the wind mass loss rate is about 5 × 10−4M⊙.

Supernova 1986J is almost the same age as SN 1987A, but was Type IIn, and likely had a massive progenitor. Located at 10 Mpc in NGC 891, it is one of the few supernovae whose radio emission can be resolved using VLBI. We present a new 5-GHz global-VLBI image of SN 1986J from 2014 as well as broadband VLA flux-density measurements. SN 1986J is unusual in that a compact synchrotron radio-emitting component appeared in the centre of the expanding shell of ejecta ~14 yr after the explosion, which now dominates the VLBI image. The central component is stationary to within the uncertainties (<570 km s−1), and it has a marginally resolved HWHM radius of (6.7−3.7+0.7) × 1016 cm. The shell has expanded with average v ≃ 5400 km s−1. The central component’s 5-GHz flux density is still increasing with time, and at present it has a 5-GHz νLν luminosity of ~4 × 1035 erg s−1, ~20 times that of the Crab Nebula. The central component may be due to a newly formed pulsar wind nebula, or an accreting black hole, or it may be due to interaction of the supernova shock with a highly structured environment left over from a progenitor which was in a close binary system. We discuss the newest observations and the constraints on its nature.

A large fraction of core-collapse supernovae are thought to result in the birth of a rotation-powered pulsar, which is later observable as a radio pulsar up to great ages. The birth properties of these pulsars, and in particular the distribution of their initial rotation periods, are however difficult to infer from studies of the radio pulsar population in our Galaxy. Yet the distributions of their birth properties is an important assumption for scenarios in which ultra-high-energy cosmic rays (UHECRs) originate in very young, extragalactic pulsars with short birth periods and/or high magnetic fields.

Using a model of the very young pulsar wind nebula’s dynamical and spectral evolution, with pulsar wind and accelerated particle parameters assumed similar to those inferred from modeling young pulsar wind nebulae (PWNe) in our Galaxy, we show that X-ray observations of supernovae, a few years to decades after the explosion, constitute a favored window to obtain meaningful constraints on the initial spin-down luminosity of the newly-formed pulsar. We examine the expected emerging PWN spectral component, taking into account the X-ray opacity of the expanding supernova ejecta, and find that it is typically best detectable in < 10 keV X-rays some years after the explosion. We use this framework to assess available X-ray observations and flux upper limits on supernovae, building on the work of Perna et al. (2008). We note that a resulting limit on spin-down luminosity corresponds univocally to a limit on the maximum magnetospheric acceleration potential, irrespective of the specific combination of magnetic field and rotation period that achieves it. We use available X-ray observations of supernovae to place constraints on the birth spin-down luminosity and period distribution of classical pulsars. We also examine the case of magnetars, born with much higher magnetic fields, and show that their much shorter initial spin-down time implies that any plausible signature of young magnetar wind nebulae can only be observed in harder X-ray or gamma-rays.

The radio non-detection of two Type Ia supernovae (SNe) SN 2011fe and SN 2014J has been modeled considering synchrotron radiation from shock accelerated electrons in the SN shock fronts. With 10% each of the bulk kinetic energy in electric and magnetic fields, a very low density of the medium around both the SNe has been estimated from the null detection of radio emission, around 1 and 4 years after the explosion of SNe 2014J and 2011fe, respectively. Keeping the fraction of energy in electrons fixed at 10%, a medium with particle density ~ 1cm−3 is found when 1% of the post shock energy is in magnetic fields. In case of a wind medium, the former predicts the mass loss rate Ṁ to be <10−9M⊙ yr−1, and the latter gives an upper limit ~10−9M⊙ yr−1, for wind velocity of 100 kms−1, for both the SNe. The tenuous media obtained from this study favor the double degenerate as well as a spin up/down model for both SNe 2011fe and 2014J.

The distributions of supernovae of different types and subtypes along the radius and in z coordinate of galaxies have been studied. We show that among SNe Ia in spiral galaxies, SNe Iax and Ia-norm have, respectively, the largest and smallest concentration to the center; the distributions of SNe Ia-91bg and Ia-91T are similar. A strong concentration of SNe Ibc to the central regions has been confirmed. In spiral galaxies, the supernovae of all types strongly concentrate to the galactic plane; the slight differences in scale height correlate with the extent to which the classes of supernovae are associated with star formation.

Supernovae (SNe) explode in environments that have been significantly modified by the SN progenitors. For core-collapse SNe, the massive progenitors ionize the ambient interstellar medium (ISM) via UV radiation and sweep the ambient ISM via fast stellar winds during the main sequence phase, replenish the surroundings with stellar material via slow winds during the luminous blue variable (LBV) or red supergiant (RSG) phase, and sweep up the circumstellar medium (CSM) via fast winds during the Wolf-Rayet (WR) phase. If a massive progenitor was in a close binary system, the binary interaction could have caused mass ejection in certain preferred directions, such as the orbital plane, and even bipolar outflow/jet. As a massive star finally explodes, the SN ejecta interacts first with the CSM that was ejected and shaped by the star itself. As the newly formed supernova remnant (SNR) expands further, it encounters interstellar structures that were shaped by the progenitor from earlier times. Therefore, the structure and evolution of a SNR is largely dependent on the initial mass and close binarity of the SN progenitor. The Large Magellanic Cloud (LMC) has an excellent sample of over 50 confirmed SNRs that are well resolved by Hubble Space Telescope, Chandra X-ray Observatory, and Spitzer Space Telescope. These multi-wavelength observations allow us to conduct stellar forensics in SNRs and understand the wide variety of morphologies and physical properties of SNRs observed.

The standard engine behind core-collapse supernovae is continuously evolving with increasingly detailed models. At this time, most simulations focus on an engine invoking turbulence above the proto-neutron star, sometimes termed the “convection-enhanced” engine. Here we review this engine and why it has become the standard for normal supernovae, focusing on a wide set of observations that provide insight into the supernova engine.

In an aspherical supernova explosion, shock emergence is not simultaneous and non-radial flows develop near the stellar surface. Oblique shock breakouts tend to be easily developed in compact progenitors like stripped-envelop core collapse supernovae. According to Matzner et al. (2013), non-spherical explosions develop non-radial flows that alters the observable emission and radiation of a supernova explosion. These flows can limit ejecta speed, change the distribution of matter and heat of the ejecta, suppress the breakout flash, and most importantly engender collisions outside the star. We construct a global numerical FLASH hydrodynamic simulation in a two dimensional spherical coordinate, focusing on the non-relativistic, adiabatic limit in a polytropic envelope to see how these fundamental differences affect the early light curve of core-collapse SNe.

To understand a wide variety of properties of young core-collapse supernova (CCSN) remnants being revealed by modern observations three-dimensional simulations of CCSNe starting from the initiation of the explosion until the expanding stellar debris transform into gaseous remnants are needed. We briefly review recent progress in modeling CCSNe on a long time scale. A current effort to model bolometric light curves based on 3D CCSN explosion models for comparison with observational data from SN 1987A is also discussed.

Core-collapse supernova explosions are driven by a central engine that converts a small fraction of the gravitational binding energy released during core collapse to outgoing kinetic energy. The suspected mode for this energy conversion is the neutrino mechanism, where a fraction of the neutrinos emitted from the newly formed protoneutron star are absorbed by and heat the matter behind the supernova shock. Accurate neutrino-matter interaction terms are crucial for simulating these explosions. In this proceedings for IAUS 331, SN 1987A, 30 years later, we explore several corrections to the neutrino-nucleon scattering opacity and demonstrate the effect on the dynamics of the core-collapse supernova central engine via two dimensional neutrino-radiation-hydrodynamics simulations. Our results reveal that the explosion properties are sensitive to corrections to the neutral-current scattering cross section at the 10-20% level, but only for densities at or above ~1012 g cm−3.

Hydrodynamical instabilities may either spin-up or down the pulsar formed in the collapse of a rotating massive star. Using numerical simulations of an idealized setup, we investigate the impact of progenitor rotation on the shock dynamics. The amplitude of the spiral mode of the Standing Accretion Shock Instability (SASI) increases with rotation only if the shock to the neutron star radii ratio is large enough. At large rotation rates, a corotation instability, also known as low-T/W, develops and leads to a more vigorous spiral mode. We estimate the range of stellar rotation rates for which pulsars are spun up or down by SASI. In the presence of a corotation instability, the spin-down efficiency is less than 30%. Given observational data, these results suggest that rapid progenitor rotation might not play a significant hydrodynamical role in the majority of core-collapse supernovae.

Extremely strong magnetic fields of the order of 1015G are required to explain the properties of magnetars, the most magnetic neutron stars. Such a strong magnetic field is expected to play an important role for the dynamics of core-collapse supernovae, and in the presence of rapid rotation may power superluminous supernovae and hypernovae associated to long gamma-ray bursts. The origin of these strong magnetic fields remains, however, obscure and most likely requires an amplification over many orders of magnitude in the protoneutron star. One of the most promising agents is the magnetorotational instability (MRI), which can in principle amplify exponentially fast a weak initial magnetic field to a dynamically relevant strength. We describe our current understanding of the MRI in protoneutron stars and show recent results on its dependence on physical conditions specific to protoneutron stars such as neutrino radiation, strong buoyancy effects and large magnetic Prandtl number.

We carried out high resolution simulations of weakly-magnetized core-collapse supernovae in two-dimensional axisymmetry in order to see the influence of the magnetic field and rotation on the explosion. We found that the magnetic field amplified by magnetorotational instability (MRI) has a great positive impact on the explosion by enhancing the neutrino heating, provided that the progenitor has large angular momentum close to the highest value found in stellar evolution calculations. We also found that even for progenitors neither involving strong magnetic flux nor large angular momentum, the magnetic field is greatly amplified by the convection aand rotation, and this leads to the boost of the explosion again by enhancing the neutrino heating.