Supernova Simulations from a 3D Progenitor Model -- Impact of Perturbations and Evolution of Explosion Properties

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.Comment: 23 pages, 22 figures, accepted for publication in MNRA

Evolution and Explosion of Very Massive Primordial Stars

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.Comment: 7 pages, including 4 figures; in. proc. MPA/ESO/MPE/USM Joint Astronomy Conference "Lighthouses of the Universe: The Most Luminous Celestial Objects and their use for Cosmology

The Supernova Channel of Super-AGB Stars

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).Comment: 13 pages, 16 figures, submitted to ApJ, uses emulateap

Nucleosynthesis in massive stars revisited

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.Comment: 4 pages, 1 figure; needs espcrc1.sty; talk at "Nuclei in the Cosmos 2000", Aarhus, Denmark, June 2000; will appear in Nucl. Phys.

Nuclear Aspects of Nucleosynthesis in Massive Stars

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.Comment: 6 pages, 2 figures; invited talk, to appear in the Proceedings of the Int. Conf. "Structure of the Nucleus at the Dawn of the Century", May 2000, Bologna, Ital

On the Maximum Mass of Accreting Primordial Supermassive Stars

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.Comment: 5 pages, 4 figures. Accepted ApJ letter

Why a Single-Star Model Cannot Explain the Bipolar Nebula of Eta Carinae

I examine the angular momentum evolution during the 1837-1856 Great Eruption of the massive star Eta Carinae. I find that the new estimate of the mass blown during that eruption implies that the envelope of Eta Car substantially spun-down during the 20 years eruption. Single-star models, most of which require the envelope to rotate close to the break-up velocity, cannot account for the bipolar nebula (the Homunculus) formed from matter expelled in that eruption. The kinetic energy and momentum of the Homunculus further constrains single-star models. I discuss how Eta Car can fit into a unified model for the formation of bipolar lobes where two oppositely ejected jets inflate two lobes (or bubbles). These jets are blown by an accretion disk, which requires stellar companions in the case of bipolar nebulae around stellar objects.Comment: ApJ, in press. New references and segments were adde

Nucleosynthesis in Massive Stars With Improved Nuclear and Stellar Physics

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.Comment: 53 pages, 17 color figures (3 as separate GIF images), slightly extended discussion and references, accepted by Ap

Long tails on thermonuclear X-ray bursts from neutron stars: a signature of inward heating?

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.Comment: Accepted for publication in Astronomy & Astrophysics, 12 pages, 14 figure

Massive Star Evolution: Nucleosynthesis and Nuclear Reaction Rate Uncertainties

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