Showing 1–50 of 225 results for author: Woosley, S

Double "acct": a distinct double-peaked supernova matching pulsational pair-instability models

Authors: C. R. Angus, S. E. Woosley, R. J. Foley, M. Nicholl, V. A. Villar, K. Taggart, M. Pursiainen, P. Ramsden, S. Srivastav, H. F. Stevance, T. Moore, K. Auchettl, W. B. Hoogendam, N. Khetan, S. K. Yadavalli, G. Dimitriadis, A. Gagliano, M. R. Siebert, A. Aamer, T. de Boer, K. C. Chambers, A. Clocchiatti, D. A. Coulter, M. R. Drout, D. Farias , et al. (13 additional authors not shown)

Abstract: We present multi-wavelength data of SN2020acct, a double-peaked stripped-envelope supernova (SN) in NGC2981 at ~150 Mpc. The two peaks are temporally distinct, with maxima separated by 58 rest-frame days, and a factor of 20 reduction in flux between. The first is luminous (M$_{r}$ = -18.00 $\pm$ 0.02 mag), blue (g - r = 0.27 $\pm$ 0.03 mag), and displays spectroscopic signatures of interaction wit… ▽ More We present multi-wavelength data of SN2020acct, a double-peaked stripped-envelope supernova (SN) in NGC2981 at ~150 Mpc. The two peaks are temporally distinct, with maxima separated by 58 rest-frame days, and a factor of 20 reduction in flux between. The first is luminous (M$_{r}$ = -18.00 $\pm$ 0.02 mag), blue (g - r = 0.27 $\pm$ 0.03 mag), and displays spectroscopic signatures of interaction with hydrogen-free circumstellar material. The second peak is fainter (M$_{r}$ = -17.29 $\pm$ 0.03 mag), and spectroscopically similar to an evolved stripped-envelope SNe, with strong blended forbidden [Ca II] and [O II] features. No other known double-peak SN exhibits a light curve similar to that of SN 2020acct. We find the likelihood of two individual SNe occurring in the same star-forming region within that time to be highly improbable, while an implausibly fine-tuned configuration would be required to produce two SNe from a single binary system. We find that the peculiar properties of SN2020acct match models of pulsational pair instability (PPI), in which the initial peak is produced by collisions of shells of ejected material, shortly followed by a terminal explosion. Pulsations from a star with a 72 M$_{\odot}$ helium core provide an excellent match to the double-peaked light curve. The local galactic environment has a metallicity of 0.4 Z$_{\odot}$, a level where massive single stars are not expected retain enough mass to encounter the PPI. However, late binary mergers or a low-metallicity pocket may allow the required core mass. We measure the rate of SN 2020acct-like events to be $<3.3\times10^{-8}$ Mpc$^{-3}$ yr$^{-1}$ at z = 0.07, or <0.1% of the total core-collapse SN rate. △ Less

Submitted 3 September, 2024; originally announced September 2024.

Comments: 31 pages, 14 figures. Submitted to ApJL, comments welcome

A cosmic formation site of silicon and sulphur revealed by a new type of supernova explosion

Authors: Steve Schulze, Avishay Gal-Yam, Luc Dessart, Adam A. Miller, Stan E. Woosley, Yi Yang, Mattia Bulla, Ofer Yaron, Jesper Sollerman, Alexei V. Filippenko, K-Ryan Hinds, Daniel A. Perley, Daichi Tsuna, Ragnhild Lunnan, Nikhil Sarin, Sean J. Brennan, Thomas G. Brink, Rachel J. Bruch, Ping Chen, Kaustav K. Das, Suhail Dhawan, Claes Fransson, Christoffer Fremling, Anjasha Gangopadhyay, Ido Irani , et al. (25 additional authors not shown)

Abstract: The cores of stars are the cosmic furnaces where light elements are fused into heavier nuclei. The fusion of hydrogen to helium initially powers all stars. The ashes of the fusion reactions are then predicted to serve as fuel in a series of stages, eventually transforming massive stars into a structure of concentric shells. These are composed of natal hydrogen on the outside, and consecutively hea… ▽ More The cores of stars are the cosmic furnaces where light elements are fused into heavier nuclei. The fusion of hydrogen to helium initially powers all stars. The ashes of the fusion reactions are then predicted to serve as fuel in a series of stages, eventually transforming massive stars into a structure of concentric shells. These are composed of natal hydrogen on the outside, and consecutively heavier compositions inside, predicted to be dominated by helium, carbon/oxygen, oxygen/neon/magnesium, and oxygen/silicon/sulphur. Silicon and sulphur are fused into inert iron, leading to the collapse of the core and either a supernova explosion or the direct formation of a black hole. Stripped stars, where the outer hydrogen layer has been removed and the internal He-rich layer (in Wolf-Rayet WN stars) or even the C/O layer below it (in Wolf-Rayet WC/WO stars) are exposed, provide evidence for this shell structure, and the cosmic element production mechanism it reflects. The types of supernova explosions that arise from stripped stars embedded in shells of circumstellar material (most notably Type Ibn supernovae from stars with outer He layers, and Type Icn supernovae from stars with outer C/O layers) confirm this scenario. However, direct evidence for the most interior shells, which are responsible for the production of elements heavier than oxygen, is lacking. Here, we report the discovery of the first-of-its-kind supernova arising from a star peculiarly stripped all the way to the silicon and sulphur-rich internal layer. Whereas the concentric shell structure of massive stars is not under debate, it is the first time that such a thick, massive silicon and sulphur-rich shell, expelled by the progenitor shortly before the SN explosion, has been directly revealed. △ Less

Submitted 3 September, 2024; originally announced September 2024.

Comments: 48 pages, 12 figures and 10 tables. Submitted to a high-impact journal. The reduced spectra and photometry will be made available via the journal webpage and the WISeREP archive after the acceptance of the paper

Nebular Nitrogen Line Emission in Stripped-Envelope Supernovae -- a New Progenitor Mass Diagnostic

Authors: Stan Barmentloo, Anders Jerkstrand, Koichi Iwamoto, Izumi Hachisu, Ken'ichi Nomoto, Jesper Sollerman, Stan Woosley

Abstract: Nitrogen is produced by CNO-cycling in massive stars, and can be ejected in significant amounts in supernova explosions. While in H-rich SNe, its [\ion{N}{II}] 6548, 6583 emission becomes obscured by strong H$α$, in explosions of He stars, this nitrogen emission becomes more visible. We here explore the formation of this line, using the \texttt{SUMO} code to compute spectra for a grid of 1D models… ▽ More Nitrogen is produced by CNO-cycling in massive stars, and can be ejected in significant amounts in supernova explosions. While in H-rich SNe, its [\ion{N}{II}] 6548, 6583 emission becomes obscured by strong H$α$, in explosions of He stars, this nitrogen emission becomes more visible. We here explore the formation of this line, using the \texttt{SUMO} code to compute spectra for a grid of 1D models with parameterized mixing informed from new 2D simulations. Because the mass fraction of nitrogen in the ejecta decreases with larger He core masses, as more of the He/N zone gets processed by shell helium burning and is lost to winds, the [\ion{N}{II}] luminosity relative to the overall optical flux probes the He core mass. By comparing to large samples of data, we find that low-mass He cores ($M_{\rm preSN}\lesssim\ 3\ M_\odot$) are exclusively associated with Type IIb SNe, with the exception of Type Ib SN 2007Y. Seeing no strong nitrogen emission in other Type Ib SNe, the implication is either an origin from low-mass stars with the He/N layer (but not the He/C) layer peeled away, or from higher-mass He cores. We also see no clear nitrogen emission in Type Ic SNe. We discuss the diagnostic potential of this new line metric, and also dependencies on mass-loss-rate and metallicity. △ Less

Submitted 13 March, 2024; originally announced March 2024.

Comments: 17 pages, 11 figures + 12 pages appendix

Modeling of the nebular-phase spectral evolution of stripped-envelope supernovae. New grids from 100 to 450 days

Authors: Luc Dessart, D. John Hillier, S. E. Woosley, Hanindyo Kuncarayakti

Abstract: We present an extended grid of multi-epoch 1D nonlocal thermodynamic equilibrium radiative transfer calculations for nebular-phase Type Ibc supernovae (SNe) from He-star explosions. Compared to Dessart+21, we study the spectral evolution from 100 to about 450d and augment the model set with progenitors that were evolved without wind mass loss. Models with the same final, preSN mass have similar yi… ▽ More We present an extended grid of multi-epoch 1D nonlocal thermodynamic equilibrium radiative transfer calculations for nebular-phase Type Ibc supernovae (SNe) from He-star explosions. Compared to Dessart+21, we study the spectral evolution from 100 to about 450d and augment the model set with progenitors that were evolved without wind mass loss. Models with the same final, preSN mass have similar yields and produce essentially the same emergent spectra. Hence, the uncertain progenitor mass loss history compromises the inference of the initial, main sequence mass. This shortcoming does not affect Type IIb SNe. However, our 1D models with a different preSN mass tend to yield widely different spectra, as seen through variations in the strong emission lines due to [NII]6548-6583, [OI]6300-6364, [CaII]7291-7323, [NiII]7378, and the forest of FeII lines below 5500A. At the lower mass end, the ejecta are He rich and at 100d cool through HeI, NII, CaII, and FeII lines, with NII and FeII dominating at 450d. These models, associated with He giants, conflict with observed SNe Ib, which typically lack strong NII emission. Instead they may lead to SNe Ibn or, because of additional stripping by a companion star, ultra-stripped SNe Ic. In contrast, for higher preSN masses, the ejecta are progressively He poor and cool at 100d through OI, CaII, and FeII lines, with OI and CaII dominating at 450d. Nonuniform, aspherical, large-scale mixing rather than composition differences likely determines the SN type at intermediate preSN masses. Variations in clumping, mixing, as well as departures from spherical symmetry would increase the spectral diversity but also introduce additional degeneracies. More robust predictions from spectral modeling require a careful attention to the initial conditions informed by physically-consistent 3D explosion models [abridged]. △ Less

Submitted 21 June, 2023; originally announced June 2023.

Comments: Accepted for publication in A&A

Journal ref: A&A 677, A7 (2023)

Supernova 2020wnt: An Atypical Superluminous Supernova with a Hidden Central Engine

Authors: Samaporn Tinyanont, Stan E. Woosley, Kirsty Taggart, Ryan J. Foley, Lin Yan, Ragnhild Lunnan, Kyle W. Davis, Charles D. Kilpatrick, Matthew R. Siebert, Steve Schulze, Chris Ashall, Ting-Wan Chen, Kishalay De, Georgios Dimitriadis, Dillon Z. Dong, Christoffer Fremling, Alexander Gagliano, Saurabh W. Jha, David O. Jones, Mansi M. Kasliwal, Hao-Yu Miao, Yen-Chen Pan, Daniel A. Perley, Vikram Ravi, César Rojas-Bravo , et al. (12 additional authors not shown)

Abstract: We present observations of a peculiar hydrogen- and helium-poor stripped-envelope (SE) supernova (SN) 2020wnt, primarily in the optical and near-infrared (near-IR). Its peak absolute bolometric magnitude of -20.9 mag and a rise time of 69~days are reminiscent of hydrogen-poor superluminous SNe (SLSNe~I), luminous transients potentially powered by spinning-down magnetars. Before the main peak, ther… ▽ More We present observations of a peculiar hydrogen- and helium-poor stripped-envelope (SE) supernova (SN) 2020wnt, primarily in the optical and near-infrared (near-IR). Its peak absolute bolometric magnitude of -20.9 mag and a rise time of 69~days are reminiscent of hydrogen-poor superluminous SNe (SLSNe~I), luminous transients potentially powered by spinning-down magnetars. Before the main peak, there is a brief peak lasting <10 days post-explosion, likely caused by interaction with circumstellar medium (CSM) ejected ~years before the SN explosion. The optical spectra near peak lack a hot continuum and OII absorptions, which are signs of heating from a central engine; they quantitatively resemble those of radioactivity-powered H/He-poor Type Ic SESNe. At ~1 year after peak, nebular spectra reveal a blue pseudo-continuum and narrow OI recombination lines associated with magnetar heating. Radio observations rule out strong CSM interactions as the dominant energy source at +266 days post peak. Near-IR observations at +200-300 day reveal carbon monoxide and dust formation, which causes a dramatic optical light curve dip. Pair-instability explosion models predict slow light curve and spectral features incompatible with observations. SN 2020wnt is best explained as a magnetar-powered core-collapse explosion of a 28 Msun pre-SN star. The explosion kinetic energy is significantly larger than the magnetar energy at peak, effectively concealing the magnetar-heated inner ejecta until well after peak. SN 2020wnt falls into a continuum between normal SNe Ic and SLSNe I and demonstrates that optical spectra at peak alone cannot rule out the presence of a central engine. △ Less

Submitted 30 November, 2022; originally announced December 2022.

Comments: Submitted to ApJ

StaNdaRT: A repository of standardized test models and outputs for supernova radiative transfer

Authors: Stéphane Blondin, Sergei Blinnikov, Fionntan P. Callan, Christine E. Collins, Luc Dessart, Wesley Even, Andreas Flörs, Andrew G. Fullard, D. John Hillier, Anders Jerkstrand, Daniel Kasen, Boaz Katz, Wolfgang Kerzendorf, Alexandra Kozyreva, Jack O'Brien, Ezequiel A. Pássaro, Nathaniel Roth, Ken J. Shen, Luke Shingles, Stuart A. Sim, Jaladh Singhal, Isaac G. Smith, Elena Sorokina, Victor P. Utrobin, Christian Vogl , et al. (4 additional authors not shown)

Abstract: We present the first results of a comprehensive supernova (SN) radiative-transfer (RT) code-comparison initiative (StaNdaRT), where the emission from the same set of standardized test models is simulated by currently-used RT codes. A total of ten codes have been run on a set of four benchmark ejecta models of Type Ia supernovae. We consider two sub-Chandrasekhar-mass ($M_\mathrm{tot} = 1.0$ M… ▽ More We present the first results of a comprehensive supernova (SN) radiative-transfer (RT) code-comparison initiative (StaNdaRT), where the emission from the same set of standardized test models is simulated by currently-used RT codes. A total of ten codes have been run on a set of four benchmark ejecta models of Type Ia supernovae. We consider two sub-Chandrasekhar-mass ($M_\mathrm{tot} = 1.0$ M$_\odot$) toy models with analytic density and composition profiles and two Chandrasekhar-mass delayed-detonation models that are outcomes of hydrodynamical simulations. We adopt spherical symmetry for all four models. The results of the different codes, including the light curves, spectra, and the evolution of several physical properties as a function of radius and time, are provided in electronic form in a standard format via a public repository. We also include the detailed test model profiles and several python scripts for accessing and presenting the input and output files. We also provide the code used to generate the toy models studied here. In this paper, we describe in detail the test models, radiative-transfer codes and output formats and provide access to the repository. We present example results of several key diagnostic features. △ Less

Submitted 15 April, 2023; v1 submitted 23 September, 2022; originally announced September 2022.

Comments: Accepted for publication in A&A. 27 pages, 12 figures (v4: updated to match published version). The ejecta models and output files from the simulations are available at https://github.com/sn-rad-trans/data1

Journal ref: A&A 668, A163 (2022)

Physics-Informed Machine Learning for Modeling Turbulence in Supernovae

Authors: Platon I. Karpov, Chengkun Huang, Iskandar Sitdikov, Chris L. Fryer, Stan Woosley, Ghanshyam Pilania

Abstract: Turbulence plays an important role in astrophysical phenomena, including core-collapse supernovae (CCSN), but current simulations must rely on subgrid models since direct numerical simulation (DNS) is too expensive. Unfortunately, existing subgrid models are not sufficiently accurate. Recently, Machine Learning (ML) has shown an impressive predictive capability for calculating turbulence closure.… ▽ More Turbulence plays an important role in astrophysical phenomena, including core-collapse supernovae (CCSN), but current simulations must rely on subgrid models since direct numerical simulation (DNS) is too expensive. Unfortunately, existing subgrid models are not sufficiently accurate. Recently, Machine Learning (ML) has shown an impressive predictive capability for calculating turbulence closure. We have developed a physics-informed convolutional neural network (CNN) to preserve the realizability condition of Reynolds stress that is necessary for accurate turbulent pressure prediction. The applicability of the ML subgrid model is tested here for magnetohydrodynamic (MHD) turbulence in both the stationary and dynamic regimes. Our future goal is to utilize this ML methodology (available on GitHub) in the CCSN framework to investigate the effects of accurately-modeled turbulence on the explosion of these stars. △ Less

Submitted 9 August, 2022; v1 submitted 17 May, 2022; originally announced May 2022.

Comments: For our ML algorithm on GitHub, see https://github.com/pikarpov-LANL/Sapsan/wiki/Estimators\#physics-informed-cnn-for-turbulence-modeling

SN 1961V: A Pulsational Pair-Instability Supernova

Authors: S. E. Woosley, Nathan Smith

Abstract: We explore a variety of models in which SN~1961V, one of the most enigmatic supernovae (SNe) ever observed, was a pulsational pair-instability supernova (PPISN). Successful models reproduce the bolometric light curve of the principal outburst and, in some cases, the emission one year before and several years afterward. All models have helium-rich ejecta, bulk hydrogenic velocities near 2000 km s… ▽ More We explore a variety of models in which SN~1961V, one of the most enigmatic supernovae (SNe) ever observed, was a pulsational pair-instability supernova (PPISN). Successful models reproduce the bolometric light curve of the principal outburst and, in some cases, the emission one year before and several years afterward. All models have helium-rich ejecta, bulk hydrogenic velocities near 2000 km s$^{-1}$, and total kinetic energies from 4 to 8 $\times 10^{50}$ erg. Each eventually leaves behind a black hole remnant. Three subclasses of PPISN models are explored, each with two different choices of carbon abundance following helium burning. Carbon is an important parameter because shell carbon burning can weaken the explosion. The three subclasses correspond to situations where SN~1961V and its immediate afterglow were: a) a single event; b) the first of two or more pulsational events separated by decades or centuries; or c) the latter stages of a complex explosion that had already been going on for a year or more. For the low carbon case, the main sequence mass for SN~1961V's progenitor would have been 100 to 115 \Msun; its pre-SN helium core mass was 45 to 52 \Msun; and the final black hole mass, 40 to 45 \Msun. For the high-carbon case, these values are increased by roughly 20 to 25\%. In some PPISN models, a $\sim10^{40}$ erg s$^{-1}$ star-like object could still be shining at the site of SN~1961V, but it has more likely been replaced by a massive accreting black hole. △ Less

Submitted 12 May, 2022; originally announced May 2022.

Comments: 18 pages, 11 figures, 3 tables; submitted to Astrophysical Journal, 5/9/22

A Carbon/Oxygen-dominated Atmosphere Days After Explosion for the "Super-Chandrasekhar" Type Ia SN 2020esm

Authors: Georgios Dimitriadis, Ryan J. Foley, Nikki Arendse, David A. Coulter, Wynn V. Jacobson-Galán, Matthew R. Siebert, Luca Izzo, David O. Jones, Charles D. Kilpatrick, Yen-Chen Pan, Kirsty Taggart, Katie Auchettl, Christa Gall, Jens Hjorth, Daniel Kasen, Anthony L. Piro, Sandra I. Raimundo, Enrico Ramirez-Ruiz, Armin Rest, Jonathan J. Swift, Stan E. Woosley

Abstract: Seeing pristine material from the donor star in a Type Ia supernova (SN Ia) explosion can reveal the nature of the binary system. In this paper, we present photometric and spectroscopic observations of SN 2020esm, one of the best-studied SNe of the class of "super-Chandrasekhar" SNe Ia (SC SNe Ia), with data obtained $-12$ to +360 days relative to peak brightness, obtained from a variety of ground… ▽ More Seeing pristine material from the donor star in a Type Ia supernova (SN Ia) explosion can reveal the nature of the binary system. In this paper, we present photometric and spectroscopic observations of SN 2020esm, one of the best-studied SNe of the class of "super-Chandrasekhar" SNe Ia (SC SNe Ia), with data obtained $-12$ to +360 days relative to peak brightness, obtained from a variety of ground- and space-based telescopes. Initially misclassified as a Type II supernova, SN 2020esm peaked at $M_{B} = -19.9$ mag, declined slowly ($Δm_{15}(B) = 0.92$ mag), and had particularly blue UV and optical colors at early times. Photometrically and spectroscopically, SN 2020esm evolved similarly to other SC SNe Ia, showing the usual low ejecta velocities, weak intermediate mass elements (IMEs), and the enhanced fading at late times, but its early spectra are unique. Our first few spectra (corresponding to a phase of $\gtrsim$10~days before peak) reveal a nearly-pure carbon/oxygen atmosphere during the first days after explosion. This composition can only be produced by pristine material, relatively unaffected by nuclear burning. The lack of H and He may further indicate that SN 2020esm is the outcome of the merger of two carbon/oxygen white dwarfs (WDs). Modeling its bolometric light curve, we find a $^{56}$Ni mass of $1.23^{+0.14}_{-0.14}$ M$_{\odot}$ and an ejecta mass of $1.75^{+0.32}_{-0.20}$ M$_{\odot}$, in excess of the Chandrasekhar mass. Finally, we discuss possible progenitor systems and explosion mechanisms of SN 2020esm and, in general, the SC SNe Ia class. △ Less

Submitted 18 December, 2021; originally announced December 2021.

Comments: 21 pages, 15 figures, accepted for publication in APJ

Pulsational pair-instability supernovae: gravitational collapse, black-hole formation, and beyond

Authors: Ninoy Rahman, Hans-Thomas Janka, Georg Stockinger, Stan Woosley

Abstract: We investigate the final collapse of rotating and non-rotating pulsational pair-instability supernova progenitors with zero-age-main-sequence masses of 60, 80, and 115$\mathrm{M}_\odot$ and iron cores between 2.37$\mathrm{M}_\odot$ and 2.72$\mathrm{M}_\odot$ by 2D hydrodynamics simulations. Using the general relativistic NADA-FLD code with energy-dependent three-flavor neutrino transport by flux-l… ▽ More We investigate the final collapse of rotating and non-rotating pulsational pair-instability supernova progenitors with zero-age-main-sequence masses of 60, 80, and 115$\mathrm{M}_\odot$ and iron cores between 2.37$\mathrm{M}_\odot$ and 2.72$\mathrm{M}_\odot$ by 2D hydrodynamics simulations. Using the general relativistic NADA-FLD code with energy-dependent three-flavor neutrino transport by flux-limited diffusion allows us to follow the evolution beyond the moment when the transiently forming neutron star (NS) collapses to a black hole (BH), which happens within 350$-$580 ms after bounce in all cases. Because of high neutrino luminosities and mean energies, neutrino heating leads to shock revival within $\lesssim$250 ms post bounce in all cases except the rapidly rotating 60$\mathrm{M}_\odot$ model. In the latter case, centrifugal effects support a 10% higher NS mass but reduce the radiated neutrino luminosities and mean energies by $\sim$20% and $\sim$10%, respectively, and the neutrino-heating rate by roughly a factor of two compared to the non-rotating counterpart. After BH formation, the neutrino luminosities drop steeply but continue on a 1$-$2 orders of magnitude lower level for several 100 ms because of aspherical accretion of neutrino and shock-heated matter, before the ultimately spherical collapse of the outer progenitor shells suppresses the neutrino emission to negligible values. In all shock-reviving models BH accretion swallows the entire neutrino-heated matter and the explosion energies decrease from maxima around 1.5$\times$10$^{51}$erg to zero within a few seconds latest. Nevertheless, the shock or a sonic pulse moves outward and may trigger mass loss, which we estimate by long-time simulations with the PROMETHEUS code. We also provide gravitational-wave signals. △ Less

Submitted 15 March, 2022; v1 submitted 17 December, 2021; originally announced December 2021.

Comments: Minor revisions for referee's comments; accepted by MNRAS

Nebular phase properties of supernova Ibc from He-star explosions

Authors: L. Dessart, D. J. Hillier, T. Sukhbold, S. E. Woosley, H. -T. Janka

Abstract: Following our recent work on Type II supernovae (SNe), we present a set of 1D nonlocal thermodynamic equilibrium radiative transfer calculations for nebular-phase Type Ibc SNe starting from state-of-the-art explosion models with detailed nucleosynthesis. Our grid of progenitor models is derived from He stars that were subsequently evolved under the influence of wind mass loss. These He stars, whic… ▽ More Following our recent work on Type II supernovae (SNe), we present a set of 1D nonlocal thermodynamic equilibrium radiative transfer calculations for nebular-phase Type Ibc SNe starting from state-of-the-art explosion models with detailed nucleosynthesis. Our grid of progenitor models is derived from He stars that were subsequently evolved under the influence of wind mass loss. These He stars, which most likely form through binary mass exchange, synthesize less oxygen than their single-star counterparts with the same zero-age main sequence (ZAMS) mass. This reduction is greater in He-star models evolved with an enhanced mass loss rate. We obtain a wide range of spectral properties at 200d. In models from He stars with an initial mass >6Msun, the [OI] 6300, 6364 is of comparable or greater strength than [CaII] 7291,7323 -- the strength of [OI] 6300, 6364 increases with He-star initial mass. In contrast, models from lower mass He stars exhibit a weak [OI] 6300, 6364, strong [CaII] 7291,7323, but also strong NII lines and FeII emission below 5500A. The ejecta density, modulated by the ejecta mass, the explosion energy, and clumping, has a critical impact on the gas ionization, line cooling, and the spectral properties. FeII dominates the emission below 5500A and is stronger at earlier nebular epochs. It ebbs as the SN ages, while the fractional flux in [OI] 6300, 6364 and [CaII] 7291,7323 increases, with a similar rate, as the ejecta recombine. Although the results depend on the adopted wind mass loss rate and pre-SN mass, we find that He stars of 6-8Msun initially (ZAMS mass of 23-28Msun) match adequately the properties of standard SNe Ibc. Our results for less massive He stars are more perplexing, since the predicted spectra are not seen in nature. They may be missed by current surveys or associated with Type Ibn SNe in which interaction dominates over decay power. [Abridged] △ Less

Submitted 25 September, 2021; originally announced September 2021.

Comments: Accepted for publication in A&A

Journal ref: A&A 656, A61 (2021)

The explosion of 9$-$29$M_\odot$ stars as Type II supernovae : results from radiative-transfer modeling at one year after explosion

Authors: Luc Dessart, D. John Hillier, Tuguldur Sukhbold, Stan Woosley, H. -T. Janka

Abstract: We present a set of nonlocal thermodynamic equilibrium steady-state calculations of radiative transfer for one-year old type II supernovae (SNe) starting from state-of-the-art explosion models computed with detailed nucleosynthesis. This grid covers single-star progenitors with initial masses between 9 and 29$M_{\odot}$, all evolved with KEPLER at solar metallicity and ignoring rotation. The [OI]… ▽ More We present a set of nonlocal thermodynamic equilibrium steady-state calculations of radiative transfer for one-year old type II supernovae (SNe) starting from state-of-the-art explosion models computed with detailed nucleosynthesis. This grid covers single-star progenitors with initial masses between 9 and 29$M_{\odot}$, all evolved with KEPLER at solar metallicity and ignoring rotation. The [OI]$λλ$$6300,6364$ line flux generally grows with progenitor mass, and H$α$ exhibits an equally strong and opposite trend. The [CaII]$λλ$$7291,\,7323$ strength increases at low $^{56}$Ni mass, low explosion energy, or with clumping. This CaII doublet, which forms primarily in the explosively-produced Si/S zones, depends little on the progenitor mass, but may strengthen if Ca$^+$ dominates in the H-rich emitting zones or if Ca is abundant in the O-rich zones. Indeed, Si-O shell merging prior to core collapse may boost the CaII doublet at the expense of the OI doublet, and may thus mimic the metal line strengths of a lower mass progenitor. We find that the $^{56}$Ni bubble effect has a weak impact, probably because it is too weak to induce much of an ionization shift in the various emitting zones. Our simulations compare favorably to observed SNe II, including SN2008bk (e.g., 9$M_{\odot}$ model), SN2012aw (12$M_{\odot}$ model), SN1987A (15$M_{\odot}$ model), or SN2015bs (25$M_{\odot}$ model with no Si-O shell merging). SNe II with narrow lines and a low $^{56}$Ni mass are well matched by the weak explosion of 9$-$11$M_{\odot}$ progenitors. The nebular-phase spectra of standard SNe II can be explained with progenitors in the mass range 12$-$15$M_{\odot}$, with one notable exception for SN2015bs. In the intermediate mass range, these mass estimates may increase by a few $M_{\odot}$ with allowance for clumping of the O-rich material or CO molecular cooling. △ Less

Submitted 27 May, 2021; originally announced May 2021.

Comments: Accepted for publication in A&A

Journal ref: A&A 652, A64 (2021)

The Pair-Instability Mass Gap for Black Holes

Authors: S. E. Woosley, Alexander Heger

Abstract: Stellar evolution theory predicts a "gap" in the black hole birth function caused by the pair instability. Presupernova stars that have a core mass below some limiting value, Mlo, after all pulsational activity is finished, collapse to black holes, whereas more massive ones, up to some limiting value, Mhi, explode, promptly and completely, as pair-instability supernovae. Previous work has suggeste… ▽ More Stellar evolution theory predicts a "gap" in the black hole birth function caused by the pair instability. Presupernova stars that have a core mass below some limiting value, Mlo, after all pulsational activity is finished, collapse to black holes, whereas more massive ones, up to some limiting value, Mhi, explode, promptly and completely, as pair-instability supernovae. Previous work has suggested Mlo is approximately 50 solar masses and Mhi is approximately 130 solar masses. These calculations have been challenged by recent LIGO observations that show many black holes merging with individual masses, Mlo is least some 65 solar masses. Here we explore four factors affecting the theoretical estimates for the boundaries of this mass gap: nuclear reaction rates, evolution in detached binaries, rotation, and hyper-Eddington accretion after black hole birth. Current uncertainties in reaction rates by themselves allow Mlo to rise to 64 solar masses and Mhi as large as 161 solar masses. Rapid rotation could further increase Mlo to about 70 solar masses, depending on the treatment of magnetic torques. Evolution in detached binaries and super-Eddington accretion can, with great uncertainty, increase Mlo still further. Dimensionless Kerr parameters close to unity are allowed for the more massive black holes produced in close binaries, though they are generally smaller. △ Less

Submitted 14 March, 2021; originally announced March 2021.

Comments: 16 pages, 4 figure panels, 6 tables, submitted to ApJ

Model Light Curves for Type Ib and Ic Supernovae

Authors: Stan Woosley, Tuguldur Sukhbold, Daniel Kasen

Abstract: Using the Monte Carlo code, SEDONA, multiband photometry and spectra are calculated for supernovae derived from stripped helium stars with presupernova masses from 2.2 to 10.0 $M_\odot$. The models are representative of evolution in close binaries and have previously been exploded using a parametrized one-dimensional model for neutrino-transport. A subset, those with presupernova masses in the ran… ▽ More Using the Monte Carlo code, SEDONA, multiband photometry and spectra are calculated for supernovae derived from stripped helium stars with presupernova masses from 2.2 to 10.0 $M_\odot$. The models are representative of evolution in close binaries and have previously been exploded using a parametrized one-dimensional model for neutrino-transport. A subset, those with presupernova masses in the range 2.2 - 5.6 $M_\odot$, have many properties in common with observed Type Ib and Ic supernovae, including a median ejected mass near 2 $M_\odot$, explosion energies near $1 \times 10^{51}$ erg, typical $^{56}$Ni masses 0.07 - 0.09 $M_\odot$, peak times of about 20 days, and a narrow range for the $V$-$R$ color index 10 days post $V$-maximum near 0.3 mag. The median peak bolometric luminosity, near 10$^{42.3}$ erg s$^{-1}$, is fainter, however, than for several observational tabulations and the brightest explosion has a bolometric luminosity of only 10$^{42.50}$ erg s$^{-1}$. The brightest absolute $B$, $V$, and $R$ magnitudes at peak are $-17.2$, $-17.8$, and $-18.0$. These limits are fainter than some allegedly typical Type Ib and Ic supernovae and could reflect problems in our models or the observational analysis. Helium stars with lower and higher masses also produce interesting transients that may have been observed including fast, faint, blue transients and long, red, faint Type Ic supernovae. New models are specifically presented for SN 2007Y, SN 2007gr, SN 2009jf, LSQ13abf, SN 2008D, and SN 2010X. △ Less

Submitted 15 February, 2021; v1 submitted 15 September, 2020; originally announced September 2020.

Comments: submitted to ApJ

Production of Lithium in Primordial Supernovae

Authors: Alexander Heger, Stan Woosley

Abstract: The first generation of stars is quite unique. The absence of metals likely affects their formation, with current models suggesting a much more top-heavy initial mass fraction than what we observe today, and some of their other properties, such as rotation rates and binarity, are largely unknown or constrained by direct observations. But even non-rotation single stars of a given mass will evolve q… ▽ More The first generation of stars is quite unique. The absence of metals likely affects their formation, with current models suggesting a much more top-heavy initial mass fraction than what we observe today, and some of their other properties, such as rotation rates and binarity, are largely unknown or constrained by direct observations. But even non-rotation single stars of a given mass will evolve quite differently due to the absence of the metals: the stars will mostly remain much more compact until their death, with the hydrogen-rich later reaching down ten teems deeper in radius then in modern stars. When they explode as supernovae, the exposure to the supernova neutrino flux is much enhanced, allowing for copious production of lithium. This production will not be constant for all stars but largely vary across the mass range. Such production even more challenges the presence of the Spite Plateau. △ Less

Submitted 11 February, 2020; originally announced February 2020.

Comments: 5 pages, 8 figures, Proc. "Lithium in the Universe: to Be or not to Be"

The Birth Function for Black Holes and Neutron Stars in Close Binaries

Authors: S. Woosley, Tuguldur Sukhbold, H. -T. Janka

Abstract: The mass function for black holes and neutron stars at birth is explored for mass-losing helium stars. These should resemble, more closely than similar studies of single hydrogen-rich stars, the results of evolution in close binary systems. The effects of varying the mass-loss rate and metallicity are calculated using a simple semi-analytic approach to stellar evolution that is tuned to reproduce… ▽ More The mass function for black holes and neutron stars at birth is explored for mass-losing helium stars. These should resemble, more closely than similar studies of single hydrogen-rich stars, the results of evolution in close binary systems. The effects of varying the mass-loss rate and metallicity are calculated using a simple semi-analytic approach to stellar evolution that is tuned to reproduce detailed numerical calculations. Though the total fraction of black holes made in stellar collapse events varies considerably with metallicity, mass-loss rate, and mass cutoff, from 5$\%$ to 30$\%$, the shapes of their birth functions are very similar for all reasonable variations in these quantities. Median neutron star masses are in the range 1.32 - 1.37 $M_\odot$ regardless of metallicity. The median black hole mass for solar metallicity is typically 8 to 9 $M_\odot$ if only initial helium cores below 40 $M_\odot$ (ZAMS mass less than 80 $M_\odot$) are counted, and 9 - 13 $M_\odot$, in most cases, if helium cores with initial masses up to 150 $M_\odot$ (ZAMS mass less than 300 $M_\odot$) contribute. As long as the mass-loss rate as a function of mass exhibits no strong non-linearities, the black hole birth function from 15 to 35 $M_\odot$ has a slope that depends mostly on the initial mass function for main sequence stars. These findings imply the possibility of constraining the initial mass function and the properties of mass loss in close binaries using ongoing measurements of gravitational wave radiation. The expected rotation rates of the black holes are briefly discussed. △ Less

Submitted 5 April, 2020; v1 submitted 28 January, 2020; originally announced January 2020.

Comments: submitted to ApJ

The Explosion of Helium Stars Evolved With Mass Loss

Authors: Thomas Ertl, Stan E. Woosley, Tuguldur Sukhbold, H. -Thomas Janka

Abstract: Light curves, explosion energies, and remnant masses are calculated for a grid of supernovae resulting from massive helium stars that have been evolved including mass loss. These presupernova stars should approximate the results of binary evolution for stars in interacting systems that lose their envelopes close to the time of helium core ignition. Initial helium star masses are in the range 2.5 t… ▽ More Light curves, explosion energies, and remnant masses are calculated for a grid of supernovae resulting from massive helium stars that have been evolved including mass loss. These presupernova stars should approximate the results of binary evolution for stars in interacting systems that lose their envelopes close to the time of helium core ignition. Initial helium star masses are in the range 2.5 to 40\,\Msun, which correspond to main sequence masses of about 13 to 90\,\Msun. Common Type Ib and Ic supernovae result from stars whose final masses are approximately 2.5 to 5.6\,\Msun. For heavier stars, a large fraction of collapses lead to black holes, though there is an island of explodability for presupernova masses near 10\,\Msun. The median neutron star mass in binaries is 1.35--1.38\,\Msun \ and the median black hole mass is between 9 and 11\,\Msun. Even though black holes less massive than 5 \Msun\ are rare, they are predicted down to the maximum neutron star mass. There is no empty ``gap'', only a less populated mass range. For standard assumptions regarding the explosions and nucleosynthesis, the models predict light curves that are fainter than the brighter common Type Ib and Ic supernovae. Even with a very liberal, but physically plausible increase in $^{56}$Ni production, the highest energy models are fainter, at peak, than 10$^{42.6}$\,erg\,s$^{-1}$, and very few approach that limit. The median peak luminosity ranges from 10$^{42.0}$ to 10$^{42.3}$\,erg\,s$^{-1}$. Possible alternatives to the standard neutrino-powered and radioactive-illuminated models are explored. Magnetars are a promising alternative. Several other unusual varieties of Type I supernovae at both high and low mass are explored. △ Less

Submitted 17 December, 2019; v1 submitted 3 October, 2019; originally announced October 2019.

Comments: 50 pages, 36 figures, 10 tables; submitted to ApJ

Gas Dynamics of the Nickel-56 Decay Heating in Pair-Instability Supernovae

Authors: Ke-Jung Chen, S. E. Woosley, Daniel J. Whalen

Abstract: Very massive 140-260 Msun stars can die as highly-energetic pair-instability supernovae (PI SNe) with energies of up to 100 times those of core-collapse SNe that can completely destroy the star, leaving no compact remnant behind. These explosions can synthesize $0.1-30$ Msun of radioactive Ni56, which can cause them to rebrighten at later times when photons due to Ni56 decay diffuse out of the eje… ▽ More Very massive 140-260 Msun stars can die as highly-energetic pair-instability supernovae (PI SNe) with energies of up to 100 times those of core-collapse SNe that can completely destroy the star, leaving no compact remnant behind. These explosions can synthesize $0.1-30$ Msun of radioactive Ni56, which can cause them to rebrighten at later times when photons due to Ni56 decay diffuse out of the ejecta. However, heat from the decay of such large masses of Ni56 could also drive important dynamical effects deep in the ejecta that are capable of mixing elements and affecting the observational signatures of these events. We have now investigated the dynamical effect of Ni56 heating on PI SN ejecta with high-resolution two-dimensional hydrodynamic simulations performed with the CASTRO code. We find that expansion of the hot Ni56 bubble forms a shell at the base of the silicon layer of the ejecta about 200 days after the explosion but that no hydrodynamical instabilities develop that would mix Ni56 with the Si/O-rich ejecta. However, while the dynamical effects of Ni56 heating may be weak they could affect the observational signatures of some PI SNe by diverting decay energy into internal expansion of the ejecta at the expense of rebrightening at later times. △ Less

Submitted 16 June, 2020; v1 submitted 29 April, 2019; originally announced April 2019.

Comments: Accepted to ApJ, 14 pages

Three-Dimensional Simulations of Magnetar-Powered Superluminous Supernovae

Authors: Ke-Jung Chen, S. E. Woosley, Daniel J. Whalen

Abstract: A rapidly spinning magnetar in a young supernova (SN) can produce a superluminous transient by converting a fraction of its rotational energy into radiation. Here, we present the first three-dimensional hydrodynamical simulations ever performed of a magnetar-powered SN in the circumstellar medium formed by the ejection of the outer layers of the star prior to the blast. We find that hydrodynamical… ▽ More A rapidly spinning magnetar in a young supernova (SN) can produce a superluminous transient by converting a fraction of its rotational energy into radiation. Here, we present the first three-dimensional hydrodynamical simulations ever performed of a magnetar-powered SN in the circumstellar medium formed by the ejection of the outer layers of the star prior to the blast. We find that hydrodynamical instabilities form on two scales in the ejecta, not just one as in ordinary core-collapse SNe: in the hot bubble energized by the magnetar and in the forward shock of the SN as it plows up ambient gas. Pressure from the bubble also makes the instabilities behind the forward shock more violent and causes more mixing in the explosion than in normal SNe, with important consequences for the light curves and spectra of the event that cannot be captured by one-dimensional models. We also find that the magnetar can accelerate Ca and Si to velocities of $\sim $ 12000 km/s and account for their broadened emission lines in observations. Our simulations also reveal that energy from even weak magnetars can accelerate iron-group elements deep in the ejecta to $5000-7000$ km/s and explain the high-velocity Fe observed at early times in some core-collapse SNe such as SN 1987A. △ Less

Submitted 20 April, 2020; v1 submitted 29 April, 2019; originally announced April 2019.

Comments: Accepted to ApJ, 12 pages

Multidimensional Radiation Hydrodynamics Simulations of Pulsational Pair-Instability Supernovae

Authors: Ke-Jung Chen, Daniel J. Whalen, S. E. Woosley, Weiqun Zhang

Abstract: Stars with masses of 80 - 130 Msun can encounter the pulsational pair-instability at the end of their lives, which triggers consecutive episodes of explosive burning that eject multiple massive shells. Collisions between these shells produce bright transients known as pulsational pair-instability supernovae (PPI SNe) that may explain some extreme supernovae. In this paper, we present the first 2D… ▽ More Stars with masses of 80 - 130 Msun can encounter the pulsational pair-instability at the end of their lives, which triggers consecutive episodes of explosive burning that eject multiple massive shells. Collisions between these shells produce bright transients known as pulsational pair-instability supernovae (PPI SNe) that may explain some extreme supernovae. In this paper, we present the first 2D and 3D radiation hydrodynamics simulations of PPI SNe with the CASTRO code. Radiative cooling causes the collided shells to evolve into thin, dense structures with hot spots that can enhance the peak luminosity of the SN by factors of 2 - 3. The light curve peaks at $1.9 - 2.1 \times 10^{43}$ erg s$^{-1}$ for 50 days and then plateaus at $2 - 3 \times 10^{42}$ erg s$^{-1}$ for 200 days, depending on viewing angle. The presence of C and O and absence of Si and Fe in its spectra can uniquely identify this transient as a PPI SN in follow-up observations. Our models suggest that multidimensional radiation hydrodynamics is required to model the evolution and light curves of all shell-collision SNe such as Type IIne, not just PPI SNe. △ Less

Submitted 21 July, 2023; v1 submitted 29 April, 2019; originally announced April 2019.

Comments: 17 pages, 12 figures, Accepted for publication in ApJ

Gravity and Light: Combining Gravitational Wave and Electromagnetic Observations in the 2020s

Authors: R. J. Foley, K. D. Alexander, I. Andreoni, I. Arcavi, K. Auchettl, J. Barnes, G. Baym, E. C. Bellm, A. M. Beloborodov, N. Blagorodnova, J. P. Blakeslee, P. R. Brady, M. Branchesi, J. S. Brown, N. Butler, M. Cantiello, R. Chornock, D. O. Cook, J. Cooke, D. L. Coppejans, A. Corsi, S. M. Couch, M. W. Coughlin, D. A. Coulter, P. S. Cowperthwaite , et al. (88 additional authors not shown)

Abstract: As of today, we have directly detected exactly one source in both gravitational waves (GWs) and electromagnetic (EM) radiation, the binary neutron star merger GW170817, its associated gamma-ray burst GRB170817A, and the subsequent kilonova SSS17a/AT 2017gfo. Within ten years, we will detect hundreds of events, including new classes of events such as neutron-star-black-hole mergers, core-collapse s… ▽ More As of today, we have directly detected exactly one source in both gravitational waves (GWs) and electromagnetic (EM) radiation, the binary neutron star merger GW170817, its associated gamma-ray burst GRB170817A, and the subsequent kilonova SSS17a/AT 2017gfo. Within ten years, we will detect hundreds of events, including new classes of events such as neutron-star-black-hole mergers, core-collapse supernovae, and almost certainly something completely unexpected. As we build this sample, we will explore exotic astrophysical topics ranging from nucleosynthesis, stellar evolution, general relativity, high-energy astrophysics, nuclear matter, to cosmology. The discovery potential is extraordinary, and investments in this area will yield major scientific breakthroughs. Here we outline some of the most exciting scientific questions that can be answered by combining GW and EM observations. △ Less

Submitted 11 March, 2019; originally announced March 2019.

Comments: submitted to Astro2020

Catching Element Formation In The Act

Authors: Chris L. Fryer, Frank Timmes, Aimee L. Hungerford, Aaron Couture, Fred Adams, Wako Aoki, Almudena Arcones, David Arnett, Katie Auchettl, Melina Avila, Carles Badenes, Eddie Baron, Andreas Bauswein, John Beacom, Jeff Blackmon, Stephane Blondin, Peter Bloser, Steve Boggs, Alan Boss, Terri Brandt, Eduardo Bravo, Ed Brown, Peter Brown, Steve Bruenn. Carl Budtz-Jorgensen, Eric Burns , et al. (194 additional authors not shown)

Abstract: Gamma-ray astronomy explores the most energetic photons in nature to address some of the most pressing puzzles in contemporary astrophysics. It encompasses a wide range of objects and phenomena: stars, supernovae, novae, neutron stars, stellar-mass black holes, nucleosynthesis, the interstellar medium, cosmic rays and relativistic-particle acceleration, and the evolution of galaxies. MeV gamma-ray… ▽ More Gamma-ray astronomy explores the most energetic photons in nature to address some of the most pressing puzzles in contemporary astrophysics. It encompasses a wide range of objects and phenomena: stars, supernovae, novae, neutron stars, stellar-mass black holes, nucleosynthesis, the interstellar medium, cosmic rays and relativistic-particle acceleration, and the evolution of galaxies. MeV gamma-rays provide a unique probe of nuclear processes in astronomy, directly measuring radioactive decay, nuclear de-excitation, and positron annihilation. The substantial information carried by gamma-ray photons allows us to see deeper into these objects, the bulk of the power is often emitted at gamma-ray energies, and radioactivity provides a natural physical clock that adds unique information. New science will be driven by time-domain population studies at gamma-ray energies. This science is enabled by next-generation gamma-ray instruments with one to two orders of magnitude better sensitivity, larger sky coverage, and faster cadence than all previous gamma-ray instruments. This transformative capability permits: (a) the accurate identification of the gamma-ray emitting objects and correlations with observations taken at other wavelengths and with other messengers; (b) construction of new gamma-ray maps of the Milky Way and other nearby galaxies where extended regions are distinguished from point sources; and (c) considerable serendipitous science of scarce events -- nearby neutron star mergers, for example. Advances in technology push the performance of new gamma-ray instruments to address a wide set of astrophysical questions. △ Less

Submitted 7 February, 2019; originally announced February 2019.

Comments: 14 pages including 3 figures

Report number: LA-UR-18-29748

The nature of PISN candidates: clues from nebular spectra

Authors: P. A. Mazzali, T. J. Moriya, M. Tanaka, S. E. Woosley

Abstract: A group of super-luminous supernovae (SL-SNe) characterised by broad light curves have been suggested to be Pair Instability SNe (PISNe). Nebular spectra computed using PISN models have failed to reproduce the broad emission lines observed in these SNe, casting doubts on their true nature. Here, models of both PISNe and the explosion following the collapse of the core of a very massive star (100 M… ▽ More A group of super-luminous supernovae (SL-SNe) characterised by broad light curves have been suggested to be Pair Instability SNe (PISNe). Nebular spectra computed using PISN models have failed to reproduce the broad emission lines observed in these SNe, casting doubts on their true nature. Here, models of both PISNe and the explosion following the collapse of the core of a very massive star (100 Msun) are used to compute nebular spectra, which are compared to the spectrum of the prototypical PISN candidate, SN 2007bi. PISN models are confirmed to produce synthetic spectra showing narrow emission lines, resulting from the confinement of 56Ni to the lowest velocities (<~ 2000 km/s) and in clear disagreement with the spectrum of SN 2007bi. Spectra more closely resembling SN 2007bi are obtained if the PISN models are fully mixed in abundance. Massive core-collapse models produce enough 56Ni to power the light curve of PISN candidates, but their spectra are also not adequate. The nebular spectrum of SN 2007bi can be successfully reproduced if the inner region is artificially filled with oxygen-rich, low-velocity ejecta. This most likely requires a grossly aspherical explosion. A major difference between PISN and massive collapse models is that the former emit much more strongly in the NIR. It is concluded that: a) current PISN candidates, in particular SN 2007bi, are more likely the result of the collapse and explosion of massive stars below the PI limit; b) significant asymmetry is required to reproduce the late-time spectrum of SN 2007bi. △ Less

Submitted 22 January, 2019; originally announced January 2019.

Comments: 13 pages, 9 figures, 1 table, accepted by MNRAS

The Evolution of Massive Helium Stars Including Mass Loss

Authors: S. E. Woosley

Abstract: The evolution of helium stars with initial masses in the range 1.6 to 120 Msun is studied, including the effects of mass loss by winds. These stars are assumed to form in binary systems when their expanding hydrogenic envelopes are promptly lost just after helium ignition. Significant differences are found with single star evolution, chiefly because the helium core loses mass during helium burning… ▽ More The evolution of helium stars with initial masses in the range 1.6 to 120 Msun is studied, including the effects of mass loss by winds. These stars are assumed to form in binary systems when their expanding hydrogenic envelopes are promptly lost just after helium ignition. Significant differences are found with single star evolution, chiefly because the helium core loses mass during helium burning rather than gaining it from hydrogen shell burning. Consequently presupernova stars for a given initial mass function have considerably smaller mass when they die and will be easier to explode. Even accounting for this difference, the helium stars with mass loss develop more centrally condensed cores that should explode more easily than their single-star counterparts. The production of low mass black holes may be diminished. Helium stars with initial masses below 3.2 Msun experience significant radius expansion after helium depletion, reaching blue supergiant proportions. This could trigger additional mass exchange or affect the light curve of the supernova. The most common black hole masses produced in binaries is estimated to be about 9 Msun. A new maximum mass for black holes derived from pulsational pair-instability supernovae is derived - 46 Msun, and a new potential gap at 10 - 12 Msun is noted. Models pertinent to SN 2014ft are presented and a library of presupernova models is generated. △ Less

Submitted 22 April, 2019; v1 submitted 1 January, 2019; originally announced January 2019.

Comments: 23 pages, 16 figures, 8 tables, submitted Astrophysical Journal Dec. 23, 2018; accepted April 21, 2019

Three-dimensional mixing and light curves: constraints on the progenitor of supernova 1987A

Authors: Victor Utrobin, Annop Wongwathanarat, H. -Thomas Janka, Ewald Mueller, T. Ertl, Stan Woosley

Abstract: With the same method as used previously, we investigate neutrino-driven explosions of a larger sample of blue supergiant models. The larger sample includes three new presupernova stars. The results are compared with light-curve observations of the peculiar type IIP SN 1987A. The explosions were modeled in 3D with the neutrino-hydrodynamics code PROMETHEUS-HOTB, and light-curve calculations were pe… ▽ More With the same method as used previously, we investigate neutrino-driven explosions of a larger sample of blue supergiant models. The larger sample includes three new presupernova stars. The results are compared with light-curve observations of the peculiar type IIP SN 1987A. The explosions were modeled in 3D with the neutrino-hydrodynamics code PROMETHEUS-HOTB, and light-curve calculations were performed in spherical symmetry with the radiation-hydrodynamics code CRAB. Our results confirm the basic findings of the previous work: 3D neutrino-driven explosions with SN 1987A-like energies synthesize an amount of Ni-56 that is consistent with the radioactive tail of the light curve. Moreover, the models mix hydrogen inward to minimum velocities below 400 km/s as required by spectral observations. Hydrodynamic simulations with the new progenitor models, which possess smaller radii than the older ones, show much better agreement between calculated and observed light curves in the initial luminosity peak and during the first 20 days. A set of explosions with similar energies demonstrated that a high growth factor of Rayleigh-Taylor instabilities at the (C+O)/He composition interface combined with a weak interaction of fast Rayleigh-Taylor plumes, where the reverse shock occurs below the He/H interface, provides a sufficient condition for efficient outward mixing of Ni-56 into the hydrogen envelope. This condition is realized to the required extent only in one of the older stellar models, which yielded a maximum velocity of around 3000 km/s for the bulk of ejected Ni-56, but failed to reproduce the helium-core mass of 6 Msun inferred from the absolute luminosity of the presupernova star. We conclude that none of the single-star progenitor models proposed for SN 1987A to date satisfies all constraints set by observations. (Abridged) △ Less

Submitted 19 March, 2019; v1 submitted 24 December, 2018; originally announced December 2018.

Comments: 16 pages, 9 figures, 3 tables; accepted for publication in Astron. Astrophys

Journal ref: A&A 624, A116 (2019)

Observatory science with eXTP

Authors: Jean J. M. in 't Zand, Enrico Bozzo, Jinlu Qu, Xiang-Dong Li, Lorenzo Amati, Yang Chen, Immacolata Donnarumma, Victor Doroshenko, Stephen A. Drake, Margarita Hernanz, Peter A. Jenke, Thomas J. Maccarone, Simin Mahmoodifar, Domitilla de Martino, Alessandra De Rosa, Elena M. Rossi, Antonia Rowlinson, Gloria Sala, Giulia Stratta, Thomas M. Tauris, Joern Wilms, Xuefeng Wu, Ping Zhou, Iván Agudo, Diego Altamirano , et al. (159 additional authors not shown)

Abstract: In this White Paper we present the potential of the enhanced X-ray Timing and Polarimetry (eXTP) mission for studies related to Observatory Science targets. These include flaring stars, supernova remnants, accreting white dwarfs, low and high mass X-ray binaries, radio quiet and radio loud active galactic nuclei, tidal disruption events, and gamma-ray bursts. eXTP will be excellently suited to stu… ▽ More In this White Paper we present the potential of the enhanced X-ray Timing and Polarimetry (eXTP) mission for studies related to Observatory Science targets. These include flaring stars, supernova remnants, accreting white dwarfs, low and high mass X-ray binaries, radio quiet and radio loud active galactic nuclei, tidal disruption events, and gamma-ray bursts. eXTP will be excellently suited to study one common aspect of these objects: their often transient nature. Developed by an international Consortium led by the Institute of High Energy Physics of the Chinese Academy of Science, the eXTP mission is expected to be launched in the mid 2020s. △ Less

Submitted 10 December, 2018; originally announced December 2018.

Comments: Accepted for publication on Sci. China Phys. Mech. Astron. (2019)

A UV Resonance Line Echo from a Shell Around a Hydrogen-Poor Superluminous Supernova

Authors: R. Lunnan, C. Fransson, P. M. Vreeswijk, S. E. Woosley, G. Leloudas, D. A. Perley, R. M. Quimby, Lin Yan, N. Blagorodnova, B. D. Bue, S. B. Cenko, A. De Cia, D. O. Cook, C. U. Fremling, P. Gatkine, A. Gal-Yam, M. M. Kasliwal, S. R. Kulkarni, F. J. Masci, P. E. Nugent, A. Nyholm, A. Rubin, N. Suzuki, P. Wozniak

Abstract: Hydrogen-poor superluminous supernovae (SLSN-I) are a class of rare and energetic explosions discovered in untargeted transient surveys in the past decade. The progenitor stars and the physical mechanism behind their large radiated energies ($\sim10^{51}$ erg) are both debated, with one class of models primarily requiring a large rotational energy, while the other requires very massive progenitors… ▽ More Hydrogen-poor superluminous supernovae (SLSN-I) are a class of rare and energetic explosions discovered in untargeted transient surveys in the past decade. The progenitor stars and the physical mechanism behind their large radiated energies ($\sim10^{51}$ erg) are both debated, with one class of models primarily requiring a large rotational energy, while the other requires very massive progenitors to either convert kinetic energy into radiation via interaction with circumstellar material (CSM), or engender a pair-instability explosion. Observing the structure of the CSM around SLSN-I offers a powerful test of some scenarios, though direct observations are scarce. Here, we present a series of spectroscopic observations of the SLSN-I iPTF16eh, which reveal both absorption and time- and frequency-variable emission in the Mg II resonance doublet. We show that these observations are naturally explained as a resonance scattering light echo from a circumstellar shell. Modeling the evolution of the emission, we find a shell radius of 0.1 pc and velocity of 3300 km s$^{-1}$, implying the shell was ejected three decades prior to the supernova explosion. These properties match theoretical predictions of pulsational pair-instability shell ejections, and imply the progenitor had a He core mass of $\sim 50-55~{\rm M}_{\odot}$, corresponding to an initial mass of $\sim 115~{\rm M}_{\odot}$. △ Less

Submitted 18 August, 2018; v1 submitted 14 August, 2018; originally announced August 2018.

Comments: Accepted. Fixed typo in table header, otherwise unchanged from previous version

The 30-Year Search for the Compact Object in SN 1987A

Authors: Dennis Alp, Josefin Larsson, Claes Fransson, Remy Indebetouw, Anders Jerkstrand, Antero Ahola, David Burrows, Peter Challis, Phil Cigan, Aleksandar Cikota, Robert P. Kirshner, Jacco Th. van Loon, Seppo Mattila, C. -Y. Ng, Sangwook Park, Jason Spyromilio, S. E. Woosley, Maarten Baes, Patrice Bouchet, Roger A. Chevalier, Kari A. Frank, Bryan M. Gaensler, Haley L. Gomez, H. -Thomas Janka, Bruno Leibundgut , et al. (10 additional authors not shown)

Abstract: Despite more than 30 years of searches, the compact object in Supernova (SN) 1987A has not yet been detected. We present new limits on the compact object in SN 1987A using millimeter, near-infrared, optical, ultraviolet, and X-ray observations from ALMA, VLT, HST, and Chandra. The limits are approximately 0.1 mJy ($0.1\times 10^{-26}$ erg s$^{-1}$ cm$^{-2}$ Hz$^{-1}$) at 213 GHz, 1 Lsun (… ▽ More Despite more than 30 years of searches, the compact object in Supernova (SN) 1987A has not yet been detected. We present new limits on the compact object in SN 1987A using millimeter, near-infrared, optical, ultraviolet, and X-ray observations from ALMA, VLT, HST, and Chandra. The limits are approximately 0.1 mJy ($0.1\times 10^{-26}$ erg s$^{-1}$ cm$^{-2}$ Hz$^{-1}$) at 213 GHz, 1 Lsun ($6\times 10^{-29}$ erg s$^{-1}$ cm$^{-2}$ Hz$^{-1}$) in optical if our line-of-sight is free of ejecta dust, and $10^{36}$ erg s$^{-1}$ ($2\times 10^{-30}$ erg s$^{-1}$ cm$^{-2}$ Hz$^{-1}$) in 2-10 keV X-rays. Our X-ray limits are an order of magnitude less constraining than previous limits because we use a more realistic ejecta absorption model based on three-dimensional neutrino-driven SN explosion models (presented in an accompanying article). The allowed bolometric luminosity of the compact object is 22 Lsun if our line-of-sight is free of ejecta dust, or 138 Lsun if dust-obscured. Depending on assumptions, these values limit the effective temperature of a neutron star to <4-8 MK and do not exclude models, which typically are in the range 3-4 MK. For the simplest accretion model, the accretion rate for an efficiency $η$ is limited to $< 10^{-11} η^{-1}$ Msun yr$^{-1}$, which excludes most predictions. For pulsar activity modeled by a rotating magnetic dipole in vacuum, the limit on the magnetic field strength ($B$) for a given spin period ($P$) is $B < 10^{14} P^2$ G s$^{-2}$. By combining information about radiation reprocessing and geometry, it is likely that the compact object is a dust-obscured thermally-emitting neutron star, which may appear as a region of higher-temperature ejecta dust emission. △ Less

Submitted 30 July, 2018; v1 submitted 11 May, 2018; originally announced May 2018.

Comments: 35 pages, 10 figures. Accepted for publication in ApJ. Updated to match accepted version; slightly shortened and clarified the discussion

Models for the Unusual Supernova iPTF14hls

Authors: S. E. Woosley

Abstract: Supernova iPTF14hls maintained a bright, variable luminosity for more than 600 days, while lines of hydrogen and iron in its spectrum had different speeds, but showed little evolution. Here several varieties of models are explored for iPTF14hls-like events. They are based upon circumstellar medium (CSM) interaction in an ordinary supernova, pulsational pair-instability supernovae (PPISN), and magn… ▽ More Supernova iPTF14hls maintained a bright, variable luminosity for more than 600 days, while lines of hydrogen and iron in its spectrum had different speeds, but showed little evolution. Here several varieties of models are explored for iPTF14hls-like events. They are based upon circumstellar medium (CSM) interaction in an ordinary supernova, pulsational pair-instability supernovae (PPISN), and magnetar formation. Each is able to explain the enduring emission and brightness of iPTF14hls, but has shortcomings when confronted with other observed characteristics. The PPISN model can, in some cases, produce a presupernova transient like the one observed at the site of iPTF14hls in 1954. It also offers a clear path to providing the necessary half solar mass of material at $\sim 5 \times 10^{16}$ cm for CSM interaction to work, and can give an irregular light curve without invoking additional assumptions. It explains the 4000 km s$^{-1}$ seen in the iron lines, but without additional energy input, strains to explain the nearly constant 8000 km s$^{-1}$ velocity seen in H$_α$. Magnetar models can also explain many of the observed features, but give a smooth light curve and may require an evolving magnetic field strength. Their dynamics may be difficult to reconcile with the observation of slow-moving hydrogen at late times. The various models predict different spectral characteristics and a remnant that, today, could be a black hole, magnetar, or even a star. Further observations and calculations of radiation transport will narrow the range of possibilities. △ Less

Submitted 25 January, 2018; originally announced January 2018.

Comments: submitted to Astrophysical Journal

Emission line models for the lowest-mass core collapse supernovae. I: Case study of a 9 $M_\odot$ one-dimensional neutrino-driven explosion

Authors: A. Jerkstrand, T. Ertl, H. -T. Janka, E. Müller, T. Sukhbold, S. E. Woosley

Abstract: 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 C… ▽ More 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. △ Less

Submitted 12 October, 2017; originally announced October 2017.

Comments: Resubmitted to MNRAS after referee comments

High Resolution Study of Presupernova Compactness

Authors: Tuguldur Sukhbold, Stan Woosley, Alexander Heger

Abstract: 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… ▽ More 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$ü$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. △ Less

Submitted 22 March, 2018; v1 submitted 9 October, 2017; originally announced October 2017.

Comments: 24 pages, 17 figures, 4 tables. Submitted to ApJ

Magnetar-Powered Supernovae in Two Dimensions. II. Broad-Line Supernovae Ic

Authors: Ke-Jung Chen, Takashi J. Moriya, Stan Woosley, Tuguldur Sukhbold, Daniel J. Whalen, Yudai Suwa, Volker Bromm

Abstract: Nascent neutron stars with millisecond periods and magnetic fields in excess of $10^{16}$ Gauss can drive highly energetic and asymmetric explosions known as magnetar-powered supernovae. These exotic explosions are one theoretical interpretation for supernovae Ic-BL which are sometimes associated with long gamma-ray bursts. Twisted magnetic field lines extract the rotational energy of the neutron… ▽ More Nascent neutron stars with millisecond periods and magnetic fields in excess of $10^{16}$ Gauss can drive highly energetic and asymmetric explosions known as magnetar-powered supernovae. These exotic explosions are one theoretical interpretation for supernovae Ic-BL which are sometimes associated with long gamma-ray bursts. Twisted magnetic field lines extract the rotational energy of the neutron star and release it as a disk wind or a jet with energies greater than 10$^{52}$ erg over $\sim 20$ sec. What fractions of the energy of the central engine go into the wind and the jet remain unclear. We have performed two-dimensional hydrodynamical simulations of magnetar-powered supernovae (SNe) driven by disk winds and jets with the CASTRO code to investigate the effect of the central engine on nucleosynthetic yields, mixing, and light curves. We find that these explosions synthesize less than 0.05 Msun of Ni and that this mass is not very sensitive to central engine type. The morphology of the explosion can provide a powerful diagnostic of the properties of the central engine. In the absence of a circumstellar medium these events are not very luminous, with peak bolometric magnitudes $M_b \sim -16.5 $ due to low Ni production. △ Less

Submitted 21 June, 2017; originally announced June 2017.

Comments: Accepted to ApJ

Very Deep Inside the SN 1987A Core Ejecta: Molecular Structures Seen in 3D

Authors: F. J. Abellán, R. Indebetouw, J. M. Marcaide, M. Gabler, C. Fransson, J. Spyromilio, D. N. Burrows, R. Chevalier, P. Cigan, B. M. Gaensler, H. L. Gomez, H. -Th. Janka, R. Kirshner, J. Larsson, P. Lundqvist, M. Matsuura, R. McCray, C. -Y. Ng, S. Park, P. Roche, L. Staveley-Smith, J. Th. Van Loon, J. C. Wheeler, S. E. Woosley

Abstract: Most massive stars end their lives in core-collapse supernova explosions and enrich the interstellar medium with explosively nucleosynthesized elements. Following core collapse, the explosion is subject to instabilities as the shock propagates outwards through the progenitor star. Observations of the composition and structure of the innermost regions of a core-collapse supernova provide a direct p… ▽ More Most massive stars end their lives in core-collapse supernova explosions and enrich the interstellar medium with explosively nucleosynthesized elements. Following core collapse, the explosion is subject to instabilities as the shock propagates outwards through the progenitor star. Observations of the composition and structure of the innermost regions of a core-collapse supernova provide a direct probe of the instabilities and nucleosynthetic products. SN 1987A in the Large Magellanic Cloud (LMC) is one of very few supernovae for which the inner ejecta can be spatially resolved but are not yet strongly affected by interaction with the surroundings. Our observations of SN 1987A with the Atacama Large Millimeter/submillimeter Array (ALMA) are of the highest resolution to date and reveal the detailed morphology of cold molecular gas in the innermost regions of the remnant. The 3D distributions of carbon and silicon monoxide (CO and SiO) emission differ, but both have a central deficit, or torus-like distribution, possibly a result of radioactive heating during the first weeks ("nickel heating"). The size scales of the clumpy distribution are compared quantitatively to models, demonstrating how progenitor and explosion physics can be constrained. △ Less

Submitted 14 June, 2017; originally announced June 2017.

Comments: ApJL accepted, 9 pages, 5 figures

Very Low Energy Supernovae: Light Curves and Spectra of Shock Breakout

Authors: Elizabeth Lovegrove, S. E. Woosley, Weiqun Zhang

Abstract: The brief transient emitted as a shock wave erupts through the surface of a presupernova star carries information about the stellar radius and explosion energy. Here the CASTRO code, which treats radiation transport using multigroup flux-limited diffusion, is used to simulate the light curves and spectra of shock breakout in very low-energy supernovae (VLE SNe), explosions in giant stars with fina… ▽ More The brief transient emitted as a shock wave erupts through the surface of a presupernova star carries information about the stellar radius and explosion energy. Here the CASTRO code, which treats radiation transport using multigroup flux-limited diffusion, is used to simulate the light curves and spectra of shock breakout in very low-energy supernovae (VLE SNe), explosions in giant stars with final kinetic energy much less than 10$^{51}$ erg. VLE SNe light curves, computed here with the KEPLER code, are distinctively faint, red, and long-lived, making them challenging to find with transient surveys. The accompanying shock breakouts are brighter, though briefer, and potentially easier to detect. Previous analytic work provides general guidance, but numerical simulations are challenging due to the range of conditions and lack of equilibration between color and effective temperatures. We consider previous analytic work and extend discussions of color temperature and opacity to the lower energy range explored by these events. Since this is the first application of the CASTRO code to shock breakout, test simulations of normal energy shock breakout of SN1987A are carried out and compared with the literature. A set of breakout light curves and spectra are then calculated for VLE SNe with final kinetic energies in the range $10^{47} - 10^{50}$ ergs for red supergiants with main sequence masses 15 Msun and 25 Msun. The importance of uncertainties in stellar atmosphere model, opacity, and ambient medium is discussed, as are observational prospects with current and forthcoming missions. △ Less

Submitted 7 June, 2017; originally announced June 2017.

Comments: 19 pages; submitted to Astrophysical Journal

ALMA spectral survey of Supernova 1987A --- molecular inventory, chemistry, dynamics and explosive nucleosynthesis

Authors: M. Matsuura, R. Indebetouw, S. Woosley, V. Bujarrabal, F. J. Abellan, R. McCray, J. Kamenetzky, C. Fransson, M. J. Barlow, H. L. Gomez, P. Cigan, I De Looze, J. Spyromilio, L. Staveley-Smith, G. Zanardo, P. Roche, J. Larsson, S. Viti, J. Th. van Loon, J. C. Wheeler, M. Baes, R. Chevalier, P. Lundqvist, J. M. Marcaide, E. Dwek , et al. (4 additional authors not shown)

Abstract: We report the first molecular line survey of Supernova 1987A in the millimetre wavelength range. In the ALMA 210--300 and 340--360 GHz spectra, we detected cold (20--170 K) CO, 28SiO, HCO+ and SO, with weaker lines of 29SiO from ejecta. This is the first identification of HCO+ and SO in a young supernova remnant. We find a dip in the J=6--5 and 5--4 SiO line profiles, suggesting that the ejecta mo… ▽ More We report the first molecular line survey of Supernova 1987A in the millimetre wavelength range. In the ALMA 210--300 and 340--360 GHz spectra, we detected cold (20--170 K) CO, 28SiO, HCO+ and SO, with weaker lines of 29SiO from ejecta. This is the first identification of HCO+ and SO in a young supernova remnant. We find a dip in the J=6--5 and 5--4 SiO line profiles, suggesting that the ejecta morphology is likely elongated. The difference of the CO and SiO line profiles is consistent with hydrodynamic simulations, which show that Rayleigh-Taylor instabilities cause mixing of gas, with heavier elements much more disturbed, making more elongated structure. We obtained isotopologue ratios of 28SiO/29SiO>13, 28SiO/30SiO>14, and 12CO/13CO>21, with the most likely limits of 28SiO/29SiO>128, 28SiO/30SiO>189. Low 29Si and 30Si abundances in SN 1987A are consistent with nucleosynthesis models that show inefficient formation of neutron-rich isotopes in a low metallicity environment, such as the Large Magellanic Cloud. The deduced large mass of HCO+ (~5x10^-6 Msun) and small SiS mass (<6x10^-5 Msun) might be explained by some mixing of elements immediately after the explosion. The mixing might have caused some hydrogen from the envelope to sink into carbon and oxygen-rich zones after the explosion, enabling the formation of a substantial mass of HCO+. Oxygen atoms may have penetrated into silicon and sulphur zones, suppressing formation of SiS. Our ALMA observations open up a new window to investigate chemistry, dynamics and explosive-nucleosynthesis in supernovae. △ Less

Submitted 7 April, 2017; originally announced April 2017.

Comments: Accepted by MNRAS, 16 pages

GRB 161219B-SN 2016jca: a powerful stellar collapse

Authors: C. Ashall, P. A. Mazzali, E. Pian, S. E. Woosley, E. Palazzi, S. J. Prentice, S. Kobayashi, S. Holmbo, A. Levan, D. Perley, M. D. Stritzinger, F. Bufano, A. V. Filippenko, A. Melandri, S. Oates, A. Rossi, J. Selsing, W. Zheng, A. J. Castro-Tirado, G. Chincarini, P. D'Avanzo, M. De Pasquale, S. Emery, A. S. Fruchter, K. Hurley , et al. (4 additional authors not shown)

Abstract: We report observations and analysis of the nearby gamma-ray burst GRB\,161219B (redshift $z=0.1475$) and the associated Type Ic supernova (SN) 2016jca. GRB\,161219B had an isotropic gamma-ray energy of $\sim 1.6 \times 10^{50}$\,erg. Its afterglow is likely refreshed at an epoch preceding the first photometric points (0.6\,d), which slows down the decay rates. Combined analysis of the SN light cur… ▽ More We report observations and analysis of the nearby gamma-ray burst GRB\,161219B (redshift $z=0.1475$) and the associated Type Ic supernova (SN) 2016jca. GRB\,161219B had an isotropic gamma-ray energy of $\sim 1.6 \times 10^{50}$\,erg. Its afterglow is likely refreshed at an epoch preceding the first photometric points (0.6\,d), which slows down the decay rates. Combined analysis of the SN light curve and multiwavelength observations of the afterglow suggest that the GRB jet was broad during the afterglow phase (full opening angle $\sim 42^\circ \pm 3^\circ$). Our spectral series shows broad absorption lines typical of GRB supernovae (SNe), which testify to the presence of material with velocities up to $\sim 0.25$c. The spectrum at 3.73\,d allows for the very early identification of a SN associated with a GRB. Reproducing it requires a large photospheric velocity ($35,000 \pm 7000$\,\kms). The kinetic energy of the SN is estimated through models to be \KE $\approx 4 \times 10^{52}$\,erg in spherical symmetry. The ejected mass in the explosion was \Mej $\approx 6.5 \pm 1.5$\,\Msun, much less than that of other GRB-SNe, demonstrating diversity among these events. The total amount of \Nifs\ in the explosion was $0.27 \pm 0.05$\,\Msun. The observed spectra require the presence of freshly synthesised \Nifs\ at the highest velocities, at least 3 times more than a standard GRB-SN. We also find evidence for a decreasing \Nifs\ abundance as a function of decreasing velocity. This suggests that SN\,2016jca was a highly aspherical explosion viewed close to on-axis, powered by a compact remnant. Applying a typical correction for asymmetry, the energy of SN\,2016jca was $\sim$ (1--3) $\times 10^{52}$\,erg, confirming that most of the energy produced by GRB-SNe goes into the kinetic energy of the SN ejecta. △ Less

Submitted 6 June, 2019; v1 submitted 14 February, 2017; originally announced February 2017.

Comments: Accepted for publication in MNRAS, minor changes to figure

Pulsational Pair-Instability Supernovae

Authors: S. E. Woosley

Abstract: The final evolution of stars in the mass range 70 - 140 solar masses is explored. Depending upon their mass loss history and rotation rates, these stars will end their lives as pulsational pair-instability supernovae producing a great variety of observational transients with total durations ranging from weeks to millennia and luminosities from 10$^{41}$ to over 10$^{44}$ erg s$^{-1}$. No non-rotat… ▽ More The final evolution of stars in the mass range 70 - 140 solar masses is explored. Depending upon their mass loss history and rotation rates, these stars will end their lives as pulsational pair-instability supernovae producing a great variety of observational transients with total durations ranging from weeks to millennia and luminosities from 10$^{41}$ to over 10$^{44}$ erg s$^{-1}$. No non-rotating model radiates more than $5 \times 10^{50}$ erg of light or has a kinetic energy exceeding $5 \times 10^{51}$ erg, but greater energies are possible, in principle, in magnetar-powered explosions which are explored. Many events resemble Type Ibn, Icn, and IIn supernovae, and some potential observational counterparts are mentioned. Some PPISN can exist in a dormant state for extended periods, producing explosions millennia after their first violent pulse. These dormant supernovae contain bright Wolf-Rayet stars, possibly embedded in bright x-ray and radio sources. The relevance of PPISN to supernova impostors like Eta Carinae, to super-luminous supernovae, and to sources of gravitational radiation is discussed. No black holes between 52 and 133 solar masses are expected from stellar evolution in close binaries. △ Less

Submitted 12 January, 2017; v1 submitted 31 August, 2016; originally announced August 2016.

Comments: 39 pages, 28 figures, submitted to Astrophysical Journal, Aug. 31, 2016. Revised and accepted Jan. 2, 2017

The Effect of Pair-Instability Mass Loss on Black Hole Mergers

Authors: K. Belczynski, A. Heger, W. Gladysz, A. J. Ruiter, S. Woosley, G. Wiktorowicz, H. -Y. Chen, T. Bulik, R. O'Shaughnesy, D. E. Holz, C. L. Fryer, E. Berti

Abstract: Mergers of two stellar origin black holes are a prime source of gravitational waves and are under intensive investigations. One crucial ingredient in their modeling has so far been neglected. Pair-instability pulsation supernovae with associated severe mass loss may suppress formation of massive black holes, decreasing black hole merger rates for the highest black hole masses. The mass loss asso… ▽ More Mergers of two stellar origin black holes are a prime source of gravitational waves and are under intensive investigations. One crucial ingredient in their modeling has so far been neglected. Pair-instability pulsation supernovae with associated severe mass loss may suppress formation of massive black holes, decreasing black hole merger rates for the highest black hole masses. The mass loss associated with pair-instability pulsation supernovae limits the Population I/II stellar-origin black hole mass to 50 Msun, in tension with earlier predictions that the maximum black hole mass could be as high as 100 Msun. Suppression of double black hole merger rates by pair-instability pulsation supernovae is negligible for our evolutionary channel. Our standard evolutionary model with inclusion of pair-instability pulsation supernovae and with pair-instability supernovae is fully consistent with the LIGO detections of black hole mergers: GW150914, GW151226 and LVT151012. The LIGO observations seem to exclude high (>400 km/s) BH natal kicks. We predict the detection of several up to 60 BH-BH mergers with total mass 10--150 Msun (most likely range: 20--80 Msun) in the forthcoming 60 effective days of LIGO O2 observations. △ Less

Submitted 11 July, 2016; originally announced July 2016.

Comments: A&A; comments welcome

Journal ref: A&A 594, A97 (2016)

Magnetar-Powered Supernovae in Two Dimensions. I. Superluminous Supernovae

Authors: Ke-Jung Chen, S. E. Woosley, Tuguldur Sukhbold

Abstract: Previous studies have shown that the radiation emitted by a rapidly rotating magnetar embedded in a young supernova can greatly amplify its luminosity. These one-dimensional studies have also revealed the existence of an instability arising from the piling up of radiatively accelerated matter in a thin dense shell deep inside the supernova. Here we examine the problem in two dimensions and find th… ▽ More Previous studies have shown that the radiation emitted by a rapidly rotating magnetar embedded in a young supernova can greatly amplify its luminosity. These one-dimensional studies have also revealed the existence of an instability arising from the piling up of radiatively accelerated matter in a thin dense shell deep inside the supernova. Here we examine the problem in two dimensions and find that, while instabilities cause mixing and fracture this shell into filamentary structures that reduce the density contrast, the concentration of matter in a hollow shell persists. The extent of the mixing depends upon the relative energy input by the magnetar and the kinetic energy of the inner ejecta. The light curve and spectrum of the resulting supernova will be appreciably altered, as will the appearance of the supernova remnant, which will be shellular and filamentary. A similar pile up and mixing might characterize other events where energy is input over an extended period by a centrally concentrated source, e.g. a pulsar, radioactive decay, a neutrino-powered wind, or colliding shells. The relevance of our models to the recent luminous transient ASASSN-15lh is briefly discussed. △ Less

Submitted 2 October, 2016; v1 submitted 27 April, 2016; originally announced April 2016.

Comments: Accepted to ApJ

The Progenitor of GW 150914

Authors: S. E. Woosley

Abstract: The spectacular detection of gravitational waves (GWs) from GW 150914 and its reported association with a gamma-ray burst (GRB) offer new insights into the evolution of massive stars. Here it is shown that no single star of any mass and credible metallicity is likely to produce the observed GW signal. Stars with helium cores in the mass range 35 to 133 solar masses encounter the pair instability a… ▽ More The spectacular detection of gravitational waves (GWs) from GW 150914 and its reported association with a gamma-ray burst (GRB) offer new insights into the evolution of massive stars. Here it is shown that no single star of any mass and credible metallicity is likely to produce the observed GW signal. Stars with helium cores in the mass range 35 to 133 solar masses encounter the pair instability and either explode or pulse until the core mass is less than 40 solar masses, smaller than the combined mass of the observed black holes. The rotation of more massive helium cores is either braked by interaction with a slowly rotating hydrogen envelope, if one is present, or by mass loss, if one is not. The very short interval between the GW signal and the observed onset of the putative GRB in GW 150914 is also too short to have come from a single star. A more probable model for making the gravitational radiation is the delayed merger of two black holes made by 70 and 90 solar mass stars in a binary system. The more massive component was a pulsational-pair instability supernova before making the first black hole. △ Less

Submitted 16 April, 2016; v1 submitted 1 March, 2016; originally announced March 2016.

Comments: Submitted to Astrophysical Journal Letters February 29, 2016

Inferring supernova IIb/Ib/Ic ejecta properties from light curves and spectra: Correlations from radiative-transfer models

Authors: Luc Dessart, D. John Hillier, Stan Woosley, Eli Livne, Roni Waldman, Sung-Chul Yoon, Norbert Langer

Abstract: 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… ▽ More 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Å at maximum in SNe IIb/Ib, although some scatter arises from the uncertain level of chemical mixing. The OI7772Å 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$μ$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. △ Less

Submitted 19 February, 2016; originally announced February 2016.

Comments: Accepted for publication in the MNRAS

The Most Luminous Supernovae

Authors: Tuguldur Sukhbold, Stan Woosley

Abstract: 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-millisecon… ▽ More 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. △ Less

Submitted 13 March, 2016; v1 submitted 15 February, 2016; originally announced February 2016.

Comments: 5 pages, 2 figures and 1 table. Submitted to The Astrophysical Journal Letters

Low-energy Population III supernovae and the origin of extremely metal-poor stars

Authors: Ke-Jung Chen, Alexander Heger, Daniel J. Whalen, Takashi J. Moriya, Volker Bromm, Stan Woosley

Abstract: Some ancient, dim, metal-poor stars may have formed in the ashes of the first supernovae (SNe). If their chemical abundances can be reconciled with the elemental yields of specific Population III (Pop III) explosions, they could reveal the properties of primordial stars. But multidimensional simulations of such explosions are required to predict their yields because dynamical instabilities can dre… ▽ More Some ancient, dim, metal-poor stars may have formed in the ashes of the first supernovae (SNe). If their chemical abundances can be reconciled with the elemental yields of specific Population III (Pop III) explosions, they could reveal the properties of primordial stars. But multidimensional simulations of such explosions are required to predict their yields because dynamical instabilities can dredge material up from deep in the ejecta that would otherwise be predicted to fall back on to the central remnant and be lost in one-dimensional (1D) models. We have performed two-dimensional (2D) numerical simulations of two low-energy Pop III SNe, a 12.4 Msun explosion and a 60 Msun explosion, and find that they produce elemental yields that are a good fit to those measured in the most iron-poor star discovered to date, SMSS J031300.36-670839.3 (J031300). Fallback on to the compact remnant in these weak explosions accounts for the lack of measurable iron in J031300 and its low iron-group abundances in general. Our 2D explosions produce higher abundances of heavy elements (atomic number Z > 20) than their 1D counterparts due to dredge-up by fluid instabilities. Since almost no Ni is ejected by these weak SNe, their low luminosities will prevent their detection in the near-infrared with the James Webb Space Telescope and future 30-m telescopes on the ground. The only evidence that they ever occurred will be in the fossil abundance record. △ Less

Submitted 21 June, 2017; v1 submitted 26 January, 2016; originally announced January 2016.

Comments: Accepted to MNRAS

Core-Collapse Supernovae from 9 to 120 Solar Masses Based on Neutrino-powered Explosions

Authors: Tuguldur Sukhbold, T. Ertl, S. E. Woosley, Justin M. Brown, H. -T. Janka

Abstract: Nucleosynthesis, light curves, explosion energies, and remnant masses are calculated for a grid of supernovae resulting from massive stars with solar metallicity and masses from 9.0 to 120 solar masses. The full evolution is followed using an adaptive reaction network of up to 2000 nuclei. A novel aspect of the survey is the use of a one-dimensional neutrino transport model for the explosion. This… ▽ More Nucleosynthesis, light curves, explosion energies, and remnant masses are calculated for a grid of supernovae resulting from massive stars with solar metallicity and masses from 9.0 to 120 solar masses. The full evolution is followed using an adaptive reaction network of up to 2000 nuclei. A novel aspect of the survey is the use of a one-dimensional neutrino transport model for the explosion. This explosion model has been calibrated to give the observed energy for SN 1987A, using several standard progenitors, and for the Crab supernova using a 9.6 solar mass progenitor. As a result of using a calibrated central engine, the final kinetic energy of the supernova is variable and sensitive to the structure of the presupernova star. Many progenitors with extended core structures do not explode, but become black holes, and the masses of exploding stars do not form a simply connected set. The resulting nucleosynthesis agrees reasonably well with the sun provided that a reasonable contribution from Type Ia supernovae is also allowed, but with a deficiency of light s-process isotopes. The resulting neutron star IMF has a mean gravitational mass near 1.4 solar masses. The average black hole mass is about 9 solar masses if only the helium core implodes, and 14 solar masses if the entire presupernova star collapses. Only ~10% of supernovae come from stars over 20 solar masses and some of these are Type Ib or Ic. Some useful systematics of Type IIp light curves are explored. △ Less

Submitted 13 February, 2016; v1 submitted 15 October, 2015; originally announced October 2015.

Comments: 49 pages, 39 figures and 9 tables. Accepted to The Astrophysical Journal 26/1/2016

Radiative-transfer models for supernovae IIb/Ib/Ic from binary-star progenitors

Authors: Luc Dessart, D. John Hillier, Stan Woosley, Eli Livne, Roni Waldman, Sung-Chul Yoon, Norbert Langer

Abstract: 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… ▽ More 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". △ Less

Submitted 28 July, 2015; originally announced July 2015.

Comments: accepted to MNRAS

The first direct measurement of 12C(12C,n)23Mg at stellar energies

Authors: B. Bucher, X. D. Tang, X. Fang, A. Heger, S. Almaraz-Calderon, A. Alongi, A. D. Ayangeakaa, M. Beard, A. Best, J. Browne, C. Cahillane, M. Couder, R. J. deBoer, A. Kontos, L. Lamm, Y. J. Li, A. Long, W. Lu, S. Lyons, M. Notani, D. Patel, N. Paul, M. Pignatari, A. Roberts, D. Robertson , et al. (6 additional authors not shown)

Abstract: Neutrons produced by the carbon fusion reaction 12C(12C,n)23Mg play an important role in stellar nucleosynthesis. However, past studies have shown large discrepancies between experimental data and theory, leading to an uncertain cross section extrapolation at astrophysical energies. We present the first direct measurement that extends deep into the astrophysical energy range along with a new and i… ▽ More Neutrons produced by the carbon fusion reaction 12C(12C,n)23Mg play an important role in stellar nucleosynthesis. However, past studies have shown large discrepancies between experimental data and theory, leading to an uncertain cross section extrapolation at astrophysical energies. We present the first direct measurement that extends deep into the astrophysical energy range along with a new and improved extrapolation technique based on experimental data from the mirror reaction 12C(12C,p)23Na. The new reaction rate has been determined with a well-defined uncertainty that exceeds the precision required by astrophysics models. Using our constrained rate, we find that 12C(12C,n)23Mg is crucial to the production of Na and Al in Pop-III Pair Instability Supernovae. It also plays a non-negligible role in the production of weak s-process elements as well as in the production of the important galactic gamma-ray emitter 60Fe. △ Less

Submitted 14 July, 2015; originally announced July 2015.

Comments: http://journals.aps.org/prl/abstract/10.1103/PhysRevLett.114.251102

Journal ref: Phys. Rev. Lett. 114, 251102 (2015)

The Remarkable Deaths of 9 - 11 Solar Mass Stars

Authors: S. E. Woosley, Alexander Heger

Abstract: The post-helium burning evolution of stars from 7 to 11 solar masses is complicated by the lingering effects of degeneracy and off-center ignition. Here stars in this mass range are studied using a standard set of stellar physics. Two important aspects of the study are the direct coupling of a reaction network of roughly 220 nuclei to the structure calculation at all stages and the use of a sub gr… ▽ More The post-helium burning evolution of stars from 7 to 11 solar masses is complicated by the lingering effects of degeneracy and off-center ignition. Here stars in this mass range are studied using a standard set of stellar physics. Two important aspects of the study are the direct coupling of a reaction network of roughly 220 nuclei to the structure calculation at all stages and the use of a sub grid model to describe the convective bounded flame that develops during neon and oxygen burning. Below 9.0 solar masses, degenerate oxygen-neon cores form that may become either white dwarfs or electron-capture supernovae. Above 10.3 solar masses the evolution proceeds "normally" to iron-core collapse, without composition inversions or degenerate flashes. Emphasis here is upon the stars in between which typically ignite oxygen burning off center. After oxygen burns in a convectively bounded flame, silicon burning ignites in a degenerate flash that commences closer to the stellar center and with increasing violence for stars of larger mass. In some cases the silicon flash is so violent that it could lead to the early ejection of the hydrogen envelope. This might have interesting observable consequences. For example, the death of a 10.0 solar mass star could produce two supernova-like displays, a faint low energy event due to the silicon flash, and an unusually bright supernova many months later as the low energy ejecta from core collapse collides with the previously ejected envelope. The potential relation to the Crab supernova is discussed. △ Less

Submitted 25 May, 2015; originally announced May 2015.

Comments: Submitted to Astrophysical Journal January 5, 2015; revised May 8, 2015

A two-parameter criterion for classifying the explodability of massive stars by the neutrino-driven mechanism

Authors: T. Ertl, H. -Th. Janka, S. E. Woosley, T. Sukhbold, M. Ugliano

Abstract: Thus far, judging the fate of a massive star (either a neutron star (NS) or a black hole) solely by its structure prior to core collapse has been ambiguous. Our work and previous attempts find a non-monotonic variation of successful and failed supernovae with zero-age main-sequence mass, for which no single structural parameter can serve as a good predictive measure. However, we identify two param… ▽ More Thus far, judging the fate of a massive star (either a neutron star (NS) or a black hole) solely by its structure prior to core collapse has been ambiguous. Our work and previous attempts find a non-monotonic variation of successful and failed supernovae with zero-age main-sequence mass, for which no single structural parameter can serve as a good predictive measure. However, we identify two parameters computed from the pre-collapse structure of the progenitor, which in combination allow for a clear separation of exploding and non-exploding cases with only few exceptions (~1-2.5%) in our set of 621 investigated stellar models. One parameter is M4, defining the normalized enclosed mass for a dimensionless entropy per nucleon of s=4, and the other is mu4 = d(m/M_sun)/d(r/1000 km) at s=4, being the normalized mass-derivative at this location. The two parameters mu4 and M4*mu4 can be directly linked to the mass-infall rate, Mdot, of the collapsing star and the electron-type neutrino luminosity of the accreting proto-NS, L_nue ~ M_ns*Mdot, which play a crucial role in the "critical luminosity" concept for the theoretical description of neutrino-driven explosions as runaway phenomenon of the stalled accretion shock. All models were evolved employing the approach of Ugliano et al. for simulating neutrino-driven explosions in spherical symmetry. The neutrino emission of the accretion layer is approximated by a gray transport solver, while the uncertain neutrino emission of the 1.1 M_sun proto-NS core is parametrized by an analytic model. The free parameters connected to the core-boundary prescription are calibrated to reproduce the observables of Supernova 1987A for five different progenitor models. △ Less

Submitted 8 January, 2016; v1 submitted 25 March, 2015; originally announced March 2015.

Comments: 23 pages, 12 figures; accepted by ApJ; revised version considerably enlarged (Fig. 7 and Sect.3.6 added)

The LOFT perspective on neutron star thermonuclear bursts

Authors: J. J. M. in 't Zand, D. Altamirano, D. R. Ballantyne, S. Bhattacharyya, E. F. Brown, Y. Cavecchi, D. Chakrabarty, J. Chenevez, A. Cumming, N. Degenaar, M. Falanga, D. K. Galloway, A. Heger, J. José, L. Keek, M. Linares, S. Mahmoodifar, C. M. Malone, M. Méndez, M. C. Miller, F. B. S. Paerels, J. Poutanen, A. Rózanska, H. Schatz, M. Serino , et al. (9 additional authors not shown)

Abstract: This is a White Paper in support of the mission concept of the Large Observatory for X-ray Timing (LOFT), proposed as a medium-sized ESA mission. We discuss the potential of LOFT for the study of thermonuclear X-ray bursts on accreting neutron stars. For a summary, we refer to the paper. This is a White Paper in support of the mission concept of the Large Observatory for X-ray Timing (LOFT), proposed as a medium-sized ESA mission. We discuss the potential of LOFT for the study of thermonuclear X-ray bursts on accreting neutron stars. For a summary, we refer to the paper. △ Less

Submitted 12 January, 2015; originally announced January 2015.

Comments: White Paper in Support of the Mission Concept of the Large Observatory for X-ray Timing

Pair-Instability Supernovae of Non-Zero Metallicity Stars

Authors: Ke-Jung Chen, Alexander Heger, Stan Woosley, Ann Almgren, Daniel Whalen

Abstract: Observational evidence suggests that some very massive stars in the local Universe may die as pair-instability supernovae. We present 2D simulations of the pair-instability supernova of a non-zero metallicity star. We find that very little mixing occurs in this explosion because metals in the stellar envelope drive strong winds that strip the hydrogen envelope from the star prior to death. Consequ… ▽ More Observational evidence suggests that some very massive stars in the local Universe may die as pair-instability supernovae. We present 2D simulations of the pair-instability supernova of a non-zero metallicity star. We find that very little mixing occurs in this explosion because metals in the stellar envelope drive strong winds that strip the hydrogen envelope from the star prior to death. Consequently, a reverse shock cannot form and trigger fluid instabilities during the supernova. Only weak mixing driven by nuclear burning occurs in the earliest stages of the supernova, and it is too weak to affect the observational signatures of the explosion. △ Less

Submitted 27 January, 2015; v1 submitted 28 July, 2014; originally announced July 2014.

Comments: 6 pages, 4 figures, ASP proceedings for the 9th International Conference on Numerical Modelinf of Space Plasma Flows (ASTRONUM-2014), Long Beach, CA