Search

Asteroseismic masses of globular cluster (GC) stars are invaluable to investigate stellar evolution. Previously, only two GCs have been seismically studied. We present new detections of solar-like oscillations in the clusters M9 and M19, focusing on two key areas: stellar mass loss and GC multiple populations. Using K2 photometry, we detect solar-like oscillations in stars on the red giant branch and early asymptotic giant branch. We measure an integrated mass-loss for M9 of $0.16\pm0.02$(rand)$\pm0.03$(sys)$M_{\odot}$ and M19 of $0.33\pm0.03$(rand)$^{+0.09}_{-0.07}$(sys)$M_{\odot}$. Comparing these to the mass-loss estimates from previous seismically studied clusters, we derive a preliminary relationship between stellar mass-loss and metallicity for Type I GCs. We find that the mass-loss for M19 -- a Type II GC -- is significantly larger, suggesting Type II clusters follow a different mass-loss-metallicity trend. We also examine the mass distributions in each evolutionary phase for evidence of a bimodality that could indicate mass differences between sub-populations. While no clear bimodality is observed, there is tentative evidence suggesting the presence of two mass populations. Classification through spectroscopic abundances into the sub-populations is needed to verify these findings. This study reinforces that asteroseismology of GC stars provides an excellent testbed for studying stellar evolution. However, to advance the field we need high-quality photometry of more GCs, a goal that could be realised with the upcoming Roman Telescope., Comment: 18 pages, 15 figures

Stellar mergers and accretion events have been crucial in shaping the evolution of the Milky Way (MW). These events have been dynamically identified and chemically characterised using red giants and main-sequence stars. RR Lyrae (RRL) variables can play a crucial role in tracing the early formation of the MW since they are ubiquitous, old (t$\ge$10 Gyr) low-mass stars and accurate distance indicators. We exploited Data Release 3 of the GALAH survey to identify 78 field RRLs suitable for chemical analysis. Using synthetic spectra calculations, we determined atmospheric parameters and abundances of Fe, Mg, Ca, Y, and Ba. Most of our stars exhibit halo-like chemical compositions, with an iron peak around [Fe/H]$\approx -$1.40, and enhanced Ca and Mg content. Notably, we discovered a metal-rich tail, with [Fe/H] values ranging from $-$1 to approximately solar metallicity. This sub-group includes almost ~1/4 of the sample, it is characterised by thin disc kinematics and displays sub-solar $\alpha$-element abundances, marginally consistent with the majority of the MW stars. Surprisingly, they differ distinctly from typical MW disc stars in terms of the s-process elements Y and Ba. We took advantage of similar data available in the literature and built a total sample of 535 field RRLs for which we estimated kinematical and dynamical properties. We found that metal-rich RRLs (1/3 of the sample) likely represent an old component of the MW thin disc. We also detected RRLs with retrograde orbits and provided preliminary associations with the Gaia-Sausage-Enceladus, Helmi, Sequoia, Sagittarius, and Thamnos stellar streams., Comment: Accepted for publication in MNRAS. 29 pages, 20 figures

Carbon-deficient red giants (CDGs) are a peculiar class of stars that have eluded explanation for decades. We aim to better characterise CDGs by using asteroseismology (Kepler, TESS) combined with spectroscopy (APOGEE, LAMOST), and astrometry (Gaia). We discovered 15 new CDGs in the Kepler field, and confirm that CDGs are rare, being only $0.15\%$ of our background sample. Remarkably, we find that our CDGs are almost exclusively in the red clump (RC) phase. Asteroseismic masses reveal that our CDGs are primarily low-mass stars ($M \lesssim$ 2~M$_{\odot}$), in contrast to previous studies which suggested they are intermediate mass ($M = 2.5 - 5.0~\rm M_{\odot}$) based on HR diagrams. A very high fraction of our CDGs ($50\%$) are also Li-rich giants. We observe a bimodal distribution of luminosity in our CDGs, with one group having normal RC luminosity and the other being a factor of two more luminous than expected for their masses. We find demarcations in chemical patterns and luminosities which lead us to split them into three groups: (i) normal-luminosity CDGs, (ii) over-luminous CDGs, and (iii) over-luminous highly-polluted CDGs. We conclude that a merger of a helium white dwarf with an RGB star is the most likely scenario for the two groups of over-luminous stars. Binary mass-transfer from intermediate-mass AGB stars is a possibility for the highly-polluted over-luminous group. For the normal-luminosity CDGs, we cannot distinguish between core He-flash pollution or lower-mass merger scenarios. Due to the overlap with the CDGs, Li-rich giants may have similar formation channels., Comment: 27 pages, 14 figures, accepted by ApJ

Based on the APOGEE survey we conducted a search for carbon-deficient red giants (CDGs). We found 103 new CDGs, increasing the number in the literature by more than a factor of 3. CDGs are very rare, representing $0.03$~per cent of giants. They appear as an extended tail off the normal carbon distribution. We show that they are found in all components of the Galaxy, contrary to previous findings. The location of CDGs in the Hertzsprung-Russell diagram (HRD) shows that they are primarily intermediate-mass stars ($2-4~\rm{M}_{\odot}$). Their extended distribution may indicate that CDGs can also sometimes have $M < 2.0~\rm{M}_{\odot}$. We attempted to identify the evolutionary phases of the CDGs using stellar model tracks. We found that the bulk of the CDGs are likely in the subgiant branch or red clump phase, whereas other CDGs may be in the red giant branch or early asymptotic giant branch phases. Degeneracy in the HRD makes exact identification difficult. We examined their C, N, and O compositions and confirmed previous studies showing that the envelope material has undergone extensive hydrogen burning through the CN(O) cycle. The new-CDGs have [C+N+O/Fe] that generally sum to zero, indicating that they started with scaled-solar composition. However, the previously known-CDGs generally have [C+N+O/Fe$] > 0.0$, indicating that some He-burning products were added to their envelopes. As to the site(s) in which this originally occurred, we do not find a convincing solution., Comment: 12 pages, 5 figures, accepted by MNRAS

Asteroseismology provides a new avenue for accurately measuring the masses of evolved globular cluster (GC) stars through the detection of their solar-like oscillations. We present the first detections of solar-like oscillations in 47 red giant branch (RGB) and early asymptotic giant branch (EAGB) stars in the metal-poor GC M80; only the second ever with measured seismic masses. We investigate two major areas of stellar evolution and GC science; the multiple populations and stellar mass-loss. We detected a distinct bimodality in the EAGB mass distribution. We showed that this is likely due to sub-population membership. If confirmed, it would be the first direct measurement of a mass difference between sub-populations. A mass difference was not detected between the sub-populations in our RGB sample. We instead measured an average RGB mass of $0.782\pm0.009~\msun$, which we interpret as the average between the sub-populations. Differing mass-loss rates on the RGB has been proposed as the second parameter that could explain the horizontal branch (HB) morphology variations between GCs. We calculated an integrated RGB mass-loss separately for each sub-population: $0.12\pm0.02~\msun$ (SP1) and $0.25\pm0.02~\msun$ (SP2). Thus, SP2 stars have greatly enhanced mass-loss on the RGB. Mass-loss is thought to scale with metallicity, which we confirm by comparing our results to a higher metallicity GC, M4. We also find that M80 stars have insignificant mass-loss on the HB. This is different to M4, suggesting that there is a metallicity and temperature dependence in the HB mass-loss. Finally, our study shows the robustness of the $\Delta\nu$-independent mass scaling relation in the low-metallicity (and low-surface gravity) regime., Comment: 20 pages, 11 figures

Mass loss remains a major uncertainty in stellar modelling. In low-mass stars, mass loss is most significant on the red giant branch (RGB), and will impact the star's evolutionary path and final stellar remnant. Directly measuring the mass difference of stars in various phases of evolution represents one of the best ways to quantify integrated mass loss. Globular clusters (GCs) are ideal objects for this. M4 is currently the only GC for which asteroseismic data exists for stars in multiple phases of evolution. Using K2 photometry, we report asteroseismic masses for 75 red giants in M4, the largest seismic sample in a GC to date. We find an integrated RGB mass loss of $\Delta\bar{M} = 0.17 \pm 0.01 ~\mathrm{M}_{\odot}$, equivalent to a Reimers' mass-loss coefficient of $\eta_R = 0.39$. Our results for initial mass, horizontal branch mass, $\eta_R$, and integrated RGB mass loss show remarkable agreement with previous studies, but with higher precision using asteroseismology. We also report the first detections of solar-like oscillations in early asymptotic giant branch (EAGB) stars in GCs. We find an average mass of $\bar{M}_{\text{EAGB}}=0.54 \pm 0.01 ~\mathrm{M}_{\odot}$, significantly lower than predicted by models. This suggests larger-than-expected mass loss on the horizontal branch. Alternatively, it could indicate unknown systematics in seismic scaling relations for the EAGB. We discover a tentative mass bi-modality in the RGB sample, possibly due to the multiple populations. In our red horizontal branch sample, we find a mass distribution consistent with a single value. We emphasise the importance of seismic studies of GCs since they could potentially resolve major uncertainties in stellar theory., Comment: 15 pages, 12 figures, accepted by MNRAS

We report novel observational evidence on the evolutionary status of lithium-rich giant stars by combining asteroseismic and lithium abundance data. Comparing observations and models of the asteroseismic gravity-mode period spacing $\Delta\Pi_{1}$, we find that super-Li-rich giants (SLR, A(Li)~$> 3.2$~dex) are almost exclusively young red-clump (RC) stars. Depending on the exact phase of evolution, which requires more data to refine, SLR stars are either (i) less than $\sim 2$~Myr or (ii) less than $\sim40$~Myr past the main core helium flash (CHeF). Our observations set a strong upper limit for the time of the inferred Li-enrichment phase of $< 40$~Myr post-CHeF, lending support to the idea that lithium is produced around the time of the CHeF. In contrast, the more evolved RC stars ($> 40$~Myr post-CHeF) generally have low lithium abundances (A(Li)~$<1.0$~dex). Between the young, super-Li-rich phase, and the mostly old, Li-poor RC phase, there is an average reduction of lithium by about 3 orders of magnitude. This Li-destruction may occur rapidly. We find the situation to be less clear with stars having Li abundances between the two extremes of super-Li-rich and Li-poor. This group, the `Li-rich' stars ($3.2 >$~A(Li)~$> 1.0$~dex), shows a wide range of evolutionary states., Comment: 9 pages, 5 figures, 1 table, accepted in Astrophysical Journal Letters

The vast majority of stars with mass similar to the Sun are expected to only destroy lithium over the course of their lives, via low-temperature nuclear burning. This has now been supported by observations of hundreds of thousands of red giant stars (Brown et al. 1989, Kumar et al. 2011, Deepak et al. 2019, Singh et al. 2019, Casey et al. 2019). Here we perform the first large-scale systematic investigation into the Li content of stars in the red clump phase of evolution, which directly follows the red giant branch phase. Surprisingly we find that all red clump stars have high levels of lithium for their evolutionary stage. On average the lithium content increases by a factor of 40 after the end of the red giant branch stage. This suggests that all low-mass stars undergo a lithium production phase between the tip of the red giant branch and the red clump. We demonstrate that our finding is not predicted by stellar theory, revealing a stark tension between observations and models. We also show that the heavily studied (Brown et al. 1989, Reddy et al. 2005, Kumar et al. 2011, Singh et al. 2019, Casey et al. 2019) very Li-rich giants, with A(Li) $> +1.5$ dex, represent only the extreme tail of the lithium enhancement distribution, comprising 3% of red clump stars. Our findings suggest a new definition limit for Li-richness in red clump stars, A(Li) $> -0.9$ dex, which is much lower than the limit of A(Li) $> +1.5$ dex used over many decades (Brown et al. 1989, Castilho et al. 1995, Reddy et al. 2005, Carlberg et al. 2016, Casey et al. 2019, Holanda et al. 2020)., Comment: 17 pages, 7 figures

Lithium depletion and enrichment in the cosmos is not yet well understood. To help tighten constraints on stellar and Galactic evolution models, we present the largest high-resolution analysis of Li abundances A(Li) to date, with results for over 100 000 GALAH field stars spanning effective temperatures $5900\,\mathrm{K} \lesssim \rm{T_{eff}} \lesssim7000\,\mathrm{K}$ and metallicities $-3 \lesssim \rm[Fe/H] \lesssim +0.5$. We separated these stars into two groups, on the warm and cool side of the so-called Li-dip, a localised region of the Kiel diagram wherein lithium is severely depleted. We discovered that stars in these two groups show similar trends in the A(Li)-[Fe/H] plane, but with a roughly constant offset in A(Li) of 0.4 dex, the warm group having higher Li abundances. At $\rm[Fe/H]\gtrsim-0.5$, a significant increasing in Li abundance with increasing metallicity is evident in both groups, signalling the onset of significant Galactic production. At lower metallicity, stars in the cool group sit on the Spite plateau, showing a reduced lithium of around 0.4 dex relative to the primordial value predicted from Big Bang nucleosynthesis (BBN). However, stars in the warm group between [Fe/H] = -1.0 and -0.5, form an elevated plateau that is largely consistent with the BBN prediction. This may indicate that these stars in fact preserve the primordial Li produced in the early Universe., Comment: 5 pages, 2 figures, Accepted for publication by MNRAS

The slow ($s$) and intermediate ($i$) neutron ($n$) capture processes occur both in asymptotic giant branch (AGB) stars, and in massive stars. To study the build-up of the $s$- and $i$-products at low metallicity, we investigate the abundances of Y, Ba, La, Nd, and Eu in 98 stars, at $-2.4<\text{[Fe/H]}<-0.9$, in the Sculptor dwarf spheroidal galaxy. The chemical enrichment from AGB stars becomes apparent at $\text{[Fe/H]}\approx-2$ in Sculptor, and causes [Y/Ba], [La/Ba], [Nd/Ba] and [Eu/Ba] to decrease with metallicity, reaching subsolar values at the highest $\text{[Fe/H]}\approx-1$. To investigate individual nucleosynthetic sites, we compared three $n$-rich Sculptor stars with theoretical yields. One carbon-enhanced metal-poor (CEMP-no) star with high $\text{[Sr, Y, Zr]}>+0.7$ is best fit with a model of a rapidly-rotating massive star, the second (likely CH star) with the $i$-process, while the third has no satisfactory fit. For a more general understanding of the build-up of the heavy elements, we calculate for the first time the cumulative contribution of the $s$- and $i$-processes to the chemical enrichment in Sculptor, and compare with theoretical predictions. By correcting for the $r$-process, we derive $\text{[Y/Ba]}_{s/i}=-0.85\pm0.16$, $\text{[La/Ba]}_{s/i}=-0.49\pm0.17$, and $\text{[Nd/Ba]}_{s/i}=-0.48\pm0.12$, in the overall $s$- and/or $i$-process in Sculptor. These abundance ratios are within the range of those of CEMP stars in the Milky Way, which have either $s$- or $i$-process signatures. The low $\text{[Y/Ba]}_{s/i}$ and $\text{[La/Ba]}_{s/i}$ that we measure in Sculptor are inconsistent with them arising from the $s$-process only, but are more compatible with models of the $i$-process. Thus we conclude that both the $s$- and $i$-processes were important for the build-up of $n$-capture elements in the Sculptor dwarf spheroidal galaxy., Comment: Accepted in A&A, 16 pages, 10 figures, Paper II in series

We simulate stellar convection at high Reynolds number (Re$\lesssim$7000) with causal time stepping but no explicit viscosity. We use the 3D Euler equations with shock capturing (Colella & Woodward 1984). Anomalous dissipation of turbulent kinetic energy occurs as an emergent feature of advection ("Onsager damping"), caused by the moderate shocks which terminate the turbulent kinetic energy spectrum; see also (Perry 2021). In strongly stratified stellar convection the asymptotic limit for the global damping length of turbulent kinetic energy is $\ell_d \sim \langle u^3 \rangle /\langle \epsilon \rangle$. This "dissipative anomaly" (Onsager 1949) fixes the value of the "mixing length parameter", $\alpha = \ell_{\rm MLT}/H_P =\overline{\langle\Gamma_1\rangle}$, which is $\sim\, 5/3$ for complete ionization. The estimate is numerically robust, agrees to within 10% with estimates from stellar evolution with constant $\alpha$. For weak stratification $\ell_d$ shrinks to the depth of a thin convective region. Our flows are filamentary, produce surfaces of separation at boundary layers, resolve the energy-containing eddies, and develop a turbulent cascade down to the grid scale which agrees with the $4096^3$ direct numerical simulation of Kaneda (2003). The cascade converges quickly, and satisfies a power-law velocity spectrum similar to Kolmogorov (1941). Our flows exhibit intermittency, anisotropy, and interactions between coherent structures, features missing from K41 theory. We derive a dissipation rate from Reynolds stresses which agrees with (i) our flows, (ii) experiment (Warhaft 2002), and (iii) high Re simulations of the Navier-Stokes equations (Iyer, et al. 2018)., Comment: 31 pages, 6 figures. The time evolution and the fly-through movies may be found at \url{http://www.astro.keele.ac.uk/shyne/321D/convection-and-convective-boundary-mixing/visualisations}

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

Context. Classical novae are thermonuclear explosions hosted by accreting white dwarfs in stellar binary systems. Material piles up on top of the white dwarf star under mildly degenerate conditions, driving a thermonuclear runaway. The energy released by the suite of nuclear processes operating at the envelope (mostly proton-capture reactions and beta-decays) heats the material up to peak temperatures ranging from 100 to 400 MK. In these events, about 10-3 - 10-7 Msun, enriched in CNO and, sometimes, other intermediate-mass elements (e.g., Ne, Na, Mg, Al), are ejected into the interstellar medium. Aims. To date, most of the efforts undertaken in the modeling of classical nova outbursts have focused on the early stages of the explosion and ejection, ignoring the interaction of the ejecta, first with the accretion disk orbiting the white dwarf, and ultimately with the secondary star. Methods. A suite of three-dimensional, SPH simulations of the interaction between the nova ejecta, the accretion disk, and the stellar companion have been performed to fill this gap, aimed at testing the influence of the different parameters (i.e., mass and velocity of the ejecta, mass and geometry of the accretion disk) on the dynamical and chemical properties of the system. Results. We discuss the conditions that lead to the disruption of the accretion disk and to mass loss from the binary system. In addition, we discuss the likelihood of chemical contamination of the stellar secondary induced by the impact with the nova ejecta and its potential effect on the next nova cycle., Comment: Accepted for Publication in Astronomy and Astrophysics

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

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

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

Large star-to-star variations of the abundances of proton-capture elements, such as Na and O, in globular clusters (GCs) are interpreted as the effect of internal pollution resulting from the presence of multiple stellar populations. To better constrain this scenario we investigate the abundance distribution of the heavy element rubidium (Rb) in NGC 6752, NGC 1904, and NGC 104 (47 Tuc). Combining the results from our sample with those in the literature, we found that Rb exhibits no star-to-star variations, regardless the cluster metallicity, with the possible intriguing, though very uncertain, exception of the metal-rich bulge cluster NGC 6388. If no star-to-star variations will be confirmed for all GCs, it implies that the stellar source of the proton-capture element variations must not have produced significant amounts of Rb. This element is observed to be enhanced at extremely high levels in intermediate-mass AGB (IM-AGB) stars in the Magellanic Clouds (i.e., at a metallicity similar to 47 Tuc and NGC 6388). This may present a challenge to this popular candidate polluter, unless the mass range of the observed IM-AGB stars does not participate in the formation of the second-generation stars in GCs. A number of possible solutions are available to resolve this conundrum, also given that the Magellanic Clouds observations are very uncertain and may need to be revised. The fast rotating massive stars scenario would not face this potential problem as the slow mechanical winds of these stars during their main-sequence phase do not carry any Rb enhancements; however, these candidates face even bigger issues such as the production of Li and the close over-imposition with core-collapse supernova timescales. Observations of Sr, Rb, and Zr in metal-rich clusters such as NGC 6388 and NGC 6441 are sorely needed to clarify the situation., Comment: Accepted for publication in The Astrophysical Journal

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 stellar initial mass function (IMF) plays a crucial role in determining the number of surviving stars in galaxies, the chemical composition of the interstellar medium, and the distribution of light in galaxies. A key unsolved question is whether the IMF is universal in time and space. Here we use state-of-the-art results of stellar evolution to show that the IMF of our Galaxy made a transition from an IMF dominated by massive stars to the present-day IMF at an early phase of the Galaxy formation. Updated results from stellar evolution in a wide range of metallicities have been implemented in a binary population synthesis code, and compared with the observations of carbon-enhanced metal-poor (CEMP) stars in our Galaxy. We find that applying the present-day IMF to Galactic halo stars causes serious contradictions with four observable quantities connected with the evolution of AGB stars. Furthermore, a comparison between our calculations and the observations of CEMP stars may help us to constrain the transition metallicity for the IMF which we tentatively set at [Fe/H] = -2. A novelty of the current study is the inclusion of mass loss suppression in intermediate-mass AGB stars at low-metallicity. This significantly reduces the overproduction of nitrogen-enhanced stars that was a major problem in using the high-mass star dominated IMF in previous studies. Our results also demonstrate that the use of the present day IMF for all time in chemical evolution models results in the overproduction of Type I.5 supernovae. More data on stellar abundances will help to understand how the IMF has changed and what caused such a transition., Comment: 8 pages, 2 figures, accepted by MNRAS Letter

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

Models of primordial and hyper-metal-poor stars with masses similar to the Sun experience an ingestion of protons into the hot core during the core helium flash phase at the end of their red giant branch evolution. This produces a concurrent secondary flash powered by hydrogen burning that gives rise to further nucleosynthesis in the core. We perform post-process nucleosynthesis calculations on a one-dimensional stellar evolution calculation of a star of 1 solar mass and metallicity [Fe/H] = -6.5 that suffers a proton ingestion episode. Our network includes 320 nuclear species and 2,366 reactions and treats mixing and burning simultaneously. The mixing and burning of protons into the hot convective core leads to the production of 13C, which then burns via the 13C(alpha,n)16O reaction releasing a large number of free neutrons. During the first two years of neutron production the neutron poison 14N abundance is low, allowing the prodigious production of heavy elements such as strontium, barium, and lead via slow neutron captures (the s process). These nucleosynthetic products are later mixed to the stellar surface and ejected via stellar winds. We compare our results with observations of the hyper-metal-poor halo star HE 1327-2326, which shows a strong Sr overabundance. Our model provides the possibility of self-consistently explaining the Sr overabundance in HE 1327-2326 together with its C, N, and O overabundances (all within a factor of ~4) if the material were heavily diluted, for example, via mass transfer in a wide binary system. The model produces at least 18 times too much Ba than observed, but this may be within the large modelling uncertainties. In this scenario, binary systems of low mass must have formed in the early Universe. If true then this puts constraints on the primordial initial mass function., Comment: Accepted for publication on Astronomy & Astrophysics Letters

Degenerate ignition of helium in low-mass stars at the end of the red giant branch phase leads to dynamic convection in their helium cores. One-dimensional (1D) stellar modeling of this intrinsically multi-dimensional dynamic event is likely to be inadequate. Previous hydrodynamic simulations imply that the single convection zone in the helium core of metal-rich Pop I stars grows during the flash on a dynamic timescale. This may lead to hydrogen injection into the core, and a double convection zone structure as known from one-dimensional core helium flash simulations of low-mass Pop III stars. We perform hydrodynamic simulations of the core helium flash in two and three dimensions to better constrain the nature of these events. To this end we study the hydrodynamics of convection within the helium cores of a 1.25 \Msun metal-rich Pop I star (Z=0.02), and a 0.85 \Msun metal-free Pop III star (Z=0) near the peak of the flash. These models possess single and double convection zones, respectively. We use 1D stellar models of the core helium flash computed with state-of-the-art stellar evolution codes as initial models for our multidimensional hydrodynamic study, and simulate the evolution of these models with the Riemann solver based hydrodynamics code Herakles which integrates the Euler equations coupled with source terms corresponding to gravity and nuclear burning. The hydrodynamic simulation of the Pop I model involving a single convection zone covers 27 hours of stellar evolution, while the first hydrodynamic simulations of a double convection zone, in the Pop III model, span 1.8 hours of stellar life. We find differences between the predictions of mixing length theory and our hydrodynamic simulations. The simulation of the single convection zone in the Pop I model shows a strong growth of the size of the convection zone due to turbulent entrainment. Hence we predict that for the Pop I model a hydrogen injection phase (i.e. hydrogen injection into the helium core) will commence after about 23 days, which should eventually lead to a double convection zone structure known from 1D stellar modeling of low-mass Pop III stars. Our two and three-dimensional hydrodynamic simulations of the double (Pop III) convection zone model show that the velocity field in the convection zones is different from that predicted by stellar evolutionary calculations. The simulations suggest that the double convection zone decays quickly, the flow eventually being dominated by internal gravity waves., Comment: 16 pages, 18 figures, submitted to AaA

(Abridged) We demonstrate that the amount of extra mixing required to fit the observed low C/N and 12C/13C ratios in first giant branch (FGB) stars is also sufficient to explain the C and N abundances of Galactic AGB stars. We simulate the effect of extra mixing on the FGB by setting the composition of the envelope to that observed in low-mass FGB stars, and then evolve the models to the tip of the AGB. The inclusion of FGB extra mixing compositional changes has a strong effect on the C and N abundance in our AGB models, leading to compositions consistent with those measured in Galactic C-rich stars. The composition of the models is also consistent with C abundances measured in mainstream silicon carbide grains. While our models cover the range of C abundances measured in C stars in NGC 1846, we cannot simultaneously match the composition of the O and C-rich stars. Our models only match the O isotopic composition of K and some M, MS giants, and are not able to match the O composition of C-rich AGB stars. By increasing the 16O intershell abundance (based on observational evidence) it is possible to reproduce the observed trend of increasing 16O/18O and 16O/17O ratios with evolutionary phase. We conclude 1) if extra mixing occurs during the AGB it likely only occurs efficiently in low metallicity objects, or when the stars are heavily obscured making spectroscopic observations difficult, and 2) that the intershell compositions of AGB stars needs further investigation., Comment: 24 pages, accepted for publication in ApJ

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