Researcher profile

S. E. Woosley

S. E. Woosley contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
27works
0followers
6topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

27 published item(s)

preprint2020arXiv

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

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.

preprint2020arXiv

Three-Dimensional Simulations of Magnetar-Powered Superluminous Supernovae

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.

preprint2011arXiv

From Convection to Explosion: End-to-End Simulation of Type Ia Supernovae

We present our end-to-end capability for computing the convective phase through the explosion phase of Type Ia supernovae. We compute the convective phase up to the time of ignition using our low Mach number code, MAESTRO, and the subsequent explosion phase using our compressible code, CASTRO. Both codes share the same BoxLib software framework and use finite-volume, block-structured adaptive mesh refinement (AMR) to enable high-resolution, three-dimensional full-star simulations that scale to 100,000+ cores. We present preliminary results from the first-ever simulations of convection preceding ignition using MAESTRO with AMR. We also demonstrate our ability to initialize a compressible simulation of the explosion phase in CASTRO using data obtained directly from MAESTRO just before ignition. Some care must be taken during this initialization procedure when interpreting the size and distribution of hot spots.

preprint2011arXiv

Models for Gamma-Ray Burst Progenitors and Central Engines

Most gamma-ray bursts are made during the deaths of massive stars. Here the environmental circumstances, stellar evolutionary paths, and explosion physics that might produce the bursts are reviewed. Neither of the two leading models - collapsar and millisecond magnetar - can be excluded, and both may operate in progenitor stars of different masses, metallicities, and rotation rates. Potential diagnostics are discussed and uncertainties highlighted. Both models are capable of producing a wide variety of transients whose properties vary with both stellar properties and viewing angle. Some of these are reviewed including the possibility of very long (days) low luminosity bursts, so far undiscovered, short hard bursts from massive stellar progenitors, and bursts from very massive Population III stars.

preprint2010arXiv

Discovery of the nearby long, soft GRB 100316D with an associated supernova

We report the Swift discovery of nearby long, soft gamma-ray burst GRB 100316D, and the subsequent unveiling of its low redshift host galaxy and associated supernova. We derive the redshift of the event to be z = 0.0591 +/- 0.0001 and provide accurate astrometry for the GRB-SN. We study the extremely unusual prompt emission with time-resolved gamma-ray to X-ray spectroscopy, and find that the spectrum is best modelled with a thermal component in addition to a synchrotron emission component with a low peak energy. The X-ray light curve has a remarkably shallow decay out to at least 800 s. The host is a bright, blue galaxy with a highly disturbed morphology and we use Gemini South, VLT and HST observations to measure some of the basic host galaxy properties. We compare and contrast the X-ray emission and host galaxy of GRB 100316D to a subsample of GRB-SNe. GRB 100316D is unlike the majority of GRB-SNe in its X-ray evolution, but resembles rather GRB 060218, and we find that these two events have remarkably similar high energy prompt emission properties. Comparison of the host galaxies of GRB-SNe demonstrates, however, that there is a great diversity in the environments in which GRB-SNe can be found. GRB 100316D is an important addition to the currently sparse sample of spectroscopically confirmed GRB-SNe, from which a better understanding of long GRB progenitors and the GRB--SN connection can be gleaned.

preprint2010arXiv

Electron fraction constraints based on Nuclear Statistical Equilibrium with beta equilibrium

The electron-to-nucleon ratio or electron fraction is a key parameter in many astrophysical studies. Its value is determined by weak-interaction rates that are based on theoretical calculations subject to several nuclear physics uncertainties. Consequently, it is important to have a model independent way of constraining the electron fraction value in different astrophysical environments. Here we show that nuclear statistical equilibrium combined with beta equilibrium can provide such a constraint. We test the validity of this approximation in presupernova models and give lower limits for the electron fraction in type Ia supernova and accretion-induced collapse.

preprint2010arXiv

Integrated Nucleosynthesis in Neutrino Driven Winds

Although they are but a small fraction of the mass ejected in core-collapse supernovae, neutrino-driven winds (NDWs) from nascent proto-neutron stars (PNSs) have the potential to contribute significantly to supernova nucleosynthesis. In previous works, the NDW has been implicated as a possible source of r-process and light p-process isotopes. In this paper we present time-dependent hydrodynamic calculations of nucleosynthesis in the NDW which include accurate weak interaction physics coupled to a full nuclear reaction network. Using two published models of PNS neutrino luminosities, we predict the contribution of the NDW to the integrated nucleosynthetic yield of the entire supernova. For the neutrino luminosity histories considered, no true r-process occurs in the most basic scenario. The wind driven from an older $1.4 M_\odot$ model for a PNS is moderately neutron-rich at late times however, and produces $^{87}$Rb, $^{88}$Sr, $^{89}$Y, and $^{90}$Zr in near solar proportions relative to oxygen. The wind from a more recently studied $1.27 M_\odot$ PNS is proton-rich throughout its entire evolution and does not contribute significantly to the abundance of any element. It thus seems very unlikely that the simplest model of the NDW can produce the r-process. At most, it contributes to the production of the N = 50 closed shell elements and some light p-nuclei. In doing so, it may have left a distinctive signature on the abundances in metal poor stars, but the results are sensitive to both uncertain models for the explosion and the masses of the neutron stars involved.

preprint2010arXiv

Modeling the Diversity of Type Ia Supernova Explosions

Type Ia supernovae (SNe Ia) are a prime tool in observational cosmology. A relation between their peak luminosities and the shapes of their light curves allows to infer their intrinsic luminosities and to use them as distance indicators. This relation has been established empirically. However, a theoretical understanding is necessary in order to get a handle on the systematics in SN Ia cosmology. Here, a model reproducing the observed diversity of normal SNe Ia is presented. The challenge in the numerical implementation arises from the vast range of scales involved in the physical mechanism. Simulating the supernova on scales of the exploding white dwarf requires specific models of the microphysics involved in the thermonuclear combustion process. Such techniques are discussed and results of simulations are presented.

preprint2008arXiv

Mixing in Zero and Solar Metallicity Supernovae

Two-dimensional simulations of mixing and fall back in non-rotating massive stars have been carried out using realistic initial models for the presupernova star and assuming standard spherically symmetric explosions of 1.2*10^51ergs. Stars of 15 and 25 solar masses with both primordial and solar composition were modeled. The zero metallicity supernova progenitors were compact blue stars and the amount of Rayleigh-Taylor induced mixing in them was greatly reduced compared with what was seen in the red supergiants with solar metallicity. The compact zero-metal stars also experienced more fallback than their solar metallicity counterparts. As a result, the ejected nucleosynthesis from the two populations was very different. For the simple explosion model assumed, low metallicity stars ejected ejected too little iron and intermediate mass elements even to explain the abundance patterns in the most iron-poor stars found to date, suggesting that some important ingredient is missing. Rotation is likely to alter these conclusions by producing a greater fraction of red supergiants among Population III stars. The velocities of the heavy elements in all models considered - both red and blue supergiants - were less than observed in SN 1987A, suggesting that at least occasionally, asymmetric aspects of the explosion mechanism and fallback play a major role in mixing.

preprint2008arXiv

Nucleosynthesis and Evolution of Massive Metal-Free Stars

The evolution and explosion of metal-free stars with masses 10--100 solar masses are followed, and their nucleosynthetic yields, light curves, and remnant masses determined. When the supernova yields are integrated over a Salpeter initial mass function, the resulting elemental abundance pattern is qualitatively solar, but with marked deficiencies of odd-Z elements with 7 <= Z <= 13. Neglecting the contribution of the neutrino wind from the neutron stars that they make, no appreciable abundances are made for elements heavier than germanium. The computed pattern compares favorably with what has been observed in metal-deficient stars with [Z] ~< -3. Most of the stars end their lives as blue supergiants and make supernovae with distinctive light curves resembling SN 1987A, but some produce primary nitrogen by dredge up and become red supergiants. A novel automated fitting algorithm is developed for determining optimal combinations of explosion energy, mixing, and initial mass function in the large model data base to agree with specified data sets. The model is applied to the low metallicity sample of Cayrel et al. (2004) and the two ultra-iron-poor stars HE0107-5240 and HE1327-2326. Best agreement with these low metallicity stars is achieved with very little mixing, and none of the metal-deficient data sets considered show the need for a high energy explosion component. To the contrary, explosion energies somewhat less than 1.2 B seem to be preferred in most cases. (abbreviated)

preprint2008arXiv

Nucleosynthesis in Early Neutrino Driven Winds

Nucleosynthesis in early neutrino winds is investigated. Presented is a brief overview of two recent problems of supernova nucleosynthesis. In the first part we investigate the effect of nuclear parameters on the synthesis of Mo92 and Mo94. Based on recent experimental results, we find that the proton rich winds of the model investigated here can not be the source of solar Mo92 and Mo94. In the second part we investigate the nucleosynthesis from neutron rich bubbles and show that they do not contribute to the overall nucleosynthesis.

preprint2007arXiv

Type Ia Supernova: Burning and Detonation in the Distributed Regime

A simple, semi-analytic representation is developed for nuclear burning in Type Ia supernovae in the special case where turbulent eddies completely disrupt the flame. The speed and width of the ``distributed&#39;&#39; flame front are derived. For the conditions considered, the burning front can be considered as a turbulent flame brush composed of corrugated sheets of well-mixed flames. These flames are assumed to have a quasi-steady-state structure similar to the laminar flame structure, but controlled by turbulent diffusion. Detonations cannot appear in the system as long as distributed flames are still quasi-steady-state, but this condition is violated when the distributed flame width becomes comparable to the size of largest turbulent eddies. When this happens, a transition to detonation may occur. For current best estimates of the turbulent energy, the most likely density for the transition to detonation is in the range 0.5 - 1.5 x 10^7 g cm^{-3}.

preprint2006arXiv

A Radial Velocity Survey of the Cygnus OB2 Association

We conducted a radial velocity survey of the Cygnus OB2 Association over a 6 year (1999 - 2005) time interval to search for massive close binaries. During this time we obtained 1139 spectra on 146 OB stars to measure mean systemic radial velocities and radial velocity variations. We spectroscopically identify 73 new OB stars for the first time, the majority of which are likely to be Association members. Spectroscopic evidence is also presented for a B3Iae classification and temperature class variation (B3 - B8) on the order of 1 year for Cygnus OB2 No. 12. Calculations of the intial mass function with the current spectroscopic sample yield Gamma = -2.2 +/- 0.1. Of the 120 stars with the most reliable data, 36 are probable and 9 are possible single-lined spectroscopic binaries. We also identify 3 new and 8 candidate double-lined spectroscopic binaries. These data imply a lower limit on the massive binary fraction of 30% - 42%. The calculated velocity dispersion for Cygnus OB2 is 2.44 +/- km/s, which is typical of open clusters. No runaway OB stars were found.

preprint2006arXiv

Constraining deflagration models of Type Ia supernovae through intermediate-mass elements

The physical structure of a nuclear flame is a basic ingredient of the theory of Type Ia supernovae (SNIa). Assuming an exponential density reduction with several characteristic times we have followed the evolution of a planar nuclear flame in an expanding background from an initial density 6.6 10^7 g/cm3 down to 2 10^6 g/cm3. The total amount of synthesized intermediate-mass elements (IME), from silicon to calcium, was monitored during the calculation. We have made use of the computed mass fractions, X_IME, of these elements to give an estimation of the total amount of IME synthesized during the deflagration of a massive white dwarf. Using X_IME and adopting the usual hypothesis that turbulence decouples the effective burning velocity from the laminar flame speed, so that the relevant flame speed is actually the turbulent speed on the integral length-scale, we have built a simple geometrical approach to model the region where IME are thought to be produced. It turns out that a healthy production of IME involves the combination of not too short expansion times, t_c > 0.2 s, and high turbulent intensities. According to our results it could be difficult to produce much more than 0.2 solar masses of intermediate-mass elements within the deflagrative paradigma. The calculations also suggest that the mass of IME scales with the mass of Fe-peak elements, making it difficult to conciliate energetic explosions with low ejected nickel masses, as in the well observed SN1991bg or in SN1998de. Thus a large production of Si-peak elements, especially in combination with a low or a moderate production of iron, could be better addressed by either the delayed detonation route in standard Chandrasekhar-mass models or, perhaps, by the off-center helium detonation in the sub Chandrasekhar-mass scenario.

preprint2006arXiv

Off-center ignition in type Ia supernova: I. Initial evolution and implications for delayed detonation

The explosion of a carbon-oxygen white dwarf as a Type Ia supernova is known to be sensitive to the manner in which the burning is ignited. Studies of the pre-supernova evolution suggest asymmetric, off-center ignition, and here we explore its consequences in two- and three-dimensional simulations. Compared with centrally ignited models, one-sided ignitions initially burn less and release less energy. For the distributions of ignition points studied, ignition within two hemispheres typically leads to the unbinding of the white dwarf, while ignition within a small fraction of one hemisphere does not. We also examine the spreading of the blast over the surface of the white dwarf that occurs as the first plumes of burning erupt from the star. In particular, our studies test whether the collision of strong compressional waves can trigger a detonation on the far side of the star as has been suggested by Plewa et al. (2004). The maximum temperature reached in these collisions is sensitive to how much burning and expansion has already gone on, and to the dimensionality of the calculation. Though detonations are sometimes observed in 2D models, none ever happens in the corresponding 3D calculations. Collisions between the expansion fronts of multiple bubbles also seem, in the usual case, unable to ignite a detonation. &#34;Gravitationally confined detonation&#34; is therefore not a robust mechanism for the explosion. Detonation may still be possible in these models however, either following a pulsation or by spontaneous detonation if the turbulent energy is high enough.

preprint2006arXiv

Prompt and Afterglow Emission Properties of Gamma-Ray Bursts with Spectroscopically Identified Supernovae

We present a detailed spectral analysis of the prompt and afterglow emission of four nearby long-soft gamma-ray bursts (GRBs 980425, 030329, 031203, and 060218) that were spectroscopically found to be associated with type Ic supernovae, and compare them to the general GRB population. For each event, we investigate the spectral and luminosity evolution, and estimate the total energy budget based upon broadband observations. The observational inventory for these events has become rich enough to allow estimates of their energy content in relativistic and sub-relativistic form. The result is a global portrait of the effects of the physical processes responsible for producing long-soft GRBs. In particular, we find that the values of the energy released in mildly relativistic outflows appears to have a significantly smaller scatter than those found in highly relativistic ejecta. This is consistent with a picture in which the energy released inside the progenitor star is roughly standard, while the fraction of that energy that ends up in highly relativistic ejecta outside the star can vary dramatically between different events.

preprint2006arXiv

Surface detonation in type Ia supernova explosions?

We explore the evolution of thermonuclear supernova explosions when the progenitor white dwarf star ignites asymmetrically off-center. Several numerical simulations are carried out in two and three dimensions to test the consequences of different initial flame configurations such as spherical bubbles displaced from the center, more complex deformed configurations, and teardrop-shaped ignitions. The burning bubbles float towards the surface while releasing energy due to the nuclear reactions. If the energy release is too small to gravitationally unbind the star, the ash sweeps around it, once the burning bubble approaches the surface. Collisions in the fuel on the opposite side increase its temperature and density and may -- in some cases -- initiate a detonation wave which will then propagate inward burning the core of the star and leading to a strong explosion. However, for initial setups in two dimensions that seem realistic from pre-ignition evolution, as well as for all three-dimensional simulations the collimation of the surface material is found to be too weak to trigger a detonation.

preprint2006arXiv

Swift detects a remarkable gamma-ray burst, GRB 060614, that introduces a new classification scheme

Gamma ray bursts (GRBs) are known to come in two duration classes, separated at ~2 s. Long bursts originate from star forming regions in galaxies, have accompanying supernovae (SNe) when near enough to observe and are likely caused by massive-star collapsars. Recent observations show that short bursts originate in regions within their host galaxies with lower star formation rates consistent with binary neutron star (NS) or NS - black hole (BH) mergers. Moreover, although their hosts are predominantly nearby galaxies, no SNe have been so far associated with short GRBs. We report here on the bright, nearby GRB 060614 that does not fit in either class. Its ~102 s duration groups it with long GRBs, while its temporal lag and peak luminosity fall entirely within the short GRB subclass. Moreover, very deep optical observations exclude an accompanying supernova, similar to short GRBs. This combination of a long duration event without accompanying SN poses a challenge to both a collapsar and merging NS interpretation and opens the door on a new GRB classification scheme that straddles both long and short bursts.

preprint2006arXiv

Type Ia Supernova Light Curves

The diversity of Type Ia supernova (SN Ia) photometry is explored using a grid of 130 one-dimensional models. It is shown that the observable properties of SNe Ia resulting from Chandrasekhar-mass explosions are chiefly determined by their final composition and some measure of ``mixing&#39;&#39; in the explosion. A grid of final compositions is explored including essentially all combinations of 56Ni, stable ``iron&#39;&#39;, and intermediate mass elements that result in an unbound white dwarf. Light curves (and in some cases spectra) are calculated for each model using two different approaches to the radiation transport problem. Within the resulting templates are models that provide good photometric matches to essentially the entire range of observed SNe Ia. On the whole, the grid of models spans a wide range in B-band peak magnitudes and decline rates, and does not obey a Phillips relation. In particular, models with the same mass of 56Ni show large variations in their light curve decline rates. We identify the physical parameters responsible for this dispersion, and consider physically motivated ``cuts&#39;&#39; of the models that agree better with the Phillips relation. For example, models that produce a constant total mass of burned material of 1.1 +/- Msun do give a crude Phillips relation, albeit with much scatter. The scatter is further reduced if one restricts that set to models that make 0.1 to 0.3 Msun of stable iron and nickel isotopes, and then mix the ejecta strongly between the center and 0.8 Msun. We conclude that the supernovae that occur most frequently in nature are highly constrained by the Phillips relation and that a large part of the currently observed scatter in the relation is likely a consequence of the intrinsic diversity of these objects.

preprint2005arXiv

Discovery of the short gamma-ray burst GRB 050709

Gamma-Ray Bursts (GRBs) fall into two classes: short-hard and long-soft bursts. The latter are now known to have X-ray and optical afterglows, to occur at cosmological distances in star-forming galaxies, and to be associated with the explosion of massive stars. In contrast, the distance scale, the energy scale, and the progenitors of short bursts have remained a mystery. Here we report the discovery of a short-hard burst whose accurate localization has led to follow-up observations that have identified the X-ray afterglow and (for the first time) the optical afterglow of a short-hard burst. These, in turn, have led to identification of the host galaxy of the burst as a late-type galaxy at z=0.16 showing that at least some short-hard bursts occur at cosmological distances in the outskirts of galaxies, and are likely to be due to the merging of compact binaries.

preprint2003arXiv

Neutrino Nucleosynthesis

We study neutrino process nucleosynthesis in massive stars using newly calculated cross sections, an expanded reaction network, and complete and self-consistent models of the progenitor star. We reevaluate the production of light isotopes from abundant progenitors as well as that of rare, heavy, proton-rich isotopes. In particular, new results are given for B11, F19, La138, and Ta180. The production of these isotopes places limits on neutrino spectrum and oscialltions.

preprint2000arXiv

Pair-Instability Supernovae, Gravity Waves, and Gamma-Ray Transients

Growing theoretical evidence suggests that the first generation of stars may have been quite massive (~100-300 solar masses). If they retain their high mass until death, such stars will, after about 3Myr, make pair-instability supernovae. We consider the complete evolution of two zero-metallicity stars of 250 and 300 solar masses. Explosive oxygen and silicon burning cause the 130 solar mass helium core to explode, but explosive burning is unable to drive an explosion in the 300 solar mass star and it collapses to a black hole. For this star, the calculated angular momentum in the presupernova model is sufficient to delay black hole formation and the star initially forms a 50 solar mass, 1000km core within which neutrinos are trapped. Although the star does not become dynamically unstable, the calculated growth time of secular rotational instabilities is shorter than the black hole formation time, and such instabilities may develop. We estimate the energy and amplitude of the gravitational waves emitted during this collapse. After the black hole forms, accretion continues through a disk. Although the disk is far too large and cool to transport energy efficiently to the rotational axis by neutrino annihilation, it has ample potential energy to produce a 1e54erg jet driven by magnetic fields. The interaction of this jet with surrounding circumstellar gas may produce an energetic gamma-ray transient, but given the redshift and time scale, this is probably not a model for typical gamma-ray bursts.

preprint2000arXiv

Prediction of Astrophysical Reaction Rates: Methods, Data Needs, and Consequences for Nucleosynthesis Studies

The majority of nuclear reactions in astrophysics involve unstable nuclei which are not fully accessible by experiments yet. Therefore, there is high demand for reliable predictions of cross sections and reaction rates by theoretical means. The majority of reactions can be treated in the framework of the statistical model (Hauser-Feshbach). The global parametrizations of the nuclear properties needed for predictions far off stability probe our understanding of the strong force and take it to its limit. The sensitivity of astrophysical scenarios to nuclear inputs is illustrated in the framework of a detailed nucleosynthesis study in type II supernovae. Abundances resulting from calculations in the same explosion model with two different sets of reaction rates are compared. Key reactions and required nuclear information are identified.

preprint1999arXiv

Hyper-Accreting Black Holes and Gamma-Ray Bursts

A variety of current models for gamma-ray bursts (GRBs) suggest a common engine - a black hole of several solar masses accreting matter from a disk at a rate 0.01 to 10 solar masses per second. Using a numerical model for relativistic disk accretion, we have studied steady-state accretion at these high rates. Inside a radius ~ 10**8 cm, for accretion rates greater than about 0.01 solar masses per second, a global state of balanced power comes to exist between neutrino losses, chiefly pair capture on nucleons, and dissipation. Energy emitted in neutrinos is less, and in the case of low accretion rates, very much less, than the maximum efficiency factor for black hole accretion (0.057 for no rotation; 0.42 for extreme Kerr rotation) times Mdot c**2. The efficiency for producing a pair fireball along the rotational axis by neutrino annihilation is calculated and found to be highly variable and very sensitive to the accretion rate. For some of the higher accretion rates studied, it can be several per cent or more; for accretion rates less than 0.05 solar masses per second, it is essentially zero. The efficiency of the Blandford-Znajek mechanism in extracting rotational energy from the black hole is also estimated. In light of these results, the viability of various gamma-ray burst models is discussed and the sensitivity of the results to disk viscosity, black hole rotation rate, and black hole mass explored. A diverse range of GRB energies seems unavoidable and neutrino annihilation in hyper-accreting black hole systems can explain bursts up to 10**52 erg. Larger energies may be inferred for beaming systems.

preprint1996arXiv

Nucleosynthesis in neutrino-driven winds: II. Implications for heavy element synthesis

During the first 20 seconds of its life, the enormous neutrino luminosity of a neutron star drives appreciable mass loss from its surface. This neutrino-driven wind has been previously identified as a likely site for the r-process. Qian & Woosley (1996) have derived, both analytically and numerically, the physical conditions relevant for heavy element synthesis in the wind. These conditions include the entropy (S), the electron fraction (Ye), the dynamic time scale, and the mass loss rate. Here we explore the implications of these conditions for nucleosynthesis. We find that the standard wind models derived in that paper are inadequate to make the r-process, though they do produce some rare species above the iron group. We further determine the general restrictions on the entropy, the electron fraction, and the dynamic time scale that are required to make the r-process. In particular, we derive from nuclear reaction network calculations the conditions required to give a sufficient neutron-to-seed ratio for production of the platinum peak. These conditions range from Ye = 0.2 and S < 100 per baryon for reasonable dynamic time scales of 0.001-0.1 s, to Ye = 0.4-0.495 and S > 400 per baryon for a dynamic time scale of 0.1 s. These conditions are also derived analytically to illustrate the physics determining the neutron-to-seed ratio.

preprint1995arXiv

The Neutron Star and Black Hole Initial Mass Function

Using recently calculated models for massive stellar evolution and supernovae coupled to a model for Galactic chemical evolution, neutron star and black hole birth functions (number of neutron stars and black holes as a function of their mass) are determined for the Milky Way Galaxy. For those stars that explode as Type II supernovae, the models give birth functions that are bimodal with peaks at 1.27 and 1.76 M\sun \ and average masses within those peaks of 1.28 and 1.73 M\sun. For those stars that explode as Type Ib there is a narrower spread of remnant masses, the average being 1.32 M\sun, and less evidence for bimodality. These values will be increased, especially in the more massive Type II supernovae, if significant accretion continues during the initial launching of the shock, and the number of heavier neutron stars could be depleted by black hole formation. The principal reason for the dichotomy in remnant masses for Type II is the difference in the presupernova structure of stars above and below 19 M\sun, the mass separating stars that burn carbon convectively from those that produce less carbon and burn radiatively. The Type Ib&#39;s and the lower mass group of the Type II&#39;s compare favorably with measured neutron star masses, and in particular to the Thorsett et al. (1993) determination of the average neutron star mass in 17 systems; 1.35 $\pm$ 0.27 M\sun. Variations in the exponent of a Salpeter initial mass function are shown not to affect the locations of the two peaks in the distribution function, but do affect their relative amplitudes. Sources of uncertainty, in particular placement of the mass cut and sensitivity to the explosion energy, are discussed, and estimates of the total number of neutron stars and black holes in the Galaxy are given. Accretion induced collapse should give a unique gravitational mass of 1.27