BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

S. E. Woosley

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

ResearcherUnclaimed source record
astro-ph.HEastro-phastro-ph.SRastro-ph.COnucl-thastro-ph.IMhep-phnucl-ex

Catalog footprint

What is connected

56works
8topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

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

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Progenitor of GW 150914

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Remarkable Deaths of 9 - 11 Solar Mass Stars

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Deaths of Very Massive Stars

The theory underlying the evolution and death of stars heavier than 10 Msun on the main sequence is reviewed with an emphasis upon stars much heavier than 30 Msun. These are stars that, in the absence of substantial mass loss, are expected to either produce black holes when they die, or, for helium cores heavier than about 35 Msun, encounter the pair instability. A wide variety of outcomes is possible depending upon the initial composition of the star, its rotation rate, and the physics used to model its evolution. These heavier stars can produce some of the brightest supernovae in the universe, but also some of the faintest. They can make gamma-ray bursts or collapse without a whimper. Their nucleosynthesis can range from just CNO to a broad range of elements up to the iron group. Though rare nowadays, they probably played a disproportionate role in shaping the evolution of the universe following the formation of its first stars.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Carbon Deflagration in Type Ia Supernova: I. Centrally Ignited Models

A leading model for Type Ia supernovae (SNe Ia) begins with a white dwarf near the Chandrasekhar mass that ignites a degenerate thermonuclear runaway close to its center and explodes. In a series of papers, we shall explore the consequences of ignition at several locations within such dwarfs. Here we assume central ignition, which has been explored before, however, the problem is worth revisiting, if only to validate those previous studies and to further elucidate the relevant physics for future work. A perturbed sphere of hot iron ash with a radius of ~100 km is initialized at the middle of the star. The subsequent explosion is followed in several simulations using a thickened flame model in which the flame speed is either fixed --- within the range expected from turbulent combustion --- or based on the local turbulent intensity. Global results, including the explosion energy and bulk nucleosynthesis (e.g. 56Ni of 0.48--0.56 $\Msun$) turn out to be insensitive to this speed. In all completed runs, the energy released by the nuclear burning is adequate to unbind the star, but not enough to give the energy and brightness of typical SNe Ia. As found previously, the chemical stratification observed in typical events is not reproduced. These models produce a large amount of unburned carbon and oxygen in central low velocity regions, which is inconsistent with spectroscopic observations, and the intermediate mass elements and iron group elements are strongly mixed during the explosion.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Finding the First Cosmic Explosions I: Pair-Instability Supernovae

The first stars are the key to the formation of primitive galaxies, early cosmological reionization and chemical enrichment, and the origin of supermassive black holes. Unfortunately, in spite of their extreme luminosities, individual Population III stars will likely remain beyond the reach of direct observation for decades to come. However, their properties could be revealed by their supernova explosions, which may soon be detected by a new generation of NIR observatories such as JWST and WFIRST. We present light curves and spectra for Pop III pair-instability supernovae calculated with the Los Alamos radiation hydrodynamics code RAGE. Our numerical simulations account for the interaction of the blast with realistic circumstellar envelopes, the opacity of the envelope, and Lyman absorption by the neutral IGM at high redshift, all of which are crucial to computing the NIR signatures of the first cosmic explosions. We find that JWST will detect pair-instability supernovae out to z > 30, WFIRST will detect them in all-sky surveys out to z ~ 15 - 20 and LSST and Pan-STARRS will find them at z ~ 7 - 8. The discovery of these ancient explosions will probe the first stellar populations and reveal the existence of primitive galaxies that might not otherwise have been detected.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Finding the First Cosmic Explosions. III. Pulsational Pair-Instability Supernovae

Population III supernovae have been the focus of growing attention because of their potential to directly probe the properties of the first stars, particularly the most energetic events that can be seen at the edge of the observable universe. But until now pair-pulsation supernovae, in which explosive thermonuclear burning in massive stars fails to unbind them but can eject their outer layers into space, have been overlooked as cosmic beacons at the earliest redshifts. These shells can later collide and, like Type IIn supernovae, produce superluminous events in the UV at high redshifts that could be detected in the near infrared today. We present numerical simulations of a 110 M$_{\odot}$ pair-pulsation explosion done with the Los Alamos radiation hydrodynamics code RAGE. We find that collisions between consecutive pair pulsations are visible in the near infrared out to z $\sim$ 15 - 20 and can probe the earliest stellar populations at cosmic dawn.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Nucleosynthetic Constraints on the Mass of the Heaviest Supernovae

Assuming a Salpeter initial mass function and taking the solar abundances as a representative sample, we explore the sensitivity of nucleosynthesis in massive stars to the truncation of supernova explosions above a certain mass. It is assumed that stars of all masses contribute to nucleosynthesis by their pre-explosive winds, but above a certain limiting main sequence mass, the presupernova star becomes a black hole and ejects nothing more. The solar abundances from oxygen to atomic mass 90 are fit quite well assuming no cut-off at all, i.e., by assuming all stars up to 120 solar masses make successful supernovae. Little degradation in the fit occurs if the upper limit is reduced to 25 solar masses. The limit can be further reduced, but the required event rate of supernovae in the remaining range rises rapidly to compensate for the lost nucleosynthesis of the more massive stars. The nucleosynthesis of the s-process declines precipitously and the production of species made in the winds, e.g., carbon, becomes unacceptably large compared with elements made in the explosion, e.g., silicon and oxygen. However, by varying uncertain physics, especially the mass loss rate for massive stars and the rate for the neon-22 to magnesium-25 reaction rate, acceptable nucleosynthesis might still be achieved with a cutoff as low as 18 solar masses. This would require a supernova frequency three times greater than the fiducial value obtained when all stars explode in order to produce the required oxygen-16. The nucleosynthesis of iron-60 and aluminum-26 is also examined.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Deflagration Stage of Chandrasekhar Mass Models For Type Ia Supernovae: I. Early Evolution

We present high-resolution, full-star simulations of the post-ignition phase of Type Ia supernovae using the compressible hydrodynamics code Castro. Initial conditions, including the turbulent velocity field and ignition site, are imported directly from a simulation of the last few hours of presupernova convection using a low Mach number code, Maestro. Adaptive mesh refinement allows the initial burning front to be modeled with an effective resolution of 36,864^3 zones (~136 m/zone). The initial rise and expansion of the deflagration front are tracked until burning reaches the star's edge and the role of the background turbulence on the flame is investigated. The effect of artificially moving the ignition location closer to the star's center is explored. The degree to which turbulence affects the burning front decreases with increasing ignition radius since the buoyancy force is stronger at larger radii. Even central ignition --- in the presence of a background convective flow field --- is rapidly carried off-center as the flame is carried by the flow field. We compare our results to analytic models for burning thermals, and find that they reproduce the general trends of the bubble's size and mass, but underpredict the amount of buoyant acceleration due to simplifying assumptions of the bubble's properties. Overall, we find that the amount of mass that burns prior to flame break out is small, consistent with a "gravitationally confined detonation" occurring at a later epoch, but additional burning will occur following breakout that may modify this conclusion.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Seeing the First Supernovae at the Edge of the Universe with JWST

The first stars ended the cosmic Dark Ages and created the first heavy elements necessary for the formation of planets and life. The properties of these stars remain uncertain, and it may be decades before individual Pop III stars are directly observed. Their masses, however, can be inferred from their supernova explosions, which may soon be found in both deep-field surveys by JWST and in all-sky surveys by WFIRST. We have performed radiation hydrodynamical simulations of the near infrared signals of Pop III pair-instability supernovae in realistic circumstellar environments with Lyman absorption by the neutral intergalactic medium. We find that JWST and WFIRST will detect these explosions out to z ~ 30 and 20, respectively, unveiling the first generation of stars in the universe.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

High-Resolution Simulations of Convection Preceding Ignition in Type Ia Supernovae Using Adaptive Mesh Refinement

We extend our previous three-dimensional, full-star simulations of the final hours of convection preceding ignition in Type Ia supernovae to higher resolution using the adaptive mesh refinement capability of our low Mach number code, MAESTRO. We report the statistics of the ignition of the first flame at an effective 4.34 km resolution, and general flow field properties at an effective 2.17 km resolution. We find that off-center ignition is likely, with radius of 50 km most favored and a likely range of 40 to 75 km. This is consistent with our previous coarser (8.68 km resolution) simulations, implying that we have achieved sufficient resolution in our determination of likely ignition radii. The dynamics of the last few hot spots preceding ignition suggest that a multiple ignition scenario is not likely. With improved resolution, we can more clearly see the general flow pattern in the convective region, characterized by a strong outward plume with a lower speed recirculation. We show that the convective core is turbulent with a Kolmogorov spectrum and has a lower turbulent intensity and larger integral length scale than previously thought (on the order of 16 km s$^{-1}$ and 200 km, respectively), and we discuss the potential consequences for the first flames.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Long Gamma-Ray Transients from Collapsars

In the collapsar model for common gamma-ray bursts, the formation of a centrifugally supported disk occurs during the first $\sim$10 seconds following the collapse of the iron core in a massive star. This only occurs in a small fraction of massive stellar deaths, however, and requires unusual conditions. A much more frequent occurrence could be the death of a star that makes a black hole and a weak or absent outgoing shock, but in a progenitor that only has enough angular momentum in its outermost layers to make a disk. We consider several cases where this is likely to occur - blue supergiants with low mass loss rates, tidally-interacting binaries involving either helium stars or giant stars, and the collapse to a black hole of very massive pair-instability supernovae. These events have in common the accretion of a solar mass or so of material through a disk over a period much longer than the duration of a common gamma-ray burst. A broad range of powers is possible, $10^{47}$ to $10^{50}\,$erg s$^{-1}$, and this brightness could be enhanced by beaming. Such events were probably more frequent in the early universe where mass loss rates were lower. Indeed this could be one of the most common forms of gamma-ray transients in the universe and could be used to study first generation stars. Several events could be active in the sky at any one time. A recent example of this sort of event may have been the SWIFT transient Sw-1644+57.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Pair Instability Supernovae: Light Curves, Spectra, and Shock Breakout

For the initial mass range (140 < M < 260 Msun) stars die in a thermonuclear runaway triggered by the pair-production instability. The supernovae they make can be remarkably energetic (up to ~10^53 ergs) and synthesize considerable amounts of radioactive isotopes. Here we model the evolution, explosion, and observational signatures of representative pair-instability supernovae (PI SNe) spanning a range of initial masses and envelope structures. The predicted light curves last for hundreds of days and range in luminosity, from very dim to extremely bright, L ~ 10^44 ergs/s. The most massive events are bright enough to be seen at high redshift, but the extended light curve duration (~1 year) -- prolonged by cosmological time-dilation -- may make it difficult to detect them as transients. An alternative approach may be to search for the brief and luminous outbreak occurring when the explosion shock wave reaches the stellar surface. Using a multi-wavelength radiation-hydrodynamics code we calculate that, in the rest-frame, the shock breakout transients of PI SNe reach luminosities of 10^45-10^46 ergs/s, peak at wavelengths ~30-170 Angstroms, and last for several hours. We explore the detectability of PI SNe emission at high redshift, and discuss how observations of the light curves, spectra, and breakout emission can be used to constrain the mass, radius, and metallicity of the progenitor.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Sub-Chandrasekhar Mass Models For Type Ia Supernovae

For carbon-oxygen white dwarfs accreting hydrogen or helium at rates in the range ~1-10 x 10^(-8) Msun/y, a variety of explosive outcomes is possible well before the star reaches the Chandrasekhar mass. These outcomes are surveyed for a range of white dwarf masses (0.7 - 1.1 Msun), accretion rates (1 - 7 x 10^(-8) Msun/y), and initial white dwarf temperatures (0.01 and 1 Lsun). The results are particularly sensitive to the convection that goes on during the last few minutes before the explosion. Unless this convection maintains a shallow temperature gradient, and unless the density is sufficiently high, the accreted helium does not detonate. Below a critical helium ignition density, which we estimate to be 5 - 10 x 10^5 g cm^(-3), either helium novae or helium deflagrations result. The hydrodynamics, nucleosynthesis, light curves, and spectra of a representative sample of detonating and deflagrating models are explored. Some can be quite faint indeed, powered at peak for a few days by the decay of 48Cr and 48V. Only the hottest, most massive white dwarfs considered with the smallest helium layers, show reasonable agreement with the light curves and spectra of common Type Ia supernovae. For the other models, especially those involving lighter white dwarfs, the helium shell mass exceeds 0.05 Msun and the mass of the 56Ni that is synthesized exceeds 0.01 Msun. These explosions do not look like ordinary Type Ia supernovae, or any other frequently observed transient.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

XRF 100316D / SN 2010bh and the Nature of Gamma Ray Burst Supernovae

We present ground-based and HST optical and infrared observations of XRF 100316D / SN 2010bh. It is seen that the optical light curves of SN 2010bh evolve at a faster rate than the archetype GRB-SN 1998bw, but at a similar rate to SN 2006aj, a supernova that was spectroscopically linked with XRF 060218, and at a similar rate to non-GRB associated type Ic SN 1994I. We estimate the rest-frame extinction of this event from our optical data to be E(B-V)=0.18 +/- 0.08 mag. We find the V-band absolute magnitude of SN 2010bh to be M_{V}=-18.62 +/- 0.08, which is the faintest peak V-band magnitude observed to-date for a spectroscopically-confirmed GRB-SNe. When we investigate the origin of the flux at t-t_{o}=0.598 days, it is shown that the light is not synchrotron in origin, but is likely coming from the supernova shock break-out. We then use our optical and infrared data to create a quasi-bolometric light curve of SN 2010bh which we model with a simple analytical formula. The results of our modeling imply that SN 2010bh synthesized a nickel mass of M_{Ni} \approx 0.10 M_{sun}, ejected M_{ej} \approx 2.2 M_{sun} and has an explosion energy of E_{k} \approx 1.4 x 10^{52} erg. Finally, for a sample 22 GRB-SNe we check for a correlation between the stretch factors and luminosity factors in the R band and conclude that no statistically-significant correlation exists.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

3D simulations of Rayleigh-Taylor mixing in core-collapse SNe with CASTRO

We present multidimensional simulations of the post-explosion hydrodynamics in three different 15 solar mass supernova models with zero, 10^{-4} solar metallicity, and solar metallicities. We follow the growth of the Rayleigh-Taylor instability that mixes together the stellar layers in the wake of the explosion. Models are initialized with spherically symmetric explosions and perturbations are seeded by the grid. Calculations are performed in two-dimensional axisymmetric and three-dimensional Cartesian coordinates using the new Eulerian hydrodynamics code, CASTRO. We find as in previous work, that Rayleigh-Taylor perturbations initially grow faster in 3D than in 2D. As the Rayleigh-Taylor fingers interact with one another, mixing proceeds to a greater degree in 3D than in 2D, reducing the local Atwood number and slowing the growth rate of the instability in 3D relative to 2D. By the time mixing has stopped, the width of the mixed region is similar in 2D and 3D simulations provided the Rayleigh-Taylor fingers show significant interaction. Our results imply that 2D simulations of light curves and nucleosynthesis in supernovae (SNe) that die as red giants may capture the features of an initially spherically symmetric explosion in far less computational time than required by a full 3D simulation. However, capturing large departures from spherical symmetry requires a significantly perturbed explosion. Large scale asymmetries cannot develop through an inverse cascade of merging Rayleigh-Taylor structures; they must arise from asymmetries in the initial explosion.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A tale of two GRB-SNe at a common redshift of z = 0.54

We present ground-based and HST optical observations of the optical transients (OTs) of long-duration Gamma Ray Bursts (GRBs) 060729 and 090618, both at a redshift of z = 0.54. For GRB 060729, bumps are seen in the optical light curves (LCs), and the late-time broadband spectral energy distributions (SEDs) of the OT resemble those of local type Ic supernovae (SNe). For GRB 090618, the dense sampling of our optical observations has allowed us to detect well-defined bumps in the optical LCs, as well as a change in colour, that are indicative of light coming from a core-collapse SN. The accompanying SNe for both events are individually compared with SN1998bw, a known GRB-supernova, and SN1994I, a typical type Ic supernova without a known GRB counterpart, and in both cases the brightness and temporal evolution more closely resemble SN1998bw. We also exploit our extensive optical and radio data for GRB 090618, as well as the publicly-available SWIFT -XRT data, and discuss the properties of the afterglow at early times. In the context of a simple jet-like model, the afterglow of GRB 090618 is best explained by the presence of a jet-break at t-to > 0.5 days. We then compare the rest-frame, peak V -band absolute magnitudes of all of the GRB and X-Ray Flash (XRF)-associated SNe with a large sample of local type Ibc SNe, concluding that, when host extinction is considered, the peak magnitudes of the GRB/XRF-SNe cannot be distinguished from the peak magnitudes of non-GRB/XRF SNe.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Spectral Lag Relations in GRB Pulses Detected with HETE-2

Using a pulse-fit method, we investigate the spectral lags between the traditional gamma-ray band (50-400 keV) and the X-ray band (6-25 keV) for 8 GRBs with known redshifts (GRB 010921, GRB 020124, GRB 020127, GRB 021211, GRB 030528, GRB 040924, GRB 041006, GRB 050408) detected with the WXM and FREGATE instruments aboard the HETE-2 satellite. We find several relations for the individual GRB pulses between the spectral lag and other observables, such as the luminosity, pulse duration, and peak energy (Epeak). The obtained results are consistent with those for BATSE, indicating that the BATSE correlations are still valid at lower energies (6-25 keV). Furthermore, we find that the photon energy dependence for the spectral lags can reconcile the simple curvature effect model. We discuss the implication of these results from various points of view.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Nucleosynthetic Imprint of 15-40 Solar Mass Primordial Supernovae on Metal-Poor Stars

The inclusion of rotationally-induced mixing in stellar evolution can alter the structure and composition of presupernova stars. We survey the effects of progenitor rotation on nucleosynthetic yields in Population III and II supernovae using the new adaptive mesh refinement (AMR) code CASTRO. We examine spherical explosions in 15, 25 and 40 solar mass stars at Z = 0 and 10^-4 solar metallicity with three explosion energies and two rotation rates. Rotation in the Z = 0 models resulted in primary nitrogen production and a stronger hydrogen burning shell which led all models to die as red supergiants. On the other hand, the Z=10^-4 solar metallicity models that included rotation ended their lives as compact blue stars. Because of their extended structure, the hydrodynamics favors more mixing and less fallback in the metal free stars than the Z = 10^-4 models. As expected, higher energy explosions produce more enrichment and less fallback than do lower energy explosions, and less massive stars produce more enrichment and leave behind smaller remnants than do more massive stars. We compare our nucleosynthetic yields to the chemical abundances in the three most iron-poor stars yet found and reproduce the abundance pattern of one, HE 0557-4840, with a zero metallicity 15 solar mass, 2.4 x 10^51 erg supernova. A Salpeter IMF averaged integration of our yields for Z=0 models with explosion energies of 2.4x10^51 ergs or less is in good agreement with the abundances observed in larger samples of extremely metal-poor stars, provided 15 solar mass stars are included. Since the abundance patterns of extremely metal-poor stars likely arise from a representative sample of progenitors, our yields suggest that low-mass supernovae contributed the bulk of the metals to the early universe.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Bright Supernovae from Magnetar Birth

Following an initial explosion that might be launched either by magnetic interactions or neutrinos, a rotating magnetar radiating according to the classic dipole formula could power a very luminous supernova. While some 56Ni might be produced in the initial explosion, the peak of the light curve in a Type I supernova would not be directly related to its mass. In fact, the peak luminosity would be most sensitive to the dipole field strength of the magnetar. The tail of the light curve could resemble radioactive decay for some time but, assuming complete trapping of the pulsar emission, would eventually be brighter. Depending on the initial explosion energy, both high and moderate velocities could accompany a very luminous light curve.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Low Mach Number Modeling of Type Ia Supernovae. IV. White Dwarf Convection

We present the first three-dimensional, full-star simulations of convection in a white dwarf preceding a Type Ia supernova, specifically the last few hours before ignition. For these long-time calculations we use our low Mach number hydrodynamics code, MAESTRO, which we have further developed to treat spherical stars centered in a three-dimensional Cartesian geometry. The main change required is a procedure to map the one-dimensional radial base state to and from the Cartesian grid. Our models recover the dipole structure of the flow seen in previous calculations, but our long-time integration shows that the orientation of the dipole changes with time. Furthermore, we show the development of gravity waves in the outer, stable portion of the star. Finally, we evolve several calculations to the point of ignition and discuss the range of ignition radii.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The Diversity of Type Ia Supernovae from Broken Symmetries

Type Ia supernovae result when carbon-oxygen white dwarfs in binary systems accrete mass from companion stars, reach a critical mass, and explode. The near uniformity of their light curves makes these supernovae good standard candles for measuring cosmic expansion, but a correction must be applied to account for the fact that the brighter supernovae have broader light curves. One-dimensional modelling, with a certain choice of parameters, can reproduce this general trend in the width-luminosity relation, but the processes of ignition and detonation have recently been shown to be intrinsically asymmetric. Here we report on multi-dimensional modelling of the explosion physics and radiative transfer that reveals that the breaking of spherical symmetry is a critical factor in determining both the width luminosity relation and the observed scatter about it. The deviation from sphericity can also explain the finite polarization detected in the light from some supernovae. The slope and normalization of the width-luminosity relation has a weak dependence on certain properties of the white dwarf progenitor, in particular the trace abundances of elements other than carbon and oxygen. Failing to correct for this effect could lead to systematic overestimates of up to 2% in the distance to remote supernovae.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Type II Supernovae: Model Light Curves and Standard Candle Relationships

A survey of Type II supernovae explosion models has been carried out to determine how their light curves and spectra vary with their mass, metallicity, and explosion energy. The presupernova models are taken from a recent survey of massive stellar evolution at solar metallicity supplemented by new calculations at subsolar metallicity. Explosions are simulated by the motion of a piston near the edge of the iron core and the resulting light curves and spectra are calculated using full multi-wavelength radiation transport. Formulae are developed that describe approximately how the model observables (light curve luminosity and duration) scale with the progenitor mass, explosion energy, and radioactive nucleosynthesis. Comparison with observational data shows that the explosion energy of typical supernovae (as measured by kinetic energy at infinity) varies by nearly an order of magnitude -- from 0.5 to 4.0 x 10^51 ergs, with a typical value of ~0.9 x 10^51 ergs. Despite the large variation, the models exhibit a tight relationship between luminosity and expansion velocity, similar to that previously employed empirically to make SNe IIP standardized candles. This relation is explained by the simple behavior of hydrogen recombination in the supernova envelope, but we find a sensitivity to progenitor metallicity and mass that could lead to systematic errors. Additional correlations between light curve luminosity, duration, and color might enable the use of SNe IIP to obtain distances accurate to ~20% using only photometric data.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Detailed Abundances for 28 Metal-poor Stars: Stellar Relics in the Milky Way

We present the results of an abundance analysis for a sample of stars with $-4<$[Fe/H]$<-2$. The data were obtained with the HIRES spectrograph at Keck Observatory. The set includes 28 stars, with effective temperature ranging from 4800 to 6600 K. For 13 stars with [Fe/H]$<-2.6$, including nine with [Fe/H]$<-3.0$ and one with [Fe/H]$=-4.0$, these are the first reported detailed abundances. For the most metal-poor star in our sample, CS 30336-049, we measure an abundance pattern that is very similar to stars in the range [Fe/H]$\sim-3.5$, including a normal C+N abundance. We also find that it has very low but measurable Sr and Ba, indicating some neutron-capture activity even at this low of a metallicity. We explore this issue further by examining other very neutron-capture-deficient stars, and find that at the lowest levels, [Ba/Sr] exhibits the ratio of the main r-process. We also report on a new r-process-enhanced star, CS 31078-018. This star has [Fe/H]$=-2.85$, [Eu/Fe]$=1.23$, and [Ba/Eu]$=-0.51$. CS 31078-018 exhibits an ``actinide boost'', i.e. much higher [Th/Eu] than expected and at a similar level to CS 31082-001. Our spectra allow us to further constrain the abundance scatter at low metallicities, which we then use to fit to the zero-metallicity Type II supernova yields of Heger & Woosley (2008). We find that supernovae with progenitor masses between 10 and 20 M$_{\odot}$ provide the best matches to our abundances.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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)

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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'' 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}.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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. "Gravitationally confined detonation" 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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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'' in the explosion. A grid of final compositions is explored including essentially all combinations of 56Ni, stable ``iron'', 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'' 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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

Massive Star Evolution: Nucleosynthesis and Nuclear Reaction Rate Uncertainties

We present a nucleosynthesis calculation of a 25 solar mass star of solar composition that includes all relevant isotopes up to polonium. In particular, all stable isotopes and necessary nuclear reaction rates are covered. We follow the stellar evolution from hydrogen burning till iron core collapse and simulate the explosion using a ``piston'' approach. We discuss the influence of two key nuclear reaction rates, C12(a,g) and Ne22(a,n), on stellar evolution and nucleosynthesis. The former significantly influences the resulting core sizes (iron, silicon, oxygen) and the overall presupernova structure of the star. It thus has significant consequences for the supernova explosion itself and the compact remnant formed. The later rate considerably affects the s-process in massive stars and we demonstrate the changes that different currently suggested values for this rate cause.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint1998arXiv

Gamma-Ray Bursts From Neutron Star Phase Transitions

The phase-transition induced collapse of a neutron star to a more compact configuration (typically a ``strange'' star) and the subsequent core bounce is often invoked as a model for gamma-ray bursts. We present the results of numerical simulations of this kind of event using realistic neutrino physics and a high density equation of state. The nature of the collapse itself is represented by the arbitrary motion of a piston deep within the star, but if any shock is to develop, the transition, or at least its final stages, must occur in less than a sonic time. Fine surface zoning is employed to adequately represent the acceleration of the shock to relativistic speeds and to determine the amount and energy of the ejecta. We find that these explosions are far too baryon-rich (ejected Mass > 0.01 solar masses) and have much too low an energy to explain gamma-ray bursts. The total energy of the ejecta having relativistic lorentz factors > 40 is less than 10^46 erg even in our most optimistic models (deep bounce, no neutrino losses or photodisintegration). However, the total energy of all the ejecta, mostly mildly relativistic, is roughly 10^51 erg and, if they occur, these events might be observed. They would also contribute to Galactic nucleosynthesis, especially the r-process, even though the most energetic layers are composed of helium and nucleons, not heavy elements.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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's and the lower mass group of the Type II'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

0 citations0 reviewsNo code linkedOpen access