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

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

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

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

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

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

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

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

It has been suggested that the anomalous Na, Mg and Al observed in Globular Cluster Red Giant stars could be the result of a thermally unstable hydrogen shell. Currently accepted reaction rates indicate that temperatures of approximately 70-75 million K are required to produce the observed enhancements in Na and Al along with depletions in Mg. The work presented here attempts to model the H shell instability by a simple mechanism of altering the energy production in the region of the H shell. We show that even extreme cases only give rise to small intermittent temperature increases that have minimal affect on the surface abundances. Full evolutionary modelling incorporating this technique simply accelerates the evolution of the RGB phase producing the same surface abundances as other models but at an earlier time. We conclude that, unless hydrogen shell instabilities manifest themselves quite differently, they are unlikely to lead to the required temperatures and alternative explanations of the abundance anomalies are more promising., Comment: 16 pages, accepted for publication in MNRAS

The Asymptotic Giant Branch (AGB) phase is very short but its importance is seen in its nucleosynthesis. A revolution in stellar modelling has taken place in the last 20 years, inspired jointly by this rich nucleosynthesis and partly by new data. For example, the isotopic data coming from pre-solar grains (see this volume) forces theorists to include species that were previously ignored, species which are energetically of no importance (i.e. they play no role in determining the stellar structure) but which can be used to constrain the models. Nucleosynthesis is now important as a tracer of temperature and mixing, and not simply a by-product of energy generation. But along with these advances come more quantitative demands. It is now increasingly important to know what is known well and what is less sure. This is the goal of this paper., Comment: 9 pages, review presented at Ringberg Workshop on "Astronomy with Radioactivites", to appear in New Astronomy Reviews

We investigate the gravonuclear instabilities reported by Bono et al. (1997a,b) during the onset of helium-shell burning at the end of the horizontal-branch (HB) phase. These instabilities are characterized by relaxation oscillations within the helium shell which lead to loops in the evolutionary tracks. We find the occurrence of these instabilities depends critically on how the breathing pulses are suppressed near the end of the HB phase. If they are suppressed by omitting the gravitational energy term in the stellar structure equations, then the helium profile within the core at the end of the HB phase will contain a broad region of varying helium abundance. The helium-burning shell which forms in this region is too thick to be unstable, and gravonuclear instabilities do not occur. If, on the other hand, the breathing pulses are suppressed by prohibiting any increase in the central helium abundance, then the final helium profile can exhibit a large discontinuity at the edge of the helium-exhausted core. The helium shell which forms just exterior to this discontinuity is then much thinner and can be thermally unstable. Even in this case, however, the gravonuclear instabilities disappear as soon as the nuclear burning broadens the helium shell into its characteristic S-shape. We conclude that the gravonuclear instabilities found by Bono et al. are a consequence of the ad hoc procedure used to suppress the breathing pulses., Comment: 6 pages, 5 figures, to appear in "The Galactic Halo: From Globular Clusters to Field Stars", Proceedings of the 35th Liege International Astrophysical Colloquium

We summarise the main properties of Asymptotic Giant Branch stars, including their structure, evolution and nucleosynthesis. The main physical mechanisms are outlined, as are the uncertainties. In keeping with the multi-disciplinary nature of this meeting, this paper is designed for those who are not experts in stellar structure., Comment: 7 pages, 3 figures, Invited Review presented at Nuclei in the Cosmos V, Volos, July 1998. For non-experts in stellar strucure! For related papers please see http://www.maths.monash.edu.au/~johnl/

We present a review of the main phases of stellar evolution with particular emphasis on the nucleosynthesis and mixing mechanisms in low- and intermediate-mass stars. In addition to explicit studies of the effects of the first, second and third dredge-up, we also discuss cool bottom processing and hot bottom burning., Comment: 30 pages, latex, 18 figures, uses style files aipproc.cls aipproc.sty epsf.sty ; to be published in (refereed) conference proceedings "Astrophysical Implications of the Laboratory Study of Presolar Materials", ed. T. Bernatowitz and E. Zinner (AIP: Sunnyside, NY), in press; also available at http://www.maths.monash.edu.au/~boothroy/