Search

The observation of X-rays during quiescence from transiently accreting neutron stars provides unique clues about the nature of dense matter. This, however, requires extensive modeling of the crusts and matching the results to observations. The pycnonuclear fusion reaction rates implemented in these models are theoretically calculated by extending phenomenological expressions and have large uncertainties spanning many orders of magnitude. We present the first sensitivity studies of these pycnonuclear fusion reactions in realistic network calculations. We also couple the reaction network with the thermal evolution code dStar to further study their impact on the neutron star cooling curves in quiescence. Varying the pycnonuclear fusion reaction rates alters the depth at which nuclear heat is deposited although the total heating remains constant. The enhancement of the pycnonuclear fusion reaction rates leads to an overall shallower deposition of nuclear heat. The impurity factors are also altered depending on the type of ashes deposited on the crust. These total changes correspond to a variation of up to 9 eV in the modeled cooling curves. While this is not sufficient to explain the shallow heat source, it is comparable to the observational uncertainties and can still be important for modeling the neutron star crust., Comment: AASTeX63, 11 pages with 9 figures

Nuclear Astrophysics is a field at the intersection of nuclear physics and astrophysics, which seeks to understand the nuclear engines of astronomical objects and the origin of the chemical elements. This white paper summarizes progress and status of the field, the new open questions that have emerged, and the tremendous scientific opportunities that have opened up with major advances in capabilities across an ever growing number of disciplines and subfields that need to be integrated. We take a holistic view of the field discussing the unique challenges and opportunities in nuclear astrophysics in regards to science, diversity, education, and the interdisciplinarity and breadth of the field. Clearly nuclear astrophysics is a dynamic field with a bright future that is entering a new era of discovery opportunities., Comment: 96 pages. Submitted to Journal of Physics G

Evidence for the presence of short-lived radioactive isotopes when the Solar System formed is preserved in meteorites, providing insights into the conditions at the birth of our Sun. A low-mass core-collapse supernova had been postulated as a candidate for the origin of $^{10}$Be, reinforcing the idea that a supernova triggered the formation of the Solar System. We present a detailed study of the production of $^{10}$Be by the $\nu$ process in supernovae, which is very sensitive to the reaction rate of the major destruction channel, $^{10}$Be(p,$\alpha$)$^{7}$Li. With data from recent nuclear experiments that show the presence of a resonant state in $^{11}$B at $\approx$~193 keV, we derive new values for the $^{10}$Be(p,$\alpha$)$^{7}$Li reaction rate which are significantly higher than previous estimates. We show that, with the new $^{10}$Be(p,$\alpha$)$^{7}$Li reaction rate, a low mass CCSN is unlikely to produce enough $^{10}$Be to explain the observed $^{10}$Be/ $^{9}$Be ratio in meteorites, even for a wide range of neutrino spectra considered in our models. These findings point towards spallation reactions induced by solar energetic particles in the early Solar System as the origin of $^{10}$Be., Comment: 12 pages, 9 Figures, 4 Tables. Submitted to Physical Review C

Nuclear reactions heat and cool the crust of accreting neutron stars and need to be understood to interpret observations of X-ray bursts and of long-term cooling in transiently accreting systems. It was recently suggested that previously neglected neutron transfer reactions may play a significant role in the nuclear processes. We present results from full nuclear network calculations that now include these reactions and determine their impact on crust composition, crust impurity, heating, and cooling. We find that a large number of neutron transfer reactions indeed occur and impact crust models. In particular, we identify a new type of reaction cycle that brings a pair of nuclei across the nuclear chart into equilibrium via alternating neutron capture and neutron release, interspersed with a neutron transfer. While neutron transfer reactions lead to changes in crust model predictions, and need to be considered in future studies, previous conclusions concerning heating, cooling, and compositional evolution are remarkably robust., Comment: 8 pages, 10 figures, accepted for Astrophysical Journal

Cataclysmic Variables (CVs) and Symbiotic Binaries are close (or not so close) binary star systems which contain both a white dwarf (WD) primary and a larger cooler secondary star that typically fills its Roche Lobe. The cooler star is losing mass through the inner Lagrangian point of the binary and a fraction of this material is accreted by the WD. Here we report on our hydrodynamic studies of the thermonuclear runaway (TNR) in the accreted material that ends in a Classical Nova explosion. We have followed the evolution of the TNRs on both carbon-oxygen (CO) and oxygen-neon (ONe) WDs. We report on 3 studies in this paper. First, simulations in which we accrete only solar matter using NOVA (our 1-D, fully implicit, hydro code). Second, we use MESA for similar studies in which we accrete only Solar matter and compare the results. Third, we accrete solar matter until the TNR is ongoing and then switch the composition in the accreted layers to a mixed composition: either 25% WD and 75% solar or 50% WD and 50% Solar matter. We find that the amount of accreted material is inversely proportional to the initial 12C abundance (as expected). Thus, accreting solar matter results in a larger amount of accreted material to fuel the outburst; much larger than in earlier studies where a mixed composition was assumed from the beginning of the simulation. Our most important result is that all these simulations eject significantly less mass than accreted and, therefore, the WD is growing in mass toward the Chandrasekhar Limit., Comment: To appear in "The Golden Age of Cataclysmic Variables and Related Objects V (GOLDEN2019)" 2-7 September 2019, Palermo, Italy. Proceedings of Science, in press

X-ray observations of transiently accreting neutron stars during quiescence provide information about the structure of neutron star crusts and the properties of dense matter. Interpretation of the observational data requires an understanding of the nuclear reactions that heat and cool the crust during accretion, and define its nonequilibrium composition. We identify here in detail the typical nuclear reaction sequences down to a depth in the inner crust where the mass density is 2E12 g/cm^3 using a full nuclear reaction network for a range of initial compositions. The reaction sequences differ substantially from previous work. We find a robust reduction of crust impurity at the transition to the inner crust regardless of initial composition, though shell effects can delay the formation of a pure crust somewhat to densities beyond 2E12 g/cm^3. This naturally explains the small inner crust impurity inferred from observations of a broad range of systems. The exception are initial compositions with A >= 102 nuclei, where the inner crust remains impure with an impurity parameter of Qimp~20 due to the N = 82 shell closure. In agreement with previous work we find that nuclear heating is relatively robust and independent of initial composition, while cooling via nuclear Urca cycles in the outer crust depends strongly on initial composition. This work forms a basis for future studies of the sensitivity of crust models to nuclear physics and provides profiles of composition for realistic crust models., Comment: 25 Pages, accepted for publication in Ap. J

In this work we report briefly on the gravitational wave (GW) signal computed in the context of a self-consistent, 3D simulation of a core-collapse supernova (CCSN) explosion of a 15M$_\odot$ progenitor star. We present a short overview of the GW signal, including signal amplitude, frequency distribution, and the energy emitted in the form of GWs for each phase of explosion, along with neutrino luminosities, and discuss correlations between them., Comment: 5 pages, 2 figures, submitted in the proceedings of the 52nd Rencontres de Moriond

Nova explosions occur on the white dwarf component of a Cataclysmic Variable binary stellar system that is accreting matter lost by its companion. When sufficient material has been accreted by the white dwarf, a thermonuclear runaway occurs and ejects material in what is observed as a Classical Nova explosion. We describe both the recent advances in our understanding of the progress of the outburst and outline some of the puzzles that are still outstanding. We report on the effects of improving both the nuclear reaction rate library and including a modern nuclear reaction network in our one-dimensional, fully implicit, hydrodynamic computer code. In addition, there has been progress in observational studies of Supernovae Ia with implications about the progenitors and we discuss that in this review., Comment: 18 pages, 7 figures

Core-collapse supernovae, the culmination of massive stellar evolution, are spectacular astronomical events and the principle actors in the story of our elemental origins. Our understanding of these events, while still incomplete, centers around a neutrino-driven central engine that is highly hydrodynamically unstable. Increasingly sophisticated simulations reveal a shock that stalls for hundreds of milliseconds before reviving. Though brought back to life by neutrino heating, the development of the supernova explosion is inextricably linked to multi-dimensional fluid flows. In this paper, the outcomes of three-dimensional simulations that include sophisticated nuclear physics and spectral neutrino transport are juxtaposed to learn about the nature of the three dimensional fluid flow that shapes the explosion. Comparison is also made between the results of simulations in spherical symmetry from several groups, to give ourselves confidence in the understanding derived from this juxtaposition., Comment: 13 pages, 4 figures, to appear in proceedings of XXXIV Mazurian Lakes Conference on Physics to be published in Acta Physica Polonica B

Simulations of r-process nucleosynthesis require nuclear physics information for thousands of neutron-rich nuclear species from the line of stability to the neutron drip line. While arguably the most important pieces of nuclear data for the r-process are the masses and beta decay rates, individual neutron capture rates can also be of key importance in setting the final r-process abundance pattern. Here we consider the influence of neutron capture rates in forming the A~80 and rare earth peaks., Comment: 10 pages, 5 figures, appears in the Proceedings of the 14th International Symposium on Capture Gamma-Ray Spectroscopy and Related Topics

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., Comment: 11 pages 10 figures in the 2011 Asiago Symbiotic Workshop

The Spallation Neutron Source (SNS) at Oak Ridge National Laboratory, Tennessee, provides an intense flux of neutrinos in the few tens-of-MeV range, with a sharply-pulsed timing structure that is beneficial for background rejection. In this document, the product of a workshop at the SNS in May 2012, we describe this free, high-quality stopped-pion neutrino source and outline various physics that could be done using it. We describe without prioritization some specific experimental configurations that could address these physics topics., Comment: White paper associated with the Workshop on Neutrinos at the Spallation Neutron Source, May 2012

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., Comment: 25 Pages, Bulletin of the Astronomical Society of India (BASI) in press

In two preceding papers we have shown that, when reaction networks are well-removed from equilibrium, explicit asymptotic and quasi-steady-state approximations can give algebraically-stabilized integration schemes that rival standard implicit methods in accuracy and speed for extremely stiff systems. However, we also showed that these explicit methods remain accurate but are no longer competitive in speed as the network approaches equilibrium. In this paper we analyze this failure and show that it is associated with the presence of fast equilibration timescales that neither asymptotic nor quasi-steady-state approximations are able to remove efficiently from the numerical integration. Based on this understanding, we develop a partial equilibrium method to deal effectively with the approach to equilibrium and show that explicit asymptotic methods, combined with the new partial equilibrium methods, give an integration scheme that plausibly can deal with the stiffest networks, even in the approach to equilibrium, with accuracy and speed competitive with that of implicit methods. Thus we demonstrate that such explicit methods may offer alternatives to implicit integration of even extremely stiff systems, and that these methods may permit integration of much larger networks than have been possible before in a number of fields., Comment: Updated references

We show that, even for extremely stiff systems, explicit integration may compete in both accuracy and speed with implicit methods if algebraic methods are used to stabilize the numerical integration. The required stabilizing algebra depends on whether the system is well-removed from equilibrium or near equilibrium. This paper introduces a quantitative distinction between these two regimes and addresses the former case in depth, presenting explicit asymptotic methods appropriate when the system is extremely stiff but only weakly equilibrated. A second paper examines quasi-steady-state methods as an alternative to asymptotic methods in systems well away from equilibrium and a third paper extends these methods to equilibrium conditions in extremely stiff systems using partial equilibrium methods. All three papers present systematic evidence for timesteps competitive with implicit methods. Because an explicit method can execute a timestep faster than an implicit method, algebraically-stabilized explicit algorithms might permit integration of larger networks than have been feasible before in various disciplines., Comment: Updated references

The location of electron capture heat sources in the crust of accreting neutron stars depends on the masses of extremely neutron-rich nuclei. We present first results from a new implementation of the time-of-flight technique to measure nuclear masses of rare isotopes at the National Superconducting Cyclotron Laboratory. The masses of 16 neutron-rich nuclei in the scandium -- nickel range were determined simultaneously, improving the accuracy compared to previous data in 12 cases. The masses of $^{61}${V}, $^{63}${Cr}, $^{66}${Mn}, and $^{74}${Ni} were measured for the first time with mass excesses of $-30.510(890)$ MeV, $-35.280(650)$ MeV, $-36.900(790)$ MeV, and $-49.210(990)$ MeV, respectively. With the measurement of the $^{66}$Mn mass, the locations of the two dominant electron capture heat sources in the outer crust of accreting neutron stars that exhibit superbursts are now experimentally constrained. We find that the location of the $^{66}$Fe$\rightarrow^{66}$Mn electron capture transition occurs significantly closer to the surface than previously assumed because our new experimental Q-value is 2.1 MeV (2.6$\sigma$) smaller than predicted by the FRDM mass model., Comment: 4 pages, 3 figures

Unraveling the mechanism for core-collapse supernova explosions is an outstanding computational challenge and the problem remains essentially unsolved despite more than four decades of effort. However, much progress in realistic modeling has occurred recently through the availability of multi-teraflop machines and the increasing sophistication of supernova codes. These improvements have led to some key insights which may clarify the picture in the not too distant future. Here we briefly review the current status of the three explosion mechanisms (acoustic, MHD, and neutrino heating) that are currently under active investigation, concentrating on the neutrino heating mechanism as the one most likely responsible for producing explosions from progenitors in the mass range ~10 to ~25 solar masses. We then briefly describe the CHIMERA code, a supernova code we have developed to simulate core-collapse supernovae in 1, 2, and 3 spatial dimensions. We finally describe the results of an ongoing suite of 2D simulations initiated from a 12, 15, 20, and 25 solar mass progenitor. These have all exhibited explosions and are currently in the expanding phase with the shock at between 5,000 and 10,000 km. We finally very briefly describe an ongoing simulation in 3 spatial dimensions initiated from the 15 solar mass progenitor., Comment: 10 pages, 4 figures

Much progress in realistic modeling of core-collapse supernovae has occurred recently through the availability of multi-teraflop machines and the increasing sophistication of supernova codes. These improvements are enabling simulations with enough realism that the explosion mechanism, long a mystery, may soon be delineated. We briefly describe the CHIMERA code, a supernova code we have developed to simulate core-collapse supernovae in 1, 2, and 3 spatial dimensions. We then describe the results of an ongoing suite of 2D simulations initiated from a 12, 15, 20, and 25 solar mass progenitor. These have all exhibited explosions and are currently in the expanding phase with the shock at between 5,000 and 20,000 km. We also briefly describe an ongoing simulation in 3 spatial dimensions initiated from the 15 solar mass progenitor., Comment: 5 pages, 3 figures

Electron captures on nuclei play an important role in the dynamics of the collapsing core of a massive star that leads to a supernova explosion. Recent calculations of these capture rates were based on microscopic models which account for relevant degrees of freedom. Due to computational restrictions such calculations were limited to a modest number of nuclei, mainly in the mass range A=45-110. Recent supernova simulations show that this pool of nuclei, however, omits the very neutron-rich and heavy nuclei which dominate the nuclear composition during the last phase of the collapse before neutrino trapping. Assuming that the composition is given by Nuclear Statistical Equilibrium we present here electron capture rates for collapse conditions derived from individual rates for roughly 2700 individual nuclei. For those nuclei which dominate in the early stage of the collapse, the individual rates are derived within the framework of microscopic models, while for the nuclei which dominate at high densities we have derived the rates based on the Random Phase Approximation with a global parametrization of the single particle occupation numbers. In addition, we have improved previous rate evaluations by properly including screening corrections to the reaction rates into account., Comment: 32 pages, 13 figures, 1 table; elsart; to appear in Nuclear Physics A

The r-process constitutes one of the major challenges in nuclear astrophysics. Its astrophysical site has not yet been identified but there is observational evidence suggesting that at least two possible sites should contribute to the solar system abundance of r-process elements and that the r-process responsible for the production of elements heavier than Z=56 operates quite robustly producing always the same relative abundances. From the nuclear-physics point of view the r-process requires the knowledge of a large number of reaction rates involving exotic nuclei. These include neutron capture rates, beta-decays and fission rates, the latter for the heavier nuclei produced in the r-process. We have developed for the first time a complete database of reaction rates that in addition to neutron-capture rates and beta-decay half-lives includes all possible reactions that can induce fission (neutron-capture, beta-decay and spontaneous fission) and the corresponding fission yields. In addition, we have implemented these reaction rates in a fully implicit reaction network. We have performed r-process calculations for the neutrino-driven wind scenario to explore whether or not fission can contribute to provide a robust r-process pattern.

We discuss the results of 3D simulations of tidal disruptions of white dwarfs by moderate-mass black holes as they may exist in the cores of globular clusters or dwarf galaxies. Our simulations follow self-consistently the hydrodynamic and nuclear evolution from the initial parabolic orbit over the disruption to the build-up of an accretion disk around the black hole. For strong enough encounters (pericentre distances smaller than about 1/3 of the tidal radius) the tidal compression is reversed by a shock and finally results in a thermonuclear explosion. These explosions are not restricted to progenitor masses close to the Chandrasekhar limit, we find exploding examples throughout the whole white dwarf mass range. There is, however, a restriction on the masses of the involved black holes: black holes more massive than $2\times 10^5$ M$_\odot$ swallow a typical 0.6 M$_\odot$ dwarf before their tidal forces can overwhelm the star's self-gravity. Therefore, this mechanism is characteristic for black holes of moderate masses. The material that remains bound to the black hole settles into an accretion disk and produces an X-ray flare close to the Eddington limit of $L_{\rm Edd} \simeq 10^{41} {\rm erg/s} M_{\rm bh}/1000 M$_\odot$), typically lasting for a few months. The combination of a peculiar thermonuclear supernova together with an X-ray flare thus whistle-blows the existence of such moderate-mass black holes. The next generation of wide field space-based instruments should be able to detect such events., Comment: 8 pages, 2 figures, EuroWD08

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., Comment: Accepted to the Astrophysical Journal

In this paper we discuss opportunities for a neutrino program at the Spallation Neutrons Source (SNS) being commissioning at ORNL. Possible investigations can include study of neutrino-nuclear cross sections in the energy rage important for supernova dynamics and neutrino nucleosynthesis, search for neutrino-nucleus coherent scattering, and various tests of the standard model of electro-weak interactions., Comment: Paper submitted for the 2008 Carolina International Symposium on Neutrino Physics

We examine the role of neutron capture on 130Sn during r-process freeze-out in the neutrino-driven wind environment of the core-collapse supernova. We find that the global r-process abundance pattern is sensitive to the magnitude of the neutron capture cross section of 130Sn. The changes to the abundance pattern include not only a relative decrease in the abundance of 130Sn and an increase in the abundance of 131Sn, but also a shift in the distribution of material in the rare earth and third peak regions., Comment: 12 pages, 9 figures

We investigate the impact of neutron capture rates near the A=130 peak on the $r$-process abundance pattern. We show that these capture rates can alter the abundances of individual nuclear species, not only in the region of A=130 peak, but also throughout the abundance pattern. We discuss the nonequilibrium processes that produce these abundance changes and determine which capture rates have the most significant impact., Comment: 8 pages, 4 figures, submitted to Phys. Rev

We consider hot accretion disk outflows from black hole - neutron star mergers in the context of the nucleosynthesis they produce. We begin with a three dimensional numerical model of a black hole - neutron star merger and calculate the neutrino and antineutrino fluxes emitted from the resulting accretion disk. We then follow the element synthesis in material outflowing the disk along parameterized trajectories. We find that at least a weak r-process is produced, and in some cases a main r-process as well. The neutron-rich conditions required for this production of r-process nuclei stem directly from the interactions of the neutrinos emitted by the disk with the free neutrons and protons in the outflow., Comment: 10 pages, 4 figures, one table and additional references added

The existence of supermassive black holes lurking in the centers of galaxies and of stellar binary systems containing a black hole with a few solar masses has been established beyond reasonable doubt. The idea that black holes of intermediate masses ($\sim 1000$ \msun) may exist in globular star clusters has gained credence over recent years but no conclusive evidence has been established yet. An attractive feature of this hypothesis is the potential to not only disrupt solar-type stars but also compact white dwarf stars. In close encounters the white dwarfs can be sufficiently compressed to thermonuclearly explode. The detection of an underluminous thermonuclear explosion accompanied by a soft, transient X-ray signal would be compelling evidence for the presence of intermediate mass black holes in stellar clusters. In this paper we focus on the numerical techniques used to simulate the entire disruption process from the initial parabolic orbit, over the nuclear energy release during tidal compression, the subsequent ejection of freshly synthesized material and the formation process of an accretion disk around the black hole., Comment: 9 pages, 4 figures, Computer Physics Communications, accepted; movie can be found at http://www.faculty.iu-bremen.de/srosswog/; reference corrected

Suggestive evidence has accumulated that intermediate mass black holes (IMBH) exist in some globular clusters. As stars diffuse in the cluster, some will inevitable wander sufficiently close to the hole that they suffer tidal disruption. An attractive feature of the IMBH hypothesis is its potential to disrupt not only solar-type stars but also compact white dwarf stars. Attention is given to the fate of white dwarfs that approach the hole close enough to be disrupted and compressed to such extent that explosive nuclear burning may be triggered. Precise modeling of the dynamics of the encounter coupled with a nuclear network allow for a realistic determination of the explosive energy release, and it is argued that ignition is a natural outcome for white dwarfs of all varieties passing well within the tidal radius. Although event rates are estimated to be significantly less than the rate of normal Type Ia supernovae, such encounters may be frequent enough in globular clusters harboring an IMBH to warrant a search for this new class of supernova., Comment: 13 pages, 4 figures, ApJ, accepted

Here we present the results from two sets of simulations, in two and three spatial dimensions. In two dimensions, the simulations include multifrequency flux-limited diffusion neutrino transport in the "ray-by-ray-plus" approximation, two-dimensional self gravity in the Newtonian limit, and nuclear burning through a 14-isotope alpha network. The three-dimensional simulations are model simulations constructed to reflect the post stellar core bounce conditions during neutrino shock reheating at the onset of explosion. They are hydrodynamics-only models that focus on critical aspects of the shock stability and dynamics and their impact on the supernova mechanism and explosion. In two dimensions, we obtain explosions (although in one case weak) for two progenitors (11 and 15 Solar mass models). Moreover, in both cases the explosion is initiated when the inner edge of the oxygen layer accretes through the shock. Thus, the shock is not revived while in the iron core, as previously discussed in the literature. The three-dimensional studies of the development of the stationary accretion shock instability (SASI) demonstrate the fundamentally new dynamics allowed when simulations are performed in three spatial dimensions. The predominant l=1 SASI mode gives way to a stable m=1 mode, which in turn has significant ramifications for the distribution of angular momentum in the region between the shock and proto-neutron star and, ultimately, for the spin of the remnant neutron star. Moreover, the three-dimensional simulations make clear, given the increased number of degrees of freedom, that two-dimensional models are severely limited by artificially imposed symmetries., Comment: 9 pages, 3 figures

The overwhelming evidence that the core collapse supernova mechanism is inherently multidimensional, the complexity of the physical processes involved, and the increasing evidence from simulations that the explosion is marginal presents great computational challenges for the realistic modeling of this event, particularly in 3 spatial dimensions. We have developed a code which is scalable to computations in 3 dimensions which couples PPM Lagrangian with remap hydrodynamics [1], multigroup, flux-limited diffusion neutrino transport [2], with many improvements), and a nuclear network [3]. The neutrino transport is performed in a ray-by-ray plus approximation wherein all the lateral effects of neutrinos are included (e.g., pressure, velocity corrections, advection) except the transport. A moving radial grid option permits the evolution to be carried out from initial core collapse with only modest demands on the number of radial zones. The inner part of the core is evolved after collapse along with the rest of the core and mantle by subcycling the lateral evolution near the center as demanded by the small Courant times. We present results of 2-D simulations of a symmetric and an asymmetric collapse of both a 15 and an 11 M progenitor. In each of these simulations we have discovered that once the oxygen rich material reaches the shock there is a synergistic interplay between the reduced ram pressure, the energy released by the burning of the shock heated oxygen rich material, and the neutrino energy deposition which leads to a revival of the shock and an explosion., Comment: 10 pages, 3 figures

Recent halo star abundance observations exhibit an important feature of consequence to the r-process: the presence of a main r-process between the second and third peaks which is consistent among halo stars. We explore fission cycling and steady-beta flow as the driving mechanisms behind this feature. The presence of fission cycling during the r-process can account for nucleosynthesis yields between the second and third peaks, whereas the presence of steady-beta flow can account for consistent r-process patterns, robust under small variations in astrophysical conditions. We employ the neutrino-driven wind of the core-collapse supernova to examine fission cycling and steady-beta flow in the r-process. As the traditional neutrino-driven wind model does not produce the required very neutron-rich conditions for these mechanisms, we examine changes to the neutrino physics necessary for fission cycling to occur in the neutrino-driven wind environment, and we explore under what conditions steady-beta flow is obtained., Comment: 9 pages, 8 figures

Based on the shell model for Gamow-Teller and the Random Phase Approximation for forbidden transitions, we have calculated reaction rates for inelastic neutrino-nucleus scattering (INNS) under supernova (SN) conditions, assuming a matter composition given by Nuclear Statistical Equilibrium. The rates have been incorporated into state-of-the-art stellar core-collapse simulations with detailed energy-dependent neutrino transport. While no significant effect on the SN dynamics is observed, INNS increases the neutrino opacities noticeably and strongly reduces the high-energy tail of the neutrino spectrum emitted in the neutrino burst at shock breakout. Relatedly the expected event rates for the observation of such neutrinos by earthbound detectors are reduced by up to about 60%., Comment: 4 pages, 2 figures, 1 table

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., Comment: 9 pages, 2 figures, as it appeared in the Proceedings of the Tours 2006 Symposium on Nuclear Physics

It has long been suggested that fission cycling may play an important role in the r-process. Fission cycling can only occur in a very neutron rich environment. In traditional calculations of the neutrino driven wind of the core-collapse supernova, the environment is not sufficiently neutron rich to produce the r-process elements. However, we show that with a reduction of the electron neutrino flux coming from the supernova, fission cycling does occur and furthermore it produces an abundance pattern which is consistent with observed r-process abundance pattern in halo stars. Such a reduction can be caused by active-sterile neutrino oscillations or other new physics., Comment: Typos corrected. Presented at NIC-IX, International Symposium on Nuclear Astrophysics - Nuclei in the Cosmos - IX, CERN, Geneva, Switzerland, 25-30 June, 2006

We investigate nucleosynthesis in the supernovae post-core bounce neutrino-driven wind environment in the presence of active-sterile neutrino transformation. We consider active-sterile neutrino oscillations for a range of mixing parameters: vacuum mass-squared differences of 0.1 eV^2 < dm^2 < 100 eV^2, and vacuum mixing angles of sin^2(2 theta_v) > 10^-4. We find a consistent r-process pattern for a large range of mixing parameters that is in rough agreement with the halo star CS 22892-052 abundances and the pattern shape is determined by fission cycling. We find that the allowed region for the formation of the r-process peaks overlaps the LSND and NSBL (3+1) allowed region., Comment: 11 pages, 7 figures, Corrected Typos

Core collapse supernovae are the leading actor in the story of the cosmic origin of the chemical elements. Existing models, which generally assume spherical symmetry and parameterize the explosion, have been able to broadly replicate the observed elemental pattern. However, inclusion of neutrino interactions produces noticeable improvement in the composition of the ejecta when compared to observations. Neutrino interactions may also provide a supernova source for light p-process nuclei., Comment: 7 pages, 2 figures, in proceedings of Astronomy with Radioactivities V, Clemson University, September 5-9, 2005, to appear in New Astronomy Reviews

We present a new nucleosynthesis process, that we denote nu p-process, which occurs in supernovae (and possibly gamma-ray bursts) when strong neutrino fluxes create proton-rich ejecta. In this process, antineutrino absorptions in the proton-rich environment produce neutrons that are immediately captured by neutron-deficient nuclei. This allows for the nucleosynthesis of nuclei with mass numbers A >64. Making this process a possible candidate to explain the origin of the solar abundances of 92,94Mo and 96,98Ru. This process also offers a natural explanation for the large abundance of Sr seen in an hyper-metal-poor star., Comment: 5 pages, 3 figures, submitted to Physical Review Letters

Over the billions of years since the Big Bang, the lives, deaths and afterlives of stars have enriched the Universe in the heavy elements that make up so much of ourselves and our world. This review summarizes the methods used to evolve these nuclear abundances within astrophysical simulations. These methods fall into 2 categories; evolution via rate equations and via equilibria. Because the rate equations in nucleosynthetic applications involve a wide range of timescales, implicit methods have proven mandatory, leading to the need to solve matrix equations. Efforts to improve the performance of such rate equation methods are focused on efficient solution of these matrix equations, in particular by making best use of the sparseness of these matrices, and finding methods that require less frequent matrix solutions. Recent work to produce hybrid schemes which use local equilibria to reduce the computational cost of the rate equations is also discussed. Such schemes offer significant improvements in the speed of reaction networks and are accurate under circumstances where calculations which assume complete equilibrium fail., Comment: 27 pages, 2 figures, a review for a special issue of Nuclear Physics A

We examine the nucleosynthesis products that are produced in the outflow from rapidly accreting disks. We find that the type of element synthesis varies dramatically with the degree of neutrino trapping in the disk and therefore the accretion rate of the disk. Disks with relatively high accretion rates such as 10 M_solar/s can produce very neutron rich nuclei that are found in the r process. Disks with more moderate accretion rates can produce copious amounts of Nickel as well as the light elements such as Lithium and Boron. Disks with lower accretion rates such as 0.1 M_solar/s produce large amounts of Nickel as well as some unusual nuclei such as Ti-49, Sc-45, Zn-64, and Mo-92. This wide array of potential nucleosynthesis products is due to the varying influence of electron neutrinos and antineutrinos emitted from the disk on the neutron-to-proton ratio in the outflow. We use a parameterization for the outflow and discuss our results in terms of entropy and outflow acceleration., Comment: 12 pages, 12 figures; submitted to ApJ

With presently known input physics and computer simulations in 1D, a self-consistent treatment of core collapse supernovae does not yet lead to successful explosions, while 2D models show some promise. Thus, there are strong indications that the delayed neutrino mechanism works combined with a multi-D convection treatment for unstable layers. On the other hand there is a need to provide correct nucleosynthesis abundances for the progressing field of galactic evolution and observations of low metallicity stars. The innermost ejecta is directly affected by the explosion mechanism, i.e. most strongly the yields of Fe-group nuclei for which an induced piston or thermal bomb treatment will not provide the correct yields because the effect of neutrino interactions is not included. We apply parameterized variations to the neutrino scattering cross sections and alternatively, parameterized variations are applied to the neutrino absorption cross sections on nucleons in the ``gain region''. We find that both measures lead to similar results, causing explosions and a Ye>0.5 in the innermost ejected layers, due to the combined effect of a short weak interaction time scale and a negligible electron degeneracy, unveiling the proton-neutron mass difference. We include all weak interactions (electron and positron capture, beta-decay, neutrino and antineutrino capture on nuclei, and neutrino and antineutrino capture on nucleons) and present first nucleosynthesis results for these innermost ejected layers to discuss how they improve predictions for Fe-group nuclei. The proton-rich environment results in enhanced abundances of 45Sc, 49Ti, and 64Zn as requested by chemical evolution studies and observations of low metallicity stars as well as appreciable production of nuclei in the mass range up to A=80., Comment: 13 pages, 8 figures. Final version

We ensure successful explosions (of otherwise non-explosive models) by enhancing the neutrino luminosity via reducing the neutrino scattering cross sections or by increasing the heating efficiency via enhancing the neutrino absorption cross sections in the heating region. Our investigations show that the resulting electron fraction Ye in the innermost ejecta is close to 0.5, in some areas even exceeding 0.5. We present the effects of the resulting values for Ye on the nucleosynthesis yields of the innermost zones of core collapse supernovae., Comment: 4pages, 2figures; contribution to Nuclei In The Cosmos VIII, to appear in Nucl. Phys. A

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., Comment: Accepted by APJL, 4 pages, 1 figure, LaTex (uses emulateapj.sty)

We study neutral current neutrino-nucleus reaction cross-sections for Mn, Fe, Co and Ni isotopes. An earlier study for a few selected nuclei has shown that in the supernova environment the cross sections are increased for low energy neutrinos due to finite-temperature effects. Our work supports this finding for a much larger set of nuclei. Furthermore we extend previous work to higher neutrino energies considering allowed and forbidden multipole contributions to the cross sections. The allowed contributions are derived from large-scale shell model calculations of the Gamow-Teller strength, while the other multipole contributions are calculated within the Random Phase Approximation. We present the cross sections as functions of initial and final neutrino energies and for a range of supernova-relevant temperatures. These cross sections will allow improved estimates of inelastic neutrino reactions on nuclei in supernova simulations., Comment: 26 pages, 15 figures, submitted to Nuclear Physics A

The most important weak nuclear interaction to the dynamics of stellar core collapse is electron capture, primarily on nuclei with masses larger than 60. In prior simulations of core collapse, electron capture on these nuclei has been treated in a highly parameterized fashion, if not ignored. With realistic treatment of electron capture on heavy nuclei come significant changes in the hydrodynamics of core collapse and bounce. We discuss these as well as the ramifications for the post-bounce evolution in core collapse supernovae., Comment: Accepted by PRL, 5 pages, 2 figures

New facilities to measure neutrino-nucleus cross sections, such as those possible in conjunction with spallation neutron sources, could provide an experimental foundation for the many neutrino-nucleus weak interaction rates needed in supernova models. This would enable more realistic supernova models and provide a greatly improved ability to understand the physics fundamental to supernovae by comparison of these models with detailed observations. Charged- and neutral-current neutrino interactions on nuclei in the stellar core play a central role in supernova dynamics and nucleosynthesis as well as being important for supernova neutrino detection. Measurements of these reactions on judiciously chosen targets would provide an invaluable test of the complex theoretical models used to compute the large number of neutrino-nucleus cross sections that are needed., Comment: Review for special section on Neutrino Physics at Spallation Neutron Sources, 24 pages, 8 figures, uses iopart.cls

Proton capture by 17F plays an important role in the synthesis of nuclei in nova explosions. A revised rate for this reaction, based on a measurement of the 1H(17F,p)17F excitation function using a radioactive 17F beam at ORNL's Holifield Radioactive Ion Beam Facility, is used to calculate the nucleosynthesis in nova outbursts on the surfaces of 1.25 and 1.35 solar mass ONeMg white dwarfs and a 1.00 solar mass CO white dwarf. We find that the new 17F(p,gamma)18Ne reaction rate changes the abundances of some nuclides (e.g., 17O) synthesized in the hottest zones of an explosion on a 1.35 solar mass white dwarf by more than a factor of 10,000 compared to calculations using some previous estimates for this reaction rate, and by more than a factor of 3 when the entire exploding envelope is considered. In a 1.25 solar mass white dwarf nova explosion, this new rate changes the abundances of some nuclides synthesized in the hottest zones by more than a factor of 600, and by more than a factor of 2 when the entire exploding envelope is considered. Calculations for the 1.00 solar mass white dwarf nova show that this new rate changes the abundance of 18Ne by 21%, but has negligible effect on all other nuclides. Comparison of model predictions with observations is also discussed., Comment: 20 pages, 6 figures, accepted for publication in ApJ

Supernova simulations to date have assumed that during core collapse electron captures occur dominantly on free protons, while captures on heavy nuclei are Pauli-blocked and are ignored. We have calculated rates for electron capture on nuclei with mass numbers A=65-112 for the temperatures and densities appropriate for core collapse. We find that these rates are large enough so that, in contrast to previous assumptions, electron capture on nuclei dominates over capture on free protons. This leads to significant changes in core collapse simulations., Comment: 5 pages, 4 figures, REVTeX 4

General relativistic multi-group and multi-flavor Boltzmann neutrino transport in spherical symmetry adds a new level of detail to the numerical bridge between microscopic nuclear and weak interaction physics and the macroscopic evolution of the astrophysical object. Although no supernova explosions are obtained, we investigate the neutrino luminosities in various phases of the postbounce evolution for a wide range of progenitor stars between 13 and 40 solar masses. The signal probes the dynamics of material layered in and around the protoneutron star and is, within narrow limits, sensitive to improvements in the weak interaction physics. Only changes that dramatically exceed physical limitations allow experiments with exploding models. We discuss the differences in the neutrino signal and find the electron fraction in the innermost ejecta to exceed 0.5 as a consequence of thermal balance and weak equilibrium at the masscut., Comment: 8 pages, 4 figures. Proceedings of the Nuclear Physics in Astrophysics Conference, Debrecen, Hungary, 2002, to appear in Nuc. Phys. A. Color figures added and reference actualized

Based on general relativistic three-flavor Boltzmann neutrino transport, we compare simulations of the postbounce evolution for five different stellar progenitors with 13 and 20 solar masses (Nomoto & Hashimoto '88), and 15, 25, and 40 solar masses (Woosley & Weaver '95)., Comment: 4 pages, 3 figures, to appear in W. Hillebrandt and E. Mueller, eds., Proceedings of the 11th Workshop on "Nuclear Astrophysics" held at Rinberg Castle, February 11-16, 2002

With the exception of the Big Bang, responsible for 1,2H, 3,4He, and 7Li, stars act as sources for the composition of the interstellar medium. Cosmic rays are related to the latter and very probably due to acceleration of the mixed interstellar medium by shock waves from supernova remnants. Thus, the understanding of the abundance evolution in the interstellar medium and especially the enrichment of heavy elements, as a function of space and time, is essential. It reflects the history of star formation and the lifetimes of the diverse contributing stellar objects. Therefore, the understanding of the endpoints of stellar evolution is essential as well. These are mainly planetary nebulae and type II/Ib/Ic supernovae as evolutionary endpoints of single stars, but also events in binary systems can contribute, like e.g. supernovae of type Ia, novae and possibly X-ray bursts and neutron star or neutron star - black hole mergers. Despite many efforts, a full and self-consistent understanding of supernovae (the main contributors to nucleosynthesis in galaxies) is not existing, yet. Their fingerprints, however, seen either in spectra, lightcurves, radioactivities/decay gamma-rays or in galactic evolution, can help to constrain the composition of their ejecta and related model uncertainties., Comment: Review, 20 pages, uses style icrc.sty, to be published in Proc. ICRC 2001