Observational evidence from archival, pre-explosion images, suggests that progenitors of type-IIP SNe (SNe-IIP) have 8 ⩽ M P ⩽ 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).

Supernova remnants (SNRs) are powerful particle accelerators. As a supernova (SN) blast wave propagates through the circumstellar medium (CSM), electrons and protons scatter across the shock and gain energy by entrapment in the magnetic field. The accelerated particles generate further magnetic field fluctuations and local amplification, leading to cosmic ray production. The wealth of data from Supernova 1987A is providing a template of the SN-CSM interaction, and an important guide to the radio detection and identification of core-collapse SNe based on their spectral properties. Thirty years after the explosion, radio observations of SNR 1987A span from 70 MHz to 700 GHz. We review extensive observing campaigns with the Australia Telescope Compact Array (ATCA) and the Atacama Large Millimeter/submillimeter Array (ALMA), and follow-ups with other radio telescopes. Observations across the radio spectrum indicate rapid changes in the remnant morphology, while current ATCA and ALMA observations show that the SNR has entered a new evolutionary phase.

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.

Supernova remnants (SNRs) are diffuse extended sources characterized by a complex morphology and a non-uniform distribution of ejecta. Such a morphology reflects pristine structures and features of the progenitor supernova (SN) and the early interaction of the SN blast wave with the inhomogeneous circumstellar medium (CSM). Deciphering the observations of SNRs might open the possibility to investigate the physical properties of both the interacting ejecta and the shocked CSM. This requires accurate numerical models which describe the evolution from the SN explosion to the remnant development and which connect the emission properties of the remnants to the progenitor SNe. Here we show how multi-dimensional SN-SNR hydrodynamic models have been very effective in deciphering observations of SNR Cassiopeia A and SN 1987A, thus unveiling the structure of ejecta in the immediate aftermath of the SN explosion and constraining the 3D pre-supernova structure and geometry of the environment surrounding the progenitor SN.

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.

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 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.

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.

In this contribution, we present results from fully general-relativistic three-dimensional (3D) simulations of a non-rotating 15M ⊙ star using different nuclear equations of state (EOSs). We show that the SASI (standing-accretion-shock-instability) activity occurs much more vigorously in models with softer EOS. By performing detailed analysis of the gravitational-wave (GW) emission, we find a new GW signature that is produced predominantly by the SASI-induced downflows to the proto-neutron star. We discuss the detectability of the GW signal by performing a coherent network analysis where multiple detectors including LIGO Hanford, LIGO Livingston, VIRGO, and KAGRA are considered. We point out that the GW signal, whose typical frequency is in the best sensitivity range of the laser-interferometers, could potentially provide the live broadcast that pictures how the supernova shock is dancing in the core. The detection horizon of the signal is estimated as 2~3 kpc for the current generation detectors, which can extend up to ~100 kpc for the third generation detectors like Cosmic Explorer. We furthermore perform a correlation analysis between the SASI-modulated GW and neutrino signals. Our results show that the time correlation of the two signals becomes highest when we take into account the travel timescale of adverting material from the (average) neutrino-sphere to the proto-neutron star surface.

SN 1987A has been observed with the Chandra X-ray Observatory over the entire course of the mission. We have re-analyzed the archival data by constructing an empirical point spread function and reconstructing high-resolution images using a Bayesian multi-scale image reconstruction algorithm. We are able to resolve structure in the equatorial ring of SN 1987A with unprecedented detail, at scales of $\approx \frac{1}{4}$ arcsec. We describe how the point spread function is constructed, and the reconstruction method, and explore the evolution of the inner ring at different epochs and passbands.

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 have observed the oxygen-rich SNR 1E 0102.2-7219 with the integral field spectrograph WiFeS at Siding Spring Observatory and discovered sulfur-rich ejecta for the first time. Follow-up deep DDT observations with MUSE on the VLT (8100 s on source) reaching down to a noise level of ~5 × 10−20ergs−1cm−2Å−1spaxel−1 have led to the additional discovery of fast-moving hydrogen as well as argon-rich and chlorine-rich material. The detection of fast-moving hydrogen knots challenges the interpretation that the progenitor of 1E 0102 was a compact core of a Wolf-Rayet star that had shed its entire envelope. In addition to the detection of hydrogen and the products of oxygen-burning, this unprecedented sharp (0.2″ spaxel size at ~0.7″ seeing) and deep MUSE view of an oxygen-rich SNR in the Magellanic Clouds reveals further exciting discoveries, including [Fe xiv]λ5303 and [Fe xi]λ7892 emission, which we associate with the forward shock. We present this exciting data set and discuss some of its implications for the explosion mechanism and nucleosynthesis of the associated supernova.

As supernova remnants (SNRs) age, they become efficient cosmic ray accelerators at their outer shell shocks. The current paradigm for shock acceleration theory favors turbulent field environs in the proximity of these shocks, turbulence driven by current instabilities involving energetic ions. With the imminent prospect of dedicated X-ray polarimeters becoming a reality, the possibility looms of probing turbulence on scales that couple to the super-TeV electrons that emit X-rays. This paper presents model X-ray polarization signatures from energetic electrons moving in simulated MHD turbulence of varying levels of “chaos.” The emission volumes are finite slabs that represent the active regions of young SNR shells. We find that the turbulent field energy must be quite limited relative to that of the total field in order for the X-ray polarization degree to be as strong as the radio measures obtained in some remnants. Results presented are pertinent to the planned IXPE and XIPE polarimeters.

Fostered by the possibilities of multi-dimensional computational modeling, in particular the advent of three-dimensional (3D) simulations, our understanding of the neutrino-driven explosion mechanism of core-collapse supernovae (SNe) has experienced remarkable progress over the past decade. First self-consistent, first-principle models have shown successful explosions in 3D, and even failed cases may be cured by moderate changes of the microphysics inside the neutron star (NS), better grid resolution, or more detailed progenitor conditions at the onset of core collapse, in particular large-scale perturbations in the convective Si and O burning shells. 3D simulations have also achieved to follow neutrino-driven explosions continuously from the initiation of the blast wave, through the shock breakout from the progenitor surface, into the radioactively powered evolution of the SN, and towards the free expansion phase of the emerging remnant. Here we present results from such simulations, which form the basis for direct comparisons with observations of SNe and SN remnants in order to derive constraints on the still disputed explosion mechanism. It is shown that predictions based on hydrodynamic instabilities and mixing processes associated with neutrino-driven explosions yield good agreement with measured NS kicks, light-curve properties of SN 1987A and asymmetries of iron and 44Ti distributions observed in SN 1987A and Cassiopeia A.

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.

HESS J1731−347 a.k.a. SNR G353.6−0.7 is one of the five known very high energy (VHE, Energy > 0.1 TeV) shell-type supernova remnants. We carried out Giant Metrewave Radio Telescope (GMRT) observations of this TeV SNR in 1390, 610 and 325 MHz bands. We detected the 325 and 610 MHz radio counterparts of the SNR G353.6−0.7 (Nayana et al. 2017). We also determined the spectral indices of individual filaments and our values are consistent with the non-thermal radio emission. We compared the radio morphology with that of VHE emission. The peak in radio emission corresponds to the faintest feature in the VHE emission. We explain this anti-correlated emission in a possible leptonic origin of the VHE γ-rays.

The data we are receiving from galactic cosmic rays are reaching an unprecedented precision, over very wide energy ranges. Nevertheless, many problems are still open, while new ones seem to appear when data happen to be redundant. We will discuss some paths to possible progress in the theoretical modeling and experimental exploration of the galactic cosmic radiation.

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.

The dust produced by supernovae is an important topic for understanding supernova physics and the chemical evolution of galaxies. Recent ALMA observations of SN 1987A have allowed us to peer into the inner ejecta to the cool dust, with spatial resolution from 0.″3 at ~300 GHz down to 0.″09 at ~680 GHz – an improvement over the previous 300 GHz Cycle 0 observations at 0.″69. Comparison of the dust location and morphology with other multiwavelength emission presents an interesting picture of the role dust plays in the ejecta. The mm-FIR SED is compared to radiative models to study the dust composition 30 years after the initial explosion. Fits to the ring emission also probe the drift of the center of the system over time.