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.

We present a detailed near-infrared chemical abundance analysis of 10 red giant members of the Galactic open cluster NGC 752. High-resolution (R$\simeq$45000) near-infrared spectral data were gathered with the Immersion Grating Infrared Spectrograph (IGRINS), providing simultaneous coverage of the complete H and K bands. We derived the abundances of H-burning (C, N, O), $\alpha$ (Mg, Si, S, Ca), light odd-Z (Na, Al, P, K), Fe-group (Sc, Ti, Cr, Fe, Co, Ni) and neutron-capture (Ce, Nd, Yb) elements. We report the abundances of S, P, K, Ce, and Yb in NGC 752 for the first time. Our analysis yields solar metallicity and solar abundance ratios for almost all of the elements heavier than the CNO group in NGC 752. O and N abundances were measured from a number of OH and CN features in the $H$ band, and C abundances were determined mainly from CO molecular lines in the K band. High excitation \ion{C}{i} lines present in both near-infrared and optical spectra were also included in the C abundance determinations. Carbon isotopic ratios were derived from the R-branch band heads of first overtone (2-0) and (3$-$1) $^{12}$CO and (2-0) $^{13}$CO lines near 23440 \AA and (3-1) $^{13}$CO lines at about 23730 \AA. The CNO abundances and $^{12}$C/$^{13}$C ratios are all consistent with our giants having completed "first dredge-up" envelope mixing of CN-cyle products. We independently assessed NGC 752 stellar membership from Gaia astrometry, leading to a new color-magnitude diagram for this cluster. Applications of Victoria isochrones and MESA models to these data yield an updated NGC 752 cluster age (1.52 Gyr) and evolutionary stage indications for the program stars. The photometric evidence and spectroscopic light element abundances all suggest that the most, perhaps all of the program stars are members of the helium-burning red clump in this cluster.

We present an analysis of the relative age distribution of the Milky Way halo, based on samples of blue horizontal-branch (BHB) stars obtained from the Panoramic Survey Telescope and Rapid Response System and \textit{Galaxy Evolution Explorer} photometry, as well a Sloan Digital Sky Survey spectroscopic sample. A machine-learning approach to the selection of BHB stars is developed, using support vector classification, with which we produce chronographic age maps of the Milky Way halo out to 40\,kpc from the Galactic center. We identify a characteristic break in the relative age profiles of our BHB samples, corresponding to a Galactocentric radius of $R_{\rm{GC}} \sim 14$\,kpc. Within the break radius, we find an age gradient of $-63.4 \pm 8.2$ Myr kpc$^{-1}$, which is significantly steeper than obtained by previous studies that did not discern between the inner- and outer-halo regions. The gradient in the relative age profile and the break radius signatures persist after correcting for the influence of metallicity on our spectroscopic calibration sample. We conclude that neither are due to the previously recognized metallicity gradient in the halo, as one passes from the inner-halo to the outer-halo region. Our results are consistent with a dissipational formation of the inner-halo population, involving a few relatively massive progenitor satellites, such as those proposed to account for the assembly of \textit{Gaia}-Enceladus, which then merged with the inner halo of the Milky Way.

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.

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.

We present near-infrared spectroscopic analysis of 12 red giant members of the Galactic open cluster NGC 6940. High-resolution (R$\simeq$45000) and high signal-to-noise ratio (S/N > 100) near-infrared H and K band spectra were gathered with the Immersion Grating Infrared Spectrograph (IGRINS) on the 2.7m Smith Telescope at McDonald Observatory. We obtained abundances of H-burning (C, N, O), ${\alpha}$ (Mg, Si, S, Ca), light odd-Z (Na, Al, P, K), Fe-group (Sc, Ti, Cr, Fe, Co, Ni) and neutron-capture (Ce, Nd, Yb) elements. We report the abundances of S, P, K, Ce, and Yb in NGC 6940 for the first time. Many OH and CN features in the H band were used to obtain O and N abundances. C abundances were measured from four different features: CO molecular lines in the K band, high excitation C I lines present in both near-infrared and optical, CH and $C_2$ bands in the optical region. We have also determined $^{12}C/^{13}C$ ratios from the R-branch band heads of first overtone (2-0) and (3-1) $^{12}CO$ (2-0) $^{13}CO$ lines near 23440 $\overset{\lower.5em\circ}{\mathrm{A}}$ and (3-1) $^{13}CO$ lines at about 23730 $\overset{\lower.5em\circ}{\mathrm{A}}$. We have also investigated the HF feature at 23358.3 $\overset{\lower.5em\circ}{\mathrm{A}}$, finding solar fluorine abundances without ruling out a slight enhancement. For some elements (such as the ${\alpha}$ group), IGRINS data yield more internally self-consistent abundances. We also revisited the CMD of NGC 6940 by determining the most probable cluster members using Gaia DR2. Finally, we applied Victoria isochrones and MESA models in order to refine our estimates of the evolutionary stages of our targets.

We selected blue horizontal branch (BHB) stars within the expected distance range and sky position of the Hermus Stream from Data Release 10 of the Sloan Digital Sky Survey. We identify a moving group of $19$ BHB stars that are concentrated within two degrees of the Hermus Stream, between $10$ and $14$ kpc from the Sun. The concentration in velocity is inconsistent with a Gaussian distribution with 98% confidence (2.33 sigma).The stars in the moving group have line-of-sight velocities of v$_{\rm gsr} \sim 50$ km s$^{-1}$, a velocity dispersion of $\sigma_v\lesssim11$ km s$^{-1}$, a line-of-sight depth of $\sim 1$ kpc, and a metallicity of [Fe/H] $=-2.1 \pm 0.4$. The best-fit orbit has a perigalacticon of $\sim4$ kpc, apogalacticon of $\sim 17$ kpc, orbital period of $\sim247$ Myr, eccentricity $e = 0.62$, and inclination $i \sim75^\circ$ from $b=90^\circ$. The BHB stars in the stream are estimated to be $12$ Gyr old. An N-body simulation of a mass-follows-light ultrafaint dwarf galaxy with mass $10^6M_{\odot}$ and radius $40$ pc is consistent with the observed properties. The properties of the identified moving group of 19 BHB stars are close enough to those of the Hermus Stream (which is traced predominantly in turnoff stars) that we find it likely that they are associated. If that is the case, then our orbit fit would imply that there is no relationship between the Hermus and Phoenix streams, as previously proposed.

M55 (NGC 6809) and NGC 6362 are among the few globular clusters for which masses and radii have been derived to high precision for member binary stars. They also contain RR Lyrae variables which, together with their non-variable horizontal-branch (HB) populations, provide tight constraints on the cluster reddenings and distance moduli through fits of stellar models to their pulsational and evolutionary properties. Reliable (m-M)_V estimates yield M_V and M_bol values of comparable accuracy for binary stars because the V-band bolometric corrections applicable to them have no more than a weak dependence on effective temperature (T_eff) and [Fe/H]. Chemical abundances derived from the binary mass-M_V relations are independent of determinations based on their spectra. The temperatures of the binaries, which are calculated directly from their luminosities and the measured radii, completely rule out the low T_eff scale that has been determined for metal-deficient stars in some recent spectroscopic and interferometric studies. If [alpha/Fe] = 0.4 and [O/Fe] = 0.5 +/- 0.1, we find that M55 has (m-M)_V = 13.95 +/- 0.05, [Fe/H] = -1.85 +/- 0.1, and an age of 12.9 +/- 0.8 Gyr, whereas NGC 6362 has (m-M)_V = 14.56 +/- 0.05, [Fe/H] = -0.90 +/- 0.1, and an age of 12.4 +/- 0.8 Gyr. The HB of NGC 6362 shows clear evidence for multiple stellar populations. Constraints from the RR Lyrae standard candle and from local subdwarfs (with Gaia DR2 parallaxes) are briefly discussed.

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)

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.

We present a new set of horizontal-branch (HB) models computed with the MESA stellar evolution code. The models adopt $\alpha$-enhanced \cite{ags09} metals mixtures and include the gravitational settling of He. They are used in our HB population synthesis tool to generate theoretical distributions of HB stars in order to describe the multiple stellar populations in the globular clusters 47Tuc, M3, and M13. The observed HB in 47Tuc is reproduced very well by our simulations for [Fe/H] $= -0.70$ and [$\alpha$/Fe] $= +0.4$ if the initial helium mass fraction varies by $\Delta Y_0 \sim 0.03$ and approximately 21%, 37%, and 42% of the stars have $Y_0 = 0.257$, 0.270, and 0.287, respectively. These simulations yield $(m-M)_V = 13.27$, implying an age near 13.0 Gyr. In the case of M3 and M13, our synthetic HBs for [Fe/H] $= -1.55$ and [$\alpha$/Fe] $= 0.4$ match the observed ones quite well if M3 has $\Delta Y_0 \sim 0.01$ and $(m-M)_V = 15.02$, resulting in an age of 12.6 Gyr, whereas M13 has $\Delta Y_0 \sim 0.08$ and $(m-M)_V = 14.42$, implying an age of 12.9 Gyr. Mass loss during giant-branch evolution and $\Delta Y_0$ appear to be the primary second parameters for M3 and M13. New observations for 7 of the 9 known RR Lyrae in M13 are also reported. Surprisingly, periods predicted for the $c$-type variables tend to be too high (by up to $\sim 0.1$~d).

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 present a new, high-resolution chronographic (age) map of the Milky Way's halo, based on the inferred ages of ~130,000 field blue horizontal-branch (BHB) stars with photometry from the Sloan Digital Sky Survey. Our map exhibits a strong central concentration of BHB stars with ages greater than 12 Gyr, extending up to ~15 kpc from the Galactic center (reaching close to the solar vicinity), and a decrease in the mean ages of field stars with distance by 1-1.5 Gyr out to ~45-50 kpc, along with an apparent increase of the dispersion of stellar ages, and numerous known (and previously unknown) resolved over-densities and debris streams, including the Sagittarius Stream. These results agree with expectations from modern LambdaCDM cosmological simulations, and support the existence of a dual (inner/outer) halo system, punctuated by the presence of over-densities and debris streams that have not yet completely phase-space mixed.

Up-to-date isochrones, zero-age horizontal branch (ZAHB) loci, and evolutionary tracks for core He-burning stars are applied to the color-magnitude diagrams (CMDs) of M3, M15, and M92, focusing in particular on their RR Lyrae populations. Periods for ab- and c-type variables are calculated using the latest theoretical calibrations of log P_ab and log P_c as a function of luminosity, mass, effective temperature, and metallicity. Our models are generally able to reproduce the measured periods to well within the uncertainties implied by the stellar properties on which pulsation periods depend, as well as the mean periods and cluster-to-cluster differences in <P_ab> and <P_c>, on the assumption of well-supported values of E(B-V), (m-M)_V, and [Fe/H]. While many of the RR Lyrae in M3 lie close to the same ZAHB that fits the faintest HB stars at bluer or redder colors, the M92 variables are all significantly evolved stars from ZAHB locations on the blue side of the instability strip. M15 appears to contain a similar population of HB stars as M92, along with additional helium-enhanced populations not present in the latter which comprise most of its RR Lyrae stars. The large number of variables in M15 and the similarity of the observed values of <P_ab> and <P_c> in M15 and M92 can be explained by HB models that allow for variations in Y. Similar ages (~12.5 Gyr) are found for all three clusters, making them significantly younger than the field halo subgiant HD 140283. Our analysis suggests a preference for stellar models that take diffusive processes into account.

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.

In a pioneering effort, Preston et al. reported that the colors of blue horizontal-branch (BHB) stars in the halo of the Galaxy shift with distance, from regions near the Galactic center to about 12 kpc away, and interpreted this as a correlated variation in the ages of halo stars, from older to younger, spanning a range of a few Gyrs. We have applied this approach to a sample of some 4700 spectroscopically confirmed BHB stars selected from the Sloan Digital Sky Survey to produce the first "chronographic map" of the halo of the Galaxy. We demonstrate that the mean de-reddened g$-$r color, <(g$-$r)o>, increases outward in the Galaxy from $-$0.22 to $-$0.08 (over a color window spanning [$-$0.3:0.0]) from regions close to the Galactic center to ~40 kpc, independent of the metallicity of the stars. Models of the expected shift in the color of the field BHB stars based on modern stellar evolutionary codes confirm that this color gradient can be associated with an age difference of roughly 2-2.5 Gyrs, with the oldest stars concentrated in the central ~15 kpc of the Galaxy. Within this central region, the age difference spans a mean color range of about 0.05 mag (~0.8 Gyrs). Furthermore, we show that chronographic maps can be used to identify individual substructures, such as the Sagittarius Stream, and overdensities in the direction of Virgo and Monoceros, based on the observed contrast in their mean BHB colors with respect to the foreground/background field population.

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.

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.

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.

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)

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.

The accretion of hydrogen-rich material onto carbon-oxygen white dwarfs (CO WDs) is crucial for understanding type Ia supernova (SN Ia) from the single-degenerate model, but this process has not been well understood due to the numerical difficulties in treating H and He flashes during the accretion. For the CO WD masses from 0.5 to $1.378\,{M}_\odot$ and accretion rates in the range from $10^{-8}$ to $10^{-5}\,{M}_\odot\,\mbox{yr}^{-1}$, we simulated the accretion of solar-composition material onto CO WDs using the state-of-the-art stellar evolution code of {\sc MESA}. For comparison with the steady-state models (e.g \citet{nskh07}), we firstly ignored the contribution from nuclear burning to the luminosity when determining the Eddington accretion rate and found that the properties of H burning in our accreting CO WD models are similar to those from the steady-state models, except that the critical accretion rates at which the WDs turn into red giants or H-shell flashes occur on their surfaces are slightly higher than those from the steady-state models. However, the super-Eddington wind is triggered at much lower accretion rates, than previously thought, when the contribution of nuclear burning to the total luminosity is included. This super-Eddington wind naturally prevents the CO WDs with high accretion rates from becoming red giants, thus presenting an alternative to the optically thick wind proposed by \cite{hkn96}. Furthermore, the super-Eddington wind works in low-metallicity environments, which may explain SNe Ia observed at high redshifts.

We propose that the abundance anomalies of proton-capture elements in globular clusters, such as the C-N, Na-O, Mg-Al and Na-F anti-correlations, were produced by super-massive stars with M ~ 10,000 Msun. Such stars could form in the runaway collisions of massive stars that sank to the cluster center as a result of dynamical friction, or via the direct monolithic collapse of the low-metallicity gas cloud from which the cluster formed. To explain the observed abundance anomalies, we assume that the super-massive stars had lost significant parts of their initial masses when only a small mass fraction of hydrogen, Delta X ~ 0.15, was transformed into helium. We speculate that the required mass loss might be caused by the super-Eddington radiation continuum-driven stellar wind or by the diffusive mode of the Jeans instability.

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.

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.

These notes provide a tutorial for those who want to use the state-of-the-art stellar evolution code MESA and post-processing nucleosyntheis tools of NuGrid. As an example, an application of MESA and NuGrid tools for simulations of a nova outburst and associated nucleosynthesis occurring on a 1.3 Msun ONe white dwarf is presented.

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 show that the evolutionary track of a low-mass red giant should make an extended zigzag on the Hertzsprung-Russel diagram just after the bump luminosity, if fast internal rotation and enhanced extra mixing in the radiative zone bring the temperature gradient close to the adiabatic one. This can explain both the location and peculiar surface chemical composition of Li-rich K giants studied by Kumar, Reddy, & Lambert (2011). We also discuss a striking resemblance between the photometric and composition peculiarities of these stars and giant components of RS CVn binaries. We demonstrate that the observationally constrained values of the temperature gradient in the Li-rich K giants agree with the required rate of extra mixing only if the turbulence which is believed to be responsible for this extra mixing is highly anisotropic, with its associated transport coefficients in the horizontal direction strongly dominating over those in the vertical direction.

First results of our 3D numerical simulations of thermohaline convection driven by 3He burning in a low-mass RGB star at the bump luminosity are presented. They confirm our previous conclusion that this convection has a mixing rate which is a factor of 50 lower than the observationally constrained rate of RGB extra-mixing. It is also shown that the large-scale instabilities of salt-fingering mean field (those of the Boussinesq and advection-diffusion equations averaged over length and time scales of many salt fingers), which have been observed to increase the rate of oceanic thermohaline mixing up to one order of magnitude, do not enhance the RGB thermohaline mixing. We speculate on possible alternative solutions of the problem of RGB extra-mixing, among which the most promising one that is related to thermohaline mixing is going to take advantage of the shifting of salt-finger spectrum towards larger diameters by toroidal magnetic field.

Low-mass stars are known to experience extra-mixing in their radiative zones on the red-giant branch (RGB) above the bump luminosity. To determine if the salt-fingering transport of chemical composition driven by 3He burning is efficient enough to produce RGB extra-mixing, 2D numerical simulations of thermohaline convection for physical conditions corresponding to the RGB case have been carried out. We have found that the effective ratio of a salt-finger's length to its diameter a_eff < 0.5 is more than ten times smaller than the value needed to reproduce observations (a_obs > 7). On the other hand, using the thermohaline diffusion coefficient from linear stability analysis together with a=a_obs is able to describe the RGB extra-mixing at all metallicities so well that it is tempting to believe that it may represent the true mechanism. In view of these results, follow-up 3D numerical simulations of thermohaline convection for the RGB case are clearly needed.

The model of magnetic braking of solar rotation considered by Charbonneau & MacGregor (1993) has been modified so that it is able to reproduce for the first time the rotational evolution of both the fastest and slowest rotators among solar-type stars in open clusters of different ages, without coming into conflict with other observational constraints, such as the time evolution of the atmospheric Li abundance in solar twins and the thinness of the solar tachocline. This new model assumes that rotation-driven turbulent diffusion, which is thought to amplify the viscosity and magnetic diffusivity in stellar radiative zones, is strongly anisotropic with the horizontal components of the transport coefficients strongly dominating over those in the vertical direction. Also taken into account is the poloidal field decay that helps to confine the width of the tachocline at the solar age. The model's properties are investigated by numerically solving the azimuthal components of the coupled momentum and magnetic induction equations in two dimensions using a finite element method

We critically examine the constraints on internal angular momentum transport which can be inferred from the spin down of open cluster stars. The rotation distribution inferred from rotation velocities and periods are consistent for larger and more recent samples, but smaller samples of rotation periods appear biased relative to vsini studies. We therefore focus on whether the rotation period distributions observed in star forming regions can be evolved into the observed ones in the Pleiades, NGC2516, M34, M35, M37, and M50 with plausible assumptions about star-disk coupling and angular momentum loss from magnetized solar-like winds. Solid body models are consistent with the data for low mass fully convective stars but highly inconsistent for higher mass stars where the surface convection zone can decouple for angular momentum purposes from the radiative interior. The Tayler-Spruit magnetic angular momentum transport mechanism, commonly employed in models of high mass stars, predicts solid-body rotation on extremely short timescales and is therefore unlikely to operate in solar-type pre-MS and MS stars at the predicted rate. Models with core-envelope decoupling can explain the spin down of 1.0 and 0.8 solar mass slow rotators with characteristic coupling timescales of 55+-25 Myr and 175+-25 Myr respectively. The upper envelope of the rotation distribution is more strongly coupled than the lower envelope of the rotation distribution, in accord with theoretical predictions that the angular momentum transport timescale should be shorter for more rapidly rotating stars. Constraints imposed by the solar rotation curve are also discussed (Abridged)

We revise a magnetic buoyancy model that has recently been proposed as a mechanism for extra mixing in the radiative zones of low-mass red giants. The most important revision is our accounting of the heat exchange between rising magnetic flux rings and their surrounding medium. This increases the buoyant rising time by five orders of magnitude, therefore the number of magnetic flux rings participating in the mixing has to be increased correspondingly. On the other hand, our revised model takes advantage of the fact that the mean molecular weight of the rings formed in the vicinity of the hydrogen burning shell has been reduced by 3He burning. This increases their thermohaline buoyancy (hence, decreases the total ring number) considerably, making it equivalent to the pure magnetic buoyancy produced by a frozen-in toroidal field with B_phi ~ 10 MG. We emphasize that some toroidal field is still needed for the rings to remain cohesive while rising. Besides, this field prevents the horizontal turbulent diffusion from eroding the mu contrast between the rings and their surrounding medium. We propose that the necessary toroidal magnetic field is generated by differential rotation of the radiative zone, that stretches a pre-existing poloidal field around the rotation axis, and that magnetic flux rings are formed as a result of its buoyancy-related instability.

We examine the stability and observational consequences of mixing induced by 3He burning in the envelopes of first ascent red giants. We demonstrate that there are two unstable modes: a rapid, nearly adiabatic mode that we cannot identify with an underlying physical mechanism, and a slow, nearly radiative mode that can be identified with thermohaline convection. We present observational constraints that make the operation of the rapid mode unlikely to occur in real stars. Thermohaline convection turns out to be fast enough only if fluid elements have finger-like structures with a length to diameter ratio l/d > 10. We identify some potentially serious obstacles for thermohaline convection as the predominant mixing mechanism for giants. We show that rotation-induced horizontal turbulent diffusion may suppress the 3He-driven thermohaline convection. Another potentially serious problem for it is to explain observational evidence of enhanced extra mixing. The 3He exhaustion in stars approaching the red giant branch (RGB) tip should make the 3He mixing inefficient on the asymptotic giant branch (AGB). In spite of this, there are observational data indicating the presence of extra mixing in low-mass AGB stars similar to that operating on the RGB. Overmixing may also occur in carbon-enhanced metal-poor stars.

Internal gravity waves (IGWs) are naturally produced by convection in stellar envelopes, and they could be an important mechanism for transporting angular momentum in the radiative interiors of stars. Prior work has established that they could operate over a short enough time scale to explain the internal solar rotation as a function of depth. We demonstrate that the natural action of IGWs is to produce large scale oscillations in the solar rotation as a function of depth, which is in marked contrast to the nearly uniform rotation in the outer radiative envelope of the Sun. An additional angular momentum transport mechanism is therefore required, and neither molecular nor shear-induced turbulent viscosity is sufficient to smooth out the profile. Magnetic processes, such as the Tayler-Spruit dynamo, could flatten the rotation profile. We therefore conclude that IGWs must operate in conjunction with magnetic angular momentum transport processes if they operate at all. Furthermore, both classes of mechanisms must be inhibited to some degree by mean molecular weight gradients in order to explain the recent evidence for a rapidly rotating embedded core in the Sun.

We critically examine the constraints imposed by carbon-enhanced metal-poor (CEMP) stars on the mixing mechanisms that operate in red giants. CEMP stars are created when the surface layers of a metal-poor dwarf are enriched with He-burning products via mass transfer from an evolved donor. The difference between main sequence (MS) and red giant CEMP abundances can be used as a diagnostic of the timescale for the mixing of the processed material into stellar interiors on the MS. Abundance trends with luminosity among red giant CEMP stars test theories of canonical extra mixing for low mass giants with a high bulk metallicity. We find a significant dilution in CN enrichment in giant CEMP stars relative to their MS precursors, and take this as evidence that thermohaline mixing induced by mean molecular weight inversions is ineffective in CEMP stars. This contradicts models that rely on efficient thermohaline mixing induced by small mu gradients in red giants, because such models would predict that MS CEMP stars with large mu inversions would be homogenized on a very short timescale. The data do not rule out slower MS thermohaline mixing comparable to previously published estimates. We also find that canonical extra mixing is strongly suppressed in CEMP giants relative to stars with the same iron abundance. A likely cause is the shift in the location of non-equilibrium CN processing relative to the steep mu gradient in the hydrogen burning shell, which also occurs in solar metallicity RGB stars. Implications for the mass accreted by CEMP stars and the mechanism for canonical extra mixing are discussed.

Angular momentum transport by internal magnetic fields is an important ingredient for stellar interior models. In this paper we critically examine the basic heuristic assumptions in the model of the Tayler-Spruit dynamo, which describes how a pinch-type instability of a toroidal magnetic field in differentially rotating stellar radiative zones may result in large-scale fluid motion. We agree with prior published work both on the existence of the instability and its nearly horizontal geometry for perturbations. However, the approximations in the original Acheson dispersion relation are valid only for small length scales, and we disagree that the dispersion relation can be extrapolated to horizontal length scales of order the radius of the star. We contend that dynamical effects, in particular angular momentum conservation, limit the maximum horizontal length scale. We therefore present transport coefficients for chemical mixing and angular momentum redistribution by magnetic torques that are significantly different from previous published values. The new magnetic viscosity is reduced by 2 to 3 orders of magnitude compared to the old one, and we find that magnetic angular momentum transport by this mechanism is very sensitive to gradients in the mean molecular weight. The revised coefficients are more compatible with empirical constraints on the timescale of core-envelope coupling in young stars than the previous ones. However, solar models including only this mechanism possess a rapidly rotating core, in contradiction with helioseismic data. (abridged)

We show that enhanced extra mixing in low-mass red giants can result in a fluorine abundance that is correlated with abundance variations of other elements participating in H burning, such as C, N, O and Na. This finding is used to explain the fluorine abundance variations recently found in bright red giants of the globular cluster M4.

Stellar models which incorporate simple diffusion or shear induced mixing are used to describe canonical extra mixing in low mass red giants of low and solar metallicity. These models are able to simultaneously explain the observed Li and CN abundance changes along upper red giant branch (RGB) in field low-metallicity stars and match photometry, rotation and carbon isotopic ratios for stars in the old open cluster M67. The shear mixing model requires that main sequence (MS) progenitors of upper RGB stars possessed rapidly rotating radiative cores and that specific angular momentum was conserved in each of their mass shells during their evolution. We surmise that solar-type stars will not experience canonical extra mixing on the RGB because their more efficient MS spin-down resulted in solid-body rotation, as revealed by helioseismological data for the Sun. Thus, RGB stars in the old, high metallicity cluster NGC 6791 should show no evidence for mixing in their carbon isotopic ratios. We develop the idea that canonical extra mixing in a giant component of a binary system may be switched to its enhanced mode with much faster and somewhat deeper mixing as a result of the giant's tidal spin-up. This scenario can explain photometric and composition peculiarities of RS CVn binaries. The tidally enforced enhanced extra mixing might contribute to the star-to-star abundance variations of O, Na and Al in globular clusters. This idea may be tested with observations of carbon isotopic ratios and CN abundances in RS CVn binaries.

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 modify our combined ("deep mixing" plus "primordial") scenario explaining the star-to-star abundance variations in globular cluster red giants in such a way that it confirms with new experimental data: (i) the new and better constrained (NACRE) thermonuclear reaction rates and (ii) a discovery of the O-Na anticorrelation in stars below the main-sequence turn-off in the cluster NGC6752. For the latter we propose that some main-sequence dwarfs in globular clusters accrete material lost by red giant primary components of close binaries during a common envelope event. As a consequence of the new reaction rates, we are drawn to the conclusion that the anomalies in [Al/Fe] in globular cluster red giants are in fact manifestations of 26Al (instead of the stable isotope 27Al) abundance variations produced by deep mixing.

The recent discovery of low-mass red giants with enhanced atmospherical Li abundance in galactic low-metallicity stellar clusters adds to the mysterious phenomenon of the Li-rich giants. We propose a combined scenario for the Li-enrichment: engulfing a giant planet (or brown dwarf) by a red giant (external source) activates inside the giant the "7Be-mechanism" producing Li internally. This episodical Li-production can happen at any time on the red giant branch and is naturally followed by Li-depletion as is observed. Limitations of our scenario are discussed as well.

We have calculated the evolution of low metallicity red giant stars under the assumption of deep mixing between the convective envelope and the hydrogen burning shell. We find that the extent of the observed abundance anomalies, and in particular the universal O-Na anticorrelation, can be totally explained by mixing which does not lead to significant helium enrichment of the envelope. On the other hand, models with extremely deep mixing and strong helium enrichment predict anomalies of sodium and oxygen, which are much larger than the observed ones. This latter result depends solely on the nucleosynthesis inside the hydrogen burning shell, but not on the details of the mixing descriptions. These, however, influence the evolution of surface abundances with brightness, which we compare with the limited observational material available. Our models allow, nevertheless, to infer details on the depth and speed of the mixing process in several clusters. Models with strong helium enrichment evolve to high luminosities and show an increased mass loss. However, under peculiar assumptions, red giants reach very high luminosities even without extreme helium mixing. Due to the consequently increased mass loss, such models could be candidates for blue horizontal branch stars, and, at the same time, would be consistent with the observed abundance anomalies.

We present the evolution of low mass and low metallicity RGB models computed under the assumption of deep mixing between the convective envelope and the hydrogen burning shell. We discuss the helium enrichment of the envelope which is necessary or allowed to achieve or to keep consistency with the observations, and conclude on the possible connection between chemical anomalies in red giants and the horizontal branch morphology in globular clusters.

Star-to-star abundance variations of C, N, O, Na and Al in globular-cluster red giants have been recently supplemented by the finding that [Mg/Fe] is depleted in stars with extremely large [Al/Fe] (Shetrone 1996a). This and other new spectroscopic results allow us to test current models of stellar evolution and nucleosynthesis, as well as those of the formation and chemical enrichment of globular clusters. In an effort to explain self-consistently these observations we have considered two possibilities: (1) a deep mixing scenario which assumes that in red giants some kind of (extra)mixing transports products of nuclear reactions from the hydrogen burning shell (HBS) to the base of the convective envelope; and (2) a combination of primordial and deep mixing scenarios. It is shown that (1) cannot account for the anticorrelation of [Mg/Fe] vs. [Al/Fe] without additional ad hoc assumptions, among which we identify a strong but still undetected low energy resonance in the reaction 24Mg(p,gamma)25Al, and episodical increases of the HBS temperature up to the value T approx. 74x10^6 K. In (2) intermediate mass AGB stars are assumed to produce the decreased 24Mg and increased 25Mg initial abundances in some globular-cluster low mass stars and Al is synthesized at the expense of 25Mg in the HBS and transported to the surface of the red giant by extramixing. We discuss advantages and deficiencies of both scenarios and propose some observational tests.

We follow the evolution of 3He in stars of mass 0.8 -- 10 M (solar units) for Pop. I and II compositions from the main sequence until advanced stages on the AGB. Under standard assumptions we confirm earlier results of more restricted investigations that low-mass stars up to 5 M are net producers of 3He. We show that the inclusion of additional mixing due to diffusion beneath the convective envelope simultaneously leads to observed carbon isotope anomalies observed in globular cluster Red Giants and to a strong reduction of 3He, such that stars exhibiting such anomalies will have destroyed 3He contrary to the standard picture.