Source author record

Stan Woosley

Stan Woosley appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

18works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

18 published item(s)

preprint2022arXiv

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

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.

preprint2020arXiv

Production of Lithium in Primordial Supernovae

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.

preprint2016arXiv

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

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.

preprint2016arXiv

The Most Luminous Supernovae

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.

preprint2015arXiv

Cosmological Impact of Population III Binaries

We present the results of the stellar feedback from Pop III binaries by employing improved, more realistic Pop III evolutionary stellar models. To facilitate a meaningful comparison, we consider a fixed mass of 60 solar masses (Msun) incorporated in Pop III stars, either contained in a single star, or split up in binary stars of 30 Msun each or an asymmetric case of one 45 Msun and one 15 Msun star. Whereas the sizes of the resulting HII regions are comparable across all cases, the HeIII regions around binary stars are significantly smaller than that of the single star. Consequently, the He$^{+}$ 1640 angstrom recombination line is expected to become much weaker. Supernova feedback exhibits great variety due to the uncertainty in possible explosion pathways. If at least one of the component stars dies as a hypernova about ten times more energetic than conventional core-collapse supernovae, the gas inside the host minihalo is effectively blown out, chemically enriching the intergalactic medium (IGM) to an average metallicity of $10^{-4}-10^{-3}$ solar metallicity (Zsun), out to $\sim 2$ kpc. The single star, however, is more likely to collapse into a black hole, accompanied by at most very weak explosions. The effectiveness of early chemical enrichment would thus be significantly reduced, in difference from the lower mass binary stars, where at least one component is likely to contribute to heavy element production and dispersal. Important new feedback physics is also introduced if close binaries can form high-mass x-ray binaries, leading to the pre-heating and -ionization of the IGM beyond the extent of the stellar HII regions.

preprint2015arXiv

Pair-Instability Supernovae of Non-Zero Metallicity Stars

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.

preprint2015arXiv

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

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

preprint2014arXiv

General Relativistic Instability Supernova of a Supermassive Population III Star

The formation of supermassive Population III stars with masses $\gtrsim$ 10,000 Msun in primeval galaxies in strong UV backgrounds at $z \sim$ 15 may be the most viable pathway to the formation of supermassive black holes by $z \sim$ 7. Most of these stars are expected to live for short times and then directly collapse to black holes, with little or no mass loss over their lives. But we have now discovered that non-rotating primordial stars with masses close to 55,000 Msun can instead die as highly energetic thermonuclear supernovae powered by explosive helium burning, releasing up to 10$ ^{55}$ erg, or about 10,000 times the energy of a Type Ia supernova. The explosion is triggered by the general relativistic contribution of thermal photons to gravity in the core of the star, which causes the core to contract and explosively burn. The energy release completely unbinds the star, leaving no compact remnant, and about half of the mass of the star is ejected into the early cosmos in the form of heavy elements. The explosion would be visible in the near infrared at $z \lesssim$ 20 to {\it Euclid} and the Wide-Field Infrared Survey Telescope (WFIRST), perhaps signaling the birth of supermassive black hole seeds and the first quasars.

preprint2014arXiv

Pair Instability Supernovae of Very Massive Population III Stars

Numerical studies of primordial star formation suggest that the first stars in the universe may have been very massive. Stellar models indicate that non-rotating Population III stars with initial masses of 140-260 Msun die as highly energetic pair-instability supernovae. We present new two-dimensional simulations of primordial pair-instability supernovae done with the CASTRO code. Our simulations begin at earlier times than previous multidimensional models, at the onset of core collapse, to capture any dynamical instabilities that may be seeded by collapse and explosive burning. Such instabilities could enhance explosive yields by mixing hot ash with fuel, thereby accelerating nuclear burning, and affect the spectra of the supernova by dredging up heavy elements from greater depths in the star at early times. Our grid of models includes both blue supergiants and red supergiants over the range in progenitor mass expected for these events. We find that fluid instabilities driven by oxygen and helium burning arise at the upper and lower boundaries of the oxygen shell $\sim$ 20 - 100 seconds after core bounce. Instabilities driven by burning freeze out after the SN shock exits the helium core. As the shock later propagates through the hydrogen envelope, a strong reverse shock forms that drives the growth of Rayleigh--Taylor instabilities. In red supergiant progenitors, the amplitudes of these instabilities are sufficient to mix the supernova ejecta.

preprint2014arXiv

Two-Dimensional Simulations of Pulsational Pair-Instability Supernovae

Massive stars that end their lives with helium cores in the range of 35 to 65 Msun are known to produce repeated thermonuclear outbursts due to a recurring pair-instability. In some of these events, solar masses of material are ejected in repeated outbursts of several times 10$^{50}$ erg each. Collisions between these shells can sometimes produce very luminous transients that are visible from the edge of the observable universe. Previous 1D studies of these events produce thin, high-density shells as one ejection plows into another. Here, in the first multidimensional simulations of these collisions, we show that the development of a Rayleigh-Taylor instability truncates the growth of the high density spike and drives mixing between the shells. The progenitor is a 110 Msun solar-metallicity star that was shown in earlier work to produce a superluminous supernova. The light curve of this more realistic model has a peak luminosity and duration that are similar to those of 1D models but a structure that is smoother.

preprint2013arXiv

The Compactness of Presupernova Stellar Cores

The success or failure of the neutrino-transport mechanism for producing a supernova in an evolved massive star is known to be sensitive not only to the mass of the iron core that collapses, but also to the density gradient in the silicon and oxygen shells surrounding that core. Here we study the systematics of a presupernova core's "compactness" (O'Connor & Ott 2011) as a function of the mass of the star and the physics used in its calculation. Fine-meshed surveys of presupernova evolution are calculated for stars from 15 to 65 Msun. The metallicity and the efficiency of semiconvection and overshoot mixing are both varied and bare carbon-oxygen cores are explored as well as full hydrogenic stars. Two different codes, KEPLER and MESA, are used for the study. A complex interplay of carbon and oxygen burning, especially in shells, can cause rapid variations in the compactness for stars of very nearly the same mass. On larger scales, the distribution of compactness with main sequence mass is found to be robustly non-monotonic, implying islands of "explodability", particularly around 8 to 20 Msun and 25 to 30 Msun. The carbon-oxygen (CO) core mass of a presupernova star is a better, though still ambiguous discriminant of its core structure than the main sequence mass.

preprint2013arXiv

Type Ia Supernovae from Merging White Dwarfs I. Prompt Detonations

Merging white dwarfs are a possible progenitor of Type Ia supernovae (SNe Ia). While it is not entirely clear if and when an explosion is triggered in such systems, numerical models suggest that a detonation might be initiated before the stars have coalesced to form a single compact object. Here we study such "peri-merger" detonations by means of numerical simulations, modeling the disruption and nucleosynthesis of the stars until the ejecta reach the coasting phase. Synthetic light curves and spectra are generated for comparison with observations. Three models are considered with primary masses 0.96 Msun, 1.06 Msun, and 1.20 Msun. Of these, the 0.96 Msun dwarf merging with an 0.81 Msun companion, with a Ni56 yield of 0.58 Msun, is the most promising candidate for reproducing common SNe Ia. The more massive mergers produce unusually luminous SNe Ia with peak luminosities approaching those attributed to "super-Chandrasekhar" mass SNe Ia. While the synthetic light curves and spectra of some of the models resemble observed SNe Ia, the significant asymmetry of the ejecta leads to large orientation effects. The peak bolometric luminosity varies by more than a factor of 2 with the viewing angle, and the velocities of the spectral absorption features are lower when observed from angles where the light curve is brightest. The largest orientation effects are seen in the ultraviolet, where the flux varies by more than an order of magnitude. Despite the large variation with viewing angle, the set of three models roughly obeys a width-luminosity relation, with the brighter light curves declining more slowly in the B-band. Spectral features due to unburned carbon from the secondary star are also seen in some cases.

preprint2013arXiv

Type Ia Supernovae from Merging White Dwarfs II) Post-Merger Detonations

Merging carbon-oxygen (CO) white dwarfs are a promising progenitor system for Type Ia supernovae (SN Ia), but the underlying physics and timing of the detonation are still debated. If an explosion occurs after the secondary star is fully disrupted, the exploding primary will expand into a dense CO medium that may still have a disk-like structure. This interaction will decelerate and distort the ejecta. Here we carry out multi-dimensional simulations of ``tamped" SN Ia models, using both particle and grid-based codes to study the merger and explosion dynamics, and a radiative transfer code to calculate synthetic spectra and light curves. We find that post-merger explosions exhibit an hourglass-shaped asymmetry, leading to strong variations in the light curves with viewing angle. The two most important factors affecting the outcome are the scale-height of the disk, which depends sensitively on the binary mass ratio, and the total ${}^{56}$Ni yield, which is governed by the central density of the remnant core. The synthetic broadband light curves rise and decline very slowly, and the spectra generally look peculiar, with weak features from intermediate mass elements but relatively strong carbon absorption. We also consider the effects of the viscous evolution of the remnant, and show that a longer time delay between merger and explosion probably leads to larger ${}^{56}$Ni yields and more symmetrical remnants. We discuss the relevance of this class of aspherical ``tamped" SN Ia for explaining the class of ``super-Chandrasekhar'' SN Ia.

preprint2013arXiv

Very Low Energy Supernovae from Neutrino Mass Loss

The continuing difficulty of achieving a reliable explosion in simulations of core-collapse supernovae, especially for more massive stars, has led to speculation concerning the observable transients that might be produced if such a supernova fails. Even if a prompt outgoing shock fails to form in a collapsing presupernova star, one must still consider the hydrodynamic response of the star to the abrupt loss of mass via neutrinos as the core forms a protoneutron star. Following a suggestion by Nadezhin (1980), we calculate the hydrodynamical responses of typical supernova progenitor stars to the rapid loss of approximately 0.2 to 0.5 M_sun of gravitational mass from their centers. In a red supergiant star, a very weak supernova with total kinetic energy ~ 10^47 erg results. The binding energy of a large fraction of the hydrogen envelope before the explosion is of the same order and, depending upon assumptions regarding the neutrino loss rates, most of it is ejected. Ejection speeds are ~ 100 km/s and luminosities ~ 10^39 erg/s are maintained for about a year. A significant part of the energy comes from the recombination of hydrogen. The color of the explosion is extremely red and the events bear some similarity to "luminous red novae," but have much lower speeds.

preprint2012arXiv

Evidence for Type Ia Supernova Diversity from Ultraviolet Observations with the Hubble Space Telescope

We present ultraviolet (UV) spectroscopy and photometry of four Type Ia supernovae (SNe 2004dt, 2004ef, 2005M, and 2005cf) obtained with the UV prism of the Advanced Camera for Surveys on the Hubble Space Telescope. This dataset provides unique spectral time series down to 2000 Angstrom. Significant diversity is seen in the near maximum-light spectra (~ 2000--3500 Angstrom) for this small sample. The corresponding photometric data, together with archival data from Swift Ultraviolet/Optical Telescope observations, provide further evidence of increased dispersion in the UV emission with respect to the optical. The peak luminosities measured in uvw1/F250W are found to correlate with the B-band light-curve shape parameter dm15(B), but with much larger scatter relative to the correlation in the broad-band B band (e.g., ~0.4 mag versus ~0.2 mag for those with 0.8 < dm15 < 1.7 mag). SN 2004dt is found as an outlier of this correlation (at > 3 sigma), being brighter than normal SNe Ia such as SN 2005cf by ~0.9 mag and ~2.0 mag in the uvw1/F250W and uvm2/F220W filters, respectively. We show that different progenitor metallicity or line-expansion velocities alone cannot explain such a large discrepancy. Viewing-angle effects, such as due to an asymmetric explosion, may have a significant influence on the flux emitted in the UV region. Detailed modeling is needed to disentangle and quantify the above effects.

preprint2012arXiv

On the nature of supernovae Ib and Ic

Using non-LTE time-dependent radiative-transfer calculations, we study the impact of mixing and non-thermal processes associated with radioactive decay on SN IIb/Ib/Ic light curves (LCs) and spectra. Starting with short-period binary models of \leq5Msun He-rich stars (18-25Msun on the main-sequence), we produce 1.2B ejecta which we artificially mix to alter the chemical stratification. While the total 56Ni mass influences the LC peak, the spatial distribution of 56Ni, controlled by mixing processes, impacts both the multi-band LCs and spectra. With enhanced mixing, our synthetic LCs start their post-breakout re-brightening phase earlier, follow a more gradual rise to peak, appear redder, and fade faster after peak due to enhanced gamma-ray escape. Non-thermal electrons, crucial for the production of HeI lines, deposit a dominant fraction of their energy as heat. Because energy deposition is generally local well after the LC peak, the broad HeI lines characteristic of maximum-light SN IIb/Ib spectra require mixing that places 56Ni and helium nuclei to within a gamma-ray mean-free-path. This requirement indicates that SNe IIb and Ib most likely arise from the explosion of stripped-envelope massive stars (main-sequence masses \leq25Msun) that have evolved through mass-transfer in a binary system, rather than from more massive single WR stars. In contrast, the lack of HeI lines in SNe Ic may result from a variety of causes: A genuine helium deficiency; strongly-asymmetric mixing; weak mixing; or a more massive, perhaps single, progenitor characterized by a larger oxygen-rich core. Our models, subject to different mixing magnitudes, can produce a variety of SN types, including IIb, IIc, Ib, and Ic. As it is poorly constrained by explosion models, mixing challenges our ability to infer the progenitor and explosion properties of SNe IIb/Ib/Ic.

preprint2011arXiv

Core-collapse explosions of Wolf-Rayet stars and the connection to type IIb/Ib/Ic supernovae

We present non-LTE time-dependent radiative-transfer simulations of supernova (SN) IIb/Ib/Ic spectra and light curves, based on ~1B-energy piston-driven ejecta, with and without 56Ni, produced from single and binary Wolf-Rayet (W-R) stars evolved at solar and sub-solar metallicities. Our bolometric light curves show a 10-day long post-breakout plateau with a luminosity of 1-5x10^7Lsun. In our 56Ni-rich models, with ~3Msun ejecta masses, this plateau precedes a 20-30-day long re-brightening phase initiated by the outward-diffusing heat wave powered by radioactive decay at depth. In low ejecta-mass models with moderate mixing, Gamma-ray leakage starts as early as ~50d after explosion and causes the nebular luminosity to steeply decline by ~0.02mag/d. Such signatures, which are observed in standard SNe IIb/Ib/Ic, are consistent with low-mass progenitors derived from a binary-star population. We propose that the majority of stars with an initial mass ~<20Msun yield SNe II-P if 'effectively" single, SNe IIb/Ib/Ic if part of a close binary system, and SN-less black holes if more massive. Our ejecta, with outer hydrogen mass fractions as low as ~>0.01 and a total hydrogen mass of ~>0.001Msun, yield the characteristic SN IIb spectral morphology at early times. However, by ~15d after the explosion, only Halpha may remain as a weak absorption feature. Our binary models, characterised by helium surface mass fractions of ~>0.85, systematically show HeI lines during the post-breakout plateau, irrespective of the 56Ni abundance. Synthetic spectra show a strong sensitivity to metallicity, which offers the possibility to constrain it directly from SN spectroscopic modelling.

preprint2006arXiv

The Supernova Gamma-Ray Burst Connection

The chief distinction between ordinary supernovae and long-soft gamma-ray bursts (GRBs) is the degree of differential rotation in the inner several solar masses when a massive star dies, and GRBs are rare mainly because of the difficulty achieving the necessary high rotation rate. Models that do provide the necessary angular momentum are discussed, with emphasis on a new single star model whose rapid rotation leads to complete mixing on the main sequence and avoids red giant formation. This channel of progenitor evolution also gives a broader range of masses than previous models, and allows the copious production of bursts outside of binaries and at high redshifts. However, even the production of a bare helium core rotating nearly at break up is not, by itself, a sufficient condition to make a gamma-ray burst. Wolf-Rayet mass loss must be low, and will be low in regions of low metallicity. This suggests that bursts at high redshift (low metallicity) will, on the average, be more energetic, have more time structure, and last longer than bursts nearby. Every burst consists of three components: a polar jet (~0.1 radian), high energy, subrelativistic mass ejection (~1 radian), and low velocity equatorial mass that can fall back after the initial explosion. The relative proportions of these three components can give a diverse assortment of supernovae and high energy transients whose properties may vary with redshift.