We present the SuperNova Explosion Code (SNEC), an open-source Lagrangian code for the hydrodynamics and equilibrium-diffusion radiation transport in the expanding envelopes of supernovae. Given a model of a progenitor star, an explosion energy, and an amount and distribution of radioactive nickel, SNEC generates the bolometric light curve, as well as the light curves in different broad bands assuming black body emission. As a first application of SNEC, we consider the explosions of a grid of 15 M (at zero-age main sequence) stars whose hydrogen envelopes are stripped to different extents and at different points in their evolution. The resulting light curves exhibit plateaus with durations of ∼20 − 100 days if 1.5 − 2 M of hydrogen-rich material is left and no plateau if less hydrogen-rich material is left. If these shorter plateau lengths are not seen for Type IIP supernovae in nature, it suggests that, at least for zero-age main sequence masses 20 M , hydrogen mass loss occurs as an all or nothing process. This perhaps points to the important role binary interactions play in generating the observed mass-stripped supernovae (i.e., Type Ib/c events). These light curves are also unlike what is typically seen for Type IIL supernovae, arguing that simply varying the amount of mass loss cannot explain these events. The most stripped models begin to show double-peaked light curves similar to what is often seen for Type IIb supernovae, confirming previous work that these supernovae can come from progenitors that have a small amount of hydrogen and a radius of ∼ 500 R .

Recent modeling of hydrogen-rich Type II supernova (SN II) light curves suggests the presence of dense circumstellar material (CSM) surrounding the exploding progenitor stars. This has important implications for the activity and structure of massive stars near the end of their lives. Since previous work focused on just a few events, here we expand to a larger sample of twenty well-observed SNe II. For each event we are able to constrain the progenitor zero-age main-sequence (ZAMS) mass, explosion energy, and the mass and radial extent of dense CSM. We then study the distribution of each of these properties across the full sample of SNe. The inferred ZAMS masses are found to be largely consistent with a Salpeter distribution with minimum and maximum masses of 10.4 and 22.9 M , respectively. We also compare the individual ZAMS masses we measure with specific SNe II that have pre-explosion imaging to check their consistency. Our masses are generally comparable to or larger than the pre-explosion imaging masses, potentially helping ease the red supergiant problem. The explosion energies vary from (0.1 − 1.3) × 10 51 erg, and for ∼ 70% of the SNe we obtain CSM masses in the range between 0.18 − 0.83 M . We see a potential correlation between the CSM mass and explosion energy, which suggests that pre-explosion activity has a strong impact on the structure of the star. This may be important to take into account in future studies of the ability of the neutrino mechanism to explode stars. We also see a possible correlation between the CSM's radial extent and ZAMS mass, which could be related to the time with respect to explosion when the CSM is first generated.

A longstanding problem in the study of supernovae (SNe) has been the relationship between the TypeIIP and TypeIIL subclasses. Whether they come from distinct progenitors or they are from similar stars with some property that smoothly transitions from one class to another has been the subject of much debate. Here,using onedimensional radiation-hydrodynamic SN models, we showthat the multi-band light curves of SNe IIL are well fit by ordinary red supergiants surrounded by dense circumstellar material (CSM). The inferred extent of this material, coupled with a typical wind velocity of 10 100 km s 1 --, suggests enhanced activity by these stars during the last months to ∼years of their lives, which may be connected with advanced stages of nuclear burning. Furthermore, we find that,even for more plateau-like SNe,dense CSM provides a better fit to the first 20 days of their light curves, indicating that the presence of such material may be more widespread than previously appreciated. Here we choose to model the CSM with a wind-like density profile, but it is unclear whether this just generally represents some other mass distribution, such as a recent mass ejection, thick disk, or even inflated envelope material. Better understanding the exact geometry and density distribution of this material will be an important question for future studies.

We study gravitational waves (GWs) from a set of 2D multigroup neutrino radiation hydrodynamic simulations of core-collapse supernovae (CCSNe). Our goal is to systematize the current knowledge about the post-bounce CCSN GW signal and recognize the templatable features that could be used by the ground-based laser interferometers. We demonstrate that, starting from ∼400 ms after core bounce, the dominant GW signal represents the fundamental quadrupole (l = 2) oscillation mode (f-mode) of the proto–neutron star (PNS), which can be accurately reproduced by a linear perturbation analysis of the angle-averaged PNS profile. Before that, in the time interval between ∼200 and ∼400 ms after bounce, the dominant mode has two radial nodes and represents a g-mode. We associate the high-frequency noise in the GW spectrograms above the main signal with p-modes, while below the dominant frequency there is a region with very little power. The collection of models presented here summarizes the dependence of the CCSN GW signal on the progenitor mass, equation of state, many-body corrections to the neutrino opacity, and rotation. Weak dependence of the dominant GW frequency on the progenitor mass motivates us to provide a simple fit for it as a function of time, which can be used as a prior when looking for CCSN candidates in the LIGO data.