We explore the potential of current and future liquid scintillator neutrino detectors of O(10) kt mass to localize a presupernova neutrino signal in the sky. In the hours preceding the core collapse of a nearby star (at distance D < 1 kpc), tens to hundreds of inverse beta decay events will be recorded, and their reconstructed topology in the detector can be used to estimate the direction to the star. Although the directionality of inverse beta decay is weak (~8% forward-backward asymmetry for currently available liquid scintillators), we find that for a fiducial signal of 200 events (which is realistic for Betelgeuse), a positional error of ~60 degrees can be achieved, resulting in the possibility to narrow the list of potential stellar candidates to less than ten, typically. For a configuration with improved forward-backward asymmetry (~40%, as expected for a lithium-loaded liquid scintillator), the angular sensitivity improves to ~15 degrees, and - when a distance upper limit is obtained from the overall event rate - it is in principle possible to uniquely identify the progenitor star. Any localization information accompanying an early supernova alert will be useful to multi-messenger observations and to particle physics tests using collapsing stars.

We explore the evolution of a select grid of solar metallicity stellar models from their pre-main sequence phase to near their final fates in a neutrino Hertzsprung-Russell diagram, where the neutrino luminosity replaces the traditional photon luminosity. Using a calibrated \MESA\ solar model for the solar neutrino luminosity ($L_{\nu,\odot}$ = 0.02398 $\cdot$ $L_{\gamma,\odot}$ = 9.1795 $\times$ 10$^{31}$ erg s$^{-1}$) as a normalization, we identify $\simeq$ 0.3 MeV electron neutrino emission from helium burning during the helium flash (peak $L_{\nu} / L_{\nu,\odot} \simeq$ 10$^4$, flux $\Phi_{\nu, {\rm He \ flash}} \simeq$ 170 (10 pc/$d$)$^{2}$ cm$^{-2}$ s$^{-1}$ for a star located at a distance of $d$ parsec, timescale $\simeq$ 3 days) and the thermal pulse (peak $L_{\nu} / L_{\nu,\odot} \simeq$ 10$^9$, flux $\Phi_{\nu, {\rm TP}} \simeq$ 1.7$\times$10$^7$ (10 pc/$d$)$^{2}$ cm$^{-2}$ s$^{-1}$, timescale $\simeq$ 0.1 yr) phases of evolution in low mass stars as potential probes for stellar neutrino astronomy. We also delineate the contribution of neutrinos from nuclear reactions and thermal processes to the total neutrino loss along the stellar tracks in a neutrino Hertzsprung-Russell diagram. We find, broadly but with exceptions, that neutrinos from nuclear reactions dominate whenever hydrogen and helium burn, and that neutrinos from thermal processes dominate otherwise.

The collapse of degenerate oxygen-neon cores (i.e., electron-capture supernovae or accretion-induced collapse) proceeds through a phase in which a deflagration wave ("flame") forms at or near the center and propagates through the star. In models, the assumed speed of this flame influences whether this process leads to an explosion or to the formation of a neutron star. We calculate the laminar flame speeds in degenerate oxygen-neon mixtures with compositions motivated by detailed stellar evolution models. These mixtures include trace amounts of carbon and have a lower electron fraction than those considered in previous work. We find that trace carbon has little effect on the flame speeds, but that material with electron fraction $Y_e \approx 0.48-0.49$ has laminar flame speeds that are $\approx 2$ times faster than those at $Y_e = 0.5$. We provide tabulated flame speeds and a corresponding fitting function so that the impact of this difference can be assessed via full star hydrodynamical simulations of the collapse process.

Astronomy, as a field, has long encouraged the development of free, open digital infrastructure (e.g., National Research Council 2010, 2011). Examples range from simple scripts that enable individual scientific research, through software instruments for entire communities, to data reduction pipelines for telescope operations at national facilities. As with the digital infrastructure of our larger society today (e.g., Eghbal 2016), nearly all astronomical research relies on free, open source software (FOSS) written and maintained by a small number of developers. And like the physical infrastructure of roads or bridges, digital infrastructure needs regular upkeep and maintenance (e.g., Eghbal 2016). In astronomy, financial support for maintaining existing digital infrastructure is generally much harder to secure than funding for developing new digital infrastructures that promise new science. Sustaining astronomy's digital infrastructure is a new topic for many, the sustainability challenges are not always widely known...

We update the capabilities of the open-knowledge software instrument Modules for Experiments in Stellar Astrophysics (MESA). RSP is a new functionality in MESAstar that models the non-linear radial stellar pulsations that characterize RR Lyrae, Cepheids, and other classes of variable stars. We significantly enhance numerical energy conservation capabilities, including during mass changes. For example, this enables calculations through the He flash that conserve energy to better than 0.001 %. To improve the modeling of rotating stars in MESA, we introduce a new approach to modifying the pressure and temperature equations of stellar structure, and a formulation of the projection effects of gravity darkening. A new scheme for tracking convective boundaries yields reliable values of the convective-core mass, and allows the natural emergence of adiabatic semiconvection regions during both core hydrogen- and helium-burning phases. We quantify the parallel performance of MESA on current generation multicore architectures and demonstrate improvements in the computational efficiency of radiative levitation. We report updates to the equation of state and nuclear reaction physics modules. We briefly discuss the current treatment of fallback in core-collapse supernova models and the thermodynamic evolution of supernova explosions. We close by discussing the new MESA Testhub software infrastructure to enhance source-code development.

Recent near-IR power-spectra and panchromatic Extragalactic Background Light measurements provide upper limits on the near-IR surface brightness (SB>31 mag/arcsec^2) that may come from Pop III stars and accretion disks around resulting stellar-mass black holes (BHs) in the epoch of First Light (z=7-17). Physical parameters for zero metallicity Pop III stars at z>7 can be estimated from MESA stellar evolution models through helium-depletion, and for BH accretion disks from quasar microlensing results and multicolor accretion models. Second-generation stars can form at higher multiplicity, so that BH accretion disks may be fed by Roche-lobe overflow from lower-mass companions in their AGB stage. The near-IR SB constraints can be used to calculate the number of caustic transits behind lensing clusters that JWST and the 25~39 m ground-based telescopes may detect for both Pop III stars and stellar mass BH accretion disks. Because Pop III stars and stellar mass BH accretion disks have sizes of a few x 10^-11 arcsec at z>7, typical caustic magnifications can be mu=10^4~10^5, with rise times of hours and decline times of < 1 year for cluster transverse velocities of v<1000 km/s. Microlensing by intracluster medium objects can modify transit magnifications, and lengthen visibility times. Depending on BH masses, accretion-disk radii and feeding efficiencies, stellar-mass BH accretion-disk caustic transits could outnumber those from Pop III stars. To observe Pop III caustic transits directly may require monitoring 3~30 lensing clusters to AB< 29 mag over a decade or more. Such a program must be started with JWST in Cycle 1, and -- depending on the role of microlensing in the Intra Cluster Light -- should be continued for decades with the GMT and TMT, where JWST and the ground-based telescopes each will play a unique and strongly complementary role.

The identification of the electromagnetic counterpart of the gravitational wave event, GW170817, and discovery of neutrinos and gamma-rays from TXS 0506+056 heralded the new era of multi-messenger astrophysics. As the number of multi-messenger events rapidly grow over the next decade, the cyberinfrastructure requirements to handle the increase in data rates, data volume, need for event follow up, and analysis across the different messengers will also explosively grow. The cyberinfrastructure requirements to enhance multi-messenger astrophysics will both be a major challenge and opportunity for astronomers, physicists, computer scientists and cyberinfrastructure specialists. Here we outline some of these requirements and argue for a distributed cyberinfrastructure institute for multi-messenger astrophysics to meet these challenges.

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.

About two generations ago, a large part of AMO science was dominated by experimental high energy collision studies and perturbative theoretical methods. Since then, AMO science has undergone a transition and is now dominated by quantum, ultracold, and ultrafast studies. But in the process, the field has passed over the complexity that lies between these two extremes. Most of the Universe resides in this intermediate region. We put forward that the next frontier for AMO science is to explore the AMO complexity that describes most of the Cosmos.

The discovery via gravitational waves of binary black hole systems with total masses greater than $60M_\odot$ has raised interesting questions for stellar evolution theory. Among the most promising formation channels for these systems is one involving a common envelope binary containing a low metallicity, core helium burning star with mass $\sim 80-90M_\odot$ and a black hole with mass $\sim 30-40M_\odot$. For this channel to be viable, the common envelope binary must eject more than half the giant star's mass and reduce its orbital separation by as much as a factor of 80. We discuss issues faced in numerically simulating the common envelope evolution of such systems and present a 3D AMR simulation of the dynamical inspiral of a low-metallicity red supergiant with a massive black hole companion.

KIC 08626021 is a pulsating DB white dwarf of considerable recent interest, and first of its class to be extensively monitored by Kepler for its pulsation properties. Fitting the observed oscillation frequencies of KIC 08626021 to a model can yield insights into its otherwise-hidden internal structure. Template-based white dwarf models choose a luminosity profile where the luminosity is proportional to the enclosed mass, $L_r \propto M_r$, independent of the effective temperature $T_{\rm eff}$. Evolutionary models of young white dwarfs with $T_{\rm eff} \gtrsim$ 25,000 K suggest neutrino emission gives rise to luminosity profiles with $L_r$ $\not\propto$ $M_r$. We explore this contrast by comparing the oscillation frequencies between two nearly identical white dwarf models: one with an enforced $L_r \propto M_r$ luminosity profile and the other with a luminosity profile determined by the star's previous evolution history. We find the low order g-mode frequencies differ by up to $\simeq$ 70 $\mu$Hz over the range of Kepler observations for KIC 08626021. This suggests that by neglecting the proper thermal structure of the star (e.g., accounting for the effect of plasmon neutrino losses), the model frequencies calculated by using an $L_r \propto M_r$ profile may have uncorrected, effectively-random errors at the level of tens of $\mu$Hz. A mean frequency difference of 30 $\mu$Hz, based on linearly extrapolating published results, suggests a template model uncertainty in the fit precision of $\simeq$ 12% in white dwarf mass, $\simeq$ 9% in the radius, and $\simeq$ 3% in the central oxygen mass fraction.

Accurately reproducing the physics behind the detonations of Type Ia supernovae and the resultant nucleosynthetic yields is important for interpreting observations of spectra and remnants. The scales of the processes involved span orders of magnitudes, making the problem computationally impossible to ever fully resolve in full star simulations in the present and near future. In the lower density regions of the star, the curvature of the detonation front will slow the detonation, affecting the production of intermediate mass elements. We find that shock strengthening due to the density gradient present in the outer layers of the progenitor is essential for understanding the nucleosynthesis there, with burning extending well below the density at which a steady-state detonation is extinct. We show that a complete reaction network is not sufficient to obtain physical detonations at high densities and modest resolution due to numerical mixing at the unresolved reaction front. At low densities, below 6$\times$10$^{5}$ g cm$^{-3}$, it is possible to achieve high enough resolution to separate the shock and the reaction region,and the abundance structure predicted by fully resolved quasi-steady-state calculations is obtained. For our best current benchmark yields, we utilize a method in which the unresolved portion of Lagrangian histories are reconstructed based on fully resolved quasi-steady-state detonation calculations. These computations demonstrate that under-resolved simulations agree approximately, $\sim$10\% in post-shock values of temperature, pressure, density, and abundances, with expected detonation structures sufficiently far from the under-resolved region, but that there is still room for some improvement in the treatment of subgrid reactions in the hydrodynamics to before better than 1$\%$ can be achieved at all densities.

We summarize panchromatic Extragalactic Background Light data to place upper limits on the integrated near-infrared surface brightness (SB) that may come from Population III stars and possible accretion disks around their stellar-mass black holes (BHs) in the epoch of First Light, broadly taken from z$\simeq$7-17. Theoretical predictions and recent near-infrared power-spectra provide tighter constraints on their sky-signal. We outline the physical properties of zero metallicity Population III stars from MESA stellar evolution models through helium-depletion and of BH accretion disks at z$\gtrsim$7. We assume that second-generation non-zero metallicity stars can form at higher multiplicity, so that BH accretion disks may be fed by Roche-lobe overflow from lower-mass companions. We use these near-infrared SB constraints to calculate the number of caustic transits behind lensing clusters that the James Webb Space Telescope and the next generation ground-based telescopes may observe for both Population III stars and their BH accretion disks. Typical caustic magnifications can be $\mu$$\simeq$10$^4$-10$^5$, with rise times of hours and decline times of $\lesssim$1 year for cluster transverse velocities of $v_{T}$$\lesssim$1000 km s$^{-1}$. Microlensing by intracluster medium objects can modify transit magnifications, but lengthen visibility times. Depending on BH masses, accretion-disk radii and feeding efficiencies, stellar-mass BH accretion-disk caustic transits could outnumber those from Population III stars. To observe Population III caustic transits directly may require to monitor 3-30 lensing clusters to AB$\lesssim$29 mag over a decade.

We update the capabilities of the software instrument Modules for Experiments in Stellar Astrophysics (MESA) and enhance its ease of use and availability. Our new approach to locating convective boundaries is consistent with the physics of convection, and yields reliable values of the convective core mass during both hydrogen and helium burning phases. Stars with $M<8\,{\rm M_\odot}$ become white dwarfs and cool to the point where the electrons are degenerate and the ions are strongly coupled, a realm now available to study with MESA due to improved treatments of element diffusion, latent heat release, and blending of equations of state. Studies of the final fates of massive stars are extended in MESA by our addition of an approximate Riemann solver that captures shocks and conserves energy to high accuracy during dynamic epochs. We also introduce a 1D capability for modeling the effects of Rayleigh-Taylor instabilities that, in combination with the coupling to a public version of the STELLA radiation transfer instrument, creates new avenues for exploring Type II supernovae properties. These capabilities are exhibited with exploratory models of pair-instability supernova, pulsational pair-instability supernova, and the formation of stellar mass black holes. The applicability of MESA is now widened by the capability of importing multi-dimensional hydrodynamic models into MESA. We close by introducing software modules for handling floating point exceptions and stellar model optimization, and four new software tools -- MESAWeb, MESA-Docker, pyMESA, and mesastar.org -- to enhance MESA's education and research impact.

We explore properties of core-collapse supernova progenitors with respect to the composite uncertainties in the thermonuclear reaction rates by coupling the reaction rate probability density functions provided by the STARLIB reaction rate library with $\texttt{MESA}$ stellar models. We evolve 1000 15 $M_{\odot}$ models from the pre main-sequence to core O-depletion at solar and subsolar metallicities for a total of 2000 Monte Carlo stellar models. For each stellar model, we independently and simultaneously sample 665 thermonuclear reaction rates and use them in a $\texttt{MESA}$ in situ reaction network that follows 127 isotopes from $^{1}$H to $^{64}$Zn. With this framework we survey the core mass, burning lifetime, composition, and structural properties at five different evolutionary epochs. At each epoch we measure the probability distribution function of the variations of each property and calculate Spearman Rank-Order Correlation coefficients for each sampled reaction rate to identify which reaction rate has the largest impact on the variations on each property. We find that uncertainties in $^{14}$N$(p,\gamma)^{15}$O, triple-$\alpha$, $^{12}$C$(\alpha,\gamma)^{16}$O, $^{12}$C($^{12}$C,$p$)$^{23}$Na, $^{12}$C($^{16}$O,$p$)$^{27}$Al, $^{16}$O($^{16}$O,$n$)$^{31}$S, $^{16}$O($^{16}$O,$p$)$^{31}$P, and $^{16}$O($^{16}$O,$\alpha$)$^{28}$Si reaction rates dominate the variations of the properties surveyed. We find that variations induced by uncertainties in nuclear reaction rates grow with each passing phase of evolution, and at core H-, He-depletion are of comparable magnitude to the variations induced by choices of mass resolution and network resolution. However, at core C-, Ne-, and O-depletion, the reaction rate uncertainties can dominate the variation causing uncertainty in various properties of the stellar model in the evolution towards iron core-collapse.

We present a new calculation of the neutrino flux received at Earth from a massive star in the $\sim 24$ hours of evolution prior to its explosion as a supernova (presupernova). Using the stellar evolution code MESA, the neutrino emissivity in each flavor is calculated at many radial zones and time steps. In addition to thermal processes, neutrino production via beta processes is modeled in detail, using a network of 204 isotopes. We find that the total produced $\nu_{e}$ flux has a high energy spectrum tail, at $E \gtrsim 3 - 4$ MeV, which is mostly due to decay and electron capture on isotopes with $A = 50 - 60$. In a tentative window of observability of $E \gtrsim 0.5$ MeV and $t < 2$ hours pre-collapse, the contribution of beta processes to the $\nu_{e}$ flux is at the level of $\sim90\%$ . For a star at $D=1$ kpc distance, a 17 kt liquid scintillator detector would typically observe several tens of events from a presupernova, of which up to $\sim 30\%$ due to beta processes. These processes dominate the signal at a liquid argon detector, thus greatly enhancing its sensitivity to a presupernova.

The physical process whereby a carbon--oxygen white dwarf explodes as a Type Ia supernova (SN Ia) remains highly uncertain. The degree of neutronization in SN Ia ejecta holds clues to this process because it depends on the mass and the metallicity of the stellar progenitor, and on the thermodynamic history prior to the explosion. We report on a new method to determine ejecta neutronization using Ca and S lines in the X-ray spectra of Type Ia supernova remnants (SNRs). Applying this method to \textit{Suzaku} data of Tycho, Kepler, 3C 397 and G337.2$-$0.7 in the Milky Way, and N103B in the Large Magellanic Cloud, we find that the neutronization of the ejecta in N103B is comparable to that of Tycho and Kepler, which suggests that progenitor metallicity is not the only source of neutronization in SNe Ia. We then use a grid of SN Ia explosion models to infer the metallicities of the stellar progenitors of our SNRs. The implied metallicities of 3C 397, G337.2$-$0.7, and N103B are major outliers compared to the local stellar metallicity distribution functions, indicating that progenitor metallicity can be ruled out as the origin of neutronization for these SNRs. Although the relationship between ejecta neutronization and equivalent progenitor metallicity is subject to uncertainties stemming from the $^{12}$C$\,$+$^{16}$O reaction rate, which affects the Ca/S mass ratio, our main results are not sensitive to these details.

Recent Suzaku X-ray spectra of SNR 3C 397 indicate enhanced stable iron-group element abundances of Ni, Mn, Cr, and Fe. Seeking to address key questions about the progenitor and explosion mechanism of 3C 397, we compute nucleosynthetic yields from a suite of multidimensional hydrodynamics models in the near-Chandrasekhar mass, single-degenerate paradigm for supernova Type Ia. Varying the progenitor white dwarf internal structure, composition, ignition, and explosion mechanism, we find the best match to the observed iron-peak elements of 3C 397 are dense (central density $\ge$ 6$\times$10$^{9}$ g cm$^{-3}$), low-carbon white dwarfs that undergo a weak, centrally-ignited deflagration, followed by a subsequent detonation. The amount of $^{56}$Ni produced is consistent with a normal or bright normal supernova Type Ia. A pure deflagration of a centrally-ignited, low central density ($\simeq$ 2$\times$10$^{9}$ g cm$^{-3}$) progenitor white dwarf, frequently considered in the literature, is also found to produce good agreement with 3C 397 nucleosynthetic yields, but leads to a subluminous SN Ia event, in conflict with X-ray linewidth data. Additionally, in contrast to prior work which suggested a large super-solar metallicity for the white dwarf progenitor for SNR 3C 397, we find satisfactory agreement for solar and sub-solar metallicity progenitors. We discuss a range of implications our results have for the single-degenerate channel.

We substantially update the capabilities of the open-source software instrument Modules for Experiments in Stellar Astrophysics (MESA). MESA can now simultaneously evolve an interacting pair of differentially rotating stars undergoing transfer and loss of mass and angular momentum, greatly enhancing the prior ability to model binary evolution. New MESA capabilities in fully coupled calculation of nuclear networks with hundreds of isotopes now allow MESA to accurately simulate advanced burning stages needed to construct supernova progenitor models. Implicit hydrodynamics with shocks can now be treated with MESA, enabling modeling of the entire massive star lifecycle, from pre-main sequence evolution to the onset of core collapse and nucleosynthesis from the resulting explosion. Coupling of the GYRE non-adiabatic pulsation instrument with MESA allows for new explorations of the instability strips for massive stars while also accelerating the astrophysical use of asteroseismology data. We improve treatment of mass accretion, giving more accurate and robust near-surface profiles. A new MESA capability to calculate weak reaction rates "on-the-fly" from input nuclear data allows better simulation of accretion induced collapse of massive white dwarfs and the fate of some massive stars. We discuss the ongoing challenge of chemical diffusion in the strongly coupled plasma regime, and exhibit improvements in MESA that now allow for the simulation of radiative levitation of heavy elements in hot stars. We close by noting that the MESA software infrastructure provides bit-for-bit consistency for all results across all the supported platforms, a profound enabling capability for accelerating MESA's development.

We explore the variation in single star 15-30 $\rm{M}_{\odot}$, non-rotating, solar metallicity, pre-supernova MESA models due to changes in the number of isotopes in a fully-coupled nuclear reaction network and adjustments in the mass resolution. Within this two-dimensional plane we quantitatively detail the range of core masses at various stages of evolution, mass locations of the main nuclear burning shells, electron fraction profiles, mass fraction profiles, burning lifetimes, stellar lifetimes, and compactness parameter at core-collapse for models with and without mass loss. Up to carbon burning we generally find mass resolution has a larger impact on the variations than the number of isotopes, while the number of isotopes plays a more significant role in determining the span of the variations for neon, oxygen and silicon burning. Choice of mass resolution dominates the variations in the structure of the intermediate convection zone and secondary convection zone during core and shell hydrogen burning respectively, where we find a minimum mass resolution of $\approx$ 0.01 $\rm{M}_{\odot}$ is necessary to achieve convergence in the helium core mass at the $\approx$5% level. On the other hand, at the onset of core-collapse we find $\approx$30% variations in the central electron fraction and mass locations of the main nuclear burning shells, a minimum of $\approx$127 isotopes is needed to attain convergence of these values at the $\approx$10% level.

It has been proposed that mixing induced by convective overshoot can disrupt the inward propagation of carbon deflagrations in super-asymptotic giant branch stars. To test this theory, we study an idealized model of convectively bounded carbon flames with 3D hydrodynamic simulations of the Boussinesq equations using the pseudospectral code Dedalus. Because the flame propagation timescale is much longer than the convection timescale, we approximate the flame as fixed in space, and only consider its effects on the buoyancy of the fluid. By evolving a passive scalar field, we derive a {\it turbulent} chemical diffusivity produced by the convection as a function of height, $D_{\rm t}(z)$. Convection can stall a flame if the chemical mixing timescale, set by the turbulent chemical diffusivity, $D_{\rm t}$, is shorter than the flame propagation timescale, set by the thermal diffusivity, $\kappa$, i.e., when $D_{\rm t}>\kappa$. However, we find $D_{\rm t}<\kappa$ for most of the flame because convective plumes are not dense enough to penetrate into the flame. Extrapolating to realistic stellar conditions, this implies that convective mixing cannot stall a carbon flame and that "hybrid carbon-oxygen-neon" white dwarfs are not a typical product of stellar evolution.

The problem of estimating non-resonant astrophysical S-factors and thermonuclear reaction rates, based on measured nuclear cross sections, is of major interest for nuclear energy generation, neutrino physics, and element synthesis. Many different methods have been applied in the past to this problem, almost all of them based on traditional statistics. Bayesian methods, on the other hand, are now in widespread use in the physical sciences. In astronomy, for example, Bayesian statistics is applied to the observation of extra-solar planets, gravitational waves, and type Ia supernovae. However, nuclear physics, in particular, has been slow to adopt Bayesian methods. We present astrophysical S-factors and reaction rates based on Bayesian statistics. We develop a framework that incorporates robust parameter estimation, systematic effects, and non-Gaussian uncertainties in a consistent manner. The method is applied to the d(p,$\gamma$)$^3$He, $^3$He($^3$He,2p)$^4$He, and $^3$He($\alpha$,$\gamma$)$^7$Be reactions, important for deuterium burning, solar neutrinos, and big bang nucleosynthesis.

Computational skills are required across all astronomy disciplines. Many students enter degree programs without sufficient skills to solve computational problems in their core classes or contribute immediately to research. We recommend advocacy for computational literacy, familiarity with fundamental software carpentry skills, and mastery of basic numerical methods by the completion of an undergraduate degree in Astronomy. We recommend the AAS Education Task Force advocate for a significant increase in computational literacy. We encourage the AAS to modestly fund efforts aimed at providing Open Education Resources (OER) that will significantly impact computational literacy in astronomy education.

We refine our previously introduced parameterized model for explosive carbon-oxygen fusion during thermonuclear supernovae (SN Ia) by adding corrections to post-processing of recorded Lagrangian fluid element histories to obtain more accurate isotopic yields. Deflagration and detonation products are verified for propagation in a uniform density medium. A new method is introduced for reconstructing the temperature-density history within the artificially thick model deflagration front. We obtain better than 5\% consistency between the electron capture computed by the burning model and yields from post-processing. For detonations, we compare to a benchmark calculation of the structure of driven steady-state planar detonations performed with a large nuclear reaction network and error-controlled integration. We verify that, for steady-state planar detonations down to a density of 5x10^6 g/cc, our post processing matches the major abundances in the benchmark solution typically to better than 10% for times greater than 0.01 s after the shock front passage. As a test case to demonstrate the method, presented here with post-processing for the first time, we perform a two dimensional simulation of a SN Ia in the Chandrasekhar-mass deflagration-detonation transition (DDT) scenario. We find that reconstruction of deflagration tracks leads to slightly more complete silicon burning than without reconstruction. The resulting abundance structure of the ejecta is consistent with inferences from spectroscopic studies of observed SNe Ia. We confirm the absence of a central region of stable Fe-group material for the multi-dimensional DDT scenario. Detailed isotopic yields are tabulated and only change modestly when using deflagration reconstruction.

In Type Ia Supernovae (\sneia), the relative abundances of chemical elements are affected by the neutron excess in the composition of the progenitor white dwarf. Since these products leave signatures in the spectra near maximum light, spectral features may be used to constrain the composition of the progenitor. We calculate the nucleosynthetic yields for three \snia simulations, assuming single degenerate, Chandrasekhar mass progenitors, for a wide range of progenitor metallicities, and calculate synthetic light curves and spectra to explore correlations between progenitor metallicity and the strength of spectral features. We use two 2D simulations of the deflagration-detonation-transition scenario with different $^{56}$Ni yields and the W7 simulation to control for differences between explosion models and total yields. While the overall yields of intermediate mass elements (16 $<$ A $\leq$ 40) differ between the three cases, trends in the yields are similar. With increasing metallicity, $^{28}$Si yields remain nearly constant, $^{40}$Ca yields decline, and Ti and $^{54}$Fe yields increase. In the synthetic spectra, we identify two features at 30 days post explosion that appear to deepen with progenitor metallicity: a Ti feature around 4200\,\AA\ and a Fe feature around 5200\,\AA\@. In all three simulations, their pseudo equivalent widths show a systematic trend with progenitor metallicity. This suggests that these two features may allow differentiation among progenitor metallicities of observed \sneia and potentially help reduce the intrinsic Hubble scatter.

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.

This white paper informs the nuclear astrophysics community and funding agencies about the scientific directions and priorities of the field and provides input from this community for the 2015 Nuclear Science Long Range Plan. It summarizes the outcome of the nuclear astrophysics town meeting that was held on August 21-23, 2014 in College Station at the campus of Texas A&M University in preparation of the NSAC Nuclear Science Long Range Plan. It also reflects the outcome of an earlier town meeting of the nuclear astrophysics community organized by the Joint Institute for Nuclear Astrophysics (JINA) on October 9- 10, 2012 Detroit, Michigan, with the purpose of developing a vision for nuclear astrophysics in light of the recent NRC decadal surveys in nuclear physics (NP2010) and astronomy (ASTRO2010). The white paper is furthermore informed by the town meeting of the Association of Research at University Nuclear Accelerators (ARUNA) that took place at the University of Notre Dame on June 12-13, 2014. In summary we find that nuclear astrophysics is a modern and vibrant field addressing fundamental science questions at the intersection of nuclear physics and astrophysics. These questions relate to the origin of the elements, the nuclear engines that drive life and death of stars, and the properties of dense matter. A broad range of nuclear accelerator facilities, astronomical observatories, theory efforts, and computational capabilities are needed. With the developments outlined in this white paper, answers to long standing key questions are well within reach in the coming decade.

We investigate properties of carbon-oxygen white dwarfs with respect to the composite uncertainties in the reaction rates using the stellar evolution toolkit, Modules for Experiments in Stellar Astrophysics (MESA) and the probability density functions in the reaction rate library STARLIB. These are the first Monte Carlo stellar evolution studies that use complete stellar models. Focusing on 3 M$_{\odot}$ models evolved from the pre main-sequence to the first thermal pulse, we survey the remnant core mass, composition, and structure properties as a function of 26 STARLIB reaction rates covering hydrogen and helium burning using a Principal Component Analysis and Spearman Rank-Order Correlation. Relative to the arithmetic mean value, we find the width of the 95\% confidence interval to be $\Delta M_{{\rm 1TP}}$ $\approx$ 0.019 M$_{\odot}$ for the core mass at the first thermal pulse, $\Delta$$t_{\rm{1TP}}$ $\approx$ 12.50 Myr for the age, $\Delta \log(T_{{\rm c}}/{\rm K}) \approx$ 0.013 for the central temperature, $\Delta \log(\rho_{{\rm c}}/{\rm g \ cm}^{-3}) \approx$ 0.060 for the central density, $\Delta Y_{\rm{e,c}} \approx$ 2.6$\times$10$^{-5}$ for the central electron fraction, $\Delta X_{\rm c}(^{22}\rm{Ne}) \approx$ 5.8$\times$10$^{-4}$, $\Delta X_{\rm c}(^{12}\rm{C}) \approx$ 0.392, and $\Delta X_{\rm c}(^{16}\rm{O}) \approx$ 0.392. Uncertainties in the experimental $^{12}$C($\alpha,\gamma)^{16}\rm{O}$, triple-$\alpha$, and $^{14}$N($p,\gamma)^{15}\rm{O}$ reaction rates dominate these variations. We also consider a grid of 1 to 6 M$_{\odot}$ models evolved from the pre main-sequence to the final white dwarf to probe the sensitivity of the initial-final mass relation to experimental uncertainties in the hydrogen and helium reaction rates.

We explore the effects of rotation on convective carbon, oxygen, and silicon shell burning during the late stages of evolution in a 20Msun star. Using the Modules for Experiments in Stellar Astrophysics (MESA) we construct 1D stellar models both with no rotation and with an initial rigid rotation of 50% of critical. At different points during the evolution, we map the 1D models into 2D and follow the multidimensional evolution using the FLASH compressible hydrodynamics code for many convective turnover times until a quasi-steady state is reached. We characterize the strength and scale of convective motions via decomposition of the momentum density into vector spherical harmonics. We find that rotation influences the total power in solenoidal modes, with a slightly larger impact for carbon and oxygen shell burning than for silicon shell burning. Including rotation in one-dimensional (1D) stellar evolution models alters the structure of the star in a manner that has a significant impact on the character of multidimensional convection. Adding modest amounts of rotation to a stellar model that ignores rotation during the evolutionary stage, however, has little impact on the character of resulting convection. Since the spatial scale and strength of convection present at the point of core collapse directly influence the supernova mechanism, our results suggest that rotation could play an important role in setting the stage for massive stellar explosions.

Our understanding of stars and their fates is based on coupling observations to theoretical models. Unlike laboratory physicists, we cannot perform experiments on stars, but rather must patiently take what nature allows us to observe. Simulation offers a means of virtual experimentation, enabling a detailed understanding of the most violent ongoing explosions in the Universe---the deaths of stars.

Carbon-rich grains with isotopic anomalies compared to the Sun are found in primitive meteorites. They were made by stars, and carry the original stellar nucleosynthesis signature. Silicon carbide grains of Type X and C, and low-density graphites condensed in the ejecta of core-collapse supernovae. We present a new set of models for the explosive He shell and compare them with the grains showing 12C/13C and 14N/15N ratios lower than solar. In the stellar progenitor H was ingested into the He shell and not fully destroyed before the explosion. Different explosion energies and H concentrations are considered. If the SN shock hits the He-shell region with some H still present, the models can reproduce the C and N isotopic signatures in C-rich grains. Hot-CNO cycle isotopic signatures are obtained, including a large production of 13C and 15N. The short-lived radionuclides 22Na and 26Al are increased by orders of magnitude. The production of radiogenic 22Ne from the decay of 22Na in the He shell might solve the puzzle of the Ne-E(L) component in low-density graphite grains. This scenario is attractive for the SiC grains of type AB with 14N/15N ratios lower than solar, and provides an alternative solution for SiC grains originally classified as nova grains. Finally, this process may contribute to the production of 14N and 15N in the Galaxy, helping to produce the 14N/15N ratio in the solar system.

We present the first three dimensional (3D) simulation of the final minutes of iron core growth in a massive star, up to and including the point of core gravitational instability and collapse. We self-consistently capture the development of strong convection driven by violent Si burning in the shell surrounding the iron core. This convective burning builds the iron core to its critical (Chandrasekhar) mass and collapse ensues, driven by electron capture and photodisintegration. The non-spherical structure and motion (turbulent fluctuations) generated by 3D convection is substantial at the point of collapse. We examine the impact of such physically-realistic 3D initial conditions on the core-collapse supernova mechanism using 3D simulations including multispecies neutrino leakage. We conclude that non-spherical progenitor structure should not be ignored, and has a significant and favorable impact on the likelihood for neutrino-driven explosions.

We explore the detailed and broad properties of carbon burning in Super Asymptotic Giant Branch (SAGB) stars with 2755 MESA stellar evolution models. The location of first carbon ignition, quenching location of the carbon burning flames and flashes, angular frequency of the carbon core, and carbon core mass are studied as a function of the ZAMS mass, initial rotation rate, and mixing parameters such as convective overshoot, semiconvection, thermohaline and angular momentum transport. In general terms, we find these properties of carbon burning in SAGB models are not a strong function of the initial rotation profile, but are a sensitive function of the overshoot parameter. We quasi-analytically derive an approximate ignition density, $\rho_{ign} \approx 2.1 \times 10^6$ g cm$^{-3}$, to predict the location of first carbon ignition in models that ignite carbon off-center. We also find that overshoot moves the ZAMS mass boundaries where off-center carbon ignition occurs at a nearly uniform rate of $\Delta M_{\rm ZAMS}$/$\Delta f_{\rm{ov}}\approx$ 1.6 $M_{\odot}$. For zero overshoot, $f_{\rm{ov}}$=0.0, our models in the ZAMS mass range $\approx$ 8.9 to 11 $M_{\odot}$ show off-center carbon ignition. For canonical amounts of overshooting, $f_{\rm{ov}}$=0.016, the off-center carbon ignition range shifts to $\approx$ 7.2 to 8.8 $M_{\odot}$. Only systems with $f_{\rm{ov}}$ $\geq 0.01$ and ZAMS mass $\approx$ 7.2-8.0 $M_{\odot}$ show carbon burning is quenched a significant distance from the center. These results suggest a careful assessment of overshoot modeling approximations on claims that carbon burning quenches an appreciable distance from the center of the carbon core.

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.

Rigorous statistical methods for estimating thermonuclear reaction rates and nucleosynthesis are becoming increasingly established in nuclear astrophysics. The main challenge being faced is that experimental reaction rates are highly complex quantities derived from a multitude of different measured nuclear parameters (e.g., astrophysical S-factors, resonance energies and strengths, particle and gamma-ray partial widths). We discuss the application of the Monte Carlo method to two distinct, but related, questions. First, given a set of measured nuclear parameters, how can one best estimate the resulting thermonuclear reaction rates and associated uncertainties? Second, given a set of appropriate reaction rates, how can one best estimate the abundances from nucleosynthesis (i.e., reaction network) calculations? The techniques described here provide probability density functions that can be used to derive statistically meaningful reaction rates and final abundances for any desired coverage probability. Examples are given for applications to s-process neutron sources, core-collapse supernovae, classical novae, and big bang nucleosynthesis.

We revisit the evidence for the contribution of the long-lived radioactive nuclides 44Ti, 55Fe, 56Co, 57Co, and 60Co to the UVOIR light curve of SN 1987A. We show that the V-band luminosity constitutes a roughly constant fraction of the bolometric luminosity between 900 and 1900 days, and we obtain an approximate bolometric light curve out to 4334 days by scaling the late time V-band data by a constant factor where no bolometric light curve data is available. Considering the five most relevant decay chains starting at 44Ti, 55Co, 56Ni, 57Ni, and 60Co, we perform a least squares fit to the constructed composite bolometric light curve. For the nickel isotopes, we obtain best fit values of M(56Ni) = (7.1 +- 0.3) x 10^{-2} Msun and M(57Ni) = (4.1 +- 1.8) x 10^{-3} Msun. Our best fit 44Ti mass is M(44Ti) = (0.55 +- 0.17) x 10^{-4} Msun, which is in disagreement with the much higher (3.1 +- 0.8) x 10^{-4} Msun recently derived from INTEGRAL observations. The associated uncertainties far exceed the best fit values for 55Co and 60Co and, as a result, we only give upper limits on the production masses of M(55Co) < 7.2 x 10^{-3} Msun and M(60Co) < 1.7 x 10^{-4} Msun. Furthermore, we find that the leptonic channels in the decay of 57Co (internal conversion and Auger electrons) are a significant contribution and constitute up to 15.5% of the total luminosity. Consideration of the kinetic energy of these electrons is essential in lowering our best fit nickel isotope production ratio to [57Ni/56Ni]=2.5+-1.1, which is still somewhat high but is in agreement with gamma-ray observations and model predictions.

We compile spectroscopic abundance data from 84 literature sources for 50 elements across 3058 stars in the solar neighborhood, within 150 pc of the Sun, to produce the Hypatia Catalog. We evaluate the variability of the spread in abundance measurements reported for the same star by different surveys. We also explore the likely association of the star within the Galactic disk, the corresponding observation and abundance determination methods for all catalogs in Hypatia, the influence of specific catalogs on the overall abundance trends, and the effect of normalizing all abundances to the same solar scale. The resulting large number of stellar abundance determinations in the Hypatia Catalog are analyzed only for thin-disk stars with observations that are consistent between literature sources. As a result of our large dataset, we find that the stars in the solar neighborhood may be reveal an asymmetric abundance distribution, such that a [Fe/H]-rich group near to the mid-plane is deficient in Mg, Si, S, Ca, Sc II, Cr II, and Ni as compared to stars further from the plane. The Hypatia Catalog has a wide number of applications, including exoplanet hosts, thick and thin disk stars, or stars with different kinematic properties.

There is compelling evidence that the peak brightness of a Type Ia supernova is affected by the electron fraction Ye at the time of the explosion. The electron fraction is set by the aboriginal composition of the white dwarf and the reactions that occur during the pre explosive convective burning. To date, determining the makeup of the white dwarf progenitor has relied on indirect proxies, such as the average metallicity of the host stellar population. In this paper, we present analytical calculations supporting the idea that the electron fraction of the progenitor systematically influences the nucleosynthesis of silicon group ejecta in Type Ia supernovae. In particular, we suggest the abundances generated in quasi nuclear statistical equilibrium are preserved during the subsequent freezeout. This allows one to potential recovery of Ye at explosion from the abundances recovered from an observed spectra. We show that measurement of 28Si, 32S, 40Ca, and 54Fe abundances can be used to construct Ye in the silicon rich regions of the supernovae. If these four abundances are determined exactly, they are sufficient to recover Ye to 6 percent. This is because these isotopes dominate the composition of silicon-rich material and iron rich material in quasi nuclear statistical equilibrium. Analytical analysis shows that the 28Si abundance is insensitive to Ye, the 32S abundance has a nearly linear trend with Ye, and the 40Ca abundance has a nearly quadratic trend with Ye. We verify these trends with post-processing of 1D models and show that these trends are reflected in model synthetic spectra.

Carbon-rich grains are observed to condense in the ejecta of recent core-collapse supernovae, within a year after the explosion. Silicon carbide grains of type X are C-rich grains with isotpic signatures of explosive supernova nucleosynthesis have been found in primitive meteorites. Much rarer silicon carbide grains of type C are a special sub-group of SiC grains from supernovae. They show peculiar abundance signatures for Si and S, isotopically heavy Si and isotopically light S, which appear to to be in disagreement with model predictions. We propose that C grains are formed mostly from C-rich stellar material exposed to lower SN shock temperatures than the more common type X grains. In this scenario, extreme $^{32}$S enrichments observed in C grains may be explained by the presence of short-lived $^{32}$Si ($\tau$$_{1/2}$ = 153 years) in the ejecta, produced by neutron capture processes starting from the stable Si isotopes. No mixing from deeper Si-rich material and/or fractionation of Si from S due to molecular chemistry is needed to explain the $^{32}$S enrichments. The abundance of $^{32}$Si in the grains can provide constraints on the neutron density reached during the supernova explosion in the C-rich He shell material. The impact of the large uncertainty of the neutron capture cross sections in the $^{32}$Si region is discussed.

We substantially update the capabilities of the open source software package Modules for Experiments in Stellar Astrophysics (MESA), and its one-dimensional stellar evolution module, MESA Star. Improvements in MESA Star's ability to model the evolution of giant planets now extends its applicability down to masses as low as one-tenth that of Jupiter. The dramatic improvement in asteroseismology enabled by the space-based Kepler and CoRoT missions motivates our full coupling of the ADIPLS adiabatic pulsation code with MESA Star. This also motivates a numerical recasting of the Ledoux criterion that is more easily implemented when many nuclei are present at non-negligible abundances. This impacts the way in which MESA Star calculates semi-convective and thermohaline mixing. We exhibit the evolution of 3-8 Msun stars through the end of core He burning, the onset of He thermal pulses, and arrival on the white dwarf cooling sequence. We implement diffusion of angular momentum and chemical abundances that enable calculations of rotating-star models, which we compare thoroughly with earlier work. We introduce a new treatment of radiation-dominated envelopes that allows the uninterrupted evolution of massive stars to core collapse. This enables the generation of new sets of supernovae, long gamma-ray burst, and pair-instability progenitor models. We substantially modify the way in which MESA Star solves the fully coupled stellar structure and composition equations, and we show how this has improved MESA's performance scaling on multi-core processors. Updates to the modules for equation of state, opacity, nuclear reaction rates, and atmospheric boundary conditions are also provided. We describe the MESA Software Development Kit (SDK) that packages all the required components needed to form a unified and maintained build environment for MESA. [Abridged]

The stellar mass range 8<M/Mo<12 corresponds to the most massive AGB stars and the most numerous massive stars. It is host to a variety of supernova progenitors and is therefore very important for galactic chemical evolution and stellar population studies. In this paper, we study the transition from super-AGB star to massive star and find that a propagating neon-oxygen burning shell is common to both the most massive electron capture supernova (EC-SN) progenitors and the lowest mass iron-core collapse supernova (FeCCSN) progenitors. Of the models that ignite neon burning off-center, the 9.5Mo model would evolve to an FeCCSN after the neon-burning shell propagates to the center, as in previous studies. The neon-burning shell in the 8.8Mo model, however, fails to reach the center as the URCA process and an extended (0.6 Mo) region of low Ye (0.48) in the outer part of the core begin to dominate the late evolution; the model evolves to an EC-SN. This is the first study to follow the most massive EC-SN progenitors to collapse, representing an evolutionary path to EC-SN in addition to that from SAGB stars undergoing thermal pulses. We also present models of an 8.75Mo super-AGB star through its entire thermal pulse phase until electron captures on 20Ne begin at its center and of a 12Mo star up to the iron core collapse. We discuss key uncertainties and how the different pathways to collapse affect the pre-supernova structure. Finally, we compare our results to the observed neutron star mass distribution.

STARLIB is a next-generation, all-purpose nuclear reaction-rate library. For the first time, this library provides the rate probability density at all temperature grid points for convenient implementation in models of stellar phenomena. The recommended rate and its associated uncertainties are also included. Currently, uncertainties are absent from all other rate libraries, and, although estimates have been attempted in previous evaluations and compilations, these are generally not based on rigorous statistical definitions. A common standard for deriving uncertainties is clearly warranted. STARLIB represents a first step in addressing this deficiency by providing a tabular, up-to-date database that supplies not only the rate and its uncertainty but also its distribution. Because a majority of rates are lognormally distributed, this allows the construction of rate probability densities from the columns of STARLIB. This structure is based on a recently suggested Monte Carlo method to calculate reaction rates, where uncertainties are rigorously defined. In STARLIB, experimental rates are supplemented with: (i) theoretical TALYS rates for reactions for which no experimental input is available, and (ii) laboratory and theoretical weak rates. STARLIB includes all types of reactions of astrophysical interest to Z = 83, such as (p,g), (p,a), (a,n), and corresponding reverse rates. Strong rates account for thermal target excitations. Here, we summarize our Monte Carlo formalism, introduce the library, compare methods of correcting rates for stellar environments, and discuss how to implement our library in Monte Carlo nucleosynthesis studies. We also present a method for accessing STARLIB on the Internet and outline updated Monte Carlo-based rates.

Novae are cataclysmic variable binary systems in which a white dwarf primary is accreting material from a low mass companion. The importance of this accretion takes on added significance if the WD can increase its mass to reach the Chandrasekhar limit thus exploding as a Type Ia supernova. In this study we accrete material of Solar composition onto carbon-oxygen white dwarfs of 0.70, 1.00 and 1.35 Msun with accretion rates from 1.6e-10 to 1.6e-6 Msun per yr. We have utilized the MESA stellar evolution code for our modeling and evolve them for many nova cycles or, in some cases, evolution to a red giant stage. Differing behaviors occur as a function of both the WD mass and the accretion rate. For the lower WD masses, the models undergo recurrent hydrogen flashes at low accretion rates; for higher accretion rates, steady-burning of hydrogen occurs and eventually gives way to recurrent hydrogen flashes. At the highest accretion rates, these models go through a steady-burning phase but eventually transition into red giants. For the highest white dwarf mass recurrent hydrogen flashes occur at lower accretion rates but for higher rates the models exhibit steady-burning interspersed with helium flashes. We find that for all our models that undergo recurrent hydrogen flashes, as well as the steady-burning models that exhibit helium flashes, the mass of the WD continues to grow toward the Chandrasekhar limit. These results suggest that the accretion of Solar abundance material onto carbon-oxygen white dwarfs in cataclysmic variable systems, the single degenerate scenario, is a viable channel for progenitors of Type Ia supernova explosions.

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.

Despite much effort in the past decades, the C-burning reaction rate is uncertain by several orders of magnitude, and the relative strength between the different channels 12C(12C,alpha)20Ne, 12C(12C,p)23Na and 12C(12C,n)23Mg is poorly determined. Additionally, in C-burning conditions a high 12C+12C rate may lead to lower central C-burning temperatures and to 13C(alpha,n)16O emerging as a more dominant neutron source than 22Ne(alpha,n)25Mg, increasing significantly the s-process production. This is due to the rapid decrease of the 13N(gamma,p)12C with decreasing temperature, causing the 13C production via 13N(beta+)13C. Presented here is the impact of the 12C+12C reaction uncertainties on the s-process and on explosive p-process nucleosynthesis in massive stars, including also fast rotating massive stars at low metallicity. Using various 12C+12C rates, in particular an upper and lower rate limit of ~ 50000 higher and ~ 20 lower than the standard rate at 5*10^8 K, five 25 Msun stellar models are calculated. The enhanced s-process signature due to 13C(alpha,n)16O activation is considered, taking into account the impact of the uncertainty of all three C-burning reaction branches. Consequently, we show that the p-process abundances have an average production factor increased up to about a factor of 8 compared to the standard case, efficiently producing the elusive Mo and Ru proton-rich isotopes. We also show that an s-process being driven by 13C(alpha,n)16O is a secondary process, even though the abundance of 13C does not depend on the initial metal content. Finally, implications for the Sr-peak elements inventory in the Solar System and at low metallicity are discussed.

Symbiotic binaries are systems containing white dwarfs (WDs) and red giants. Symbiotic novae are those systems in which thermonuclear eruptions occur on the WD components. These are to be distinguished from events driven by accretion disk instabilities analogous to dwarf novae eruptions in cataclysmic variable outbursts. Another class of symbiotic systems are those in which the WD is extremely luminous and it seems likely that quiescent nuclear burning is ongoing on the accreting WD. A fundamental question is the secular evolution of the WD. Do the repeated outbursts or quiescent burning in these accreting systems cause the WD to gain or lose mass? If it is gaining mass, can it eventually reach the Chandrasekhar Limit and become a supernova (a SN Ia if it can hide the hydrogen and helium in the system)? In order to better understand these systems, we have begun a new study of the evolution of Thermonuclear Runaways (TNRs) in the accreted envelopes of WDs using a variety of initial WD masses, luminosities and mass accretion rates. We use our 1-D hydro code, NOVA, which includes the new convective algorithm of Arnett, Meakin and Young, the Hix and Thielemann nuclear reaction solver, the Iliadis reaction rate library, the Timmes equation of state, and the OPAL opacities. We assume a solar composition (Lodders abundance distribution) and do not allow any mixing of accreted material with core material. This assumption strongly influences our results. We report here (1) that the WD grows in mass for all simulations so that canonical `steady burning' does not occur, and (2) that only a small fraction of the accreted matter is ejected in some (but not all) simulations. We also find that the accreting systems, before thermonuclear runaway, are too cool to be seen in X-ray searches for SN Ia progenitors.

We present a brief summary of the Single Degenerate Scenario for the progenitors of Type Ia Supernovae in which it is assumed that a low mass carbon-oxygen white dwarf is growing in mass as a result of accretion from a secondary star in a close binary system. Recent hydrodynamic simulations of accretion of solar material onto white dwarfs without mixing always produce a thermonuclear runaway and steady burning does not occur. For a broad range in WD mass (0.4 Solar masses to 1.35 Solar Masses), the maximum ejected material occurs for the 1.25 Solar Mass sequences and then decreases as the white dwarf mass decreases. Therefore, the white dwarfs are growing in mass as a consequence of the accretion of solar material and as long as there is no mixing of accreted material with core material. In contrast, a thermonuclear runaway in the accreted hydrogen-rich layers on the low luminosity WDs in close binary systems where mixing of core matter with accreted material has occurred is the outburst mechanism for Classical, Recurrent, and Symbiotic novae. The differences in characteristics of these systems is likely the WD mass and mass accretion rate. The high levels of enrichment of CN ejecta in elements ranging from carbon to sulfur confirm that there is dredge-up of matter from the core of the WD and enable them to contribute to the chemical enrichment of the interstellar medium. Therefore, studies of CNe can lead to an improved understanding of Galactic nucleosynthesis, some sources of pre-solar grains, and the Extragalactic distance scale. The characteristics of the outburst depend on the white dwarf mass, luminosity, mass accretion rate, and the chemical composition of both the accreting material and WD material. The properties of the outburst also depends on when, how, and if the accreted layers are mixed with the WD core and the mixing mechanism is still unknown.

We systematically explore zero impact parameter collisions of white dwarfs with the Eulerian adaptive grid code FLASH for 0.64+0.64 M$_{\odot}$ and 0.81+0.81 M$_{\odot}$ mass pairings. Our models span a range of effective linear spatial resolutions from 5.2$\times10^{7}$ to 1.2$\times10^{7}$ cm. However, even the highest resolution models do not quite achieve strict numerical convergence, due to the challenge of properly resolving small-scale burning and energy transport. The lack of strict numerical convergence from these idealized configurations suggest that quantitative predictions of the ejected elemental abundances that are generated by binary white dwarf collision and merger simulations should be viewed with caution. Nevertheless, the convergence trends do allow some patterns to be discerned. We find that the 0.64+0.64 M$_{\odot}$ head-on collision model produces 0.32 M$_{\odot}$ of \nickel[56] and 0.38 M$_{\odot}$ of \silicon[28], while the 0.81+0.81 M$_{\odot}$ head-on collision model produces 0.39 M$_{\odot}$ of \nickel[56] and 0.55 M$_{\odot}$ of \silicon[28] at the highest spatial resolutions. Both mass pairings produce $\sim$0.2 M$_{\odot}$ of unburned \carbon[12]+\oxygen[16]. We also find the 0.64+0.64 M$_{\odot}$ head-on collision begins carbon burning in the central region of the stalled shock between the two white dwarfs, while the more energetic \hbox{0.81+0.81 M$_{\odot}$} head-on collision raises the initial post-shock temperature enough to burn the entire stalled shock region to nuclear statistical equilibrium.

We present a study exploring a systematic effect on the brightness of type Ia supernovae using numerical models that assume the single-degenerate paradigm. Our investigation varied the central density of the progenitor white dwarf at flame ignition, and considered its impact on the explosion yield, particularly the production and distribution of radioactive Ni-56, which powers the light curve. We performed a suite of two-dimensional simulations with randomized initial conditions, allowing us to characterize the statistical trends that we present. The simulations indicate that production of Fe-group material is statistically independent of progenitor central density, but the mass of stable Fe-group isotopes is tightly correlated with central density, with a decrease in the production of Ni-56 at higher central densities. These results imply progenitors with higher central densities produce dimmer events. We provide details of the post-explosion distribution of Ni-56 in the models, including the lack of a consistent centrally-located deficit of Ni-56, which may be compared to observed remnants. By performing a self-consistent extrapolation of our model yields and considering the main-sequence lifetime of the progenitor star and the elapsed time between the formation of the white dwarf and the onset of accretion, we develop a brightness-age relation that improves our prediction of the expected trend for single degenerates and we compare this relation with observations.

Several lines of evidence, from isotopic analyses of meteorites to studies of the Sun's elemental and isotopic composition, indicate that the solar system was contaminated early in its evolution by ejecta from a nearby supernova (SN). Previous models have invoked SN material being injected into an extant protoplanetary disk, or isotropically expanding ejecta sweeping over a distant (>10 pc) cloud core, simultaneously enriching it and triggering its collapse. Here we consider a new astrophysical setting: the injection of clumpy SN ejecta, as observed in the Cas A SN remnant, into the molecular gas at the periphery of an HII region created by the SN's progenitor star. To track these interactions we have conducted a suite of high-resolution (1500^3 effective) 3D simulations that follow the evolution of individual clumps as they move into molecular gas. Even at these high resolutions, our simulations do not quite achieve numerical convergence, due to the challenge of properly resolving the small-scale mixing of ejecta and molecular gas, although they do allow some robust conclusions to be drawn. Isotropically exploding ejecta do not penetrate into the molecular cloud, but, if cooling is properly accounted for, clumpy ejecta penetrate to distances ~10^18 cm and mix effectively with star-forming molecular gas. The ~2 M_\odot high-metallicity ejecta from a core-collapse SN is likely to mix with ~2 \times 10^4 M_\odot of molecular gas. Thus all stars forming late (~5 Myr) in the evolution of an HII region may be contaminated by SN ejecta at a level ~10^-4. This level of contamination is consistent with the abundances of short-lived radionuclides and possibly some stable isotopic shifts in the early solar system, and is potentially consistent with the observed variability in stellar elemental abundances. SN contamination of forming planetary systems may be a common, universal process.

Type Ia supernovae are bright stellar explosions distinguished by standardizable light curves that allow for their use as distance indicators for cosmological studies. Despite their highly successful use in this capacity, the progenitors of these events are incompletely understood. We describe simulating type Ia supernovae in the paradigm of a thermonuclear runaway occurring in a massive white dwarf star. We describe the multi-scale physical processes that realistic models must incorporate and the numerical models for these that we employ. In particular, we describe a flame-capturing scheme that addresses the problem of turbulent thermonuclear combustion on unresolved scales. We present the results of our study of the systematics of type Ia supernovae including trends in brightness following from properties of the host galaxy that agree with observations. We also present performance results from simulations on leadership-class architectures.

[Shortened] The 12C + 12C fusion reaction has been the subject of considerable experimental efforts to constrain uncertainties at temperatures relevant for stellar nucleosynthesis. In order to investigate the effect of an enhanced carbon burning rate on massive star structure and nucleosynthesis, new stellar evolution models and their yields are presented exploring the impact of three different 12C + 12C reaction rates. Non-rotating stellar models were generated using the Geneva Stellar Evolution Code and were later post-processed with the NuGrid Multi-zone Post-Processing Network tool. The enhanced rate causes core carbon burning to be ignited more promptly and at lower temperature. This reduces the neutrino losses, which increases the core carbon burning lifetime. An increased carbon burning rate also increases the upper initial mass limit for which a star exhibits a convective carbon core. Carbon shell burning is also affected, with fewer convective-shell episodes and convection zones that tend to be larger in mass. Consequently, the chance of an overlap between the ashes of carbon core burning and the following carbon shell convection zones is increased, which can cause a portion of the ashes of carbon core burning to be included in the carbon shell. Therefore, during the supernova explosion, the ejecta will be enriched by s-process nuclides synthesized from the carbon core s process. The yields were used to estimate the weak s-process component in order to compare with the solar system abundance distribution. The enhanced rate models were found to produce a significant proportion of Kr, Sr, Y, Zr, Mo, Ru, Pd and Cd in the weak component, which is primarily the signature of the carbon-core s process. Consequently, it is shown that the production of isotopes in the Kr-Sr region can be used to constrain the 12C + 12C rate using the current branching ratio for a- and p-exit channels.

We carry out a comprehensive smooth particle hydrodynamics simulation survey of double-degenerate white dwarf binary mergers of varying mass combinations in order to establish correspondence between initial conditions and remnant configurations. We find that all but one of our simulation remnants share general properties such as a cold, degenerate core surrounded by a hot disk, while our least massive pair of stars forms only a hot disk. We characterize our remnant configurations by the core mass, the rotational velocity of the core, and the half-mass radius of the disk. We also find that some of our simulations with very massive constituent stars exhibit helium detonations on the surface of the primary star before complete disruption of the secondary. However, these helium detonations are insufficiently energetic to ignite carbon, and so do not lead to prompt carbon detonations.

We explore the nucleosynthesis trends from two mechanisms during freeze-out expansions in core-collapse supernovae. The first mechanism is related to the convection and instabilities within homogeneous stellar progenitor matter that is accreted through the supernova shock. The second mechanism is related to the impact of the supersonic wind termination shock (reverse shock) within the tumultuous inner regions of the ejecta above the proto-neutron star. Our results suggest that isotopes in the mass range 12<=A<=122 that are produced during the freeze-out expansions may be classified in two families. The isotopes of the first family manifest a common mass fraction evolutionary profile, whose specific shape per isotope depends on the characteristic transition between two equilibrium states (equilibrium state transition) during each type of freeze-out expansion. The first family includes the majority of isotopes in this mass range. The second family is limited to magic nuclei and isotopes in their locality, which do not sustain any transition, become nuclear flow hubs, and dominate the final composition. We use exponential and power-law adiabatic profiles to identify dynamic large-scale and small-scale equilibrium patterns among nuclear reactions. (truncated abstract) In addition, we introduce non-monotonic parameterized profiles to probe the impact of the reverse shock and multi-dimensional explosion asymmetries on nucleosynthesis. Cases are shown in which the non-monotonic profiles favor the production of radioactivities. Non-monotonic freeze-out profiles involve longer non-equilibrium nucleosynthesis intervals compared with the exponential and power-law profiles, resulting in mass fraction trends and yield distributions that may not be achieved by the monotonic freeze-out profiles.

The contribution by massive stars (M > 9 solar masses) to the weak s-process component of the solar system abundances is primarily due to the 22Ne neutron source, which is activated near the end of helium-core burning. The residual 22Ne left over from helium-core burning is then reignited during carbon burning, initiating further s-processing that modifies the isotopic distribution. This modification is sensitive to the stellar structure and the carbon burning reaction rate. Recent work on the 12C + 12C reaction suggests that resonances located within the Gamow peak may exist, causing a strong increase in the astrophysical S-factor and consequently the reaction rate. To investigate the effect of an increased rate, 25 solar mass stellar models with three different carbon burning rates, at solar metallicity, were generated using the Geneva Stellar Evolution Code (GENEC) with nucleosynthesis post-processing calculated using the NuGrid Multi-zone Post-Processing Network code (MPPNP). The strongest rate caused carbon burning to occur in a large convective core rather than a radiative one. The presence of this large convective core leads to an overlap with the subsequent convective carbon-shell, significantly altering the initial composition of the carbon-shell. In addition, an enhanced rate causes carbon-shell burning episodes to ignite earlier in the evolution of the star, igniting the 22Ne source at lower temperatures and reducing the neutron density.

Type Ia supernovae are bright stellar explosions thought to occur when a thermonuclear runaway consumes roughly a solar mass of degenerate stellar material. These events produce and disseminate iron-peak elements, and properties of their light curves allow for standardization and subsequent use as cosmological distance indicators. The explosion mechanism of these events remains, however, only partially understood. Many models posit the explosion beginning with a deflagration born near the center of a white dwarf that has gained mass from a stellar companion. In order to match observations, models of this single-degenerate scenario typically invoke a subsequent transition of the (subsonic) deflagration to a (supersonic) detonation that rapidly consumes the star. We present an investigation into the systematics of thermonuclear supernovae assuming this paradigm. We utilize a statistical framework for a controlled study of two-dimensional simulations of these events from randomized initial conditions. We investigate the effect of the composition and thermal history of the progenitor on the radioactive yield, and thus brightness, of an event. Our results offer an explanation for some observed trends of mean brightness with properties of the host galaxy.

We compare the yields of Ti44 and Ni56 produced from post-processing the thermodynamic trajectories from three different core-collapse models -- a Cassiopeia A progenitor, a double shock hypernova progenitor, and a rotating 2D explosion -- with the yields from exponential and power-law trajectories. The peak temperatures and densities achieved in these core-collapse models span several of the distinct nucleosynthesis regions we identify, resulting in different trends in the Ti44 and Ni56 yields for different mass elements. The Ti44 and Ni56 mass fraction profiles from the exponential and power-law profiles generally explain the tendencies of the post-processed yields, depending on which regions are traversed by the model. We find integrated yields of Ti44 and Ni56 from the exponential and power-law trajectories are generally within a factor 2 or less of the post-process yields. We also analyze the influence of specific nuclear reactions on the Ti44 and Ni56 abundance evolution. Reactions that affect all yields globally are the 3a, p(e-,nu)n and n(e+,nubar)p. The rest of the reactions are ranked according to their degree of impact on the synthesis of Ti44. The primary ones include Ti44(a,p)V47, Ca40(a,g)Ti44, V45(p,g)Cr46, Ca40(a,p)Sc43, F17(a,p)Ne20, Na21(a,p)Mg24, Sc41(p,g)Ti42, Sc43(p,g)Ti44, Ti44(p,g)V45, and Ni57(p,g)Cu58, along with numerous weak reactions. Our analysis suggests that not all Ti44 need be produced in an a-rich freeze-out in core-collapse events, and that reaction rate equilibria in combination with timescale effects for the expansion profile may account for the paucity of Ti44 observed in supernovae remnants.

We present a comprehensive study of white dwarf collisions as an avenue for creating type Ia supernovae. Using a smooth particle hydrodynamics code with a 13-isotope, {\alpha}-chain nuclear network, we examine the resulting 56Ni yield as a function of total mass, mass ratio, and impact parameter. We show that several combinations of white dwarf masses and impact parameters are able to produce sufficient quantities of 56Ni to be observable at cosmological distances. We find the 56Ni production in double-degenerate white dwarf collisions ranges from sub-luminous to the super-luminous, depending on the parameters of the collision. For all mass pairs, collisions with small impact parameters have the highest likelihood of detonating, but 56Ni production is insensitive to this parameter in high-mass combinations, which significantly increases their likelihood of detection. We also find that the 56Ni dependence on total mass and mass ratio is not linear, with larger mass primaries producing disproportionately more 56Ni than their lower mass secondary counterparts, and symmetric pairs of masses producing more 56Ni than asymmetric pairs.

Stellar physics and evolution calculations enable a broad range of research in astrophysics. Modules for Experiments in Stellar Astrophysics (MESA) is a suite of open source libraries for a wide range of applications in computational stellar astrophysics. A newly designed 1-D stellar evolution module, MESA star, combines many of the numerical and physics modules for simulations of a wide range of stellar evolution scenarios ranging from very-low mass to massive stars, including advanced evolutionary phases. MESA star solves the fully coupled structure and composition equations simultaneously. It uses adaptive mesh refinement and sophisticated timestep controls, and supports shared memory parallelism based on OpenMP. Independently usable modules provide equation of state, opacity, nuclear reaction rates, and atmosphere boundary conditions. Each module is constructed as a separate Fortran 95 library with its own public interface. Examples include comparisons to other codes and show evolutionary tracks of very low mass stars, brown dwarfs, and gas giant planets; the complete evolution of a 1 Msun star from the pre-main sequence to a cooling white dwarf; the Solar sound speed profile; the evolution of intermediate mass stars through the thermal pulses on the He-shell burning AGB phase; the interior structure of slowly pulsating B Stars and Beta Cepheids; evolutionary tracks of massive stars from the pre-main sequence to the onset of core collapse; stars undergoing Roche lobe overflow; and accretion onto a neutron star. Instructions for downloading and installing MESA can be found on the project web site (http://mesa.sourceforge.net/).

We explore the effects of the deflagration to detonation transition (DDT) density on the production of Ni-56 in thermonuclear supernova explosions (type Ia supernovae). Within the DDT paradigm, the transition density sets the amount of expansion during the deflagration phase of the explosion and therefore the amount of nuclear statistical equilibrium (NSE) material produced. We employ a theoretical framework for a well-controlled statistical study of two-dimensional simulations of thermonuclear supernovae with randomized initial conditions that can, with a particular choice of transition density, produce a similar average and range of Ni-56 masses to those inferred from observations. Within this framework, we utilize a more realistic "simmered" white dwarf progenitor model with a flame model and energetics scheme to calculate the amount of Ni-56 and NSE material synthesized for a suite of simulated explosions in which the transition density is varied in the range 1-3x10^7 g/cc. We find a quadratic dependence of the NSE yield on the log of the transition density, which is determined by the competition between plume rise and stellar expansion. By considering the effect of metallicity on the transition density, we find the NSE yield decreases by 0.055 +/- 0.004 solar masses for a 1 solar metallicity increase evaluated about solar metallicity. For the same change in metallicity, this result translates to a 0.067 +/- 0.004 solar mass decrease in the Ni-56 yield, slightly stronger than that due to the variation in electron fraction from the initial composition. Observations testing the dependence of the yield on metallicity remain somewhat ambiguous, but the dependence we find is comparable to that inferred from some studies.

Recent observational studies of type Ia supernovae (SNeIa) suggest correlations between the peak brightness of an event and the age of the progenitor stellar population. This trend likely follows from properties of the progenitor white dwarf (WD), such as central density, that follow from properties of the host stellar population. We present a statistically well-controlled, systematic study utilizing a suite of multi-dimensional SNeIa simulations investigating the influence of central density of the progenitor WD on the production of Fe-group material, particularly radioactive Ni-56, which powers the light curve. We find that on average, as the progenitor's central density increases, production of Fe-group material does not change but production of Ni-56 decreases. We attribute this result to a higher rate of neutronization at higher density. The central density of the progenitor is determined by the mass of the WD and the cooling time prior to the onset of mass transfer from the companion, as well as the subsequent accretion heating and neutrino losses. The dependence of this density on cooling time, combined with the result of our central density study, offers an explanation for the observed age-luminosity correlation: a longer cooling time raises the central density at ignition thereby producing less Ni-56 and thus a dimmer event. While our ensemble of results demonstrates a significant trend, we find considerable variation between realizations, indicating the necessity for averaging over an ensemble of simulations to demonstrate a statistically significant result.

The merger of two white dwarfs (a.k.a. double degenerate merger) has often been cited as a potential progenitor of type Ia supernovae. Here we combine population synthesis, merger and explosion models with radiation-hydrodynamics light-curve models to study the implications of such a progenitor scenario on the observed type Ia supernova population. Our standard model, assuming double degenerate mergers do produce thermonuclear explosions, produces supernova light-curves that are broader than the observed type Ia sample. In addition, we discuss how the shock breakout and spectral features of these double degenerate progenitors will differ from the canonical bare Chandrasekhar-massed explosion models. We conclude with a discussion of how one might reconcile these differences with current observations.

The slow neutron capture process in massive stars (the weak s-process) produces most of the s-only isotopes in the mass region 60 < A < 90. The nuclear reaction rates used in simulations of this process have a profound effect on the final s-process yields. We generated 1D stellar models of a 25 solar mass star varying the 12C + 12C rate by a factor of 10 and calculated full nucleosynthesis using the post-processing code PPN. Increasing or decreasing the rate by a factor of 10 affects the convective history and nucleosynthesis, and consequently the final yields.

We explore collisions between two white dwarfs as a pathway for making Type Ia Supernovae (SNIa). White dwarf number densities in globular clusters allow 10-100 redshift <1 collisions per year, and observations by (Chomiuk et al. 2008) of globular clusters in the nearby S0 galaxy NGC 7457 have detected what is likely to be a SNIa remnant. We carry out simulations of the collision between two 0.6 solar mass white dwarfs at various impact parameters and mass resolutions. For impact parameters less than half the radius of the white dwarf, we find such collisions produce approximately 0.4 solar masses of Ni56, making such events potential candidates for underluminous SNIa or a new class of transients between Novae and SNIa.

We present a theoretical framework for formal study of systematic effects in Supernovae Type Ia (SN Ia) that utilizes 2-d simulations to implement a form of the deflagration-detonation transition (DDT) explosion scenario. The framework is developed from a randomized initial condition that leads to a sample of simulated SN Ia whose Ni56 masses have a similar average and range to those observed, and have many other modestly realistic features such as the velocity extent of intermediate mass elements. The intended purpose is to enable statistically well-defined studies of both physical and theoretical parameters of the SN Ia explosion simulation. We present here a thorough description of the outcome of the SN Ia explosions produced by our current simulations. A first application of this framework is utilized to study the dependence of the SN Ia on the Ne22 content, which is known to be directly influenced by the progenitor stellar population's metallicity. Our study is very specifically tailored to measure how the Ne22 content influences the competition between the rise of plumes of burned material and the expansion of the star before these plumes reach DDT conditions. This competition controls the amount of material in nuclear statistical equilibrium (NSE) and therefore Ni56 produced by setting the density at which nucleosynthesis takes place during the detonation phase of the explosion. Although the outcome following from any particular ignition condition can change dramatically with Ne22 content, with a sample of 20 ignition conditions we find that the systematic change in the expansion of the star prior to detonation is not large enough to compete with the dependence on initial neutron excess discussed by Timmes, Brown & Truran (2003). (Abridged)

This white paper, directed to the Stars and Stellar Evolution panel, has three objectives: 1) to provide the Astro2010 Decadal Survey with a vista into the goals of the nuclear physics and nuclear astrophysics community; 2) to alert the astronomical community of joint opportunities for discoveries at the interface between nuclear physics and astronomy; and 3) to delineate efforts in nuclear physics and describe the observational and theoretical advances in astrophysics necessary to make progress towards answering the following questions in the Nuclear Science 2007 Long Range Plan: 1) What is the origin and distribution of the elements? 2) What are the nuclear reactions that power stars and stellar explosions? 3) What is the nature of dense matter? The scope of this white paper concerns the specific area of "low energy" nuclear astrophysics. We define this as the area of overlap between astrophysics and the study of nuclear structure and reactions. Of the questions listed above, two -- What is the origin of the elements? and What is the nature of dense matter? -- were specifically listed in the National Academies Study, "Connecting Quarks with the Cosmos".

The s-process production in massive stars at very low metallicities is expected to be negligible due to the low abundance of the neutron source 22Ne, to primary neutron poisons and decreasing iron seed abundances. However, recent models of massive stars including the effects of rotation show that a strong production of 22Ne is possible in the helium core, as a consequence of the primary nitrogen production (observed in halo metal poor stars). Using the PPN post-processing code, we studied the impact of this primary 22Ne on the s process. We find a large production of s elements between strontium and barium, starting with the amount of primary 22Ne predicted by stellar models. There are several key reaction rate uncertainties influencing the s-process efficiency. Among them, 17O(alpha,gamma) may play a crucial role strongly influencing the s process efficiency, or it may play a negligible role, according to the rate used in the calculations. We also report on the development of a new parallel (MPI) post-processing code (MPPNP) designed to follow the complete nucleosynthesis in stars on highly resolved grids. We present here the first post-processing run from the ZAMS up to the end of helium burning for a 15 solar mass model.

We present preliminary results from recent high-resolution double-degenerate merger simulations with the Smooth Particle Hydrodynamics (SPH) technique. We put particular emphasis on verification and validation in our effort and show the importance of details in the initial condition setup for the final outcome of the simulation. We also stress the dynamical importance of including shocks in the simulations. These results represent a first step toward a suite of simulations that will shed light on the question whether double-degenerate mergers are a viable path toward type 1a supernovae. In future simulations, we will make use of the capabilities of the NuGrid collaboration in post-processing SPH particle trajectories with a complete nuclear network to follow the detailed nuclear reactions during the dynamic merger phase.

We examine observed heavy element abundances in the Cassiopeia A supernova remnant as a constraint on the nature of the Cas A supernova. We compare bulk abundances from 1D and 3D explosion models and spatial distribution of elements in 3D models with those derived from X-ray observations. We also examine the cospatial production of 26Al with other species. We find that the most reliable indicator of the presence of 26Al in unmixed ejecta is a very low S/Si ratio (~0.05). Production of N in O/S/Si-rich regions is also indicative. The biologically important element P is produced at its highest abundance in the same regions. Proxies should be detectable in supernova ejecta with high spatial resolution multiwavelength observations.

Simulations of nucleosynthesis in astrophysical environments are at the intersection of nuclear physics reaction rate research and astrophysical applications, for example in the area of galactic chemical evolution or near-field cosmology. Unfortunately, at present the available yields for such applications are based on heterogeneous assumptions between the various contributing nuclear production sites, both in terms of modeling the thermodynamic environment itself as well as the choice of specifc nuclear reaction rates and compilations. On the other side, new nuclear reaction rate determinations are often taking a long time to be included in astrophysical applications. The NuGrid project addresses these issues by providing a set of codes and a framework in which these codes interact. In this contribution we describe the motivation, goals and first results of the NuGrid project. At the core is a new and evolving post-processing nuclesoynthesis code (PPN) that can follow quiescent and explosive nucleosynthesis following multi-zone 1D-stellar evolution as well as multi-zone hydrodynamic input, including explosions. First results are available in the areas of AGB and massive stars.

We investigate the physical conditions where 44Ti and 56Ni are created in core-collapse supernovae. In this preliminary work we use a series of post-processing network calculations with parametrized expansion profiles that are representative of the wide range of temperatures, densities and electron-to-baryon ratios found in 3D supernova simulations. Critical flows that affect the final yields of 44Ti and 56Ni are assessed.

The "collapsar" engine for gamma-ray bursts invokes as its energy source the failure of a normal supernova and the formation of a black hole. Here we present the results of the first three-dimensional simulation of the collapse of a massive star down to a black hole, including the subsequent accretion and explosion. The explosion differs significantly from the axisymmetric scenario obtained in two-dimensional simulations; this has important consequences for the nucleosynthetic yields. We compare the nucleosynthetic yields to those of hypernovae. Calculating yields from three-dimensional explosions requires new strategies in post-process nucleosynthesis; we discuss NuGrid's plan for three-dimensional yields.

We review some of the uncertainties in calculating nucleosynthetic yields, focusing on the explosion mechanism. Current yield calculations tend to either use a piston, energy injection, or enhancement of neutrino opacities to drive an explosion. We show that the energy injection, or more accurately, an entropy injection mechanism is best-suited to mimic our current understanding of the convection-enhanced supernova engine. The enhanced neutrino-opacity technique is in qualitative disagreement with simulations of core-collapse supernovae and will likely produce errors in the yields. But piston-driven explosions are the most discrepant. Piston-driven explosion severely underestimate the amount of fallback, leading to order-of-magnitude errors in the yields of heavy elements. To obtain yields accurate to the factor of a few level, we must use entropy or energy injection and this has become the NuGrid collaboration approach.

Many nucleosynthesis and mixing processes of low-mass stars as they evolve from the Main Sequence to the thermal-pulse Asymptotic Giant Branch phase (TP-AGB) are well understood (although of course important physics components, e.g. rotation, magnetic fields, gravity wave mixing, remain poorly known). Nevertheless, in the last years presolar grain measurements with high resolution have presented new puzzling problems and strong constraints on nucleosynthesis processes in stars. The goal of the NuGrid collaboration is to present uniform yields for a large range of masses and metallicities, including low$-$mass stars and massive stars and their explosions. Here we present the first calculations of stellar evolution and high-resolution, post-processing simulations of an AGB star with an initial mass of 2 M_sun and solar-like metallicity (Z=0.01), based on the post-processing code PPN. In particular, we analyze the formation and evolution of the radiative 13C-pocket between the 17th TP and the 18th TP. The s-process nucleosynthesis profile of a sample of heavy isotopes is also discussed, before the next convective TP occurrence.

The nucleosynthetic yield from a supernova explosion depends upon a variety of effects: progenitor evolution, explosion process, details of the nuclear network, and nuclear rates. Especially in studies of integrated stellar yields, simplifications reduce these uncertainties. But nature is much more complex, and to actually study nuclear rates, we will have to understand the full, complex set of processes involved in nucleosynthesis. Here we discuss a few of these complexities and detail how the NuGrid collaboration will address them.

We have continued our studies of the Classical Nova outburst by evolving TNRs on 1.25Msun and 1.35Msun WDs (ONeMg composition) under conditions which produce mass ejection and a rapid increase in the emitted light, by examining the effects of changes in the nuclear reaction rates on both the observable features and the nucleosynthesis during the outburst. In order to improve our calculations over previous work, we have incorporated a modern nuclear reaction network into our hydrodynamic computer code. We find that the updates in the nuclear reaction rate libraries change the amount of ejected mass, peak luminosity, and the resulting nucleosynthesis. In addition, as a result of our improvements, we discovered that the pep reaction was not included in our previous studies of CN explosions. Although the energy production from this reaction is not important in the Sun, the densities in WD envelopes can exceed $10^4$ gm cm$^{-3}$ and the presence of this reaction increases the energy generation during the time that the p-p chain is operating. The effect of the increased energy generation is to reduce the evolution time to the peak of the TNR and, thereby, the accreted mass as compared to the evolutionary sequences done without this reaction included. As expected from our previous work, the reduction in accreted mass has important consequences on the characteristics of the resulting TNR and is discussed in this paper.

Proton-rich material in a state of nuclear statistical equilibrium (NSE) is one of the least studied regimes of nucleosynthesis. One reason for this is that after hydrogen burning, stellar evolution proceeds at conditions of equal number of neutrons and protons or at a slight degree of neutron-richness. Proton-rich nucleosynthesis in stars tends to occur only when hydrogen-rich material that accretes onto a white dwarf or neutron star explodes, or when neutrino interactions in the winds from a nascent proto-neutron star or collapsar-disk drive the matter proton-rich prior to or during the nucleosynthesis. In this paper we solve the NSE equations for a range of proton-rich thermodynamic conditions. We show that cold proton-rich NSE is qualitatively different from neutron-rich NSE. Instead of being dominated by the Fe-peak nuclei with the largest binding energy per nucleon that have a proton to nucleon ratio close to the prescribed electron fraction, NSE for proton-rich material near freeze-out temperature is mainly composed of Ni56 and free protons. Previous results of nuclear reaction network calculations rely on this non-intuitive high proton abundance, which this paper will explain. We show how the differences and especially the large fraction of free protons arises from the minimization of the free energy as a result of a delicate competition between the entropy and the nuclear binding energy.

Prior to the explosion of a carbon-oxygen white dwarf in a Type Ia supernova there is a long "simmering," during which the 12C + 12C reaction gradually heats the white dwarf on a long (~ 1000 yr) timescale. Piro & Bildsten showed that weak reactions during this simmering set a maximum electron abundance Ye at the time of the explosion. We investigate the nuclear reactions during this simmering with a series of self-heating, at constant pressure, reaction network calculations. Unlike in AGB stars, proton captures onto 22Ne and heavier trace nuclei do not play a significant role. The 12C abundance is sufficiently high that the neutrons preferentially capture onto 12C, rather than iron group nuclei. As an aid to hydrodynamical simulations of the simmering phase, we present fits to the rates of heating, electron capture, change in mean atomic mass, and consumption of 12C in terms of the screened thermally averaged cross section for 12C + 12C. Our evaluation of the net heating rate includes contributions from electron captures into the 3.68 MeV excited state of 13C. This results in a slightly larger energy release, per 12C consumed, than that found by Piro & Bildsten, but less than that released for a burn to only 20Ne and 23Na. We compare our one-zone results to more accurate integrations over the white dwarf structure to estimate the amount of 12C that must be consumed to raise the white dwarf temperature, and hence to determine the net reduction of Ye during simmering.

We review the current state of modeling convective mixing in AGB stars. The focus is on results obtained through multi-dimensional hydrodynamic simulations of AGB convection, both in the envelope and the unstable He-shell. Using two different codes and a wide range of resolutions and modeling assumptions we find that mixing across convective boundaries is significant for He-shell flash convection. We present a preliminary quantitative analysis of this convectively induced extra mixing, based on a sub-set of our simulations. Other non-standard mixing will be discussed briefly.

Classical novae participate in the cycle of Galactic chemical evolution in which grains and metal enriched gas in their ejecta, supplementing those of supernovae, AGB stars, and Wolf-Rayet stars, are a source of heavy elements for the ISM. Once in the diffuse gas, this material is mixed with the existing gases and then incorporated into young stars and planetary systems during star formation. Infrared observations have confirmed the presence of carbon, SiC, hydrocarbons, and oxygen-rich silicate grains in nova ejecta, suggesting that some fraction of the pre-solar grains identified in meteoritic material come from novae. The mean mass returned by a nova outburst to the ISM probably exceeds ~2 x 10^{-4} Solar Masses. Using the observed nova rate of 35 per year in our Galaxy, it follows that novae introduce more than ~7 x 10^{-3} Solar Masses per year of processed matter into the ISM. Novae are expected to be the major source of 15N and 17O in the Galaxy and to contribute to the abundances of other isotopes in this atomic mass range. Here, we report on how changes in the nuclear reaction rates affect the properties of the outburst and alter the predictions of the contributions of novae to Galactic chemical evolution. We also discuss the necessity of including the pep reaction in studies of thermonuclear runaways in material accreted onto white dwarfs.

Carbon-oxygen white dwarfs contain neon-22 formed from alpha-captures onto nitrogen during core He burning in the progenitor star. In a white dwarf (type Ia) supernova, the neon-22 abundance determines, in part, the neutron-to-proton ratio and hence the abundance of radioactive nickel-56 that powers the lightcurve. The neon-22 abundance also changes the burning rate and hence the laminar flame speed. We tabulate the flame speedup for different initial carbon and neon-22 abundances and for a range of densities. This increase in the laminar flame speed--about 30% for a neon-22 mass fraction of 6%--affects the deflagration just after ignition near the center of the white dwarf, where the laminar speed of the flame dominates over the buoyant rise, and in regions of lower density ~ 10^7 g/cm3 where a transition to distributed burning is conjectured to occur. The increase in flame speed will decrease the density of any transition to distributed burning.

The conditions for the leading r-process site candidate, neutrino-driven winds, can not be reproduced self-consistently in current supernova models. For that reason, we investigate an alternate model involving the mass ejected by fallback in a supernova explosion, through hydrodynamic and nucleosynthesis calculations. The nucleosynthetic products of this ejected material produces r-process elements, including those in the vicinity of the elusive 3rd peak at mass number 195. Trans-iron element production beyond the second peak is made possible by a rapid (<1ms) freezeout of alpha particles which leaves behind a large nucleon (including protons!) to r-process seed ratio. This rapid phase is followed by a relatively long (>15 ms) simmering phase at ~2e9 K, which is the thermodynamic consequence of the hydrodynamic trajectory of the turbulent flows in the fallback outburst. During the slow phase high mass elements beyond the second peak are first made through rapid capture of both protons and neutrons. The flow stays close to valley of stability during this phase. After freeze-out of protons the remaining neutrons cause a shift out to short-lived isotopes as is typical for the r-process. A low electron fraction isn't required in this model, however, the detailed final distribution is sensitive to the electron fraction. Our simulations suggest that supernova fallback is a viable alternative scenario for the r-process.

We present the first hydrodynamic, multi-dimensional simulations of He-shell flash convection. Specifically, we investigate the properties of shell convection at a time immediately before the He- luminosity peak during the 15th thermal pulse of a stellar evolution track with initially two solar masses and metallicity Z=0.01. This choice is a representative example of a low-mass asymptotic giant branch thermal pulse. We construct the initial vertical stratification with a set of polytropes to resemble the stellar evolution structure. Convection is driven by a constant volume heating in a thin layer at the bottom of the unstable layer. We calculate a grid of 2D simulations with different resolutions and heating rates. Our set of simulations includes one low-resolution 3D run. The computational domain includes 11.4 pressure scale heights. He-shell flash convection is dominated by large convective cells that are centered in the lower half of the convection zone. Convective rolls have an almost circular appearance because focusing mechanisms exist in the form of the density stratification for downdrafts and the heating of localized eddies that generate upflows. Nevertheless, downdrafts appear to be somewhat more focused. The He-shell flash convection generates a rich spectrum of gravity waves in both stable layers above and beneath the convective shell. The magnitude of the convective velocities from our 1D mixing-length theory model and the rms-averaged vertical velocities from the hydrodynamic model are consistent within a factor of a few. However, the velocity profile in the hydrodynamic simulation is more asymmetric, and decays exponentially inside the convection zone. [abbreviated]

The details of ignition of Type Ia supernovae remain fuzzy, despite the importance of this input for any large-scale model of the final explosion. Here, we begin a process of understanding the ignition of these hotspots by examining the burning of one zone of material, and then investigate the ignition of a detonation due to rapid heating at single point. We numerically measure the ignition delay time for onset of burning in mixtures of degenerate material and provide fitting formula for conditions of relevance in the Type Ia problem. Using the neon abundance as a proxy for the white dwarf metallicity, we then find that ignition times can decrease by ~20% with addition of even 5% of neon by mass. When temperature fluctuations that successfully kindle a region are very rare, such a reduction in ignition time can increase the probability of ignition by orders of magnitude. If the neon comes largely at the expense of carbon, a similar increase in the ignition time can occur. We then consider the ignition of a detonation by an explosive energy input in one localized zone, eg a Sedov blast wave leading to a shock-ignited detonation. Building on previous work on curved detonations, we find that surprisingly large inputs of energy are required to successfully launch a detonation, leading to required matchheads of ~4500 detonation thicknesses - tens of centimeters to hundreds of meters - which is orders of magnitude larger than naive considerations might suggest. This is a very difficult constraint to meet for some pictures of a deflagration-to-detonation transition, such as a Zel'dovich gradient mechanism ignition in the distributed burning regime.

We compare a suite of 3D explosion calculations and stellar models incorporating advanced physics with observational constraints on the progenitor of Cassiopeia A. We consider binary and single stars from 16 to 40 solar masses with a range of explosion energies and geometries. The parameter space allowed by observations of nitrogen rich high velocity ejecta, ejecta mass, compact remnant mass, and 44Ti and 56Ni abundances individually and as an ensemble is considered. A progenitor of 15-25 solar masses which loses its hydrogen envelope to a binary interaction and undergoes an energetic explosion can match all the observational constraints.

In recent years proper motion measurements have been added to the set of observational constraints on the current properties of Galactic X-ray binaries. We develop an analysis that allows us to consider all this available information and reconstruct the full evolutionary history of X-ray binaries back to the time of compact object formation. This analysis accounts for mass transfer through the ongoing X-ray phase, tidal circularization before the onset of Roche-lobe overflow, motion through the Galactic potential after the formation of the compact object, and binary orbital dynamics and hydrodynamic modeling of the core collapse. We apply the analysis to the soft X-ray transient GRO J1655-40 and, for the first time, use its full 3D peculiar velocity constraints right after core collapse instead of lower limits on the current space velocity given by the present-day radial velocity. We find that the system has remained within 200pc from the Galactic plane throughout its entire life time and that the mass loss and a kick possibly associated with the black hole (BH) formation imparted a kick velocity of 45-115 km/s to the binary's center of mass. Right after BH formation, the system consists of a 3.5-6.3 Msun BH and a 2.3-4 Msun main-sequence star. At the onset of the X-ray phase the donor is still on the main sequence. We find that a symmetric BH formation event cannot be formally excluded, but that the associated system parameters are only marginally consistent with the currently observed binary properties. BH formation mechanisms involving an asymmetric supernova explosion with associated BH kick velocities of a few tens of km/s, on the other hand, satisfy the constraints much more comfortably. We also derive an upper limit on the BH kick magnitude of 210 km/s. (abridged)

Nova explosions occur on the white dwarf (WD) component of a Cataclysmic Variable stellar system which is accreting matter lost by a companion. A Type Ia supernova explosion is thought to result when a WD, in a similar binary configuration, grows in mass to the Chandrasekhar Limit. Here, we present calculations of accretion of Solar matter, at a variety of mass accretion rates, onto hot ($2.3 \times 10^{5}$K), luminous (30L$_\odot$), massive (1.25M$_\odot$, 1.35M$_\odot$) Carbon-Oxygen WDs. In contrast to our nova simulations where the WD has a low initial luminosity and a thermonuclear runaway (TNR) occurs and ejects material, these simulations do not eject material (or only a small fraction of the accreted material) and the WD grows in mass. A hydrogen TNR does not occur because hydrogen fuses to helium in the surface layers, and we call this process Surface Hydrogen Burning (SHB). As the helium layer grows in mass, it gradually fuses either to carbon and oxygen or to more massive nuclei depending on the WD mass and mass accretion rate. If such a WD were to explode in a SN Ia event, therefore, it would show neither hydrogen nor helium in its spectrum as is observed. Moreover, the luminosities and effective temperatures of our simulations agree with the observations of some of the Super Soft X-ray Binary Sources and, therefore, our results strengthen previous speculation that some of them (CAL 83 and CAL 87 for example) are probably progenitors of SN Ia explosions. Finally, we have achieved SHB for values of the mass accretion rate that almost span the observed values of the Cataclysmic Variables.

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.

Critically understanding the `standard candle'-like behavior of Type Ia supernovae requires understanding their explosion mechanism. One family of models for Type Ia Supernovae begins with a deflagration in a Carbon-Oxygen white dwarf which greatly accelerates through wrinkling and flame instabilities. While the planar speed and behavior of astrophysically-relevant flames is increasingly well understood, more complex behavior, such as the flame's response to stretch and curvature, has not been extensively explored in the astrophysical literature; this behavior can greatly enhance or suppress instabilities and local flame-wrinkling, which in turn can increase or decrease the bulk burning rate. In this paper, we explore the effects of curvature on both nuclear flames and simpler model flames to understand the effect of curvature on the flame structure and speed.

We explore the idea that the observed variations in the peak luminosities of Type Ia supernovae originate in part from a scatter in metallicity of the main-sequence stars that become white dwarfs. Previous, numerical, studies have not self-consistently explored metallicities greater than solar. One-dimensional Chandrasekhar mass models of SNe Ia produce most of their 56Ni in a burn to nuclear statistical equilibrium between the mass shells 0.2 and 0.8 solar masses, for which the electron to nucleon ratio is constant during the burn. We show analytically that, under these conditions, charge and mass conservation constrain the mass of 56Ni produced to depend linearly on the original metallicity of the white dwarf progenitor. Detailed post-processing of W7-like models confirms this linear dependence. The effect that we identify is most evident at metallicities larger than solar, and is in agreement with previous self-consistent calculations over the metallicity range common to both calculations. The observed scatter in the metallicity (1/3--3 times solar) of the solar neighborhood is enough to induce a 25% variation in the mass of 56Ni ejected by Type Ia supernovae. This is sufficient to vary the peak V-band brightness by approximately 0.2. This scatter in metallicity is present out to the limiting redshifts of current observations (z < 1). Sedimentation of 22Ne can possibly amplify the variation in 56Ni mass up to 50%. Further numerical studies can determine if other metallicity-induced effects, such as a change in the mass of the 56Ni-producing region, offset or enhance this variation.

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

We look in detail at the process of mapping an astrophysical initial model from a stellar evolution code onto the computational grid of an explicit, Godunov type code while maintaining hydrostatic equilibrium. This mapping process is common in astrophysical simulations, when it is necessary to follow short-timescale dynamics after a period of long timescale buildup. We look at the effects of spatial resolution, boundary conditions, the treatment of the gravitational source terms in the hydrodynamics solver, and the initialization process itself. We conclude with a summary detailing the mapping process that yields the lowest ambient velocities in the mapped model.

Observed novae abundances and explosion energies estimated from observations indicate that there must be significant mixing of the heavier material of the white dwarf (C+O) into the lighter accreted material (H+He). Accordingly, nova models must incorporate a mechanism that will dredge up the heavier white dwarf material, and fluid motions from an early convection phase is one proposed mechanism. We present results from two-dimensional simulations of classical nova precursor models that demonstrate the beginning of a convective phase during the `simmering' of a Nova precursor. We use a new hydrostatic equilibrium hydrodynamics module recently developed for the adaptive-mesh code FLASH. The two-dimensional models are based on the one-dimensional models of Ami Glasner, and were evolved with FLASH from a pre-convective state to the onset of convection.

We present the results of a simulation of a wind-driven non-linear gravity wave breaking on the surface of a white dwarf. The ``wind'' consists of H/He from an accreted envelope, and the simulation demonstrates that this breaking wave mechanism can produce a well-mixed layer of H/He with C/O from the white dwarf above the surface. Material from this mixed layer may then be transported throughout the accreted envelope by convection, which would enrich the C/O abundance of the envelope as is expected from observations of novae.

We present a case study of validating an astrophysical simulation code. Our study focuses on validating FLASH, a parallel, adaptive-mesh hydrodynamics code for studying the compressible, reactive flows found in many astrophysical environments. We describe the astrophysics problems of interest and the challenges associated with simulating these problems. We describe methodology and discuss solutions to difficulties encountered in verification and validation. We describe verification tests regularly administered to the code, present the results of new verification tests, and outline a method for testing general equations of state. We present the results of two validation tests in which we compared simulations to experimental data. The first is of a laser-driven shock propagating through a multi-layer target, a configuration subject to both Rayleigh-Taylor and Richtmyer-Meshkov instabilities. The second test is a classic Rayleigh-Taylor instability, where a heavy fluid is supported against the force of gravity by a light fluid. Our simulations of the multi-layer target experiments showed good agreement with the experimental results, but our simulations of the Rayleigh-Taylor instability did not agree well with the experimental results. We discuss our findings and present results of additional simulations undertaken to further investigate the Rayleigh-Taylor instability.

We present direct numerical simulations of mixing at Rayleigh-Taylor unstable interfaces performed with the FLASH code, developed at the ASCI/Alliances Center for Astrophysical Thermonuclear Flashes at the University of Chicago. We present initial results of single-mode studies in two and three dimensions. Our results indicate that three-dimensional instabilities grow significantly faster than two-dimensional instabilities and that grid resolution can have a significant effect on instability growth rates. We also find that unphysical diffusive mixing occurs at the fluid interface, particularly in poorly resolved simulations.

We discuss some of the computational challenges encountered in simulating the evolution of clusters of galaxies. Eulerian adaptive mesh refinement (AMR) techniques can successfully address these challenges but are currently being used by only a few groups. We describe our publicly available AMR code, FLASH, which uses an object-oriented framework to manage its AMR library, physics modules, and automated verification. We outline the development of the FLASH framework to include collisionless particles, permitting it to be used for cluster simulation.

The burning regimes encountered by laminar deflagrations and ZND detonations propagating through helium-rich compositions in the presence of buoyancy-driven turbulence are analyzed. Particular attention is given to models of X-ray bursts which start with a thermonuclear runaway on the surface of a neutron star, and the thin shell helium instability of intermediate-mass stars. In the X-ray burst case, turbulent deflagrations propagating in the lateral or radial directions encounter a transition from the distributed regime to the flamlet regime at a density of 10^8 g cm^{-3}. In the radial direction, the purely laminar deflagration width is larger than the pressure scale height for densities smaller than 10^6 g cm^{-3}. Self-sustained laminar deflagrations travelling in the radial direction cannot exist below this density. Similarily, the planar ZND detonation width becomes larger than the pressure scale height at 10^7 g cm^{-3}, suggesting that a steady-state, self-sustained detonations cannot come into existance in the radial direction. In the thin helium shell case, turbulent deflagrations travelling in the lateral or radial directions encounter the distributed regime at densities below 10^7 g cm^{-3}, and the flamelet regime at larger densities. In the radial direction, the purely laminar deflagration width is larger than the pressure scale height for densities smaller than 10^4 g cm^{-3}, indicating that steady-state laminar deflagrations cannot form below this density. The planar ZND detonation width becomes larger than the pressure scale height at 5 10^4 g cm^{-3}, suggesting that steady-state, self-sustained detonations cannot come into existance in the radial direction.

Gamma-rays from the decay of $^{26}$Al offer a stringent constraint on the Galaxy's global star formation rate over the past million years, supplementing other methods for quantifying the recent Galactic star formation rate, such as equivalent widths of H$\alpha$ emission. Advantages and disadvantages of using $^{26}$Al gamma-ray measurements as a tracer of the massive star formation rate are analyzed. Estimates of the Galactic $^{26}$Al mass derived from COMPTEL measurements are coupled with a simple, analytical model of the $^{26}$Al injection rate from massive stars and restrict the Galaxy's recent star formation rate to \hbox{5 $\pm$ 4 M\sun yr$^{-1}$}. In addition, we show that the derived $^{26}$Al mass implies a present day \hbox{Type II + Ib} supernovae rate of 3.4 $\pm$ 2.8 per century, which seems consistent with other independent estimates of the Galactic core collapse supernova rate. If some independent measure of the massive star initial mass function or star formation rate or \hbox{Type II + Ib} supernovae rate were to become available (perhaps through estimates of the Galactic $^{60}$Fe mass), then a convenient way to restrain, or possibly determine, the other parameters is presented.

We examine the consequences and implications of the possibility that the best-fit $m$=4/3 line of the silicon isotopic ratios measured in mainstream SiC grains is identical or parallel to to the mean ISM evolution line of the silicon isotopes. Even though the mean ISM evolution proceeds along a line of unity slope when deviations are expressed in terms of the native representation (the mean ISM), the evolution line can become a slope 4/3 line in the solar representation, provided that the solar composition is displaced from the mean ISM evolution. During the course of this analysis, we introduce new methods for relating the solar composition to that of the mean ISM at the time of solar birth. These new developments offer a unique view on the meaning of the mainstream SiC particles, and affords a new way of quantitatively answering the question whether the sun has a special composition relative to the mean ISM at solar birth. If the correlation slope of the silicon isotopes in the mean ISM could be decisively established, then its value would quantify the difference between the solar and mean ISM silicon abundances. Our formalism details the transformations between the two representations, and applies not only to $^{29}$Si and $^{30}$Si, but to any two purely secondary isotopes of any element (O, Ne, Mg, and perhaps S). Both the advantages and disadvantages of this technique are critically reviewed.

Plausible ratios of deuterium to hydrogen D/H as a function of metallicity, time, and redshift are investigated. Guided by the heavy element abundance patterns observed locally in Galactic dwarf stars and at large redshift in quasi-stellar object absorption line systems, empirical evolution of the relative abundance ratios Li/D, B/D, N/D, O/D, and F/D for QSO absorption line systems are given for the possible evolutionary patterns in D/H. Shortened abstract.

We calculate and discuss the chemical evolution of the isotopic silicon abundances in the interstellar medium at distances and times appropriate to the birth of the solar system. This has several objectives, some of which are related to anomalous silicon isotope ratios within presolar grains extracted from meteorites; namely: (1) What is the relative importance for silicon isotopic compositions in the bulk ISM of Type II supernovae, Type Ia supernovae, and AGB stars? (2) Are $^{29}$Si and $^{30}$Si primary or secondary nucleosynthesis products? (3) In what isotopic direction in a three-isotope plot do core-collapse supernovae of different mass move the silicon isotopic composition? (4) Why do present calculations not reproduce the solar ratios for silicon isotopes, and what does that impose upon studies of anomalous Si isotopes in meteoritic silicon carbide grains? (5) Are chemical-evolution features recorded in the anomalous SiC grains? Our answers are formulated on the basis of the Woosley \& Weaver (1995) supernova yield survey. Renormalization with the calculated interstellar medium silicon isotopic composition and solar composition is as an important and recurring concept of this paper. Possible interpretations of the silicon isotope anomalies measured in single SiC grains extracted from carbonaceous meteorites are then presented. The calculations suggest that the temporal evolution of the isotopic silicon abundances in the interstellar medium may be recorded in these grains.

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

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

Using recently calculated models for massive stellar evolution and supernovae coupled to a model for Galactic chemical evolution, neutron star and black hole birth functions (number of neutron stars and black holes as a function of their mass) are determined for the Milky Way Galaxy. For those stars that explode as Type II supernovae, the models give birth functions that are bimodal with peaks at 1.27 and 1.76 M\sun \ and average masses within those peaks of 1.28 and 1.73 M\sun. For those stars that explode as Type Ib there is a narrower spread of remnant masses, the average being 1.32 M\sun, and less evidence for bimodality. These values will be increased, especially in the more massive Type II supernovae, if significant accretion continues during the initial launching of the shock, and the number of heavier neutron stars could be depleted by black hole formation. The principal reason for the dichotomy in remnant masses for Type II is the difference in the presupernova structure of stars above and below 19 M\sun, the mass separating stars that burn carbon convectively from those that produce less carbon and burn radiatively. The Type Ib's and the lower mass group of the Type II's compare favorably with measured neutron star masses, and in particular to the Thorsett et al. (1993) determination of the average neutron star mass in 17 systems; 1.35 $\pm$ 0.27 M\sun. Variations in the exponent of a Salpeter initial mass function are shown not to affect the locations of the two peaks in the distribution function, but do affect their relative amplitudes. Sources of uncertainty, in particular placement of the mass cut and sensitivity to the explosion energy, are discussed, and estimates of the total number of neutron stars and black holes in the Galaxy are given. Accretion induced collapse should give a unique gravitational mass of 1.27

Abundance histories for QSO absorption line systems as a function of redshift are presented for all metals lighter than gallium. Coupling various conventional cosmological models with a simple, first--order model for the chemical evolution of the QSO absorption line systems allows transformation of the observed abundance trends in the Galaxy to abundance histories for the gas in QSO absorption line systems. Comparison of the transformed abundance patterns with the zinc to hydrogen [Zn/H] and silicon to hydrogen [Si/H] ratios observed in damped Lyman-$\alpha$ systems finds agreement over more than two orders of magnitude in abundance with a $\Lambda$=0, \ $\Omega$=0.2 -- 1.0, \ and $\tau_{\rm {delay}}$=3 Gyr cosmology, where $\tau_{\rm {delay}}$ is the time between the start of the Big Bang and the beginning of galaxy formation. Alternative meanings for the required time delay are explored, along with extensions to the simple first--order model. Predictions are made for the abundance histories of elements not yet observed in QSO absorption line systems.

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