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.
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 re
alism 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.
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.
We present a numerical investigation of the tidal disruption of white dwarfs by moderately massive black holes, with particular reference to the centers of dwarf galaxies and globular clusters. Special attention is given to the fate of white dwarfs o
f all masses that approach the black hole close enough to be disrupted and severely compressed to such extent that explosive nuclear burning can be triggered. Consistent modeling of the gas dynamics together with the nuclear reactions allows for a realistic determination of the explosive energy release. In the most favorable cases, the nuclear energy release may be comparable to that of typical type Ia supernovae. Although the explosion will increase the mass fraction escaping on hyperbolic orbits, a good fraction of the debris remains to be swallowed by the hole, causing a bright soft X-ray flare lasting for about a year. Such transient signatures, if detected, would be a compelling testimony for the presence of a moderately mass black hole (below $10^5 M_odot$).
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.
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 captu
re 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.
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.
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 th
e 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.
