Currently there is great interest in the outcomes and astrophysical implications of mergers of double degenerate binaries. In a commonly adopted approximation, the components of such binaries are represented by polytropes with an index n=3/2. We present detailed comparisons of stellar mass-transfer and merger simulations of polytropic binaries that have been carried out using two very different numerical algorithms --- a finite-volume "grid" code and a smoothed-particle hydrodynamics (SPH) code. We find that there is agreement in both the ultimate outcomes of the evolutions and the intermediate stages if the initial conditions for each code are chosen to match as closely as possible. We find that even with closely matching initial setups, the time it takes to reach a concordant evolution differs between the two codes because the initial depth of contact cannot be matched exactly. There is a general tendency for SPH to yield higher mass transfer rates and faster evolution to the final outcome. We also present comparisons of simulations calculated from two different energy equations: in one series we assume a polytropic equation of state and in the other series an ideal gas equation of state. In the latter series of simulations an atmosphere forms around the accretor which can exchange angular momentum and cause a more rapid loss of orbital angular momentum. In the simulations presented here, the effect of the ideal equation of state is to de-stabilize the binary in both SPH and grid simulations, but the effect is more pronounced in the grid code.

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.

We review existing SPH setup methods and outline their advantages, limitations and drawbacks. We present a new method for constructing initial conditions for smoothed particle hydrodynamics (SPH) simulations, which may also be of interest for N-body simulations, and demonstrate this method on a number of applications. This new method is inspired by adaptive binning techniques using weighted Voronoi tesselations. Particles are placed and iteratively moved based on their proximity to neighboring particles and the desired spatial resolution. This new method can satisfy arbitrarily complex spatial resolution requirements.

[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 use three-dimensional hydrodynamical simulations to study the rapid infall phase of the common envelope interaction of a red giant branch star of mass equal to 0.88 \msun and a companion star of mass ranging from 0.9 down to 0.1 \msun. We first compare the results obtained using two different numerical techniques with different resolutions, and find overall very good agreement. We then compare the outcomes of those simulations with observed systems thought to have gone through a common envelope. The simulations fail to reproduce those systems in the sense that most of the envelope of the donor remains bound at the end of the simulations and the final orbital separations between the donor's remnant and the companion, ranging from 26.8 down to 5.9 \rsun, are larger than the ones observed. We suggest that this discrepancy vouches for recombination playing an essential role in the ejection of the envelope and/or significant shrinkage of the orbit happening in the subsequent phase.

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.

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.

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 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.

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 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.

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 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.

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.

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.

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 present results from an extensive survey of 64 cavities in the X-ray halos of clusters, groups and normal elliptical galaxies. We show that the evolution of the size of the cavities as they rise in the X-ray atmosphere is inconsistent with the standard model of adiabatic expansion of purely hydrodynamic models. We also note that the majority of the observed bubbles should have already been shredded apart by Rayleigh-Taylor and Richtmyer-Meshkov instabilities if they were of purely hydrodynamic nature. Instead we find that the data agrees much better with a model where the cavities are magnetically dominated and inflated by a current-dominated magneto-hydrodynamic jet model, recently developed by Li et al. (2006) and Nakamura et al. (2006). We conduct complex Monte-Carlo simulations of the cavity detection process including incompleteness effects to reproduce the cavity sample's characteristics. We find that the current-dominated model agrees within 1sigma, whereas the other models can be excluded at >5sigma confidence. To bring hydrodynamic models into better agreement, cavities would have to be continuously inflated. However, these assessments are dependent on our correct understanding of the detectability of cavities in X-ray atmospheres, and will await confirmation when automated cavity detection tools become available in the future. Our results have considerable impact on the energy budget associated with active galactic nucleus feedback.

This is the third paper in a series analyzing X-ray emission from the hot interstellar medium in a sample of 54 normal elliptical galaxies observed by Chandra, focusing on 36 galaxies with sufficient signal to compute radial temperature profiles. We distinguish four qualitatively different types of profile: positive gradient (outwardly rising), negative gradients (falling), quasi-isothermal (flat) and hybrid (falling at small radii, then rising). We measure the mean logarithmic temperature gradients in two radial regions: from 0--2 $J$-band effective radii $R_J$ (excluding the central point source), and from 2--$4R_J$. We find the outer gradient to be uncorrelated with intrinsic host galaxy properties, but strongly influenced by the environment: galaxies in low-density environments tend to show negative outer gradients, while those in high-density environments show positive outer gradients, suggesting influence of circumgalactic hot gas. The inner temperature gradient is unaffected by the environment but strongly correlated with intrinsic host galaxy characteristics: negative inner gradients are more common for smaller, optically faint, low radio-luminosity galaxies, whereas positive gradients are found in bright galaxies with stronger radio sources. There is no evidence for bimodality in the distribution of inner or outer gradients. We propose three scenarios to explain the inner temperature gradients: (1) Weak AGN heat the ISM locally, higher-luminosity AGN heat the system globally through jets inflating cavities at larger radii; (2) The onset of negative inner gradients indicates a declining importance of AGN heating relative to other sources, such as compressional heating or supernovae; (3) The variety of temperature profiles are snapshots of different stages of a time-dependent flow.

We report on the discovery of a new quantitative relationship between X-ray gas morphology and radio and X-ray AGN luminosities in normal elliptical galaxies. This is the second paper in a series using data on 54 objects from the Chandra public archive and builds on the findings of Paper I, which demonstrated that hydrostatic equilibrium in elliptical galaxies holds, at best, only approximately and that the shape of the X-ray isophotes is unrelated to the shape of the gravitational potential. Instead, the gas is almost always asymmetrically disturbed. In this paper, we quantify the amount of asymmetry and study its correlation with other galaxy properties. We also determine radio powers and derive X-ray AGN luminosities for our galaxy sample. We find that the amount of asymmetry in the gas is correlated with both measures of AGN activity, in the sense that the hot gas is more disturbed in galaxies with higher radio and X-ray AGN luminosities. We find no evidence that galaxy density has significant effects on gas morphology. We do however find evidence for a correlation between gas asymmetry and the presence of hot ambient gas, which we interpret as a signature of hydrodynamic interactions with an external ambient medium. Surprisingly, the AGN-morphology connection persists all the way down to the weakest AGN luminosities in rather X-ray faint galaxies. This is strong morphological evidence that supports the general importance of AGN feedback, even in normal elliptical galaxies.

We review the current state of the art in double degenerate merger simulations to better understand the role this phenomenon plays in type Ia progenitors. Because the fate of a merged system may well depend on the exact evolution of the matter temperature (as well as mixing of the merged system), precision simulations are required to determine the true fate of these systems. Unfortunately, if we compare the results of current simulations, we find many-order of magnitude differences in quantities like mass-transfer rates in the merger process. We discuss these differences and outline an approach using verification and validation that should allow us to achieve a level of precision sufficient to determine the true fate (thermonuclear vs. collapse) of double degenerate mergers. Understanding the fate of lower-mass systems (e.g. R Coronae Borealis stars) may be key in our final testing phase.

We present preliminary results from Smooth Particle Hydrodynamics (SPH) simulations of common envelope evolution. We qualitatively compare the interaction between a 0.9M red giant with two different companion masses: a 0.05M brown dwarf and a 0.25M white dwarf companion.

We present results from magneto-hydrodynamic jet simulations in the magnetically dominated regime and find that they inflate cavities into the intracluster medium. Mock Chandra X-ray observations of our simulations show a strong resemblance to real cavities observed in X-ray images of clusters, both in terms of their morphology and thermodynamic structure. An analysis of the evolution of bubble sizes in the multi-cavity system Hydra A, as well as in a large sample of 64 cavities in 32 clusters shows that bubbles tend to expand much faster than expected in the purely adiabatic regime. Instead, we find that the bubbles follow more closely a trend predicted by our current-carrying jet models.

We argue the behavior of Poynting flux-dominated outflows from AGN in the galactic cluster systems by performing three-dimensional MHD simulations within the framework of the "magnetic tower" model. Of particular interests are the structure of MHD waves, the cylindrical radial force balance, the (de)collimation, and the stability properties of magnetic tower jets. Transition between the jet/lobe and the formation of wiggling jet by growing current-driven instability are discussed.

We present an X-ray analysis of 54 normal elliptical galaxies in the Chandra archive and isolate their hot gas component from the contaminating point source emission, allowing us to conduct, for the first time, a morphological analysis on the gas alone. A comparison with optical images and photometry shows that the hot gas morphology has surprisingly little in common with the shape of the stellar distribution. We observe no correlation between optical and X-ray ellipticities in the inner regions where stellar mass dominates over dark matter. A shallow correlation would be expected if the gas had settled into hydrostatic equilibrium with the gravitational potential. Instead, observed X-ray ellipticities exceed optical ellipticities in many cases. We exclude rotation as the dominant factor to produce the gas ellipticities. The gas appears disturbed, and hydrostatic equilibrium is the exception rather than the rule. Nearly all hydrostatic models can be ruled out at 99% confidence, based of their inability to reproduce the optical-X-ray correlation and large X-ray ellipticities. Hydrostatic models not excluded are those in which dark matter either dominates over stellar mass inside the inner half-light radius or has a prominently cigar-shaped distribution, both of which can be ruled out on other grounds. We conclude that, even for rather X-ray faint elliptical galaxies, the gas is at least so far out of equilibrium that it does not retain any information about the shape of the potential, and that X-ray derived radial mass profiles may be in error by factors of order unity.

We present a technique to adaptively bin sparse X-ray data using weighted Voronoi tesselations (WVTs). WVT binning is a generalisation of Cappellari & Copin's (2001) Voronoi binning algorithm, developed for integral field spectroscopy. WVT binning is applicable to many types of data and creates unbiased binning structures with compact bins that do not lead the eye. We apply the algorithm to simulated data, as well as several X-ray data sets, to create adaptively binned intensity images, hardness ratio maps and temperature maps with constant signal-to-noise ratio per bin. We also illustrate the separation of diffuse gas emission from contributions of unresolved point sources in elliptical galaxies. We compare the performance of WVT binning with other adaptive binning and adaptive smoothing techniques. We find that the CIAO tool csmooth creates serious artefacts and advise against its use to interpret diffuse X-ray emission.

We report on the discovery of a new correlation between global parameters of the hot interstellar gas in elliptical galaxies. We reanalyze archival Chandra data for 30 normal early-type systems, removing the contributions of resolved and unresolved point sources to reveal the X-ray morphology of the hot gas. We determine the half-light radius, R_X, and the mean surface brightness, I_X, from the gas surface brightness profiles. A spectral analysis determines the temperature, T_X, of the gas within 3 optical effective radii. We find that the galaxies lie on an X-ray Gas Fundamental Plane (XGFP) of the form T_X ~ R_X^{0.28} I_X^{0.22}. This is close to, but distinct from, a simple luminosity-temperature relation. The intrinsic width of the XGFP is only 0.07dex, nearly identical to that of the stellar (optical) fundamental plane (SFP). This is surprising since X-ray gas masses are typically ~10^{-2} of the stellar masses. We show that the XGFP is not a simple consequence of the virial theorem or hydrostatic equilibrium, and that it is essentially independent of the SFP. The XGFP thus represents a genuinely new constraint on the hydrodynamical evolution of elliptical galaxies.