Search

Here we present the discovery of 895 s-process-rich candidates from 454,180 giant stars observed by the Large Sky Area Multi-Object Fibre Spectroscopic Telescope (LAMOST) using a data-driven approach. This sample constitutes the largest number of s-process enhanced stars ever discovered. Our sample includes 187 s-process-rich candidates that are enhanced in both barium and strontium, 49 stars with significant barium enhancement only and 659 stars that show only a strontium enhancement. Most of the stars in our sample are in the range of effective temperature and log g typical of red giant branch (RGB) populations, which is consistent with our observational selection bias towards finding RGB stars. We estimate that only a small fraction (0.5 per cent) of binary configurations are favourable for s-process enriched stars. The majority of our s-process-rich candidates (95 per cent) show strong carbon enhancements, whereas only five candidates (less than 3 percent) show evidence of sodium enhancement. Our kinematic analysis reveals that 97 percent of our sample are disc stars, with the other 3 percent showing velocities consistent with the Galactic halo. The scaleheight of the disc is estimated to be zh=0.634kpc, comparable with values in the literature. A comparison with the yields from asymptotic giant branch (AGB) models suggests that the main neutron source responsible for the Ba and Sr enhancements is the 13C(alpha,n)16O reaction. We conclude that the s-process-rich candidates may have received their overabundances via mass transfer from a previous ABG companion with an initial mass in the range 1-3Msun., Comment: 9 pages, 6 figures

Chemical tagging seeks to identify unique star formation sites from present-day stellar abundances. Previous techniques have treated each abundance dimension as being statistically independent, despite theoretical expectations that many elements can be produced by more than one nucleosynthetic process. In this work we introduce a data-driven model of nucleosynthesis where a set of latent factors (e.g., nucleosynthetic yields) contribute to all stars with different scores, and clustering (e.g., chemical tagging) is modelled by a mixture of multivariate Gaussians in a lower-dimensional latent space. We use an exact method to simultaneously estimate the factor scores for each star, the partial assignment of each star to each cluster, and the latent factors common to all stars, even in the presence of missing data entries. We use an information-theoretic Bayesian principle to estimate the number of latent factors and clusters. Using the second Galah data release we find that six latent factors are preferred to explain N = 2,566 stars with 17 chemical abundances. We identify the rapid- and slow-neutron capture processes, as well as latent factors consistent with Fe-peak and \alpha-element production, and another where K and Zn dominate. When we consider N ~ 160,000 stars with missing abundances we find another 7 factors, as well as 16 components in latent space. Despite these components showing separation in chemistry that is explained through different yield contributions, none show significant structure in their positions or motions. We argue that more data, and joint priors on cluster membership that are constrained by dynamical models, are necessary to realise chemical tagging at a galactic-scale. We release software that allows for model parameters to be optimised in seconds given a fixed number of latent factors, components, and $10^7$ abundance measurements., Comment: Accepted to ApJ

Theoretical models of stellar evolution predict that most of the lithium inside a star is destroyed as the star becomes a red giant. However, observations reveal that about 1% of red giants are peculiarly rich in lithium, often exceeding the amount in the interstellar medium or predicted from the Big Bang. With only about 150 lithium-rich giants discovered in the past four decades, and no distinguishing properties other than lithium enhancement, the origin of lithium-rich giant stars is one of the oldest problems in stellar astrophysics. Here we report the discovery of 2,330 low-mass (1 to 3$\,M_\odot$) lithium-rich giant stars, which we argue are consistent with internal lithium production that is driven by tidal spin-up by a binary companion. Our sample reveals that most lithium-rich giants have helium-burning cores ($80^{+7}_{-6}\%$), and that the frequency of lithium-rich giants rises with increasing stellar metallicity. We find that while planet accretion may explain some lithium-rich giants, it cannot account for the majority that have helium-burning cores. We rule out most other proposed explanations as the primary mechanism for lithium-rich giants, including all stages related to single star evolution. Our analysis shows that giants remain lithium-rich for only about two million years. A prediction from this lithium depletion timescale is that most lithium-rich giants with a helium-burning core have a binary companion., Comment: ApJ accepted. Catalog of 2,330 lithium-rich giants available at https://zenodo.org/record/3515777

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., Comment: 14 pages including 3 figures

Abridged: Getting a better understanding of the evolution and nucleosynthetic yields of the most metal-poor stars (appr. Z<=10^-5) is critical because they are part of the big picture of the history of the primitive Universe. Yet many of the remaining unknowns of stellar evolution lie in the birth, life, and death of these objects. We review stellar evolution of intermediate-mass (IMS) Z<=10-5 models existing in the literature, with a focus on the problem of their final fates. The depth and efficiency of mixing episodes are critical to determine the mass limits for the formation of electron-capture supernovae, but our knowledge of these phenomena is not complete because they are strongly affected by the choice of input physics. We also consider the alternative SNI1/2 channel to form SNe out of the most metal-poor IMS. In this case, it is critical to understand the thermally-pulsing AGB evolution until the late stages. Efficient second dredge-up and, later, third dredge-up episodes could be able to pollute stellar envelopes enough for the stars to undergo thermal pulses in a way very similar to that of higher initial Z objects. Inefficient 2nd and/or 3rd dredge-up may leave an almost pristine envelope, unable to sustain strong stellar winds. This may allow the H-exhausted core to grow to M_Ch before the envelope is lost, and thus let the star explode as a SNI1/2. After reviewing the information available on these two possible channels for the formation of SNe, we discuss existing nucleosynthetic yields of stars of metallicity Z<=10^-5, and present an example of nucleosynthetic calculations for a thermally-pulsing Super-AGB star of Z=10^-5. We compare theoretical predictions with observations of the lowest [Fe/H] objects detected. The review closes by discussing current open questions as well as possible fruitful avenues for future research., Comment: 36 pages, 13 figures, 3 tables, PASA review, accepted

Extremely metal-poor (EMP) stars are an integral piece in the puzzle that is the early Universe, and although anomolous subclasses of EMP stars such as carbon-enhanced metal-poor (CEMP) stars are well-studied, they make up less than half of all EMP stars with [Fe/H] $\sim -3.0$. The amount of carbon depletion occurring on the red giant branch (carbon offset) is used to determine the evolutionary status of EMP stars, and this offset will differ between CEMP and normal EMP stars. The depletion mechanism employed in stellar models (from which carbon offfsets are derived) is very important, however the only widely available carbon offsets in the literature are derived from stellar models using a thermohaline mixing mechanism that cannot simultaneously match carbon and lithium abundances to observations for a single diffusion coeffcient. Our stellar evolution models utilise a modified thermohaline mixing model that can match carbon and lithium in the metal-poor globular cluster NGC 6397. We compare our models to the bulk of the EMP star sample at [Fe/H] $= -3$ and show that our modified models follow the trend of the observations and deplete less carbon compared to the standard thermohaline mixing theory. We conclude that stellar models that employ the standard thermohaline mixing formalism overestimate carbon offsets and hence CEMP star frequencies, particularly at metallicities where carbon-normal stars dominate the EMP star population., Comment: 8 pages, 4 figures, 1 table

Stars with unusual elemental abundances offer clues about rare astrophysical events or nucleosynthetic pathways. Stars with significantly depleted magnesium and enhanced potassium ([Mg/Fe] < -0.5; [K/Fe] > 1) have to date only been found in the massive globular cluster NGC 2419 and, to a lesser extent, NGC 2808. The origin of this abundance signature remains unknown, as does the reason for its apparent exclusivity to these two globular clusters. Here we present 112 field stars, identified from 454,180 LAMOST giants, that show significantly enhanced [K/Fe] and possibly depleted [Mg/Fe] abundance ratios. Our sample spans a wide range of metallicities (-1.5 < [Fe/H] < 0.3), yet none show abundance ratios of [K/Fe] or [Mg/Fe] that are as extreme as those observed in NGC 2419. If confirmed, the identified sample of stars represents evidence that the nucleosynthetic process producing the anomalous abundances ratios of [K/Fe] and [Mg/Fe] probably occurs at a wide range of metallicities. This would suggest that pollution scenarios that are limited to early epochs (such as Population III supernovae) are an unlikely explanation, although they cannot be ruled out entirely. This sample is expected to help guide modelling attempts to explain the origin of the Mg-K abundance signature.

A series of astronomical observations obtained over the period 1986 to 2018 supports the idea that life is a cosmic rather than a purely terrestrial or planetary phenomenon. These include (1) the detection of biologically relevant molecules in interstellar clouds and in comets, (2) mid-infrared spectra of interstellar grains and the dust from comets, (3) a diverse set of data from comets including the Rosetta mission showing consistency with biology and (4) the frequency of Earth-like or habitable planets in the Galaxy. We argue that the conjunction of all the available data suggests the operation of cometary biology and interstellar panspermia rather than the much weaker hypothesis of comets being only the source of the chemical building blocks of life. We conclude with specific predictions on the properties expected of extra-terrestrial life if it is discovered on Enceladus, Europa or beyond. A radically different biochemistry elsewhere can be considered as a falsification of the theory of interstellar panspermia., Comment: 13 pages, 5 figures, submitted to PASA

Theoretical predictions for the core helium burning phase of stellar evolution are highly sensitive to the uncertain treatment of mixing at convective boundaries. In the last few years, interest in constraining the uncertain structure of their deep interiors has been renewed by insights from asteroseismology. Recently, Spruit (2015) proposed a limit for the rate of growth of helium-burning convective cores based on the higher buoyancy of material ingested from outside the convective core. In this paper we test the implications of such a limit for stellar models with a range of initial mass and metallicity. We find that the constraint on mixing beyond the Schwarzschild boundary has a significant effect on the evolution late in core helium burning, when core breathing pulses occur and the ingestion rate of helium is fastest. Ordinarily, core breathing pulses prolong the core helium burning lifetime to such an extent that models are at odds with observations of globular cluster populations. Across a wide range of initial stellar masses ($0.83 \leq M/\text{M}_\odot \leq 5$), applying the Spruit constraint reduces the core helium burning lifetime because core breathing pulses are either avoided or their number and severity reduced. The constraint suggested by Spruit therefore helps to resolve significant discrepancies between observations and theoretical predictions. Specifically, we find improved agreement for $R_2$, the observed ratio of asymptotic giant branch to horizontal branch stars in globular clusters; the luminosity difference between these two groups; and in asteroseismology, the mixed-mode period spacing detected in red clump stars in the \textit{Kepler} field., Comment: Accepted for publication in MNRAS; 11 pages, 6 figures

Thermohaline mixing is a favoured mechanism for the so-called "extra mixing" on the red giant branch of low-mass stars. The mixing is triggered by the molecular weight inversion created above the hydrogen shell during first dredge-up when helium-3 burns via 3He(3He,2p)4He. The standard 1D diffusive mixing scheme cannot simultaneously match carbon and lithium abundances to NGC 6397 red giants. We investigate two modifications to the standard scheme: (1) an advective two stream mixing algorithm, and (2) modifications to the standard 1D thermohaline mixing formalism. We cannot simultaneously match carbon and lithium abundances using our two stream mixing approach. However we develop a modified diffusive scheme with an explicit temperature dependence that can simultaneously fit carbon and lithium abundances to NGC 6397 stars. Our modified diffusive scheme induces mixing that is faster than the standard theory predicts in the hotter part of the thermohaline region and mixing that is slower in the cooler part. Our results infer that the extra mixing mechanism needs further investigation and more observations are required, particularly for stars in different clusters spanning a range in metallicity., Comment: 15 pages, 15 figures

We review the lives, deaths and nucleosynthetic signatures of intermediate mass stars in the range approximately 6.5-12 Msun, which form super-AGB stars near the end of their lives. We examine the critical mass boundaries both between different types of massive white dwarfs (CO, CO-Ne, ONe) and between white dwarfs and supernovae and discuss the relative fraction of super-AGB stars that end life as either an ONe white dwarf or as a neutron star (or an ONeFe white dwarf), after undergoing an electron capture supernova. We also discuss the contribution of the other potential single-star channels to electron-capture supernovae, that of the failed massive stars. We describe the factors that influence these different final fates and mass limits, such as composition, the efficiency of convection, rotation, nuclear reaction rates, mass loss rates, and third dredge-up efficiency. We stress the importance of the binary evolution channels for producing electron-capture supernovae. We discuss recent nucleosynthesis calculations and elemental yield results and present a new set of s-process heavy element yield predictions. We assess the contribution from super-AGB star nucleosynthesis in a Galactic perspective, and consider the (super-)AGB scenario in the context of the multiple stellar populations seen in globular clusters. A brief summary of recent works on dust production is included. Lastly we conclude with a discussion of the observational constraints and potential future advances for study into these stars on the low mass/high mass star boundary., Comment: 28 pages, 11 figures. Invited review for Publications of the Astronomical Society of Australia, to be published in special issue on "Electron Capture Supernovae". Submitted

We summarise the evolution and nucleosynthesis in AGB and Super-AGB stars. We then examine the major sources of uncertainty, especially mass-loss., Comment: 8 pages, no figures. Invited review presented at The 11th Pacific Rim Conference on Stellar Astrophysics "Physics and Chemistry of the Late Stages of Stellar Evolution"

Globular clusters are of paramount importance for testing theories of stellar evolution and early galaxy formation. Strong evidence for multiple populations of stars in globular clusters derives from observed abundance anomalies. A puzzling example is the recently detected Mg-K anticorrelation in NGC 2419. We perform Monte Carlo nuclear reaction network calculations to constrain the temperature-density conditions that gave rise to the elemental abundances observed in this elusive cluster. We find a correlation between stellar temperature and density values that provide a satisfactory match between simulated and observed abundances in NGC 2419 for all relevant elements (Mg, Si, K, Ca, Sc, Ti, and V). Except at the highest densities ($\rho \gtrsim 10^8$~g/cm$^3$), the acceptable conditions range from $\approx$ $100$~MK at $\approx$ $10^8$~g/cm$^3$ to $\approx$ $200$~MK at $\approx$ $10^{-4}$~g/cm$^3$. This result accounts for uncertainties in nuclear reaction rates and variations in the assumed initial composition. We review hydrogen burning sites and find that low-mass stars, AGB stars, massive stars, or supermassive stars cannot account for the observed abundance anomalies in NGC 2419. Super-AGB stars could be viable candidates for the polluter stars if stellar model parameters can be fine-tuned to produce higher temperatures. Novae, either involving CO or ONe white dwarfs, could be interesting polluter candidates, but a current lack of low-metallicity nova models precludes firmer conclusions. We also discuss if additional constraints for the first-generation polluters can be obtained by future measurements of oxygen, or by evolving models of second-generation low-mass stars with a non-canonical initial composition., Comment: 24 pages, 9 Figures

The treatment of convective boundaries during core helium burning is a fundamental problem in stellar evolution calculations. In Paper~I we showed that new asteroseismic observations of these stars imply they have either very large convective cores or semiconvection/partially mixed zones that trap g-modes. We probe this mixing by inferring the relative lifetimes of asymptotic giant branch (AGB) and horizontal branch (HB) from $R_2$, the observed ratio of these stars in recent HST photometry of 48 Galactic globular clusters. Our new determinations of $R_2$ are more self-consistent than those of previous studies and our overall calculation of $R_2 = 0.117 \pm 0.005$ is the most statistically robust now available. We also establish that the luminosity difference between the HB and the AGB clump is $\Delta \log{L}_\text{HB}^\text{AGB} = 0.455 \pm 0.012$. Our results accord with earlier findings that standard models predict a lower $R_2$ than is observed. We demonstrate that the dominant sources of uncertainty in models are the prescription for mixing and the stochastic effects that can result from its numerical treatment. The luminosity probability density functions that we derive from observations feature a sharp peak near the AGB clump. This constitutes a strong new argument against core breathing pulses, which broaden the predicted width of the peak. We conclude that the two mixing schemes that can match the asteroseismology are capable of matching globular cluster observations, but only if (i) core breathing pulses are avoided in models with a semiconvection/partially mixed zone, or (ii) that models with large convective cores have a particular depth of mixing beneath the Schwarzschild boundary during subsequent early-AGB `gravonuclear' convection., Comment: 22 pages, 14 figures. Accepted for publication in MNRAS

There is now strong evidence that some stars have been born with He mass fractions as high as $Y \approx 0.40$ (e.g., in $\omega$ Centauri). However, the advanced evolution, chemical yields, and final fates of He-rich stars are largely unexplored. We investigate the consequences of He-enhancement on the evolution and nucleosynthesis of intermediate-mass asymptotic giant branch (AGB) models of 3, 4, 5, and 6 M$_\odot$ with a metallicity of $Z = 0.0006$ ([Fe/H] $\approx -1.4$). We compare models with He-enhanced compositions ($Y=0.30, 0.35, 0.40$) to those with primordial He ($Y=0.24$). We find that the minimum initial mass for C burning and super-AGB stars with CO(Ne) or ONe cores decreases from above our highest mass of 6 M$_\odot$ to $\sim$ 4-5 M$_\odot$ with $Y=0.40$. We also model the production of trans-Fe elements via the slow neutron-capture process (s-process). He-enhancement substantially reduces the third dredge-up efficiency and the stellar yields of s-process elements (e.g., 90% less Ba for 6 M$_\odot$, $Y=0.40$). An exception occurs for 3 M$_\odot$, where the near-doubling in the number of thermal pulses with $Y=0.40$ leads to $\sim$ 50% higher yields of Ba-peak elements and Pb if the $^{13}$C neutron source is included. However, the thinner intershell and increased temperatures at the base of the convective envelope with $Y=0.40$ probably inhibit the $^{13}$C neutron source at this mass. Future chemical evolution models with our yields might explain the evolution of s-process elements among He-rich stars in $\omega$ Centauri., Comment: 21 pages, 16 figures, accepted for publication by MNRAS. Stellar yields included as online data tables

The detection of mixed oscillation modes offers a unique insight into the internal structure of core helium burning (CHeB) stars. The stellar structure during CHeB is very uncertain because the growth of the convective core, and/or the development of a semiconvection zone, is critically dependent on the treatment of convective boundaries. In this study we calculate a suite of stellar structure models and their non-radial pulsations to investigate why the predicted asymptotic g-mode $\ell = 1$ period spacing $\Delta\Pi_1$ is systematically lower than is inferred from Kepler field stars. We find that only models with large convective cores, such as those calculated with our newly proposed "maximal-overshoot" scheme, can match the average $\Delta\Pi_1$ reported. However, we also find another possible solution that is related to the method used to determine $\Delta\Pi_1$: mode trapping can raise the observationally inferred $\Delta\Pi_1$ well above its true value. Even after accounting for these two proposed resolutions to the discrepancy in average $\Delta\Pi_1$, models still predict more CHeB stars with low $\Delta\Pi_1$ ($ < 270$ s) than are observed. We establish two possible remedies for this: i) there may be a difficulty in determining $\Delta\Pi_1$ for early CHeB stars (when $\Delta\Pi_1$ is lowest) because of the effect that the sharp composition profile at the hydrogen burning shell has on the pulsations, or ii) the mass of the helium core at the flash is higher than predicted. Our conclusions highlight the need for the reporting of selection effects in asteroseismic population studies in order to safely use this information to constrain stellar evolution theory., Comment: 24 pages. 24 figures. Published in MNRAS

We conduct a series of comparisons between spectroscopic and photometric observations of globular clusters and stellar models to examine their predictive power. Data from medium-to-high resolution spectroscopic surveys of lithium allow us to investigate first dredge-up and extra mixing in two clusters well separated in metallicity. Abundances at first dredge-up are satisfactorily reproduced but there is preliminary evidence to suggest that the models overestimate the luminosity at which the surface composition first changes in the lowest-metallicity system. Our models also begin extra mixing at luminosities that are too high, demonstrating a significant discrepancy with observations at low metallicity. We model the abundance changes during extra mixing as a thermohaline process and determine that the usual diffusive form of this mechanism cannot simultaneously reproduce both the carbon and lithium observations. Hubble Space Telescope photometry provides turnoff and bump magnitudes in a large number of globular clusters and offers the opportunity to better test stellar modelling as function of metallicity. We directly compare the predicted main-sequence turn-off and bump magnitudes as well as the distance-independent parameter $\Delta M_V ~^{\rm{MSTO}}_{\rm{bump}}$. We require 15 Gyr isochrones to match the main-sequence turn-off magnitude in some clusters and cannot match the bump in low-metallicity systems. Changes to the distance modulus, metallicity scale and bolometric corrections may impact on the direct comparisons but $\Delta M_V ~^{\rm{MSTO}}_{\rm{bump}}$, which is also underestimated from the models, can only be improved through changes to the input physics. Overshooting at the base of the convective envelope with an efficiency that is metallicity dependent is required to reproduce the empirically determined value of $\Delta M_V ~^{\rm{MSTO}}_{\rm{bump}}$., Comment: 20 pages, 11 Figures, 4 Tables, Accepted for publication in MNRAS

We examine the physical basis for algorithms to replace mixing-length theory (MLT) in stellar evolutionary computations. Our 321D procedure is based on numerical solutions of the Navier-Stokes equations. These implicit large eddy simulations (ILES) are three-dimensional (3D), time-dependent, and turbulent, including the Kolmogorov cascade. We use the Reynolds-averaged Navier-Stokes (RANS) formulation to make concise the 3D simulation data, and use the 3D simulations to give closure for the RANS equations. We further analyze this data set with a simple analytical model, which is non-local and time-dependent, and which contains both MLT and the Lorenz convective roll as particular subsets of solutions. A characteristic length (the damping length) again emerges in the simulations; it is determined by an observed balance between (1) the large-scale driving, and (2) small-scale damping. The nature of mixing and convective boundaries is analyzed, including dynamic, thermal and compositional effects, and compared to a simple model. We find that (1) braking regions (boundary layers in which mixing occurs) automatically appear {\it beyond} the edges of convection as defined by the Schwarzschild criterion, (2) dynamic (non-local) terms imply a non-zero turbulent kinetic energy flux (unlike MLT), (3) the effects of composition gradients on flow can be comparable to thermal effects, and (4) convective boundaries in neutrino-cooled stages differ in nature from those in photon-cooled stages (different P\'eclet numbers). The algorithms are based upon ILES solutions to the Navier-Stokes equations, so that, unlike MLT, they do not require any calibration to astronomical systems in order to predict stellar properties. Implications for solar abundances, helioseismology, asteroseismology, nucleosynthesis yields, supernova progenitors and core collapse are indicated., Comment: 22 pages, 4 figures, 2 tables; significantly re-written, critique of Pasetto, et al. model added, accepted for publication by ApJ

We have serendipitously identified the first lithium-rich giant star located close to the red giant branch bump in a globular cluster. Through intermediate-resolution FLAMES spectra we derived a lithium abundance of A(Li)=2.55 (assuming local thermodynamical equilibrium), which is extremely high considering the star's evolutionary stage. Kinematic and photometric analysis confirm the object as a member of the globular cluster NGC 362. This is the fourth Li-rich giant discovered in a globular cluster but the only one known to exist at a luminosity close to the bump magnitude. The three previous detections are clearly more evolved, located close to, or beyond the tip of their red giant branch. Our observations are able to discard the accretion of planets/brown dwarfs, as well as an enhanced mass-loss mechanism as a formation channel for this rare object. Whilst the star sits just above the cluster bump luminosity, its temperature places it towards the blue side of the giant branch in the colour-magnitude diagram. We require further dedicated observations to unambiguously identify the star as a red giant: we are currently unable to confirm whether Li production has occurred at the bump of the luminosity function or if the star is on the pre zero-age horizontal branch. The latter scenario provides the opportunity for the star to have synthesised Li rapidly during the core helium flash or gradually during its red giant branch ascent via some extra mixing process., Comment: Accepted for publication in The Astrophysical Journal Letters

In recent years much interest has been shown in the process of thermohaline mixing in red giants. In low and intermediate mass stars this mechanism first activates at the position of the bump in the luminosity function, and has been identified as a likely candidate for driving the slow mixing inferred to occur in these stars. One particularly important consequence of this process, which is driven by a molecular weight inversion, is the destruction of lithium. We show that the degree of lithium destruction, or in some cases production, is extremely sensitive to the numerical details of the stellar models. Within the standard 1D diffusion approximation to thermohaline mixing, we find that dfferent evolution codes, with their default numerical schemes, can produce lithium abundances that differ from one another by many orders of magnitude. This disagreement is worse for faster mixing. We perform experiments with four independent stellar evolution codes, and derive conditions for the spatial and temporal resolution required for a converged numerical solution. The results are extremely sensitive to the timesteps used. We find that predicted lithium abundances published in the literature until now should be treated with caution., Comment: 19 pages, 18 figures, accepted for publication by MNRAS; revised version has corrected Fig 17

We explore the final fates of massive intermediate-mass stars by computing detailed stellar models from the zero age main sequence until near the end of the thermally pulsing phase. These super-AGB and massive AGB star models are in the mass range between 5.0 and 10.0 Msun for metallicities spanning the range Z=0.02-0.0001. We probe the mass limits M_up, M_n and M_mass, the minimum masses for the onset of carbon burning, the formation of a neutron star, and the iron core-collapse supernovae respectively, to constrain the white dwarf/electron-capture supernova boundary. We provide a theoretical initial to final mass relation for the massive and ultra-massive white dwarfs and specify the mass range for the occurrence of hybrid CO(Ne) white dwarfs. We predict electron-capture supernova (EC-SN) rates for lower metallicities which are significantly lower than existing values from parametric studies in the literature. We conclude the EC-SN channel (for single stars and with the critical assumption being the choice of mass-loss rate) is very narrow in initial mass, at most approximately 0.2 Msun. This implies that between ~ 2-5 per cent of all gravitational collapse supernova are EC-SNe in the metallicity range Z=0.02 to 0.0001. With our choice for mass-loss prescription and computed core growth rates we find, within our metallicity range, that CO cores cannot grow sufficiently massive to undergo a Type 1.5 SN explosion., Comment: 15 pages, 7 figures, accepted for publication in MNRAS

Among the short-lived radioactive nuclei inferred to be present in the early solar system via meteoritic analyses, there are several heavier than iron whose stellar origin has been poorly understood. In particular, the abundances inferred for 182Hf (half-life = 8.9 million years) and 129I (half-life = 15.7 million years) are in disagreement with each other if both nuclei are produced by the rapid neutron-capture process. Here, we demonstrate that contrary to previous assumption, the slow neutron-capture process in asymptotic giant branch stars produces 182Hf. This has allowed us to date the last rapid and slow neutron-capture events that contaminated the solar system material at roughly 100 million years and 30 million years, respectively, before the formation of the Sun., Comment: 24 pages, 5 figures, published on Science Vol. 345, pp. 650-653

Convergent lines of evidence suggest that globular clusters host multiple stellar populations. It appears that they experience at least two episodes of star formation whereby a fraction of first-generation stars contribute astrated ejecta to form the second generation(s). To identify the polluting progenitors we require distinguishing chemical signatures such as that provided by lithium. Theoretical models predict that lithium can be synthesised in AGB stars, whereas no net Li production is expected from other candidates. It has been shown that in order to reproduce the abundance pattern found in M4, Li production must occur within the polluters, favouring the AGB scenario. Here we present Li and Al abundances for a large sample of RGB stars in M12 and M5. These clusters have a very similar metallicity, whilst demonstrating differences in several cluster properties. Our results indicate that the first-generation and second-generation stars share the same Li content in M12; we recover an abundance pattern similar to that observed in M4. In M5 we find a higher degree of complexity and a simple dilution model fails in reproducing the majority of the stellar population. In both clusters we require Li production across the different stellar generations, but production seems to have occurred to different extents. We suggest that such a difference might be related to the cluster mass with the Li production being more efficient in less-massive clusters. This is the first time a statistically significant correlation between the Li spread within a GC and its luminosity has been demonstrated. Finally, although Li-producing polluters are required to account for the observed pattern, other mechanisms, such as MS depletion, might have played a role in contributing to the Li internal variation, though at relatively low level., Comment: Accepted for publication in The Astrophysical Journal. 15 pages, 14 figures

The chemical evolution of the Universe is governed by the chemical yields from stars, which in turn is determined primarily by the initial stellar mass. Even stars as low as 0.9Msun can, at low metallicity, contribute to the chemical evolution of elements. Stars less massive than about 10Msun experience recurrent mixing events that can significantly change the surface composition of the envelope, with observed enrichments in carbon, nitrogen, fluorine, and heavy elements synthesized by the slow neutron capture process (the s-process). Low and intermediate mass stars release their nucleosynthesis products through stellar outflows or winds, in contrast to massive stars that explode as core-collapse supernovae. Here we review the stellar evolution and nucleosynthesis for single stars up to ~10Msun from the main sequence through to the tip of the asymptotic giant branch (AGB). We include a discussion of the main uncertainties that affect theoretical calculations and review the latest observational data, which are used to constrain uncertain details of the stellar models. We finish with a review of the stellar yields available for stars less massive than about 10Msun and discuss efforts by various groups to address these issues and provide homogeneous yields for low and intermediate-mass stars covering a broad range of metallicities., Comment: Invited Review accepted by Publications of the Astronomical Society of Australia; 70 pages

We present a new grid of stellar models and nucleosynthetic yields for super-AGB stars with metallicities Z=0.001 and 0.0001, applicable for use within galactic chemical evolution models. Contrary to more metal rich stars where hot bottom burning is the main driver of the surface composition, in these lower metallicity models the effect of third dredge-up and corrosive second dredge-up also have a strong impact on the yields. These metal-poor and very metal-poor super-AGB stars create large amounts of He4, C13 and N14, as well as the heavy magnesium isotopes Mg25 and Mg26. There is a transition in yield trends at metallicity Z approximately 0.001, below which we find positive yields of C12, O16, N15, Al27 and Si28, which is not the case for higher metallicities. We explore the large uncertainties derived from wind prescriptions in super-AGB stars, finding approximately 2 orders of magnitude difference in yields of Ne22, Na23, Mg24,25,26, Al27 and our s-process proxy isotope g. We find inclusion of variable composition low temperature molecular opacities is only critical for super-AGB stars of metallicities below Z approximately 0.001. We analyze our results, and those in the literature, to address the question: Are super-AGB stars the polluters responsible for extreme population in the globular cluster NGC 2808? Our results, as well as those from previous studies, seem unable to satisfactorily match the extreme population in this globular cluster., Comment: 19 pages, 9 figures, accepted for publication in MNRAS

The vital importance of composition-dependent low-temperature opacity in low-mass (M < 3Msun) asymptotic giant branch (AGB) stellar models of metallicity Z > 0.001 has recently been demonstrated (e.g. Marigo 2002; Ventura & Marigo 2010). Its significance to more metal-poor, intermediate mass (M > 2.5Msun) models has yet to be investigated. We show that its inclusion in lower-metallicity models ([Fe/H] < -2) is essential, and that there exists no threshold metallicity below which composition-dependent molecular opacity may be neglected. We find it to be crucial in all intermediate-mass models investigated ([Fe/H] < -2 and 2.5 < M/Msun < 5), because of the evolution of the surface chemistry, including the orders of magnitude increase in the abundance of molecule-forming species. Its effect on these models mirrors that previously reported for higher-metallicity models - increase in radius, decrease in Teff, faster mass loss, shorter thermally pulsing AGB lifetime, reduced enrichment in third dredge-up products (by a factor of three to ten), and an increase in the mass limit for hot bottom burning. We show that the evolution of low-metallicity models with composition-dependent low-temperature opacity is relatively independent of initial metal abundance because its contribution to the opacity is far outweighed by changes due to dredge-up. Our results imply a significant reduction in the expected number of nitrogen-enhanced metal-poor stars, which may help explain their observed paucity. We note that these findings are partially a product of the macrophysics adopted in our models, in particular the Vassiliadis & Wood (1993) mass loss rate which is strongly dependent on radius., Comment: 13 pages, 13 figures, 2 tables; accepted for publication in The Astrophysical Journal

We have computed detailed evolution and nucleosynthesis models for super and massive AGB stars over the mass range 6.5-9.0 Msun in divisions of 0.5 Msun with metallicities Z=0.02, 0.008 and 0.004. These calculations, in which we find third dredge-up and hot bottom burning, fill the gap between existing low and intermediate-mass AGB star models and high mass star models that become supernovae. For the considered metallicities, the composition of the yields is largely dominated by the thermodynamic conditions at the base of the convective envelope rather than by the pollution arising from third dredge up. We investigate the effects of various uncertainties, related to the mass-loss rate, mixing length parameter, and the treatment of evolution after the envelope instability that develops near the end of the (Super)AGB phase. Varying these parameters alter the yields mainly because of their impact on the amount of third dredge up enrichment, and to a lesser extent on the hot bottom burning conditions. Our models produce significant amounts of He4, Li7 (depending on the mass-loss formulation) C13, N14, O17, Na23, Mg25, as well the radioactive isotope Al26 in agreement with previous investigation. In addition our results show enrichment of Ne22, Mg26 and Fe60, as well as a substantial increase in our proxy neutron capture species representing all species heavier than iron. These stars may provide important contributions to the Galaxy's inventory of the heavier Mg isotopes, N14, Li7 and Al27., Comment: 22 pages, 15 figures, accepted for publication in MNRAS

The asymptotic giant branch (AGB) phase is the final stage of nuclear burning for low-mass stars. Although Milky Way globular clusters are now known to harbour (at least) two generations of stars they still provide relatively homogeneous samples of stars that are used to constrain stellar evolution theory. It is predicted by stellar models that the majority of cluster stars with masses around the current turn-off mass (that is, the mass of the stars that are currently leaving the main sequence phase) will evolve through the AGB phase. Here we report that all of the second-generation stars in the globular cluster NGC 6752 -- 70 per cent of the cluster population -- fail to reach the AGB phase. Through spectroscopic abundance measurements, we found that every AGB star in our sample has a low sodium abundance, indicating that they are exclusively first-generation stars. This implies that many clusters cannot reliably be used for star counts to test stellar evolution timescales if the AGB population is included. We have no clear explanation for this observation., Comment: Published in Nature (online 29 May 2013, hard copy 13 June), 12 pages, 3 figures + supplementary information section

Globular clusters display significant variations in their light-element content, pointing to the existence of a second stellar generation formed from the ejecta of an earlier generation. The nature of these internal polluters is still a matter of debate: the two most popular scenarios indicate intermediate-mass AGB stars and fast rotating massive stars. Abundances determination for some key elements can help distinguish between these competitor candidates. We present in this paper Y abundances for a sample of 103 red giant branch stars in NGC 6121. Within measurement errors, we find that the [Y/Fe] is constant in this cluster contrary to a recent suggestion. For a subsample of six stars we also find [Rb/Fe] to be constant, consistent with previous studies showing no variation in other s-process elements. We also present a new set of stellar yields for intermediate-mass AGB stellar models of 5 and 6 solar masses, including heavy element s-process abundances. The uncertainties on the mass-loss rate, the mixing-length parameter, and the nuclear reaction rates have a major impact on the stellar abundances. Within the IM-AGB pollution scenario, the constant abundance of heavy elements inside the cluster requires a marginal s-process efficiency in IM-AGB stars. Such a constrain could still be satisfied by the present models assuming a stronger mass-loss rate. The uncertainties mentioned above are limiting the predictive power of intermediate-mass AGB models. For these reasons, at the moment we are not able to clearly rule out their role as main polluters of the second population stars in globular clusters., Comment: Accepted for publication in MNRAS, 18 pages, 21 figures

The composition of the most primitive Solar System condensates, such as calcium-aluminum-rich inclusions (CAI) and micron-sized corundum grains, show that short-lived radionuclides (SLR), e.g., 26Al, were present in the early Solar System. Their abundances require a local origin, which however is far from being understood. We present for the first time the abundances of several SLR up to 60Fe predicted from stars with initial mass in the range roughly 7-11 Msun. These stars evolve through core H, He, and C burning. After core C burning they go through a "Super"-asymptotic giant branch (Super-AGB) phase, with the H and He shells activated alternately, episodic thermal pulses in the He shell, a very hot temperature at the base of the convective envelope (~ 10^8 K), and strong stellar winds driving the H-rich envelope into the surrounding interstellar medium. The final remnants of the evolution of Super-AGB stars are mostly O-Ne white dwarfs. Our Super-AGB models produce 26Al/27Al yield ratios ~ 0.02 - 0.26. These models can account for the canonical value of the 26Al/27Al ratio using dilutions with the Solar Nebula of the order of 1 part of Super-AGB mass per several 10^2 to several 10^3 of Solar Nebula mass, resulting in associated changes in the O composition in the range 3 permil to 20 permil. This is in agreement with observations of the O isotopic ratios in primitive Solar System condensated, which do not carry the signature of a stellar pollutor. The radionuclides 41Ca and 60Fe are produced by neutron captures in Super-AGB stars and their meteoritic abundances are also matched by some of our models, depending on the nuclear and stellar physics uncertainties as well as the meteoritic experimental data. We also expect and are currently investigating Super-AGB production of SLR heavier than iron, such as 107Pd., Comment: 23 pages, 5 figures, 6 tables, accepted for publication on Meteoritics and Planetary Science

The literature is rich in analysis and results related to thermally pulsing-asymptotic giant branch (TP-AGB) stars, but the problem of the instabilities that arise and cause the divergence of models during the late stages of their evolution is rarely addressed. We investigate the physical conditions, causes and consequences of the interruption in the calculations of massive AGB stars in the late thermally-pulsing AGB phase. We have thoroughly analysed the physical structure of a solar metallicity 8.5 solar mass star and described the physical conditions at the base of the convective envelope (BCE) just prior to divergence. We find that the local opacity maximum caused by M-shell electrons of Fe-group elements lead to the accumulation of an energy excess, to the departure of thermal equilibrium conditions at the base of the convective envelope and, eventually, to the divergence of the computed models. For the 8.5 solar mass case we present in this work the divergence occurs when the envelope mass is about 2 solar mass. The remaining envelope masses range between somewhat less than 1 and more than 2 solar mass for stars with initial masses between 7 and 10 solar mass and, therefore, our results are relevant for the evolution and yields of super-AGB stars. If the envelope is ejected as a consequence of the instability we are considering, the occurrence of electron-capture supernovae would be avoided at solar metallicity., Comment: 19 pages, 9 figures, accepted for A&A

It is now widely accepted that globular cluster red giant branch stars owe their strange abundance patterns to a combination of pollution from progenitor stars and in situ extra mixing. In this hybrid theory a first generation of stars imprint abundance patterns into the gas from which a second generation forms. The hybrid theory suggests that extra mixing is operating in both populations and we use the variation of [C/Fe] with luminosity to examine how efficient this mixing is. We investigate the observed red giant branches of M3, M13, M92, M15 and NGC 5466 as a means to test a theory of thermohaline mixing. The second parameter pair M3 and M13 are of intermediate metallicity and our models are able to account for the evolution of carbon along the RGB in both clusters. Although, in order to fit the most carbon-depleted main-sequence stars in M13 we require a model whose initial [C/Fe] abundance leads to a carbon abundance lower than is observed. Furthermore our results suggest that stars in M13 formed with some primary nitrogen (higher C+N+O than stars in M3). In the metal-poor regime only NGC 5466 can be tentatively explained by thermohaline mixing operating in multiple populations. We find thermohaline mixing unable to model the depletion of [C/Fe] with magnitude in M92 and M15. It appears as if extra mixing is occurring before the luminosity function bump in these clusters. To reconcile the data with the models would require first dredge-up to be deeper than found in extant models., Comment: 13 Pages, 3 figures. Accepted for publication in the Astrophysical Journal

We use the 3D stellar structure code DJEHUTY to model the ingestion of protons into the intershell convection zone of a 1 solar mass asymptotic giant branch star of metallicity Z=10^-4. We have run two simulations: a low resolution one of around 300,000 zones, and a high resolution one consisting of 2,000,000 zones. Both simulations have been evolved for about 4 hours of stellar time. We observe the existence of fast, downward flowing plumes that are able to transport hydrogen into close proximity to the helium burning shell before burning takes place. The intershell in the 3D model is richer in protons than the 1D model by several orders of magnitude and so we obtain substantially higher hydrogen-burning luminosities -- over 10^8 solar luminosities in the high resolution simulation -- than are found in the 1D model. Convective velocities in these simulations are over 10 times greater than the predictions of mixing length theory, though the 3D simulations have greater energy generation due to the enhanced hydrogen burning. We find no evidence of the convective zone splitting into two, though this could be as a result of insufficient spatial resolution or because the models have not been evolved for long enough. We suggest that the 1D mixing length theory and particularly the use of a diffusion algorithm for mixing do not give an accurate picture of these events. An advective mixing scheme may give a better representation of the transport processes seen in the 3D models., Comment: 24 pages, 10 figures. Accepted for publication in ApJ. Movies to accompany the manuscript can be found at http://www.mso.anu.edu.au/~rjs/djehutyPIE.html

We review the observational evidence for extra mixing in stars on the red giant branch (RGB) and discuss why thermohaline mixing is a strong candidate mechanism. We recall the simple phenomenological description of thermohaline mixing, and aspects of mixing in stars in general. We use observations of M3 to constrain the form of the thermohaline diffusion coefficient and any associated free parameters. This is done by matching [C/Fe] and [N/Fe] along the RGB of M3. After taking into account a presumed initial primordial bimodality of [C/Fe] in the CN-weak and CN-strong stars our thermohaline mixing models can explain the full spread of [C/Fe]. Thermohaline mixing can produce a significant change in [N/Fe] as a function of absolute magnitude on the RGB for initially CN-weak stars, but not for initially CN-strong stars, which have so much nitrogen to begin with that any extra mixing does not significantly affect the surface nitrogen composition., Comment: 33 Pages, 10 Figures. Accepted for publication in The Astrophysical Journal

There is an apparent dichotomy between the metal-poor ([Fe/H] <= -2) yet carbon-normal giants and their carbon-rich counterparts. The former undergo significant depletion of carbon on the red giant branch after they have undergone first dredge-up, whereas the latter do not appear to experience significant depletion. We investigate this in the context that the extra mixing occurs via the thermohaline instability that arises due to the burning of helium-3. We present the evolution of [C/Fe], [N/Fe] and carbon-12/carbon-13 for three models: a carbon-normal metal-poor star, and two stars that have accreted material from a 1.5 solar mass AGB companion, one having received 0.01 solar masses of material and the other having received 0.1 solar masses. We find the behaviour of the carbon-normal metal-poor stars is well reproduced by this mechanism. In addition, our models also show that the efficiency of carbon-depletion is significantly reduced in carbon-rich stars. This extra-mixing mechanism is able to reproduce the observed properties of both carbon-normal and carbon-rich stars., Comment: Accepted for publication in MNRAS. 7 pages, 8 figures

Abundances of C, N, and O are determined in four bright red giants that span the known abundance range for light (Na and Al) and s-process (Zr and La) elements in the globular cluster NGC 1851. The abundance sum C+N+O exhibits a range of 0.6 dex, a factor of 4, in contrast to other clusters in which no significant C+N+O spread is found. Such an abundance range offers support for the Cassisi et al. (2008) scenario in which the double subgiant branch populations are coeval but with different mixtures of C+N+O abundances. Further, the Na, Al, Zr, and La abundances are correlated with C+N+O, and therefore, NGC 1851 is the first cluster to provide strong support for the scenario in which AGB stars are responsible for the globular cluster light element abundance variations., Comment: Accepted for publication in ApJ Letters

We revisit the problem of why stars become red giants. We modify the physics of a standard stellar evolution code in order to determine what does and what does not contribute to a star becoming a red giant. In particular, we have run tests to try to separate the effects of changes in the mean molecular weight and in the energy generation. The implications for why stars become red giants are discussed. We find that while a change in the mean molecular weight is necessary (but not sufficient) for a 1 solar mass star to become a red giant, this is not the case in a star of 5 solar masses. It therefore seems that there may be more than one way to make a giant., Comment: 7 pages, 9 figures. Accepted for publication in PASA

We study the post-main sequence evolution of products of collisions between main sequence stars (blue stragglers), with particular interest paid to the horizontal branch and asymptotic giant branch phases. We found that the blue straggler progeny populate the colour-magnitude diagram slightly blueward of the red giant branch and between 0.2 and 1 magnitudes brighter than the horizontal branch. We also found that the lifetimes of collision products on the horizontal branch is consistent with the numbers of so-called "evolved blue straggler stars" (E-BSS) identified by various authors in a number of globular clusters, and is almost independent of mass or initial composition profile. The observed ratio of the number of E-BSS to blue stragglers points to a main sequence lifetime for blue stragglers of approximately 1-2 Gyr on average., Comment: accepted by the Astrophysical Journal

Recent work has proposed that a merger event between a red-giant and a He white dwarf may be responsible for the production of R-stars (Izzard et al, 2007). We investigate the proposed evolution and nucleosynthesis of such a model. We simulate the hypothesized late ignition of the core flash by increasing the neutrino losses until the ignition occurs sufficiently far from the centre that the subsequent evolution produces dredge-up of carbon to the extent that the post-flash object is a carbon star. Detailed nucleosynthesis is performed within this approximation, and we show that the overall properties are broadly consistent with the observations. Details will depend on the dynamics of the merger event., Comment: 8 pages, 6 figures, to be published in PASA

We report on an investigation into stellar evolution and nucleosynthesis in the low and extremely low metallicity regime, including models of stars with a pure Big Bang composition (i.e. Z=0). The metallicity range of the extremely metal poor (EMP) models we have calculated is -6.5 < [Fe/H] < -3.0, whilst our models are in the mass range 0.85 < M < 3.0 Msun. Many of the EMP and Z=0 models experience violent evolutionary episodes not seen at higher metallicities. We refer to these events as `Dual Flashes' since they are characterised by peaks in the hydrogen and helium burning luminosities occurring at roughly the same time. Some of the material processed by these events is later dredged up by the convective envelope, causing signifcant surface pollution. These events have been reported by previous studies, so our results reaffirm their occurrence -- at least in 1D stellar models. The novelty of this study is that we have calculated the entire evolution of the Z=0 and EMP models, from the ZAMS to the end of the TPAGB, including detailed nucleosynthesis. We have also calculated the nucleosynthetic yields, which will soon be available in electronic format. Although subject to many uncertainties these are, as far as we are aware, the only yields available in this mass and metallicity range. In this paper we briefy describe some of the results in the context of abundance observations of EMP halo stars. This work formed part of SWC's PhD thesis (completed in March 2007) and a series of subsequent papers will describe the results of the study in more detail., Comment: 5 pages, 6 figures. To appear in conference proceedings "First Stars III", ed. B. O'Shea, A. Heger & T.Abel

We present stellar yields calculated from detailed models of low and intermediate-mass asymptotic giant branch (AGB) stars. We evolve models with a range of mass from 1 to 6Msun, and initial metallicities from solar to 1/200th of the solar metallicity. Each model was evolved from the zero age main sequence to near the end of the thermally-pulsing AGB phase, and through all intermediate phases including the core He-flash for stars initially less massive than 2.5Msun. For each mass and metallicity, we provide tables containing structural details of the stellar models during the TP-AGB phase, and tables of the stellar yields for 74 species from hydrogen through to sulphur, and for a small number of iron-group nuclei. All tables are available for download. Our results have many applications including use in population synthesis studies and the chemical evolution of galaxies and stellar systems, and for comparison to the composition of AGB and post-AGB stars and planetary nebulae., Comment: 26 pages; to appear in Publications of the Astronomical Society of Australia (PASA); typos fixed in the text and in Tables 4 and 6

Three-dimensional stellar modeling has enabled us to identify a deep-mixing mechanism that must operate in all low mass giants. This mixing process is not optional, and is driven by a molecular weight inversion created by the 3He(3He,2p)4He reaction. In this paper we characterize the behavior of this mixing, and study its impact on the envelope abundances. It not only eliminates the problem of 3He overproduction, reconciling stellar and big bang nucleosynthesis with observations, but solves the discrepancy between observed and calculated CNO isotope ratios in low mass giants, a problem of more than 3 decades' standing. This mixing mechanism, which we call `$\delta\mu$-mixing', operates rapidly (relative to the nuclear timescale of overall evolution, ~ 10^8 yrs) once the hydrogen burning shell approaches the material homogenized by the surface convection zone. In agreement with observations, Pop I stars between 0.8 and 2.0$\Msun$ develop 12C/13C ratios of 14.5 +/- 1.5, while Pop II stars process the carbon to ratios of 4.0 +/- 0.5. In stars less than 1.25$\Msun$, this mechanism also destroys 90% to 95% of the 3He produced on the main sequence., Comment: Final accepted version (submitted to Astrophys J in Jan 2007...)

The early-R stars are carbon-rich K-type giants. They are enhanced in C12, C13 and N14, have approximately solar oxygen, magnesium isotopes, s-process and iron abundances, have the luminosity of core-helium burning stars, are not rapid rotators, are members of the Galactic thick disk and, most peculiarly of all, are all single stars. Conventional single-star stellar evolutionary models cannot explain such stars, but mergers in binary systems have been proposed to explain their origin. We have synthesized binary star populations to calculate the number of merged stars with helium cores which could be early-R stars. We find many possible evolutionary channels. The most common of which is the merger of a helium white dwarf with a hydrogen-burning red giant branch star during a common envelope phase followed by a helium flash in a rotating core which mixes carbon to the surface. All the channels together give ten times more early-R stars than we require to match recent Hipparcos observations - we discuss which channels are likely to be the true early-R stars and which are not. For the first time we have constructed a viable model of the early-R stars with which we can test some of our ideas regarding common envelope evolution in giants, stellar mergers, rotation, the helium flash and the origin of the early-R stars., Comment: 17 Pages, accepted for publication in A&A

We continue our study of the core helium flash using the three dimensional hydrodynamics code Djehuty. Continuing from earlier calculations, we now take relaxed 3D configurations and add various amounts of rotation. We find that rotation periods consistent with those observed in white dwarfs produce negligible changes in the structure and evolution of the core flash, at least for the very small timescales we have yet been able to investigate. There is no sign of any extra mixing due to the rotation. There is some inconclusive evidence for a slight change in the luminosity, at the 1% level., Comment: Accepted for publication in the (refereed) Proceedings of the Nuclei In The Cosmos IX meeting

Low-mass stars, ~1-2 solar masses, near the Main Sequence are efficient at producing He-3, which they mix into the convective envelope on the giant branch and should distribute into the Galaxy by way of envelope loss. This process is so efficient that it is difficult to reconcile the low observed cosmic abundance of He-3 with the predictions of both stellar and Big Bang nucleosynthesis. In this paper we find, by modeling a red giant with a fully three-dimensional hydrodynamic code and a full nucleosynthetic network, that mixing arises in the supposedly stable and radiative zone between the hydrogen-burning shell and the base of the convective envelope. This mixing is due to Rayleigh-Taylor instability within a zone just above the hydrogen-burning shell, where a nuclear reaction lowers the mean molecular weight slightly. Thus we are able to remove the threat that He-3 production in low-mass stars poses to the Big Bang nucleosynthesis of He-3., Comment: Accepted by Science, and available from Science Express online

We investigate the suggestion that there are stellar populations in some globular clusters with enhanced helium (Y from 0.28 to 0.40) compared to the primordial value. We assume that a previous generation of massive Asymptotic Giant Branch (AGB) stars have polluted the cluster. Two independent sets of AGB yields are used to follow the evolution of helium and CNO using a Salpeter initial mass function (IMF) and two top-heavy IMFs. In no case are we able to produce the postulated large Y ~ 0.35 without violating the observational constraint that the CNO content is nearly constant., Comment: accepted for publication in ApJ

Using a chemical evolution model we investigate the intriguing suggestion that there are populations of stars in some globular clusters (e.g. NGC 2808, omega Centauri) with enhanced levels of helium (Y from about 0.28 to 0.40) compared to the majority of the population that presumably have a primordial helium abundance. We assume that a previous generation of massive low-metallicity Asymptotic Giant Branch (AGB) stars has polluted the cluster gas via a slow stellar wind. We use two independent sets of AGB yields computed from detailed models to follow the evolution of helium, carbon, nitrogen and oxygen in the cluster gas using a Salpeter initial mass function (IMF) and a number of top-heavy IMFs. In no case were we able to fit the observational constraints, Y > 0.30 and C+N+O approximately constant. Depending on the shape of the IMF and the yields, we either obtained Y approximately greater than 0.30 and large increases in C+N+O or Y < 0.30 and C+N+O approximately constant. These results suggest that either AGB stars alone are not responsible for the large helium enrichment or that any dredge-up from this generation of stars was less than predicted by standard models., Comment: To appear in Proceedings of The VIII Torino Workshop on Nucleosynthesis in AGB stars: Constraints on AGB Nucleosynthesis from Observations, edited by C. Abia and I. Dominguez

We model the core helium flash in a low-mass red giant using Djehuty, a fully three-dimensional (3D) code. The 3D structures were generated from converged models obtained during the 1D evolutionary calculation of a 1$\Msun$ star. Independently of which starting point we adopted, we found that after some transient relaxation the 3D model settled down with a briskly convecting He-burning shell that was not very different from what the 1D model predicted., Comment: 32 pages, 13 figures, accepted by ApJ. Figures reduced from high-res eps to low-res gif for rapid download. Contact Lattanzio for a copy of the higher res images (included, inplace, in the manuscript)

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., Comment: resubmitted to Phys Rev C

We present nucleosynthesis calculations and the resulting 19F stellar yields for a large set of models with different masses and metallicity. We find that the production of fluorine depends on the temperature of the convective pulses, the amount of primary 12C mixed into the envelope by third dredge up and the extent of the partial mixing zone. Then we perform a detailed analysis of the reaction rates involved in the production of 19F and the effects of their uncertainties. We find that the major uncertainties are associated with the 14C(alpha,gamma)18O and the 19F(alpha,p)22Ne reaction rates. For these two reactions we present new estimates of the rates and their uncertainties. The importance of the partial mixing zone is reduced when using our estimate for the 14C(alpha,gamma)18O rate. Taking into account both the uncertainties related to the partial mixing zone and those related to nuclear reactions, the highest values of 19F enhancements observed in AGB stars are not matched by the models. This is a problem that will have to be revised by providing a better understanding of the formation and nucleosynthesis in the partial mixing zone, also in relation to reducing the uncertainties of the 14C(alpha,gamma)18O reaction rate. At the same time the possible effect of Cool Bottom Processing at the base of the convective envelope should be included in the computation of AGB nucleosynthesis. This process could in principle help matching the highest 19F abundances observed by decreasing the C/O ratio at the surface of the star, while leaving the 19F abundance unchanged., Comment: 40 pages, 8 figures, accepted for publication on the Astrophysical Journal