Active Carbon and Oxygen Shell Burning Hydrodynamics

We have simulated 2.5$\times10^3$ s of the late evolution of a $23 \rm M_\odot$ star with full hydrodynamic behavior. We present the first simulations of a multiple-shell burning epoch, including the concurrent evolution and interaction of an oxygen and carbon burning shell. In addition, we have evolved a 3D model of the oxygen burning shell to sufficiently long times (300 s) to begin to assess the adequacy of the 2D approximation. We summarize striking new results: (1) strong interactions occur between active carbon and oxygen burning shells, (2) hydrodynamic wave motions in nonconvective regions, generated at the convective-radiative boundaries, are energetically important in both 2D and 3D with important consequences for compositional mixing, and (3) a spectrum of mixed p- and g-modes are unambiguously identified with corresponding adiabatic waves in these computational domains. We find that 2D convective motions are exaggerated relative to 3D because of vortex instability in 3D. We discuss the implications for supernova progenitor evolution and symmetry breaking in core collapse.Comment: 5 pages, 4 figures in emulateapj format. Accepted for publication in ApJ Letters. High resolution figure version available at http://spinach.as.arizona.ed

Collapsars - Gamma-Ray Bursts and Explosions in "Failed Supernovae"

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

On relative supernova rates and nucleosynthesis roles

It is shown that the Ni-56-Fe-56 observed in SN 1987A argues that core collapse supernovae may be responsible for more that 50 percent of the iron in the galaxy. Furthermore it is argued that the time averaged rate of thermonuclear driven Type I supernovae may be at least an order of magnitude lower than the average rate of core collapse supernovae. The present low rate of Type II supernovae (below their time averaged rate of approx. 1/10 yr) is either because the past rate was much higher because many core collapse supernovae are dim like SN 1987A. However, even in this latter case they are only an order of magnitude dimmer that normal Type II's due to the contribution of Ni-56 decay to the light curve

The Impact of Hydrodynamic Mixing on Supernova Progenitors

Recent multidimensional hydrodynamic simulations have demonstrated the importance of hydrodynamic motions in the convective boundary and radiative regions of stars to transport of energy, momentum, and composition. The impact of these processes increases with stellar mass. Stellar models which approximate this physics have been tested on several classes of observational problems. In this paper we examine the implications of the improved treatment on supernova progenitors. The improved models predict substantially different interior structures. We present pre-supernova conditions and simple explosion calculations from stellar models with and without the improved mixing treatment at 23 solar masses. The results differ substantially.Comment: 12 pages, 2 figures, accepted for publication in the Astrophysical Journal Letter

2D Multi-Angle, Multi-Group Neutrino Radiation-Hydrodynamic Simulations of Postbounce Supernova Cores

We perform axisymmetric (2D) multi-angle, multi-group neutrino radiation-hydrodynamic calculations of the postbounce phase of core-collapse supernovae using a genuinely 2D discrete-ordinate (S_n) method. We follow the long-term postbounce evolution of the cores of one nonrotating and one rapidly-rotating 20-solar-mass stellar model for ~400 milliseconds from 160 ms to ~550 ms after bounce. We present a multi-D analysis of the multi-angle neutrino radiation fields and compare in detail with counterpart simulations carried out in the 2D multi-group flux-limited diffusion (MGFLD) approximation to neutrino transport. We find that 2D multi-angle transport is superior in capturing the global and local radiation-field variations associated with rotation-induced and SASI-induced aspherical hydrodynamic configurations. In the rotating model, multi-angle transport predicts much larger asymptotic neutrino flux asymmetries with pole to equator ratios of up to ~2.5, while MGFLD tends to sphericize the radiation fields already in the optically semi-transparent postshock regions. Along the poles, the multi-angle calculation predicts a dramatic enhancement of the neutrino heating by up to a factor of 3, which alters the postbounce evolution and results in greater polar shock radii and an earlier onset of the initially rotationally weakened SASI. In the nonrotating model, differences between multi-angle and MGFLD calculations remain small at early times when the postshock region does not depart significantly from spherical symmetry. At later times, however, the growing SASI leads to large-scale asymmetries and the multi-angle calculation predicts up to 30% higher average integral neutrino energy deposition rates than MGFLD.Comment: 20 pages, 21 figures. Minor revisions. Accepted for publication in ApJ. A version with high-resolution figures may be obtained from http://www.stellarcollapse.org/papers/Ott_et_al2008_multi_angle.pd

Massive Stars in the Range 13−25M⊙\rm 13-25 M_\odot: Evolution and Nucleosynthesis. II. the Solar Metallicity Models

We present the evolutionary properties of a set of massive stellar models (namely 13, 15, 20 and 25 $\rm M_\odot$) from the main sequence phase up to the onset of the iron core collapse. All these models have initial solar chemical composition, i.e. Y=0.285 and Z=0.02. A 179 isotope network, extending from neutron up to $\rm ^{68}Zn$ and fully coupled to the evolutionary code has been adopted from the Carbon burning onward. Our results are compared, whenever possible, to similar computations available in literature.Comment: 42 pages, 18 figures, 26 tables, accepted for publicatin in ApJ

A dynamical model of surrogate reactions

A new dynamical model is developed to describe the whole process of surrogate reactions; transfer of several nucleons at an initial stage, thermal equilibration of residues leading to washing out of shell effects and decay of populated compound nuclei are treated in a unified framework. Multi-dimensional Langevin equations are employed to describe time-evolution of collective coordinates with a time-dependent potential energy surface corresponding to different stages of surrogate reactions. The new model is capable of calculating spin distributions of the compound nuclei, one of the most important quantity in the surrogate technique. Furthermore, various observables of surrogate reactions can be calculated, e.g., energy and angular distribution of ejectile, and mass distributions of fission fragments. These features are important to assess validity of the proposed model itself, to understand mechanisms of the surrogate reactions and to determine unknown parameters of the model. It is found that spin distributions of compound nuclei produced in $^{18}$O+$^{238}$U $\rightarrow ^{16}$O+$^{240*}$U and $^{18}$O+$^{236}$U $\rightarrow ^{16}$O+$^{238*}$U reactions are equivalent and much less than 10$\hbar$, therefore satisfy conditions proposed by Chiba and Iwamoto (PRC 81, 044604(2010)) if they are used as a pair in the surrogate ratio method.Comment: 17 pages, 5 figure

Observational Tests and Predictive Stellar Evolution II: Non-standard Models

We examine contributions of second order physical processes to results of stellar evolution calculations amenable to direct observational testing. In the first paper in the series (Young et al. 2001) we established baseline results using only physics which are common to modern stellar evolution codes. In the current paper we establish how much of the discrepancy between observations and baseline models is due to particular elements of new physics. We then consider the impact of the observational uncertainties on the maximum predictive accuracy achievable by a stellar evolution code. The sun is an optimal case because of the precise and abundant observations and the relative simplicity of the underlying stellar physics. The Standard Model is capable of matching the structure of the sun as determined by helioseismology and gross surface observables to better than a percent. Given an initial mass and surface composition within the observational errors, and no additional constraints for which the models can be optimized, it is not possible to predict the sun's current state to better than ~7%. Convectively induced mixing in radiative regions, seen in multidimensional hydrodynamic simulations, dramatically improves the predictions for radii, luminosity, and apsidal motions of eclipsing binaries while simultaneously maintaining consistency with observed light element depletion and turnoff ages in young clusters (Young et al. 2003). Systematic errors in core size for models of massive binaries disappear with more complete mixing physics, and acceptable fits are achieved for all of the binaries without calibration of free parameters. The lack of accurate abundance determinations for binaries is now the main obstacle to improving stellar models using this type of test.Comment: 33 pages, 8 figures, accepted for publication in the Astrophysical Journa

Nucleosynthesis in Type Ia Supernovae

Among the major uncertainties involved in the Chandrasekhar mass models for Type Ia supernovae are the companion star of the accreting white dwarf (or the accretion rate that determines the carbon ignition density) and the flame speed after ignition. We present nucleosynthesis results from relatively slow deflagration (1.5 - 3 % of the sound speed) to constrain the rate of accretion from the companion star. Because of electron capture, a significant amount of neutron-rich species such as ^{54}Cr, ^{50}Ti, ^{58}Fe, ^{62}Ni, etc. are synthesized in the central region. To avoid the too large ratios of ^{54}Cr/^{56}Fe and ^{50}Ti/^{56}Fe, the central density of the white dwarf at thermonuclear runaway must be as low as \ltsim 2 \e9 \gmc. Such a low central density can be realized by the accretion as fast as \dot M \gtsim 1 \times 10^{-7} M_\odot yr^{-1}. These rapidly accreting white dwarfs might correspond to the super-soft X-ray sources.Comment: 10 page LaTeX, 7 PostScript figures, to appear in Nuclear Physics A, Vol. A621 (1997

Collisional Dark Matter and the Origin of Massive Black Holes

If the cosmological dark matter is primarily in the form of an elementary particle which has cross section and mass for self-interaction having a ratio similar to that of ordinary nuclear matter, then seed black holes (formed in stellar collapse) will grow in a Hubble time, due to accretion of the dark matter, to a mass range 10^6 - 10^9 solar masses. Furthermore, the dependence of the final black hole mass on the galaxy velocity dispersion will be approximately as observed and the growth rate will show a time dependence consistent with observations. Other astrophysical consequences of collisional dark matter and tests of the idea are noted.Comment: 7 pages, no figures, LaTeX2e, Accepted for publication in Phys. Rev. Lett. Changed conten