In Pop III stellar models convection-induced mixing between H- and He-rich burning layers can induce a burst of nuclear energy and thereby substantially alter the subsequent evolution and nucleosynthesis in the first massive stars. We investigate H-He shell and core interactions in 26 stellar evolution simulations with masses $15 - 140 \,\mathrm{M}_{\odot}$, using five sets of mixing assumptions. In 22 cases, H-He interactions induce local nuclear energy release in the range $ \sim 10^{9} - 10^{13.5}\, \mathrm{L}_{\odot}$. The luminosities on the upper end of this range amount to a substantial fraction of the layer's internal energy over a convective advection timescale, indicating a dynamic stellar response that would violate 1D stellar evolution modelling assumptions. We distinguish four types of H-He interactions depending on the evolutionary phase and convective stability of the He-rich material. H-burning conditions during H-He interactions give $^{12}\mathrm{C}/^{13}\mathrm{C}$ ratios between $\approx 1.5$ to $\sim1000$ and [C/N] ratios from $\approx -1.7 $ to $\approx 3.2$ with a correlation that agrees well with observations of CEMP-no stars. We also explore Ca production from hot CNO breakout and find the simulations presented here likely cannot explain the observed Ca abundance in the most Ca-poor CEMP-no star. We describe the evolution leading to H-He interactions, which occur during or shortly after core-contraction phases. Three simulations without a H-He interaction are computed to Fe-core infall and a $140\,\mathrm{M}_{\odot}$ simulation becomes pair-unstable. We also discuss present modelling limitations and the need for 3D hydrodynamic models to fully understand these stellar evolutionary phases.

We present two mixing models for post-processing of 3D hydrodynamic simulations applied to convective-reactive i-process nucleosynthesis in a rapidly accreting white dwarf (RAWD) with [Fe/H] = -2.6, in which H is ingested into a convective He shell during a He flash. A 1D advective two-stream model is formulated with physically motivated radial and horizontal mixing coefficients constrained by 3D hydrodynamic simulations. A more traditional approach uses diffusion coefficients calculated from the same simulations. All 3D simulations include the energy feedback of the $^{12}$C(p,$\gamma)^{13}$N reaction from the entrainment of stably stratified H. Global oscillations of shell H ingestion in two of the RAWD simulations cause bursts of entrainment of H and energy feedback into the flow. With the same nuclear network as in the 3D simulations, the 1D advective two-stream model reproduces the rate and location of the H burning within the He shell closely matching the 3D simulation predictions, as well as qualitatively displaying the asymmetry of the $X_{\mathrm{H}}$ profiles between the up- and downstream. With a full i-process network the advective mixing model captures the difference in the n-capture nucleosynthesis in the up- and downstream. For example, $^{89}$Kr and $^{90}$Kr with half-lives of 3.18 min and 32.3 s differ by a factor 2 - 10 in the two streams. In this particular application the diffusion approach provides globally the same abundance distribution as the advective two-stream mixing model. The resulting i-process yields from the diffusive and advective post-processing models are compared with observations of the exemplary CEMP-r/s star CS31062-050.

Several anomalous elemental abundance ratios have been observed in the metal-poor star HD94028. We assume that its high [As/Ge] ratio is a product of a weak intermediate (i) neutron-capture process. Given that observational errors are usually smaller than predicted nuclear physics uncertainties, we have first set up a benchmark one-zone i-process nucleosynthesis simulation results of which provide the best fit to the observed abundances. We have then performed Monte Carlo simulations in which 113 relevant (n,$\gamma$) reaction rates of unstable species were randomly varied within Hauser-Feshbach model uncertainty ranges for each reaction to estimate the impact on the predicted stellar abundances. One of the interesting results of these simulations is a double-peaked distribution of the As abundance, which is caused by the variation of the $^{75}$Ga (n,$\gamma$) cross section. This variation strongly anti-correlates with the predicted As abundance, confirming the necessity for improved theoretical or experimental bounds on this cross section. The $^{66}$Ni (n,$\gamma$) reaction is found to behave as a major bottleneck for the i-process nucleosynthesis. Our analysis finds the Pearson product-moment correlation coefficient $r_\mathrm{P} > 0.2$ for all of the i-process elements with $32 \leq Z \leq 42$, with significant changes in their predicted abundances showing up when the rate of this reaction is reduced to its theoretically constrained lower bound. Our results are applicable to any other stellar nucleosynthesis site with the similar i-process conditions, such as Sakurai's object (V4334 Sagittarii) or rapidly-accreting white dwarfs.

Modeling the evolution of the elements in the Milky Way is a multidisciplinary and challenging task. In addition to simulating the 13 billion years evolution of our Galaxy, chemical evolution simulations must keep track of the elements synthesized and ejected from every astrophysical site of interest (e.g., supernova, compact binary merger). The elemental abundances of such ejecta, which are a fundamental input for chemical evolution codes, are usually taken from theoretical nucleosynthesis calculations performed by the nuclear astrophysics community. Therefore, almost all chemical evolution predictions rely on the nuclear physics behind those calculations. In this proceedings, we highlight the impact of nuclear physics uncertainties on galactic chemical evolution predictions. We demonstrate that nuclear physics and galactic evolution uncertainties both have a significant impact on interpreting the origin of neutron-capture elements in our Solar System. Those results serve as a motivation to create and maintain collaborations between the fields of nuclear astrophysics and galaxy evolution.

The origin of many elements of the periodic table remains an unsolved problem. While many nucleosynthetic channels are broadly understood, significant uncertainties remain regarding certain groups of elements such as the intermediate and rapid neutron-capture processes, the p-process, or the origin of odd-Z elements in the most metal-poor stars. Canada has a long tradition of leadership in nuclear astrophysics, dating back to the work of Alastair Cameron in the 1950s. Recent faculty hires have further boosted activity in the field, including transient observation and theory, survey science on galactic nucleosynthesis, and nuclear experiments. This white paper contains a brief overview of recent activity in the community, highlighting strengths in each sub-field, and provides recommendations to improve interdisciplinary collaboration. Sustaining Canadian leadership in the next decade will require, on the observational side, access to transient and non-transient surveys like LSST, SKA, or MSE, support for target-of-opportunity observing in current and future Canadian telescopes, and participation in next-generation X-ray telescopes such as ATHENA. State-of-the-art theoretical predictions will require an ambitious succession plan for the Niagara supercomputer to support large parallel jobs. We propose a funding instrument for postdoctoral training that reflects the interdisciplinary nature of nuclear astrophysics research, and the creation of a national collaborative funding program that allows for joint projects and workshop organization.

Interactions between convective shells in evolved massive stars have been linked to supernova impostors, to the production of the odd-Z elements Cl, K, and Sc, and they might also help generate the large-scale asphericities that are known to facilitate shock revival in supernova explosion models. We investigate the process of ingestion of C-shell material into a convective O-burning shell, including the hydrodynamic feedback from the nuclear burning of the ingested material. Our 3D hydrodynamic simulations span almost 3 dex in the total luminosity $L_\mathrm{tot}$. All but one of the simulations reach a quasi-stationary state with the entrainment rate and convective velocity proportional to $L_\mathrm{tot}$ and $L_\mathrm{tot}^{1/3}$, respectively. Carbon burning provides $14\,$--$\,33\%$ of the total luminosity, depending on the set of reactions considered. Equivalent simulations done on $768^3$ and $1152^3$ grids are in excellent quantitative agreement. The flow is dominated by a few large-scale convective cells. An instability leading to large-scale oscillations with Mach numbers in excess of $0.2$ develops in an experimental run with the energy yield from C burning increased by a factor of 10. This run represents most closely the conditions expected in a violent O-C shell merger, which is a potential production site for odd-Z elements such as K and Sc and which may seed asymmetries in the supernova progenitor. 1D simulations may underestimate the energy generation from the burning of ingested material by as much as a factor two owing to their missing the effect of clumpiness of entrained material on the nuclear reaction rate.

Aims. We investigate the s-process during the AGB phase of stellar models whose cores are enforced to rotate at rates consistent with asteroseismology observations of their progenitors and successors. Methods. We calculated new 2M$_{\odot}$, Z=0.01 models, rotating at 0, 125, and 250 km/s at the start of main sequence. An artificial, additional viscosity was added to enhance the transport of angular momentum in order to reduce the core rotation rates to be in agreement with asteroseismology observations. We compared rotation rates of our models with observed rotation rates during the MS up to the end of core He burning, and the white dwarf phase. Results. We present nucleosynthesis calculations for these rotating AGB models that were enforced to match the asteroseismic constraints on rotation rates of MS, RGB, He-burning, and WD stars. In particular, we calculated one model that matches the upper limit of observed rotation rates of core He-burning stars and we also included a model that rotates one order of magnitude faster than the upper limit of the observations. The s-process production in both of these models is comparable to that of non-rotating models. Conclusions. Slowing down the core rotation rate in stars to match the above mentioned asteroseismic constraints reduces the rotationally induced mixing processes to the point that they have no effect on the s-process nucleosynthesis. This result is independent of the initial rotation rate of the stellar evolution model. However, there are uncertainties remaining in the treatment of rotation in stellar evolution, which need to be reduced in order to confirm our conclusions, including the physical nature of our approach to reduce the core rotation rates of our models, and magnetic processes.

We have modelled the multicycle evolution of rapidly-accreting CO white dwarfs (RAWDs) with stable H burning intermittent with strong He-shell flashes on their surfaces for $0.7\leq M_\mathrm{RAWD}/M_\odot\leq 0.75$ and [Fe/H] ranging from $0$ to $-2.6$. We have also computed the i-process nucleosynthesis yields for these models. The i process occurs when convection driven by the He-shell flash ingests protons from the accreted H-rich surface layer, which results in maximum neutron densities $N_\mathrm{n,max}\approx 10^{13}$-$10^{15}\ \mathrm{cm}^{-3}$. The H-ingestion rate and the convective boundary mixing (CBM) parameter $f_\mathrm{top}$ adopted in the one-dimensional nucleosynthesis and stellar evolution models are constrained through 3D hydrodynamic simulations. The mass ingestion rate and, for the first time, the scaling laws for the CBM parameter $f_\mathrm{top}$ have been determined from 3D hydrodynamic simulations. We confirm our previous result that the high-metallicity RAWDs have a low mass retention efficiency ($\eta < 10\%$). A new result is that RAWDs with [Fe/H]$< -2$ have $\eta > 20\%$, therefore their masses may reach the Chandrasekhar limit and they may eventually explode as SNeIa. This result and the good fits of the i-process yields from the metal-poor RAWDs to the observed chemical composition of the CEMP-r/s stars suggest that some of the present-day CEMP-r/s stars could be former distant members of triple systems, orbiting close binary systems with RAWDs that may have later exploded as SNeIa.

We analyze the production of the element Cr in galactic chemical evolution (GCE) models using the NuGrid nucleosynthesis yields set. We show that the unusually large [Cr/Fe] abundance at [Fe/H] $\approx 0$ reported by previous studies using those yields and predicted by our Milky Way model originates from the merging of convective Si-burning and C-burning shells in a 20 $M_\odot$ model at metallicity $Z=0.01$, about an hour before the star explodes. This merger mixes the incomplete burning material in the Si shell, including $^{51}$V and $^{52}$Cr, out to the edge of the carbon/oxygen (CO) core. The adopted supernova model ejects the outer 2 $M_\odot$ of the CO core, which includes a significant fraction of the Cr-rich material. When including this 20 $M_\odot$ model at $Z=0.01$ in the yields interpolation scheme of our GCE model for stars in between 15 and 25 $M_\odot$, we overestimate [Cr/Fe] by an order of magnitude at [Fe/H] $\approx$ 0 relative to observations in the Galactic disk. This raises a number of questions regarding the occurrence of Si-C shell mergers in nature, the accuracy of different simulation approaches, and the impact of such mergers on the pre-supernova structure and explosion dynamics. According to the conditions in this 1D stellar model, the substantial penetration of C-shell material into the Si-shell could launch a convective-reactive global oscillation, if a merger does take place. In any case, GCE provides stringent constraints on the outcome of this stellar evolution phase.

The production of the neutron-capture isotopes beyond iron that we observe today in the solar system is the result of the combined contribution of the r-process, the s- process and possibly the i-process. Low-mass AGB (2 < M/Msun < 3) and massive (M >10 Msun ) stars have been identified as the sites of the s-process. In this work we consider the evolution and nucleosynthesis of low-mass AGB stars. We provide an update of the NuGrid Set models, adopting the same general physics assumptions but using an updated convective-boundary mixing model accounting for the contribution from internal gravity waves. The combined data set includes the initial masses Mzams/Msun = 2, 3 for Z = 0.03, 0.02, 0.01. These models are computed with the MESA stellar code and the evolution is followed up to the end of the AGB phase. The nucleosynthesis was calculated for all isotopes in post-processing with the NuGrid mppnp code. The convective boundary mixing model leads to the formation of a 13C-pocket three times wider compared to the one obtained in the previous set of models, bringing the simulation results now in closer agreement with observations. We also discuss the potential impact of other processes inducing mixing, like rotation, adopting parametric models compatible with theory and observations. Complete yield data tables, derived data products and online analytic data access are provided.

Gamma-ray astronomy explores the most energetic photons in nature to address some of the most pressing puzzles in contemporary astrophysics. It encompasses a wide range of objects and phenomena: stars, supernovae, novae, neutron stars, stellar-mass black holes, nucleosynthesis, the interstellar medium, cosmic rays and relativistic-particle acceleration, and the evolution of galaxies. MeV gamma-rays provide a unique probe of nuclear processes in astronomy, directly measuring radioactive decay, nuclear de-excitation, and positron annihilation. The substantial information carried by gamma-ray photons allows us to see deeper into these objects, the bulk of the power is often emitted at gamma-ray energies, and radioactivity provides a natural physical clock that adds unique information. New science will be driven by time-domain population studies at gamma-ray energies. This science is enabled by next-generation gamma-ray instruments with one to two orders of magnitude better sensitivity, larger sky coverage, and faster cadence than all previous gamma-ray instruments. This transformative capability permits: (a) the accurate identification of the gamma-ray emitting objects and correlations with observations taken at other wavelengths and with other messengers; (b) construction of new gamma-ray maps of the Milky Way and other nearby galaxies where extended regions are distinguished from point sources; and (c) considerable serendipitous science of scarce events -- nearby neutron star mergers, for example. Advances in technology push the performance of new gamma-ray instruments to address a wide set of astrophysical questions.

Massive-star binaries can undergo a phase where one of the two stars expands during its advanced evolutionary stage as a giant and envelops its companion, ejecting the hydrogen envelope and tightening its orbit. Such a common envelope phase is required to tighten the binary orbit in the formation of many of the observed X-ray binaries and merging compact binary systems. In the formation scenario for neutron star binaries, the system might pass through a phase where a neutron star spirals into the envelope of its giant star companion. These phases lead to mass accretion onto the neutron star. Accretion onto these common-envelope-phase neutron stars can eject matter that has undergone burning near to the neutron star surface. This paper presents nucleosynthetic yields of this ejected matter, using population synthesis models to study the importance of these nucleosynthetic yields in a galactic chemical evolution context. Depending on the extreme conditions in temperature and density found in the accreted material, both proton-rich and neutron-rich nucleosynthesis can be obtained, with efficient production of neutron rich isotopes of low Z material at the most extreme conditions, and proton rich isotopes, again at low Z, in lower density models. Final yields are found to be extremely sensitive to the physical modeling of the accretion phase. We show that neutron stars accreting in binary common envelopes might be a new relevant site for galactic chemical evolution, and therefore more comprehensive studies are needed to better constrain nucleosynthesis in these objects.

The special computational challenges of simulating 3-D hydrodynamics in deep stellar interiors are discussed, and numerical algorithmic responses described. Results of recent simulations carried out at scale on the NSF's Blue Waters machine at the University of Illinois are presented, with a special focus on the computational challenges they address. Prospects for future work using GPU-accelerated nodes such as those on the DoE's new Summit machine at Oak Ridge National Laboratory are described, with a focus on numerical algorithmic accommodations that we believe will be necessary.

We report on our ongoing investigation into the nucleosynthetic and hydrodynamic nature of mixing at the interface between the H- and He-convection zones in massive Pop III stars. Studying a grid of 26 1D stellar evolution simulations with different mixing assumptions, we find that H-He interactions occur in 23/26 cases. We demonstrate the nucleosynthesis expected in a H-He interaction in an 80M$_\odot$. Finally, we describe our progress in simulating a Pop III double convection zone in the PPMStar hydrodynamics code.

The stellar yields for galactic modeling applications (SYGMA) code is an open-source module that models the chemical ejecta and feedback of simple stellar populations (SSPs). It is intended for use in hydrodynamical simulations and semi-analytic models of galactic chemical evolution. The module includes the enrichment from asymptotic giant branch (AGB) stars, massive stars, SNIa and neutron-star mergers. An extensive and extendable stellar yields library includes the NuGrid yields with all elements and many isotopes up to Bi. Stellar feedback from mechanic and frequency-dependent radiative luminosities are computed based on NuGrid stellar models and their synthetic spectra. The module further allows for customizable initial-mass functions and supernova Ia (SNIa) delay-time distributions to calculate time-dependent ejecta based on stellar yield input. A variety of r-process sites can be included. A comparison of SSP ejecta based on NuGrid yields with those from Portinari et al. (1998) and Marigo (2001) reveals up to a factor of 3.5 and 4.8 less C and N enrichment from AGB stars at low metallicity, a result we attribute to NuGrid's modeling of hot-bottom burning. Different core-collapse supernova explosion and fallback prescriptions may lead to substantial variations for the accumulated ejecta of C, O and Si in the first $10^7\, \mathrm{yr}$ at $Z=0.001$. An online interface of the open-source SYGMA module enables interactive simulations, analysis and data extraction of the evolution of all species formed by the evolution of simple stellar populations.

We provide here a significant extension of the NuGrid Set 1 models in mass coverage and toward lower metallicity, adopting the same physics assumptions. The combined data set now includes the initial masses M/Msun = 1, 1.65, 2, 3, 4, 5, 6, 7, 12, 15, 20, 25 for Z = 0.02, 0.01, 0.006, 0.001, 0.0001 with alpha-enhanced composition for the lowest three metallicities. These models are computed with the MESA stellar evolution code and are evolved up to the AGB, the white dwarf stage, or until core collapse. The nucleosynthesis was calculated for all isotopes in post-processing with the NuGrid mppnp code. Explosive nucleosynthesis is based on semi-analytic 1D shock models. Metallicity-dependent mass loss, convective boundary mixing in low- and intermediate mass models and H and He core burning massive star models is included. Convective O-C shell mergers in some stellar models lead to the strong production of odd-Z elements P, Cl, K and Sc. In AGB models with hot dredge-up the convective boundary mixing efficiency is reduced to accommodate for its energetic feedback. In both low-mass and massive star models at the lowest metallicity H-ingestion events are observed and lead to i-process nucleosynthesis and substantial N-15 production. Complete yield data tables, derived data products and online analytic data access are provided.

The first peak s-process elements Rb, Sr, Y and Zr in the post-AGB star Sakurai's object (V4334 Sagittarii) have been proposed to be the result of i-process nucleosynthesis in a post-AGB very-late thermal pulse event. We estimate the nuclear physics uncertainties in the i-process model predictions to determine whether the remaining discrepancies with observations are significant and point to potential issues with the underlying astrophysical model. We find that the dominant source in the nuclear physics uncertainties are predictions of neutron capture rates on unstable neutron rich nuclei, which can have uncertainties of more than a factor 20 in the band of the i-process. We use a Monte Carlo variation of 52 neutron capture rates and a 1D multi-zone post-processing model for the i-process in Sakurai's object to determine the cumulative effect of these uncertainties on the final elemental abundance predictions. We find that the nuclear physics uncertainties are large and comparable to observational errors. Within these uncertainties the model predictions are consistent with observations. A correlation analysis of the results of our MC simulations reveals that the strongest impact on the predicted abundances of Rb, Sr, Y and Zr is made by the uncertainties in the (n,$\gamma$) reaction rates of $^{85}$Br, $^{86}$Br, $^{87}$Kr, $^{88}$Kr, $^{89}$Kr, $^{89}$Rb, $^{89}$Sr, and $^{92}$Sr. (abridged)

Collaborations in astronomy and astrophysics are faced with numerous cyber infrastructure challenges, such as large data sets, the need to combine heterogeneous data sets, and the challenge to effectively collaborate on those large, heterogeneous data sets with significant processing requirements and complex science software tools. The cyberhubs system is an easy-to-deploy package for small to medium-sized collaborations based on the Jupyter and Docker technology, that allows web-browser enabled, remote, interactive analytic access to shared data. It offers an initial step to address these challenges. The features and deployment steps of the system are described, as well as the requirements collection through an account of the different approaches to data structuring, handling and available analytic tools for the NuGrid and PPMstar collaborations. NuGrid is an international collaboration that creates stellar evolution and explosion physics and nucleosynthesis simulation data. The PPMstar collaboration performs large-scale 3D stellar hydrodynamics simulation of interior convection in the late phases of stellar evolution. Examples of science that is presently performed on cyberhubs, in the areas 3D stellar hydrodynamic simulations, stellar evolution and nucleosynthesis and Galactic chemical evolution, are presented.

Rapidly accreting white dwarfs (RAWDs) have been proposed as contributors to the chemical evolution of heavy elements in the Galaxy. Here, we test this scenario for the first time and determine the contribution of RAWDs to the solar composition of first-peak neutron-capture elements. We add the metallicity-dependent contribution of RAWDs to the one-zone galactic chemical evolution code OMEGA according to RAWD rates from binary stellar population models combined with metallicity-dependent $i$-process stellar yields calculated following the models of Denissenkov et al. (2017). With this approach we find that the contribution of RAWDs to the evolution of heavy elements in the Galaxy could be responsible for a significant fraction of the solar composition of Kr, Rb, Sr, Y, Zr, Nb, and Mo ranging from $2$ to $45\%$ depending on the element, the enrichment history of the Galactic gas, and the total mass ejected per RAWD. This contribution could explain the missing solar Lighter Element Primary Process for some elements (e.g., Sr, Y, and Zr). We do not overproduce any isotope relative to the solar composition, but $^{96}$Zr is produced in a similar amount. The $i$ process produces efficiently the Mo stable isotopes $^{95}$Mo and $^{97}$Mo. When nuclear reaction rate uncertainties are combined with our GCE uncertainties, the upper limits for the predicted RAWD contribution increase by a factor of $1.5-2$ for Rb, Sr, Y, and Zr, and by 3.8 and 2.4 for Nb and Mo, respectively. We discuss the implication of the RAWD stellar evolution properties on the single degenerate Type Ia supernova scenario.

Convective boundary mixing (CBM) in the advanced evolutionary stages of massive stars is not well understood. Structural changes caused by convection have an impact on the evolution as well as the subsequent supernova, or lack thereof. The effects of convectively driven mixing across convective boundaries during the post He core burning evolution of $25\mathrm{M}_{\odot}$, solar-metallicity, non-rotating stellar models is studied using the MESA stellar evolution code. CBM is modelled using the exponentially decaying diffusion coefficient equation, the free parameter of which, $f_{\mathrm{cbm}}$, is varied systematically throughout the course of the stellar model's evolution with values of $(0.002, 0.012, 0.022, 0.032)$. The effect of varying this parameter produces mass ranges at collapse in the ONe, Si, Fe cores of $(1.82\mathrm{M}_{\odot}, 4.36\mathrm{M}_{\odot})$, $(1.67\mathrm{M}_{\odot}, 1.99\mathrm{M}_{\odot})$ and $(1.46\mathrm{M}_{\odot}, 1.70\mathrm{M}_{\odot})$ respectively, with percent differences from the model with minimal CBM as large as 86.3%. At the presupernova stage, the compactness of the stellar cores from O'Connor & Ott (2011), $\xi_M$, exhibit a range of $(0.120, 0.354)$, suggesting that the extent of CBM in the advanced burning stages of massive stars is an important consideration for the explodability and type of compact remnant. The nucleosynthetic yields from the models, most notably C, O, Ne, Mg and Si are also significantly affected by the CBM assumptions, showing non-linear trends with increased mixing. The simulations show that interactions between convective C, Ne and O shells produce significant non-linear changes in the evolution, whereas from the end of Si burning, the structural changes attributed to the CBM are dominated by the growth of the convective C shell. Progenitor structures for all the models are available from HERE (link and DOI to appear).

We have investigated a highly energetic H-ingestion event during shell He burning leading to H-burning luminosities of $\log(L_\mathrm{H}/\mathrm{\Lsun}) \sim 13$ in a $45 \Msun$ Pop III massive stellar model. In order to track the nucleosynthesis which may occur in such an event, we run a series of single-zone nucleosynthesis models for typical conditions found in the stellar evolution model. Such nucleosynthesis conditions may lead to $\textit{i}-$process neutron densities of up to $ \sim 10^{13} \mathrm{cm}^{-3}$. The resulting simulation abundance pattern, where Mg comes from He burning and Ca from the $\textit{i}$ process, agrees with the general observed pattern of the most iron-poor star currently known, SMSS J031300.36-670839.3. However, Na is also efficiently produced in these $i$ process conditions, and the prediction exceeds observations by $\sim 2.5 \unitstyle{dex}$. While this probably rules out this model for SMSS J031300.36-670839.3, the typical $\textit{i}-$process signature of combined He burning and $\textit{i}$ process of higher than solar [Na/Mg], [Mg/Al] and low [Ca/Mg] is reproducing abundance features of the two next most iron-poor stars HE\,1017-5240 and HE\,1327-2326 very well. The $i$ process does not reach Fe which would have to come from a low level of additional enrichment. $\textit{i}$ process in hyper-metal poor or Pop III massive stars may be able to explain certain abundance patterns observed in some of the most-metal poor CEMP-no stars.

We address the deficiency of odd-Z elements P, Cl, K and Sc in galactic chemical evolution models through an investigation of the nucleosynthesis of interacting convective O- and C shells in massive stars. 3D hydrodynamic simulations of O-shell convection with moderate C-ingestion rates show no dramatic deviation from spherical symmetry. We derive a spherically averaged diffusion coefficient for 1D nucleosynthesis simulations which show that such convective-reactive ingestion events can be a production site for P, Cl, K and Sc. An entrainment rate of $10^{-3}{{\rm M}_\odot}$/s features overproduction factors $OP_\mathrm{s} \approx 7$. Full O-C shell mergers in our 1D stellar evolution massive star models have overproduction factors $OP_\mathrm{m}>1 \mathrm{dex}$ but for such cases 3D hydrodynamic simulations suggest deviations from spherical symmetry. p-process species can be produced with overproduction factors of $OP_\mathrm{m}>1 \mathrm{dex}$, e.g. for $^{130,132}Ba$. Using the uncertain prediction of the $15{{\rm M}_\odot}$, $Z=0.02$ massive star model ($OP_\mathrm{m} \approx 15$) as representative for merger or entrainment convective-reactive events involving O- and C-burning shells, and assume that such events occur in more than 50% of all stars, our chemical evolution models reproduce the observed Galactic trends of the odd-Z elements.

Atmospheric parameters and chemical compositions for ten stars with metallicities in the region of -2.2< [Fe/H] <-0.6 were precisely determined using high resolution, high signal to noise, spectra. For each star the abundances, for 14 to 27 elements, were derived using both LTE and NLTE approaches. In particular, differences by assuming LTE or NLTE are about 0.10 dex; depending on [Fe/H], Teff, gravity and element lines used in the analysis. We find that the O abundance has the largest error, ranging from 0.10 and 0.2 dex. The best measured elements are Cr, Fe, and Mn; with errors etween 0.03 and 0.11 dex. The stars in our sample were included in previous different observational work. We provide a consistent data analysis. The data dispersion introduced in the literature by different techniques and assumptions used by the different authors is within the observational errors, excepting for HD103095. We compare these results with stellar observations from different data sets and a number of theoretical galactic chemical evolution (GCE) simulations. We find a large scatter in the GCE results, used to study the origin of the elements. Within this scatter as found in previous GCE simulations, we cannot reproduce the evolution of the elemental ratios [Sc/Fe], [Ti/Fe], and [V/Fe] at different metallicities. The stellar yields from core collapse supernovae (CCSN) are likely primarily responsible for this discrepancy. Possible solutions and open problems are discussed.

V4334 Sgr (a.k.a. Sakurai's object) is the central star of an old planetary nebula that underwent a very late thermal pulse a few years before its discovery in 1996. We have been monitoring the evolution of the optical emission line spectrum since 2001. The goal is to improve the evolutionary models by constraining them with the temporal evolution of the central star temperature. In addition the high resolution spectral observations obtained by X-shooter and ALMA show the temporal evolution of the different morphological components.

Based on stellar evolution simulations, we demonstrate that rapidly accreting white dwarfs in close binary systems are an astrophysical site for the intermediate neutron-capture process. During recurrent and very strong He-shell flashes in the stable H-burning accretion regime H-rich material enters the He-shell flash convection zone. $^{12}$C(p,$\gamma)^{13}$N reactions release enough energy to potentially impact convection, and i process is activated through the $^{13}$C($\alpha$,n)$^{16}$O reaction. The H-ingestion flash may not cause a split of the convection zone as it was seen in simulations of He-shell flashes in post-AGB and low-Z AGB stars. We estimate that for the production of first-peak heavy elements this site can be of similar importance for galactic chemical evolution as the s-process production by low-mass AGB stars. The He-shell flashes result in the expansion and, ultimately, ejection of the accreted and then i-process enriched material, via super-Eddington luminosity winds or Roche-lobe overflow. The white dwarf models do not retain any significant amount of the accreted mass, with a He retention efficiency of $\leq 10\%$ depending on mass and convective boundary mixing assumptions. This makes the evolutionary path of such systems to supernova Ia explosion highly unlikely.

We use the OMEGA galactic chemical evolution code to investigate how the assumptions used for the treatment of galactic inflows and outflows impact numerical predictions. The goal is to determine how our capacity to reproduce the chemical evolution trends of a galaxy is affected by the choice of implementation used to include those physical processes. In pursuit of this goal, we experiment with three different prescriptions for galactic inflows and outflows and use OMEGA within a Markov Chain Monte Carlo code to recover the set of input parameters that best reproduces the chemical evolution of nine elements in the dwarf spheroidal galaxy Sculptor. This provides a consistent framework for comparing the best-fit solutions generated by our different models. Despite their different degrees of intended physical realism, we found that all three prescriptions can reproduce in an almost identical way the stellar abundance trends observed in Sculptor. While the three models have the same capacity to fit the data, the best values recovered for the parameters controlling the number of Type Ia supernovae and the strength of galactic outflows, are substantially different and in fact mutually exclusive from one model to another. For the purpose of understanding how a galaxy evolves, we conclude that only reproducing the evolution of a limited number of elements is insufficient and can lead to misleading conclusions. More elements or additional constraints such as the galaxy's star formation efficiency and the gas fraction are needed in order to break the degeneracy between the different modeling assumptions. Our results show that the successes and failures of chemical evolution models are predominantly driven by the input stellar yields, rather than by the complexity of the galaxy model itself. Simple models such as OMEGA are therefore sufficient to test and validate stellar yields.

This work investigates the properties of convection in stars with particular emphasis on entrainment across the upper convective boundary (CB). Idealised simulations of turbulent convection in the O-burning shell of a massive star are performed in $4\pi$ geometry on $768^3$ and $1536^3$ grids, driven by a representative heating rate. A heating series is also performed on the $768^3$ grid. The $1536^3$ simulation exhibits an entrainment rate at the upper CB of $1.33\times10^{-6}~M_\odot~\mathrm{s}^{-1}$. The $768^3$ simulation with the same heating rate agrees within 17 per cent. The entrainment rate at the upper convective boundary is found to scale linearly with the driving luminosity and with the cube of the shear velocity at the upper boundary, while the radial RMS fluid velocity scales with the cube root of the driving luminosity, as expected. The mixing is analysed in a 1D diffusion framework, resulting in a simple model for CB mixing. The analysis confirms previous findings that limiting the MLT mixing length to the distance to the CB in 1D simulations better represents the spherically-averaged radial velocity profiles from the 3D simulations and provides an improved determination of the reference diffusion coefficient $D_0$ for the exponential diffusion CB mixing model in 1D. From the 3D simulation data we adopt as the convective boundary the location of the maximum gradient in the horizontal velocity component which has $2\sigma$ spatial fluctuations of $\approx0.17 H_P$ . The exponentially decaying diffusion CB mixing model with $f = 0.03$ reproduces the spherically-averaged 3D abundance profiles.

Galactic chemical evolution is a topic that involves nuclear physics, stellar evolution, galaxy evolution, observation, and cosmology. Continuous communication and feedback between these fields is a key component in improving our understanding of how galaxies form and how elements are created and recycled in galaxies and intergalactic space. In this proceedings, we present the current state of the JINA-NuGrid chemical evolution pipeline. It is designed to probe the impact of nuclear astrophysics uncertainties on galactic chemical evolution, to improve our knowledges regarding the origin of the elements in a cosmological context, and to create the required interdisciplinary connections.

We present hydrodynamic simulations of the common envelope binary interaction between a giant star and a compact companion carried out with the adaptive mesh refinement code ENZO and the smooth particle hydrodynamics code PHANTOM. These simulations mimic the parameters of one of the simulations by Passy et al., but assess the impact of a larger, more realistic initial orbital separation on the simulation outcome. We conclude that for both codes the post-common envelope separation is somewhat larger and the amount of unbound mass slightly greater when the initial separation is wide enough that the giant does not yet overflow or just overflows its Roche lobe. PHANTOM has been adapted to the common envelope problem here for the first time and a full comparison with ENZO is presented, including an investigation of convergence as well as energy and angular momentum conservation. We also set our simulations in the context of past simulations. This comparison reveals that it is the expansion of the giant before rapid in-spiral and not spinning up of the star that causes a larger final separation. We also suggest that the large range in unbound mass for different simulations is difficult to explain and may have something to do with simulations that are not fully converged.

We used a one-zone chemical evolution model to address the question of how many masses and metallicities are required in grids of massive stellar models in order to ensure reliable galactic chemical evolution predictions. We used a set of yields that includes seven masses between 13 and 30 Msun, 15 metallicities between 0 and 0.03 in mass fraction, and two different remnant mass prescriptions. We ran several simulations where we sampled subsets of stellar models to explore the impact of different grid resolutions. Stellar yields from low- and intermediate-mass stars and from Type Ia supernovae have been included in our simulations, but with a fixed grid resolution. We compared our results with the stellar abundances observed in the Milky Way for O, Na, Mg, Si, Ca, Ti, and Mn. Our results suggest that the range of metallicity considered is more important than the number of metallicities within that range, which only affects our numerical predictions by about 0.1 dex. We found that our predictions at [Fe/H] < -2 are very sensitive to the metallicity range and the mass sampling used for the lowest metallicity included in the set of yields. Variations between results can be as high as 0.8 dex, for any remnant mass prescription. At higher [Fe/H], we found that the required number of masses depends on the element of interest and on the remnant mass prescription. With a monotonic remnant mass prescription where every model explodes as a core-collapse supernova, the mass resolution induces variations of 0.2 dex on average. But with a remnant mass prescription that includes islands of non-explodability, the mass resolution can cause variations of about 0.2 to 0.7 dex depending on the choice of metallicity range. With such a prescription, explosive or non-explosive models can be missed if not enough masses are selected, resulting in over- or under-estimations of the mass ejected by massive stars.

MSE is an 11.25m aperture observatory with a 1.5 square degree field of view that will be fully dedicated to multi-object spectroscopy. More than 3200 fibres will feed spectrographs operating at low (R ~ 2000 - 3500) and moderate (R ~ 6000) spectral resolution, and approximately 1000 fibers will feed spectrographs operating at high (R ~ 40000) resolution. MSE is designed to enable transformational science in areas as diverse as tomographic mapping of the interstellar and intergalactic media; the in-situ chemical tagging of thick disk and halo stars; connecting galaxies to their large scale structure; measuring the mass functions of cold dark matter sub-halos in galaxy and cluster-scale hosts; reverberation mapping of supermassive black holes in quasars; next generation cosmological surveys using redshift space distortions and peculiar velocities. MSE is an essential follow-up facility to current and next generations of multi-wavelength imaging surveys, including LSST, Gaia, Euclid, WFIRST, PLATO, and the SKA, and is designed to complement and go beyond the science goals of other planned and current spectroscopic capabilities like VISTA/4MOST, WHT/WEAVE, AAT/HERMES and Subaru/PFS. It is an ideal feeder facility for E-ELT, TMT and GMT, and provides the missing link between wide field imaging and small field precision astronomy. MSE is optimized for high throughput, high signal-to-noise observations of the faintest sources in the Universe with high quality calibration and stability being ensured through the dedicated operational mode of the observatory. (abridged)

The $s$-process nucleosynthesis in Asymptotic Giant Branch (AGB) stars depends on the modeling of convective boundaries. We present models and s-process simulations that adopt a treatment of convective boundaries based on the results of hydrodynamic simulations and on the theory of mixing due to gravity waves in the vicinity of convective boundaries. Hydrodynamics simulations suggest the presence of convective boundary mixing (CBM) at the bottom of the thermal pulse-driven convective zone. Similarly, convection-induced mixing processes are proposed for the mixing below the convective envelope during third dredge-up where the 13C pocket for the s process in AGB stars forms. In this work we apply a CBM model motivated by simulations and theory to models with initial mass $M = 2$ and $M = 3M_\odot$, and with initial metal content Z = 0.01 and Z = 0.02. As reported previously, the He-intershell abundance of 12C and 16O are increased by CBM at the bottom of pulse-driven convection zone. This mixing is affecting the $^{22}Ne(\alpha,n)^{25}Mg$ activation and the s-process effciency in the 13C-pocket. In our model CBM at the bottom of the convective envelope during the third dredgeup represents gravity wave mixing. We take further into account that hydrodynamic simulations indicate a declining mixing efficiency already about a pressure scale height from the convective boundaries, compared to mixing-length theory. We obtain the formation of the 13C-pocket with a mass of $\approx 10^{-4}M_\odot$. The final $s$-process abundances are characterized by 0.36 < [s=Fe] < 0.78 and the heavy-to-light s-process ratio is 0.23 < [hs=ls] < 0.45. Finally, we compare our results with stellar observations, pre-solar grain measurements and previous work.

We use a simple one-zone galactic chemical evolution model to quantify the uncertainties generated by the input parameters in numerical predictions, for a galaxy with properties similar to those of the Milky Way. We compiled several studies from the literature to gather the current constraints for our simulations regarding the typical value and uncertainty of seven basic parameters, which are: the lower and upper mass limits of the stellar initial mass function (IMF), the slope of the high-mass end of the stellar IMF, the slope of the delay-time distribution function of Type Ia supernovae (SNe Ia), the number of SNeIa per Msun formed, the total stellar mass formed, and the final mass of gas. We derived a probability distribution function to express the range of likely values for every parameter, which were then included in a Monte Carlo code to run several hundred simulations with randomly selected input parameters. This approach enables us to analyze the predicted chemical evolution of 16 elements in a statistical way by identifying the most probable solutions along with their 68% and 95% confidence levels. Our results show that the overall uncertainties are shaped by several input parameters that individually contribute at different metallicities, and thus at different galactic ages. The level of uncertainty then depends on the metallicity and is different from one element to another. Among the seven input parameters considered in this work, the slope of the IMF and the number of SNe Ia are currently the two main sources of uncertainty.

We provide a set of stellar evolution and nucleosynthesis calculations that applies established physics assumptions simultaneously to low- and intermediate-mass and massive star models. Our goal is to provide an internally consistent and comprehensive nuclear production and yield data base for applications in areas such as pre-solar grain studies. Our non-rotating models assume convective boundary mixing where it has been adopted before. We include 8 (12) initial masses for $Z = 0.01$ ($0.02$). Models are followed either until the end of the asymptotic giant branch phase or the end of Si burning, complemented by a simple analytic core-collapse supernova models with two options for fallback and shock velocities. The explosions show which pre-supernova yields will most strongly be effected by the explosive nucleosynthesis. We discuss how these two explosion parameters impacts the light elements and the $s$ and $p$ process. For low- and intermediate-mass models our stellar yields from H to Bi include the effect of convective boundary mixing at the He-intershell boundaries and the stellar evolution feedback of the mixing process that produces the $^{13}$C pocket. All post-processing nucleosynthesis calculations use the same nuclear reaction rate network and nuclear physics input. We provide a discussion of the nuclear production across the entire mass range organized by element group. All our stellar nucleosynthesis profile and time evolution output is available electronically, and tools to explore the data on the NuGrid VOspace hosted by the Canadian Astronomical Data Centre are introduced.

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

We present a detailed analysis of the composition and nucleosynthetic origins of the heavy elements in the metal-poor ([Fe/H]=-1.62+/-0.09) star HD94028. Previous studies revealed that this star is mildly enhanced in elements produced by the slow neutron-capture process (s-process; e.g., [Pb/Fe]=+0.79+/-0.32) and rapid neutron-capture process (r-process; e.g., [Eu/Fe]=+0.22+/-0.12), including unusually large molybdenum ([Mo/Fe]=+0.97+/-0.16) and ruthenium ([Ru/Fe]=+0.69+/-0.17) enhancements. However, this star is not enhanced in carbon ([C/Fe]=-0.06+/-0.19). We analyze an archival near-ultraviolet spectrum of HD94028, collected using the Space Telescope Imaging Spectrograph on board the Hubble Space Telescope, and other archival optical spectra collected from ground-based telescopes. We report abundances or upper limits derived from 64 species of 56 elements. We compare these observations with s-process yields from low-metallicity AGB evolution and nucleosynthesis models. No combination of s- and r-process patterns can adequately reproduce the observed abundances, including the super-solar [As/Ge] ratio (+0.99+/-0.23) and the enhanced [Mo/Fe] and [Ru/Fe] ratios. We can fit these features when including an additional contribution from the intermediate neutron-capture process (i process), which perhaps operated by the ingestion of H in He-burning convective regions in massive stars, super-AGB stars, or low-mass AGB stars. Currently, only the i process appears capable of consistently producing the super-solar [As/Ge] ratios and ratios among neighboring heavy elements found in HD94028. Other metal-poor stars also show enhanced [As/Ge] ratios, hinting that operation of the i process may have been common in the early Galaxy.

Motivated by recent results in stellar evolution that predict the existence of hybrid white dwarf (WD) stars with a C-O core inside an O-Ne shell, we simulate thermonuclear (Type Ia) supernovae from these hybrid progenitors. We use the FLASH code to perform multidimensional simulations in the deflagration to detonation transition (DDT) explosion paradigm. Our hybrid progenitor models were produced with the MESA stellar evolution code and include the effects of the Urca process, and we map the progenitor model to the FLASH grid. We performed a suite of DDT simulations over a range of ignition conditions consistent with the progenitor's thermal and convective structure assuming multiple ignition points. To compare the results from these hybrid WD stars to previous results from C-O white dwarfs, we construct a set of C-O WD models with similar properties and similarly simulate a suite of explosions. We find that despite significant variability within each suite, trends distinguishing the explosions are apparent in their $^{56}$Ni yields and the kinetic properties of the ejecta. We comment on the feasibility of these explosions as the source of some classes of observed subluminous events.

We investigate the evolution of super-AGB thermal pulse (TP) stars for a range of metallicities (Z) and explore the effect of convective boundary mixing (CBM). With decreasing metallicity and evolution along the TP phase, the He-shell flash and the third dredge-up (TDU) occur closer together in time. After some time (depending upon the CBM parameterisation), efficient TDU begins while the pulse-driven convection zone (PDCZ) is still present, causing a convective exchange of material between the PDCZ and the convective envelope. This results in the ingestion of protons into the convective He-burning pulse. Even small amounts of CBM encourage the interaction of the convection zones leading to transport of protons from the convective envelope into the He layer. H-burning luminosities exceed $10^9$ (in some cases $10^{10}$) $\mathrm{L}_\odot$. We also calculate models of dredge-out in the most massive super-AGB stars and show that the dredge-out phenomenon is another likely site of convective-reactive H-$^{12}$C combustion. We discuss the substantial uncertainties of stellar evolution models under these conditions. Nevertheless, the simulations suggest that in the convective-reactive H-combustion regime of H ingestion the star may encounter conditions for the intermediate neutron capture process (i process). We speculate that some CEMP-s/r stars could originate in i-process conditions in the H-ingestion phases of low-Z SAGB stars. This scenario would however suggest a very low electron-capture supernova rate from super-AGB stars. We also simulate potential outbursts triggered by such H-ingestion events, present their light curves and briefly discuss their transient properties.

We investigate whether the anomalous elemental abundance patterns in some of the C-enhanced metal-poor-s+r (CEMP-r/s) stars are consistent with predictions of nucleosynthesis yields from the i-process, a neutron-capture regime at neutron densities intermediate between those typical for the slow (s) and rapid (r) processes. Conditions necessary for the i-process are expected to be met at multiple stellar sites, such as the He-core and He-shell flashes in low-metallicity low-mass stars, super-AGB and post-AGB stars, as well as low-metallicity massive stars. We have found that single-exposure one-zone simulations of the i-process reproduce the abundance patterns in some of the CEMP-r/s stars much better than the model that assumes a superposition of yields from s- and r-process sources. Our previous study of nuclear data uncertainties relevant to the i-process revealed that they could have a significant impact on the i-process yields obtained in our idealized one-zone calculations, leading, for example, to ~0.7dex uncertainty in our predicted [Ba/La] ratio. Recent 3D hydrodynamic simulations of convection driven by a He-shell flash in post-AGB Sakurai's object have discovered a new mode of non-radial instabilities: the Global Oscillation of Shell H-ingestion. This has demonstrated that spherically symmetric stellar evolution simulations cannot be used to accurately model physical conditions for the i-process.

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

We demonstrate that among the potential sources of the primordial abundance variations of the proton-capture elements in globular-cluster stars proposed so far, such as the hot-bottom burning in massive AGB stars and H burning in the convective cores of supermassive and fast-rotating massive MS stars, only the supermassive MS stars with M > 10,000 Msun can explain all the observed abundance correlations without any fine-tuning of model parameters. We use our assumed chemical composition for the pristine gas in M13 (NGC6205) and its mixtures with 50% and 90% of the material partially processed in H burning in the 60,000 Msun MS model star as the initial compositions for the normal, intermediate and extreme populations of low-mass stars in this globular cluster, as suggested by its O-Na anti-correlation. We evolve these stars from the zero-age MS to the RGB tip with the thermohaline and parametric prescriptions for the RGB extra mixing. We find that the 3He-driven thermohaline convection cannot explain the evolutionary decline of [C/Fe] in M13 RGB stars, which, on the other hand, is well reproduced with the universal values for the mixing depth and rate calibrated using the observed decrease of [C/Fe] with MV in the globular cluster NGC5466 that does not have the primordial abundance variations.

Massive stars are key sources of radiative, kinetic, and chemical feedback in the universe. Grids of massive star models computed by different groups each using their own codes, input physics choices and numerical approximations, however, lead to inconsistent results for the same stars. We use three of these 1D codes---GENEC, KEPLER and MESA---to compute non-rotating stellar models of $15~\mathrm{M}_\odot$, $20~\mathrm{M}_\odot$, and $25~\mathrm{M}_\odot$ and compare their nucleosynthesis. We follow the evolution from the main sequence until the end of core helium burning. The GENEC and KEPLER models hold physics assumptions used in large grids of published models. The MESA code was set up to use convective core overshooting such that the CO core masses are consistent with those obtained by GENEC. For all models, full nucleosynthesis is computed using the NuGrid post-processing tool MPPNP. We find that the surface abundances predicted by the models are in reasonable agreement. In the helium core, the standard deviation of the elemental overproduction factors for Fe to Mo is less than $30\,\%$---smaller than the impact of the present nuclear physics uncertainties. For our three initial masses, the three stellar evolution codes yield consistent results. Differences in key properties of the models, e.g., helium and CO core masses and the time spent as a red supergiant, are traced back to the treatment of convection and, to a lesser extent, mass loss. The mixing processes in stars remain the key uncertainty in stellar modelling. Better constrained prescriptions are thus necessary to improve the predictive power of stellar evolution models.

When carbon is ignited off-centre in a CO core of a super-AGB star, its burning in a convective shell tends to propagate to the centre. Whether the C flame will actually be able to reach the centre depends on the efficiency of extra mixing beneath the C convective shell. Whereas thermohaline mixing is too inefficient to interfere with the C-flame propagation, convective boundary mixing can prevent the C burning from reaching the centre. As a result, a C-O-Ne white dwarf (WD) is formed, after the star has lost its envelope. Such a "hybrid" WD has a small CO core surrounded by a thick ONe zone. In our 1D stellar evolution computations, the hybrid WD is allowed to accrete C-rich material, as if it were in a close binary system and accreted H-rich material from its companion with a sufficiently high rate at which the accreted H would be processed into He under stationary conditions, assuming that He could then be transformed into C. When the mass of the accreting WD approaches the Chandrasekhar limit, we find a series of convective Urca shell flashes associated with high abundances of 23Na and 25Mg. They are followed by off-centre C ignition leading to convection that occupies almost the entire star. To model the Urca processes, we use the most recent well-resolved data for their reaction and neutrino-energy loss rates. Because of the emphasized uncertainty of the convective Urca process in our hybrid WD models of SN Ia progenitors, we consider a number of their potentially possible alternative instances for different mixing assumptions, all of which reach a phase of explosive C ignition, either off or in the centre. Our hybrid SN Ia progenitor models have much lower C to O abundance ratios at the moment of the explosive C ignition than their pure CO counterparts, which may explain the observed diversity of the SNe Ia.

When carbon is ignited off-center in a CO core of a super-AGB star, its burning in a convective shell tends to propagate to the center. Whether the C flame will actually be able to reach the center depends on the efficiency of extra mixing beneath the C convective shell. Whereas thermohaline mixing is too inefficient to interfere with the C-flame propagation, convective boundary mixing can prevent the C burning from reaching the center. As a result, a C-O-Ne white dwarf (WD) is formed, after the star has lost its envelope. Such a "hybrid" WD has a small CO core surrounded by a thick ONe zone. In our 1D stellar evolution computations the hybrid WD is allowed to accrete C-rich material, as if it were in a close binary system and accreted H-rich material from its companion with a sufficiently high rate at which the accreted H would be processed into He under stationary conditions, assuming that He could then be transformed into C. When the mass of the accreting WD approaches the Chandrasekhar limit, we find a series of convective Urca shell flashes associated with high abundances of 23Na and 25Mg. They are followed by off-center C ignition leading to convection that occupies almost the entire star. To model the Urca processes, we use the most recent well-resolved data for their reaction and neutrino-energy loss rates. Because of the emphasized uncertainty of the convective Urca process in our hybrid WD models of SN Ia progenitors, we consider a number of their potentially possible alternative instances for different mixing assumptions, all of which reach a phase of explosive C ignition, either off or in the center. Our hybrid SN Ia progenitor models have much lower C to O abundance ratios at the moment of the explosive C ignition than their pure CO counterparts, which may explain the observed diversity of the SNe Ia.

Recently an increasing number of studies were devoted to measure the abundances of neutron-capture elements heavier than iron in stars belonging to Galactic Open Clusters (OCs). OCs span a sizeable range in metallicity -0.6<[Fe/H]<+0.4), and they show abundances of light elements similar to disk stars of the same age. A different pattern is observed for heavy elements. A large scatter is observed for Ba, with most OCs showing [Ba/Fe] and [Ba/La] overabundant with respect to the Sun. The origin of this overabundance is not clearly understood. With the goal of providing new observational insights we determined radial velocities, atmospheric parameters and chemical composition of 27 giant stars members of five OCs: Cr 110, Cr 261, NGC 2477, NGC 2506 and NGC 5822. We used high-resolution spectra obtained with the UVES spectrograph at ESO Paranal. We perform a detailed spectroscopic analysis of these stars to measure the abundance of up to 22 elements per star. We study the dependence of element abundance on metallicity and age with unprecedented detail, complementing our analysis with data culled from the literature. We confirm the trend of Ba overabundance in OCs, and show its large dispersion for clusters younger than ~4 Gyr. Finally, the implications of our results for stellar nucleosynthesis are discussed. We show in this work that the Ba enrichment compared to other neutron-capture elements in OCs cannot be explained by the contributions from the slow neutron-capture process and the rapid neutron-capture process. Instead, we argue that this anomalous signature can be explained by assuming an additional contribution by the intermediate neutron-capture process.

V4334 Sgr (Sakurai's object) is an enigmatic evolved star that underwent a very late thermal pulse a few years before its discovery in 1996. It ejected a new, hydrogen-deficient nebula in the process. Emission lines from the newly ejected gas were first discovered in 1998 (He I 1083 nm) and 2001 (optical). We have monitored the optical emission spectrum since. From 2001 through 2007 the optical spectrum showed an exponential decline in flux, consistent with a shock that occurred around 1998 and started cooling soon after that. In this paper we show that since 2008 the line fluxes have been continuously rising again. Our preliminary interpretation is that this emission comes from a region close to the central star, and is excited by a second shock. This shock may have been induced by an increase in the stellar mass loss and wind velocity associated with a rise in the stellar temperature.

We performed 3-D simulations of proton-rich material entrainment into \czw-rich He-shell flash convection and the subsequent H-ingestion flash that took place in the post-AGB star Sakurai's object. Observations of the transient nature and anomalous abundance features are available to validate our method and assumptions, with the aim to apply them to very low metallicity stars in the future. We include nuclear energy feedback from H burning and cover the full $4\pi$ geometry of the shell. Runs on $768^3$ and $1536^3$ grids agree well with each other and have been followed for $1500\mathrm{min}$ and $1200\mathrm{min}$. After a $850\mathrm{min}$ long quiescent entrainment phase the simulations enter into a global non-spherical oscillation that is launched and sustained by individual ignition events of H-rich fluid pockets. Fast circumferential flows collide at the antipode and cause the formation and localized ignition of the next H-overabundant pocket. The cycle repeats for more than a dozen times while its amplitude decreases. During the global oscillation the entrainment rate increases temporarily by a factor $\approx 100$. Entrained entropy quenches convective motions in the upper layer until the burning of entrained H establishes a separate convection zone. The lower-resolution run hints at the possibility that another global oscillation, perhaps even more violent will follow. The location of the H-burning convection zone agrees with a 1-D model in which the mixing efficiency is calibrated to reproduce the light curve. The simulations have been performed at the NSF Blue Waters supercomputer at NCSA.

We present the first 3-dimensional, fully compressible gas-dynamics simulations in $4\pi$ geometry of He-shell flash convection with proton-rich fuel entrainment at the upper boundary. This work is motivated by the insufficiently understood observed consequences of the H-ingestion flash in post-AGB stars (Sakurai's object) and metal-poor AGB stars. Our investigation is focused on the entrainment process at the top convection boundary and on the subsequent advection of H-rich material into deeper layers, and we therefore ignore the burning of the proton-rich fuel in this study. We find that, for our deep convection zone, coherent convective motions of near global scale appear to dominate the flow. At the top boundary convective shear flows are stable against Kelvin-Helmholtz instabilities. However, such shear instabilities are induced by the boundary-layer separation in large-scale, opposing flows. This links the global nature of thick shell convection with the entrainment process. We establish the quantitative dependence of the entrainment rate on grid resolution. With our numerical technique simulations with $1024^3$ cells or more are required to reach a numerical fidelity appropriate for this problem. However, only the result from the $1536^3$ simulation provides a clear indication that we approach convergence with regard to the entrainment rate. Our results demonstrate that our method, which is described in detail, can provide quantitative results related to entrainment and convective boundary mixing in deep stellar interior environments with veryvstiff convective boundaries. For the representative case we study in detail, we find an entrainment rate of $4.38 \pm 1.48 \times 10^{-13}M_\odot \mathrm{/s}$.

Classical novae are the result of thermonuclear flashes of hydrogen accreted by CO or ONe white dwarfs, leading eventually to the dynamic ejection of the surface layers. These are observationally known to be enriched in heavy elements, such as C, O and Ne that must originate in layers below the H-flash convection zone. Building on our previous work, we now present stellar evolution simulations of ONe novae and provide a comprehensive comparison of our models with published ones. Some of our models include exponential convective boundary mixing to account for the observed enrichment of the nova ejecta even when accreted material has a solar abundance distribution. Our models produce maximum temperature evolution profiles and nucleosynthesis yields in good agreement with models that generate enriched ejecta by assuming that the accreted material was pre-mixed. We confirm for ONe novae the result we reported previously, i.e.\ we found that $^3$He could be produced {\it in situ} in solar-composition envelopes accreted with slow rates ($\dot{M} < 10^{-10}\,M_\odot/\mbox{yr}$) by cold ($T_{\rm WD} < 10^7$ K) CO WDs, and that convection was triggered by $^3$He burning before the nova outburst in that case. In addition, we now find that the interplay between the $^3$He production and destruction in the solar-composition envelope accreted with an intermediate rate, e.g.\ $\dot{M} = 10^{-10}\,M_\odot/\mbox{yr}$, by the $1.15\,M_\odot$ ONe WD with a relatively high initial central temperature, e.g.\ $T_{\rm WD} = 15\times 10^6$ K, leads to the formation of a thick radiative buffer zone that separates the bottom of the convective envelope from the WD surface. (Abridged)

We present barium, carbon, and silicon isotopic compositions of 38 acid-cleaned presolar SiC grains from Murchison. Comparison with previous data shows that acid washing is highly effective in removing barium contamination. Strong depletions in $\delta$($^{138}$Ba/$^{136}$Ba) values are found, down to $-$400 permil, which can only be modeled with a flatter $^{13}$C profile within the $^{13}$C pocket than is normally used. The dependence of $\delta$($^{138}$Ba/$^{136}$Ba) predictions on the distribution of $^{13}$C within the pocket in AGB models allows us to probe the $^{13}$C profile within the $^{13}$C pocket and the pocket mass in asymptotic giant branch (AGB) stars. In addition, we provide constraints on the $^{22}$Ne$(\alpha,n)^{25}$Mg rate in the stellar temperature regime relevant to AGB stars, based on $\delta$($^{134}$Ba/$^{136}$Ba) values of mainstream grains. We found two nominally mainstream grains with strongly negative $\delta$($^{134}$Ba/$^{136}$Ba) values that cannot be explained by any of the current AGB model calculations. Instead, such negative values are consistent with the intermediate neutron capture process ($i$-process), which is activated by the Very Late Thermal Pulse (VLTP) during the post-AGB phase and characterized by a neutron density much higher than the $s$-process. These two grains may have condensed around post-AGB stars. Finally, we report abundances of two $p$-process isotopes, $^{130}$Ba and $^{132}$Ba, in single SiC grains. These isotopes are destroyed in the $s$-process in AGB stars. By comparing their abundances with respect to that of $^{135}$Ba, we conclude that there is no measurable decay of $^{135}$Cs ($t_{1/2}$= 2.3 Ma) to $^{135}$Ba in individual SiC grains, indicating condensation of barium, but not cesium into SiC grains before $^{135}$Cs decayed.

Evolution of a supernova type Ia progenitor requires formation of a CO white dwarf, which implies a dependence on the C burning rate (CBR). It can also be affected by the recently identified possibility of C flame quenching by convective boundary mixing. We present first results of our study of the combined effect of these two potential sources of uncertainty on the SN Ia progenitor evolution. We consider the possibility that the CBR is higher than its currently recommended value by as much as a factor of 1000 if unidentified resonances are important, or that it is significantly lower because of the hindrance effect. For stellar models that assume the Schwarzschild boundary for convection, the maximum initial mass for the formation of CO WDs increases from M_i ~ 5.5 Msun for the CBR factor of 1000 to M_i > 7.0 Msun for the CBR factor of 0.01. For C-flame quenching models, hybrid C-O-Ne WDs form for a range of initial mass of Delta M_i ~ 1 Msun, which increases a fraction of stars that form WDs capable of igniting C in a thermonuclear runaway. The most extreme case is found for the CBR factor of 0.1 that is supported by the hindrance model. This nuclear physics assumption, combined with C flame quenching, leads to the formation of a hybrid C-O-Ne WD with a mass of 1.3 Msun. Such WDs do not need to accrete much mass to reach the Chandrasekhar limit.

Nuclear astrophysics simulations aiming to study the origin of the elements in stars require a multitude of nuclear physics input. Both systematic model dependent and statistically correlated uncertainties need to be considered. An application where realistic uncertainty assessments are especially important is the intermediate neutron capture process or i process: a neutron capture regime with neutron densities intermediate between the slow and rapid processes. Accordingly, the main network flux proceeds on the neutron-rich unstable isotopes up to 4-5 species off the valley of stability. The i process has been clearly identified to be active in post-AGB stars during the Very Late Thermal Pulse H-ingestion event, and a recent work infers about its important role in early generations of stars. Here we demonstrate the effect of propagating systematic nuclear uncertainties from different theoretical models to final abundances for a region around the 2$^\mathrm{nd}$ peak at $A-Z=80$ for elemental ratio predictions involving Ba, La and Eu in i-process conditions. These elements are used to distinguish different n-capture contributions observed in low-metallicity stars. For the simple 1-zone model adopted here, predictions vary as much as a factor of 22 in on possible observational plane ([La/Eu] vs. [Ba/La]). To consider statistically correlated uncertainties, we similarly perform a full nuclear physics uncertainty study within a given Hauser-Feshbach model and demonstrate the role of correlations on the final stellar abundance uncertainties. We show that in i-process conditions the main result of neglecting correlations is to underestimate the impact of nuclear uncertainties on the final nucleosynthesis yields by as much as two orders of magnitude. (continued in document)

Graphite is one of the many presolar circumstellar condensate species found in primitive meteorites. While the isotopic compositions of low-density graphite grains indicate an origin in core-collapse supernovae, some high-density grains have extreme isotopic anomalies in C, Ca and Ti, which cannot be explained by envelope predictions of asymptotic giant branch (AGB) stars or theoretical supernova models. The Ca and Ti isotopic anomalies, however, match the predictions of He-shell abundances in AGB stars. In this study, we show that the C, Ca, and Ti isotopic anomalies are consistent with nucleosynthesis predictions of the H-ingestion phase during a very late thermal pulse (VLTP) event in post-AGB stars. The low $^{12}$C/$^{13}$C isotopic ratios in these grains are a result of abundant $^{12}$C efficiently capturing the protons that are being ingested during the VLTP. Very high neutron densities of $\sim 10^{15}$ cm$^{-3}$, typical of the $i$-process, are achieved during this phase in post-AGB stars. The large $^{42,43,44}$Ca excesses in some graphite grains are indicative of neutron capture nucleosynthesis during VLTP. The comparison of VLTP nucleosynthesis calculations to the graphite data also indicate that apparent anomalies in the Ti isotopic ratios are due to large contributions from $^{46,48}$Ca, which cannot be resolved from the isobars $^{46,48}$Ti during the measurements. We conclude that presolar graphite grains with moderate to extreme Ca and Ti isotopic anomalies originate in post-AGB stars that suffer a very late thermal pulse.

We report C, Si, and S isotope measurements on 34 presolar silicon carbide grains of Type AB, characterized by 12C/13C < 10. Nitrogen, Mg-Al-, and Ca-Ti-isotopic compositions were measured on a subset of these grains. Three grains show large 32S excesses, a signature that has been previously observed for grains from supernovae (SNe). Enrichments in 32S may be due to contributions from the Si/S zone and the result of S molecule chemistry in still unmixed SN ejecta or due to incorporation of radioactive 32Si from C-rich explosive He shell ejecta. However, a SN origin remains unlikely for the three AB grains considered here, because of missing evidence for 44Ti, relatively low 26Al/27Al ratios (a few times 10-3), and radiogenic 32S along with low 12C/13C ratios. Instead, we show that born-again asymptotic giant branch (AGB) stars that have undergone a very-late thermal pulse (VLTP), known to have low 12C/13C ratios and enhanced abundances of the light s-process elements, can produce 32Si, which makes such stars attractive sources for AB grains with 32S excesses. This lends support to the proposal that at least some AB grains originate from born-again AGB stars, although uncertainties in the born-again AGB star models and possible variations of initial S-isotopic compositions in the parent stars of AB grains make it difficult to draw a definitive conclusion.

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

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

After off-center C ignition in the cores of super-AGB stars the C flame propagates all the way down to the center, trailing behind it the C-shell convective zone, and thus building a degenerate ONe core. This standard picture is obtained in stellar evolution simulations if the bottom C-shell convection boundary is assumed to be a discontinuity associated with a strict interpretation of the Schwarzschild condition for convective instability. However, this boundary is prone to additional mixing processes, such as thermohaline convection and convective boundary mixing. Using hydrodynamic simulations we show that, contrary to previous results, thermohaline mixing is too inefficient to interfere with the C-flame propagation. However, even a small amount of convective boundary mixing removes the physical conditions required for the C-flame propagation all the way to the center. This result holds even if we allow for some turbulent heat transport in the CBM region. As a result, super AGB stars build in their interiors hybrid C-O-Ne degenerate cores composed of a relatively large CO core (M_CO ~ 0.2 M_sun) surrounded by a thick ONe zone (M_ONe ~ 0.85 M_sun) with another thin CO layer above. If exposed by mass loss, these cores will become hybrid C-O-Ne white dwarfs. Otherwise, the ignition of C-rich material in the central core, surrounded by the thick ONe zone, may trigger a thermonuclear supernova explosion. The quenching of the C-flame may have implications for the ignition mechanism of SN Ia in the double-degenerate merger scenario.

The R Coronae Borealis (RCB) stars are hydrogen-deficient, variable stars that are most likely the result of He-CO WD mergers. They display extremely low oxygen isotopic ratios, 16O/18O ~ 1 - 10, 12C/13C>=100, and enhancements up to 2.6dex in F and in s-process elements from Zn to La, compared to solar. These abundances provide stringent constraints on the physical processes during and after the double-degenerate merger. As shown before O-isotopic ratios observed in RCB stars cannot result from the dynamic double-degenerate merger phase, and we investigate now the role of the long-term 1D spherical post-merger evolution and nucleosynthesis based on realistic hydrodynamic merger progenitor models. We adopt a model for extra envelope mixing to represent processes driven by rotation originating in the dynamical merger. Comprehensive nucleosynthesis post-processing simulations for these stellar evolution models reproduce, for the first time, the full range of the observed abundances for almost all the elements measured in RCB stars: 16O/18O ratios between 9 and 15, C-isotopic ratios above 100, and ~1.4 - 2.35dex F enhancements, along with enrichments in s-process elements. The nucleosynthesis processes in our models constrain the length and temperature in the dynamic merger shell-of-fire feature as well as the envelope mixing in the post-merger phase. s-process elements originate either in the shell-of-fire merger feature or during the post-merger evolution, but the contribution from the AGB progenitors is negligible. The post-merger envelope mixing must eventually cease ~ 10^6yr after the dynamic merger phase, before the star enters the RCB phase.

Classical novae are the results of surface thermonuclear explosions of hydrogen accreted by white dwarfs (WDs) from their low-mass main-sequence or red-giant binary companions. Chemical composition analysis of their ejecta shows that nova outbursts occur on both carbon-oxygen (CO) and more massive oxygen-neon (ONe) WDs, and that there is cross-boundary mixing between the accreted envelope and underlying WD. We demonstrate that the state-of-the-art stellar evolution code MESA and post-processing nucleosynthesis tools of NuGrid can successfully be used for modeling of CO and ONe nova outbursts and nucleosynthesis. The convective boundary mixing (CBM) in our 1D numerical simulations is implemented using a diffusion coefficient that is exponentially decreasing with a distance below the bottom of the convective envelope. We show that this prescription produces maximum temperature evolution profiles and nucleosynthesis yields in good agreement with those obtained using the commonly adopted 1D nova model in which the CBM is mimicked by assuming that the accreted envelope has been pre-mixed with WD's material. In a previous paper, we have found that 3He can be produced in situ in solar-composition envelopes accreted with slow rates (dM/dt < 1e-10 M_sun/yr) by cold (T_WD < 1d7 K) CO WDs, and that convection is triggered by 3He burning before the nova outburst in this case. Here, we confirm this result for ONe novae. Additionally, we find that the interplay between the 3He production and destruction in the solar-composition envelope accreted with an intermediate rate, e.g. dM/dt = 1e-10 M_sun/yr, by the 1.15 M_sun ONe WD with a relatively high initial central temperature, e.g. T_WD = 15e6 K, leads to the formation of a thick radiative buffer zone that separates the bottom of the convective envelope from the WD surface.

About a year after core collapse supernova, dust starts to condense in the ejecta. In meteorites, a fraction of C-rich presolar grains (e.g., silicon carbide (SiC) grains of Type-X and low density graphites) are identified as relics of these events, according to the anomalous isotopic abundances. Several features of these abundances remain unexplained and challenge the understanding of core-collapse supernovae explosions and nucleosynthesis. We show, for the first time, that most of the measured C-rich grain abundances can be accounted for in the C-rich material from explosive He burning in core-collapse supernovae with high shock velocities and consequent high temperatures. The inefficiency of the $^{12}$C($\alpha$,$\gamma$)$^{16}$O reaction relative to the rest of the $\alpha$-capture chain at $T > 3.5\times10^8 \mathrm{K}$ causes the deepest He-shell material to be carbon rich and silicon rich, and depleted in oxygen. The isotopic ratio predictions in part of this material, defined here as the C/Si zone, are in agreement with the grain data. The high-temperature explosive conditions that our models reach at the bottom of the He shell, can also be representative of the nucleosynthesis in hypernovae or in the high-temperature tail of a distribution of conditions in asymmetric supernovae. Finally, our predictions are consistent with the observation of large $^{44}$Ca/$^{40}$Ca observed in the grains. This is due to the production of $^{44}$Ti together with $^{40}$Ca in the C/Si zone, and/or to the strong depletion of $^{40}$Ca by neutron captures.

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

Classical novae are the results of surface thermonuclear explosions of H-rich material accreted by white dwarfs (WDs) from their low-mass main-sequence or red-giant binary companions. Chemical composition analysis of their ejecta shows that nova outbursts occur on both carbon-oxygen (CO) and oxygen-neon (ONe) WDs, and that there is cross-boundary mixing between the accreted envelope and underlying WD. We have combined the stellar evolution code MESA and post-processing nucleosynthesis tools of NuGrid into a framework that allows to produce up-to-date models of nova outbursts and compute detailed nucleosynthesis in novae occurring on CO and ONe WDs. The convective boundary mixing (CBM) in our 1D numerical simulations is implemented using a diffusion coefficient that is exponentially decreasing with a distance below the bottom of the convective envelope. This MESA CBM prescription is based on the findings in 3D hydrodynamic simulations that the velocity field, and along with it the mixing expressed in terms of a diffusion coefficient, decays exponentially in the stable layer adjacent to a convective boundary. The framework can also use the commonly adopted 1D nova model in which the CBM is mimicked by assuming that the accreted envelope has been pre-enriched with WD's material. In this preliminary report, we present the most interesting new results related to CO and ONe nova outbursts that have been obtained with the Nova Framework.

Novae are cataclysmic variables driven by accretion of H-rich material onto a white-dwarf (WD) star from its low-mass main-sequence binary companion. New time-domain observational capabilities, such as the Palomar Transient Factory and Pan-STARRS, have revealed a diversity of their behaviour that should be theoretically addressed. Nova outbursts depend sensitively on nuclear physics data, and more readily available nova simulations are needed in order to effectively prioritize experimental effort in nuclear astrophysics. In this paper we use the MESA stellar evolution code to construct multicycle nova evolution sequences with CO WD cores. We explore a range of WD masses and accretion rates as well as the effect of different cooling times before the onset of accretion. In addition, we study the dependence on the elemental abundance distribution of accreted material and convective boundary mixing at the core-envelope interface. Models with such convective boundary mixing display an enrichment of the accreted envelope with C and O from the underlying white dwarf that is commensurate with observations. We compare our results with the previous work and investigate a new scenario for novae with the 3He-triggered convection.

We study the response of giant stars to mass loss. One-dimensional simulations of red and asymptotic giant branch stars with mass loss rates from $10^{-3}$ up to a few \msun/yr show in no case any significant radius increase. The largest radius increase of 0.2% was found in the case with the lowest mass loss rate. For dynamical-timescale mass loss rates, that may be encountered during a common envelope phase, the evolution is not adiabatic. The superadiabatic outer layer of the giant's envelope has a local thermal timescale comparable to the dynamical timescale. Therefore, this layer has enough time to readjust thermally. Moreover, the giant star is driven out of hydrostatic equilibrium and evolves dynamically. In these cases no increase of the stellar radius with respect to its initial value is found. If the mass loss rate is high enough, the superadiabaticity of the outer layer is lost progressively and a radiative zone forms due to a combination of thermal and dynamical readjustment. Conditions for unstable mass transfer based on adiabatic mass loss models that predict a significant radius increase, may need to be re-evaluated.

[Abbreviated] We have investigated the color-magnitude diagram of Omega Centauri and find that the blue main sequence (bMS) can be reproduced only by models that have a of helium abundance in the range Y=0.35-$0.40. To explain the faint subgiant branch of the reddest stars ("MS-a/RG-a" sequence), isochrones for the observed metallicity ([Fe/H]\approx0.7) appear to require both a high age (~13Gyr) and enhanced CNO abundances ([CNO/Fe]\approx0.9$). Y~0.35 must also be assumed in order to counteract the effects of high CNO on turnoff colors, and thereby to obtain a good fit to the relatively blue turnoff of this stellar population. This suggest a short chemical evolution period of time (<1Gyr) for Omega Cen. Our intermediate-mass (super-)AGB models are able to reproduce the high helium abundances, along with [N/Fe]~2 and substantial O depletions if uncertainties in the treatment of convection are fully taken into account. These abundance features distinguish the bMS stars from the dominant [Fe/H] $\approx1.7$ population. The most massive super-AGB stellar models (M_zams>=6.8M_sun, M_He,core>=1.245M_sun) predict too large N-enhancements, which limits their role in contributing to the extreme populations. We show quantitatively that highly He- and N-enriched AGB ejecta have particularly efficient cooling properties. Based on these results and on the reconstruction of the orbit of Omega Cen with respect to the Milky Way we propose the galactic plane passage gas purging scenario for the chemical evolution of this cluster. Our model addresses the formation and properties of the bMS population (including their central location in the cluster). We follow our model descriptively through four passage events, which could explain not only some key properties of the bMS, but also of the MS-a/RGB-a and the s-enriched stars.

A leading formation scenario for R Coronae Borealis (RCB) stars invokes the merger of degenerate He and CO white dwarfs (WD) in a binary. The observed ratio of 16O/18O for RCB stars is in the range of 0.3-20 much smaller than the solar value of ~500. In this paper, we investigate whether such a low ratio can be obtained in simulations of the merger of a CO and a He white dwarf. We present the results of five 3-dimensional hydrodynamic simulations of the merger of a double white dwarf system where the total mass is 0.9 Mdot and the initial mass ratio (q) varies between 0.5 and 0.99. We identify in simulations with $q\lesssim0.7$ a feature around the merged stars where the temperatures and densities are suitable for forming 18O. However, more 16O is being dredged-up from the C- and O-rich accretor during the merger than the amount of 18O that is produced. Therefore, on a dynamical time scale over which our hydrodynamics simulation runs, a 16O/18O ratio of ~2000 in the "best" case is found. If the conditions found in the hydrodynamic simulations persist for 10^6 seconds the oxygen ratio drops to 16 in one case studied, while in a hundred years it drops to ~4 in another case studied, consistent with the observed values in RCB stars. Therefore, the merger of two white dwarfs remains a strong candidate for the formation of these enigmatic stars.

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

We report on oxygen abundances determined from medium-resolution near-IR spectroscopy for a sample of 57 carbon-enhanced metal-poor (CEMP) stars selected from the Hamburg/ESO survey. The majority of our program stars exhibit oxygen-to-iron ratios in the range +0.5 < [O/Fe]< +2.0. The [O/Fe] values for this sample are statistically compared to available high-resolution estimates for known CEMP stars, as well as to high-resolution estimates for a set of carbon-normal metal-poor stars. Carbon, nitrogen, and oxygen abundance patterns for a sub-sample of these stars are compared to yield predictions for very metal-poor asymptotic giant-branch abundances in the recent literature. We find that the majority of our sample exhibit patterns that are consistent with previously studied CEMP stars having s-process-element enhancements, and thus have very likely been polluted by carbon- and oxygen-enhanced material transferred from a metal-poor asymptotic giant-branch companion.

Depending on mass and metallicity as well as evolutionary phase, stars occasionally experience convective-reactive nucleosynthesis episodes. We specifically investigate the situation when nucleosynthetically unprocessed, H-rich material is convectively mixed with a He-burning zone, for example in convectively unstable shell on top of electron-degenerate cores in AGB stars, young white dwarfs or X-ray bursting neutron stars. Such episodes are frequently encountered in stellar evolution models of stars of extremely low or zero metal content [...] We focus on the convective-reactive episode in the very-late thermal pulse star Sakurai's object (V4334 Sagittarii). Asplund etal. (1999) determined the abundances of 28 elements, many of which are highly non-solar, ranging from H, He and Li all the way to Ba and La, plus the C isotopic ratio. Our simulations show that the mixing evolution according to standard, one-dimensional stellar evolution models implies neutron densities in the He that are too low to obtain a significant neutron capture nucleosynthesis on the heavy elements. We have carried out 3D hydrodynamic He-shell flash convection [...] we assume that the ingestion process of H into the He-shell convection zone leads only after some delay time to a sufficient entropy barrier that splits the convection zone [...] we obtain significantly higher neutron densities (~few 10^15 1/cm^3) and reproduce the key observed abundance trends found in Sakurai's object. These include an overproduction of Rb, Sr and Y by about 2 orders of magnitude higher than the overproduction of Ba and La. Such a peculiar nucleosynthesis signature is impossible to obtain with the mixing predictions in our one-dimensional stellar evolution models. [...] We determine how our results depend on uncertainties of nuclear reaction rates, for example for the C13(\alpha, n)O16 reaction.

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 the N, O, F and Na abundance and 12C/13C isotopic ratio measurements or upper limits for a sample of 10 C-rich, metal-poor giant stars, eight enhanced in s-process (CEMP-s) elements and two poor in n-capture elements (CEMP-no). The abundances are derived from IR, K-band, high-resolution CRIRES@VLT spectra obtained. The metallicity of our sample ranges from [Fe/H]=-3.4 to -1.3. F abundance could be measured only in two CEMP-s stars. With [F/Fe]=0.64, one is mildly F-overabundant, while the other is F-rich, at [F/Fe]=1.44. For the remaining eight objects, including both CEMP-no in our sample, only upper limits on F abundance could be placed. Our measurements and upper limits show that there is a spread in [F/C+N] ratio in CEMP-s stars as predicted by theory. Predictions from nucleosynthetic models for low-mass, low-metallicity Asymptotic Giant Branch stars, account for the derived F abundances, while the upper limits on F content derived for most of the stars are lower than the predicted values. The measured Na content is accounted for by AGB models in the 1.25 to 1.7 Msun range, confirming that the stars responsible for the peculiar abundance pattern observed in CEMP-s stars are low-mass, low-metallicity AGB stars, in agreement with the most accepted astrophysical scenario. We conclude that the mechanism of F production in current state-of-the-art low-metallicity low-mass AGB models needs further scrutiny and that F measurements in a larger number of metal-poor stars are needed to better constraint the models.

The {\alpha}-formalism is a common way to parametrize the common envelope interaction between a giant star and a more compact companion. The {\alpha} parameter describes the fraction of orbital energy released by the companion that is available to eject the giant star's envelope. By using new, detailed stellar evolutionary calculations we derive a user-friendly prescription for the {\lambda} parameter and an improved approximation for the envelope binding energy, thus revising the {\alpha} equation. We then determine {\alpha} both from simulations and observations in a self consistent manner. By using our own stellar structure models as well as population considerations to reconstruct the primary's parameters at the time of the common envelope interaction, we gain a deeper understanding of the uncertainties. We find that systems with very low values of q (the ratio of the companion's mass to the mass of the primary at the time of the common envelope interaction) have higher values of {\alpha}. A fit to the data suggests that lower mass companions are left at comparable or larger orbital separations to more massive companions. We conjecture that lower mass companions take longer than a stellar dynamical time to spiral in to the giant's core, and that this is key to allowing the giant to use its own thermal energy to help unbind its envelope. As a result, although systems with light companions might not have enough orbital energy to unbind the common envelope, they might stimulate a stellar reaction that results in the common envelope ejection.

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

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.

The He-shell flash convection in AGB stars is the site for the high-temperature component of the s-process in low- and intermediate mass giants, driven by the Ne22 neutron source. [...] The upper convection boundary plays a critical role during the H-ingestion episode that may lead to neutron-bursts in the most metal-poor AGB stars. We address these problems through global 3-dimensional hydrodynamic simulations including the entire spherical He-shell flash convection zone (as oposed to the 3D box-in-a-star simulations). An important aspect of our current effort is to establish the feasibility of our appoach. We explain why we favour the explicit treatment over the anelastic approximation for this problem. The simulations presented in this paper use a Cartesian grid of 512^3 cells and have been run on four 8-core workstations for four days to simulate ~5000s, which corresponds to almost ten convective turn-over times. The convection layer extends radially at the simulated point in the flash evolution over 7 H_p pressure scale-heights and exceeds the size of the underlying core. Convection is dominated by large convective cells that fill more than an entire octant. [...]

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.

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

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

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

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

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

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

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

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.

The goal of this paper is to analyze the impact of a primary neutron source on the s-process nucleosynthesis in massive stars at halo metallicity. Recent stellar models including rotation at very low metallicity predict a strong production of primary N14. Part of the nitrogen produced in the H-burning shell diffuses by rotational mixing into the He core where it is converted to Ne22 providing additional neutrons for the s process. We present nucleosynthesis calculations for a 25 Msun star at [Fe/H] = -3, -4, where in the convective core He-burning about 0.8 % in mass is made of primary Ne22. The usual weak s-process shape is changed by the additional neutron source with a peak between Sr and Ba, where the s-process yields increase by orders of magnitude with respect to the yields obtained without rotation. Iron seeds are fully consumed and the maximum production of Sr, Y and Zr is reached. On the other hand, the s-process efficiency beyond Sr and the ratio Sr/Ba are strongly affected by the amount of Ne22 and by nuclear uncertainties, first of all by the Ne22(alpha,n)Mg25 reaction. Finally, assuming that Ne22 is primary in the considered metallicity range, the s-process efficiency decreases with metallicity due to the effect of the major neutron poisons Mg25 and Ne22. This work represents a first step towards the study of primary neutron source effect in fast rotating massive stars, and its implications are discussed in the light of spectroscopic observations of heavy elements at halo metallicity.

While in the past spherodicity was assumed, and still is used in modeling of most nebulae, we know now that only a small number of planetary nebulae (PNe) are really spherical or at least nearly round. Round planetary nebulae are the minority of objects. In the case of those objects that underwent a very late helium flash (called VLTP objects or ``born-again'' PNe) it seems to be different. The first, hydrogen-rich PN, is more or less round. The ejecta from the VLTP event, in contrast, are extremely asymmetrical.

The prediction of the spins of the compact remnants is a fundamental goal of the theory of stellar evolution. Here, we confront the predictions for white dwarf spins from evolutionary models including rotation with observational constraints. We perform stellar evolution calculations for stars in the mass range 1... 3$\mso$, including the physics of rotation, from the zero age main sequence into the TP-AGB stage. We calculate two sets of model sequences, with and without inclusion of magnetic fields. From the final computed models of each sequence, we deduce the angular momenta and rotational velocities of the emerging white dwarfs. While models including magnetic torques predict white dwarf rotational velocities between 2 and 10 km s$^{-1}$, those from the non-magnetic sequences are found to be one to two orders of magnitude larger, well above empirical upper limits. We find the situation analogous to that in the neutron star progenitor mass range, and conclude that magnetic torques may be required in order to understand the slow rotation of compact stellar remnants in general.

In ESO period 65 (April-September 2000) the large programme 165.N-0276, led by Roger Cayrel, began making use of UVES at the Kueyen VLT telescope. Known within the Team and outside as ``First Stars'', it was aimed at obtaining high resolution, high signal-to-noise ratio spectra in the range 320 nm -- 1000 nm for a large sample of extremely metal-poor (EMP) stars identified from the HK objective prism survey.In this contribution we highlight the main results of the large programme.

The evolution of the first generations of stars at zero or extremly low metallicity, and especially some crucial properties like the primary N14 production, is charactarized by convective-reactive mixing events that are mostly absent from similar evolution phases at solar-like metallicity. These episodes occur when unprocessed H-rich material is mixed accross a convective boundary into C12 rich He-burning material, as for example in He-shell flashes of extremely-low metallicity AGB stars. In this paper we describe the astrophysical context of such convective-reactive events, including the difficulty of current one-dimensional stellar evolution models to correctly simulate these evolutionary phases. We then describe the requirements and current state of modeling convective-reactive processes in the first stars environment. We demonstrate some of the new concepts that we are applying to this problem, i.e. the highly accurate PPB advection scheme in the framework of PPM hydrodynamic simulations of mixing accross a very stiff convective boundary. We show initial results of such simulations that address the first non-reactive step of this problem, which is the entrainment of H at the top boundary of the He-shell flash convection zone.

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

The ejecta of the first probably very massive stars polluted the Big Bang primordial element mix with the first heavier elements. The resulting ultra metal-poor abundance distribution provided the initial conditions for the second stars of a wide range of initial masses reaching down to intermediate and low masses. The importance of these second stars for understanding the origin of the elements in the early universe are manifold. While the massive first stars have long vanished the second stars are still around and currently observed. They are the carriers of the information about the first stars, but they are also capable of nuclear production themselves. For example, in order to use ultra or extremely metal-poor stars as a probe for the r-process in the early universe a reliable model of the s-process in the second stars is needed. Eventually, the second stars may provide us with important clues on questions ranging from structure formation to how the stars actually make the elements, not only in the early but also in the present universe. In particular the C-rich extremely metal-poor stars, most of which show the s-process signature, are thought to be associated with chemical yields from the evolved giant phase of intermediate mass stars. Models of such AGB stars at extremely low metallicity now exist, and comparison with observation show important discrepancies, for example with regard to the synthesis of nitrogen. This may hint at burning and mixing aspects of extremely metal-poor evolved stars that are not yet included in the standard picture of evolution, as for example the hydrogen-ingestion flash. The second stars of intermediate mass may have also played an important role in the formation of heavy elements that form through slow neutron capture reaction chains (s-process). Comparison of models with observations reveal which aspects of the physics input and assumptions need to be improved. The s-process is a particularly useful diagnostic tool for probing the physical processes that are responsible for the creation of elements in stars, like for example rotation. As new observational techniques and strategies continue to penetrate the field, for example the multi-object spectroscopy, or the future spectroscopic surveys, the extremely metal-poor stars will play an increasingly important role to address some of the most fundamental and challenging, current questions of astronomy.

We investigate the reheating of the very late thermal pulse (VLTP) object V4334 Sgr (Sakurai's Object) using radio observations from the Very Large Array, and optical spectra obtained with the Very Large Telescope. We find a sudden rise of the radio flux at 5 and 8 GHz - from <= 90 micro-Jy and 80 +/- 30 micro-Jy in February 2005 to 320 micro-Jy and 280 micro-Jy in June 2006. Optical line emission is also evolving, but the emission lines are fading. The optical line emission and early radio flux are attributed to a fast shock (and not photoionization as was reported earlier) which occurred around 1998. The fading is due to post-shock cooling and recombination. The recent rapid increase in radio flux is evidence for the onset of photoionization of carbon starting around 2005. The current results indicate an increase in the stellar temperature to 12 kK in 2006. The mass ejected in the VLTP eruption is M_ej >= 1e-4 Msol, but could be as high as 1e-2 Msol, depending mainly on the distance and the clumping factor of the outflow. We derive a distance between 1.8 and 5 kpc. A high mass loss could expose the helium layer and yield abundances compatible with those of [WC] and PG1159 stars.

The mass loss properties of carbon AGB stars are not very well constrained at present. A variety of empirical or theoretical formulae with different parameterisations are available in the literature and the agreement between them is anything but good. These simple mass loss prescriptions are nontheless used in many models of stellar evolution without much consideration of their applicability in various cases. We present here an on-going project aiming at a better description of the mass loss, that could be used to improve stellar evolution models -- especially the evolution during the TP-AGB phase. As a first step, we have considered the mass loss evolution during a He-shell flash. Using stellar parameters sampled from a stellar evolutionary track, we have computed the time evolution of the atmospheric layers and wind acceleration region during a flash event with detailed frequency-dependent radiation-hydrodynamical models including dust formation. Based on these results, we have also simulated the subsequent long-term dynamical evolution of the circumstellar envelope (CSE), including the formation of a detached shell. We find that existing simple mass loss prescriptions imply different mass loss evolutions compared with our model. The second step of the project deals with the dependence of mass loss on the basic stellar parameters. At the moment we are computing a large grid of wind models for C-rich AGB stars. Preliminary results show that simple parameterisations are difficult to obtain in certain regions of the parameter space considered, due to strong non-linearities in the wind mechanism.

We study the late evolution of solar metallicity stars in the transition region between white dwarf formation and core collapse. This includes the super-asymptotic giant branch (super-AGB, SAGB) stars, which have massive enough cores to ignite carbon burning and form an oxygen-neon (ONe) core. The most massive SAGB stars have cores that may grow to the Chandrasekhar mass because of continued shell-burning. Their cores collapse, triggering a so called electron capture supernovae (ECSN). From stellar evolution models we find that the initial mass range for SAGB evolution is 7.5 ... 9.25\msun. We perform calculations with three different stellar evolution codes to investigate the sensitivity of this mass range to some of the uncertainties in current stellar models. The mass range significantly depends on the treatment of semiconvective mixing and convective overshooting. To consider the effect of a large number of thermal pulses, as expected in SAGB stars, we construct synthetic SAGB models that include a semi-analytical treatment of dredge-up, hot-bottom burning, and thermal pulse properties. This synthetic model enables us to compute the evolution of the main properties of SAGB stars from the onset of thermal pulses until the core reaches the Chandrasekhar mass or is uncovered by the stellar wind. Thereby, we determine the stellar initial mass ranges that produce ONe-white dwarfs and electron-capture supernovae. The latter is found to be 9.0 ... 9.25\msun for our fiducial model, implying that electron-capture supernovae would constitute about 4% of all supernovae in the local universe. Our synthetic approach allows us to explore the uncertainty of this number imposed by uncertainties in the third dredge-up efficiency and ABG mass loss rate. We find for ECSNe a upper limit of ~20% of all supernovae (abridged).

The origin of the so called 'detached shells' around AGB stars is not fully understood, but two common hypotheses state that these shells form either through the interaction of distinct wind phases or an eruptive mass loss associated with a He-shell flash. We present a model of the formation of detached shells around thermal pulse asymptotic giant branch (TP-AGB) stars, based on detailed modelling of mass loss and stellar evolution, leading to a combination of eruptive mass loss and wind interaction. The purpose of this paper is first of all to connect stellar evolution with wind and mass loss evolution and demonstrate its consistency with observations, but also to show how thin detached shells around TP-AGB stars can be formed. Previous attempts to link mass loss evolution with the formation of detached shells were based on approximate prescriptions for the mass loss and have not included detailed modelling of the wind formation as we do here. (abridged)

We have found that at least seven hydrogen-deficient carbon (HdC) and R Coronae Borealis (RCB) stars, have 16O/18O ratios close to and in some cases less than unity, values that are orders of magnitude lower than measured in other stars (the Solar value is 500). Greatly enhanced 18O is evident in every HdC and RCB we have measured that is cool enough to have detectable CO bands. The three HdC stars measured have 16O/18O < 1, lower values than any of the RCB stars. These discoveries are important clues in determining the evolutionary pathways of HdC and RCB stars, for which two models have been proposed: the double degenerate (white dwarf (WD) merger), and the final helium-shell flash (FF). No overproduction of 18O is expected in the FF scenario. We have quantitatively explored the idea that HdC and RCB stars originate in the mergers of CO- and He-WDs. The merger process is estimated to take only a few days, with accretion rates of 150 Msun/ yr producing temperatures at the base of the accreted envelope of 1.2 - 1.9 x 10^8 K. Analysis of a simplified one-zone calculation shows that nucleosynthesis in the dynamically accreting material may provide a suitable environment for a significant production of 18O, leading to very low values of 16O/18O, similar to those observed. We also find qualitative agreement with observed values of 12C/13C and with the CNO elemental ratios. H-admixture during the accretion process from the small H-rich C/O WD envelope may play an important role in producing the observed abundances. Overall our analysis shows that WD mergers may very well be the progenitors of O18-rich RCB and HdC stars, and that more detailed simulations and modeling are justified.

PG1159-035 is the prototype of the PG1159 spectral class which consists of extremely hot hydrogen-deficient (pre-) white dwarfs. It is also the prototype of the GW Vir variables, which are non-radial g-mode pulsators. The study of PG1159 stars reveals insight into stellar evolution and nucleosynthesis during AGB and post-AGB phases. We perform a quantitative spectral analysis of PG1159-035 focusing on the abundance determination of trace elements. We have taken high-resolution ultraviolet spectra of PG1159-035 with the Hubble Space Telescope and the Far Ultraviolet Spectroscopic Explorer. They are analysed with non-LTE line blanketed model atmospheres. We confirm the high effective temperature with high precision (Teff=140,000+/-5000 K) and the surface gravity of logg=7. For the first time we assess the abundances of silicon, phosphorus, sulfur, and iron. Silicon is about solar. For phosphorus we find an upper limit of solar abundance. A surprisingly strong depletion of sulfur (2% solar) is discovered. Iron is not detected, suggesting an upper limit of 30% solar. This coincides with the Fe deficiency found in other PG1159 stars. We redetermine the nitrogen abundance and find it to be lower by one dex compared to previous analyses. The sulfur depletion is in contradiction with current models of AGB star intershell nucleosynthesis. The iron deficiency confirms similar results for other PG1159 stars and is explained by the conversion of iron into heavier elements by n-capture in the s-processing environment of the precursor AGB star. However, the extent of the iron depletion is stronger than predicted by evolutionary models. The relatively low nitrogen abundance compared to other pulsating PG1159 stars weakens the role of nitrogen as a distinctive feature of pulsators and non-pulsators in the GW Vir instability strip.

Theoretical models of very metal-poor intermediate-mass Asymptotic Giant Branch (AGB) stars predict a large overabundance of primary nitrogen. The very metal-poor, carbon-enhanced, s-process-rich stars, which are thought to be the polluted companions of now-extinct AGB stars, provide direct tests of the predictions of these models. Recent studies of the carbon and nitrogen abundances in metal-poor stars have focused on the most carbon-rich stars, leading to a potential selection bias against stars that have been polluted by AGB stars that produced large amounts of nitrogen, and hence have small [C/N] ratios. We call these stars Nitrogen-Enhanced Metal-Poor (NEMP) stars, and define them as having [N/Fe] > +0.5 and [C/N] < -0.5. In this paper, we report on the [C/N] abundances of a sample of 21 carbon-enhanced stars, all but three of which have [C/Fe] < +2.0. If NEMP stars were made as easily as Carbon-Enhanced Metal-Poor (CEMP) stars, then we expected to find between two and seven NEMP stars. Instead, we found no NEMP stars in our sample. Therefore, this observational bias is not an important contributor to the apparent dearth of N-rich stars. Our [C/N] values are in the same range as values reported previously in the literature (-0.5 to +2.0), and all stars are in disagreement with the predicted [C/N] ratios for both low-mass and high-mass AGB stars. We suggest that the decrease in [C/N] from the low-mass AGB models is due to enhanced extra-mixing, while the lack of NEMP stars may be caused by unfavorable mass ratios in binaries or the difficulty of mass transfer in binary systems with large mass ratios.

From high-resolution VLT/UVES spectra (R~43,000), we determine abundances or upper limits for Li, C, N, O, and other important elements, as well as 12C/13C isotopic ratio for three Carbon enhanced metal poor stars. All three stars have -3.30 <= [Fe/H] <= -2.85 and moderate to high CNO abundances. CS 22958-042 is one of the most carbon-rich CEMP stars known ([C/Fe] = +3.2), while CS 29528-041 (one of the few N-enhanced metal-poor stars known) is one of the most nitrogen rich ([N/Fe] = +3.0). Oxygen is very high in CS 31080-095 ([O/Fe] = +2.35) and in CS 22958-042 ([O/Fe] = +1.35). All three stars exhibit [Sr/Fe] < 0; Ba is not detected in CS 22958-042 ([Ba/Fe] < -0.53), but it is moderately enhanced ([Ba/Fe] ~ 1) in the other two stars. CS 22958-042 displays one of the largest sodium overabundances yet found in CEMP stars ([Na/Fe] = +2.8). CS 22958-042 has 12C/13C = 9, similar to most other CEMP stars without enhanced neutron-capture elements, while 12C/13C <= 40 in CS 31080-095. CS 31080-095 and CS 29528-041 have A(Li) ~ 1.7, below the Spite Plateau, while Li is not detected in CS 22958-042. CS 22958-042 is a CEMP-no star, but the other two stars are in no known class of CEMP star and thus either constitute a new class or are a link between the CEMP-no and CEMP-s classes, adding complexity to the abundance patterns for CEMP stars. We interpret the abundance patterns in our stars to imply that current models for the presumed AGB binary progenitors lack an extra-mixing process, similar to those apparently operating in RGB stars.

New Spitzer imaging observations have revealed the structure around the Mira variable star R Hya to be a one-sided parabolic arc 100 arcsec to the West stretching from North to South. We successfully model R Hya and its surroundings in terms of an interaction of the stellar wind from an asymptotic giant branch (AGB) star with the interstellar medium (ISM) the star moves through. Our three-dimensional hydrodynamic simulation reproduces the structure as a bow shock into the oncoming ISM. We propose this as another explanation of detached shells around such stars which should be considered alongside current theories of internal origin. The simulation predicts the existence of a tail of ram-pressure-stripped AGB material stretching downstream. Indications for such a tail behind R Hya are seen in IRAS maps.

We present the first results of the MIRIAD (MIPS [Multiband Imaging Photometer for Spitzer] Infra-Red Imaging of AGB [asymptotic giant branch] Dustshells) project using the Spitzer Space Telescope. The primary aim of the project is to probe the material distribution in the extended circumstellar envelopes (CSE) of evolved stars and recover the fossil record of their mass loss history. Hence, we must map the whole of the CSEs plus the surrounding sky for background subtraction, while avoiding the central star that is brighter than the detector saturation limit. With our unique mapping strategy, we have achieved better than one MJy/sr sensitivity in three hours of integration and successfully detected a faint (< 5 MJy/sr), extended (~400 arcsec) far-infrared nebula around the AGB star R Hya. Based on the parabolic structure of the nebula, the direction of the space motion of the star with respect to the nebula shape, and the presence of extended H alpha emission co-spatial to the nebula, we suggest that the detected far-IR nebula is due to a bow shock at the interface of the interstellar medium and the AGB wind of this moving star. This is the first detection of the stellar-wind bow-shock interaction for an AGB star and exemplifies the potential of Spitzer as a tool to examine the detailed structure of extended far-IR nebulae around bright central sources. \

The evolution of central stars of planetary nebulae can proceed in several distinct ways, either leading to H-deficiency or to H-normal surface composition. Several new simulations of the evolution channels that lead to H-deficiency are now available, mainly the born-again scenarios that are triggered by a He-shell flash during the hot pre-white dwarf evolution phase. A realistic AGB progenitor evolution is important for correct HRD tracks, that allow mass determinations. New hydrodynamic simulations of He-shell flash convection including cases with H-ingestion are now performed, and allow a determination of the convective extra-mixing efficiency. This has direct consequences for the intershell abundance distribution of AGB stars that can be observed in the H-deficient CSPN.

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

We report on the progress of our on-going campaign to monitor the evolution of the VLTP objects V4334 Sgr and V605 Aql, as well as the suspected (V)LTP object CK Vul. V4334 Sgr does not show signs of increased ionization compared to our previous observations in 2004. We obtained the first radio detection of V605 Aql, indicating a strong increase in radio flux since 1987. We also present the first radio detection of CK Vul and discuss the expansion of the material ejected during the 1670 event.

The hydrogen-deficiency in extremely hot post-AGB stars of spectral class PG1159 is probably caused by a (very) late helium-shell flash or a AGB final thermal pulse that consumes the hydrogen envelope, exposing the usually-hidden intershell region. Thus, the photospheric element abundances of these stars allow to draw conclusions about details of nuclear burning and mixing processes in the precursor AGB stars. We compare predicted element abundances to those determined by quantitative spectral analyses performed with advanced non-LTE model atmospheres. A good qualitative and quantitative agreement is found for many species (He, C, N, O, Ne, F, Si) but discrepancies for others (P, S, Fe) point at shortcomings in stellar evolution models for AGB stars.

GW Vir variables are the pulsating members in the spectroscopic class of PG 1159 stars. In order to understand the characteristic differences between pulsating and non-pulsating PG 1159 stars, we analyse FUSE spectra of eleven objects, of which six are pulsating, by means of state-of-the-art NLTE model atmospheres. The numerous metal lines in the FUV spectra of these stars allow a precise determination of the photospheric parameters. We present here preliminary results of our analysis.

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

We review the observed properties of extremely hot hydrogen-deficient post-AGB stars of spectral type [WC] and PG1159. Their H-deficiency is probably caused by a (very) late helium-shell flash or a AGB final thermal pulse, laying bare interior stellar regions which are usually kept hidden below the hydrogen envelope. Thus, the photospheric element abundances of these stars allow to draw conclusions about details of nuclear burning and mixing processes in the precursor AGB stars. We summarize the state-of-the-art of stellar evolution models which simulate AGB evolution and the occurrence of a late He-shell flash. We compare predicted element abundances to those determined by quantitative spectral analyses performed with advanced non-LTE model atmospheres. A good qualitative and quantitative agreement is found. Future work can contribute to an even more complete picture of the nuclear processes in AGB stars.

Calculations that demonstrate the influence of three key nuclear reaction rates on the evolution of Asymptotic Giant Branch stars have been carried out. We study the case of a star with an initial mass of 2Msun and a metallicity of Z=0.01, somewhat less than the solar metallicity. The dredge-up of nuclear processed material from the interior of the star, and the yield predictions for carbon, are sensitive to the rate of the N14(p,gamma)O15 and triple-alpha reactions. These reactions dominate the H- and He-burning shells of stars in this late evolutionary phase. Published uncertainty estimates for each of these two rates propagated through stellar evolution calculations cause uncertainties in carbon enrichment and yield predictions of about a factor of two. The other important He-burning reaction C12(alpha,gamma)O16, although associated with the largest uncertainty in our study, does not have a significant influence on the abundance evolution compared to other modelling uncertainties. This finding remains valid when the entire evolution from the main-sequence to the tip of the AGB is considered. We discuss the experimental sources of the rate uncertainties addressed here, and give some outlook for future work.

We report the discovery of a uniquely large excess of 18O in the hydrogen-deficient carbon (HdC) star, HD 137613, based on a spectrum of the first overtone bands of CO at 2.3-2.4 micron in which three strong absorption bands of 12C18O are clearly present. Bands of 12C16O also are present but no bands of 13C16O or 12C17O are seen. We estimate an isotopic ratio 16O/18O ~1. The Solar value of this ratio is ~500. Neither He-core burning nor He-shell flash burning can produce the isotopic ratios of oxygen and carbon observed in HD 137613. However, a remarkable similarity exists between the observed abundances and those found in the outer layers of the broad He-shell of early-AGB stars, soon after the end of He-core burning. It is not known how the outer envelope down to the He-shell could be lost but some mechanism of enhanced mass loss must be involved. HD 137613 may be a post-early-AGB star with the outer layers of the former He-burning shell as its photosphere. The unusual elemental abundances of the HdC stars resemble those of the R Coronae Borealis (RCB) stars, but HdC stars do not produce clouds of dust that produce declines in brightness. None of the other RCB or HdC stars observed shows significant 18O.

95% of all stars end their lives as white dwarfs. About 20% of the hot post-AGB stars are hydrogen deficient. Most of these are the result of a late helium-shell flash, but the evolutionary status of a fraction of about 10-20% of the hottest hydrogen-deficient stars, namely four O(He) stars, is as yet unexplained. They could be the long-searched hot successors of RCrB stars, which have not been identified up to now. If this turns out to be true, then a third post-AGB evolutionary sequence is revealed, which is probably the result of a double degenerate merging process. More generally, understanding details of merging double degenerate stars is of interest in the context of SN Ia events and hence cosmology.

We have studied how the third dredge-up and the carbon star formation in low-mass Asymptotic Giant Branch stars depends on certain key nuclear reaction rates. We find from a set of complete stellar evolution calculations of a 2Msun model with Z=0.01 including mass loss, that varying either the N14(p,g)O15 or the 3-alpha reaction rate within their uncertainties as given in the NACRE compilation results in dredge-up and yields that differ by a factor of 2. Model tracks with a higher rate for the 3-alpha rate and a lower rate for the N14(p,g)O15 reaction both show more efficient third dredge-up. New experimental results for the N14(p,g)O15 reaction rates are surveyed, yielding a rate which is about 40% lower than the tabulated NACRE rate, and smaller than NACRE's lower limit. We discuss the possible implications of the revised nuclear reaction stellar evolution calculations that aim to reproduce the observed carbon star formation at low mass, which requires efficient third dredge-up.

Intermediate mass stellar evolution tracks from the main sequence to the tip of the AGB for five initial masses (2 to 6Msun) and metallicity Z=0.0001 have been computed. The detailed 1D structure and evolution models include exponential overshooting, mass loss and a detailed nucleosynthesis network with updated nuclear reaction rates. The network includes a two-particle heavy neutron sink for approximating neutron density in the He-shell flash. It is shown how the neutron-capture nucleosynthesis is important in models of very low metallicity for the formation of light neutron-heavy species, like sodium or the heavy neon and magnesium isotopes. The models have high resolution, as required for modeling the third dredge-up. All sequences have been followed from the pre-main sequence to the end of the AGB when all envelope mass is lost. Detailed structural and chemical model properties as well as yields are presented. This set of stellar models is based on standard assumptions and updated input physics. It can be confronted with observations of extremely-metal poor stars and may be used to assess the role of AGB stars in the origin of abundance anomalies of some Globular Cluster members of correspondingly low metallicity.

(abbreviated) The evolution of intermediate mass stars at very low metallicity during their final thermal pulse asymptotic giant branch phase is studied in detail. As representative examples models with initial masses of 4Msun and 5Msun with a metallicity of Z=0.0001 ([Fe/H] ~ -2.3) are discussed. The 1D stellar structure and evolution model includes time- and depth dependent overshooting motivated by hydrodynamical simulations, as well as a full nuclear network and time-dependent mixing. Particular attention is given to high time and space resolution to avoid numerical artefacts related to third dredge-up and hot-bottom burning predictions. The model calculations predict very efficient third dredge-up which mixes the envelope with the entire intershell layer or a large fraction thereof, and in some cases penetrates into the C/O core below the He-shell. In all cases primary oxygen is mixed into the envelope. The models predict efficient envelope burning during the interpulse phase. Depending on the envelope burning temperature, oxygen is destroyed to varying degrees. The combined effect of dredge-up and envelope burning does not lead to any significant oxygen depletion in any of the cases considered in this study. The large dredge-up efficiency in our model is closely related to the particular properties of the H-shell during the dredge-up phase in low-metallicity very metal poor stars, which is followed here over many thermal pulses. During the dredge-up phase, the temperature just below the convective boundary is large enough for protons to burn vigorously when they are brought into the C-rich environment below the convection boundary by the time- and depth dependent overshooting. H-burning luminosities of 10^5 to ~2* 10^6L_sun are generated. [...]

While the basic properties of AGB stellar evolution are well established, comprehensive observational studies of late phases of intermediate mass stars continue to generate puzzles for current stellar models. Here, I review current techniques to model AGB stars, and I discuss important aspects of current research of AGB (and post-AGB) stellar evolution with a particular focus on how these interrelate.

(abridged) We model the nucleosynthesis during the thermal pulse phase of a rotating, solar metallicity AGB star of 3M_sun. Rotationally induced mixing during the thermal pulses produces a layer (~2E-5M_sun) on top of the CO-core where large amounts of protons and C12 co-exist. We follow the abundance evolution in this layer, in particular that of the neutron source C13 and of the neutron poison N14. In our AGB model mixing persists during the entire interpulse phase due to the steep angular velocity gradient at the core-envelope interface. We follow the neutron production during the interpulse phase, and find a resulting maximum neutron exposure of tau_max =0.04 mbarn^-1, which is too small to produce any significant s-process. In parametric models, we then investigate the combined effects of diffusive overshooting from the convective envelope and rotationally induced mixing. Models with overshoot and weaker interpulse mixing - as perhaps expected from more slowly rotating stars - yield larger neutron exposures. We conclude that the incorporation of rotationally induce mixing processes has important consequences for the production of heavy elements in AGB stars. Through a distribution of initial rotation rates it may lead to a natural spread in the neutron exposures obtained in AGB stars of a given mass - as appears to be required by observations. Our results suggest that both processes, diffusive overshoot and rotational mixing, may be required to obtain a consistent description of the s-process in AGB stars which fulfils all observational constraints. Finally, we find that mixing due to rotation within our current framework does increase the production of N15 in the partial mixing zone, however still falling short of what seems required by observations.

We present full stellar evolution and parametric models of the surface abundance evolution of O16, Ne22, Na23 and the magnesium isotopes in an extremely metal-poor intermediate mass star M_ZAMS=5M_sun, Z=0.0001. O16 and Ne22 are injected into the envelope by the third dredge-up following thermal pulses on the asymptotic giant branch. These species and the initially present Mg24 are depleted by hot bottom burning (HBB) during the interpulse phase. As a result, Na23, Mg25 and Mg26 are enhanced. If the HBB temperatures are sufficiently high for this process to deplete oxygen efficiently, Na23 is first produced and then depleted during the interpulse phase. Although the simultaneous depletion of O16 and enhancement of Na23 is possible, the required fine tuning of the dredge-up and HBB casts some doubt on the robustness of this process as the origin of the O-Na anti-correlation observed in globular cluster stars. However, a very robust prediction of our models are low Mg24/Mg25 and Mg24/Mg26 ratios whenever significant O16 depletion can be achieved. This seems to be in stark contrast with recent observations of the magnesium isotopic ratios in the globular cluster NGC6752.

We model the nucleosynthesis during a radiative interpulse phase of a rotating 3Msun Asymptotic Giant Branch (AGB) star. We find an enhanced production of the neutron source species C13 compared to non-rotating models due to shear mixing of protons and C12 at the core-envelope interface. We estimate that the resulting total production of heavy elements by slow neutron capture s-process is too low to account for most observations. This due to the fact that rotationally induced mixing during the interpulse phase causes a pollution of the C13 pocket layer with the neutron poison N14. As a result we find a maximum neutron exposure of tau_max=0.04mbarn^{-1} in the s-process layer of our solar metallicity model with rotation. This is about a factor of 5 to 10 less than required to reproduce the observed stellar s-process abundance patterns. We compare our results with models that include hydrodynamic overshooting mixing, and with simple parametric models including the combined effects of overshooting and mixing in the interpulse. Within the parametric model a range of mixing efficiencies during the interpulse phase correlates with a spread in the s-process efficiency. Such a spread is observed in AGB and post-AGB stars as well as in pre-solar SiC grains.

The oldest stars hide information about the chemical evolution of the early universe. In this contect new models of the AGB stellar evolution phase at very low and zero metallicity are presented. Due to the deficiency or absence of CNO catalytic material hydrogen burning operates at significantly higher temperatures than in stars of higher metallicity. As a result convective mixing plays a very important role since the nuclear burning shells are not well separated by an entropy barrier. These stars can host flash-like burning events induced by mixing of C12 and protons at the core-envelope interface. Three different mixing events with proton-capture nucleosynthesis can be encountered. One of them, the hot dredge-up, is reported here for the first time. All these specific modes of nucleosynthesis are relevant for the production of CNO material in the first intermediate mass stars.

[abridged] We study the s-process in AGB stars using three different stellar evolutionary models computed for a 3Msun and solar metallicity star. First we investigate the formation and the efficiency of the main neutron source. We parametrically vary the number of protons mixed from the envelope into the C12 rich core. For p/C12 > 0.3, mainly N14 is produced, which represent a major neutron poison. The amount of C12 in the He intershell and the maximum value of the time-integrated neutron flux are proportional. Then we generate detailed s-process calculations on the basis of stellar evolutionary models constructed with three different codes. One code considers convective hydrodynamic overshoot that depends on a free parameter f, and results in partial mixing beyond convective boundaries, the most efficient third dredge up and the formation of the C13 pocket. For the other two codes an identical C13 pocket is introduced in the post-processing nucleosynthesis calculations. The models generally reproduce the spectroscopically observed s-process enhancements. The results of the cases without overshoot are remarkably similar. The code including hydrodynamic overshoot produces a He intershell composition near to that observed in H-deficient central stars of planetary nebulae. As a result of this intershell dredge up the neutron fluxes have a higher efficiency, both during the interpulse periods and within thermal pulses. The s-element distribution is pushed toward the heavier s-process elements and large abundances of neutron-rich isotopes fed by branching points in the s-process path are produced. Several observational constraints are better matched by the models without overshoot. Our study need to be extended to different masses and metallicities and in the space of the free overshoot parameter f.

We describe a general yet simple method to analyse the propagation of nuclear reaction rate uncertainties in a stellar nucleosynthesis and mixing context. The method combines post-processing nucleosynthesis and mixing calculations with a Monte Carlo scheme. With this approach we reanalyze the dependence of theoretical oxygen isotopic ratio predictions in first dredge-up red giant branch stars in a systematic way. Such predictions are important to the interpretation of pre-solar Al_2 O_3 grains from meteorites. The reaction rates with uncertainties were taken from the NACRE compilation (Angulo etal., 1999). We include seven reaction rates in our systematic analysis of stellar models with initial masses from 1 to 3 M_sun. We find that the uncertainty of reaction rate for reaction O18(p,alpha)N15 typically causes an error in the theoretical O16/O18 ratio of about +20/-5 per cent. The error of the O16/O17 prediction is +-10 to 40 per cent depending on the stellar mass, and is persistently dominated by the comparatively small uncertainty of the O16(p,gamma)F17 reaction. With the new estimates on reaction rate uncertainties by the NACRE compilation, the p-capture reactions O17(p,alpha)N14 and O17(p,gamma)F18 have virtually no impact on theoretical predictions for stellar mass <= 1.5 M_sun. However, this is not the case for masses > 1.5 M_sun, where core mixing and subsequent envelope mixing interact. In these cases where core mixing complicates post-dredge-up surface abundances, uncertainty in other reactions have a secondary but noticeable effect on surface abundances.

We investigate the evolutionary status of four stars: V348 Sgr, DY Cen and MV Sgr in the Galaxy and HV 2671 in the LMC. These stars have in common random deep declines in visual brightness which are characteristic for R Coronae Borealis (RCB) stars. RCB stars are typically cool, hydrogen deficient supergiants. The four stars studied in this paper are hotter (T$_{\rm eff}$ = 15-20 kK) than the majority of RCB stars (T$_{\rm eff}$ = 5000-7000 K). Although these are commonly grouped together as the \emph{hot RCB stars} they do not necessarily share a common evolutionary history. We present new observational data and an extensive collection of archival and previously-published data which is reassessed to ensure internal consistency. We find temporal variations of various properties on different time scales which will eventually help us to uncover the evolutionary history of these objects. DY Cen and MV Sgr have typical RCB helium abundances which excludes any currently known post-AGB evolutionary models. Moreover, their carbon and nitrogen abundances present us with further problems for their interpretation. V348 Sgr and HV 2671 are in general agreement with a born-again post-AGB evolution and their abundances are similar to Wolf-Rayet central stars of PN. The three Galactic stars in the sample have circumstellar nebulae which produce forbidden line radiation (for HV 2671 we have no information). V348 Sgr and DY Cen have low density, low expansion velocity nebulae (resolved in the case of V348 Sgr), while MV Sgr has a higher density, higher expansion velocity nebula.

The carbon and nitrogen isotopic ratios of pre-solar SiC grains of type A+B suggest a proton-limited nucleosynthetic process as encountered, for instance, during the very late thermal pulse of post-AGB stars. We study the nuclear processes during this phase and find carbon and nitrogen isotopic ratios which can reproduce those of A+B grains. These results are still preliminary because they depend on uncertain factors such as the details of mixing during the post-AGB thermal pulse, the rates of some nuclear reactions, and the assumptions on mixing during the progenitor AGB phase.

H-deficient post-AGB objects, e.g. PG1159 type star K1-16 and born-again AGB star Sakurai's object, have been reported to be significantly iron-deficient. We find that the iron deficiencies expected due to neutron-capture nucleosynthesis during either the progenitor AGB evolution and/or the neutron burst that occurs as a result of the rapid burning of protons during a post-AGB He-flash are generally in line with observations.

We present a stellar evolution model sequence of the past and future evolution of the post-AGB born again star Sakurai's object (V4334 Sgr). In order to match the short evolutionary time scale we have to assume that the convective ingestion of hydrogen-rich envelope material into the He-flash convection zone proceeds slower than predicted by the mixing length theory. For the future we predict a swift second evolution through the central star region before a second born-again evolution occurs.

We discuss the occurrence of the s-process during the radiative interpulse phase of rotating AGB stars. Due to differential rotation, protons are mixed into C12-rich layers after thermal pulses, in the course of the so called third dredge up episode. We follow the time evolution of key isotope abundances in the relevant layers with a post-processing code which includes time dependant mixing and nucleosynthesis. In rotating AGB models, the mixing persists during the entire interpulse phase due to the steep gradient of angular velocity at the envelope-core interface. As the layers containing protons and C12, which are formed this way, become hotter, a C13-pocket is formed in a natural way. However, in this situation also N14 is formed and spread over the entire C13-pocket. We include the neutron consuming N14(n,p) reaction in our network and determine to what extent it reduces the production of trans-iron elements. We propose that rotation may be responsible for the spread of efficiencies of the C13 neutron source as required by observations.

Sakurai's object (V4334 Sgr) is a born again AGB star following a very late thermal pulse. So far no stellar evolution models have been able to explain the extremely fast evolution of this star, which has taken it from the pre-white dwarf stage to its current appearance as a giant within only a few years. A very high stellar mass can be ruled out as the cause of the fast evolution. Instead the evolution time scale is reproduced in stellar models by making the assumption that the efficiency for element mixing in the He-flash convection zone during the very late thermal pulse is smaller than predicted by the mixing-length theory. As a result the main energy generation from fast proton capture occurs closer to the surface and the expansion to the giant state is accelerated to a few years. Assuming a mass of V4334 Sgr of 0.604Msun -- which is consistent with a distance of 4kpc -- a reduction of the mixing length theory mixing efficiency by a factor of ~ 100 is required to match its evolutionary time scale. This value decreases if V4334 Sgr has a smaller mass and accordingly a smaller distance. However, the effect does not disappear for the smallest possible masses. These findings may present a semi-empirical constraint on the element mixing in convective zones of the stellar interior.

Sakurai's object is a born again AGB star of the very late thermal pulse flavor. In this contribution I will discuss new models of stellar evolution and nucleosynthesis models of this phase. Two most intriguing properties of Sakurai's objects have so far not been understood theoretically: the peculiar chemical appearance, in particular the high lithium abundance and the short time scale of only a few years on which the transition from the dwarf configuration into the born again giant appearance has occurred. A new nucleosynthesis mode of hot hydrogen-deficient He3 burning can explain the extraordinary lithium abundance. During the thermal pulse He3 is ingested from the envelope together with the protons into the hot He-flash convection zone. The first network calculations show that due to the large C12 abundance protons are rather captured by carbon than destroy newly formed Be7 and ultimately Li7. Moreover, the short evolution time scale has been reproduced by making the assumption that the convective efficiency for element mixing is smaller by two to three orders of magnitude than predicted by the mixing-length theory. As a result the main energy generation from fast convective proton capture will occur at a larger mass coordinate, closer to the surface and the expansion to the giant state is accelerated to a few years in agreement with Sakurai's behavior. This result may represent an independent empirical constraint on the poorly known efficiency of element mixing in convective zones of the stellar interior.

I discuss recent new models of post-Asymptotic Giant Branch stellar evolution. These models aim to clarify the evolutionary origin and status of a variety of hydrogen-deficient post-AGB stars such as central stars of planetary nebulae of Wolf-Rayet spectral type, PG1159 stars or Sakurai's object. Starting with AGB models with overshoot such stars can evolve through one of four distinct channels. Each of these channels has typical abundance patterns depending on the relative timing of the departure from the AGB and the occurrence of the last thermal pulse. I discuss the responsible mechanisms and observational counterparts.

A thermal pulse during the post-AGB phase of stellar evolution may lead to a unique mode of light element nucleosynthesis. The stage is set by the ingestion of the unprocessed envelope material into the hot He-flash convection zone below. If the temperature is sufficiently large and the C12 abundance high enough (e.g. T_8 > 0.8, X(C12) ~ 0.4 and X(H) ~ 1E-3) protons react faster with C12 and form C13 than destroying Be7 The latter forms by alpha-capture of He3 after an initial reduction of the He3 abundance to about 3E-5 X(He4) by the ppI reaction He3(He3,2p)He4 (for T_8 ~ 1). All He3 is burned within minutes to weeks depending on the temperature. Be7 is now present at about the previously mentioned level of He3. Its further fate is determined by the reactions Be7(e-,nu)Li7 and the alpha-capture reactions of Be7 and Li7. These captures lead to the production of Be11 which in turn is finally destroyed by Be11(alpha,n)N14$. The details of this mechanism of light element production in real stars is expected to be fairly dependent on the description of mixing.

Within the framework of the current models for the slow neutron capture (s) process in Asymptotic Giant Branch (AGB) stars we explore the uncertainties introduced by the assumptions made on stellar modeling. On the basis of star models constructed with three different evolutionary codes we generate detailed neutron-capture nucleosynthesis post-processing models. The main difference among the codes is that one of them includes an overshooting mechanism. As a result, the neutron fluxes are stronger both during the interpulse periods and within thermal pulses. The efficiency of the C13 source is also studied and we find that a linear relationship exists between the initial C12 in the intershell and the maximum number of neutrons produced.

(abridged) Observations clearly show that low-mass AGB stars can provide a nucleosynthesis site of the s-process. Recent stellar evolution models indicate that radiative burning of C13 between thermal pulses in low-mass AGB stars may indeed provide the needed neutrons. Some mixing between the proton-rich envelope and the carbon-rich core may lead to the production of C13. However, the responsible physical mechanism is not yet unambiguously identified. We present stellar model calculations with overshoot and rotation. Overshoot, with a time-dependent and exponentially decaying efficiency, leads to a partial mixture of protons and C12 during the third dredge-up. According to the depth-dependent ratio of protons and C12, a small C13-pocket forms underneath a N14-rich layer. Overshoot does not allow for any mixing during the interpulse phase. Rotation introduces mixing driven by large angular velocity gradients which form at the envelope-core interface in AGB stars, in particular after a thermal pulse. This leads to partial mixing after a pulse, as in the case of overshoot. However, rotation continues to mix the region of the C13-pocket with a diffusion coefficient of log D ~ 2...3 [cm**2 / s]. This does not only spread the C13-pocket, but also the more massive N14-rich layer, and finally leads to mixture of both layers. By the time when the temperature there has risen to about 9*10E7K and neutron production sets in, the N14 abundance exceeds the C13 abundance by a factor of 5...10. We analyze the role of N14 as a neutron poison by considering the recycling of neutrons via N14(n,p)C14 and C12(p,gamma)N13(beta+)C13 qualitatively. We find that as long as X(N14) << X(C12), the s-process will still be possible to occur under radiative conditions.

The influence of extended convective mixing (overshoot) on asymptotic giant branch stellar evolution is investigated in detail. The extended mixing is treated time-dependently, and the efficiency declines exponentially with the geometric distance from the convective boundary. It has been considered at all convective boundaries, including the He-flash convection zone in the intershell region which forms during the thermal pulses. Both the structural and the chemical evolution are affected by the inclusion of overshoot. The main results include a very efficient third dredge-up which leads to the formation of carbon stars of low mass and luminosity. A C13 pocket which may serve as a neutron source for the s-process can form after the third dredge-up has reached into the C12 rich intershell. Overshoot applied to the pulse-driven convective zone during the He-flash leads to a deeper penetration of the bottom of this convective zone into the C/O core below the He-burning shell. This in turn causes He4 to be less abundant in the intershell while C12 and O16 are more abundant compared to calculations without overshoot. We show that overshoot at the He-flash convection zone as well as at the base of the envelope convection enhance the efficiency of the third dredge-up. Characteristic properties of the structural and chemical evolution of AGB stars are presented.

We calculated models of massive AGB stars with a self-consistent coupling of time-dependent mixing and nuclear burning for 30 isotopes and 74 reactions. Overshoot with an exponentially declining velocity field was considered and applied during all stages of evolution and in all convective regions. Very efficient 3rd dredge-up was found even overcompensating the growth of the hydrogen-exhausted core after a few thermal pulses. Hot bottom burning occurs for M>4Msol within the sequences with overshoot. Carbon star formation in these more massive AGB stars is delayed or even prevented by hot bottom burning despite the very efficient dredge-up. With the simultaneous treatment of mixing and burning the formation of Li-rich AGB stars due to the Cameron-Fowler mechanism was followed. For a 6Msol model the maximum Li abundance was found to be epsilon(Li7)=4.4.

Recent advances in constructing stellar evolution models of hydrogen-deficient post-asymptotic giant branch (AGB) stars are presented. Hydrogen-deficiency can originate from mixing and subsequent convective burning of protons in the deeper layers during a thermal pulse on the post-AGB (VLTP). Dredge-up alone may also be responsible for hydrogen-deficiency of post-AGB stars. Models of the last thermal pulse on the AGB with very small envelope masses have shown efficient third dredge-up. The hydrogen content of the envelope is diluted sufficiently to produce H-deficient post-AGB stars (AFTP). Moreover, dredge-up alone may also cause H-deficiency during the Born-again phase (LTP). During the second AGB phase a convective envelope develops. A previously unknown lithium enrichment at the surface of Born-again stellar models may be used to distinguish between objects with different post-AGB evolution. The observed abundance ratios of C, O and He can be reproduced by all scenarios if an AGB starting model with inclusion of overshoot is used for the post-AGB model sequence. An appendix is devoted to the numerical methods for models of proton capture nucleosynthesis in the He-flash convection zone during a thermal pulse.

We give a brief review on the properties of asymptotic giant branch models with overshoot. Then we describe new model calculations with overshoot. Initial masses are ranging from 1 to 3Msun and metallicities are Z=0.02, 0.01 and 0.001. Third dredge-up occurs efficiently for low masses and carbon stars are formed, with some at core masses as low as 0.58Msun. After the thermal pulse at which stars become C-rich the luminosities are in the range of the observed C-star luminosity function during the whole interpulse phase and for all C-star models. The dredge-up evolution depends mainly on the core mass at the first thermal pulse and on the metallicity. The Z=0.001 models of the 2 and 3Msun sequence become C-rich almost instantaneously after the onset of the first thermal pulses. For the 2Msun case the C/O ratio initially exceeds 4. During following dredge-up episodes the C/O ratio decreases.

The concept of overshoot has already been considered for numerous cases in stellar evolution calculations. We explore the consequences of overshoot at the convection zone which forms during the He-flash (thermal pulse) in AGB stars. We find dramatic changes for the abundances within the intershell region as well as for the mechanism of the 3rd dredge-up. That means that both the predicted evolution of structural quantities as well as the chemical evolution at the surface will be different, if overshoot is considered. We also present evidence for the presence of overshoot during the He-flash from detailed model calculations of the post-AGB phase and the comparison with observations. The new evolutionary models show that a very late thermal pulse during the post-AGB evolution can bring the intershell material up to the surface with minimal modification due to convective H-burning. The good agreement of the surface abundances of these models with observed surface abundances of H-deficient post-AGB stars is interpreted as a strong support for the presence of overshoot during the thermal pulses of AGB stars.

We calculated a grid of evolutionary models for white dwarfs with helium cores (He-WDs) and investigated the occurrence of hydrogen-shell flashes due to unstable hydrogen burning via CNO cycling. Our calculations show that such thermal instabilities are restricted to a certain mass range (M=0.21...0.30Msun), consistent with earlier studies. Models within this mass range undergo the more hydrogen shell flashes the less massive they are. This is caused by the strong dependence of the envelope mass on the white dwarf core mass. The maximum luminosities from hydrogen burning during the flashes are of the order of 10^5 Lsun. Because of the development of a pulse-driven convection zone whose upper boundary temporarily reaches the surface layers, the envelope's hydrogen content decreases by Delta(X)=0.06 per flash. Our study further shows that an additional high mass-loss episode during a flash-driven Roche lobe overflow to the white dwarf's companion does not affect the final cooling behaviour of the models. Independent of hydrogen shell flashes the evolution along the final white dwarf cooling branch is determined by hydrogen burning via pp-reactions down to effective temperatures as low as 8000 K.

We present an evolutionary sequence of a low mass star from the Asymptotic Giant Branch (AGB) through its post-AGB stage, during which its surface chemical composition changes from hydrogen-rich to strongly hydrogen-deficient as consequence of a very late thermal pulse, following the so-called born-again scenario. The internal structure and abundance changes during this pulse are computed with a numerical method which allows the physically consistent calculation of stellar layers where thermonuclear and mixing time scale are comparable - a situation which occurs when the helium flash driven convection zone extends to the hydrogen-rich surface layers during the pulse peak. The final surface mass fractions are [He/C/O]=[0.38/0.36/0.22], where the high oxygen abundance is due to diffusive overshoot employed during the AGB evolution. These models are the first to achieve general agreement with the surface abundance pattern observed in hydrogen-deficient post-AGB stars - e.g. the PG 1159 stars or the WR-type central stars of planetary nebulae -, confirming the born-again scenario with a physically consistent calculation and supporting the occurrence of convective overshooting in thermally pulsing AGB stars.

We present the first evolutionary models of intermediate mass stars up to their thermal pulses which include effects of rotation on the stellar structure as well as rotationally induced mixing of chemical species and angular momentum. We find a significant angular momentum transport from the core to the hydrogen-rich envelope and obtain a white dwarf rotation rate comparable to current observational upper limits of 50 km/s. Large angular momentum gradients at the bottom of the convective envelope and the tip of the pulse driven convective shell are shown to produce marked chemical mixing between the proton-rich and the 12C-rich layers during the so called third dredge-up. This leads to a subsequent production of 13C which is followed by neutron production through 13C(alpha,n) in radiative layers in between thermal pulses. Although uncertainties in the efficiency of rotational mixing processes persist, we conclude that rotation is capable of producing a 13C-rich layer as required for the occurrence of the s-process in TP-AGB stars.

We present seven evolutionary tracks for low-mass white dwarfs with helium cores, ranging in mass from 0.179 to 0.414 Msol. We generated the pre-white dwarf models from a 1 Msol sequence extending up to the tip of its red-giant branch by applying high mass-loss rates at appropriate positions, and we followed their evolution across the Hertzsprung-Russell diagram and down the cooling path. We discuss the internal structures and cooling properties of these new models and compare them with those of recently published models for low-mass white dwarfs which are based on simplified initial configurations. We also demonstrate that our new models seem to remove the apparent discrepancies between the characteristic ages of millisecond pulsars and the cooling ages of their white dwarf companions.

We investigate the influence of convective overshoot on stellar evolution models of the thermal pulse AGB phase with M_ZAMS = 3 Msol. An exponential diffusive overshoot algorithm is applied to all convective boundaries during all evolutionary stages. We demonstrate that overshooting at the bottom of the pulse-driven convective zone, which forms in the intershell during the He-shell flash, leads to more efficient third dredge-up. Some overshoot at the bottom of the convective envelope removes the He-H discontinuity, which would otherwise prohibit the occurrence of the third dredge-up for this stellar mass. However, no correlation between the amount of envelope overshoot and the efficiency of the third dredge-up has been found. Increasingly efficient third dredge-up eventually leads to a carbon star model. Due to the partial mixing efficiency in the overshoot region a C13-pocket can form after the third dredge-up event which may be crucial for n-capture nucleosynthesis.

We investigate the validity of the core mass - luminosity relation (CMLR), originally described by Paczynski (1970), for asymptotic giant branch stars under the presence of third dredge-up events. We find, that models with efficient third dredge-up with less massive cores than those associated with hot bottom burning (Bloecker and Schoenberner 1991) do not obey the linear CMLR. Complete evolutionary calculations of thermal pulse stellar models which consider overshoot according to an exponential diffusive algorithm show systematically larger third dredge-up for lower core masses (0.55 Msol < M_H < 0.8 Msol) than any other existing models. We present and discuss the luminosity evolution of these models.

We present a grid of evolutionary tracks for low-mass white dwarfs with helium cores in the mass range from 0.179 to 0.414 Msol. The lower mass limit is well-suited for comparison with white dwarf companions of millisecond pulsars. The tracks are based on a 1 Msol model sequence extending from the pre-main sequence stage up to the tip of the red-giant branch. Applying large mass loss rates at appropriate positions forced the models to move off the giant branch. The further evolution was then followed across the Hertzsprung-Russell diagram and down the cooling branch. At maximum effective temperature the envelope masses above the helium cores increase from 0.6 to 5.4 x 10^{-3} Msol for decreasing mass. We carefully checked for the occurrence of thermal instabilities of the hydrogen shell by adjusting the computational time steps accordingly. Hydrogen flashes have been found to take place only in the mass interval 0.21 < M/Msol < 0.3. The models show that hydrogen shell burning contributes significantly to the luminosity budget of white dwarfs with helium cores. For very low masses the hydrogen shell luminosity remains to be dominant even down to effective temperatures well below 10000K. Accordingly, the corresponding cooling ages are significantly larger than those gained from model calculations which neglect nuclear burning or the white dwarf progenitor evolution. Using the atmospheric parameters of the white dwarf in the PSR J1012+5307 system we determined a mass of M=0.19 +/- 0.02 Msol and a cooling age of 6 +/- 1 Gyr, in good agreement with the spin-down age, 7 Gyr, of the pulsar.

Based on radiation hydrodynamics modeling of stellar convection zones, a diffusion scheme has been devised describing the downward penetration of convective motions beyond the Schwarzschild boundary (overshoot) into the radiative interior. This scheme of exponential diffusive overshoot has already been successfully applied to AGB stars. Here we present an application to the Sun in order to determine the time scale and depth extent of this additional mixing, i.e. diffusive overshoot at the base of the convective envelope. We calculated the associated destruction of lithium during the evolution towards and on the main-sequence. We found that the slow-mixing processes induced by the diffusive overshoot may lead to a substantial depletion of lithium during the Sun's main-sequence evolution.

The focus of this study is on the treatment of those stellar regions immediately adjacent to convective zones. The results of hydrodynamical simulations by Freytag et al. (1996, A&A313,497) show that the motion of convective elements extends well beyond the boundary of the convectively unstable region. We have applied their parametrized description of the corresponding velocities to the treatment of overshoot in stellar evolution calculations up to the AGB (Pop.I, M_zams=3M_sun).