Gravitational waves from coalescing neutron stars encode information about nuclear matter at extreme densities, inaccessible by laboratory experiments. The late inspiral is influenced by the presence of tides, which depend on the neutron star equation of state. Neutron star mergers are expected to often produce rapidly-rotating remnant neutron stars that emit gravitational waves. These will provide clues to the extremely hot post-merger environment. This signature of nuclear matter in gravitational waves contains most information in the 2-4 kHz frequency band, which is outside of the most sensitive band of current detectors. We present the design concept and science case for a neutron star extreme matter observatory (NEMO): a gravitational-wave interferometer optimized to study nuclear physics with merging neutron stars. The concept uses high circulating laser power, quantum squeezing and a detector topology specifically designed to achieve the high-frequency sensitivity necessary to probe nuclear matter using gravitational waves. Above one kHz, the proposed strain sensitivity is comparable to full third-generation detectors at a fraction of the cost. Such sensitivity changes expected event rates for detection of post-merger remnants from approximately one per few decades with two A+ detectors to a few per year, and potentially allows for the first gravitational-wave observations of supernovae, isolated neutron stars, and other exotica.

To date no metal-free stars have been identified by direct observations. The most common method of constraining their properties is searching the spectra of the most metal-poor stars for the chemical elements created in the first stars and their supernova. In this approach, modelled supernova yields are compared to the observed abundance patterns in extremely metal-poor stars. The method typically only uses the abundance ratios, i.e., the yields are diluted to the observed level. Following the usual assumption of spherical symmetry we compute a simple lower limit of the mass a supernova can mix with and find that it is consistent with all published simulations of early chemical enrichment in the interstellar medium. For three different cases, we demonstrate that this dilution limit can change the conclusions from the abundance fitting. There is a large discrepancy between the dilution found in simulations of SN explosions in minihaloes and the dilution assumed in many abundance fits. Limiting the dilution can significantly alter the likelihood of which supernovae are possible progenitors of observed CEMP-no stars. In particular, some of the faint, very low-yield SNe, which have been suggested as models for the abundance pattern of SMSS0313-6708, cannot explain the measured metal abundances, as their predicted metal yields are too small by two orders of magnitude. Altogether, the new dilution model presented here emphasizes the need to better understand the mixing and dilution behaviour of aspherical SNe.

We present 3D full-sphere supernova simulations of non-rotating low-mass (~9 Msun) progenitors, covering the entire evolution from core collapse through bounce and shock revival, through shock breakout from the stellar surface, until fallback is completed several days later. We obtain low-energy explosions [~(0.5-1.0)x 10^{50} erg] of iron-core progenitors at the low-mass end of the core-collapse supernova (LMCCSN) domain and compare to a super-AGB (sAGB) progenitor with an oxygen-neon-magnesium core that collapses and explodes as electron-capture supernova (ECSN). The onset of the explosion in the LMCCSN models is modelled self-consistently using the Vertex-Prometheus code, whereas the ECSN explosion is modelled using parametric neutrino transport in the Prometheus-HOTB code, choosing different explosion energies in the range of previous self-consistent models. The sAGB and LMCCSN progenitors that share structural similarities have almost spherical explosions with little metal mixing into the hydrogen envelope. A LMCCSN with less 2nd dredge-up results in a highly asymmetric explosion. It shows efficient mixing and dramatic shock deceleration in the extended hydrogen envelope. Both properties allow fast nickel plumes to catch up with the shock, leading to extreme shock deformation and aspherical shock breakout. Fallback masses of <~5x10^{-3} Msun have no significant effects on the neutron star (NS) masses and kicks. The anisotropic fallback carries considerable angular momentum, however, and determines the spin of the newly-born NS. The LMCCSNe model with less 2nd dredge-up results in a hydrodynamic and neutrino-induced NS kick of >40 km/s and a NS spin period of ~30 ms, both not largely different from those of the Crab pulsar at birth.

We calculate the mean evolution of the iron-peak abundance ratios [(Cr,Mn,Co,Zn)/Fe] in the Galaxy, using modern supernova and hypernova chemical yields and a Galactic Chemical Evolution code that assumes homogeneous chemical evolution. We investigate a range of hypernova occurrence rates and are able to produce a chemical composition that is a reasonable fit to the observed values in metal-poor stars. This requires a hypernova occurence rate that is large (50%) in the early Universe, decreasing throughout evolution to a value that is within present day observational constraints (>~ 1%). A large hypernova occurence rate is beneficial to matching the high [Zn/Fe] observed in the most metal-poor stars, although including hypernovae with progenitor mass >= 60 solar masses is detrimental to matching the observed [(Mn,Co)/Fe] evolution at low [Fe/H]. A significant contribution from HNe seems to be critical for producing supersolar [(Co,Zn)/Fe] at low metallicity, though more work will need to be done in order to match the most extreme values. We also emphasise the need to update models for the enrichment sources at higher metallicity, as the satisfactory recovery of the solar values of [(Cr,Mn,Co,Zn)/Fe] still presents a challenge.

Fallback in core-collapse supernovae plays a crucial role in determining the properties of the compact remnants and of the ejecta composition. We perform three-dimensional simulations of mixing and fallback for selected non-rotating supernova models to study how explosion energy and asymmetries correlate with the remnant mass, remnant kick, and remnant spin. We find that the strongest kick and spin is imparted by partial fallback in an asymmetric explosion. Black hole (BH) kicks of several hundred $\mathrm{km}\,\mathrm{s}^{-1}$ and spin parameters of $\mathord{\sim}0.25$ can be obtained in this scenario. If the initial explosion energy barely exceeds the envelope binding energy, stronger fallback results, and the remnant kick and spin remain small. If the explosion energy is high with respect to the envelope binding energy, there is little fallback with a small effect on the remnant kick, but the spin-up by fallback can be substantial. For a non-rotating $12\,\mathrm{M}_\odot$ progenitor, we find that the neutron star (NS) is spun up to millisecond periods. The high specific angular momentum of the fallback material can also lead to disk formation around black holes. Fallback may thus be a pathway towards millisecond-magnetar or collapsar-type engines for hypernovae and gamma-ray bursts that does not require rapid progenitor rotation. Within our small set of simulations, none reproduced the peculiar layered fallback necessary to explain the metal-rich iron-poor composition of many carbon-enhanced metal-poor (CEMP) stars. Models with different explosion energy and different realisations of asymmetries may, however, be compatible with CEMP abundance patterns.

Understanding the evolution of the N/O ratio in the interstellar medium (ISM) of galaxies is essential if we are to complete our picture of the chemical evolution of galaxies at high redshift, since most observational calibrations of O/H implicitly depend upon the intrinsic N/O ratio. The observed N/O ratio, however, shows large scatter at low O/H, and is strongly dependent on galactic environment. We show that several heretofore unexplained features of the N/O distribution at low O/H can be explained by the N seen in metal-poor galaxies being mostly primary nitrogen that is returned to the ISM via pre-supernova winds from rapidly rotating massive stars ($M \gtrsim 10$ M$_\odot$, $v/v_{\rm crit} \gtrsim 0.4$). This mechanism naturally produces the observed N/O plateau at low O/H. We show that the large scatter in N/O at low O/H also arises naturally from variations in star-formation efficiency. By contrast, models in which the N and O come primarily from supernovae provide a very poor fit to the observed abundance distribution. We propose that the peculiar abundance patterns we observe at low O/H are a signature that dwarf galaxies retain little of their SN ejecta, leaving them with abundance patterns typical of winds.

The evolutionary paths taken by massive stars with $M \gtrsim 60 \, \mathrm{M}_\odot$ remain substantially uncertain. They begin their lives as main sequence (MS) O-stars. Depending on their masses, rotation rates, and metallicities, they can then encounter a wide range of evolutionary states with an equally broad set of possible surface compositions and spectral classifications. We present a new grid of calculations for the evolution of such stars that covers a broad range in mass, M/M$_\odot = 60$ to $150$, rotation rate, $v \, / \, v_{\rm crit} = 0$ to $0.6$, metallicity, $[\mathrm{Fe}/\mathrm{H}] = -4$ to $0$, and $\alpha$-element enhancement, $[\alpha/\mathrm{Fe}] = 0$ to $0.4$. We show that rotating stars undergo rotationally-induced dredge-up of nucleosynthetic products, mostly He and N, to their surfaces while still on the MS. Non-rotating metal-rich stars also reveal the products of nucleosynthesis on their surfaces because even modest amounts of mass loss expose their "fossil" convective cores: regions that are no longer convective, but which were part of the convective core at an early stage in the star's evolution. Thus surface enhancement of He and N is expected for rotating stars at all metallicities, and for non-rotating stars if they are relatively metal-rich. We calculate a stellar atmosphere for a representative model from our grid, properly accounting for He- and N-enhancement, and show that the resulting spectrum provides a good match to observed WNL stars, strongly suggesting that the physical mechanisms we have identified are the ultimate cause of the WNL phase.

Type-I X-ray bursts are recurring thermonuclear explosions on the surface of accreting neutron stars. Matching observed bursts to computational models can help to constrain system properties, such as the neutron star mass and radius, crustal heating rates, and the accreted fuel composition, but systematic parameter studies to date have been limited. We apply Markov chain Monte Carlo methods to 1D burst models for the first time, and obtain system parameter estimations for the `Clocked Burster', GS 1826$-$238, by fitting multiple observed epochs simultaneously. We explore multiple parameters which are often held constant, including the neutron star mass, crustal heating rate, and hydrogen composition. To improve the computational efficiency, we precompute a grid of 3840 KEPLER models - the largest set of 1D burst simulations to date - and by interpolating over the model grid, we can rapidly sample burst predictions. We obtain estimates for a CNO metallicity of $Z_\mathrm{CNO} = 0.010^{+0.005}_{-0.004}$, a hydrogen fraction of $X_0 = 0.74^{+0.02}_{-0.03}$, a distance of $d \sqrt{\xi_\mathrm{b}} = 6.5^{+0.4}_{-0.6}\, \mathrm{kpc}$, and a system inclination of $i = {69^{+2}_{-3}}^{\circ}$.

Comparison of theoretical line profiles to observations provides important tests for supernova explosion models. We study the shapes of radioactive decay lines predicted by current 3D core-collapse explosion simulations, and compare these to observations of SN 1987A and Cas A. Both the widths and shifts of decay lines vary by several thousand kilometers per second depending on viewing angle. The line profiles can be complex with multiple peaks. By combining observational constraints from 56Co decay lines, 44Ti decay lines, and Fe IR lines, we delineate a picture of the morphology of the explosive burning ashes in SN 1987A. For M_ZAMS=15-20 Msun progenitors exploding with ~1.5 *10^51 erg, ejecta structures suitable to reproduce the observations involve a bulk asymmetry of the 56Ni of at least ~400 km/s and a bulk velocity of at least ~1500 km/s. By adding constraints to reproduce the UVOIR bolometric light curve of SN 1987A up to 600d, an ejecta mass around 14 Msun is favoured. We also investigate whether observed decay lines can constrain the neutron star (NS) kick velocity. The model grid provides a constraint V_NS > V_redshift, and applying this to SN 1987A gives a NS kick of at least 500 km/s. For Cas A, our single model provides a satisfactory fit to the NuSTAR observations and reinforces the result that current neutrino-driven core-collapse SN models can achieve enough bulk asymmetry in the explosive burning material. Finally, we investigate the internal gamma-ray field and energy deposition, and compare the 3D models to 1D approximations.

Supermassive stars have been proposed as the progenitors of the massive ($\sim 10^{9}\,\rm{M}_{\odot}$) quasars observed at $z\sim7$. Prospects for directly detecting supermassive stars with next-generation facilities depend critically on their intrinsic lifetimes, as well as their formation rates. We use the 1D stellar evolution code Kepler to explore the theoretical limiting case of zero-metallicity, non-rotating stars, formed monolithically with initial masses between $10\,\rm{kM}_{\odot}$ and $190\,\rm{kM}_{\odot}$. We find that stars born with masses between $\sim60\,\rm{kM}_{\odot}$ and $\sim150\,\rm{kM}_{\odot}$ collapse at the end of the main sequence, burning stably for $\sim1.5\,\rm{Myr}$. More massive stars collapse directly through the general relativistic instability after only a thermal timescale of $\sim3\,\rm{kyr}$--$4\,\rm{kyr}$. The expected difficulty in producing such massive, thermally-relaxed objects, together with recent results for currently preferred rapidly-accreting formation models, suggests that such ``truly direct'' or ``dark'' collapses may not be typical for supermassive objects in the early Universe. We close by discussing the evolution of supermassive stars in the broader context of massive primordial stellar evolution and the possibility of supermassive stellar explosions.

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.

A 70Msun BH was discovered in Milky Way disk in a long period and almost circular detached binary system (LB-1) with a high metallicity 8Msun B star companion. Current consensus on the formation of BHs from high metallicity stars limits the black hole mass to be below 20Msun. Using simple evolutionary model, we show that the formation of a 70Msun BH in high metallicity environment is possible if stellar wind mass loss rates are reduced by factor of 5. As observations indicate, a fraction of massive stars have surface magnetic fields which may quench the wind mass-loss, independently of stellar mass and metallicity. We also computed detailed stellar evolution models and we confirm such a scenario. A non-rotating 85Msun model at Z=0.014 with decreased winds ends up as a 71Msun star prior core-collapse with a 32Msun helium core and a 28Msun CO core. Such star avoids pair-instability pulsation supernova mass loss and may form a 70Msun BH in the direct collapse. Stars that can form such BHs expand to significant size with radius of R>600Rsun, exceeding the size of LB-1 orbit. Therefore, we can explain the formation of BHs upto 70Msun at high metallicity and this result is independent from LB-1. However, if LB-1 hosts a massive BH we are unable to explain how such a binary star system could have formed without invoking some exotic scenarios.

The isotopes $^{60}$Fe and $^{26}$Al originate from massive stars and their supernovae, reflecting ongoing nucleosynthesis in the Galaxy. We studied the gamma-ray emission from these isotopes at characteristic energies 1173, 1332, and 1809 keV with over 15 years of SPI data, finding a line flux in $^{60}$Fe combined lines of $(0.31\pm 0.06) \times 10^{-3}$ ph cm$^{-2}$ s$^{-1}$ and the $^{26}$Al line flux of $(16.8\pm 0.7) \times 10^{-4}$ ph cm$^{-2}$ s$^{-1}$ above the background and continuum emission for the whole sky. Based on the exponential-disk grid maps, we characterise the emission extent of $^{26}$Al to find scale parameters $R_0 =7.0^{+1.5}_{-1.0}$ kpc and $z_0=0.8^{+0.3}_{-0.2}$ kpc, however the $^{60}$Fe lines are too weak to spatially constrain the emission. Based on a point source model test across the Galactic plane, the $^{60}$Fe emission would not be consistent with a single strong point source in the Galactic center or somewhere else, providing a hint for a diffuse nature. We carried out comparisons of emission morphology maps using different candidate-source tracers for both $^{26}$Al and $^{60}$Fe emissions, and suggests that the $^{60}$Fe emission is more likely to be concentrated towards the Galactic plane. We determine the $^{60}$Fe/$^{26}$Al $\gamma$-ray flux ratio at $(18.4\pm4.2)\,\%$ , when using a parameterized spatial morphology model. Across the range of plausible morphologies, it appears possible that $^{26}$Al and $^{60}$Fe are distributed differently in the Galaxy. Using the best fitting maps for each of the elements, we constrain flux ratios in the range 0.2--0.4. We discuss its implications for massive star models and their nucleosynthesis.

We present a seven-minute long $4\pi$-3D simulation of a shell merger event in a non-rotating $18.88\, M_\odot$ supernova progenitor before the onset of gravitational collapse. The key motivation is to capture the large-scale mixing and asymmetries in the wake of the shell merger before collapse using a self-consistent approach. The $4\pi$ geometry is crucial as it allows us to follow the growth and evolution of convective modes on the largest possible scales. We find significant differences between the kinematic, thermodynamic and chemical evolution of the 3D and the 1D model. The 3D model shows vigorous convection leading to more efficient mixing of nuclear species. In the 3D case the entire oxygen shell attains convective Mach numbers of $\mathord{\approx}\, 0.1$, whereas in the 1D model, the convective velocities are much lower and there is negligible overshooting across convective boundaries. In the 3D case, the convective eddies entrain nuclear species from the neon (and carbon) layers into the deeper part of the oxygen burning shell, where they burn and power a violent convection phase with outflows. This is a prototypical model of a convective-reactive system. Due to the strong convection and the resulting efficient mixing, the interface between the neon layer and the silicon-enriched oxygen layer disappears during the evolution, and silicon is mixed far out into merged oxygen/neon shell. Neon entrained inwards by convective downdrafts burns, resulting in lower neon mass in the 3D model compared to the 1D model at time of collapse. In addition, the 3D model develops remarkable large-scale, large-amplitude asymmetries, which may have important implications for the impending gravitational collapse and the subsequent explosion.

Major mergers of gas-rich galaxies provide promising conditions for the formation of supermassive black holes (SMBHs; $\gtrsim10^5$ M$_\odot$) by direct collapse because they can trigger mass inflows as high as $10^4-10^5$ M$_\odot$ yr$^{-1}$ on sub-parsec scales. However, the channel of SMBH formation in this case, either dark collapse (direct collapse without prior stellar phase) or supermassive star (SMS; $\gtrsim10^4$ M$_\odot$), remains unknown. Here, we investigate the limit in accretion rate up to which stars can maintain hydrostatic equilibrium. We compute hydrostatic models of SMSs accreting at $1-1000$ M$_\odot$ yr$^{-1}$, and estimate the departures from equilibrium a posteriori by taking into account the finite speed of sound. We find that stars accreting above the atomic cooling limit ($\gtrsim10$ M$_\odot$ yr$^{-1}$) can only maintain hydrostatic equilibrium once they are supermassive. In this case, they evolve adiabatically with a hylotropic structure, that is, entropy is locally conserved and scales with the square root of the mass coordinate. Our results imply that stars can only become supermassive by accretion at the rates of atomically cooled haloes ($\sim0.1-10$ M$_\odot$ yr$^{-1}$). Once they are supermassive, larger rates are possible.

All stellar mass black holes have hitherto been identified by X-rays emitted by gas that is accreting onto the black hole from a companion star. These systems are all binaries with black holes below 30 M$_{\odot}$$^{1-4}$. Theory predicts, however, that X-ray emitting systems form a minority of the total population of star-black hole binaries$^{5,6}$. When the black hole is not accreting gas, it can be found through radial velocity measurements of the motion of the companion star. Here we report radial velocity measurements of a Galactic star, LB-1, which is a B-type star, taken over two years. We find that the motion of the B-star and an accompanying H$\alpha$ emission line require the presence of a dark companion with a mass of $68^{+11}_{-13}$ M$_{\odot}$, which can only be a black hole. The long orbital period of 78.9 days shows that this is a wide binary system. The gravitational wave experiments have detected similarly massive black holes$^{7,8}$, but forming such massive ones in a high-metallicity environment would be extremely challenging to current stellar evolution theories$^{9-11}$.

We report a new mechanism for the \textsl{s} in rotating massive metal-poor stars. Our models show that above a critical rotation speed, such stars evolve in a quasi-chemically-homogeneous fashion, which gives rise to a prolific \textsl{s}-process. Rotation-induced mixing results in primary production of $^{13}$C, which subsequently makes neutrons via $^{13}\mathrm{C}(\alpha,\mathrm{n})^{16}\mathrm{O}$ during core He burning. Neutron capture can last up to $\sim 10^{13}\,\mathrm{s}$ ($\sim 3\times 10^{5}$~yr) with the peak central neutron density ranging from $\sim10^7$ to $10^{8}\,\mathrm{cm}^{-3}$. Depending on the rotation speed and the mass loss rate, a strong \textsl{s} can occur with production of elements up to Bi for progenitors with initial metallicities of $[Z]\lesssim -1.5$. This result suggests that rapidly-rotating massive metal-poor stars are likely the first site for the main \textsl{s}-process. We find that these stars can potentially explain the early onset of the \textsl{s}-process and some of the carbon-enhanced metal-poor (CEMP-\textsl{s} and CEMP-\textsl{r/s}) stars with strong enrichment attributed to the \textsl{s} or a mixture of the \textsl{r}-process and the \textsl{s}-process.

We present a new method of matching observations of Type I (thermonuclear) X-ray bursts with models, comparing the predictions of a semi-analytic ignition model with X-ray observations of the accretion-powered millisecond pulsar SAX J1808.4--3658 in outburst. We used a Bayesian analysis approach to marginalise over the parameters of interest and determine parameters such as fuel composition, distance/anisotropy factors, neutron star mass and neutron star radius. Our study includes a treatment of the system inclination effects, inferring that the rotation axis of the system is inclined $\left(69^{+4}_{-2}\right)^\circ$ from the observers line of sight, assuming a flat disc model. This method can be applied to any accreting source that exhibits Type I X-ray bursts. We find a hydrogen mass fraction of $0.57^{+0.13}_{-0.14}$ and CNO metallicity of $0.013^{+0.006}_{-0.004}$ for the accreted fuel is required by the model to match the observed burst energies, for a distance to the source of $3.3^{+0.3}_{-0.2}\,\mathrm{kpc}$. We infer a neutron star mass of $1.5^{+0.6}_{-0.3}\,\mathrm{M}_{\odot}$ and radius of $11.8^{+1.3}_{-0.9}\,\mathrm{km}$ for a surface gravity of $1.9^{+0.7}_{-0.4}\times10^{14}\,\mathrm{cm}\,\mathrm{s}^{-2}$ for SAX J1808.4--3658.

Many distinct classes of high-energy variability have been observed in astrophysical sources, on a range of timescales. The widest range (spanning microseconds-decades) is found in accreting, stellar-mass compact objects, including neutron stars and black holes. Neutron stars are of particular observational interest, as they exhibit surface effects giving rise to phenomena (thermonuclear bursts and pulsations) not seen in black holes. Here we briefly review the present understanding of thermonuclear (type-I) X-ray bursts. These events are powered by an extensive chain of nuclear reactions, which are in many cases unique to these environments. Thermonuclear bursts have been exploited over the last few years as an avenue to measure the neutron star mass and radius, although the contribution of systematic errors to these measurements remains contentious. We describe recent efforts to better match burst models to observations, with a view to resolving some of the astrophysical uncertainties related to these events. These efforts have good prospects for providing complementary information to nuclear experiments.

During the first few hundred days after the explosion, core-collapse supernovae (SNe) emit down-scattered X-rays and gamma-rays originating from radioactive line emissions, primarily from the $^{56}$Ni $\rightarrow$ $^{56}$Co $\rightarrow$ $^{56}$Fe chain. We use SN models based on three-dimensional neutrino-driven explosion simulations of single stars and mergers to compute this emission and compare the predictions with observations of SN 1987A. A number of models are clearly excluded, showing that high-energy emission is a powerful way of discriminating between models. The best models are almost consistent with the observations, but differences that cannot be matched by a suitable choice of viewing angle are evident. Therefore, our self-consistent models suggest that neutrino-driven explosions are able to produce, in principle, sufficient mixing, although remaining discrepancies may require small changes to the progenitor structures. The soft X-ray cutoff is primarily determined by the metallicity of the progenitor envelope. The main effect of asymmetries is to vary the flux level by a factor of ${\sim}$3. For the more asymmetric models, the shapes of the light curves also change. In addition to the models of SN 1987A, we investigate two models of Type II-P SNe and one model of a stripped-envelope Type IIb SN. The Type II-P models have similar observables as the models of SN 1987A, but the stripped-envelope SN model is significantly more luminous and evolves faster. Finally, we make simple predictions for future observations of nearby SNe.

Supermassive primordial stars in hot, atomically-cooling haloes at $z \sim$ 15 - 20 may have given birth to the first quasars in the universe. Most simulations of these rapidly accreting stars suggest that they are red, cool hypergiants, but more recent models indicate that some may have been bluer and hotter, with surface temperatures of 20,000 - 40,000 K. These stars have spectral features that are quite distinct from those of cooler stars and may have different detection limits in the near infrared (NIR) today. Here, we present spectra and AB magnitudes for hot, blue supermassive primordial stars calculated with the TLUSTY and CLOUDY codes. We find that photometric detections of these stars by the James Webb Space Telescope (JWST) will be limited to $z \lesssim$ 10 - 12, lower redshifts than those at which red stars can be found, because of quenching by their accretion envelopes. With moderate gravitational lensing, Euclid and the Wide-Field Infrared Space Telescope (WFIRST) could detect blue supermassive stars out to similar redshifts in wide-field surveys.

The neutrino process that occurs in the outer stellar shells during a supernova explosion and involves neutrino-nucleus reactions produces a range of rare, stable and radioactive isotopes. We improve previous $\nu$-process studies by using, for the first time, the time-dependent neutrino emission spectra, as predicted from supernova simulations, rather than a simplified parametric description modeled after the neutron-star cooling phase. In particular, our calculations use time-dependent neutrino spectra for all neutrino species, consider their deviation from a Fermi-Dirac distribution and account for the neutrino emission from the neutrino burst and accretion phases. We find that the time-dependent treatment of the neutrino emission spectra results in higher yields for the selected nuclei produced by the $\nu$~process as compared to previous studies and also compared to the approximation of assuming constant neutrino energies corresponding to the time-averaged mean energy radiated in each species. The effect is largest for nuclides produced by charged-current reactions. Our results reflect the dynamical competition between neutrino-induced reactions and the effect of the shock passage through the star. By varying the neutrino burst luminosity and the duration of the accretion phase, we study the impact of these early emission phases and their uncertainties on the $\nu$-process nucleosynthesis. We find that the deviation of the neutrino spectra from a Fermi-Dirac distribution calculated in supernova simulations has a negligible effect on the $\nu$-process yields.

We assess the potential of using presupernova neutrino signals at the Jiangmen Underground Neutrino Observatory (JUNO) to probe the yet-unknown neutrino mass hierarchy. Using models for stars of 12, 15, 20, and 25 solar masses, we find that if the electron antineutrino signals from such a star can be predicted precisely and the star is within ~440-880 pc, the number of events of electron antineutrino captures on protons detected within one day of its explosion allows to determine the hierarchy at the > ~95% confidence level. For determination at this level using such signals from Betelgeuse, which is at a distance of ~222 pc, the uncertainty in the predicted number of signals needs to be < ~14-30%. In view of more realistic uncertainties, we discuss and advocate a model-independent determination using both electron neutrino and antineutrino signals from Betelgeuse. This method is feasible if the cosmogenic background for neutrino-electron scattering events can be reduced by a factor of ~2.5-10 from the current estimate. Such reduction might be achieved by using coincidence of the background events, the exploration of which for JUNO is highly desirable.

We present synthetic light curves of fallback-powered supernovae based on a neutrino-driven explosion of a 40 Msun zero-metallicity star with significant fallback accretion onto a black hole that was previously simulated by Chan et al. (2018) until shock breakout. Here, we investigate the light curve properties of the explosion after shock breakout for various fallback models. Without extra power from fallback accretion, the light curve is that of a Type IIP supernova with a plateau magnitude of around -14 mag and a plateau duration of 40 days. With extra power for the light curve from fallback accretion, however, we find that the transient can be significantly more luminous. The light-curve shape can be SN 1987A-like or Type IIP-like, depending on the efficiency of the fallback engine. If the accretion disk forms soon after the collapse and more than 1% of the accretion energy is released as the central engine, fallback accretion powered supernovae become as luminous as superluminous supernovae. We suggest that Type II superluminous supernovae with broad hydrogen features could be related to such hydrogen-rich supernovae powered by fallback accretion. In the future, such superluminous supernovae powered by fallback accretion might be found among the supernovae from the first stars in addition to pair-instability supernovae and pulstational pair-instability supernovae.

In recent years, the discovery of massive quasars at z~7 has provided a striking challenge to our understanding of the origin and growth of supermassive black holes in the early Universe. Mounting observational and theoretical evidence indicates the viability of massive seeds, formed by the collapse of supermassive stars, as a progenitor model for such early, massive accreting black holes. Although considerable progress has been made in our theoretical understanding, many questions remain regarding how (and how often) such objects may form, how they live and die, and how next generation observatories may yield new insight into the origin of these primordial titans. This review focusses on our present understanding of this remarkable formation scenario, based on discussions held at the Monash Prato Centre from November 20--24, 2017, during the workshop "Titans of the Early Universe: The Origin of the First Supermassive Black Holes."

The recent detection of 60Fe in the cosmic rays provides conclusive evidence that there is a recently synthesized component (few MY) in the GCRs (Binns et al. 2016). In addition, these nuclei must have been synthesized and accelerated in supernovae near the solar system, probably in the Sco-Cen OB association subgroups, which are about 100 pc distant from the Sun. Recent theoretical work on the production of r-process nuclei appears to indicate that it is difficult for SNe to produce the solar system abundances relative to iron of r-process elements with high atomic number (Z), including the actinides (Th, U, Np, Pu, and Cm). Instead, it is believed by many that the heaviest r-process nuclei, or perhaps even all r-process nuclei, are produced in binary neutron star mergers. Since we now know that there is at least a component of the GCRs that has been recently synthesized and accelerated, models of r-process production by SNe and BNSM can be tested by measuring the relative abundances of these ultra-heavy r-process nuclei, and especially the actinides, since they are radioactive and provide clocks that give the time interval from nucleosynthesis to detection at Earth. Since BNSM are believed to be much less frequent in our galaxy than SNe (roughly 1000 times less frequent, the ratios of the actinides, each with their own half-life, will enable a clear determination of whether the heaviest r-process nuclei are synthesized in SNe or in BNSM. In addition, the r-process nuclei for the charge range from 34 to 82 can be used to constrain models of r-process production in BNSM and SNe. Thus, GCRs become a multi-messenger component in the study of BNSM and SNe.

A major scientific goal of JWST is to probe the epoch of re-ionization of the Universe at z above 6, and up to 20 and beyond. At these redshifts, galaxies are just beginning to form and the observable objects are early black holes, supernovae, and cosmic infrared background. The JWST has the necessary sensitivity to observe these targets individually, but a public deep and wide science enabling survey in the wavelength range from 2-5 $\mu$m is needed to discover these black holes and supernovae and to cover the area large enough for cosmic infrared background to be reliably studied. This enabling survey will also discover a large number of other transients and enable sciences such as supernova cosmology up to z $\sim$ 5, star formation history at high redshift through supernova explosions, faint stellar objects in the Milky Way, and galaxy evolution up to z approaching 10. The results of this survey will also serve as an invaluable target feeder for the upcoming era of ELT and SKA.

There is now solid experimental evidence of at least one supernova explosion within 100 pc of Earth within the last few million years, from measurements of the short-lived isotope 60Fe in widespread deep-ocean samples, as well as in the lunar regolith and cosmic rays. This is the first established example of a specific dated astrophysical event outside the Solar System having a measurable impact on the Earth, offering new probes of stellar evolution, nuclear astrophysics, the astrophysics of the solar neighborhood, cosmic-ray sources and acceleration, multi-messenger astronomy, and astrobiology. Interdisciplinary connections reach broadly to include heliophysics, geology, and evolutionary biology. Objectives for the future include pinning down the nature and location of the established near-Earth supernova explosions, seeking evidence for others, and searching for other short-lived isotopes such as 26Al and 244Pu. The unique information provided by geological and lunar detections of radioactive 60Fe to assess nearby supernova explosions make now a compelling time for the astronomy community to advocate for supporting multi-disciplinary, cross-cutting research programs.

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.

Type I X-ray bursts are thermonuclear explosions on the surface of accreting neutron stars. Hydrogen rich X-ray bursts burn protons far from the line of stability and can release energy in the form of neutrinos from $\beta$-decays. We have estimated, for the first time, the neutrino fluxes of Type I bursts for a range of initial conditions based on the predictions of a 1D implicit hydrodynamics code, KEPLER, which calculates the complete nuclear reaction network. We find that neutrino losses are between $6.7 \times 10^{-5}$ and $0.14$ of the total energy per nucleon, Q$_{nuc}$, depending upon the hydrogen fraction in the fuel. These values are significantly below the $35\,\%$ value for neutrino losses often adopted in recent literature for the rp-process. The discrepancy arises because it is only at $\beta$-decays that $\approx35\,\%$ of energy is lost due to neutrino emission, whereas there are no neutrino losses in $(p,\gamma)$ and $(\alpha,p)$ reactions. Using the total measured burst energies from KEPLER for a range of initial conditions, we have determined an approximation formula for the total energy per nucleon released during an X-ray burst, Q$_{nuc}$=1.31+6.95$\bar{X} - 1.92\bar{X}^2 $MeV/nucleon, where $\bar{X}$ is the average hydrogen mass fraction of the ignition column, with an RMS error of $0.052\,$Mev/nucleon. We provide a detailed analysis of the nuclear energy output of a burst and find an incomplete extraction of mass excess in the burst fuel, with $14\,\%$ of the mass excess in the fuel not being extracted.

The collapse of supermassive primordial stars in hot, atomically-cooled halos may have given birth to the first quasars at $z \sim$ 15 - 20. Recent numerical simulations of these rapidly accreting stars reveal that they are cool, red hypergiants shrouded by dense envelopes of pristine atomically-cooled gas at 6,000 - 8,000 K, with luminosities $L$ $\gtrsim$ 10$^{10}$ L$_{\odot}$. Could such luminous but cool objects be detected as the first stage of quasar formation in future near infrared (NIR) surveys? We have now calculated the spectra of supermassive primordial stars in their birth envelopes with the Cloudy code. We find that some of these stars will be visible to JWST at $z \lesssim$ 20 and that with modest gravitational lensing Euclid and WFIRST could detect them out to $z \sim$ 10 - 12. Rather than obscuring the star, its accretion envelope enhances its visibility in the NIR today by reprocessing its short-wavelength flux into photons that are just redward of the Lyman limit in the rest frame of the star.

In this White Paper we present the potential of the Enhanced X-ray Timing and Polarimetry (eXTP) mission for determining the nature of dense matter; neutron star cores host an extreme density regime which cannot be replicated in a terrestrial laboratory. The tightest statistical constraints on the dense matter equation of state will come from pulse profile modelling of accretion-powered pulsars, burst oscillation sources, and rotation-powered pulsars. Additional constraints will derive from spin measurements, burst spectra, and properties of the accretion flows in the vicinity of the neutron star. Under development 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.

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.

In this paper we describe the potential of the enhanced X-ray Timing and Polarimetry (eXTP) mission for studies related to accretion flows in the strong field gravity regime around both stellar-mass and supermassive black-holes. eXTP has the unique capability of using advanced 'spectral-timing-polarimetry' techniques to analyze the rapid variations with three orthogonal diagnostics of the flow and its geometry, yielding unprecedented insight into the inner accreting regions, the effects of strong field gravity on the material within them and the powerful outflows which are driven by the accretion process.

We present a suite of seven 3D supernova simulations of non-rotating low-mass progenitors using multi-group neutrino transport. Our simulations cover single star progenitors with zero-age main sequence masses between $9.6 M_\odot$ and $12.5 M_\odot$ and (ultra)stripped-envelope progenitors with initial helium core masses between $2.8 M_\odot$ and $3.5 M_\odot$. We find explosion energies between $0.1\,\mathrm{Bethe}$ and $0.4\,\mathrm{Bethe}$, which are still rising by the end of the simulations. Although less energetic than typical events, our models are compatible with observations of less energetic explosions of low-mass progenitors. In six of our models, the mass outflow rate already exceeds the accretion rate onto the proto-neutron star, and the mass and angular momentum of the compact remnant have closely approached their final value, barring the possibility of later fallback. While the proto-neutron star is still accelerated by the gravitational tug of the asymmetric ejecta, the acceleration can be extrapolated to obtain estimates for the final kick velocity. We obtain gravitational neutron star masses between $1.22 M_\odot$ and $1.44 M_\odot$, kick velocities between $11\, \mathrm{km}\, \mathrm{s}^{-1}$ and $695\, \mathrm{km}\, \mathrm{s}^{-1}$, and spin periods from $20\, \mathrm{ms}$ to $2.7\,\mathrm{s}$, which suggests that typical neutron star birth properties can be naturally obtained in the neutrino-driven paradigm. We find a loose correlation between the explosion energy and the kick velocity. There is no indication of spin-kick alignment, but a correlation between the kick velocity and the neutron star angular momentum, which needs to be investigated further as a potential point of tension between models and observations.

Based on the work of Menon & Heger (2017), we present the bolometric light curvesand spectra of the explosions of blue supergiant progenitors from binary mergers. We study SN 1987A and two other peculiar Type IIP supernovae: SN 1998A and SN 2006V. The progenitor models were produced using the stellar evolution code Keplerand then exploded using the 1D radiation hydrodynamic code Crab. The explosions of binary merger models exhibit an overall better fit to the light curve of SN 1987A than previous single star models, because of their lower helium-core masses, larger envelope masses, and smaller radii. The merger model that best matches the observational constraints of the progenitor of SN 1987A and the light curve is a model with a radius of 37 solar radii, an ejecta mass of 20.6 solar masses, an explosion energy of 1.7 Bethe, a nickel mass of 0.073 solar masses, and a nickel mixing velocity of 3,000 km/s. This model also works for SN 1998A and is comparable with earlier estimates from semi-analytic models. In the case of SN 2006V, however, a model with a radius of 150 solar radii and ejecta mass of 19.1 solar masses matches the light curve. These parameters are significantly higher than predictions from semi-analytic models for the progenitor of this supernova.

We propose a new neutron-capture site in early metal-poor and metal-free stars of $\sim 20$--$30 \,\mathrm{M}_\odot$ that results from proton ingestion in the He shell during late stages of the stars' lives. Most of the neutron capture occurs in the first $\lesssim 10^6\,$s following proton ingestion when $^{13}{\rm C}(\alpha,\mathrm{n})^{16}\mathrm{O}$ produces neutron densities typical of the intermediate neutron-capture process. This phase may be followed by another lasting $\gtrsim 10^7\,$s with $^{17}\mathrm{O}(\alpha,\mathrm{n})^{20}\mathrm{Ne}$ producing much lower neutron densities typical of the slow neutron-capture process. We explore the dependence of the proposed neutron-capture nucleosynthesis on the amount and time of proton ingestion, the initial metallicity, and the ensuing supernova shock. We obtain a range of heavy element abundance patterns including those attributed to the slow neutron-capture process or a combination of the slow and rapid neutron-capture processes. Our results can account for the observed ubiquity of heavy elements such as Sr and Ba in the early Galaxy and explain puzzling abundance patterns of these elements in at least some very metal-poor (VMP) stars including those of the carbon-enhanced varieties. In the latter case, the explanation by the single site proposed here differs from the existing paradigm that attributes various classes of VMP stars to enrichment by multiple different sites.

Observations suggest that the interstellar medium (ISM) might have been highly enriched in carbon at very early times. We explore nucleosynthesis in massive carbon-enhanced metal-poor (CEMP) stars of $12$--$40\,\mathrm{M}_\odot$ formed from such an ISM with ${\rm [Fe]}\leq-2$. We find substantial production of elements heavier than Fe, mostly up to Sr, by the weak \textsl{s}-process in stars with initial abundances of ${\rm [CNO]}\gtrsim-1.5$. Even heavier elements, up to Ba, can be produced for ${\rm [CNO]}\gtrsim-0.5$. The efficiency of this \textsl{s}-process is sensitive to the initial enhancement of C (or more generally, CNO) and mass of the star, with the yield increasing approximately linearly with the initial Fe abundance. The \textsl{s}-process in CEMP stars of $\gtrsim 20 \,\mathrm{M}_\odot$ with ${\rm [CNO]}\gtrsim -1.5$ can be an important source for heavy elements in the Early Galaxy.

We investigate the relationship between explosion energy and nucleosynthesis in Population III supernovae and provide nucleosynthetic results for the explosions of stars with progenitor masses of $15\,\mathrm{M}_\odot$, $20\,\mathrm{M}_\odot$, $30\,\mathrm{M}_\odot$, $40\,\mathrm{M}_\odot$, $60\,\mathrm{M}_\odot$, and $80\,\mathrm{M}_\odot$, and explosion energies between approximately $10^{50}$ erg and $10^{53}$ erg. We find that the typical abundance pattern observed in metal-poor stars are best matched by supernovae with progenitor mass in the range $15\,\mathrm{M}_\odot$ - $30\,\mathrm{M}_\odot$, and explosion energy of $\sim$ ($5$ - $10$) $\times$ $10^{51}$ erg. In these models, a reverse shock caused by jumps in density between shells of different composition serves to decrease synthesis of chromium and manganese, which is favourable to matching the observed abundances in metal-poor stars. Spherically symmetric explosions of our models with progenitor mass $\geq$ $40\,\mathrm{M}_\odot$ do not provide yields that are compatible with the iron-peak abundances that are typically observed in metal-poor stars, however, by approximating the yields that we might expect from these models in highly aspherical explosions, we find indications that explosions of stars $40\,\mathrm{M}_\odot$ - $80\,\mathrm{M}_\odot$ with bipolar jets may be good candidates for the enrichment sources of metal-poor stars with enhanced carbon abundances.

We study the neutrino-induced production of nuclides in explosive supernova nucleosynthesis for progenitor stars with solar metallicity including neutrino nucleus reactions for all nuclei with charge numbers $Z < 76$ with average neutrino energies in agreement with modern Supernova simulations. Considering progenitors with initial main sequence masses between 13~M$_\odot$ and 30~M$_\odot$, we find a significant production of $^{11}$B, $^{138}$La, and $^{180}$Ta by neutrino nucleosynthesis, despite the significantly reduced neutrino energies. The production of $^{19}$F turns out to be more sensitive to the progenitor mass and structure than to the $\nu$ process. With our complete set of cross sections we have identified effects of the $\nu$~process on several stable nuclei including $^{33}$S, $^{40}$Ar, $^{41}$K, $^{59}$Co, and $^{113}$In at the 10\% level. Neutrino-induced reactions contribute to a similar extent to the production of radioactive ${}^{26}$Al and increase the yield of $^{22}$Na by 50\%. Future $\gamma$~ray astronomy missions may reach the precision at which the contribution from the $\nu$~process becomes relevant. We find that the production of $^{22}$Na by the $\nu$~process could explain the Ne-E(L) component of meteoritic graphite grains. The $\nu$~process enhances the yield of $^{36}$Cl and we point out that the resulting $^{36}$Cl/$^{35}$Cl ratio is in agreement with the values infrerred for the early solar system. Our extended set of neutrino-nucleus interactions also allows us to exclude any further effects of the $\nu$ process on stable nuclei and to quantify the effects on numerous, hitherto unconsidered radioactive nuclei, e.g., $^{36}$Cl, $^{72}$As, $^{84}$Rb, and $^{88}$Y.

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.

Modelling of thermonuclear X-ray bursts on accreting neutron stars has to date focused on stable accretion rates. However, bursts are also observed during episodes of transient accretion. During such events, the accretion rate can evolve significantly between bursts, and this regime provides a unique test for burst models. The accretion-powered millisecond pulsar SAX J1808.4-3658 exhibits accretion outbursts every 2-3 years. During the well-sampled month-long outburst of 2002 October, four helium-rich X-ray bursts were observed. Using this event as a test case, we present the first multi-zone simulations of X-ray bursts under a time-dependent accretion rate. We investigate the effect of using a time-dependent accretion rate in comparison to constant, averaged rates. Initial results suggest that using a constant, average accretion rate between bursts may underestimate the recurrence time when the accretion rate is decreasing, and overestimate it when the accretion rate is increasing. Our model, with an accreted hydrogen fraction of $X=0.44$ and a CNO metallicity of $Z_\mathrm{CNO}=0.02$, reproduces the observed burst arrival times and fluences with root mean square (RMS) errors of $2.8\,\mathrm{h}$, and $0.11\times 10^{-6}\,\mathrm{erg\, cm^{-2}}$, respectively. Our results support previous modelling that predicted two unobserved bursts, and indicate that additional bursts were also missed by observations.

The recent discoveries of many double neutron star systems and their detection as LIGO-Virgo merger events call for a detailed understanding of their origin. Explosions of ultra-stripped stars in binary systems have been shown to play a key role in this context and have also generated interest as a potential explanation for rapidly evolving hydrogen-free transients. Here we present the first attempt to model such explosions based on binary evolution calculations that follow the mass transfer to the companion to obtain a consistent core-envelope structure as needed for reliable predictions of the supernova transient. We simulate the explosion in 2D and 3D, and confirm the modest explosion energies ~10^50erg and small kick velocities reported earlier in 2D models based on bare carbon-oxygen cores. The spin-up of the neutron star by asymmetric accretion is small in 3D with no indication of spin-kick alignment. Simulations up to shock breakout show the mixing of sizeable amounts of iron group material into the helium envelope. In view of recent ideas for a mixing-length treatment (MLT) of Rayleigh-Taylor instabilities in supernovae, we perform a detailed analysis of the mixing, which reveals evidence for buoyancy-drag balance, but otherwise does not support the MLT approximation. The mixing may have implications for the spectroscopic signatures of ultra-stripped supernovae that need to be investigated in the future. Our stellar evolution calculation also predicts presupernova mass loss due to an off-centre silicon deflagration flash, which suggests that supernovae from extremely stripped cores may show signs of interactions with circumstellar material.

The production of the heavy stable proton-rich isotopes between 74Se and 196Hg -- the p nuclides -- is due to the contribution from different nucleosynthesis processes, activated in different types of stars. Whereas these processes have been subject to various studies, their relative contributions to Galactic Chemical Evolution are still a matter of debate. We investigate for the first time the nucleosynthesis of p nuclides in GCE by including metallicity and progenitor mass-dependent yields of core-collapse supernovae (ccSNe) into a chemical evolution model. We used a grid of metallicities and progenitor masses from two different sets of stellar yields and followed the contribution of ccSNe to the Galactic abundances as function of time. In combination with previous studies on p-nucleus production in thermonuclear supernovae (SNIa), and using the same GCE description, this allows us to compare the respective roles of SNeIa and ccSNe in the production of p-nuclei in the Galaxy. The gamma process in ccSN is very efficient for a wide range of progenitor masses (13M_sun_- 25M_sun) at solar metallicity. Since it is a secondary process with its efficiency depending on the initial abundance of heavy elements, its contribution is strongly reduced below solar metallicity. This makes it challenging to explain the inventory of the p nuclides in the Solar System by the contribution from ccSNe alone. In particular, we find that ccSNe contribute less than 10% of the solar p nuclide abundances, with only a few exceptions. Due to the uncertain contribution from other nucleosynthesis sites in ccSNe, such as neutrino winds or $\alpha$-rich freeze out, we conclude that the light p-nuclides 74$Se,78$Kr, 84Sr, and 92Mo may either still be completely or only partially produced in ccSNe.

Core-collapse supernova science is now entering an era where engine models are beginning to make both qualitative and quantitative predictions. Although the evidence in support of the convective engine for core-collapse supernova continues to grow, it is difficult to place quantitative constraints on this engine. Some studies have made specific predictions for the remnant distribution from the convective engine, but the results differ between different groups. Here we use a broad parameterization for the supernova engine to understand the differences between distinct studies. With this broader set of models, we place error bars on the remnant mass and basic yields from the uncertainties in the explosive engine. We find that, even with only 3 progenitors and a narrow range of explosion energies, we can produce a wide range of remnant masses and nucleosynthetic yields.

Long-duration gamma-ray bursts (LGRBs) and ultra-LGRBs (ULGRBs) originate from collapsars, in the center of which a newborn rotating stellar-mass black hole (BH) surrounded by a massive accretion disk may form. In the scenario of BH hyperaccretion inflow-outflow model and Blandford-Znajek (BZ) mechanism to trigger gamma-ray bursts (GRBs), the real accretion rate to power a BZ jet is far lower than the mass supply rate from the progenitor star. The characteristics of the progenitor stars can be constrained by GRB luminosity observations, and the results exceed usual expectations. LGRBs lasting from several seconds to tens of seconds in the rest frame may originate from solar-metallicity ($Z \sim 1~ Z_\odot$, where $Z$ and $Z_\odot$ are the metallicities of progenitor stars and the Sun), massive ($M \gtrsim 34 ~M_\odot$, where $M$ and $M_\odot$ are the masses of progenitor stars and the Sun) stars or some zero-metallicity ($Z \sim 0$) stars. A fraction of low-metallicity ($Z \lesssim 10^{-2}~Z_\odot$) stars, including Population III stars, can produce ULGRBs such as GRB 111209A. The fraction of LGRBs lasting less than tens of seconds in the rest frame is more than 40$\%$, which cannot conform to the fraction of the demanded type of progenitor star. It possibly implies that the activity timescale of central engine may be much longer than the observed timescale of prompt emission phase, as indicated by X-ray late-time activities. Alternatively, LGRBs and ULGRBs may be powered by a millisecond magnetar central engine.

Supermassive stars born from pristine gas in atomically-cooled haloes are thought to be the progenitors of supermassive black holes at high redshifts. However, the way they accrete their mass is still an unsolved problem. In particular, for accretion to proceed, a large amount of angular momentum has to be extracted from the collapsing gas. Here, we investigate the constraints stellar evolution imposes on this angular momentum problem. We present an evolution model of a supermassive Population III star including simultaneously accretion and rotation. We find that, for supermassive stars to form by accretion, the accreted angular momentum has to be about 1% of the Keplerian angular momentum. This tight constraint comes from the $\Omega\Gamma$-limit, at which the combination of radiation pressure and centrifugal force cancels gravity. It implies that supermassive stars are slow rotators, with a surface velocity less than 10-20% of their first critical velocity, at which the centrifugal force alone cancels gravity. At such low velocities, the deformation of the star due to rotation is negligible.

We examine the nucleosynthesis in the innermost, neutrino-processed ejecta (a few $10^{-3}\,M_\odot$) of self-consistent, two-dimensional explosion models of core-collapse supernovae for six progenitor stars with different initial masses. Three models have initial masses near the low-mass end of the supernova range, $8.8\,M_\odot$ (e8.8; electron-capture supernova), $9.6\,M_\odot$ (z9.6), and $8.1\,M_\odot$ (u8.1), with initial metallicities of 1, 0, and $10^{-4}$ times the solar metallicity, respectively. The other three are solar-metallicity models with initial masses of $11.2\,M_\odot$ (s11), $15\,M_\odot$ (s15), and $27\,M_\odot$ (s27). The low-mass models e8.8, z9.6, and u8.1 exhibit high production factors (nucleosynthetic abundances relative to the solar ones) of 100--200 for light trans-iron elements from Zn to Zr. This is associated with appreciable ejection of neutron-rich matter in these models. Remarkably, the nucleosynthetic outcomes for progenitors e8.8 and z9.6 are almost identical, including interesting productions of $^{48}$Ca and $^{60}$Fe, irrespective of their quite different (O-Ne-Mg and Fe) cores prior to collapse. In the more massive models s11, s15, and s27, several proton-rich isotopes of light trans-iron elements, including the $p$-isotope $^{92}$Mo (for s27) are made, up to production factors of $\sim$30. Both electron-capture and core-collapse supernovae near the low-mass end can therefore be dominant contributors to the Galactic inventory of light trans-iron elements from Zn to Zr and probably $^{48}$Ca and live $^{60}$Fe. The innermost ejecta of more massive supernovae may have only sub-dominant contributions to the chemical enrichment of the Galaxy except for $^{92}$Mo.

JWST was conceived and built to answer one of the most fundamental questions that humans can address empirically: "How did the Universe make its first stars?". Our First Lights At REionization (FLARE) project transforms the quest for the epoch of reionization from the static to the time domain. It targets the complementary question: "What happened to those first stars?". It will be answered by observations of the most luminous events: supernovae and accretion on to black holes formed by direct collapse from the primordial gas clouds. These transients provide direct constraints on star-formation rates (SFRs) and the truly initial Initial Mass Function (IMF), and they may identify possible stellar seeds of supermassive black holes (SMBHs). Furthermore, our knowledge of the physics of these events at ultra-low metallicity will be much expanded. JWST's unique capabilities will detect these most luminous and earliest cosmic messengers easily in fairly shallow observations. However, these events are very rare at the dawn of cosmic structure formation and so require large area coverage. Time domain astronomy can be advanced to an unprecedented depth by means of a shallow field of JWST reaching 27 mag AB in 2 and 4.4 microns over a field as large as 0.1 square degree visited multiple times each year. Such a survey may set strong constraints or detect massive Pop III SNe at redshifts beyond 10, pinpointing the redshift of the first stars, or at least their death. Based on our current knowledge of superluminous supernovae (SLSNe), such a survey will find one or more SLSNe at redshifts above 6 in five years and possibly several direct collapse black holes. Although JWST is not designed as a wide field survey telescope, we show that such a wide field survey is possible with JWST and is critical in addressing several of its key scientific goals.

The evolutionary pathways of core-collapse supernova progenitors at the low-mass end of the spectrum are beset with major uncertainties. In recent years, a variety of evolutionary channels has been discovered in addition to the classical electron capture supernova channel of super-AGB stars. The few available progenitor models at the low-mass end have been studied with great success in supernova simulations as the peculiar density structure makes for robust neutrino-driven explosions in this mass range. Detailed nucleosynthesis calculations have been conducted both for models of electron capture supernovae and low-mass iron core supernovae and revealed an interesting production of the lighter trans-iron elements (such as Zn, Sr, Y, Zr) as well as rare isotopes like Ca-48 and Fe-60. We stress the need to explore the low-mass end of the supernova spectrum further and link various observables to understand the diversity of explosions in this regime.

Fallback in core-collapse supernovae is considered a major ingredient for explaining abundance anomalies in metal-poor stars and the natal kicks and spins of black holes (BHs). We present a first 3D simulation of BH formation and fallback in an "aborted" neutrino-driven explosion of a $40$ solar mass zero-metallicity progenitor from collapse to shock breakout. We follow the phase up to BH formation using the relativistic CoCoNuT-FMT code. For the subsequent evolution to shock breakout we apply the moving-mesh code Arepo to core-collapse supernovae for the first time. Our simulation shows that despite early BH formation, neutrino-heated bubbles can survive for tens of seconds before being accreted, leaving them sufficient time to transfer part of their energy to sustain the shock wave as is propagates through the envelope. Although the initial net energy ($\sim 2$ Bethe) of the neutrino-heated ejecta barely equals the binding energy of the envelope, $11\,\mathrm{M}_\odot$ of hydrogen are still expelled with an energy of $0.23$ Bethe. We find no significant mixing and only a modest BH kick and spin, but speculate that stronger effects could occur for slightly more energetic explosions or progenitors with less tightly bound envelopes.

Recent three-dimensional simulations have suggested that convective seed perturbations from shell burning can play an important role in triggering neutrino-driven supernova explosions. Since isolated simulations cannot determine whether this perturbation-aided mechanism is of general relevance across the progenitor mass range, we here investigate the pertinent properties of convective oxygen and silicon burning shells in a broad range of presupernova stellar evolution models. We find that conditions for perturbation-aided explosions are most favourable in the extended oxygen shells of progenitors between about 16 and 26 solar masses, which exhibit large-scale convective overturn with high convective Mach numbers. Although the highest convective Mach numbers of up to 0.3 are reached in the oxygen shells of low-mass progenitors, convection is typically dominated by small-scale modes in these shells, which implies a more modest role of initial perturbations in the explosion mechanism. Convective silicon burning rarely provides the high Mach numbers and large-scale perturbations required for perturbation-aided explosions. We also find that about $40\%$ of progenitors between 16 and 26 solar masses exhibit simultaneous oxygen and neon burning in the same convection zone as a result of a shell merger shortly before collapse.

We study the impact of large-scale perturbations from convective shell burning on the core-collapse supernova explosion mechanism using three-dimensional (3D) multi-group neutrino hydrodynamics simulations of an 18 solar mass progenitor. Seed asphericities in the O shell, obtained from a recent 3D model of O shell burning, help trigger a neutrino-driven explosion 330ms after bounce whereas the shock is not revived in a model based on a spherically symmetric progenitor for at least another 300ms. We tentatively infer a reduction of the critical luminosity for shock revival by ~20% due to pre-collapse perturbations. This indicates that convective seed perturbations play an important role in the explosion mechanism in some progenitors. We follow the evolution of the 18 solar mass model into the explosion phase for more than 2s and find that the cycle of accretion and mass ejection is still ongoing at this stage. With a preliminary value of 0.77 Bethe for the diagnostic explosion energy, a baryonic neutron star mass of 1.85 solar masses, a neutron star kick of ~600km/s and a neutron star spin period of ~20ms at the end of the simulation, the explosion and remnant properties are slightly atypical, but still lie comfortably within the observed distribution. Although more refined simulations and a larger survey of progenitors are still called for, this suggests that a solution to the problem of shock revival and explosion energies in the ballpark of observations are within reach for neutrino-driven explosions in 3D.

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.

Supermassive primordial stars are suspected to be the progenitors of the most massive quasars at z~6. Previous studies of such stars were either unable to resolve hydrodynamical timescales or considered stars in isolation, not in the extreme accretion flows in which they actually form. Therefore, they could not self-consistently predict their final masses at collapse, or those of the resulting supermassive black hole seeds, but rather invoked comparison to simple polytropic models. Here, we systematically examine the birth, evolution and collapse of accreting non-rotating supermassive stars under accretion rates of 0.01-10 solar masses per year, using the stellar evolution code KEPLER. Our approach includes post-Newtonian corrections to the stellar structure and an adaptive nuclear network, and can transition to following the hydrodynamic evolution of supermassive stars after they encounter the general relativistic instability. We find that this instability triggers the collapse of the star at masses of 150,000-330,000 solar masses for accretion rates of 0.1-10 solar masses per year, and that the final mass of the star scales roughly logarithmically with the rate. The structure of the star, and thus its stability against collapse, is sensitive to the treatment of convection, and the heat content of the outer accreted envelope. Comparison with other codes suggests differences here may lead to small deviations in the evolutionary state of the star as a function of time, that worsen with accretion rate. Since the general relativistic instability leads to the immediate death of these stars, our models place an upper limit on the masses of the first quasars at birth.

Thermonuclear flashes of hydrogen and helium accreted onto neutron stars produce the frequently observed Type I X-ray bursts. It is the current paradigm that almost all material burns in a burst, after which it takes hours to accumulate fresh fuel for the next burst. In rare cases, however, bursts are observed with recurrence times as short as minutes. We present the first one-dimensional multi-zone simulations that reproduce this phenomenon. Bursts that ignite in a relatively hot neutron star envelope leave a substantial fraction of the fuel unburned at shallow depths. In the wake of the burst, convective mixing events driven by opacity bring this fuel down to the ignition depth on the observed timescale of minutes. There, unburned hydrogen mixes with the metal-rich ashes, igniting to produce a subsequent burst. We find burst pairs and triplets, similar to the observed instances. Our simulations reproduce the observed fraction of bursts with short waiting times of ~30%, and demonstrate that short recurrence time bursts are typically less bright and of shorter duration.

Supermassive primordial stars forming in atomically-cooled halos at $z \sim15-20$ are currently thought to be the progenitors of the earliest quasars in the Universe. In this picture, the star evolves under accretion rates of $0.1 - 1$ $M_\odot$ yr$^{-1}$ until the general relativistic instability triggers its collapse to a black hole at masses of $\sim10^5$ $M_\odot$. However, the ability of the accretion flow to sustain such high rates depends crucially on the photospheric properties of the accreting star, because its ionising radiation could reduce or even halt accretion. Here we present new models of supermassive Population III protostars accreting at rates $0.001 - 10$ $M_\odot$ yr$^{-1}$, computed with the GENEVA stellar evolution code including general relativistic corrections to the internal structure. We use the polytropic stability criterion to estimate the mass at which the collapse occurs, which has been shown to give a lower limit of the actual mass at collapse in recent hydrodynamic simulations. We find that at accretion rates higher than $0.001$ $M_\odot$ yr$^{-1}$ the stars evolve as red, cool supergiants with surface temperatures below $10^4$ K towards masses $>10^5$ $M_\odot$, and become blue and hot, with surface temperatures above $10^5$ K, only for rates $\lesssim0.001$ $M_\odot$ yr$^{-1}$. Compared to previous studies, our results extend the range of masses and accretion rates at which the ionising feedback remains weak, reinforcing the case for direct collapse as the origin of the first quasars.

We have used effective reaction rates (ERR) for the helium burning reactions to predict the yield of the gamma-emitting nuclei 26Al, 44Ti, and 60Fe in core col- lapse supernovae. The variations in the predicted yields for values of the reaction rates allowed by the ERR are much smaller than obtained previously, and smaller than other uncertainties. A "filter" for supernova nucleosynthesis yields based on pre-supernova structure was used to estimate the effect of failed supernovae on the initial mass function-averaged yields; this substantially reduced the yields of all these isotopes, but the predicted yield ratio 60Fe/26Al was little affected. The robustness of this ratio is promising for comparison with data, but it is larger than observed in nature; possible causes for this discrepancy are discussed.

Extremely metal-poor stars are uniquely informative on the nature of massive Population III stars. Modulo a few elements that vary with stellar evolution, the present-day photospheric abundances observed in extremely metal-poor stars are representative of their natal gas cloud composition. For this reason, the chemistry of extremely metal-poor stars closely reflects the nucleosynthetic yields of supernovae from massive Population III stars. Here we collate detailed abundances of 53 extremely metal-poor stars from the literature and infer the masses of their Population III progenitors. We fit a simple initial mass function to a subset of 29 of theinferred Population III star masses, and find that the mass distribution is well-represented by a power law IMF with exponent $\alpha = 2.35^{+0.29}_{-0.24}$. The inferred maximum progenitor mass for supernovae from massive Population III stars is $M_{\rm{max}} = 87^{+13}_{-33}$ M$_\odot$, and we find no evidence in our sample for a contribution from stars with masses above $\sim$120 M$_\odot$. The minimum mass is strongly consistent with the theoretical lower mass limit for Population III supernovae. We conclude that the IMF for massive Population III stars is consistent with the initial mass function of present-day massive stars and there may well have formed stars much below the supernova mass limit that could have survived to the present day.

We present the results of a detailed, systematic stellar evolution study of binary mergers for blue supergiant (BSG) progenitors of Type II supernovae. In particular, these are the first evolutionary models that can simultaneously reproduce nearly all observational aspects of the progenitor of SN 1987A, $\text{Sk}-69\,^{\circ}202$, such as its position in the HR diagram, the enrichment of helium and nitrogen in the triple-ring nebula, and its lifetime before its explosion. The merger model, based on the one proposed by Podsiadlowski 1992 et al. and Podsiadlowski 2007 et al., consists of a main sequence secondary star that dissolves completely in the common envelope of the primary red supergiant at the end of their merger. We empirically explore a large initial parameter space, such as primary masses ($15\,\text{M}_{\odot}$, $16\,\text{M}_{\odot}$, and $17\,\text{M}_{\odot}$), secondary masses ($2\,\text{M}_{\odot}$, $3\,\text{M}_{\odot}$, ..., $8\,\text{M}_{\odot}$) and different depths up to which the secondary penetrates the He core of the primary during the merger. The evolution of the merged star is continued until just before iron-core collapse and the surface properties of the 84 pre-supernova models ($16\,\text{M}_{\odot}-23\,\mathrm{M}_{\odot}$) computed have been made available in this work. Within the parameter space studied, the majority of the pre-supernova models are compact, hot BSGs with effective temperature $>12\,\text{kK}$ and radii of $30\,\text{R}_{\odot}-70\,\mathrm{R}_{\odot}$ of which six match nearly all the observational properties of $\text{Sk}-69\,^{\circ}202$.

About 4.6 billion years ago, some event disturbed a cloud of gas and dust, triggering the gravitational collapse that led to the formation of the solar system. A core-collapse supernova, whose shock wave is capable of compressing such a cloud, is an obvious candidate for the initiating event. This hypothesis can be tested because supernovae also produce telltale patterns of short-lived radionuclides, which would be preserved today as isotopic anomalies. Previous studies of the forensic evidence have been inconclusive, finding a pattern of isotopes differing from that produced in conventional supernova models. Here we argue that these difficulties either do not arise or are mitigated if the initiating supernova was a special type, low in mass and explosion energy. Key to our conclusion is the demonstration that short-lived Be-10 can be readily synthesized in such supernovae by neutrino interactions, while anomalies in stable isotopes are suppressed.

Modern numerical simulations of the formation of the first stars predict that the first stars formed in multiples. In those cases, the chemical yields of multiple supernova explosions may have contributed to the formation of a next generation star. We match the chemical abundances of the oldest observed stars in the universe to a database of theoretical supernova models, to show that it is likely that the first stars formed from the ashes of two or more progenitors.

We present the first 3D simulation of the last minutes of oxygen shell burning in an 18 solar mass supernova progenitor up to the onset of core collapse. A moving inner boundary is used to accurately model the contraction of the silicon and iron core according to a 1D stellar evolution model with a self-consistent treatment of core deleptonization and nuclear quasi-equilibrium. The simulation covers the full solid angle to allow the emergence of large-scale convective modes. Due to core contraction and the concomitant acceleration of nuclear burning, the convective Mach number increases to ~0.1 at collapse, and an l=2 mode emerges shortly before the end of the simulation. Aside from a growth of the oxygen shell from 0.51 to 0.56 solar masses due to entrainment from the carbon shell, the convective flow is reasonably well described by mixing length theory, and the dominant scales are compatible with estimates from linear stability analysis. We deduce that artificial changes in the physics, such as accelerated core contraction, can have precarious consequences for the state of convection at collapse. We argue that scaling laws for the convective velocities and eddy sizes furnish good estimates for the state of shell convection at collapse and develop a simple analytic theory for the impact of convective seed perturbations on shock revival in the ensuing supernova. We predict a reduction of the critical luminosity for explosion by 12--24% due to seed asphericities for our 3D progenitor model relative to the case without large seed perturbations.

Proton-rich nuclei, the so-called p-nuclei, are made in photodisintegration processes in outer shells of massive stars in the course of the final supernova explosion. Nuclear uncertainties in the production of these nuclei have been quantified in a Monte Carlo procedure. Bespoke temperature-dependent uncertainties were assigned to different types of reactions involving nuclei from Fe to Bi and all rates were varied randomly within the uncertainties. The resulting total production uncertainties of the p-nuclei are below a factor of two, with few exceptions. Key reactions dominating the final uncertainties have been identified in an automated procedure using correlations between rate and abundance uncertainties. Our results are compared to those of a previous study manually varying reaction rates.

With an increasing number of superluminous supernovae (SLSNe) discovered the question of their origin remains open and causes heated debates in the supernova community. Currently, there are three proposed mechanisms for SLSNe: (1) pair-instability supernovae (PISN), (2) magnetar-driven supernovae, and (3) models in which the supernova ejecta interacts with a circumstellar material ejected before the explosion. Based on current observations of SLSNe, the PISN origin has been disfavoured for a number of reasons. Many PISN models provide overly broad light curves and too reddened spectra, because of massive ejecta and a high amount of nickel. In the current study we re-examine PISN properties using progenitor models computed with the GENEC code. We calculate supernova explosions with FLASH and light curve evolution with the radiation hydrodynamics code STELLA. We find that high-mass models (200 and 250 solar masses) at relatively high metallicity (Z=0.001) do not retain hydrogen in the outer layers and produce relatively fast evolving PISNe Type I and might be suitable to explain some SLSNe. We also investigate uncertainties in light curve modelling due to codes, opacities, the nickel-bubble effect and progenitor structure and composition.

We report on the first high-resolution spectroscopic analysis of HE0020-1741, a bright (V=12.9), ultra metal-poor ([Fe/H] = -4.1), carbon-enhanced ([C/Fe] = +1.7) star selected from the Hamburg/ESO Survey. This star exhibits low abundances of neutron-capture elements ([Ba/Fe] = -1.1), and an absolute carbon abundance A(C) = 6.1; based on either criterion, HE0020-1741 is sub-classified as a CEMP-no star. We show that the light-element abundance pattern of HE0020-1741 is consistent with predicted yields from a massive (M = 21.5 Mo), primordial composition, supernova (SN) progenitor. We also compare the abundance patterns of other ultra metal-poor stars from the literature with available measures of C, N, Na, Mg, and Fe abundances with an extensive grid of SN models (covering the mass range 10 Mo - 100 Mo), in order to probe the nature of their likely stellar progenitors. Our results suggest that at least two classes of progenitors are required at [Fe/H] < -4.0, as the abundance patterns for more than half of the sample studied in this work (7 out of 12 stars) cannot be easily reproduced by the predicted yields.

Nuclear uncertainties in the production of $p$ nuclei in massive stars have been quantified in a Monte Carlo procedure. Bespoke temperature-dependent uncertainties were assigned to different types of reactions involving nuclei from Fe to Bi. Their simultaneous impact was studied in postprocessing explosive trajectories for three different stellar models. It was found that the grid of mass zones in the model of a 25 $M_\odot$ star, which is widely used for investigations of $p$ nucleosynthesis, is too crude to properly resolve the detailed temperature changes required for describing the production of $p$ nuclei. Using models with finer grids for 15 $M_\odot$ and 25 $M_\odot$ stars with initial solar metallicity, it was found that most of the production uncertainties introduced by nuclear reaction uncertainties are smaller than a factor of two. Since a large number of rates were varied at the same time in the Monte Carlo procedure, possible cancellation effects of several uncertainties could be taken into account. Key rates were identified for each $p$ nucleus, which provide the dominant contribution to the production uncertainty. These key rates were found by examining correlations between rate variations and resulting abundance changes. This method is superior to studying flow patterns, especially when the flows are complex, and to individual, sequential variation of a few rates.

We used a one-zone chemical evolution model to address the question of how many masses and metallicities are required in grids of massive stellar models in order to ensure reliable galactic chemical evolution predictions. We used a set of yields that includes seven masses between 13 and 30 Msun, 15 metallicities between 0 and 0.03 in mass fraction, and two different remnant mass prescriptions. We ran several simulations where we sampled subsets of stellar models to explore the impact of different grid resolutions. Stellar yields from low- and intermediate-mass stars and from Type Ia supernovae have been included in our simulations, but with a fixed grid resolution. We compared our results with the stellar abundances observed in the Milky Way for O, Na, Mg, Si, Ca, Ti, and Mn. Our results suggest that the range of metallicity considered is more important than the number of metallicities within that range, which only affects our numerical predictions by about 0.1 dex. We found that our predictions at [Fe/H] < -2 are very sensitive to the metallicity range and the mass sampling used for the lowest metallicity included in the set of yields. Variations between results can be as high as 0.8 dex, for any remnant mass prescription. At higher [Fe/H], we found that the required number of masses depends on the element of interest and on the remnant mass prescription. With a monotonic remnant mass prescription where every model explodes as a core-collapse supernova, the mass resolution induces variations of 0.2 dex on average. But with a remnant mass prescription that includes islands of non-explodability, the mass resolution can cause variations of about 0.2 to 0.7 dex depending on the choice of metallicity range. With such a prescription, explosive or non-explosive models can be missed if not enough masses are selected, resulting in over- or under-estimations of the mass ejected by massive stars.

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.

X-ray bursts are thermonuclear flashes on the surface of accreting neutron stars and reliable burst models are needed to interpret observations in terms of properties of the neutron star and the binary system. We investigate the dependence of X-ray burst models on uncertainties in (p,$\gamma$), ($\alpha$,$\gamma$), and ($\alpha$,p) nuclear reaction rates using fully self-consistent burst models that account for the feedbacks between changes in nuclear energy generation and changes in astrophysical conditions. A two-step approach first identified sensitive nuclear reaction rates in a single-zone model with ignition conditions chosen to match calculations with a state-of-the-art 1D multi-zone model based on the {\Kepler} stellar evolution code. All relevant reaction rates on neutron deficient isotopes up to mass 106 were individually varied by a factor of 100 up and down. Calculations of the 84 highest impact reaction rate changes were then repeated in the 1D multi-zone model. We find a number of uncertain reaction rates that affect predictions of light curves and burst ashes significantly. The results provide insights into the nuclear processes that shape X-ray burst observables and guidance for future nuclear physics work to reduce nuclear uncertainties in X-ray burst models.

We present a new approach to understand the landscape of supernova explosion energies, ejected nickel masses, and neutron star birth masses. In contrast to other recent parametric approaches, our model predicts the properties of neutrino-driven explosions based on the pre-collapse stellar structure without the need for hydrodynamic simulations. The model is based on physically motivated scaling laws and simple differential equations describing the shock propagation, the contraction of the neutron star, the neutrino emission, the heating conditions, and the explosion energetics. Using model parameters compatible with multi-D simulations and a fine grid of thousands of supernova progenitors, we obtain a variegated landscape of neutron star and black hole formation similar to other parameterised approaches and find good agreement with semi-empirical measures for the "explodability" of massive stars. Our predicted explosion properties largely conform to observed correlations between the nickel mass and explosion energy. Accounting for the coexistence of outflows and downflows during the explosion phase, we naturally obtain a positive correlation between explosion energy and ejecta mass. These correlations are relatively robust against parameter variations, but our results suggest that there is considerable leeway in parametric models to widen or narrow the mass ranges for black hole and neutron star formation and to scale explosion energies up or down. Our model is currently limited to an all-or-nothing treatment of fallback and there remain some minor discrepancies between model predictions and observational constraints.

We provide a set of stellar evolution and nucleosynthesis calculations that applies established physics assumptions simultaneously to low- and intermediate-mass and massive star models. Our goal is to provide an internally consistent and comprehensive nuclear production and yield data base for applications in areas such as pre-solar grain studies. Our non-rotating models assume convective boundary mixing where it has been adopted before. We include 8 (12) initial masses for $Z = 0.01$ ($0.02$). Models are followed either until the end of the asymptotic giant branch phase or the end of Si burning, complemented by a simple analytic core-collapse supernova models with two options for fallback and shock velocities. The explosions show which pre-supernova yields will most strongly be effected by the explosive nucleosynthesis. We discuss how these two explosion parameters impacts the light elements and the $s$ and $p$ process. For low- and intermediate-mass models our stellar yields from H to Bi include the effect of convective boundary mixing at the He-intershell boundaries and the stellar evolution feedback of the mixing process that produces the $^{13}$C pocket. All post-processing nucleosynthesis calculations use the same nuclear reaction rate network and nuclear physics input. We provide a discussion of the nuclear production across the entire mass range organized by element group. All our stellar nucleosynthesis profile and time evolution output is available electronically, and tools to explore the data on the NuGrid VOspace hosted by the Canadian Astronomical Data Centre are introduced.

Using the KEPLER hydrodynamics code, 464 models of thermonuclear X-ray bursters were performed across a range of accretion rates and compositions. We present the library of simulated burst profiles from this sample, and examine variations in the simulated lightcurve for different model conditions. We find that the recurrence time varies as a power law against accretion rate, and measure its slope while mixed H/He burning is occurring for a range of metallicities, finding the power law gradient to vary from $\eta = 1.1$ to $1.24$. We also identify the accretion rates at which mixed H/He burning stops and a transition occurs to different burning regimes. We also explore how varying the accretion rate and metallicity affects burst morphology in both the rise and tail.

Accreting neutron stars exhibit Type I X-ray bursts from both frequent hydrogen/helium flashes as well as rare carbon flashes. The latter (superbursts) ignite in the ashes of the former. Hydrogen/helium bursts, however, are thought to produce insufficient carbon to power superbursts. Stable burning could create the required carbon, but this was predicted to only occur at much larger accretion rates than where superbursts are observed. We present models of a new steady-state regime of stable hydrogen and helium burning that produces pure carbon ashes. Hot CNO burning of hydrogen heats the neutron star envelope and causes helium to burn before the conditions of a helium flash are reached. This takes place when the mass accretion rate is around 10% of the Eddington limit: close to the rate where most superbursts occur. We find that increased heating at the base of the envelope sustains steady-state burning by steepening the temperature profile, which increases the amount of helium that burns before a runaway can ensue.

Core-collapse supernovae are considered to be important contributors to the primitive dust enrichment of the interstellar medium in the high-redshift universe. Theoretical models of dust formation in stellar explosions have so far provided controversial results and a generally poor fit to the observations of dust formation in local supernovae. We present a new methodology for the calculation of carbonaceous dust formation in young supernova remnants. Our new technique uses both the nucleation theory and a chemical reaction network to allow us to compute the dust growth beyond the molecular level as well as to consider chemical erosion of the forming grains. We find that carbonaceous dust forms efficiently in the core of the ejecta, but takes several years to condensate, longer than previously estimated. It forms unevenly and remains concentrated in the inner part of the remnant. These results support the role of core-collapse supernovae as dust factories and provide new insight on the observations of SN 1987A, in which large amounts of dust have been detected to form on a timescale of years after core collapse.

We summarize our studies on neutrino-driven nucleosynthesis in He shells of early core-collapse supernovae with metallicities of $Z\lesssim 10^{-3}\,Z_\odot$. We find that for progenitors of $\sim 11$--$15\,{\rm M}_\odot$, the neutrons released by $^4{\rm He}(\bar{\nu}_e,e^+n)^3{\rm H}$ in He shells can be captured to produce nuclei with mass numbers up to $A \sim 200$. This mechanism is sensitive to neutrino emission spectra and flavor oscillations. In addition, we find two new primary mechanisms for neutrino-induced production of $^{9}$Be in He shells. The first mechanism produces $^9$Be via $^7{\rm Li}(n,\gamma)^8{\rm Li}(n,\gamma)^9{\rm Li}(e^-\bar{\nu}_e)^9{\rm Be}$ and relies on a low explosion energy for its survival. The second mechanism operates in progenitors of $\sim 8\,{\rm M}_\odot$, where $^9$Be can be produced directly via $^7{\rm Li}(^3{\rm H},n_0)^9{\rm Be}$ during the rapid expansion of the shocked He-shell material. The light nuclei $^7$Li and $^3$H involved in these mechanisms are produced by neutrino interactions with $^4$He. We discuss the implications of neutrino-induced nucleosynthesis in He shells for interpreting the elemental abundances in metal-poor stars.

If very massive stars (M >~ 100 Msun) can form and avoid too strong mass loss during their evolution, they are predicted to explode as pair-instability supernovae (PISNe). One critical test for candidate events is whether their nucleosynthesis yields and internal ejecta structure, being revealed through nebular-phase spectra at t >~ 1 yr, match those of model predictions. Here we compute theoretical spectra based on model PISN ejecta at 1-3 years post-explosion to allow quantitative comparison with observations. The high column densities of PISNe lead to complete gamma-ray trapping for t >~ 2 years which, combined with fulfilled conditions of steady state, leads to bolometric supernova luminosities matching the 56Co decay. Most of the gamma-rays are absorbed by the deep-lying iron and silicon/sulphur layers. The ionization balance shows a predominantly neutral gas state, which leads to emission lines of Fe I, Si I, and S I. For low-mass PISNe the metal core expands slowly enough to produce a forest of distinct lines, whereas high-mass PISNe expand faster and produce more featureless spectra. Line blocking is complete below ~5000 A for several years, and the model spectra are red. The strongest line is typically [Ca II] 7291,7323, one of few lines from ionized species. We compare our models with proposed PISN candidates SN 2007bi and PTF12dam, finding discrepancies for several key observables and thus no support for a PISN interpretation. We discuss distinct spectral features predicted by the models, and the possibility of detecting pair-instability explosions among non-superluminous supernovae.

The existence of black holes (BHs) of mass ~ 10^{9} M_sun at z > 6 is a big puzzle in astrophysics because even optimistic estimates of the accretion time are insufficient for stellar mass BHs of ~ 10 M_sun to grow into such supermassive BHs. A resolution of this puzzle might be the direct collapse of supermassive stars with mass M ~ 10^{5} M_sun into massive seed BHs. We find that if a jet is launched from the accretion disk around the central BH, the jet can break out the star because of the structure of the radiation pressure-dominated envelope. Such ultra-long gamma-ray bursts with duration of ~ 10^{4} - 10^{6} s and flux of 10^{-11} - 10^{-8} erg s^{-1} cm^{-2} could be detectable by Swift. We estimate an event rate of < 1yr^{-1}. The total explosion energy is > 10^{55} - 10^{56} erg. The resulting negative feedback delays the growth of the remnant BH by about 70 Myr or evacuates the host galaxy completely.

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.

We report the discovery of one extremely metal-poor (EMP; [Fe/H]<-3) and one ultra metal-poor (UMP; [Fe/H]<-4) star selected from the SDSS/SEGUE survey. These stars were identified as EMP candidates based on their medium-resolution (R~2,000) spectra, and were followed-up with high-resolution (R~35,000) spectroscopy with the Magellan-Clay Telescope. Their derived chemical abundances exhibit good agreement with those of stars with similar metallicities. We also provide new insights on the formation of the UMP stars, based on comparison with a new set of theoretical models of supernovae nucleosynthesis. The models were matched with 20 UMP stars found in the literature, together with one of the program stars (SDSS J1204+1201), with [Fe/H]=-4.34. From fitting their abundances, we find that the supernovae progenitors, for stars where carbon and nitrogen are measured, had masses ranging from 20.5 M_sun to 28 M_sun and explosion energies from 0.3 to 0.9x10^51 erg. These results are highly sensitive to the carbon and nitrogen abundance determinations, which is one of the main drivers for future high-resolution follow-up of UMP candidates. In addition, we are able to reproduce the different CNO abundance patterns found in UMP stars with a single progenitor type, by varying its mass and explosion energy.

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.

Measurements of explosive nucleosynthesis yields in core-collapse supernovae provide tests for explosion models. We investigate constraints on explosive conditions derivable from measured amounts of nickel and iron after radioactive decays using nucleosynthesis networks with parameterized thermodynamic trajectories. The Ni/Fe ratio is for most regimes dominated by the production ratio of 58Ni/(54Fe + 56Ni), which tends to grow with higher neutron excess and with higher entropy. For SN 2012ec, a supernova that produced a Ni/Fe ratio of $3.4\pm1.2$ times solar, we find that burning of a fuel with neutron excess $\eta \approx 6\times 10^{-3}$ is required. Unless the progenitor metallicity is over 5 times solar, the only layer in the progenitor with such a neutron excess is the silicon shell. Supernovae producing large amounts of stable nickel thus suggest that this deep-lying layer can be, at least partially, ejected in the explosion. We find that common spherically symmetric models of $M_{\rm ZAMS} \lesssim 13$ Msun stars exploding with a delay time of less than one second ($M_{\rm cut} < 1.5$ Msun) are able to achieve such silicon-shell ejection. Supernovae that produce solar or sub-solar Ni/Fe ratios, such as SN 1987A, must instead have burnt and ejected only oxygen-shell material, which allows a lower limit to the mass cut to be set. Finally, we find that the extreme Ni/Fe value of 60-75 times solar derived for the Crab cannot be reproduced by any realistic-entropy burning outside the iron core, and neutrino-neutronization obtained in electron-capture models remains the only viable explanation.

SMSS J031300.36-670839.3 (hereafter SM0313-6708) is a sub-giant halo star, with no detectable Fe lines and large overabundances of C and Mg relative to Ca. We obtained VLT-UVES spectra extending to 3060 Angstroms showing strong OH A-X band lines enabling an oxygen abundance to be derived. The OH A-X band lines in SM0313-6708 are much stronger than the CH C-X band lines. Spectrum synthesis fits indicate an [O/C] ratio of 0.02 +- 0.175. Our high S/N UVES data also enabled us to lower the Fe abundance limit to [Fe/H]{3D},NLTE < -7.52 (3 sigma). These data support our previous suggestion that the star formed from the iron-poor ejecta of a single massive star Population III supernova.

We study the neutrino-induced production of nuclides in explosive supernova nucleosynthesis for progenitor stars with solar metallicity and initial main sequence masses between 15 M$_\odot$ and 40 M$_\odot$. We improve previous investigations i) by using a global set of partial differential cross sections for neutrino-induced charged- and neutral-current reactions on nuclei with charge numbers $Z < 76 $ and ii) by considering modern supernova neutrino spectra which have substantially lower average energies compared to those previously adopted in neutrino nucleosynthesis studies. We confirm the production of $^7$Li, $^{11}$B, $^{138}$La, and $^{180}$Ta by neutrino nucleosynthesis, albeit at slightly smaller abundances due to the changed neutrino spectra. We find that for stars with a mass smaller than 20 M$_\odot$, $^{19}$F is produced mainly by explosive nucleosynthesis while for higher mass stars it is produced by the $\nu$ process. We also find that neutrino-induced reactions, either directly or indirectly by providing an enhanced abundance of light particles, noticeably contribute to the production of the radioactive nuclides $^{22}$Na and $^{26}$Al. Both nuclei are prime candidates for gamma-ray astronomy. Other prime targets, $^{44}$Ti and $^{60}$Fe, however, are insignificantly produced by neutrino-induced reactions. We also find a large increase in the production of the long-lived nuclei $^{92}$Nb and $^{98}$Tc due to charged-current neutrino capture.

The GALAH survey is a large high-resolution spectroscopic survey using the newly commissioned HERMES spectrograph on the Anglo-Australian Telescope. The HERMES spectrograph provides high-resolution (R ~28,000) spectra in four passbands for 392 stars simultaneously over a 2 degree field of view. The goal of the survey is to unravel the formation and evolutionary history of the Milky Way, using fossil remnants of ancient star formation events which have been disrupted and are now dispersed throughout the Galaxy. Chemical tagging seeks to identify such dispersed remnants solely from their common and unique chemical signatures; these groups are unidentifiable from their spatial, photometric or kinematic properties. To carry out chemical tagging, the GALAH survey will acquire spectra for a million stars down to V~14. The HERMES spectra of FGK stars contain absorption lines from 29 elements including light proton-capture elements, alpha-elements, odd-Z elements, iron-peak elements and n-capture elements from the light and heavy s-process and the r-process. This paper describes the motivation and planned execution of the GALAH survey, and presents some results on the first-light performance of HERMES.

The recent discovery that the Fe-K line luminosities and energy centroids observed in nearby SNRs are a strong discriminant of both progenitor type and circumstellar environment has implications for our understanding of supernova progenitor evolution. Using models for the chemical composition of core-collapse supernova ejecta, we model the dynamics and thermal X-ray emission from shocked ejecta and circumstellar material, modeled as an $r^{-2}$ wind, to ages of 3000 years. We compare the X-ray spectra expected from these models to observations made with the Suzaku satellite. We also model the dynamics and X-ray emission from Type Ia progenitor models. We find a clear distinction in Fe-K line energy centroid between core-collapse and Type Ia models. The core-collapse supernova models predict higher Fe-K line centroid energies than the Type Ia models, in agreement with observations. We argue that the higher line centroids are a consequence of the increased densities found in the circumstellar environment created by the expansion of the slow-moving wind from the massive progenitors.

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.

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.

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.

Massive stars are key sources of radiative, kinetic, and chemical feedback in the universe. Grids of massive star models computed by different groups each using their own codes, input physics choices and numerical approximations, however, lead to inconsistent results for the same stars. We use three of these 1D codes---GENEC, KEPLER and MESA---to compute non-rotating stellar models of $15~\mathrm{M}_\odot$, $20~\mathrm{M}_\odot$, and $25~\mathrm{M}_\odot$ and compare their nucleosynthesis. We follow the evolution from the main sequence until the end of core helium burning. The GENEC and KEPLER models hold physics assumptions used in large grids of published models. The MESA code was set up to use convective core overshooting such that the CO core masses are consistent with those obtained by GENEC. For all models, full nucleosynthesis is computed using the NuGrid post-processing tool MPPNP. We find that the surface abundances predicted by the models are in reasonable agreement. In the helium core, the standard deviation of the elemental overproduction factors for Fe to Mo is less than $30\,\%$---smaller than the impact of the present nuclear physics uncertainties. For our three initial masses, the three stellar evolution codes yield consistent results. Differences in key properties of the models, e.g., helium and CO core masses and the time spent as a red supergiant, are traced back to the treatment of convection and, to a lesser extent, mass loss. The mixing processes in stars remain the key uncertainty in stellar modelling. Better constrained prescriptions are thus necessary to improve the predictive power of stellar evolution models.

The discovery of 150 - 300 M$_{\odot}$ stars in the Local Group and pair-instability supernova candidates at low redshifts has excited interest in this exotic explosion mechanism. Realistic light curves for pair-instability supernovae at near-solar metallicities are key to identifying and properly interpreting these events as more are found. We have modeled pair-instability supernovae of 150 - 500 M$_{\odot}$ Z $\sim$ 0.1 - 0.4 Z$_{\odot}$ stars. These stars lose up to 80% of their mass to strong line-driven winds and explode as bare He cores. We find that their light curves and spectra are quite different from those of Population III pair-instability explosions, which therefore cannot be used as templates for low-redshift events. Although non-zero metallicity pair-instability supernovae are generally dimmer than their Population III counterparts, in some cases they will be bright enough to be detected at the earliest epochs at which they can occur, the formation of the first galaxies at $z \sim$ 10 - 15. Others can masquerade as dim, short duration supernovae that are only visible in the local universe and that under the right conditions could be hidden in a wide variety of supernova classes. We also report for the first time that some pair-instability explosions can create black holes with masses of $\sim$ 100 M$_{\odot}$.

The Large Observatory For x-ray Timing (LOFT) was studied within ESA M3 Cosmic Vision framework and participated in the final down-selection for a launch slot in 2022-2024. Thanks to the unprecedented combination of effective area and spectral resolution of its main instrument, LOFT will study the behaviour of matter under extreme conditions, such as the strong gravitational field in the innermost regions of accretion flows close to black holes and neutron stars, and the supra-nuclear densities in the interior of neutron stars. The science payload is based on a Large Area Detector (LAD, 10 m 2 effective area, 2-30 keV, 240 eV spectral resolution, 1 deg collimated field of view) and a WideField Monitor (WFM, 2-50 keV, 4 steradian field of view, 1 arcmin source location accuracy, 300 eV spectral resolution). The WFM is equipped with an on-board system for bright events (e.g. GRB) localization. The trigger time and position of these events are broadcast to the ground within 30 s from discovery. In this paper we present the status of the mission at the end of its Phase A study.

Among the short-lived radioactive nuclei inferred to be present in the early solar system via meteoritic analyses, there are several heavier than iron whose stellar origin has been poorly understood. In particular, the abundances inferred for 182Hf (half-life = 8.9 million years) and 129I (half-life = 15.7 million years) are in disagreement with each other if both nuclei are produced by the rapid neutron-capture process. Here, we demonstrate that contrary to previous assumption, the slow neutron-capture process in asymptotic giant branch stars produces 182Hf. This has allowed us to date the last rapid and slow neutron-capture events that contaminated the solar system material at roughly 100 million years and 30 million years, respectively, before the formation of the Sun.

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.

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

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 stability of thermonuclear burning of hydrogen and helium accreted onto neutron stars is strongly dependent on the mass accretion rate. The burning behavior is observed to change from Type I X-ray bursts to stable burning, with oscillatory burning occurring at the transition. Simulations predict the transition at a ten times higher mass accretion rate than observed. Using numerical models we investigate how the transition depends on the hydrogen, helium, and CNO mass fractions of the accreted material, as well as on the nuclear reaction rates of triple alpha and the hot-CNO breakout reactions 15O(a,g)19Ne and 18Ne(a,p)21Na. For a lower hydrogen content the transition is at higher accretion rates. Furthermore, most experimentally allowed reaction rate variations change the transition accretion rate by at most 10%. A factor ten decrease of the 15O(a,g)19Ne rate, however, produces an increase of the transition accretion rate of 35%. None of our models reproduce the transition at the observed rate, and depending on the true 15O(a,g)19Ne reaction rate, the actual discrepancy may be substantially larger. We find that the width of the interval of accretion rates with marginally stable burning depends strongly on both composition and reaction rates. Furthermore, close to the stability transition, our models predict that X-ray bursts have extended tails where freshly accreted fuel prolongs nuclear burning.

Effective values for the key helium burning reaction rates, triple-alpha and 12C(alpha,gamma)16O, are obtained by adjusting their strengths so as to obtain the best match with the solar abundance pattern of isotopes uniquely or predominately made in core collapse supernovae. These effective rates are then used to determine the production of the neutrino isotopes. The use of effective rates considerably reduces the uncertainties in the production factors arising from uncertainties in the helium burning rates, and improves our ability to use the production of 11B to constrain the neutrino emission from supernovae.

We introduce numerical algorithms for initializing multidimensional simulations of stellar explosions with 1D stellar evolution models. The initial mapping from 1D profiles onto multidimensional grids can generate severe numerical artifacts, one of the most severe of which is the violation of conservation laws for physical quantities. We introduce a numerical scheme for mapping 1D spherically-symmetric data onto multidimensional meshes so that these physical quantities are conserved. We verify our scheme by porting a realistic 1D Lagrangian stellar profile to the new multidimensional Eulerian hydro code CASTRO. Our results show that all important features in the profiles are reproduced on the new grid and that conservation laws are enforced at all resolutions after mapping. We also introduce a numerical scheme for initializing multidimensional supernova simulations with realistic perturbations predicted by 1D stellar evolution models. Instead of seeding 3D stellar profiles with random perturbations, we imprint them with velocity perturbations that reproduce the Kolmogorov energy spectrum expected for highly turbulent convective regions in stars. Our models return Kolmogorov energy spectra and vortex structures like those in turbulent flows before the modes become nonlinear. Finally, we describe approaches to determining the resolution for simulations required to capture fluid instabilities and nuclear burning. Our algorithms are applicable to multidimensional simulations besides stellar explosions that range from astrophysics to cosmology.

The element abundance ratios of four low-mass stars with extremely low metallicities indicate that the gas out of which the stars formed was enriched in each case by at most a few, and potentially only one low-energy, supernova. Such supernovae yield large quantities of light elements such as carbon but very little iron. The dominance of low-energy supernovae is surprising, because it has been expected that the first stars were extremely massive, and that they disintegrated in pair-instability explosions that would rapidly enrich galaxies in iron. What has remained unclear is the yield of iron from the first supernovae, because hitherto no star is unambiguously interpreted as encapsulating the yield of a single supernova. Here we report the optical spectrum of SMSS J031300.36- 670839.3, which shows no evidence of iron (with an upper limit of 10^-7.1 times solar abundance). Based on a comparison of its abundance pattern with those of models, we conclude that the star was seeded with material from a single supernova with an original mass of ~60 Mo (and that the supernova left behind a black hole). Taken together with the previously mentioned low-metallicity stars, we conclude that low-energy supernovae were common in the early Universe, and that such supernovae yield light element enrichment with insignificant iron. Reduced stellar feedback both chemically and mechanically from low-energy supernovae would have enabled first-generation stars to form over an extended period. We speculate that such stars may perhaps have had an important role in the epoch of cosmic reionization and the chemical evolution of early galaxies.

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.

The existence of supermassive black holes as early as z ~ 7 is one of the great unsolved problems in cosmological structure formation. One leading theory argues that they are born during catastrophic baryon collapse in z ~ 15 protogalaxies in strong Lyman-Werner UV backgrounds. Atomic line cooling in such galaxies fragments baryons into massive clumps that are thought to directly collapse to 10^4 - 10^5 solar-mass black holes. We have now discovered that some of these fragments can instead become supermassive stars that eventually explode as pair-instability supernovae with energies of ~ 10^55 erg, the most energetic explosions in the universe. We have calculated light curves and spectra for supermassive Pop III PI SNe with the Los Alamos RAGE and SPECTRUM codes. We find that they will be visible in NIR all-sky surveys by Euclid out to z ~ 10 - 15 and by WFIRST and WISH out to z ~ 15 - 20, perhaps revealing the birthplaces of the first quasars.

One of the leading contenders for the origin of supermassive black holes at $z \gtrsim$ 7 is catastrophic baryon collapse in atomically-cooled halos at $z \sim$ 15. In this scenario, a few protogalaxies form in the presence of strong Lyman-Werner UV backgrounds that quench H$_2$ formation in their constituent halos, preventing them from forming stars or blowing heavy elements into the intergalactic medium prior to formation. At masses of 10$^ 8$ \Ms\ and virial temperatures of 10$^4$ K, gas in these halos rapidly cools by H lines, in some cases forming 10$^4$ - 10$^6$ \Ms\ Pop III stars and, a short time later, the seeds of supermassive black holes. Instead of collapsing directly to black holes some of these stars died in the most energetic thermonuclear explosions in the universe. We have modeled the explosions of such stars in the dense cores of line-cooled protogalaxies in the presence of cosmological flows. In stark contrast to the explosions in diffuse regions in previous simulations, these SNe briefly engulf the protogalaxy but then collapse back into its dark matter potential. Fallback drives turbulence that efficiently distributes metals throughout the interior of the halo and fuels the rapid growth of nascent black holes at its center. The accompanying starburst and x-ray emission from these line-cooled galaxies easily distinguish them from more slowly evolving neighbors and might reveal the birthplaces of supermassive black holes on the sky.

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.

Supermassive primordial stars are expected to form in a small fraction of massive protogalaxies in the early universe, and are generally conceived of as the progenitors of the seeds of supermassive black holes (BHs). Supermassive stars with masses of ~55,000 M_Sun, however, have been found to explode and completely disrupt in a supernova (SN) with an energy of up to ~10^55 erg instead of collapsing to a BH. Such events, ~10,000 times more energetic than typical SNe today, would be among the biggest explosions in the history of the universe. Here we present a simulation of such a SN in two stages. Using the RAGE radiation hydrodynamics code we first evolve the explosion from an early stage through the breakout of the shock from the surface of the star until the blast wave has propagated out to several parsecs from the explosion site, which lies deep within an atomic cooling dark matter (DM) halo at z ~ 15. Then, using the GADGET cosmological hydrodynamics code we evolve the explosion out to several kiloparsecs from the explosion site, far into the low-density intergalactic medium. The host DM halo, with a total mass of 4 x 10^7 M_Sun, much more massive than typical primordial star-forming halos, is completely evacuated of high density gas after < 10 Myr, although dense metal-enriched gas recollapses into the halo, where it will likely form second-generation stars with metallicities of ~ 0.05 Z_Sun after > 70 Myr. The chemical signature of supermassive star explosions may be found in such long-lived second-generation stars today.

All stellar evolution models for nucleosynthesis require an initial isotopic abundance set to use as a starting point. Generally, our knowledge of isotopic abundances of stars is fairly incomplete except for the Solar System. We present a first model for a complete average isotopic decomposition as a function of metallicity. Our model is based on the underlying nuclear astrophysics processes, and is fitted to observational data, rather than traditional forward galactic chemical evolution modeling which integrates stellar yields beginning from big bang nucleosynthesis. We first decompose the isotopic solar abundance pattern into contributions from astrophysical sources. Each contribution is then assumed to scale as a function of metallicity. The resulting total isotopic abundances are summed into elemental abundances and fitted to available halo and disk stellar data to constrain the model's free parameter values. This procedure allows us to use available elemental observational data to reconstruct and constrain both the much needed complete isotopic evolution that is not accessible to current observations, and the underlying astrophysical processes. Our model finds a best fit for Type Ia contributing \simeq 0.7 to the solar Fe abundance, and and Type Ia onset occurring at [Fe/H]\simeq -1.2.

The first primitive galaxies formed from accretion and mergers by z ~ 15, and were primarily responsible for cosmological reionization and the chemical enrichment of the early cosmos. But a few of these galaxies may have formed in the presence of strong Lyman-Werner UV fluxes that sterilized them of H_2, preventing them from forming stars or expelling heavy elements into the IGM prior to assembly. At masses of 10^8 Ms and virial temperatures of 10^4 K, these halos began to rapidly cool by atomic lines, perhaps forming 10^4 - 10^6 Ms Pop III stars and, later, the seeds of supermassive black holes. We have modeled the explosion of a supermassive Pop III star in the dense core of a line-cooled protogalaxy with the ZEUS-MP code. We find that the supernova (SN) expands to a radius of ~ 1 kpc, briefly engulfing the entire galaxy, but then collapses back into the potential well of the dark matter. Fallback fully mixes the interior of the protogalaxy with metals, igniting a violent starburst and fueling the rapid growth of a massive black hole at its center. The starburst would populate the protogalaxy with stars in greater numbers and at higher metallicities than in more slowly-evolving, nearby halos. The SN remnant becomes a strong synchrotron source that can be observed with eVLA and eMERLIN and has a unique signature that easily distinguishes it from less energetic SN remnants. Such explosions, and their attendant starbursts, may well have marked the birthplaces of supermassive black holes on the sky.

We present two new primary mechanisms for the synthesis of the rare nucleus 9Be, both triggered by nu-induced production of 3H followed by 4He(3H,gamma)7Li in the He shells of core-collapse supernovae. For progenitors of ~8M_sun, 7Li(3H,n_0)9Be occurs during the rapid expansion of the shocked He shell. Alternatively, for ultra-metal-poor progenitors of ~11-15M_sun, 7Li(n,gamma)8Li(n,gamma)9Li(e^-anti-nu_e)9Be occurs with neutrons produced by 4He(anti-nu_e,e^+n)3H, assuming a hard effective anti-nu_e spectrum from oscillations (which also leads to heavy element production through rapid neutron capture) and a weak explosion (so the 9Be survives shock passage). We discuss the associated production of 7Li and 11B, noting patterns in LiBeB production that might distinguish the new mechanisms from others.

About a year after core collapse supernova, dust starts to condense in the ejecta. In meteorites, a fraction of C-rich presolar grains (e.g., silicon carbide (SiC) grains of Type-X and low density graphites) are identified as relics of these events, according to the anomalous isotopic abundances. Several features of these abundances remain unexplained and challenge the understanding of core-collapse supernovae explosions and nucleosynthesis. We show, for the first time, that most of the measured C-rich grain abundances can be accounted for in the C-rich material from explosive He burning in core-collapse supernovae with high shock velocities and consequent high temperatures. The inefficiency of the $^{12}$C($\alpha$,$\gamma$)$^{16}$O reaction relative to the rest of the $\alpha$-capture chain at $T > 3.5\times10^8 \mathrm{K}$ causes the deepest He-shell material to be carbon rich and silicon rich, and depleted in oxygen. The isotopic ratio predictions in part of this material, defined here as the C/Si zone, are in agreement with the grain data. The high-temperature explosive conditions that our models reach at the bottom of the He shell, can also be representative of the nucleosynthesis in hypernovae or in the high-temperature tail of a distribution of conditions in asymmetric supernovae. Finally, our predictions are consistent with the observation of large $^{44}$Ca/$^{40}$Ca observed in the grains. This is due to the production of $^{44}$Ti together with $^{40}$Ca in the C/Si zone, and/or to the strong depletion of $^{40}$Ca by neutron captures.

Understanding the properties of Pop III stars is prerequisite to elucidating the nature of primeval galaxies, the chemical enrichment and reionization of the early IGM, and the origin of supermassive black holes. While the primordial IMF remains unknown, recent evidence from numerical simulations and stellar archaeology suggests that some Pop III stars may have had lower masses than previously thought, 15 - 50 \Ms in addition to 50 - 500 \Ms. The detection of Pop III supernovae by JWST, WFIRST or the TMT could directly probe the primordial IMF for the first time. We present numerical simulations of 15 - 40 \Ms Pop III core-collapse SNe done with the Los Alamos radiation hydrodynamics code RAGE. We find that they will be visible in the earliest galaxies out to z ~ 10 - 15, tracing their star formation rates and in some cases revealing their positions on the sky. Since the central engines of Pop III and solar-metallicity core-collapse SNe are quite similar, future detection of any Type II supernovae by next-generation NIR instruments will in general be limited to this epoch.

Recent studies have claimed the existence of very massive stars (VMS) up to 300 solar masses in the local Universe. As this finding may represent a paradigm shift for the canonical stellar upper-mass limit of 150 Msun, it is timely to discuss the status of the data, as well as the far-reaching implications of such objects. We held a Joint Discussion at the General Assembly in Beijing to discuss (i) the determination of the current masses of the most massive stars, (ii) the formation of VMS, (iii) their mass loss, and (iv) their evolution and final fate. The prime aim was to reach broad consensus between observers and theorists on how to identify and quantify the dominant physical processes.

We study the sensitivity of presupernova evolution and supernova nucleosynthesis yields of massive stars to variations of the helium-burning reaction rates within the range of their uncertainties. We use the current solar abundances from Lodders (2009) for the initial stellar composition. We compute a grid of 12 initial stellar masses and 176 models per stellar mass to explore the effects of independently varying the 12^C(a,g)16^O and 3a reaction rates, denoted R_a12 and R_3a, respectively. The production factors of both the intermediate-mass elements (A=16-40) and the s-only isotopes along the weak s-process path (70Ge, 76Se, 80Kr, 82Kr, 86Sr, and 87Sr) were found to be in reasonable agreement with predictions for variations of R_3a and R_a12 of +/-25%; the s-only isotopes, however, tend to favor higher values of R_3a than the intermediate-mass isotopes. The experimental uncertainty (one standard deviation) in R_3a(R_a12) is approximately +/-10%(+/-25%). The results show that a more accurate measurement of one of these rates would decrease the uncertainty in the other as inferred from the present calculations. We also observe sharp changes in production factors and standard deviations for small changes in the reaction rates, due to differences in the convection structure of the star. The compactness parameter was used to assess which models would likely explode as successful supernovae, and hence contribute explosive nucleosynthesis yields. We also provide the approximate remnant masses for each model and the carbon mass fractions at the end of core-helium burning as a key parameter for later evolution stage

Understanding the progenitors of core collapse supernovae and their population statistics is a key ingredient for many current studies in astronomy but as yet this remains elusive. Using the MESA stellar evolution code we study the dependence of the lower mass limit for making core collapse supernovae (SNe) as function of initial stellar metallicity. We find that this mass limit is smallest at approximately [Z] = -2 with a value of ~ 8.3 Msun. At [Z] = 0 the limit is ~ 9.5 Msun and continues to rise with higher metallicity. As a consequence, for a fixed initial mass function the supernova rate may be 20% to 25% higher at [Z] = -2. This affects the association of observed SN rates as a probe for the cosmological star formation rate, rate predictions for supernova surveys, and population synthesis studies.

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.

The neutrino-driven explosion mechanism for core-collapse supernovae in its modern flavor relies on the additional support of hydrodynamical instabilities in achieving shock revival. Two possible candidates, convection and the so-called standing accretion shock instability (SASI), have been proposed for this role. In this paper, we discuss new successful simulations of supernova explosions that shed light on the relative importance of these two instabilities. While convection has so far been observed to grow first in self-consistent hydrodynamical models with multi-group neutrino transport, we here present the first such simulation in which the SASI grows faster while the development of convection is initially inhibited. We illustrate the features of this SASI-dominated regime using an explosion model of a 27 solar mass progenitor, which is contrasted with a convectively-dominated model of an 8.1 solar mass progenitor with subsolar metallicity, whose early post-bounce behavior is more in line with previous 11.2 and 15 solar mass explosion models. We analyze the conditions discriminating between the two different regimes, showing that a high mass-accretion rate and a short advection time-scale are conducive for strong SASI activity. We also briefly discuss some important factors for capturing the SASI-driven regime, such as general relativity, the progenitor structure, a nuclear equation of state leading to a compact proto-neutron star, and the neutrino treatment. Finally, we evaluate possible implications of our findings for 2D and 3D supernova simulations. Our results show that a better understanding of the SASI and convection in the non-linear regime is required.

The LOFT mission concept is one of four candidates selected by ESA for the M3 launch opportunity as Medium Size missions of the Cosmic Vision programme. The launch window is currently planned for between 2022 and 2024. LOFT is designed to exploit the diagnostics of rapid X-ray flux and spectral variability that directly probe the motion of matter down to distances very close to black holes and neutron stars, as well as the physical state of ultra-dense matter. These primary science goals will be addressed by a payload composed of a Large Area Detector (LAD) and a Wide Field Monitor (WFM). The LAD is a collimated (<1 degree field of view) experiment operating in the energy range 2-50 keV, with a 10 m^2 peak effective area and an energy resolution of 260 eV at 6 keV. The WFM will operate in the same energy range as the LAD, enabling simultaneous monitoring of a few-steradian wide field of view, with an angular resolution of <5 arcmin. The LAD and WFM experiments will allow us to investigate variability from submillisecond QPO's to year-long transient outbursts. In this paper we report the current status of the project.

Superbursts are rare day-long Type I X-ray bursts due to carbon flashes on accreting neutron stars in low-mass X-ray binaries. They heat the neutron star envelope such that the burning of accreted hydrogen and helium becomes stable, and the common shorter X-ray bursts are quenched. Short bursts reappear only after the envelope cools down. We study multi-zone one-dimensional models of the neutron star envelope, in which we follow carbon burning during the superburst, and we include hydrogen and helium burning in the atmosphere above. We investigate both the case of a solar composition and a helium-rich atmosphere. This allows us to study for the first time a wide variety of thermonuclear burning behavior as well as the transitions between the different regimes in a self-consistent manner. For solar composition, burst quenching ends much sooner than previously expected. This is because of the complex interplay between the 3-alpha, hot CNO, and CNO breakout reactions. Stable burning of hydrogen and helium transitions via marginally stable burning (mHz QPOs) to less energetic bursts with short recurrence times. We find a short-lived bursting mode where weaker and stronger bursts alternate. Eventually the bursting behavior changes back to that of the pre-superburst bursts. Because of the scarcity of observations, this transition has not been directly detected after a superburst. Using the MINBAR burst catalog we identify the shortest upper limit on the quenching time for 4U 1636-536, and derive further constraints on the time scale on which bursts return.

Mapping one-dimensional stellar profiles onto multidimensional grids as initial conditions for hydrodynamics calculations can lead to numerical artifacts, one of the most severe of which is the violation of conservation laws for physical quantities such as energy and mass. Here we introduce a numerical scheme for mapping one-dimensional spherically-symmetric data onto multidimensional meshes so that these physical quantities are conserved. We validate our scheme by porting a realistic 1D Lagrangian stellar profile to the new multidimensional Eulerian hydro code CASTRO. Our results show that all important features in the profiles are reproduced on the new grid and that conservation laws are enforced at all resolutions after mapping.

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.

Superbursts are rare and energetic thermonuclear carbon flashes observed to occur on accreting neutron stars. We create the first multi-zone models of series of superbursts using a stellar evolution code. We self-consistently build up the fuel layer at different rates, spanning the entire range of observed mass accretion rates for superbursters. For all models light curves are presented. They generally exhibit a shock breakout, a precursor burst due to shock heating, and a two-component power-law decay. Shock heating alone is sufficient for a bright precursor, that follows the shock breakout on a short dynamical time scale due to the fall-back of expanded layers. Models at the highest accretion rates, however, lack a shock breakout, precursor, and the first power law decay component. The ashes of the superburst that form the outer crust are predominantly composed of iron, but a superburst leaves a silicon-rich layer behind in which the next one ignites. Comparing the model light curves to an observed superburst from 4U 1636-53, we find for our accretion composition the best agreement with a model at three times the observed accretion rate. We study the dependence on crustal heating of observables such as the recurrence time and the decay time scale. It remains difficult, however, to constrain crustal heating, if there is no good match with the observed accretion rate, as we see for 4U 1636-53.

Theoretical models suggest that the first stars in the universe could have been very massive, with typical masses $\gtrsim$ 100 \Msun. Many of them might have died as energetic thermonuclear explosions known as pair-instability supernovae (PSNe). We present multidimensional numerical simulations of PSNe with the new radiation-hydrodynamics code CASTRO. Our models capture all explosive burning and follow the explosion until the shock breaks out from the stellar surface. 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 sec after the explosion begins. Later, when the shock reaches the hydrogen envelope a strong reverse shock forms that rapidly develops additional Rayleigh-Taylor instabilities. In red supergiant progenitors, the amplitudes of these instabilities are sufficient to mix the supernova's ejecta and alter its observational signature. Our results provide useful predictions for the detection of PSNe by forthcoming telescopes.

We have studied the sensitivity to variations in the triple alpha and 12C(alpha, gamma)16O reaction rates, of the yield of the neutrino process isotopes 7Li,11B, 19F, 138La, and 180Ta in core collapse supernovae. Compared to solar abundances, less than 15% of 7Li, about 25-80% of 19F, and about half of 138La is produced in these stars. Over a range of two sigma for each helium-burning rate, 11B is overproduced and the yield varies by an amount larger than the variation caused by the effects of neutrino oscillations. The total 11B yield, however, may eventually provide constraints on supernova neutrino spectra.

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.

We have studied the sensitivity of supernova production of the gamma emitting nuclei 26Al, 44Ti and 60Fe to variations of the rates of the triple alpha and 12C(alpha, gamma) reactions. Over a range of twice their experimental uncertainties we find variations in the production of 60Fe by more than a factor of five. Smaller variations, about a factor of two to three, were observed for 26Al and 44Ti. The yields of these isotopes change significantly when the abundances of Lodders (2003) are used instead of those of Anders and Grevesse (1989). These sensitivities will limit conclusions based on a comparison of observed gamma ray intensities and stellar models until the helium burning rates are better known. Prospects for improving the helium burning rates are discussed and a new version of the Boyes rate for 12C(alpha, gamma} is presented.

We present preliminary results from two dimensional numerical studies of pair instability supernova (PSN). We study nuclear burning, hydrodynamic instabilities and explosion of very massive stars. Use a new radiation-hydrodynamics code, CASTRO.

We present preliminary results from multidimensional numerical studies of pair instability supernova (PSN), studying the fluid instabilities that occur in multiple spatial dimensions. We use the new radiation-hydrodynamics code, CASTRO, and introduce a new mapping procedure that defines the initial conditions for the multidimensional runs in such a way that conservation of physical quantities is guaranteed at any level of resolution.

Type I X-ray bursts from low-mass X-ray binaries result from a thermonuclear runaway in the material accreted onto the neutron star. Although typical recurrence times are a few hours, consistent with theoretical ignition model predictions, there are also observations of bursts occurring as promptly as ten minutes or less after the previous event. We present a comprehensive assessment of this phenomenon using a catalog of 3387 bursts observed with the BeppoSAX/WFCs and RXTE/PCA X-ray instruments. This catalog contains 136 bursts with recurrence times of less than one hour, that come in multiples of up to four events, from 15 sources. Short recurrence times are not observed from so-called ultra-compact binaries, indicating that hydrogen burning processes play a crucial role. As far as the neutron star spin frequency is known, these sources all spin fast at over 500 Hz; the rotationally induced mixing may explain burst recurrence times of the order of 10 min. Short recurrence time bursts generally occur at all mass accretion rates where normal bursts are observed, but for individual sources the short recurrence times may be restricted to a smaller interval of accretion rate. The fraction of such bursts is roughly 30%. We also report the shortest known recurrence time of 3.8 minutes.

We have studied the sensitivity to variations in the the triple alpha and 12C(a,g)16O reaction rates of the production of 26Al, 44Ti, and 60Fe in core-collapse supernovae. We used the KEPLER code to model the evolution of 15, 20, and 25 solar mass stars to the onset of core collapse and simulated the ensuing supernova explosion using a piston model for the explosion. Calculations were performed for the Anders and Grevesse (1989) and Lodders (2003) abundances. Over a range of twice the experimental uncertainty, sigma, for each helium-burning rate, the production of 26Al, 60Fe, and their ratio vary by factors of five or more. For some species, similar variations were observed for much smaller rate changes, 0.5 sigma or less. The production of 44Ti was less sensitive to changes in the helium-burning rates. Production of all three isotopes depended on the solar abundance set used for the initial stellar composition.

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.

We have studied the sensitivity of s-process nucleosynthesis in massive stars to +/- 2 sigma variations in the rates of the triple alpha and 12C(a,g)16O reactions. We simulated the evolution of massive stars from H-burning through Fe-core collapse, followed by a supernova explosion. We found that: the production factors of s-process nuclides between 58Fe and 96Zr change strongly with changes in the He burning reaction rates; using the Lodders (2003) solar abundances rather than those of Anders and Grevesse (1989) reduces s-process nucleosynthesis; later burning phases beyond core He burning and the shell C burning have a significant effect on post-explosive production factors. We also discuss the implications of the uncertainties in the helium burning rates for evidence of a new primary neutron capture process (LEPP) in massive stars.

We report the discovery of one-hour long tails on the few-minutes long X-ray bursts from the `clocked burster' GS 1826-24. We propose that the tails are due to enduring thermal radiation from the neutron star envelope. The enduring emission can be explained by cooling of deeper NS layers which were heated up through inward conduction of heat produced in the thermonuclear shell flash responsible for the burst. Similar, though somewhat shorter, tails are seen in bursts from EXO 0748-676 and 4U 1728-34. Only a small amount of cooling is detected in all these tails. This is either due to compton up scattering of the tail photons or, more likely, to a NS that is already fairly hot due to other, stable, nuclear processes.

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 explore the sensitivity of massive stars to neutrino magnetic moments. We find that the additional cooling due to the neutrino magnetic moments bring about qualitative changes to the structure and evolution of stars in the mass window 7 Msun < M < 18 Msun, rather than simply changing the time scales for the burning. We describe some of the consequences of this modified evolution: the shifts in the threshold masses for creating core-collapse supernovae and oxygen-neon-magnesium white dwarfs and the appearance of a new type of supernova in which a partial carbon-oxygen core explodes within a massive star. The resulting sensitivity to the magnetic moment is at the level of (2-4) * 10^{-11} \mu_B.

In the course of the life of a massive star, wind-wind interaction can give rise to the formation of circumstellar nebulae which are both predicted and observed in the nature. We present generic model calculations to predict the properties of such nebulae for blue supergiants. From stellar evolution calculations including rotation, we obtain the time dependence of the stellar wind properties and of the stellar radiation field. These are used as input for hydro-calculations of the circumstellar medium throughout the star's life. Here, we present the results for a rapidly rotating 12 solar masses single star. This star undergoes a blue loop during its post main sequence evolution, at the onset of which its contraction spins it up close to critical rotation. Due to the consequent anisotropic mass loss, the blue supergiant wind sweeps up the preceding slow wind into an hour glass structure. Its collision with the previously formed spherical red supergiant wind shell forms a short-lived luminous nebula consisting of two polar caps and a central inner ring. With time, the polar caps evolve into mid-latitude rings which gradually move toward the equatorial plane while the central ring is fading. These structures are reminiscent to the observed nebulae around the blue supergiant Sher 25 and the progenitor of SN 1987A. The simple model of an hour glass colliding with a spherical shell retrieves most of the intriguing nebula geometries discovered around blue supergiants, and suggests them to form an evolutionary sequence. Our results indicate that binarity is not required to obtain them.

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.

Recent analysis of five gamma-ray burst (GRB) afterglow spectra reveal the absence of molecular hydrogen absorption lines, a surprising result in light of their large neutral hydrogen column densities and the detection of H$_2$ in similar, more local star-forming regions like 30 Doradus in the Large Magellanic Cloud (LMC). Observational evidence further indicates that the bulk of the neutral hydrogen column in these sight lines lies 100 pc beyond the progenitor and that H$_2$ was absent prior to the burst, suggesting that direct flux from the star, FUV background fields, or both suppressed its formation. We present one-dimensional radiation hydrodynamical models of GRB host galaxy environments, including self-consistent radiative transfer of both ionizing and Lyman-Werner photons, nine-species primordial chemistry with dust formation of H$_2$, and dust extinction of UV photons. We find that a single GRB progenitor is sufficient to ionize neutral hydrogen to distances of 50 - 100 pc but that a galactic Lyman-Werner background is required to dissociate the molecular hydrogen in the ambient ISM. Intensities of 0.1 - 100 times the Galactic mean are necessary to destroy H$_2$ in the cloud, depending on its density and metallicity. The minimum radii at which neutral hydrogen will be found in afterglow spectra is insensitive to the mass of the progenitor or the initial mass function (IMF) of its cluster, if present.

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)

The prediction of the spins of the compact remnants is a fundamental goal of the theory of stellar evolution. Here, we confront the predictions for white dwarf spins from evolutionary models including rotation with observational constraints. We perform stellar evolution calculations for stars in the mass range 1... 3$\mso$, including the physics of rotation, from the zero age main sequence into the TP-AGB stage. We calculate two sets of model sequences, with and without inclusion of magnetic fields. From the final computed models of each sequence, we deduce the angular momenta and rotational velocities of the emerging white dwarfs. While models including magnetic torques predict white dwarf rotational velocities between 2 and 10 km s$^{-1}$, those from the non-magnetic sequences are found to be one to two orders of magnitude larger, well above empirical upper limits. We find the situation analogous to that in the neutron star progenitor mass range, and conclude that magnetic torques may be required in order to understand the slow rotation of compact stellar remnants in general.

The understanding of the formation, life, and death of Population III stars, as well as the impact that these objects had on later generations of structure formation, is one of the foremost issues in modern cosmological research and has been an active area of research during the past several years. We summarize the results presented at "First Stars III," a conference sponsored by Los Alamos National Laboratory, the Kavli Institute for Particle Astrophysics and Cosmology, and the Joint Institute for Nuclear Astrophysics. This conference, the third in a series, took place in July 2007 at the La Fonda Hotel in Santa Fe, New Mexico, U.S.A.

We explore the dependence of pre-supernova evolution and supernova nucleosynthesis yields on the uncertainties in helium burning reaction rates. Using the revised solar abundances of Lodders (2003) for the initial stellar composition, instead of those of Anders & Grevesse (1989), changes the supernova yields and limits the constraints that those yields place on the 12C(a,g)16O reaction rate. The production factors of medium-weight elements (A = 16-40) were found to be in reasonable agreement with observed solar ratios within the current experimental uncertainties in the triple alpha reaction rate. Simultaneous variations by the same amount in both reaction rates or in either of them separately, however, can induce significant changes in the central 12C abundance at core carbon ignition and in the mass of the supernova remnant. It therefore remains important to have experimental determinations of the helium burning rates so that their ratio and absolute values are known with an accuracy of 10% or better.

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.

We produce models of early WN, WC, and WO stars as a function of metallicity Z using an analytic CAK-type approach. At log(Z / Z_sun) >= -2 both WN and WC stars have the approximate dependence Mdot \propto Z^{0.5}. For a WN wind the mass-loss rate drops rapidly below log(Z / Z_sun) = -2, and no wind solution can be found for log(Z / Z_sun) < -3. For WC and WO winds the mass-loss rate plummets in the range -3 <= log(Z / Z_sun) <= -2 and tends to flatten due to the self-enrichment of C and O to around 10^{-8} M_sun yr^{-1} for log(Z / Z_sun) <= -4. No significant difference in Mdot was found for WC versus WO stars at low metallicity.

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.

We study the late evolution of solar metallicity stars in the transition region between white dwarf formation and core collapse. This includes the super-asymptotic giant branch (super-AGB, SAGB) stars, which have massive enough cores to ignite carbon burning and form an oxygen-neon (ONe) core. The most massive SAGB stars have cores that may grow to the Chandrasekhar mass because of continued shell-burning. Their cores collapse, triggering a so called electron capture supernovae (ECSN). From stellar evolution models we find that the initial mass range for SAGB evolution is 7.5 ... 9.25\msun. We perform calculations with three different stellar evolution codes to investigate the sensitivity of this mass range to some of the uncertainties in current stellar models. The mass range significantly depends on the treatment of semiconvective mixing and convective overshooting. To consider the effect of a large number of thermal pulses, as expected in SAGB stars, we construct synthetic SAGB models that include a semi-analytical treatment of dredge-up, hot-bottom burning, and thermal pulse properties. This synthetic model enables us to compute the evolution of the main properties of SAGB stars from the onset of thermal pulses until the core reaches the Chandrasekhar mass or is uncovered by the stellar wind. Thereby, we determine the stellar initial mass ranges that produce ONe-white dwarfs and electron-capture supernovae. The latter is found to be 9.0 ... 9.25\msun for our fiducial model, implying that electron-capture supernovae would constitute about 4% of all supernovae in the local universe. Our synthetic approach allows us to explore the uncertainty of this number imposed by uncertainties in the third dredge-up efficiency and ABG mass loss rate. We find for ECSNe a upper limit of ~20% of all supernovae (abridged).

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.

We study the evolution of first star (Population III) binaries. Under specific conditions, these stars may produce high redshift gamma-ray bursts (GRBs). We demonstrate that the occurrence rate of GRBs does not depend sensitively on evolutionary parameters in the population synthesis models. We show that the first binaries may form a very small group (< 1%) of fast rotating stars through binary tidal interactions that make GRBs. This finding is contrary to the Bromm & Loeb assumption that all stars in close Population III binaries will be spun up by tides and produce a GRB. We find that there is simply not enough fast rotating stars in Population III binaries to expect detection with SWIFT. Predicted detection rates, even with very optimistic assumptions on binary fraction, evolutionary parameters and GRB detection, are very small: 0.1-0.01 per year.

First steps have been taken in a more comprehensive study of the dependence of observables in Type I X-ray bursts on uncertain (p,gamma) reaction rates along the rp-process path. We use the multizone hydrodynamics code KEPLER which implicitly couples a full nuclear reaction network of more than 1000 isotopes, as needed, to follow structure and evolution of the X-ray burst layer and its ashes. This allows us to incorporate the full rp-process network, including all relevant nuclear reactions, and individually study changes in the X-ray burst light curves when modifying selected key nuclear reaction rates. In this work we considered all possible proton captures to nuclei with 10 < Z < 28 and N <= Z. When varying individual reaction rates within a symmetric full width uncertainty of a factor of 10000, early results for some rates show changes in the burst light curve as large as 10 percent of peak luminosity. This is very large compared to the current sensitivity of X-ray observations. More precise reaction rates are therefore needed to test current X-ray burst models, particularly of the burst rise, with observational data and to constrain astrophysical parameters.

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.

The 44Ti(t1/2 = 59 y) nuclide, an important signature of supernova nucleosynthesis, has recently been observed as live radioactivity by gamma-ray astronomy from the Cas A remnant. We investigate in the laboratory the major 44Ti production reaction, 40Ca(alpha,gamma)44Ti (E_cm = 0.6-1.2 MeV/u), by direct off- line counting of 44Ti nuclei. The yield, significantly higher than inferred from previous experiments, is analyzed in terms of a statistical model using microscopic nuclear inputs. The associated stellar rate has important astrophysical consequences, increasing the calculated supernova 44Ti yield by a factor ~2 over previous estimates and bringing it closer to Cas A observations.

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.

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.

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.

The collapsar model for gamma-ray bursts requires three essential ingredients: a massive core, removal of the hydrogen envelope, and enough angular momentum in the core. We study current massive star evolution models of solar metallicity to determine which massive star physics is capable of producing these ingredients. In particular, we investigate the role of hydrodynamic and magnetic internal angular momentum transport and binary mass and angular momentum transfer. We follow the evolution of rotating single stars and of binary systems that include rotational processes for both stars. Neglecting magnetic fields, we show that the cores of massive single stars can maintain a high specific angular momentum when evolved with the assumption that mean molecular weight gradients suppress rotational mixing processes. In binary systems that undergo mass transfer during core hydrogen burning the mass receiving star accretes large amounts of high angular momentum material, leading to a spin-up of the core. We find, however, that this merely compensates for the tidal angular momentum loss due to spin-orbit coupling, which leads to synchronous rotation before the mass transfer event. Therefore the resulting cores do not rotate faster than in single stars. We also present models that include magnetic fields generated by differential rotation and we consider the internal angular momentum transport by magnetic torques. We investigate the capability of magnetic torques to efficiently pump angular momentum into the cores of accreting stars. Despite our finding that this mechanism works, the magnetic coupling of core and envelope after the accreting star ends core hydrogen burning leads to slower rotation than in the non-magnetic case.

The collapsar model, the now leading model for the engine behind gamma-ray bursts and hypernovae, requires that a star collapses to form a black hole surrounded by an accretion disk of high-angular momentum material. The current best theoretical stellar models, however, do not retain enough angular momentum in the core of the star to make a centrifugally supported disk. In this paper, we present the first calculations of the helium-star/helium-star merger progenitors for the collapsar model. These progenitors invoke the merger of two helium cores during the common envelope inspiral phase of a binary system. We find that, in some cases, the merger can produce cores that are rotating 3-10 times faster than single stars. He-star/He-star gamma-ray burst progenitors have a very different redshift distribution than their single-star gamma-ray burst progenitors and we discuss how gamma-ray burst observations can constrain these progenitors.

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

The mechanism responsible for the natal kicks of neutron stars continues to be a challenging problem. Indeed, many mechanisms have been suggested, and one hydrodynamic mechanism may require large initial asymmetries in the cores of supernova progenitor stars. Goldreich, Lai, & Sahrling (1997) suggested that unstable g-modes trapped in the iron (Fe) core by the convective burning layers and excited by the $\epsilon$-mechanism may provide the requisite asymmetries. We perform a modal analysis of the last minutes before collapse of published core structures and derive eigenfrequencies and eigenfunctions, including the nonadiabatic effects of growth by nuclear burning and decay by both neutrino and acoustic losses. In general, we find two types of g-modes: inner-core g-modes, which are stabilized by neutrino losses and outer-core g-modes which are trapped near the burning shells and can be unstable. Without exception, we find at least one unstable g-mode for each progenitor in the entire mass range we consider, 11 M$_{\sun}$ to 40 M$_{\sun}$. More importantly, we find that the timescales for growth and decay are an order of magnitude or more longer than the time until the commencement of core collapse. We conclude that the $\epsilon$-mechanism may not have enough time to significantly amplify core g-modes prior to collapse.

We systematically examine how the presence in a binary affects the final core structure of a massive star and its consequences for the subsequent supernova explosion. Interactions with a companion star may change the final rate of rotation, the size of the helium core, the strength of carbon burning and the final iron core mass. Stars with initial masses larger than \sim 11\Ms that experiece core collapse will generally have smaller iron cores at the time of the explosion if they lost their envelopes due to a previous binary interaction. Stars below \sim 11\Ms, on the other hand, can end up with larger helium and metal cores if they have a close companion, since the second dredge-up phase which reduces the helium core mass dramatically in single stars does not occur once the hydrogen envelope is lost. We find that the initially more massive stars in binary systems with masses in the range 8 - 11\Ms are likely to undergo an electron-capture supernova, while single stars in the same mass range would end as ONeMg white dwarfs. We suggest that the core collapse in an electron-capture supernova (and possibly in the case of relatively small iron cores) leads to a prompt explosion rather than a delayed neutrino-driven explosion and that this naturally produces neutron stars with low-velocity kicks. This leads to a dichotomous distribution of neutron star kicks, as inferred previously, where neutron stars in relatively close binaries attain low kick velocities. We illustrate the consequences of such a dichotomous kick scenario using binary population synthesis simulations and discuss its implications. This scenario has also important consequences for the minimum initial mass of a massive star that becomes a neutron star. (Abbreviated.)

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.

We investigate the dependence of the carbon and oxygen production in stars on the 3-alpha rate by varying the energy of the 02+-state of 12C and determine the resulting yields for a selection of low-mass, intermediate-mass, and massive stars. The yields are obtained using modern stellar evolution codes that follow the entire evolution of massive stars, including the supernova explosion, and consider in detail the 3rd dredge-up process during the thermally pulsating asymptotic giant branch of low-mass and intermediate-mass stars. Our results show that the C and O production in massive stars depends strongly on the initial mass, and that it is crucial to follow the entire evolution. A rather strong C production during the He-shell flashes compared to quiescent He burning leads to a lower sensitivity of the C and O production in low-mass and intermediate-mass stars on the 3-alpha-rate than predicted in our previous work. In particular, the C production of intermediate-mass stars seems to have a maximum close to the actual value of the 02+ energy level of 12C.

Many classical nova ejecta are enriched in CNO and Ne. Rosner et al. recently suggested that the enrichment might originate in the resonant interaction between large-scale shear flows in the accreted H/He envelope and gravity waves at the interface between the envelope and the underlying C/O white dwarf. The shear flow amplifies the waves, which eventually form cusps and break. This wave breaking injects a spray of C/O into the superincumbent H/He. Using two-dimensional simulations, we formulate a quantitative expression for the amount of C/O per unit area that can be entrained into the H/He at saturation. The fraction of the envelope that is enriched depends on the horizontal distribution of shear velocity and the density contrast between the C/O white dwarf and the H/He layer but is roughly independent of the vertical shape of the shear profile. Using this parameterization for the mixed mass, we then perform several one-dimensional Lagrangian calculations of an accreting white dwarf envelope and consider two scenarios: that the wave breaking and mixing is driven by the convective flows; and that the mixing occurs prior to the onset of convection. In the absence of enrichment prior to ignition, the base of the convective zone, as calculated from mixing-length theory with the Ledoux instability criterion, does not reach the C/O interface. As a result, there is no additional mixing, and the runaway is slow. In contrast, the formation of a mixed layer during the accretion of H/He, prior to ignition, causes a more violent runaway. The envelope can be enriched by <25% of C/O by mass (consistent with that observed in some ejecta) for shear velocities, over the surface, with Mach numbers <0.4.

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.

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.

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.

We discuss the first available binary evolution models which include up-to-date rotational physics for both components, as well as angular momentum accretion and spin-orbit coupling. These models allow a self-consistent computation of the mass transfer efficiency during Roche-lobe overflow phases, and a determination of the transition from quasi-conservative to non-conservative evolution. Applications to massive binary systems lead to predictions for the spin rates of compact objects in binaries, and for the occurrence of gamma-ray bursts from collapsars in binaries. Rotational effects in accreting white dwarfs are found to stabilise the shell burning and decrease the carbon abundance in progenitor models for Chandrasekhar-mass Type Ia supernovae, and to potentially avoid a detonation of the white dwarf within the sub-Chandrasekhar mass scenario.

We present new mechanisms by which the isotopic compositions of X-type grains of presolar SiC are altered by reverse shocks in Type II supernovae. We address three epochs of reverse shocks: pressure wave from the H envelope near t = 10$^6$s; reverse shock from the presupernova wind near 10$^8-10^9$s; reverse shock from the ISM near 10$^{10}$s. Using 1-D hydrodynamics we show that the first creates a dense shell of Si and C atoms near 10$^6$s in which the SiC surely condenses. The second reverse shock causes precondensed grains to move rapidly forward through decelerated gas of different isotopic composition, during which implantation, sputtering and further condensation occur simultaneously. The third reverse shock causes only further ion implantation and sputtering, which may affect trace element isotopic compositions. Using a 25M$_{\odot}$ supernova model we propose solutions to the following unsolved questions: where does SiC condense?; why does SiC condense in preference to graphite?; why is condensed SiC $^{28}$Si-rich?; why is O richness no obstacle to SiC condensation?; how many atoms of each isotope are impacted by a grain that condenses at time t$_0$ at radial coordinate r$_0$? These many considerations are put forward as a road map for interpreting SiC X grains found in meteorites and their meaning for supernova physics.

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.

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.

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.

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.

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 present models for the complete life and death of a 60 solar mass star evolving in a close binary system, from the main sequence phase to the formation of a compact remnant and fallback of supernova debris. After core hydrogen exhaustion, the star expands, loses most of its envelope by Roche lobe overflow, and becomes a Wolf-Rayet star. We study its post-mass transfer evolution as a function of the Wolf-Rayet wind mass loss rate (which is currently not well constrained and will probably vary with initial metallicity of the star). Varying this mass loss rate by a factor 6 leads to stellar masses at collapse that range from 3.1 to 10.7 solar masses. Although the iron core masses at collapse are generally larger for stars with larger final masses, they do not depend monotonically on the final stellar mass or even the C/O-core mass. We then compute the evolution of all models through collapse and bounce. The results range from strong supernova explosions for the lower final masses to the direct collapse of the star into a black hole for the largest final mass. Correspondingly, the final remnant masses, which were computed by following the supernova evolution and fallback of material for a time scale of about one year, are between 1.2 and 10 solar masses. We discuss the remaining uncertainties of this result and outline the consequences of our results for the understanding of the progenitor evolution of X-ray binaries and gamma-ray burst models.

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.

The discrepancy in the past years of many more black-hole soft X-ray transients (SXTs), of which a dozen have now been identified, had challenged accepted wisdom in black hole evolution. The sheer number of inferred SXTs requires that many, if not most, stars of ZAMS masses $20-35\msun$ end up as black holes (Portegies Zwart et al. 1997; Ergma & van den Heuvel 1998). In this paper we show that this can be understood by challenging the accepted wisdom that the result of helium core burning in a massive star is independent of whether the core is covered by a hydrogen envelope, or "naked" while it burns. The latter case occurs in binaries when the envelope of the more massive star is transferred to the companion by Roche Lobe overflow while in either main sequence or red giant stage. For solar metallicity, whereas the helium cores which burn while naked essentially never go into high-mass black holes, those that burn while clothed do so, beginning at ZAMS mass $\sim 20\msun$, the precise mass depending on the $^{12}C(\alpha,\gamma)^{16}$O rate as we outline. In this way the SXTs can be evolved, provided that the H envelope of the massive star is removed only following the He core burning. We also show that in order to evolve a black hole of mass $\gsim 10\msun$ such as observed in Cyg X-1, even employing extremely massive progenitors of ZAMS mass $\gsim 60\msun$ for the black hole, the core must be covered by hydrogen during a substantial fraction of the core burning. In other words, the progenitor must be a WNL star. We evolve Cyg X-1 in an analogous way to which the SXTs are evolved, the difference being that the companion in Cyg X-1 is more massive than those in the SXTs, so that Cyg X-1 shines continuously.

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.

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 investigate the evolution of the surface properties of models for rotating massive stars, i.e., their luminosities, effective temperatures, surface rotational velocities, and surface abundances of all isotopes, from the zero age main sequence to the supernova stage. Our results are based on the grid of stellar models by Heger, Langer, & Woosley, which covers solar metallicity stars in the initial mass range 8 - 25 solar masses. Results are parameterized by initial mass, initial rotational velocity and major uncertainties in the treatment of the rotational mixing inside massive stars. Rotationally induced mixing processes widen the main sequence and increase the core hydrogen burning lifetime, similar to the effects of convective overshooting. It can also significantly increase the luminosity during and after core hydrogen burning, and strongly affects the evolution of the effective temperature. Our models predict surface rotational velocities for various evolutionary stages, in particular for blue supergiants, red supergiants, and for the immediate presupernova stage. We discuss the changes of the surface abundances due to rotationally induced mixing for main sequence and post main sequence stars. We single out two characteristics by which the effect of rotational mixing can be distinguished from that of massive close binary mass transfer, the only alternative process leading to non-standard chemical surface abundances in massive stars. A comparison with observed abundance anomalies in various types of massive stars supports the concept of rotational mixing in massive stars and indicates that it is responsible for most of the observed abundance anomalies.

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.

We present the results from a series of two-dimensional core-collapse simulations using a rotating progenitor star. We find that the convection in these simulations is less vigorous because a) rotation weakens the core bounce which seeds the neutrino-driven convection and b) the angular momentum profile in the rotating core stabilizes against convection. The limited convection leads to explosions which occur later and are weaker than the explosions produced from the collapse of non-rotating cores. However, because the convection is constrained to the polar regions, when the explosion occurs, it is stronger along the polar axis. This asymmetric explosion can explain the polarization measurements of core-collapse supernovae. These asymmetries also provide a natural mechanism to mix the products of nucleosynthesis out into the helium and hydrogen layers of the star. We also discuss the role the collapse of these rotating stars play on the generation of magnetic fields and neutron star kicks. Given a range of progenitor rotation periods, we predict a range of supernova energies for the same progenitor mass. The critical mass for black hole formation also depends upon the rotation speed of the progenitor.

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.

We present the first evolutionary models of intermediate mass stars up to their thermal pulses which include effects of rotation on the stellar structure as well as rotationally induced mixing of chemical species and angular momentum. We find a significant angular momentum transport from the core to the hydrogen-rich envelope and obtain a white dwarf rotation rate comparable to current observational upper limits of 50 km/s. Large angular momentum gradients at the bottom of the convective envelope and the tip of the pulse driven convective shell are shown to produce marked chemical mixing between the proton-rich and the 12C-rich layers during the so called third dredge-up. This leads to a subsequent production of 13C which is followed by neutron production through 13C(alpha,n) in radiative layers in between thermal pulses. Although uncertainties in the efficiency of rotational mixing processes persist, we conclude that rotation is capable of producing a 13C-rich layer as required for the occurrence of the s-process in TP-AGB stars.

We report on a mechanism which may lead to a spin-up of the surface of a rotating single star leaving the Hayashi line, which is much stronger than the spin-up expected from the mere contraction of the star. By analyzing rigidly rotating, convective stellar envelopes, we qualitatively work out the mechanism through which these envelopes may be spun up or down by mass loss through their lower or upper boundary, respectively. We find that the first case describes the situation in retreating convective envelopes, which tend to retain most of the angular momentum while becoming less massive, thereby increasing the specific angular momentum in the convection zone and thus in the layers close to the stellar surface. We explore the spin-up mechanism quantitatively in a stellar evolution calculation of a rotating 12 M_sun star, which is found to be spun up to critical rotation after leaving the red supergiant branch. We discuss implications of this spin-up for the circumstellar matter around several types of stars, i.e., post-AGB stars, B[e] stars, pre-main sequence stars, and, in particular, the progenitor of Supernova 1987A.

We summarize the present status of the predictions of massive star models for the evolution of their surface properties. After discussing luminosity, temperature and chemical composition, we focus on the question whether massive stars may arrive at critical rotation during their evolution, either on the main sequence or in later stages. We find both cases to be possible and briefly discuss observable consequences.

In this paper, we investigate the evolutionary status of B[e]~stars from the point of view of stellar evolution theory. We try to answer to the question of how massive hot supergiants --- i.e. evolved stars --- can be capable of producing a circumstellar disk. We find and discuss three possibilities: very massive evolved main sequence stars close to critical rotation due to their proximity to their Eddington-limit, blue supergiants which have just left the red supergiant branch, and single star merger remnants of a close binary system. While the latter process seems to be required to understand the properties of the spectroscopic binary R4 in the LMC, the other two scenarios may be capable of explaining the distribution of the B[e] stars in the HR~diagram. The three scenarios make different predictions about the duration of the B[e]~phase, the time integrated disk mass and the stellar properties during the B[e]~phase, which may ultimately allow to distinguish them observationally.

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.

We investigate the pulsational properties of RSG models --- which we evolve from ZAMS masses in the range 10 to $20 \Msun$ --- by means of linear and non-linear calculations. We find period and growth rate of the dominant fundamental mode to increase with increasing luminosity-to-mass ratio $L/M$. Our models obtain relatively large $L/M$ values due to the inclusion of rotation in the evolutionary calculations; however, the largest values are obtained at and beyond central He-exhaustion due to major internal rearrangements of the nuclear burning regions. Our non-linear calculations as well as the behavior of the linear period and growth rate of the pulsations for periods approaching the Kelvin-Helmholtz time scale of the H-rich stellar envelope point towards the possibility of large amplitude pulsations. Such properties are similar to that found in AGB stars and suggest the possibility of a ``superwind'' to occur before the RSGs explode as supernovae. We conclude that changes in global stellar properties during the last few $10^4 \yr$ before core collapse may lead to drastic changes in the pulsational and wind properties of pre-supernova stars, with marked consequences for the immediate pre-supernova structure of the star and the circumstellar medium. We compare our results with observations of long-period OH/IR variables and discuss observational evidence for our scenario from observed supernova light curves, spectra and remnants.