• The density structure surrounding the iron core of a massive star when it dies is known to have a major effect on whether or not the star explodes. Here we repeat previous surveys of presupernova evolution with some important corrections to code physics and four to ten times better mass resolution in each star. The number of presupernova masses considered is also much larger. Over 4,000 models are calculated in the range from 12 to 60 Msun with varying mass loss rates. The core structure is not greatly affected by the increased spatial resolution. The qualitative patterns of compactness measures and their extrema are the same, but with the increased number of models, the scatter seen in previous studies is replaced by several localized branches. More physics-based analyses by Ertl et al. (2016) and M$\"u$eller et al. (2016) show these branches with less scatter than the single parameter characterization of O'Connor & Ott (2011). These branches are particularly apparent for stars in the mass ranges 14 - 19 Msun and 22 - 24 Msun. The multi-valued solutions are a consequence of interference between several carbon and oxygen burning shells during the late stages of evolution. For a relevant range of masses, whether a star explodes or not may reflect more the small, almost random differences in its late evolution than its initial mass. The large number of models allow statistically meaningful statements about the radius, luminosity, and effective temperatures of presupernova stars, their core structures, and their remnant mass distribution.
  • Recent observations have revealed an amazing diversity of extremely luminous supernovae, seemingly increasing in radiant energy without bound. We consider here the physical limits of what existing models can provide for the peak luminosity and total radiated energy for non-relativistic, isotropic stellar explosions. The brightest possible supernova is a Type I explosion powered by a sub-millisecond magnetar. Such models can reach a peak luminosity of $\rm 2\times10^{46}\ erg\ s^{-1}$ and radiate a total energy of $\rm 4 \times10^{52}\ erg$. Other less luminous models are also explored, including prompt hyper-energetic explosions in red supergiants, pulsational-pair instability supernovae, and pair-instability supernovae. Approximate analytic expressions and limits are given for each case. Excluding magnetars, the peak luminosity is near $\rm 1\times10^{44}\ erg\ s^{-1}$ for the brightest models. The corresponding limits on total radiated power are $\rm3 \times 10^{51}\ erg$ (Type I) and $\rm1 \times 10^{51}\ erg$ (Type II). A magnetar-based model for the recent transient event, ASASSN-15lh is presented that strains, but does not exceed the limits of what the model can provide.
  • We present 1-D non-Local-Thermodynamic-Equilibrium time-dependent radiative-transfer simulations for a large grid of supernovae (SNe) IIb/Ib/Ic that result from the terminal explosion of the mass donor in a close-binary system. Our sample covers ejecta masses $M_{\rm e}$ of 1.7$-$5.2M$_\odot$, kinetic energies $E_{\rm kin}$ of 0.6$-$5.0$\times$10$^{51}$erg, and $^{56}$Ni masses of 0.05$-$0.30M$_\odot$. We find a strong correlation between the $^{56}$Ni mass and the photometric properties at maximum, and between the rise time to bolometric maximum and the post-maximum decline rate. We confirm the small scatter in ($V-R$) at 10d past $R$-band maximum. The quantity $V_{\rm m} \equiv \sqrt{2E_{\rm kin}/M_{\rm e}}$ is comparable to the Doppler velocity measured from HeI 5875\AA\ at maximum in SNe IIb/Ib, although some scatter arises from the uncertain level of chemical mixing. The OI7772\AA\ line may be used for SNe Ic, but the correspondence deteriorates with higher ejecta mass/energy. We identify a temporal reversal of the Doppler velocity at maximum absorption in the $\sim$1.05$\mu$m feature in all models. The reversal is due to HeI alone and could serve as a test for the presence of helium in SNe Ic. Because of variations in composition and ionisation, the ejecta opacity shows substantial variations with both velocity and time. This is in part the origin of the offset between our model light curves and the predictions from the Arnett model.
  • We present 1-D non-Local-Thermodynamic-Equilibrium time-dependent radiative-transfer simulations for supernovae (SNe) of type IIb, Ib, and Ic that result from the terminal explosion of the mass donor in a close-binary system. Here, we select three ejecta with a total kinetic energy of ~1.2e51erg, but characterised by different ejecta masses (2-5Msun), composition, and chemical mixing. The type IIb/Ib models correspond to the progenitors that have retained their He-rich shell at the time of explosion. The type Ic model arises from a progenitor that has lost its helium shell, but retains 0.32Msun of helium in a CO-rich core of 5.11Msun. We discuss their photometric and spectroscopic properties during the first 2-3 months after explosion, and connect these to their progenitor and ejecta properties including chemical stratification. For these three models, Arnett's rule overestimates the 56Ni mass by ~50% while the procedure of Katz et al., based on an energy argument, yields a more reliable estimate. The presence of strong CI lines around 9000A prior to maximum is an indicator that the pre-SN star was under-abundant in helium. As noted by others, the 1.08micron feature is a complex blend of CI, MgII, and HeI lines, which makes the identification of He uncertain in SNe Ibc unless other HeI lines can be identified. Our models show little scatter in (V-R) colour 10d after R-band maximum. We also address a number of radiative transfer properties of SNe Ibc, including the notion of a photosphere, the inference of a representative ejecta expansion rate, spectrum formation, blackbody fits and "correction factors".
  • Merging carbon-oxygen (CO) white dwarfs are a promising progenitor system for Type Ia supernovae (SN Ia), but the underlying physics and timing of the detonation are still debated. If an explosion occurs after the secondary star is fully disrupted, the exploding primary will expand into a dense CO medium that may still have a disk-like structure. This interaction will decelerate and distort the ejecta. Here we carry out multi-dimensional simulations of ``tamped" SN Ia models, using both particle and grid-based codes to study the merger and explosion dynamics, and a radiative transfer code to calculate synthetic spectra and light curves. We find that post-merger explosions exhibit an hourglass-shaped asymmetry, leading to strong variations in the light curves with viewing angle. The two most important factors affecting the outcome are the scale-height of the disk, which depends sensitively on the binary mass ratio, and the total ${}^{56}$Ni yield, which is governed by the central density of the remnant core. The synthetic broadband light curves rise and decline very slowly, and the spectra generally look peculiar, with weak features from intermediate mass elements but relatively strong carbon absorption. We also consider the effects of the viscous evolution of the remnant, and show that a longer time delay between merger and explosion probably leads to larger ${}^{56}$Ni yields and more symmetrical remnants. We discuss the relevance of this class of aspherical ``tamped" SN Ia for explaining the class of ``super-Chandrasekhar'' SN Ia.
  • The success or failure of the neutrino-transport mechanism for producing a supernova in an evolved massive star is known to be sensitive not only to the mass of the iron core that collapses, but also to the density gradient in the silicon and oxygen shells surrounding that core. Here we study the systematics of a presupernova core's "compactness" (O'Connor & Ott 2011) as a function of the mass of the star and the physics used in its calculation. Fine-meshed surveys of presupernova evolution are calculated for stars from 15 to 65 Msun. The metallicity and the efficiency of semiconvection and overshoot mixing are both varied and bare carbon-oxygen cores are explored as well as full hydrogenic stars. Two different codes, KEPLER and MESA, are used for the study. A complex interplay of carbon and oxygen burning, especially in shells, can cause rapid variations in the compactness for stars of very nearly the same mass. On larger scales, the distribution of compactness with main sequence mass is found to be robustly non-monotonic, implying islands of "explodability", particularly around 8 to 20 Msun and 25 to 30 Msun. The carbon-oxygen (CO) core mass of a presupernova star is a better, though still ambiguous discriminant of its core structure than the main sequence mass.
  • Merging white dwarfs are a possible progenitor of Type Ia supernovae (SNe Ia). While it is not entirely clear if and when an explosion is triggered in such systems, numerical models suggest that a detonation might be initiated before the stars have coalesced to form a single compact object. Here we study such "peri-merger" detonations by means of numerical simulations, modeling the disruption and nucleosynthesis of the stars until the ejecta reach the coasting phase. Synthetic light curves and spectra are generated for comparison with observations. Three models are considered with primary masses 0.96 Msun, 1.06 Msun, and 1.20 Msun. Of these, the 0.96 Msun dwarf merging with an 0.81 Msun companion, with a Ni56 yield of 0.58 Msun, is the most promising candidate for reproducing common SNe Ia. The more massive mergers produce unusually luminous SNe Ia with peak luminosities approaching those attributed to "super-Chandrasekhar" mass SNe Ia. While the synthetic light curves and spectra of some of the models resemble observed SNe Ia, the significant asymmetry of the ejecta leads to large orientation effects. The peak bolometric luminosity varies by more than a factor of 2 with the viewing angle, and the velocities of the spectral absorption features are lower when observed from angles where the light curve is brightest. The largest orientation effects are seen in the ultraviolet, where the flux varies by more than an order of magnitude. Despite the large variation with viewing angle, the set of three models roughly obeys a width-luminosity relation, with the brighter light curves declining more slowly in the B-band. Spectral features due to unburned carbon from the secondary star are also seen in some cases.
  • The continuing difficulty of achieving a reliable explosion in simulations of core-collapse supernovae, especially for more massive stars, has led to speculation concerning the observable transients that might be produced if such a supernova fails. Even if a prompt outgoing shock fails to form in a collapsing presupernova star, one must still consider the hydrodynamic response of the star to the abrupt loss of mass via neutrinos as the core forms a protoneutron star. Following a suggestion by Nadezhin (1980), we calculate the hydrodynamical responses of typical supernova progenitor stars to the rapid loss of approximately 0.2 to 0.5 M_sun of gravitational mass from their centers. In a red supergiant star, a very weak supernova with total kinetic energy ~ 10^47 erg results. The binding energy of a large fraction of the hydrogen envelope before the explosion is of the same order and, depending upon assumptions regarding the neutrino loss rates, most of it is ejected. Ejection speeds are ~ 100 km/s and luminosities ~ 10^39 erg/s are maintained for about a year. A significant part of the energy comes from the recombination of hydrogen. The color of the explosion is extremely red and the events bear some similarity to "luminous red novae," but have much lower speeds.
  • Using non-LTE time-dependent radiative-transfer calculations, we study the impact of mixing and non-thermal processes associated with radioactive decay on SN IIb/Ib/Ic light curves (LCs) and spectra. Starting with short-period binary models of \leq5Msun He-rich stars (18-25Msun on the main-sequence), we produce 1.2B ejecta which we artificially mix to alter the chemical stratification. While the total 56Ni mass influences the LC peak, the spatial distribution of 56Ni, controlled by mixing processes, impacts both the multi-band LCs and spectra. With enhanced mixing, our synthetic LCs start their post-breakout re-brightening phase earlier, follow a more gradual rise to peak, appear redder, and fade faster after peak due to enhanced gamma-ray escape. Non-thermal electrons, crucial for the production of HeI lines, deposit a dominant fraction of their energy as heat. Because energy deposition is generally local well after the LC peak, the broad HeI lines characteristic of maximum-light SN IIb/Ib spectra require mixing that places 56Ni and helium nuclei to within a gamma-ray mean-free-path. This requirement indicates that SNe IIb and Ib most likely arise from the explosion of stripped-envelope massive stars (main-sequence masses \leq25Msun) that have evolved through mass-transfer in a binary system, rather than from more massive single WR stars. In contrast, the lack of HeI lines in SNe Ic may result from a variety of causes: A genuine helium deficiency; strongly-asymmetric mixing; weak mixing; or a more massive, perhaps single, progenitor characterized by a larger oxygen-rich core. Our models, subject to different mixing magnitudes, can produce a variety of SN types, including IIb, IIc, Ib, and Ic. As it is poorly constrained by explosion models, mixing challenges our ability to infer the progenitor and explosion properties of SNe IIb/Ib/Ic.
  • Xiaofeng Wang, Lifan Wang, Alexei V. Filippenko, Eddie Baron, Markus Kromer, Dennis Jack, Tianmeng Zhang, Greg Aldering, Pierre Antilogus, David Arnett, Dietrich Baade, Brian J. Barris, Stefano Benetti, Patrice Bouchet, Adam S. Burrows, Ramon Canal, Enrico Cappellaro, Raymond Carlberg, Elisa di Carlo, Peter Challis, Arlin Crotts, John I. Danziger, Massimo Della Valle, Michael Fink, Ryan J. Foley, Claes Fransson, Avishay Gal-Yam, Peter Garnavich, Chris L. Gerardy, Gerson Goldhaber, Mario Hamuy, Wolfgang Hillebrandt, Peter A. Hoeflich, Stephen T. Holland, Daniel E. Holz, John P. Hughes, David J. Jeffery, Saurabh W. Jha, Dan Kasen, Alexei M. Khokhlov, Robert P. Kirshner, Robert Knop, Cecilia Kozma, Kevin Krisciunas, Brian C. Lee, Bruno Leibundgut, Eric J. Lentz, Douglas C. Leonard, Walter H. G. Lewin, Weidong Li, Mario Livio, Peter Lundqvist, Dan Maoz, Thomas Matheson, Paolo Mazzali, Peter Meikle, Gajus Miknaitis, Peter Milne, Stefan Mochnacki, Ken'Ichi Nomoto, Peter E. Nugent, Elaine Oran, Nino Panagia, Saul Perlmutter, Mark M. Phillips, Philip Pinto, Dovi Poznanski, Christopher J. Pritchet, Martin Reinecke, Adam Riess, Pilar Ruiz-Lapuente, Richard Scalzo, Eric M. Schlegel, Brian Schmidt, James Siegrist, Alicia M. Soderberg, Jesper Sollerman, George Sonneborn, Anthony Spadafora, Jason Spyromilio, Richard A. Sramek, Sumner G. Starrfield, Louis G. Strolger, Nicholas B. Suntzeff, Rollin Thomas, John L. Tonry, Amedeo Tornambe, James W. Truran, Massimo Turatto, Michael Turner, Schuyler D. Van Dyk, Kurt Weiler, J. Craig Wheeler, Michael Wood-Vasey, Stan Woosley, Hitoshi Yamaoka
    Feb. 6, 2012 astro-ph.CO, astro-ph.HE
    We present ultraviolet (UV) spectroscopy and photometry of four Type Ia supernovae (SNe 2004dt, 2004ef, 2005M, and 2005cf) obtained with the UV prism of the Advanced Camera for Surveys on the Hubble Space Telescope. This dataset provides unique spectral time series down to 2000 Angstrom. Significant diversity is seen in the near maximum-light spectra (~ 2000--3500 Angstrom) for this small sample. The corresponding photometric data, together with archival data from Swift Ultraviolet/Optical Telescope observations, provide further evidence of increased dispersion in the UV emission with respect to the optical. The peak luminosities measured in uvw1/F250W are found to correlate with the B-band light-curve shape parameter dm15(B), but with much larger scatter relative to the correlation in the broad-band B band (e.g., ~0.4 mag versus ~0.2 mag for those with 0.8 < dm15 < 1.7 mag). SN 2004dt is found as an outlier of this correlation (at > 3 sigma), being brighter than normal SNe Ia such as SN 2005cf by ~0.9 mag and ~2.0 mag in the uvw1/F250W and uvm2/F220W filters, respectively. We show that different progenitor metallicity or line-expansion velocities alone cannot explain such a large discrepancy. Viewing-angle effects, such as due to an asymmetric explosion, may have a significant influence on the flux emitted in the UV region. Detailed modeling is needed to disentangle and quantify the above effects.
  • We present non-LTE time-dependent radiative-transfer simulations of supernova (SN) IIb/Ib/Ic spectra and light curves, based on ~1B-energy piston-driven ejecta, with and without 56Ni, produced from single and binary Wolf-Rayet (W-R) stars evolved at solar and sub-solar metallicities. Our bolometric light curves show a 10-day long post-breakout plateau with a luminosity of 1-5x10^7Lsun. In our 56Ni-rich models, with ~3Msun ejecta masses, this plateau precedes a 20-30-day long re-brightening phase initiated by the outward-diffusing heat wave powered by radioactive decay at depth. In low ejecta-mass models with moderate mixing, Gamma-ray leakage starts as early as ~50d after explosion and causes the nebular luminosity to steeply decline by ~0.02mag/d. Such signatures, which are observed in standard SNe IIb/Ib/Ic, are consistent with low-mass progenitors derived from a binary-star population. We propose that the majority of stars with an initial mass ~<20Msun yield SNe II-P if 'effectively" single, SNe IIb/Ib/Ic if part of a close binary system, and SN-less black holes if more massive. Our ejecta, with outer hydrogen mass fractions as low as ~>0.01 and a total hydrogen mass of ~>0.001Msun, yield the characteristic SN IIb spectral morphology at early times. However, by ~15d after the explosion, only Halpha may remain as a weak absorption feature. Our binary models, characterised by helium surface mass fractions of ~>0.85, systematically show HeI lines during the post-breakout plateau, irrespective of the 56Ni abundance. Synthetic spectra show a strong sensitivity to metallicity, which offers the possibility to constrain it directly from SN spectroscopic modelling.
  • The chief distinction between ordinary supernovae and long-soft gamma-ray bursts (GRBs) is the degree of differential rotation in the inner several solar masses when a massive star dies, and GRBs are rare mainly because of the difficulty achieving the necessary high rotation rate. Models that do provide the necessary angular momentum are discussed, with emphasis on a new single star model whose rapid rotation leads to complete mixing on the main sequence and avoids red giant formation. This channel of progenitor evolution also gives a broader range of masses than previous models, and allows the copious production of bursts outside of binaries and at high redshifts. However, even the production of a bare helium core rotating nearly at break up is not, by itself, a sufficient condition to make a gamma-ray burst. Wolf-Rayet mass loss must be low, and will be low in regions of low metallicity. This suggests that bursts at high redshift (low metallicity) will, on the average, be more energetic, have more time structure, and last longer than bursts nearby. Every burst consists of three components: a polar jet (~0.1 radian), high energy, subrelativistic mass ejection (~1 radian), and low velocity equatorial mass that can fall back after the initial explosion. The relative proportions of these three components can give a diverse assortment of supernovae and high energy transients whose properties may vary with redshift.
  • We study the detectability of primordial metal-free stars. Cosmological enrichment is a local process that takes place over an extended redshift range. While the duration of this transition depends on several unknown factors, in all cases late-forming metal-free stars are found in ~ 10^8 solar mass objects, which are just large enough to cool, but small enough to not be clustered near areas of previous star formation. We discuss the observational properties of these objects, some of which may have already been detected in ongoing Lyman-alpha surveys. If metal-free stars have masses between 140 and 260 solar masses, they are expected to end their lives as pair-production supernovae (PPSNe). Using the KEPLER code, we calculate a set of light curves for these objects. The peak luminosities of typical PPSNe are only slightly greater than those of Type Ia supernovae, but they remain bright much longer (~ 1 year) and have hydrogen lines. Ongoing supernova searches are able to place stringent limits on the fraction of very massive stars out to redshift 2. The planned Joint Dark Energy Mission satellite will be able to extend these constraints out to redshift 6
  • Those massive stars that, during their deaths, give rise to gamma-ray bursts (GRBs) must be endowed with an unusually large amount of angular momentum in their inner regions, one to two orders of magnitude greater than the ones that make common pulsars. Yet the inclusion of mass loss and angular momentum transport by magnetic torques during the precollapse evolution is known to sap the core of the necessary rotation. Here we explore the evolution of very rapidly rotating, massive stars, including stripped down helium cores that might result from mergers or mass transfer in a binary, and single stars that rotate unusually rapidly on the main sequence. For the highest possible rotation rates (about 400 km/s), a novel sort of evolution is encountered in which single stars mix completely on the main sequence, never becoming red giants. Such stars, essentially massive "blue stragglers", produce helium-oxygen cores that rotate unusually rapidly. Such stars might comprise roughly 1% of all stars above 10 solar masses and can, under certain circumstances retain enough angular momentum to make GRBs. Because this possibility is very sensitive to mass loss, GRBs will be much more probable in regions of low metallicity.
  • Nonrotating, zero metallicity stars with initial masses 140 < M < 260 solar masses are expected to end their lives as pair-production supernovae (PPSNe), in which an electron-positron pair-production instability triggers explosive nuclear burning. Interest in such stars has been rekindled by recent theoretical studies that suggest primordial molecular clouds preferentially form stars with these masses. Since metal enrichment is a local process, the resulting PPSNe could occur over a broad range of redshifts, in pockets of metal-free gas. Using the implicit hydrodynamics code KEPLER, we have calculated a set of PPSN light curves that addresses the theoretical uncertainties and allows us to assess observational strategies for finding these objects at intermediate redshifts. The peak luminosities of typical PPSNe are only slightly greater than those of Type Ia, but they remain bright much longer (~ 1 year) and have hydrogen lines. Ongoing supernova searches may soon be able to limit the contribution of these very massive stars to < 1% of the total star formation rate density out to z=2 which already provides useful constraints for theoretical models. The planned Joint Dark Energy Mission satellite will be able to extend these limits out to z=6.