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.

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.

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

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

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

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

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

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

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

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

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

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

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

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

The density structure surrounding the iron core of a massive star when it dies is known to have a major effect on whether or not the star explodes. Here we repeat previous surveys of presupernova evolution with some important corrections to code physics and four to ten times better mass resolution in each star. The number of presupernova masses considered is also much larger. Over 4,000 models are calculated in the range from 12 to 60 Msun with varying mass loss rates. The core structure is not greatly affected by the increased spatial resolution. The qualitative patterns of compactness measures and their extrema are the same, but with the increased number of models, the scatter seen in previous studies is replaced by several localized branches. More physics-based analyses by Ertl et al. (2016) and M$\"u$eller et al. (2016) show these branches with less scatter than the single parameter characterization of O'Connor & Ott (2011). These branches are particularly apparent for stars in the mass ranges 14 - 19 Msun and 22 - 24 Msun. The multi-valued solutions are a consequence of interference between several carbon and oxygen burning shells during the late stages of evolution. For a relevant range of masses, whether a star explodes or not may reflect more the small, almost random differences in its late evolution than its initial mass. The large number of models allow statistically meaningful statements about the radius, luminosity, and effective temperatures of presupernova stars, their core structures, and their remnant mass distribution.

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

A large fraction of core-collapse supernovae (CCSNe), 30-50%, are expected to originate from the low-mass end of progenitors with $M_{\rm ZAMS}~= 8-12~M_\odot$. However, degeneracy effects make stellar evolution modelling of such stars challenging, and few predictions for their supernova light curves and spectra have been presented. Here we calculate synthetic nebular spectra of a 9 $M_\odot$ Fe CCSN model exploded with the neutrino mechanism. The model predicts emission lines with FWHM$\sim$1000 km/s, including signatures from each deep layer in the metal core. We compare this model to observations of the three subluminous IIP SNe with published nebular spectra; SN 1997D, SN 2005cs, and SN 2008bk. The prediction of both line profiles and luminosities are in good agreement with SN 1997D and SN 2008bk. The close fit of a model with no tuning parameters provides strong evidence for an association of these objects with low-mass Fe CCSNe. For SN 2005cs, the interpretation is less clear, as the observational coverage ended before key diagnostic lines from the core had emerged. We perform a parameterised study of the amount of explosively made stable nickel, and find that none of these three SNe show the high $^{58}$Ni/$^{56}$Ni ratio predicted by current models of electron capture SNe (ECSNe) and ECSN-like explosions. Combined with clear detection of lines from O and He shell material, these SNe rather originate from Fe core progenitors. We argue that the outcome of self-consistent explosion simulations of low-mass stars, which gives fits to many key observables, strongly suggests that the class of subluminous Type IIP SNe is the observational counterpart of the lowest mass CCSNe.

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

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

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

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

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

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

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.

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

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.

Recent observations have revealed an amazing diversity of extremely luminous supernovae, seemingly increasing in radiant energy without bound. We consider here the physical limits of what existing models can provide for the peak luminosity and total radiated energy for non-relativistic, isotropic stellar explosions. The brightest possible supernova is a Type I explosion powered by a sub-millisecond magnetar. Such models can reach a peak luminosity of $\rm 2\times10^{46}\ erg\ s^{-1}$ and radiate a total energy of $\rm 4 \times10^{52}\ erg$. Other less luminous models are also explored, including prompt hyper-energetic explosions in red supergiants, pulsational-pair instability supernovae, and pair-instability supernovae. Approximate analytic expressions and limits are given for each case. Excluding magnetars, the peak luminosity is near $\rm 1\times10^{44}\ erg\ s^{-1}$ for the brightest models. The corresponding limits on total radiated power are $\rm3 \times 10^{51}\ erg$ (Type I) and $\rm1 \times 10^{51}\ erg$ (Type II). A magnetar-based model for the recent transient event, ASASSN-15lh is presented that strains, but does not exceed the limits of what the model can provide.

We present 1-D non-Local-Thermodynamic-Equilibrium time-dependent radiative-transfer simulations for a large grid of supernovae (SNe) IIb/Ib/Ic that result from the terminal explosion of the mass donor in a close-binary system. Our sample covers ejecta masses $M_{\rm e}$ of 1.7$-$5.2M$_\odot$, kinetic energies $E_{\rm kin}$ of 0.6$-$5.0$\times$10$^{51}$erg, and $^{56}$Ni masses of 0.05$-$0.30M$_\odot$. We find a strong correlation between the $^{56}$Ni mass and the photometric properties at maximum, and between the rise time to bolometric maximum and the post-maximum decline rate. We confirm the small scatter in ($V-R$) at 10d past $R$-band maximum. The quantity $V_{\rm m} \equiv \sqrt{2E_{\rm kin}/M_{\rm e}}$ is comparable to the Doppler velocity measured from HeI 5875\AA\ at maximum in SNe IIb/Ib, although some scatter arises from the uncertain level of chemical mixing. The OI7772\AA\ line may be used for SNe Ic, but the correspondence deteriorates with higher ejecta mass/energy. We identify a temporal reversal of the Doppler velocity at maximum absorption in the $\sim$1.05$\mu$m feature in all models. The reversal is due to HeI alone and could serve as a test for the presence of helium in SNe Ic. Because of variations in composition and ionisation, the ejecta opacity shows substantial variations with both velocity and time. This is in part the origin of the offset between our model light curves and the predictions from the Arnett model.

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.

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.

We present 1-D non-Local-Thermodynamic-Equilibrium time-dependent radiative-transfer simulations for supernovae (SNe) of type IIb, Ib, and Ic that result from the terminal explosion of the mass donor in a close-binary system. Here, we select three ejecta with a total kinetic energy of ~1.2e51erg, but characterised by different ejecta masses (2-5Msun), composition, and chemical mixing. The type IIb/Ib models correspond to the progenitors that have retained their He-rich shell at the time of explosion. The type Ic model arises from a progenitor that has lost its helium shell, but retains 0.32Msun of helium in a CO-rich core of 5.11Msun. We discuss their photometric and spectroscopic properties during the first 2-3 months after explosion, and connect these to their progenitor and ejecta properties including chemical stratification. For these three models, Arnett's rule overestimates the 56Ni mass by ~50% while the procedure of Katz et al., based on an energy argument, yields a more reliable estimate. The presence of strong CI lines around 9000A prior to maximum is an indicator that the pre-SN star was under-abundant in helium. As noted by others, the 1.08micron feature is a complex blend of CI, MgII, and HeI lines, which makes the identification of He uncertain in SNe Ibc unless other HeI lines can be identified. Our models show little scatter in (V-R) colour 10d after R-band maximum. We also address a number of radiative transfer properties of SNe Ibc, including the notion of a photosphere, the inference of a representative ejecta expansion rate, spectrum formation, blackbody fits and "correction factors".

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.

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.

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

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

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

This document presents the results from the Distances subgroup of the Cosmic Frontier Community Planning Study (Snowmass 2013). We summarize the current state of the field as well as future prospects and challenges. In addition to the established probes using Type IA supernovae and baryon acoustic oscillations, we also consider prospective methods based on clusters, active galactic nuclei, gravitational wave sirens and strong lensing time delays.

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

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

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

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.

The kinetic energy of supernovae (SNe) accompanied by gamma-ray bursts (GRBs) tends to cluster near E52 erg, with 2.E52 erg an upper limit to which no compelling exceptions are found (assuming a certain degree of asphericity), and it is always significantly larger than the intrinsic energy of the GRB themselves (corrected for jet collimation). This energy is strikingly similar to the maximum rotational energy of a neutron star rotating with period 1 ms. It is therefore proposed that all GRBs associated with luminous SNe are produced by magnetars. GRBs that result from black hole formation (collapsars) may not produce luminous SNe. X-ray Flashes (XRFs), which are associated with less energetic SNe, are produced by neutron stars with weaker magnetic field or lower spin.

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.

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.

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

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

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

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

Using 2D and 3D simulation, we study the "robustness" of the double detonation scenario for Type Ia supernovae, in which a detonation in the helium shell of a carbon-oxygen white dwarf induces a secondary detonation in the underlying core. We find that a helium detonation cannot easily descend into the core unless it commences (artificially) well above the hottest layer calculated for the helium shell in current presupernova models. Compressional waves induced by the sliding helium detonation, however, robustly generate hot spots which trigger a detonation in the core. Our simulations show that this is true even for non-axisymmetric initial conditions. If the helium is ignited at multiple points, the internal waves can pass through one another or be reflected, but this added complexity does not defeat the generation of the hot spot. The ignition of very low-mass helium shells depends on whether a thermonuclear runaway can simultaneously commence in a sufficiently large region.

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.

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

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.

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.

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

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

We model neutrino emission from a newly born neutron star subsequent to a supernova explosion to study its sensitivity to the equation of state, neutrino opacities, and convective instabilities at high baryon density. We find the time period and spatial extent over which convection operates is sensitive to the behavior of the nuclear symmetry energy at and above nuclear density. When convection ends within the proto-neutron star, there is a break in the predicted neutrino emission that may be clearly observable.

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.

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.

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.

We develop a one-dimensional theoretical model for thermals burning in Type Ia supernovae based on the entrainment assumption of Morton, Taylor and Turner. Extensions of the standard model are required to account for the burning and for the expansion of the thermal due to changes in the background stratification found in the full star. The model is compared with high-resolution three-dimensional numerical simulations, both in a uniform environment, and in a full-star setting. The simulations in a uniform environment present compelling agreement with the predicted power-laws and provide model constants for the full-star model, which then provides excellent agreement with the full-star simulation. The importance of the different components in the model are compared, and are all shown to be relevant. An examination of the effect of initial conditions was then conducted using the one-dimensional model, which would have been infeasible in three dimensions. More mass was burned when the ignition kernel was larger and closer to the center of the star. The turbulent flame speed was found to be important during the early-time evolution of the thermal, but played a diminished role at later times when the evolution is dominated by the large-scale hydrodynamics responsible for entrainment. However, a higher flame speed effectively gave a larger initial ignition kernel and so resulted in more mass burned. This suggests that future studies should focus on the early-time behavior of these thermals (in particular, the transition to turbulence), and that the choice of turbulent flame speed does not play a significant role in the dynamics once the thermal has become established.

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.

In the distributed burning regime, turbulence disrupts the internal structure of the flame, and so the idea of laminar burning propagated by conduction is no longer valid. The nature of the burning depends on the turbulent Damkohler number (Da), which steadily declines from much greater than one to less that one as the density decreases to a few 10^6 g/cc. Scaling arguments predict that the turbulent flame speed s, normalized by the turbulent intensity u, follows s/u=Da^1/2 for Da<1. The flame in this regime is a single turbulently-broadened structure that moves at a steady speed, and has a width larger than the integral scale of the turbulence. The scaling is predicted to break down at Da=1, and the flame burns as a turbulently-broadened effective unity Lewis number flame. We refer to this kind of flame as a lambda-flame. The burning becomes a collection of lambda-flames spread over a region approximately the size of the integral scale. While the total burning rate continues to have a well-defined average, s_{T} ~ u, the burning is unsteady. We present a theoretical framework, supported by both 1D and 3D numerical simulations, for the burning in these two regimes. Our results indicate that the average value of s can actually be roughly twice u for Da>1, and that localized excursions to as much as five times u can occur. The lambda-flame speed and width can be predicted based on the turbulence in the star and the turbulent nuclear burning time scale of the fuel. We propose a practical method for measuring these based on the scaling relations and small-scale computationally-inexpensive simulations. This suggests that a simple turbulent flame model can be easily constructed suitable for large-scale distributed supernovae flames.

In previous studies, we examined turbulence-flame interactions in carbon-burning thermonuclear flames in Type Ia supernovae. In this study, we consider turbulence-flame interactions in the trailing oxygen flames. The two aims of the paper are to examine the response of the inductive oxygen flame to intense levels of turbulence, and to explore the possibility of transition to detonation in the oxygen flame. Scaling arguments analogous to the carbon flames are presented and then compared against three-dimensional simulations for a range of Damk\"ohler numbers ($\Da_{16}$) at a fixed Karlovitz number. The simulations suggest that turbulence does not significantly affect the oxygen flame when $\Da_{16}<1$, and the flame burns inductively some distance behind the carbon flame. However, for $\Da_{16}>1$, turbulence enhances heat transfer and drives the propagation of a flame that is {\em narrower} than the corresponding inductive flame would be. Furthermore, burning under these conditions appears to occur as part of a combined carbon-oxygen turbulent flame with complex compound structure. The simulations do not appear to support the possibility of a transition to detonation in the oxygen flame, but do not preclude it either.

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.

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.

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

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.

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.

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.

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.

We investigate the evolution of Type Ib/c supernova (SN Ib/c) progenitors in close binary systems, using new evolutionary models that include the effects of rotation, with initial masses of 12 - 25 Msun for the primary components, and of single helium stars with initial masses of 2.8 - 20 Msun. We find that, despite the impact of tidal interaction on the rotation of primary stars, the amount of angular momentum retained in the core at the presupernova stage in different binary model sequences converge to a value similar to those found in previous single star models. This amount is large enough to produce millisecond pulsars, but too small to produce magnetars or long gamma-ray bursts. We employ the most up-to-date estimate for the Wolf-Rayet mass loss rate, and its implications for SN Ib/c progenitors are discussed in detail. In terms of stellar structure, SN Ib/c progenitors in binary systems are predicted to have a wide range of final masses even up to 7 Msun, with helium envelopes of 0.16 - 1.5 Msun. Our results indicate that, if the lack of helium lines in the spectra of SNe Ic were due to small amounts of helium, the distribution of both initial and final masses of SN Ic progenitors should be bimodal. Furthermore, we find that a thin hydrogen layer (0.001 - 0.01 Msun) is expected to be present in many SN Ib progenitors at the presupernova stage. We show that the presence of hydrogen, together with a rather thick helium envelope, can lead to a significant expansion of some SN Ib/c progenitors by the time of supernova explosion. This may have important consequences for the shock break-out and supernova light curve. We also argue that some SN progenitors with thin hydrogen layers produced via Case AB/B transfer might be related to Type IIb supernova progenitors with relatively small radii of about 10 Rsun.

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.

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.

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.

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.

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.

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.

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.

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.

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.

The nature of carbon burning flames in Type Ia supernovae is explored as they interact with Kolmogorov turbulence. One-dimensional calculations using the Linear Eddy Model of Kerstein (1991) elucidate three regimes of turbulent burning. In the simplest case, large scale turbulence folds and deforms thin laminar flamelets to produce a flame brush with a total burning rate given approximately by the speed of turbulent fluctuations on the integral scale, U_L. This is the regime where the supernova explosion begins and where most of its pre-detonation burning occurs. As the density declines, turbulence starts to tear the individual flamelets, making broader structures that move faster. For a brief time, these turbulent flamelets are still narrow compared to their spacing and the concept of a flame brush moving with an overall speed of U_L remains valid. However, the typical width of the individual flamelets, which is given by the condition that their turnover time equals their burning time, continues to increase as the density declines. Eventually, mixed regions almost as large as the integral scale itself are transiently formed. At that point, a transition to detonation can occur. The conditions for such a transition are explored numerically and it is estimated that the transition will occur for densities near 1 x 10**7 g/cm**3, provided the turbulent speed on the integral scale exceeds about 15% sonic. An example calculation shows the details of a detonation actually developing.

Aims: Taking advantage of the forthcoming Catalog of the HETE-2 mission, the aim of this paper is to evaluate the main properties of HETE-2 GRBs - the E_peak, the T_90 and the E_iso - in their source frames and to derive their unbiased distribution. Methods: We first construct a complete sample containing all the bursts localized by the WXM on-board HETE-2, which are selected with a uniform criterion and whose observed parameters can be constrained. We then derive the intrinsic E_peak, T_90 and E_iso distributions using their redshift when it is available, or their pseudo-redshift otherwise. We finally compute the number of GRB (N_Vmax) within the visibility volume (V_max) of each GRB, in order to derive a weight for each detected burst accounting both for the detection significance and the star formation history of the universe. Results: The unbiased distributions obtained clearly show the predominence of X-ray flashes (XRFs) in the global GRB population. We also derive the rate of local GRBs: R0^H2 > 11 Gpc-3 yr-1, which is intermediate between the local rate obtained by considering only the high-luminosity bursts (~1 Gpc-3 yr-1) and that obtained by including the low-luminosity bursts (>200 Gpc-3 yr-1).

The large range of time and length scales involved in type Ia supernovae (SN Ia) requires the use of flame models. As a prelude to exploring various options for flame models, we consider, in this paper, high-resolution three-dimensional simulations of the small-scale dynamics of nuclear flames in the supernova environment in which the details of the flame structure are fully resolved. The range of densities examined, 1 to $8 \times 10^7$ g cm$^{-3}$, spans the transition from the laminar flamelet regime to the distributed burning regime where small scale turbulence disrupts the flame. The use of a low Mach number algorithm facilitates the accurate resolution of the thermal structure of the flame and the inviscid turbulent kinetic energy cascade, while implicitly incorporating kinetic energy dissipation at the grid-scale cutoff. For an assumed background of isotropic Kolmogorov turbulence with an energy characteristic of SN Ia, we find a transition density between 1 and $3 \times 10^7$ g cm$^{-3}$ where the nature of the burning changes qualitatively. By $1 \times 10^7$ g cm$^{-3}$, energy diffusion by conduction and radiation is exceeded, on the flame scale, by turbulent advection. As a result, the effective Lewis Number approaches unity. That is, the flame resembles a laminar flame, but is turbulently broadened with an effective diffusion coefficient, $D_T \sim u' l$, where $u'$ is the turbulent intensity and $l$ is the integral scale. For the larger integral scales characteristic of a real supernova, the flame structure is predicted to become complex and unsteady. Implications for a possible transition to detonation are discussed.

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.

We report the serendipitous detection by GALEX of fast (<1 day) rising (>1 mag) UV emission from two Type II plateau (II-P) supernovae (SNe) at z=0.185 and 0.324 discovered by the Supernova Legacy Survey. Optical photometry and VLT spectroscopy 2 weeks after the GALEX detections link the onset of UV emission to the time of shock breakout. Using radiation hydrodynamics and non-LTE radiative transfer simulations, and starting from a standard red supergiant (RSG; Type II-P SN progenitor) star evolved self-consistently from the main sequence to iron core collapse, we model the shock breakout phase and the 55 hr that follow. The small scale height of our RSG atmosphere model suggests that the breakout signature is a thermal soft X-ray burst (lambda_peak ~ 90\AA) with a duration of <~ 2000 s. Longer durations are possible but require either an extended and tenuous non-standard envelope, or an unusually dense RSG wind with \dot{M} ~ 10^(-3) Msun yr^(-1). The GALEX observations miss the peak of the luminous (M_FUV ~ -20) UV burst but unambiguously capture the rise of the emission and a subsequent 2 day long plateau. The postbreakout, UV-bright plateau is a prediction of our model in which the shift of the peak of the spectral energy distribution (SED) from ~100 to ~1000\AA and the ejecta expansion both counteract the decrease in bolometric luminosity from ~10^11 to ~10^9 L_sun over that period. Based on the observed detection efficiency of our study we make predictions for the breakout detection rate of the GALEX Time Domain Survey.

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.

GRB 041006 was detected by HETE-2 at 12:18:08 UT on 06 October 2004. This GRB displays a soft X-ray emission, a precursor before the onset of the main event, and also a soft X-ray tail after the end of the main peak. The light curves in four different energy bands display different features; At higher energy bands several peaks are seen in the light curve, while at lower energy bands a single broader bump dominates. It is expected that these different features are the result of a mixture of several components each of which has different energetics and variability. To reveal the nature of each component, we analysed the time resolved spectra and they are successfully resolved into several components. We also found that these components can be classified into two distinct classes; One is a component which has an exponential decay of $E_{p}$ with a characteristic timescale shorter than $\sim$ 30 sec, and its spectrum is well represented by a broken power law function, which is frequently observed in many prompt GRB emissions, so it should have an internal-shock origin. Another is a component whose $E_{p}$ is almost unchanged with characteristic timescale longer than $\sim$ 60 sec, and shows a very soft emission and slower variability. The spectrum of the soft component is characterized by either a broken power law or a black body spectrum. This component might originate from a relatively wider and lower velocity jet or a photosphere of the fireball. By assuming that the soft component is a thermal emission, the radiation radius is initially $4.4 \times 10^{6}$ km, which is a typical radius of a blue supergiant, and its expansion velocity is $2.4 \times 10^{5}$ km/s in the source frame.

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)

We discuss how rotation and binary interactions may be related to the diversity of type Ibc supernovae and long gamma-ray bursts. After presenting recent evolutionary models of massive single and binary stars including rotation, the Tayler-Spruit dynamo and binary interactions, we argue that the nature of SNe Ibc progenitors from binary systems may not significantly differ from that of single star progenitors in terms of rotation, and that most long GRB progenitors may be produced via the quasi-chemically homogeneous evolution at sub-solar metallicity. We also briefly discuss the possible role of magnetic fields generated in the convective core of a massive star for the transport of angular momentum, which is potentially important for future stellar evolution models of supernova and GRB progenitors.

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.

Recent nucleosynthesis calculations of Type II supernovae using advanced neutrino transport determine that the early neutrino winds are proton-rich. However, a fraction of the ejecta emitted at the same time is composed of neutron-rich pockets. In this paper we calculate the nucleosynthesis contribution from the neutron-rich pockets in the hot convective bubbles of a core-collapse supernova and show that they do not contribute significantly to the total nucleosynthesis.

The compact remnants of core collapse supernovae - neutron stars and black holes - have properties that reflect both the structure of their stellar progenitors and the physics of the explosion. In particular, the masses of these remnants are sensitive to the density structure of the presupernova star and to the explosion energy. To a considerable extent, the final mass is determined by the ``fallback'', during the explosion, of matter that initially moves outwards, yet ultimately fails to escape. We consider here the simulated explosion of a large number of massive stars (10 to 100 \Msun) of Population I (solar metallicity) and III (zero metallicity), and find systematic differences in the remnant mass distributions. As pointed out by Chevalier(1989), supernovae in more compact progenitor stars have stronger reverse shocks and experience more fallback. For Population III stars above about 25 \Msun and explosion energies less than $1.5 \times 10^{51}$ erg, black holes are a common outcome, with masses that increase monotonically with increasing main sequence mass up to a maximum hole mass of about 35 \Msun. If such stars produce primary nitrogen, however, their black holes are systematically smaller. For modern supernovae with nearly solar metallicity, black hole production is much less frequent and the typical masses, which depend sensitively on explosion energy, are smaller. We explore the neutron star initial mass function for both populations and, for reasonable assumptions about the initial mass cut of the explosion, find good agreement with the average of observed masses of neutron stars in binaries. We also find evidence for a bimodal distribution of neutron star masses with a spike around 1.2 \Msun (gravitational mass) and a broader distribution peaked around 1.4 \Msun.

The X-ray burster GS 1826-24 shows extremely regular Type I X-ray bursts whose energetics and recurrence times agree well with thermonuclear ignition models. We present calculations of sequences of burst lightcurves using multizone models which follow the rp-process nucleosynthesis with an extensive nuclear reaction network. The theoretical and observed burst lightcurves show remarkable agreement. The models naturally explain the slow ~5s rise and long ~100s tails of these bursts, as well as their dependence on mass accretion rate. This comparison provides further evidence for solar metallicity in the accreted material in this source, and constrains the distance to the source. The main difference is that the observed lightcurves do not show the distinct two-stage rise of the models. This may reflect the time for burning to spread over the stellar surface, or may indicate that our treatment of heat transport or nuclear physics needs to be revised. The trends in burst properties with accretion rate are well-reproduced by our spherically symmetric models which include chemical and thermal inertia from the ashes of previous bursts. Changes in the covering fraction of the accreted fuel are not required.

The extremely luminous supernova SN 2006gy challenges the traditional view that the collapse of a stellar core is the only mechanism by which a massive star makes a supernova, because it seems too luminous by more than a factor of ten. Here we report that the brightest supernovae in the modern Universe arise from collisions between shells of matter ejected by massive stars that undergo an interior instability arising from the production of electron-positron pairs. This "pair instability" leads to explosive burning that is insufficient to unbind the star, but ejects many solar masses of the envelope. After the first explosion, the remaining core contracts and searches for a stable burning state. When the next explosion occurs, several solar masses of material are again ejected, which collide with the earlier ejecta. This collision can radiate 1E50 erg of light, about a factor of ten more than an ordinary supernova. Our model is in good agreement with the observed light curve for SN 2006gy and also shows that some massive stars can produce more than one supernova-like outburst.

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

A long X-ray flash was detected and localized by the instruments aboard the High Energy Transient Explorer II (HETE-2) at 00:03:30 UT on 2004 September 16. The position was reported to the GRB Coordinates Network (GCN) approximately 2 hours after the burst. This burst consists of two peaks separated by 200 s, with durations of 110 s and 60 s. We have analyzed the energy spectra of the 1st and 2nd peaks observed with the Wide Field X-Ray Monitor (WXM) and the French Gamma Telescope (FREGATE). We discuss the origin of the 2nd peak in terms of flux variabilities and timescales. We find that it is most likely part of the prompt emission, and is explained by the long-acting engine model. This feature is similar to some bright X-ray flares detected in the early afterglow phase of bursts observed by the Swift satellite.

Hans Bethe contributed in many ways to our understanding of the supernovae that happen in massive stars, but, to this day, a first principles model of how the explosion is energized is lacking. Nevertheless, a quantitative theory of nucleosynthesis is possible. We present a survey of the nucleosynthesis that occurs in 32 stars of solar metallicity in the mass range 12 to 120 solar masses. The most recent set of solar abundances, opacities, mass loss rates, and current estimates of nuclear reaction rates are employed. Restrictions on the mass cut and explosion energy of the supernovae based upon nucleosynthesis, measured neutron star masses, and light curves are discussed and applied. The nucleosynthetic results, when integrated over a Salpeter initial mass function (IMF), agree quite well with what is seen in the sun. We discuss in some detail the production of the long lived radioactivities, 26Al and 60Fe, and why recent model-based estimates of the ratio 60Fe/26Al are overly large compared with what satellites have observed. A major source of the discrepancy is the uncertain nuclear cross sections for the creation and destruction of these unstable isotopes.

The observational diversity of ``gamma-ray bursts'' (GRBs) has been increasing, and the natural inclination is a proliferation of models. We explore the possibility that at least part of this diversity is a consequence of a single basic model for the central engine operating in a massive star of variable mass, differential rotation rate, and mass loss rate. Whatever that central engine may be - and here the collapsar is used as a reference point - it must be capable of generating both a narrowly collimated, highly relativistic jet to make the GRB, and a wide angle, sub-relativistic outflow responsible for exploding the star and making the supernova bright. To some extent, the two components may vary independently, so it is possible to produce a variety of jet energies and supernova luminosities. We explore, in particular, the production of low energy bursts and find a lower limit, $\sim10^{48}$ erg s$^{-1}$ to the power required for a jet to escape a massive star before that star either explodes or is accreted. Lower energy bursts and and ``suffocated'' bursts may be particularly prevalent when the metallicity is high, i.e., in the modern universe at low redshift.

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.

In this work, we present a multi-wavelength study of XRF 040912, aimed at measuring its distance scale and the intrinsic burst properties. We performed a detailed spectral and temporal analysis of both the prompt and the afterglow emission and we estimated the distance scale of the likely host galaxy. We then used the currently available sample of XRFs with known distance to discuss the connection between XRFs and classical Gamma-ray Bursts (GRBs). We found that the prompt emission properties unambiguously identify this burst as an XRF, with an observed peak energy of E_p=17+/-13 keV and a burst fluence ratio S(2-30keV)/S(30-400keV)>1. A non-fading optical source with R~24 mag and with an apparently extended morphology is spatially consistent with the X-ray afterglow, likely the host galaxy. XRF 040912 is a very dark burst since no afterglow optical counterpart is detected down to R>25 mag (3 sigma limiting magnitude) at 13.6 hours after the burst. The host galaxy spectrum detected from 3800A to 10000A, shows a single emission line at 9552A. The lack of any other strong emission lines blue-ward of the detected one and the absence of the Ly alpha cut-off down to 3800A are consistent with the hypothesis of the [OII] line at redshift z=1.563+/-0.001. The intrinsic spectral properties rank this XRF among the soft GRBs in the E_peak-E_iso diagram. Similar results were obtained for most XRFs at known redshift. Only XRF 060218 and XRF 020903 represent a good example of instrinsic XRF(i-XRF) and are possibly associated with a different progenitor population. This scenario may calls for a new definition of XRFs.

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.

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.

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.

Brighter Type Ia supernovae (SNe Ia) have broader, more slowly declining B-band light curves than dimmer SNe Ia. We study the physical origin of this width-luminosity relation (WLR) using detailed radiative transfer calculations of Chandrasekhar mass SN Ia models. We find that the luminosity dependence of the diffusion time (emphasized in previous studies) is in fact of secondary relevance in understanding the model WLR. Instead, the essential physics involves the luminosity dependence of the spectroscopic/color evolution of SNe Ia. Following maximum-light, the SN colors are increasingly affected by the development of numerous Fe II/Co II lines which blanket the B-band and, at the same time, increase the emissivity at longer wavelengths. Because dimmer SNe Ia are generally cooler, they experience an earlier onset of Fe III to Fe II recombination in the iron-group rich layers of ejecta, resulting in a more rapid evolution of the SN colors to the red. The faster B-band decline rate of dimmer SNe Ia thus reflects their faster ionization evolution.

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.

Observations show that at least some gamma-ray bursts (GRBs) happen simultaneously with core-collapse supernovae (SNe), thus linking by a common thread nature's two grandest explosions. We review here the growing evidence for and theoretical implications of this association, and conclude that most long-duration soft-spectrum GRBs are accompanied by massive stellar explosions (GRB-SNe). The kinetic energy and luminosity of well-studied GRB-SNe appear to be greater than those of ordinary SNe, but evidence exists, even in a limited sample, for considerable diversity. The existing sample also suggests that most of the energy in the explosion is contained in nonrelativistic ejecta (producing the supernova) rather than in the relativistic jets responsible for making the burst and its afterglow. Neither all SNe, nor even all SNe of Type Ibc produce GRBs. The degree of differential rotation in the collapsing iron core of massive stars when they die may be what makes the difference.

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.

Long-duration gamma-ray bursts (GRBs) are associated with type Ic supernovae that are more luminous than average and that eject material at very high velocities. Less-luminous supernovae were not hitherto known to be associated with GRBs, and therefore GRB-supernovae were thought to be rare events. Whether X-ray flashes - analogues of GRBs, but with lower luminosities and fewer gamma-rays - can also be associated with supernovae, and whether they are intrinsically 'weak' events or typical GRBs viewed off the axis of the burst, is unclear. Here we report the optical discovery and follow-up observations of the type Ic supernova SN 2006aj associated with X-ray flash XRF 060218. Supernova 2006aj is intrinsically less luminous than the GRB-supernovae, but more luminous than many supernovae not accompanied by a GRB. The ejecta velocities derived from our spectra are intermediate between these two groups, which is consistent with the weakness of both the GRB output and the supernova radio flux. Our data, combined with radio and X-ray observations, suggest that XRF 060218 is an intrinsically weak and soft event, rather than a classical GRB observed off-axis. This extends the GRB-supernova connection to X-ray flashes and fainter supernovae, implying a common origin. Events such as XRF 060218 are probably more numerous than GRB-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.

GRB 051022 was detected at 13:07:58 on 22 October 2005 by HETE-2. The location of GRB 051022 was determined immediately by the flight localization system. This burst contains multiple pulses and has a rather long duration of about 190 seconds. The detections of candidate X-ray and radio afterglows were reported, whereas no optical afterglow was found. The optical spectroscopic observations of the host galaxy revealed the redshift z = 0.8. Using the data derived by HETE-2 observation of the prompt emission, we found the absorption N_H = 8.8 -2.9/+3.1 x 10^22 cm^-2 and the visual extinction A_V = 49 -16/+17 mag in the host galaxy. If this is the case, no detection of any optical transient would be quite reasonable. The absorption derived by the Swift XRT observations of the afterglow is fully consistent with those obtained from the early HETE-2 observation of the prompt emission. Our analysis implies an interpretation that the absorbing medium could be outside external shock at R ~ 10^16 cm, which may be a dusty molecular cloud.

Type Ia supernovae (SNe Ia) are the largest thermonuclear explosions in the Universe. Their light output can be seen across great distances and has led to the discovery that the expansion rate of the Universe is accelerating. Despite the significance of SNe Ia, there are still a large number of uncertainties in current theoretical models. Computational modeling offers the promise to help answer the outstanding questions. However, even with today's supercomputers, such calculations are extremely challenging because of the wide range of length and time scales. In this paper, we discuss several new algorithms for simulations of SNe Ia and demonstrate some of their successes.

Here we report the localizations and properties of four short-duration GRBs localized by the High Energy Transient Explorer 2 satellite (HETE-2): GRBs 010326B, 040802, 051211 and 060121, all of which were detected by the French Gamma Telescope (Fregate) and localized with the Wide-field X-ray Monitor (WXM) and/or Soft X-ray Camera (SXC) instruments. We discuss eight possible criteria for determining whether these GRBs are "short population bursts" (SPBs) or "long population bursts" (LPBs). These criteria are (1) duration, (2) pulse widths, (3) spectral hardness, (4) spectral lag, (5) energy Egamma radiated in gamma rays (or equivalently, the kinetic energy E_KE of the GRB jet), (6) existence of a long, soft bump following the burst, (7) location of the burst in the host galaxy, and (8) type of host galaxy. In particular, we have developed a likelihood method for determining the probability that a burst is an SPB or a LPB on the basis of its T90 duration alone. A striking feature of the resulting probability distribution is that the T90 duration at which a burst has an equal probability of being a SPB or a LPB is T90 = 5 s, not T90 = 2 s, as is often used. All four short-duration bursts discussed in detail in this paper have T90 durations in the Fregate 30-400 keV energy band of 1.90, 2.31, 4.25, and 1.97 sec, respectively, yielding probabilities P(S|T90) = 0.97, 0.91, 0.60, and 0.95 that these bursts are SPBs on the basis of their T90 durations alone. All four bursts also have spectral lags consistent with zero. These results provide strong evidence that all four GRBs are SPBs (abstract continues).

When massive stars exhaust their fuel they collapse and often produce the extraordinarily bright explosions known as core-collapse supernovae. On occasion, this stellar collapse also powers an even more brilliant relativistic explosion known as a long-duration gamma-ray burst. One would then expect that long gamma-ray bursts and core-collapse supernovae should be found in similar galactic environments. Here we show that this expectation is wrong. We find that the long gamma-ray bursts are far more concentrated on the very brightest regions of their host galaxies than are the core-collapse supernovae. Furthermore, the host galaxies of the long gamma-ray bursts are significantly fainter and more irregular than the hosts of the core-collapse supernovae. Together these results suggest that long-duration gamma-ray bursts are associated with the most massive stars and may be restricted to galaxies of limited chemical evolution. Our results directly imply that long gamma-ray bursts are relatively rare in galaxies such as our own Milky Way.

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

We present ground-based and Hubble Space Telescope optical observations of the X-ray flash (XRF) 020903, covering 300 days. The afterglow showed a very rapid rise in the first day, followed by a relatively slow decay in the next few days. There was a clear bump in the light curve after ~25 days, accompanied by a drastic change in the spectral energy distribution. The light curve and the spectral energy distribution are naturally interpreted as the emergence -- and subsequent decay -- of a supernova (SN), similar to SN 1998bw. At peak luminosity, the SN is estimated to be 0.8 +/- 0.1 mag fainter than SN1998bw. This argues in favor of the existence of a supernova associated with this X-ray flash. A spectrum obtained 35 days after the burst shows emission lines from the host galaxy. We use this spectrum to put an upper limit on the oxygen abundance of the host at [O/H] < -0.6 dex. We also discuss a possible trend between the softness of several bursts and the early behavior of the optical afterglow, in the sense that XRFs and X-ray rich GRBs seem to have a plateau phase or even a rising light curve. This can be naturally explained in models where XRFs are similar to GRBs but seen off the jet axis.

Supernovae are nature's grandest explosions and an astrophysical laboratory in which unique conditions exist that are not achievable on Earth. They are also the furnaces in which most of the elements heavier than carbon have been forged. Scientists have argued for decades about the physical mechanism responsible for these explosions. It is clear that the ultimate energy source is gravity, but the relative roles of neutrinos, fluid instabilities, rotation and magnetic fields continue to be debated.

We present u'g'r'i'BV photometry and optical spectroscopy of the Type Ib/Ic SN 2005bf covering the first ~100 days following discovery. The u'g'BV light curves displayed a double-peaked morphology, which is among Type Ib/Ic supernovae. The bolometric light curve indicates that SN 2005bf was a remarkably luminous event. Spectroscopically, SN 2005bf underwent a unique transformation from a Type Ic-like event at early times to a typical Type Ib supernova at later phases. The initial maximum in u'g'BV was accompanied by the presence in the spectrum of high velocity absorption lines of Fe II, Ca II, and H I. The photospheric velocity derived from spectra at early epochs was unusually low compared with ordinary Type Ib supernovae. We describe one-dimensional computer simulations which attempt to account for these remarkable properties. The most favored model is that of a very energetic (2 x 10^{51} erg), asymmetric explosion of a massive (8.3 M_sun) Wolf-Rayet WN star that had lost most of its hydrogen envelope. We speculate that an unobserved relativistic jet was launched producing a two-component explosion consisting of 1) a polar explosion containing a small fraction of the total mass and moving at high velocity, and 2) the explosion of the rest of the star. At first, only the polar explosion is observed, producing the initial maximum and the high velocity absorption-line spectrum resembling a Type Ic event. At late times, this fast-moving component becomes optically-thin, revealing the slower-moving explosion of the rest of the star and transforming the observed spectrum to that of a typical Type Ib supernova. If this scenario is correct, then SN 2005bf is the best example to date of a transition object between normal Type Ib/Ic supernovae and gamma ray bursts.

We investigate marginally stable nuclear burning on the surface of accreting neutron stars as an explanation for the mHz quasi-periodic oscillations (QPOs) observed from three low mass X-ray binaries. At the boundary between unstable and stable burning, the temperature dependence of the nuclear heating rate and cooling rate almost cancel. The result is an oscillatory mode of burning, with an oscillation period close to the geometric mean of the thermal and accretion timescales for the burning layer. We describe a simple one-zone model which illustrates this basic physics, and then present detailed multizone hydrodynamical calculations of nuclear burning close to the stability boundary using the KEPLER code. Our models naturally explain the characteristic 2 minute period of the mHz QPOs, and why they are seen only in a very narrow range of X-ray luminosities. The oscillation period is sensitive to the accreted hydrogen fraction and the surface gravity, suggesting a new way to probe these parameters. A major puzzle is that the accretion rate at which the oscillations appear in the theoretical models is an order of magnitude larger than the rate implied by the X-ray luminosity when the mHz QPOs are seen. We discuss the implications for our general understanding of nuclear burning on accreting neutron stars. One possibility is that the accreted material covers only part of the neutron star surface at luminosities Lx > ~1E37 erg/s.

One of the outstanding unsolved riddles of nuclear astrophysics is the origin of the so called ``p-process'' nuclei from A = 92 to 126. Both the lighter and heavier p-process nuclei are adequately produced in the neon and oxygen shells of ordinary Type II supernovae, but the origin of these intermediate isotopes, especially 92,94Mo and 96,98Ru, has long been mysterious. Here we explore the production of these nuclei in the neutrino-driven wind from a young neutron star. We consider such early times that the wind still contains a proton excess because the rates for electron neutrino and positron captures on neutrons are faster than those for the inverse captures on protons. Following a suggestion by Frohlich et al. 2005, we also include the possibility that, in addition to the protons, alpha-particles, and heavy seed, a small flux of neutrons is maintained by the reaction p(bar(nu_e),e+)n. This flux of neutrons is critical in bridging the long waiting points along the path of the rp-process by (n,p) and (n,gamma) reactions. Using the unmodified ejecta histories from a recent two-dimensional supernova model by Janka et al. 2003, we find synthesis of p-rich nuclei up to 102Pd. However, if the entropy of these ejecta is increased by a factor of two, the synthesis extends to 120Te. Still larger increases in entropy, that might reflect the role of magnetic fields or vibrational energy input neglected in the hydrodynamical model, result in the production of numerous r-, s-, and p-process nuclei up to A approximately 170, even in winds that are proton-rich.

We present a systematic survey of the capabilities of type Ia supernova explosion models starting from a number of flame seeds distributed around the center of the white dwarf star. To this end we greatly improved the resolution of the numerical simulations in the initial stages. This novel numerical approach facilitates a detailed study of multi-spot ignition scenarios with up to hundreds of ignition sparks. Two-dimensional simulations are shown to be inappropriate to study the effects of initial flame configurations. Based on a set of three-dimensional models, we conclude that multi-spot ignition scenarios may improve type Ia supernova models towards better agreement with observations. The achievable effect reaches a maximum at a limited number of flame ignition kernels as shown by the numerical models and corroborated by a simple dimensional analysis.

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.

The thermonuclear runaway that culminates in the explosion of a Chandrasekhar mass white dwarf as a Type Ia supernova begins centuries before the star actually explodes. Here, using a 3D anelastic code, we examine numerically the convective flow during the last minute of that runaway, a time that is crucial in determining just where and how often the supernova ignites. We find that the overall convective flow is dipolar, with the higher temperature fluctuations in an outbound flow preferentially on one side of the star. Taken at face value, this suggests an asymmetric ignition that may well persist in the geometry of the final explosion. However, we also find that even a moderate amount of rotation tends to fracture this dipole flow, making ignition over a broader region more likely. Though our calculations lack the resolution to study the flow at astrophysically relevant Rayleigh numbers, we also speculate that the observed dipolar flow will become less organized as the viscosity becomes very small. Motion within the dipole flow shows evidence of turbulence, suggesting that only geometrically large fluctuations (~1 km) will persist to ignite the runaway. We also examine the probability density function for the temperature fluctuations, finding evidence for a Gaussian, rather than exponential distribution, which suggests that ignition sparks may be strongly spatially clustered.

We study the detectability of primordial metal-free stars. Cosmological enrichment is a local process that takes place over an extended redshift range. While the duration of this transition depends on several unknown factors, in all cases late-forming metal-free stars are found in ~ 10^8 solar mass objects, which are just large enough to cool, but small enough to not be clustered near areas of previous star formation. We discuss the observational properties of these objects, some of which may have already been detected in ongoing Lyman-alpha surveys. If metal-free stars have masses between 140 and 260 solar masses, they are expected to end their lives as pair-production supernovae (PPSNe). Using the KEPLER code, we calculate a set of light curves for these objects. The peak luminosities of typical PPSNe are only slightly greater than those of Type Ia supernovae, but they remain bright much longer (~ 1 year) and have hydrogen lines. Ongoing supernova searches are able to place stringent limits on the fraction of very massive stars out to redshift 2. The planned Joint Dark Energy Mission satellite will be able to extend these constraints out to redshift 6

Those massive stars that, during their deaths, give rise to gamma-ray bursts (GRBs) must be endowed with an unusually large amount of angular momentum in their inner regions, one to two orders of magnitude greater than the ones that make common pulsars. Yet the inclusion of mass loss and angular momentum transport by magnetic torques during the precollapse evolution is known to sap the core of the necessary rotation. Here we explore the evolution of very rapidly rotating, massive stars, including stripped down helium cores that might result from mergers or mass transfer in a binary, and single stars that rotate unusually rapidly on the main sequence. For the highest possible rotation rates (about 400 km/s), a novel sort of evolution is encountered in which single stars mix completely on the main sequence, never becoming red giants. Such stars, essentially massive "blue stragglers", produce helium-oxygen cores that rotate unusually rapidly. Such stars might comprise roughly 1% of all stars above 10 solar masses and can, under certain circumstances retain enough angular momentum to make GRBs. Because this possibility is very sensitive to mass loss, GRBs will be much more probable in regions of low metallicity.

Gamma-ray burst (GRB) afterglows have provided important clues to the nature of these massive explosive events, providing direct information on the nearby environment and indirect information on the central engine that powers the burst. We report the discovery of two bright X-ray flares in GRB afterglows, including a giant flare comparable in total energy to the burst itself, each peaking minutes after the burst. These strong, rapid X-ray flares imply that the central engines of the bursts have long periods of activity, with strong internal shocks continuing for hundreds of seconds after the gamma-ray emission has ended.

Nonrotating, zero metallicity stars with initial masses 140 < M < 260 solar masses are expected to end their lives as pair-production supernovae (PPSNe), in which an electron-positron pair-production instability triggers explosive nuclear burning. Interest in such stars has been rekindled by recent theoretical studies that suggest primordial molecular clouds preferentially form stars with these masses. Since metal enrichment is a local process, the resulting PPSNe could occur over a broad range of redshifts, in pockets of metal-free gas. Using the implicit hydrodynamics code KEPLER, we have calculated a set of PPSN light curves that addresses the theoretical uncertainties and allows us to assess observational strategies for finding these objects at intermediate redshifts. The peak luminosities of typical PPSNe are only slightly greater than those of Type Ia, but they remain bright much longer (~ 1 year) and have hydrogen lines. Ongoing supernova searches may soon be able to limit the contribution of these very massive stars to < 1% of the total star formation rate density out to z=2 which already provides useful constraints for theoretical models. The planned Joint Dark Energy Mission satellite will be able to extend these limits out to z=6.

We extend a low Mach number hydrodynamics method developed for terrestrial combustion, to the study of thermonuclear flames in Type Ia supernovae. We discuss the differences between 2-D and 3-D Rayleigh-Taylor unstable flame simulations, and give detailed diagnostics on the turbulence, showing that the kinetic energy power spectrum obeys Bolgiano-Obukhov statistics in 2-D, but Kolmogorov statistics in 3-D. Preliminary results from 3-D reacting bubble calculations are shown, and their implications for ignition are discussed.

(Abridged) We report on HST observations of the late-time afterglow and host galaxy of GRB021004 (z=2.33). Although this GRB is one of the best observed so far, there is substantial disagreement between different interpretations of data sets on this burst in the literature. We have observed the field of GRB021004 with the HST at multiple epochs from 3 days until almost 10 months after the burst. With STIS spectroscopy we cover the spectral region 600-1700AA in the restframe. From the limit on the flux recovery bluewards of the Lyman-limit we constrain the HI column density to be above 1x10^18 cm-2 Based on ACS and NICMOS imaging we find that the afterglow evolved achromatically within the errors. The color changes observed by other authors during the first four days must be related to a stochastic phenomenon superimposed on an afterglow component with a constant spectral shape. This achromaticity implies that the cooling break has remained on the blue side of the optical part of the spectrum for at least two weeks after the burst. The optical-to-X-ray slope is consistent with being the same at 1.4 and 52.4 days after the burst. This indicates that the cooling frequency is constant and hence, according to fireball models, that the circumburst medium has a constant density profile. The late-time slope of the lightcurve is in the range alpha_2 = 1.8-1.9, and is inconsistent with a single power-law. This could be due to a late-time flattening caused by the transition to non-relativistic expansion or due to excess emission (a lightcurve `bump') 7 days after burst. The host is a starburst galaxy with no evidence for dust and with strong Ly-alpha emission. The afterglow was located very close (~100 pc) to the center of the host suggesting that the progenitor was associated with a circumnuclear starburst.

Flame instabilities play a dominant role in accelerating the burning front to a large fraction of the speed of sound in a Type Ia supernova. We present a three-dimensional numerical simulation of a Rayleigh-Taylor unstable carbon flame, following its evolution through the transition to turbulence. A low Mach number hydrodynamics method is used, freeing us from the harsh timestep restrictions imposed by sound waves. We fully resolve the thermal structure of the flame and its reaction zone, eliminating the need for a flame model. A single density is considered, $1.5\times 10^7 \gcc$, and half carbon\slash half oxygen fuel--conditions under which the flame propagated in the flamelet regime in our related two-dimensional study. We compare to a corresponding two-dimensional simulation, and show that while fire-polishing keeps the small features suppressed in two dimensions, turbulence wrinkles the flame on far smaller scales in the three-dimensional case, suggesting that the transition to the distributed burning regime occurs at higher densities in three dimensions. Detailed turbulence diagnostics are provided. We show that the turbulence follows a Kolmogorov spectrum and is highly anisotropic on the large scales, with a much larger integral scale in the direction of gravity. Furthermore, we demonstrate that it becomes more isotropic as it cascades down to small scales. Based on the turbulent statistics and the flame properties of our simulation, we compute the Gibson scale. We show the progress of the turbulent flame through a classic combustion regime diagram, indicating that the flame just enters the distributed burning regime near the end of our simulation.

The low luminosity radio emission of the unusually faint GRB 031203 has been argued to support the idea of a class of intrinsically sub-energetic gamma-ray bursts (GRBs), currently comprising two members. While low energy GRBs probably exist, we show that the collective prompt and multiwavelength observations of the afterglow of GRB 031203 do not necessarily require a sub-energetic nature for that event. In fact, the data are more consistent with a typical, powerful GRB seen at an angle of about twice the opening angle of the central jet. The (redshift corrected) peak energy, E_p, of GRB 031203 then becomes ~ 2 MeV, similar to many other GRBs.

GRB 020124 and GRB 030323 constitute half the sample of gamma-ray bursts with a measured redshift greater than 3. This paper presents the temporal and spectral properties of these two gamma-ray bursts detected and localized with HETE-2. While they have nearly identical redshifts (z=3.20 for GRB 020124, and z=3.37 for GRB 030323), these two GRBs span about an order of magnitude in fluence, thus sampling distinct regions of the GRB luminosity function. The properties of these two bursts are compared with those of the bulk of the GRB population detected by HETE-2. We also discuss the energetics of GRB 020124 and GRB 030323 and show that they are compatible with the Epeak - Eiso relation discovered by Amati et al. (2002). Finally, we compute the maximum redshifts at which these bursts could have been detected by HETE-2 and we address various issues connected with the detection and localization of high-z GRBs.

We report High Energy Transient Explorer 2 (HETE-2) Wide Field X-ray Monitor/French Gamma Telescope observations of XRF030723 along with observations of the XRF afterglow made using the 6.5m Magellan Clay telescope and the Chandra X-ray Observatory. The observed peak energy E_pk_obs of the nu F_nu burst spectrum is found to lie within (or below) the WXM 2-25 keV passband at 98.5% confidence, and no counts are detected above 30 keV. Our best fit value is E_pk_obs=8.4+3.5/-3.4 keV. The ratio of X-ray to Gamma-ray flux for the burst follows a correlation found for GRBs observed with HETE-2, and the duration of the burst is similar to that typical of long-duration GRBs. If we require that the burst isotropic equivalent energy E_iso and E_pk_rest satisfy the relation discovered by Amati et al. (2002), a redshift of z=0.38+0.36/-0.18 can be determined, in agreement with constraints determined from optical observations. We are able to fit the X-ray afterglow spectrum and to measure its temporal fade. Although the best-fit fade is shallower than the concurrent fade in the optical, the spectral similarity between the two bands indicates that the X-ray fade may actually trace the optical fade. If this is the case, the late time rebrightening observed in the optical cannot be due to a supernova bump. We interpret the prompt and afterglow X-ray emission as arising from a jetted GRB observed off-axis and possibly viewed through a complex circumburst medium due to a progenitor wind.

An exceptionally intense gamma-ray burst, GRB030329, was detected and localized by the instruments on board the High Energy Transient Explorer satellite (HETE) at 11:37:14 UT on 29 March 2003. The burst consisted of two \~10s pulses of roughly equal brightness and an X-ray tail lasting >100s. The energy fluence in the 30-400 keV energy band was 1.08e-4 erg/cm2, making GRB030329 one of the brightest GRBs ever detected. Communication of a 2 arcmin error box 73 minutes after the burst allowed the rapid detection of a counterpart in the optical, X-ray, radio and the ensuing discovery of a supernova with most unusual characteristics. Analyses of the burst lightcurves reveal the presence of a distinct, bright, soft X-ray component underlying the main GRB: the 2-10 keV fluence of this component is ~7e-6 erg/cm2. The main pulses of GRB030329 were preceded by two soft, faint, non-thermal bumps. We present details of the HETE observations of GRB030329.

As an explosion develops in the collapsed core of a massive star, neutrino emission drives convection in a hot bubble of radiation, nucleons, and pairs just outside a proto-neutron star. Shortly thereafter, neutrinos drive a wind-like outflow from the neutron star. In both the convective bubble and the early wind, weak interactions temporarily cause a proton excess (Ye > 0.5) to develop in the ejected matter. This situation lasts for at least the first second, and the approximately 0.05 - 0.1 solar masses that is ejected has an unusual composition that may be important for nucleosynthesis. Using tracer particles to follow the conditions in a two-dimensional model of a successful supernova explosion calculated by Janka, Buras & Rampp (2003), we determine the composition of this material. Most of it is helium and 56Ni. The rest is relatively rare species produced by the decay of proton-rich isotopes unstable to positron emission. In the absence of pronounced charged-current neutrino capture, nuclear flow will be held up by long-lived waiting point nuclei in the vicinity of 64Ge. The resulting abundance pattern can be modestly rich in a few interesting rare isotopes like 45Sc, 49Ti, and 64Zn. The present calculations imply yields that, when compared with the production of major species in the rest of the supernova, are about those needed to account for the solar abundance of 45Sc and 49Ti. Since the synthesis will be nearly the same in stars of high and low metallicity, the primary production of these species may have discernible signatures in the abundances of low metallicity stars. We also discuss uncertainties in the nuclear physics and early supernova evolution to which abundances of interesting nuclei are sensitive.

As a massive star evolves through multiple stages of nuclear burning on its way to becoming a supernova, a complex, differentially rotating structure is set up. Angular momentum is transported by a variety of classic instabilities, and also by magnetic torques from fields generated by the differential rotation. We present the first stellar evolution calculations to follow the evolution of rotating massive stars including, at least approximately, all these effects, magnetic and non-magnetic, from the zero-age main sequence until the onset of iron-core collapse. The evolution and action of the magnetic fields is as described by Spruit 2002 and a range of uncertain parameters is explored. In general, we find that magnetic torques decrease the final rotation rate of the collapsing iron core by about a factor of 30 to 50 when compared with the non-magnetic counterparts. Angular momentum in that part of the presupernova star destined to become a neutron star is an increasing function of main sequence mass. That is, pulsars derived from more massive stars will rotate faster and rotation will play a more dominant role in the star's explosion. The final angular momentum of the core is determined - to within a factor of two - by the time the star ignites carbon burning. For the lighter stars studied, around 15 solar masses, we predict pulsar periods at birth near 15 ms, though a factor of two range is easily tolerated by the uncertainties. Several mechanisms for additional braking in a young neutron star, especially by fall back, are also explored.

We describe and discuss the global properties of 45 gamma-ray bursts (GRBs) observed by HETE-2 during the first three years of its mission, focusing on the properties of X-Ray Flashes (XRFs) and X-ray-rich GRBs (XRRs). We find that the numbers of XRFs, XRRs, and GRBs are comparable. We find that the durations and the sky distributions of XRFs and XRRs are similar to those of GRBs. We also find that the spectral properties of XRFs and XRRs are similar to those of GRBs, except that the values of the peak energy $E^{\rm obs}_{\rm peak}$ of the burst spectrum in $\nu F_\nu$, the peak energy flux $\Fp$, and the energy fluence $S_E$ of XRFs are much smaller -- and those of XRRs are smaller -- than those of GRBs. Finally, we find that the distributions of all three kinds of bursts form a continuum in the [$S_E$(2-30 keV),$S_E$(30-400) keV]-plane, the [$S_E$(2-400 keV), $E_{\rm peak}$]-plane, and the [$F_{\rm peak}$(50-300 keV), $E_{\rm peak}$]-plane. These results provide strong evidence that all three kinds of bursts arise from the same phenomenon.

We analyze the possible existence of a pulsational instability excited by the $\epsilon$-mechanism during the last few centuries of evolution of a Chandrasekhar mass white dwarf prior to its explosion as a Type Ia supernova. Our analysis is motivated by the temperature sensitivity of the nuclear energy generation rate ($\sim T^{23}$) in a white dwarf whose structural adiabatic index is near 4/3. Based upon a linear stability analysis, we find that the fundamental mode and higher order radial modes are indeed unstable and that the fundamental mode has the shortest growth time scale. However, the growth time scale for such instability never becomes shorter than the evolutionary timescale. Therefore, even though the star \emph{is} pulsationally unstable, we do not expect these radial modes to have time to grow and to affect the structure and explosion properties of Type Ia supernovae.

We present the discovery and monitoring of the optical transient (OT) associated with GRB 020410. The fading OT was found by Hubble Space Telescope (HST) observations taken 28 and 65 days after burst at a position consistent with the X-ray afterglow. Subsequent re-examination of early ground based observations revealed that a faint OT was present 6 hours after burst, confirming the source association with GRB 020410. A deep non-detection after one week requires that the OT re-brightened between day 7 and day 28, and further late time HST data taken approximately 100 days after burst imply that it is very red.We compare both the flux and color of the excess with supernova models and show that the data are best explained by the presence of a Type Ib/c supernova at a redshift z ~ 0.5, which occured roughly coincident with the day of GRB.

New two- and three-dimensional calculations are presented of relativistic jet propagation and break out in massive Wolf-Rayet stars. Such jets are thought responsible for gamma-ray bursts. As it erupts, the jet is surrounded by a cocoon of less energetic, but still highly relativistic ejecta that expands and becomes visible at much larger polar angles. These less energetic ejecta may be the origin of X-ray flashes and other high-energy transients which will be visible to a larger fraction of the sky, albeit to a shorter distance than common gamma-ray bursts. Jet stability is also examined in three-dimensional calculations. If the jet changes angle by more than three degrees in several seconds, it will dissipate, producing a broad beam with inadequate Lorentz factor to make a common gamma-ray burst.

The attempt to secure conclusive, spectroscopic evidence for the GRB/SN connection has been a central theme in most GRB observing time proposals since the discovery of the very unusual GRB980425 associated with the peculiar type Ib/c SN1998bw. GRB030329 provided this evidence to everybody's satisfaction. In this contribution we show the results of a spectroscopic campaign of the supernova associated with GRB030329 carried out at ESOs Very Large Telescope. We also present preliminary results from a photometric and spectroscopic campaign targeting the X-ray Flash of July 23.

Planar flames are intrinsically unstable in open domains due to the thermal expansion across the burning front--the Landau-Darrieus instability. This instability leads to wrinkling and growth of the flame surface, and corresponding acceleration of the flame, until it is stabilized by cusp formation. We look at the Landau-Darrieus instability for C/O thermonuclear flames at conditions relevant to the late stages of a Type Ia supernova explosion. Two-dimensional direct numerical simulations of both single-mode and multi-mode perturbations using a low Mach number hydrodynamics code are presented. We show the effect of the instability on the flame speed as a function of both the density and domain size, demonstrate the existence of the small scale cutoff to the growth of the instability, and look for the proposed breakdown of the non-linear stabilization at low densities. The effects of curvature on the flame as quantified through measurements of the growth rate and computation of the corresponding Markstein number. While accelerations of a few percent are observed, they are too small to have any direct outcome on the supernova explosion.

A Type Ia supernova explosion likely begins as a nuclear runaway near the center of a carbon-oxygen white dwarf. The outward propagating flame is unstable to the Landau-Darrieus, Rayleigh-Taylor, and Kelvin-Helmholtz instabilities, which serve to accelerate it to a large fraction of the speed of sound. We investigate the Rayleigh-Taylor unstable flame at the transition from the flamelet regime to the distributed-burning regime, around densities of $10^7$ g/cc, through detailed, fully resolved simulations. A low Mach number, adaptive mesh hydrodynamics code is used to achieve the necessary resolution and long time scales. As the density is varied, we see a fundamental change in the character of the burning--at the low end of the density range the Rayleigh-Taylor instability dominates the burning, whereas at the high end the burning suppresses the instability. In all cases, significant acceleration of the flame is observed, limited only by the size of the domain we are able to study. We discuss the implications of these results on the potential for a deflagration to detonation transition.

(Abrigded) We report X-ray studies of the environs of SN 1998bw and GRB 980425 using the Chandra X-Ray Observatory 1281 days after the GRB. Combining our observation of the supernova with others of the GRB afterglow, a smooth X-ray light curve, spanning ~1300 days, is obtained by assuming the burst and supernova were coincident at 35.6 Mpc. When this X-ray light curve is compared with those of the X-ray ``afterglows'' of ordinary GRBs, X-ray Flashes, and ordinary supernovae, evidence emerges for at least two classes of lightcurves, perhaps bounding a continuum. By three to ten years, all these phenomena seem to converge on a common X-ray luminosity, possibly indicative of the supernova underlying them all. This convergence strengthens the conclusion that SN 1998bw and GRB 980425 took place in the same object. One possible explanation for the two classes is a (nearly) standard GRB observed at different angles, in which case X-ray afterglows with intermediate luminosities should eventually be discovered. Finally, we comment on the contribution of GRB afterglows to the ULX source population.

The HETE-2 mission has been highly productive. It has observed more than 250 GRBs so far. It is currently localizing 25 - 30 GRBs per year, and has localized 43 GRBs to date. Twenty-one of these localizations have led to the detection of X-ray, optical, or radio afterglows, and as of now, 11 of the bursts with afterglows have known redshifts. HETE-2 has confirmed the connection between GRBs and Type Ic supernovae, a singular achievement and certainly one of the scientific highlights of the mission so far. It has provided evidence that the isotropic-equivalent energies and luminosities of GRBs are correlated with redshift, implying that GRBs and their progenitors evolve strongly with redshift. Both of these results have profound implications for the nature of GRB progenitors and for the use of GRBs as a probe of cosmology and the early universe. HETE-2 has placed severe constraints on any X-ray or optical afterglow of a short GRB. It is also solving the mystery of "optically dark' GRBs, and revealing the nature of X-ray flashes.

The HETE-2 mission has been highly productive. It has observed more than 250 GRBs so far. It is currently localizing 25 - 30 GRBs per year, and has localized 43 GRBs to date. Twenty-one of these localizations have led to the detection of X-ray, optical, or radio afterglows, and as of now, 11 of the bursts with afterglows have redshift determinations. HETE-2 has also observed more than 45 bursts from soft gamma-ray repeaters, and more than 700 X-ray bursts. HETE-2 has confirmed the connection between GRBs and Type Ic supernovae, a singular achievement and certainly one of the scientific highlights of the mission so far. It has provided evidence that the isotropic-equivalent energies and luminosities of GRBs may be correlated with redshift; such a correlation would imply that GRBs and their progenitors evolve strongly with redshift. Both of these results have profound implications for the nature of GRB progenitors and for the use of GRBs as a probe of cosmology and the early universe. HETE-2 has placed severe constraints on any X-ray or optical afterglow of a short GRB. It has made it possible to explore the previously unknown behavior optical afterglows at very early times, and has opened up the era of high-resolution spectroscopy of GRB optical afterglows. It is also solving the mystery of "optically dark" GRBs, and revealing the nature of X-ray flashes (XRFs).

We report HETE-2 WXM/FREGATE observations of the X-ray flash, XRF 020903. This event was extremely soft: the ratio log(S_X/S_gamma) = 0.7, where S_X and S_gamma are the fluences in the 2-30 and 30-400 keV energy bands, is the most extreme value observed so far by HETE-2. In addition, the spectrum has an observed peak energy E^{obs}_{peak} < 5.0 keV (99.7 % probability upper limit) and no photons were detected above ~10 keV. The burst is shorter at higher energies, which is similar to the behavior of long GRBs. We consider the possibility that the burst lies at very high redshift and that the low value of E^{obs}_{peak} is due to the cosmological redshift, and show that this is very unlikely. We find that the properties of XRF 020903 are consistent with the relation between the fluences S(7-30 keV) and S(30-400 keV) found by Barraud et al. for GRBs and X-ray-rich GRBs, and are consistent with the extension by a decade of the hardness-intensity correlation (Mallozzi et al. 1995) found by the same authors. Assuming that XRF 020903 lies at a redshift z = 0.25 as implied by the host galaxy of the candidate optical and radio afterglows of this burst, we find that the properties of XRF 020903 are consistent with an extension by a factor ~300 of the relation between the isotropic-equivalent energy E_iso and the peak E_peak of the nu F_nu spectrum (in the source frame of the burst) found by Amati et al. for GRBs. The results presented in this paper therefore provide evidence that XRFs, X-ray-rich GRBs, and GRBs form a continuum and are a single phenomenon. The results also impose strong constraints on models of XRFs and X-ray-rich GRBs.

SN 2003dh, one of the most luminous supernovae ever recorded, and the one with the highest measured velocities, accompanied gamma-ray burst 030329. Its rapid rise to maximum and equally rapid decline pose problems for any spherically symmetric model. We model the supernova here as a very energetic, polar explosion that left the equatorial portions of the star almost intact. The total progenitor mass was much greater than the mass of high-velocity ejecta, and the total mass of 56-Ni synthesized was about 0.5 solar masses. Such asymmetries and nickel masses are expected in the collapsar model. A ``composite two-dimensional'' model is calculated that agrees well with the characteristics of the observed light curve. The mass of 56-Ni required for this light curve is 0.55 solar masses and the total explosion energy, 26 x 10**51 erg.

A bright, x-ray-rich GRB was detected by HETE-2 at 11:18:34.03 UT on 11 Dec 2002. The WXM localization was to 14' and relayed to the GCN 22 s after the start of the burst. The ground SXC localization was within 2' of R.A. 08h 09m 00s, Dec 06d 44' 20" (J2000). GRB021211 consists of a single, FRED-like pulse with t90s of 2.3 s (85-400 keV) and 8.5 s (2-10 keV). The peak photon number and photon energy fluxes in the 2-400 keV band, are 34.0 +/- 1.8 ph cm^-2 s^-1 and 1.68 +/- 0.11 x 10^-6 erg cm^-2 s^-1, respectively. The energy fluences in the 2-30 keV and 30-400 kev energy bands are S_X = 1.36 +/- 0.05 x 10^-6 erg cm^-2 and S_gamma 2.17 +/- 0.15 x 10^-6 erg cm^-2, respectively. Thus GRB021211 is an X-ray-rich GRB (S_X/S_gamma = 0.63 > 0.32). The spectrum is well-fit by a Band function (alpha = -0.805, beta = -2.37, E_peak = 46.8 keV). The prompt localization allowed the detection of an optical afterglow for what would otherwise have been an ``optically dark'' GRB. GRB 021211 demonstrates that some fraction of burst afterglows are ``optically dark'' because their optical afterglows at times > 1 hr after the burst are very faint, and thus have often escaped detection. GRB 021211 shows that such ``optically dim'' bursts can have very bright afterglows at times < 20 min after the burst.

The observable properties of a Type Ia supernova are sensitive to how the nuclear runaway ignites in a Chandrasekhar mass white dwarf - at a single point at its center, off-center, or at multiple points and times. We present a simple analytic model for the runaway based upon a combination of stellar mixing-length theory and recent advances in understanding Rayleigh-Benard convection. The convective flow just prior to runaway is likely to have a strong dipolar component, though higher multipoles may contribute appreciably at the very high Rayleigh number (10$^{25}$) appropriate to the white dwarf core. A likely outcome is multi-point ignition with an exponentially increasing number of ignition points during the few tenths of a second that it takes the runaway to develop. The first sparks ignite approximately 150 - 200 km off center, followed by ignition at smaller radii. Rotation may be important to break the dipole asymmetry of the ignition and give a healthy explosion.

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.

Multi-zone models of Type I X-ray bursts are presented that use an adaptive nuclear reaction network of unprecedented size, up to 1300 isotopes. Sequences of up to 15 bursts are followed for two choices of accretion rate and metallicity. At 0.1 Eddington (and 0.02 Eddington for low metallicity), combined hydrogen-helium flashes occur. The rise times, shapes, and tails of these light curves are sensitive to the efficiency of nuclear burning at various waiting points along the rp-process path and these sensitivities are explored. The bursts show "compositional inertia", in that their properties depend on the fact that accretion occurs onto the ashes of previous bursts which contain left-over hydrogen, helium and CNO nuclei. This acts to reduce the sensitivity of burst properties to metallicity. For the accretion rates studied, only the first anomalous burst in one model produces nuclei as heavy as A=100, other bursts make chiefly nuclei with A~64. The amount of carbon remaining after hydrogen-helium bursts is typically <1% by mass, and decreases further as the ashes are periodically heated by subsequent bursts. At the lower accretion rate of 0.02 Eddington and solar metallicity, the bursts ignite in a hydrogen-free helium layer. At the base of this layer, up to 90% of the helium has already burned to carbon prior to the unstable ignition. These helium-ignited bursts have briefer, brighter light curves with shorter tails, very rapid rise times (<0.1 s), and ashes lighter than the iron group.

Over the past five years evidence has mounted that long-duration (> 2 s) gamma-ray bursts (GRBs)--the most brilliant of all astronomical explosions--signal the collapse of massive stars in our Universe. This evidence was originally based on the probable association of one unusual GRB with a supernova, but now includes the association of GRBs with regions of massive star formation in distant galaxies, the appearance of supernova-like 'bumps' in the optical afterglow light curves of several bursts and lines of freshly synthesized elements in the spectra of a few X-ray afterglows. These observations support, but do not yet conclusively demonstrate, the idea that long-duration GRBs are associated with the deaths of massive stars, presumably arising from core collapse. Here we report evidence that a very energetic supernova (a hypernova) was temporally and spatially coincident with a GRB at redshift z = 0.1685. The timing of the supernova indicates that it exploded within a few days of the GRB, strongly suggesting that core-collapse events can give rise to GRBs, thereby favouring the 'collapsar' model.

Rotation in massive stars has been studied on the main sequence and during helium burning for decades, but only recently have realistic numerical simulations followed the transport of angular momentum that occurs during more advanced stages of evolution. The results affect such interesting issues as whether rotation is important to the explosion mechanism, whether supernovae are strong sources of gravitational radiation, the star's nucleosynthesis, and the initial rotation rate of neutron stars and black holes. We find that when only hydrodynamic instabilities (shear, Eddington-Sweet, etc.) are included in the calculation, one obtains neutron stars spinning at close to critical rotation at their surface -- or even formally in excess of critical. When recent estimates of magnetic torques (Spruit 2002) are added, however, the evolved cores spin about an order of magnitude slower. This is still more angular momentum than observed in young pulsars, but too slow for the collapsar model for gamma-ray bursts.

One of the most dramatic possible consequences of stellar rotation is its influence on stellar death, particularly of massive stars. If the angular momentum of the iron core when it collapses is such as to produce a neutron star with a period of 5 ms or less, rotation will have important consequences for the supernova explosion mechanism. Still shorter periods, corresponding to a neutron star rotating at break up, are required for the progenitors of gamma-ray bursts (GRBs). Current stellar models, while providing an excess of angular momentum to pulsars, still fall short of what is needed to make GRBs. The possibility of slowing young neutron stars in ordinary supernovae by a combination of neutrino-powered winds and the propeller mechanism is discussed. The fall back of slowly moving ejecta during the first day of the supernova may be critical. GRBs, on the other hand, probably require stellar mergers for their production and perhaps less efficient mass loss and magnetic torques than estimated thus far.

We examine the propagation of 2-dimensional relativistic jets through the stellar progenitor in the collapsar model for gamma-ray bursts. In agreement with previous studies, we find that relativistic jets are collimated by their passage through the stellar mantle. Interaction of these jets with the star and their own cocoons also causes mixing that sporadically decelerates the flow. We speculate that this mixing instability is chiefly responsible for the variable Lorentz factor needed in the internal shock model and for the complex light curves seen in many GRBs. In all cases studied, the jet is shocked deep inside the star following a brief period of adiabatic expansion. The jet that finally emerges from the star thus has a moderate Lorentz factor, modulated by mixing, and a very large internal energy. In a second series of calculations, we follow the escape of that sort of jet. Because of the large ratio of internal to kinetic energy in both the jet and its cocoon, the opening angle of the final jet is significantly greater than at breakout. A small amount of material emerges at large angles, but with a Lorentz factor still sufficiently large to make a weak GRB. This leads us to propose a "unified model" in which a variety of high energy transients, ranging from x-ray flashes to "classic" GRBs, may be seen depending upon the angle at which a standard collapsar is observed. We also speculate that the breakout of a relativistic jet and its collision with the stellar wind will produce a brief transient with properties similar to the class of "short-hard" GRBs. Implications of our calculations for GRB light curves, the luminosity-variability relation, and the GRB-supernova association are also discussed. (Abridged)

We simulate shock propagation in a wide range of core-collapsing presupernovae: from compact WR stars exploding as SNe Ib/c through very extended envelopes of the narrow-line SNe IIn. We find that the same physical phenomenon of radiating shocks can produce outbursts of X-ray radiation (with photon energy $3kT \sim 1$ keV) lasting only a second in SNe Ib/c, as well as a very high flux of visual light, lasting for months, in SNe IIn.

Numerical modeling of shock breakout in compact presupernovae of type Ib/c is done by two independent methods. Peak color temperatures reach 4 to 8 million K for explosion energies 1 to 9 foe with a typical local time-scale of 0.03 seconds. In the obverver's frame the burst is smeared over few seconds due to light travel time correction R/c.

How massive stars die -- what sort of explosion and remnant each produces -- depends chiefly on the masses of their helium cores and hydrogen envelopes at death. For single stars, stellar winds are the only means of mass loss, and these are chiefly a function of the metallicity of the star. We discuss how metallicity, and a simplified prescription for its effect on mass loss, affects the evolution and final fate of massive stars. We map, as a function of mass and metallicity, where black holes and neutron stars are likely to form and where different types of supernovae are produced. Integrating over an initial mass function, we derive the relative populations as a function of metallicity. Provided single stars rotate rapidly enough at death, we speculate upon stellar populations that might produce gamma-ray bursts and jet-driven supernovae.

We look at the spreading of localized temperature perturbations in the accreted fuel layer of a non-rotating neutron star. The pressure at the base of the accreted fuel layer is large and the material is only partially degenerate. Any temperature perturbations and resulting pressure gradients will lead to enormous accelerations (both laterally and vertically) of the material in the fuel layer. If the burning of this fuel cannot proceed more rapidly than the spreading of this perturbation, then localized burning cannot take place, and it is likely that the ignition would have to proceed simultaneously throughout the envelope. We present some multidimensional simulations of the spreading of temperature perturbations in a helium atmosphere on a neutron star.

The last five years have seen growing challenges to the traditional paradigm of a core collapse supernova powered by the neutrino emission of a young proto-neutron star. Chief among these challenges are gamma-ray bursts (GRBs) and the supernovae that seem to accompany them. Here we review some recent - and not so recent - models for GRBs and supernovae in which strong magnetic fields, rotation, or accretion into a black hole play a role. The conditions for these energetic explosions are special and, at this point, there is no compelling reason to invoke them in the general case. That is, 99% of supernovae may still operate in the traditional fashion.

We review the current basic picture of the evolution of massive stars and how their evolution and structure changes as a function of initial mass. We give an overview of the fate of modern (Pop I) and primordial (Pop III) stars with emphasis on massive and very massive stars. For single stars we show how the type of explosions, the type of remnant and their frequencies changes for different initial metallicities.

After briefly describing the anelastic approximation and Glatzmaier's code, we present results from our preliminary studies of core convection during the hydrogen burning phase of a 15 solar mass star, as well as our most recent results concerning convection in the oxygen shell of a 25 solar mass star.

Growing observational evidence supports the proposition that gamma-ray bursts (GRBs) are powered by relativistic jets from massive helium stars whose cores have collapsed to black holes and an accretion disk (collapsars). We model the propagation of relativistic jets through the stellar progenitor and its wind using a two-dimensional special relativistic hydrodynamics code based on the PPM formalism. The jet emerges from the star with a plug in front and a cocoon surrounding it. During its propagation outside the star, the jet gains high Lorentz factor as its internal energy is converted into kinetic energy while the cocoon expands both outwards and sideways. External shocks between the cocoon and the stellar wind can produce $\gamma$-ray and hard x-ray transients. The interaction of the jet beam and the plug will also affect both of them substantially, and may lead to short-hard GRBs. Internal shocks in the jet itself may make long-soft GRBs.

We follow the nuclear reactions that occur in the accretion disks of stellar mass black holes that are accreting at a very high rate, 0.01 to 1 solar masses per second, as is realized in many current models for gamma-ray bursts (GRBs). The degree of neutronization in the disk is a sensitive function of the accretion rate, black hole mass, Kerr parameter, and disk viscosity. For high accretion rates and low viscosity, material arriving at the black hole will consist predominantly of neutrons. This degree of neutronization will have important implications for the dynamics of the GRB producing jet and perhaps for the synthesis of the r-process. For lower accretion rates and high viscosity, as might be appropriate for the outer disk in the collapsar model, neutron-proton equality persists allowing the possible synthesis of 56Ni in the disk wind. 56Ni must be present to make any optically bright Type Ib supernova, and in particular those associated with GRBs.

The {\it HETE-2} (hereafter \HETE) French Gamma Telescope (FREGATE) and the Wide-field X-ray Monitor (WXM) instruments detected a short ($t_{50} = 360$ msec in the FREGATE 85-300 keV energy band), hard gamma-ray burst (GRB) that occurred at 1578.72 SOD (00:26:18.72 UT) on 31 May 2002. The WXM flight localization software produced a valid location in spacecraft (relative) coordinates. However, since no on-board real-time star camera aspect was available, an absolute localization could not be disseminated. A preliminary localization was reported as a GCN Position Notice at 01:54:22 UT, 88 min after the burst. Further ground analysis produced a refined localization, which can be expressed as a 90% confidence rectangle that is 67 arcminutes in RA and 43 arcminutes in Dec (90% confidence region), centered at RA = +15$^{\rm h}$ 14$^{\rm m}$ 45$^{\rm s}$, Dec = -19$^\circ$ 21\arcmin 35\arcsec (J2000). An IPN localization of the burst was disseminated 18 hours after the GRB (Hurley et al. 2002b). A refined IPN localization was disseminated $\approx$ 5 days after the burst. This hexagonal-shaped localization error region is centered on RA = 15$^{\rm h}$ 15$^{\rm m}$ 03.57$^{\rm s}$, -19$^\circ$ 24\arcmin 51.00\arcsec (J2000), and has an area of $\approx$ 22 square arcminutes (99.7% confidence region). The prompt localization of this short, hard GRB by \HETE and the anti-Sun pointing of the \HETE instruments, coupled with the refinement of the localization by the IPN, has made possible rapid follow-up observations of the burst at radio, optical, and X-ray wavelengths.

Leading models for the "central engine" of long, soft gamma-ray bursts (GRBs) are briefly reviewed with emphasis on the collapsar model. Growing evidence supports the hypothesis that GRBs are a supernova-like phenomenon occurring in star forming regions, differing from ordinary supernovae in that a large fraction of their energy is concentrated in highly relativistic jets. The possible progenitors and physics of such explosions are discussed and the important role of the interaction of the emerging relativistic jet with the collapsing star is emphasized. This interaction may be responsible for most of the time structure seen in long, soft GRBs. What we have called "GRBs" may actually be a diverse set of phenomena with a key parameter being the angle at which the burst is observed. GRB 980425/SN 1988bw and the recently discovered hard x-ray flashes may be examples of this diversity.

We study the evolution of stars that may be the progenitors of common (long-soft) GRBs. Bare rotating helium stars, presumed to have lost their envelopes due to winds or companions, are followed from central helium ignition to iron core collapse. Including realistic estimates of angular momentum transport (Heger, Langer, & Woosley 2000) by non-magnetic processes and mass loss, one is still able to create a collapsed object at the end with sufficient angular momentum to form a centrifugally supported disk, i.e., to drive a collapsar engine. However, inclusion of current estimates of magnetic torques (Spruit 2002) results in too little angular momentum for collapsars.

A next generation Gamma Ray Burst (GRB) mission to follow the upcoming Swift mission is described. The proposed Energetic X-ray Imaging Survey Telescope, EXIST, would yield the limiting (practical) GRB trigger sensitivity, broad-band spectral and temporal response, and spatial resolution over a wide field. It would provide high resolution spectra and locations for GRBs detected at GeV energies with GLAST. Together with the next generation missions Constellation-X, NGST and LISA and optical-survey (LSST) telescopes, EXIST would enable GRBs to be used as probes of the early universe and the first generation of stars. EXIST alone would give ~10-50" positions (long or short GRBs), approximate redshifts from lags, and constrain physics of jets, orphan afterglows, neutrinos and SGRs.

We present the first calculations to follow the evolution of all stable nuclei and their radioactive progenitors in stellar models computed from the onset of central hydrogen burning through explosion as Type II supernovae. Calculations are performed for Pop I stars of 15, 19, 20, 21, and 25 M_sun using the most recently available experimental and theoretical nuclear data, revised opacity tables, neutrino losses, and weak interaction rates, and taking into account mass loss due to stellar winds. A novel ``adaptive'' reaction network is employed with a variable number of nuclei (adjusted each time step) ranging from about 700 on the main sequence to more than 2200 during the explosion. The network includes, at any given time, all relevant isotopes from hydrogen through polonium (Z=84). Even the limited grid of stellar masses studied suggests that overall good agreement can be achieved with the solar abundances of nuclei between 16O and 90Zr. Interesting discrepancies are seen in the 20 M_sun model and, so far, only in that model, that are a consequence of the merging of the oxygen, neon, and carbon shells about a day prior to core collapse. We find that, in some stars, most of the ``p-process'' nuclei can be produced in the convective oxygen burning shell moments prior to collapse; in others, they are made only in the explosion. Serious deficiencies still exist in all cases for the p-process isotopes of Ru and Mo.

We calculate Gamma-Ray Burst afterglow light-curves from a relativistic jet of initial opening angle theta_0, as seen by observers at a wide range of viewing angles, theta_obs, from the jet axis. We describe three increasingly more realistic models and compare the resulting light-curves. An observer at theta_obs < theta_0 should see a light curve very similar to that for an on-axis observer. An observer at theta_obs > theta_0 should see a rising light curve at early times, the flux peaking when the jet Lorentz factor sim 1/theta_obs. After this time the flux is not very different from that seen by an on-axis observer. A strong linear polarization (<40%) may occur near the peak in the light curve, and slowly decay with time. We show that if GRB jets have a universal energy, then orphan afterglows associated with off-axis jets should be seen up to a constant theta_obs, therefore the detection rate of orphan afterglows would be proportional to the true GRB rate. We also discuss the proposed connection between supernova 1998bw and GRB 980425.

On September 21 at 18950.56 SOD (05:15:50.56) UT the FREGATE gamma-ray instrument on the High Energy Transient Explorer (HETE) detected a bright gamma-ray burst (GRB). The burst was also seen by the X-detector on the WXM X-ray instrument and was therefore well-localized in the X direction; however, the burst was outside the fully-coded field-of-view of the WXM Y-detector, and therefore information on the Y direction of the burst was limited. Cross-correlation of the HETE and Ulysses time histories yielded an Interplanetary Network (IPN) annulus that crosses the HETE error strip at a ~45 degree angle. The intersection of the HETE error strip and the IPN annulus produces a diamond-shaped error region for the location of the burst having an area of 310 square arcminutes. Based on the FREGATE and WXM light curves, the duration of the burst is characterized by a t90 = 18.4 s in the WXM 4 - 25 keV energy range, and 23.8 s and 21.8 s in the FREGATE 6 - 40 and 32 - 400 keV energy ranges, respectively. The fluence of the burst in these same energy ranges is 4.8 10^{-6}, 5.5 10^{-6}, and 11.4 10^{-6} erg cm^{-2}, respectively. Subsequent optical and radio observations by ground-based observers have identified the afterglow of GRB010921 and determined an apparent redshift of z = 0.450.

While the modern stellar IMF shows a rapid decline with increasing mass, theoretical investigations suggest that very massive stars (>100 solar masses) may have been abundant in the early universe. Other calculations also indicate that, lacking metals, these same stars reach their late evolutionary stages without appreciable mass loss. After central helium burning, they encounter the electron-positron pair instability, collapse, and burn oxygen and silicon explosively. If sufficient energy is released by the burning, these stars explode as brilliant supernovae with energies up to 100 times that of an ordinary core collapse supernova. They also eject up to 50 solar masses of radioactive Ni56. Stars less massive than 140 solar masses or more massive than 260 solar masses should collapse into black holes instead of exploding, thus bounding the pair-creation supernovae with regions of stellar mass that are nucleosynthetically sterile. Pair-instability supernovae might be detectable in the near infrared out to redshifts of 20 or more and their ashes should leave a distinctive nucleosynthetic pattern.

Growing evidence suggests that the first generation of stars may have been quite massive (~100-300 M_sun). Could these stars have left a distinct nucleosynthetic signature? We explore the nucleosynthesis of helium cores in the mass range M_He=64 to 133 Msun, corresponding to main-sequence star masses of approximately 140 to 260 M_sun. Above M_He=133 M_sun, without rotation and using current reaction rates, a black hole is formed and no nucleosynthesis is ejected. For lighter helium core masses, ~40 to 63 M_sun, violent pulsations occur, induced by the pair instability and accompanied by supernova-like mass ejection, but the star eventually produces a large iron core in hydrostatic equilibrium. It is likely that this core, too, collapses to a black hole, thus cleanly separating the heavy element nucleosynthesis of pair instability supernovae from those of other masses, both above and below. Indeed, black hole formation is a likely outcome for all Population III stars with main sequence masses between about 25 M_sun and 140 M_sun (M_He = 9 to 63 M_sun) as well as those above 260 M_sun. Nucleosynthesis in pair-instability supernovae varies greatly with the mass of the helium core which determines the maximum temperature reached during the bounce. At the upper range of exploding core masses, a maximum of 57 M_sun of Ni56 is produced making these the most energetic and the brightest thermonuclear explosions in the universe. Integrating over a distribution of masses, we find that pair instability supernovae produce a roughly solar distribution of nuclei having even nuclear charge, but are remarkably deficient in producing elements with odd nuclear charge. Also, essentially no elements heavier than zinc are produced due to a lack of s- and r-processes.

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.

In order to study the processes creating intermediate and heavy nuclei in massive stars it is necessary to provide neutron capture cross sections and reaction rates close to stability and for moderately unstable neutron-rich nuclei. Furthermore, one has to know the efficiency of neutron-releasing reactions in the main evolutionary phases of a massive star. We present simulations of the nucleosynthesis in a 15 and 25 solar mass star, for the first time followed completely from main sequence hydrogen burning until the type II supernova explosion including all nuclides up to Bi. Theoretical reaction rates were calculated with the NON-SMOKER code, providing a complete library of Hauser-Feshbach cross sections and rates for nuclear and astrophysical applications. Experimental rates at stability were taken from different sources. The impact of uncertainties in the rates on nucleosynthesis are illustrated by two examples, the reactions 62Ni(n,gamma)63Ni and 22Ne(alpha,n)25Mg.

The stability of metal-free very massive stars ($Z$ = 0; $M = 120 - 500 \msol$) is analyzed and compared with metal-enriched stars. Such zero-metal stars are unstable to nuclear-powered radial pulsations on the main sequence, but the growth time scale for these instabilities is much longer than for their metal-rich counterparts. Since they stabilize quickly after evolving off the ZAMS, the pulsation may not have sufficient time to drive appreciable mass loss in Z = 0 stars. For reasonable assumptions regarding the efficiency of converting pulsational energy into mass loss, we find that, even for the larger masses considered, the star may die without losing a large fraction of its mass. We find a transition between the $\epsilon$- and $\kappa$-mechanisms for pulsational instability at $Z\sim 2\E{-4} - 2\E{-3}$. For the most metal-rich stars, the $\kappa$-mechanism yields much shorter $e$-folding times, indicating the presence of a strong instability. We thus stress the fundamental difference of the stability and late stages of evolution between very massive stars born in the early universe and those that might be born today.

Recent shell-model calculations of weak-interaction rates for nuclei in the mass range A = 45 - 65 have resulted in substantial revisions to the hitherto standard set of Fuller, Fowler, & Newman (FFN). In particular, key electron-capture rates, such as that for Co60 are much smaller. We consider here the effects of these revised rates on the presupernova (post-oxygen burning) evolution of massive stars in the mass range 11 to 40 M_sun. Moreover, we include, for the first time in models by our group, the effects of modern rates for beta-decay in addition to electron capture and positron emission. Values for the central electron mole number at the time of iron core collapse in the new models are typically larger, by delta Y_e = 0.005 to 0.015, than those of Woosley & Weaver 1995, with a tendency for the more massive models to display larger differences. About half of this change is a consequence of including beta-decay; the other half, result of the smaller rates for electron capture. Unlike what might be expected solely on basis of the larger Y_e's, the new iron core masses are systematically smaller owing to a decrease in the entropy in the outer iron core. The changes in iron core mass range from zero to 0.1 M_sun. We also observe, as predicted by Aufderheide et al. (1994), a tendency towards beta-equilibrium just prior to the collapse of the core, and the subsequent loss of that equilibrium as core collapse proceeds. We discuss the key weak reaction rates, both beta-decay and electron-capture, responsible for the evolution of Y_e and make suggestions for future measurements.

We investigate the evolution, final fate, and nucleosynthetic yields of rotating and non-rotating very massive stars (VMS) of zero metallicity. First we address the issue of mass loss during hydrogen burning due to vibrational instabilities. We find that these objects are much more stable than what was found in previous studies of VMS of solar composition, and expect only negligible mass loss driven by the pulsations. As these stars thus reach the end of their evolution with massive helium cores, they encounter the pair-creation instability. We find that for helium core masses of ~64...133 solar masses these stars are completely disrupted with explosion energies of up to ~1E53 erg and eject up to ~60 solar masses of Ni56 Stars with more massive helium cores collapse into black holes. We present the first calculations that follow the collapse of such a massive rotating star and predict that X-ray burst and significant gravitational wave emission could result.

We have performed the first calculations to follow the evolution of all stable nuclei and their radioactive progenitors in a finely-zoned stellar model computed from the onset of central hydrogen burning through explosion as a Type II supernova. Calculations were done for 15, 20, and 25 solar masses Pop I stars using the most recently available set of experimental and theoretical nuclear data, revised opacity tables, and taking into account mass loss due to stellar winds. Here results are presented for one 15 solar masses model.

Preliminary results of a new set of stellar evolution and nucleosynthesis calculations for massive stars are presented. These results were obtained with an extended reaction network up to Bi. The discussion focuses on the importance of nuclear rates in pre- and post-explosive nucleosynthesis. The need for further experiments to study specific reactions and nuclear properties (optical alpha+nucleus potentials) is emphasized.

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.

Improved values for stellar weak interaction rates have been recently calculated based upon a large shell model diagonalization. Using these new rates (for both beta decay and electron capture), we have examined the presupernova evolution of massive stars in the range 15-40 Msun. Comparing our new models with a standard set of presupernova models by Woosley and Weaver, we find significantly larger values for the electron-to-baryon ratio Ye at the onset of collapse and iron core masses reduced by approximately 0.1 Msun. The inclusion of beta-decay accounts for roughly half of the revisions, while the other half is a consequence of the improved nuclear physics. These changes will have important consequences for nucleosynthesis and the supernova explosion mechanism.

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.

We present the first calculations to follow the evolution of all stable isotopes (and their abundant radioactive progenitors) in a finely zoned stellar model computed from the onset of central hydrogen burning through explosion as a Type II supernova. The calculations were performed for a 15 solar mass Pop I star using the most recently available set of experimental and theoretical nuclear data, revised opacity tables, and taking into account mass loss due to stellar winds. We find the approximately solar production of proton-rich isotopes above a mass number of A=120 due to the gamma-process. We also find a weak s-process, which along with the gamma-process and explosive helium and carbon burning, produces nearly solar abundances of almost all nuclei from A=60 to 85. A few modifications of the abundances of heavy nuclei above mass 90 by the s-process are also noted and discussed. New weak rates lead to significant alteration of the properties of the presupernova core.

We present a statistical model which shows the influence of turbulence on a thermonuclear flame propagating in C+O white dwarf matter. Based on a Monte Carlo description of turbulence, it provides a method for investigating the physics in the so-called distributed burning regime. Using this method we perform numerical simulations of turbulent flames and show that in this particular regime the flamelet model for the turbulent flame velocity loses its validity. In fact, at high turbulent intensities burning in the distributed regime can lead to a deceleration of the turbulent flame and thus induces a competing process to turbulent effects that cause a higher flame speed. It is also shown that in dense C+O matter turbulent heat transport is described adequately by the Peclet number, Pe, rather than by the Reynolds number, which means that flame propagation is decoupled from small-scale turbulence. Finally, at the onset of our results we argue that the available turbulent energy in an exploding C+O white dwarf is probably too low in order to make a deflagration to detonation transition possible.

We investigate the possibility of a delayed detonation in a type Ia supernova under the assumption that the transition to detonation is triggered by turbulence only. Our discussion is based on the Zeldovich mechanism and suggests that typical turbulent velocities present during the explosion are not strong enough to allow this transition to occur. Although we are able to show that in carbon-rich matter (e.g., $X(^{12}$C$) = 0.75$) the possibility of a deflagration to detonation transition (DDT) is enhanced, even in this case the turbulent velocities needed are larger than the expected value of $u'(L) \approx 10^7 {cm s}^{-1}$ on a length-scale of $L \approx 10^6$ cm. Thus we conclude that a DDT may not be a common event during a thermonuclear explosion of a Chandrasekhar-mass white dwarf.

A variety of arguments suggest that the most common form of gamma-ray bursts (GRBs), those longer than a few seconds, involve the formation of black holes in supernova-like events. Two kinds of ``collapsar'' models are discussed, those in which the black hole forms promptly - a second or so after iron core collapse - and those in which formation occurs later, following ``fallback'' over a period of minutes to hours. In most cases, extraction of energy from a rapidly accreting disk (and a rapidly rotating black hole) is achieved by magnetohydrodynamical processes, although neutrino-powered models remain viable in cases where the accretion rate is >0.05 solar masses per second. GRBs are but one observable phenomenon accompanying black hole birth and other possibilities are discussed, some of which (long, faint GRBs and soft x-ray transients) may await discovery. Since they all involve black holes of similar mass accreting one to several M\sun, collapsars have a nearly standard total energy, around 10**52 erg, but both the fraction of that energy ejected as highly relativistic matter and the distribution of that energy with angle can be highly variable. An explanation is presented why inferred GRB luminosity might correlate inversely with time scales and arguments are given against the production of ordinary GRBs by supergiant stars.

A deep all-sky imaging hard x-ray survey and wide-field monitor is needed to extend soft (ROSAT) and medium (ABRIXAS2) x-ray surveys into the 10-100 keV band (and beyond) at comparable sensitivity (~0.05 mCrab). This would enable discovery and study of >3000 obscured AGN, which probably dominate the hard x-ray background; detailed study of spectra and variability of accreting black holes and a census of BHs in the Galaxy; Gamma-ray bursts and associated massive star formation (PopIII) at very high redshift and Soft Gamma-ray Repeaters throughout the Local Group; and a full galactic survey for obscured supernova remnants. The Energetic X-ray Imaging Survey Telescope (EXIST) is a proposed array of 8 x 1m^2 coded aperture telescopes fixed on the International Space Station (ISS) with 160deg x 40deg field of view which images the full sky each 90 min orbit. EXIST has been included in the most recent NASA Strategic Plan as a candidate mission for the next decade. An overview of the science goals and mission concept is presented.

We consider the explosion of supernovae and the possible production of a variety of high energy transients by delayed black hole formation in massive stars endowed with rotation. Following the launch of a ``successful'' shock by the usual neutrino powered mechanism, the inner layers of the star move outwards, but lack adequate energy to eject all the matter exterior to the neutron star. Over a period of minutes to hours a variable amount of mass, about 0.1 to 5 solar masses falls back into the collapsed remnant, often turning it into a black hole and establishing an accretion disk. The accretion rate, about 0.001 to 0.01 solar masses per second, is inadequate to produce a jet mediated by neutrino annihilation, but similar to that invoked in magnetohydrodynamic (MHD) models for gamma-ray bursts (GRBs). We thus consider the effect of jets formed by ``fallback'' in stars that are already in the process of exploding. We justify a parameterization of the jet power as a constant times the mass accretion rate, $\epsilon \dot {\rm M} c^2$, and explore the consequences of $\epsilon$ = 0.001 and 0.01. In supergiants, shock breakout produces bright x-ray transients that might be a diagnostic of the model. Jets produced by fallback should be more frequent than those made by the prompt formation of a black hole and may power the most common form of gamma-ray transient in the universe, although not the most common form seen so far by BATSE. Those are still attributed to prompt black hole formation, but it may be that the diverse energies observed for GRBs so far reflect chiefly the variable collimation of the jet inside the star and a consequently highly variable fraction of relativistic ejecta. Indeed, these events may all have a common total energy near 10$^{52}$ erg.

A diverse range of phenomena is possible when a black hole experiences very rapid accretion from a disk due to the incomplete explosion of a massive presupernova star endowed with rotation. In the most extreme case, the outgoing shock fails promptly in a rotating helium star, a black hole and an accretion disk form, and a strong gamma-ray burst (GRB) results. However, there may also be more frequently realized cases where the black hole forms after a delay of from several tens of seconds to several hours as approximately 0.1 to 5 solar masses falls back into the collapsed remnant following a mildly successful supernova explosion. There, the same MHD mechanisms frequently invoked to produce GRBs would also produce jets in stars already in the process of exploding. The presupernova star could be a Wolf-Rayet star or a red or blue supergiant. Depending upon its initial pressure, the collimation of the jet may also vary since ``hot'' jets will tend to diverge and share their energy with the rest of the star. From these situations, one expects diverse outcomes ranging from GRBs with a large range of energies and durations, to asymmetric, energetic supernovae with weak GRBs. SN 1998bw may have been the explosion of a star in which fall back produced a black hole and a less collimated jet than in the case of prompt black hole formation.

The evolution of rotating stars with zero-age main sequence (ZAMS) masses in the range 8 to 25 M_sun is followed through all stages of stable evolution. The initial angular momentum is chosen such that the star's equatorial rotational velocity on the ZAMS ranges from zero to ~ 70 % of break-up. Redistribution of angular momentum and chemical species are then followed as a consequence of rotationally induced circulation and instablities. The effects of the centrifugal force on the stellar structure are included. Uncertain mixing efficiencies are gauged by observations. We find, as noted in previous work, that rotation increases the helium core masses and enriches the stellar envelopes with products of hydrogen burning. We determine, for the first time, the angular momentum distribution in typical presupernova stars along with their detailed chemical structure. Angular momentum loss due to (non-magnetic) stellar winds and the redistribution of angular momentum during core hydrogen burning are of crucial importance for the specific angular momentum of the core. Neglecting magnetic fields, we find angular momentum transport from the core to the envelope to be unimportant after core helium burning. We obtain specific angular momenta for the iron core and overlaying material of 1E16...1E17 erg s. These values are insensitive to the initial angular momentum. They are small enough to avoid triaxial deformations of the iron core before it collapses, but could lead to neutron stars which rotate close to break-up. They are also in the range required for the collapsar model of gamma-ray bursts. The apparent discrepancy with the measured rotation rates of young pulsars is discussed.

While many models have been proposed for GRBs, those currently favored are all based upon the formation of and/or rapid accretion into stellar mass black holes. We present population synthesis calculations of these models using a Monte Carlo approach in which the many uncertain parameters intrinsic to such calculations are varied. We estimate the event rate for each class of model as well as the propagation distance for those having significant delay between formation and burst production, i.e., double neutron star (DNS) mergers and black hole-neutron star (BH/NS) mergers. For reasonable assumptions regarding the many uncertainties in population synthesis, we calculate a daily event rate in the universe for i) merging neutron stars: ~100/day; ii) neutron-star black hole mergers: ~450/day; iii) collapsars: ~10,000/day; iv) helium star black hole mergers: ~1000/day; and v) white dwarf black hole mergers: ~20/day. The range of uncertainty in these numbers however, is very large, typically two to three orders of magnitude. These rates must additionally be multiplied by any relevant beaming factor and sampling fraction (if the entire universal set of models is not being observed). Depending upon the mass of the host galaxy, half of the DNS and BH/NS mergers will happen within 60kpc (for a Milky-Way massed galaxy) to 5Mpc (for a galaxy with negligible mass) from the galactic center. Because of the delay time, neutron star and black hole mergers will happen at a redshift 0.5 to 0.8 times that of the other classes of models. Information is still lacking regarding the hosts of short hard bursts, but we suggest that they are due to DNS and BH/NS mergers and thus will ultimately be determined to lie outside of galaxies and at a closer mean distance than long complex bursts (which we attribute to collapsars).

We explore the sensitivity of the nucleosynthesis of intermediate mass elements (28 < A < 80) in supernovae derived from massive stars to the nuclear reaction rates employed in the model. Two standard sources of reaction rate data (Woosley et al. 1978; and Thielemann et al. 1987) are employed in pairs of calculations that are otherwise identical. Both include as a common backbone the experimental reactions rates of Caughlan & Fowler (1988). Two stellar models are calculated for each of two main sequence masses: 15 and 25 solar masses. Each star is evolved from core hydrogen burning to a presupernova state carrying an appropriately large reaction network, then exploded using a piston near the edge of the iron core as described by Woosley & Weaver (1995). The final stellar yields from the models calculated with the two rate sets are compared and found to differ in most cases by less than a factor of two over the entire range of nuclei studied. Reasons for the major discrepancies are discussed in detail along with the physics underlying the two reaction rate sets employed. The nucleosynthesis results are relatively robust and less sensitive than might be expected to uncertainties in nuclear reaction rates, though they are sensitive to the stellar model employed.

The merger of compact binaries, especially black holes and neutron stars, is frequently invoked to explain gamma-ray bursts (GRB's). In this paper, we present three dimensional hydrodynamical simulations of the relatively neglected mergers of white dwarfs and black holes. During the merger, the white dwarf is tidally disrupted and sheared into an accretion disk. Nuclear reactions are followed and the energy release is negligible. Peak accretion rates are ~0.05 Msun/s (less for lower mass white dwarfs) lasting for approximately a minute. Many of the disk parameters can be explained by a simple analytic model which we derive and compare to our simulations. This model can be used to predict accretion rates for white dwarf and black hole (or neutron star) masses which are not simulated in this paper. Although the mergers studied here create disks with larger radii, and longer accretion times than those from the merger of double neutron stars, a larger fraction of the merging star's mass becomes part of the disk. Thus the merger of a white dwarf and a black hole could produce a long duration GRB. The event rate of these mergers may be as high as 1/Myr per galaxy.

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.

Using a two-dimensional hydrodynamics code (PROMETHEUS), we study the continued evolution of rotating massive helium stars whose iron core collapse does not produce a successful outgoing shock, but instead forms a black hole. We study the formation of a disk, the associated flow patterns, and the accretion rate for disk viscosity parameter, alpha ~ 0.001 and 0.1. For the standard 14 solar mass model the average accretion rate for 15 s is 0.07 solar masses per second and the total energy deposited along the rotational axes by neutrino annihilation is (1 - 14) x 10**51 erg, depending upon the evolution of the Kerr parameter and uncertain neutrino efficiencies. Simulated deposition of this energy in the polar regions results in strong relativistic outflow - jets beamed to about 1.5% of the sky. The jets remain highly focused, and are capable of penetrating the star in 5 - 10 s. After the jet breaks through the surface of the star, highly relativistic flow can commence. Because of the sensitivity of the mass ejection and jets to accretion rate, angular momentum, and disk viscosity, and the variation of observational consequences with viewing angle, a large range of outcomes is possible ranging from bright GRBs like GRB 971214 to faint GRB-supernovae like SN 1998bw. X-ray precursors are also possible as the jet first breaks out of the star. While only a small fraction of supernovae make GRBs, we predict that all GRBs longer than a few seconds will make supernovae similar to SN 1998bw. However, hard, energetic GRBs shorter than a few seconds will be difficult to make in this model.

Recently a Type Ic supernova, SN 1998bw, was discovered coincident with a gamma-ray burst, GRB 980425. The supernova had unusual radio, optical, and spectroscopic properties. We explore here models based upon helium stars in the range 9 - 14 Msun and carbon-oxygen stars 6 - 11 Msun which experience unusually energetic explosions (kinetic energy 0.5 - 2.8 x 10**52 erg). Bolometric light curves and multi-band photometry are calculated and compared favorably with observations. Both LTE and non-LTE spectra are also calculated for the model that agrees best with the light curve, a carbon-oxygen core of 6 Msun exploded with a kinetic energy of 2.2 x 10**52 erg. We also examine potential mechanisms for producing the observed gamma-ray burst (GRB) - shock break-out and relativistic shock deceleration in circumstellar material. For spherically symmetric models, both fail to produce a GRB of even the low luminosity inferred for GRB 980425. The high explosion energies required to understand the supernova are in contrast to what is expected for such massive stars and may indicate that a new sort of explosion has been identified, possibly the consequence of a collapsar (Woosley 1993, 1996). Indeed a more likely explanation for what was seen is a highly asymmetric explosion in which the GRB was produced by a relativistic jet, perhaps viewed obliquely, and only a fraction of the total stellar mass was ejected, the remainder accreting into a black hole. The ejected mass (but not the 56Ni mass), explosion energy, and velocities may then be smaller. Other Type Ic supernova - GRB coincidences are suggested.

We present a model for gamma-ray bursts (GRB's) in which a stellar mass black hole acquires a massive accretion disk by merging with the helium core of its red giant companion. The black hole enters the helium core after it, or its neutron star progenitor, first experiences a common envelope phase that carries it inwards through the hydrogen envelope. Accretion of the last several solar masses of helium occurs on a time scale of roughly a minute and provides a neutrino luminosity of approximately 10^51 - 10^52 erg/s. Neutrino annihilation, 0.01% to 0.1% efficient, along the rotational axis then gives a baryon loaded fireball of electron-positron pairs and radiation (about 10$^{50}$ erg total) whose beaming and relativistic interaction with circumstellar material makes the GRB (e.g., Rees & Meszaros 1992). The useful energy can be greatly increased if energy can be extracted from the rotational energy of the black hole by magnetic interaction with the disk. Such events should occur at a rate comparable to that of merging neutron stars and black hole neutron star pairs and may be responsible for long complex GRB's, but not short hard ones.

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.

The light curve of Supernova (SN) 1993J is calculated using two approaches to radiation transport as exemplified by the two computer codes, STELLA and EDDINGTON. Particular attention is paid to shock breakout and the photometry in the U, B, and V bands during the first 120 days. The hydrodynamical model, the explosion of a 13 Msun star which had lost most of its hydrogenic envelope to a companion, is the same in each calculation. The comparison elucidates differences between the approaches and also serves to validate the results of both. STELLA includes implicit hydrodynamics and is able to model supernova evolution at early times, before the expansion is homologous. STELLA also employs multi-group photonics and is able to follow the radiation as it decouples from the matter. EDDINGTON uses a different algorithm for integrating the transport equation, assumes homologous expansion, and uses a finer frequency resolution. Good agreement is achieved between the two codes only when compatible physical assumptions are made about the opacity. A new result for SN 1993J is a prediction of the continuum spectrum near the shock breakout (calculated by STELLA) which is superior to the results of other standard single energy group hydrocodes such as VISPHOT or TITAN. Based on the results of our independent codes, we discuss the uncertainties involved in the current time dependent models of supernova light curves.

Ten years later, astronomers are still puzzled by the stellar evolution that produced SN 1987A --- a blue supergiant. In single star models, the new OPAL opacities make blue solutions more difficult to achieve, though still possible for certain choices of convection physics. We also consider rotation, which has the desirable effect of producing large surface enhancements of nitrogen and helium, but the undesirable effect of increasing the helium-core mass at the expense of the envelope. The latter makes blue solutions more difficult. Still, we seek a model that occurs with high probability in the LMC and for which the time-scale for making the last transition from red to blue, $\sim$ 20,000 years, has a physical interpretation --- the Kelvin-Helmholtz time of the helium core. Single star models satisfy both criteria and might yet prove to be the correct explanation for Sk -69 202, provided new rotational or convection physics can simultaneously give a blue star and explain the ring structure. Some speculations on how this might be achieved are presented and some aspects of binary models briefly discussed.

During the first 20 seconds of its life, the enormous neutrino luminosity of a neutron star drives appreciable mass loss from its surface. Previous investigations have shown that this neutrino-driven wind could be the site where the r-process occurs. The nucleosynthesis is sensitive to four physical parameters characterizing the wind: its mass loss rate, the entropy per baryon, the electron fraction, and the dynamic time scale. Different authors, using numerical models for supernovae, have arrived at qualitatively different values for these key parameters. Here we derive their values analytically and test our analytic results by numerical calculations using an implicit hydrodynamic code. Employing our analytic and numerical methods, we also investigate how various factors can affect our results. The derived entropy typically falls short, by a factor of two to three, of the value required to produce a strong r-process. Various factors that might give a higher entropy or a more rapid expansion in the wind are discussed.

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.

The flame born in the deep interior of a white dwarf that becomes a Type Ia supernova is subject to several instabilities. We briefly review these instabilities and the corresponding flame acceleration. We discuss the conditions necessary for each of the currently proposed explosion mechanisms and the attendant uncertainties. A grid of critical masses for detonation in the range $10^7$ - $2 \times 10^9$ g cm$^{-3}$ is calculated and its sensitivity to composition explored. Prompt detonations are physically improbable and appear unlikely on observational grounds. Simple deflagrations require some means of boosting the flame speed beyond what currently exists in the literature. ``Active turbulent combustion'' and multi-point ignition are presented as two plausible ways of doing this. A deflagration that moves at the ``Sharp-Wheeler'' speed, $0.1 g_{\rm eff} t$, is calculated in one dimension and shows that a healthy explosion is possible in a simple deflagration if the front moves with the speed of the fastest floating bubbles. The relevance of the transition to the ``distributed burning regime'' is discussed for delayed detonations. No model emerges without difficulties, but detonation in the distributed regime is plausible, will produce intermediate mass elements, and warrants further study.

A series of two dimensional numerical simulations of explosive nuclear burning is presented for white dwarfs near the Chandraskhar mass. We assume that the burning begins as a slow deflagration front at or near the center of the star, and continues until the density in the burning regions has declined to about 10$^7$ g cm$^{-3}$, where the flame is essentially extinguished. We employ a novel numerical representation of the turbulent flame brush based upon ideas previously developed for modelling laboratory combustion and explore in some detail the sensitivity of the outcome to the manner in which burning is initiated. In particular, we simulate 1) a centrally ignited deflagration, 2) off-center ignition at a single ``point", and 3) simultaneous off-center ignition at five ``points". We find that the amount of $^{56}$Ni that is produced and other observable properties depend sensitively upon how the fuel is ignited.

The possibility of making stellar mass black holes in supernovae that otherwise produce viable Type II and Ib supernova explosions is discussed and estimates given of their number in the Milky Way Galaxy. Observational diagnostics of stellar mass black hole formation are reviewed. While the equation of state sets the critical mass, fall back during the explosion is an equally important (and uncertain) element in determining if a black hole is formed. SN 1987A may or may not harbor a black hole, but if the critical mass for neutron stars is 1.5 - 1.6 M\sun, as Brown and Bethe suggest, it probably does. Observations alone do not yet resolve the issue. Reasons for this state of ambiguity are discussed and suggestions given as to how gamma-ray and x-ray observations in the future might help.

The production of the radioactive isotopes $^{44}$Ti and $^{60}$Co in all types of supernovae is examined and compared to observational constraints including Galactic $\gamma$--ray surveys, measurements of the diffuse 511 keV radiation, $\gamma$--ray observations of Cas A, the late time light curve of SN 1987A, and isotopic anomalies found in silicon carbide grains in meteorites. The (revised) line flux from $^{44}$Ti decay in the Cas A supernova remnant reported by COMPTEL on the Compton Gamma-Ray Observatory is near the upper bound expected from our models. The necessary concurrent ejection of $^{56}$Ni would also imply that Cas A was a brighter supernova than previously thought unless extinction in the intervening matter was very large. Thus, if confirmed, the reported amount of $^{44}$Ti in Cas A provides very interesting constraints on both the supernova environment and its mechanism. The abundances of $^{44}$Ti and $^{60}$Co ejected by Type II supernovae are such that gamma-radiation from $^{44}$Ti decay SN 1987A could be detected by a future generation of gamma-ray telescopes and that the decay of $^{60}$Co might provide an interesting contribution to the late time light curve of SN 1987A and other Type II supernovae. To produce the solar $^{44}$Ca abundance and satisfy all the observational constraints, nature may prefer at least the occasional explosion of sub-Chandrasekhar mass white dwarfs as Type Ia supernovae. Depending on the escape fraction of positrons due to $^{56}$Co made in all kinds of Type Ia supernovae, a significant fraction of the steady state diffuse 511 keV emission may arise from the annihilation of positrons produced during the decay of $^{44}$Ti to $^{44}$Ca. The Ca and Ti isotopic anomalies in pre-solar grains confirm the production of $^{44}$Ti in supernovae and that

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

Using recently calculated yields for Type II supernovae, along with models for chemical evolution and the distribution of mass in the interstellar medium, the current abundances and spatial distributions of two key gamma-ray radioactivities, $^{26}$Al and $^{60}$Fe, are determined. The estimated steady state production rates are 2.0 $\pm$ 1.0 M\sun \ Myr$^{-1}$ for $^{26}$Al and 0.75 $\pm$ 0.4 M\sun \ Myr$^{-1}$ for $^{60}$Fe. This corresponds to 2.2 $\pm$ 1.1 M\sun \ of $^{26}$Al and 1.7 $\pm$ 0.9 M\sun \ of $^{60}$Fe in the present interstellar medium. Sources of uncertainty are discussed, one of the more important being the current rate of core collapse supernovae in the Galaxy. Our simple model gives three per century, but reasonable changes in the star formation rate could easily accommodate a core collapse rate one-half as large, and thus one-half the yields. When these stellar and chemical evolution results are mapped into a three dimensional model of the Galaxy, the calculated 1809 keV gamma-ray flux map is consistent with the {\it Compton Gamma Ray Observatory} observations of a steep decline in the flux outside a longitude of $\pm$ 50$^\circ$ from the Galactic center, and the slight flux enhancements observed in the vicinity of spiral arms. Other potential stellar sources of $^{26}$Al and $^{60}$Fe are mentioned, especially the possibility of $^{60}$Fe synthesis in Type Ia supernovae. Predictions for the $^{60}$Fe mass distribution, total mass, and flux map are given.

Using the output from a grid of 60 Type II supernova models (Woosley \& Weaver 1994) of varying mass (11 $\ltaprx$ M/M\sun $\ltaprx$ 40) and metallicity (0, 10$^{-4}$, 0.01, 0.1, and 1 Z\sol), the chemical evolution of 76 stable isotopes, from hydrogen to zinc, is calculated. The chemical evolution calculation employs a simple dynamical model for the Galaxy (infall with a 4 billion year $e$-folding time scale onto a exponential disk and 1/r$^2$ bulge), and standard evolution parameters, such as a Salpeter initial mass function and a quadratic Schmidt star formation rate. The theoretical results are compared in detail with observed stellar abundances in stars with metallicities in the range -3.0 $\ltaprx$ [Fe/H] $\ltaprx$ 0.0 dex. While our discussion focuses on the solar neighborhood where there are the most observations, the supernovae rates, an intrinsically Galactic quantity, are also discussed.