Search

Theoretical physical-chemical models for the formation of planetary systems depend on data quality for the Sun's composition, that of stars in the solar neighbourhood, and of the estimated "pristine" compositions for stellar systems. The effective scatter and the observational uncertainties of elements within a few hundred parsecs from the Sun, even for the most abundant metals like carbon, oxygen and silicon, are still controversial. Here we analyse the stellar production and the chemical evolution of key elements that underpin the formation of rocky (C, O, Mg, Si) and gas/ice giant planets (C, N, O, S). We calculate 198 galactic chemical evolution (GCE) models of the solar neighbourhood to analyse the impact of different sets of stellar yields, of the upper mass limit for massive stars contributing to GCE ($M_{\rm up}$) and of supernovae from massive-star progenitors which do not eject the bulk of the iron-peak elements (faint supernovae). Even considering the GCE variation produced via different sets of stellar yields, the observed dispersion of elements reported for stars in the Milky Way disk is not reproduced. Among others, the observed range of super-solar [Mg/Si] ratios, sub-solar [S/N], and the dispersion of up to 0.5 dex for [S/Si] challenge our models. The impact of varying $M_{\rm up}$ depends on the adopted supernova yields. Thus, observations do not provide a constraint on the M$_{\rm up}$ parametrization. When including the impact of faint supernova models in GCE calculations, elemental ratios vary by up to 0.1-0.2 dex in the Milky Way disk; this modification better reproduces observations., Comment: 36 pages, 26 figures, 1 Table, 1 Appendix, Accepted for publication in MNRAS

Type Ia supernova explosions (SNIa) are fundamental sources of elements for the chemical evolution of galaxies. They efficiently produce intermediate-mass (with Z between 11 and 20) and iron group elements - for example, about 70% of the solar iron is expected to be made by SNIa. In this work, we calculate complete abundance yields for 39 models of SNIa explosions, based on three progenitors - a 1.4M deflagration detonation model, a 1.0 double detonation model and a 0.8 M double detonation model - and 13 metallicities, with 22Ne mass fractions of 0, 1x10-7, 1x10-6, 1x10-5, 1x10-4, 1x10-3, 2x10-3, 5x10-3, 1x10-2, 1.4x10-2, 5x10-2, and 0.1 respectively. Nucleosynthesis calculations are done using the NuGrid suite of codes, using a consistent nuclear reaction network between the models. Complete tables with yields and production factors are provided online at Zenodo: Yields. We discuss the main properties of our yields in the light of the present understanding of SNIa nucleosynthesis, depending on different progenitor mass and composition. Finally, we compare our results with a number of relevant models from the literature., Comment: 42 pages, 21 figures. Accepted for publication in ApJS 21-06-23

Certain types of silicon carbide (SiC) grains, e.g., SiC-X grains, and low density (LD) graphites are C-rich presolar grains that are thought to have condensed in the ejecta of core-collapse supernovae (CCSNe). In this work we compare C, N, Al, Si, and Ti isotopic abundances measured in presolar grains with the predictions of 21 CCSN models. The impact of a range of SN explosion energies is considered, with the high energy models favouring the formation of a C/Si zone enriched in $^{12}$C, $^{28}$Si, and $^{44}$Ti. Eighteen of the 21 models have H ingested into the He-shell and different abundances of H remaining from such H-ingestion. CCSN models with intermediate to low energy (that do not develop a C/Si zone) cannot reproduce the $^{28}$Si and $^{44}$Ti isotopic abundances in grains without assuming mixing with O-rich CCSN ejecta. The most $^{28}$Si-rich grains are reproduced by energetic models when material from the C/Si zone is mixed with surrounding C-rich material, and the observed trends of the $^{44}$Ti/$^{48}$Ti and $^{49}$Ti/$^{48}$Ti ratios are consistent with the C-rich C/Si zone. For the models with H-ingestion, high and intermediate explosion energies allow the production of enough $^{26}$Al to reproduce the $^{26}$Al/$^{27}$Al measurements of most SiC-X and LD graphites. In both cases, the highest $^{26}$Al/$^{27}$Al ratio is obtained with H still present at $X_H \approx 0.0024$ in He-shell material when the SN shock is passing. The existence of H in the former convective He-shell points to late H-ingestion events in the last days before massive stars explode as a supernova.

Short-lived radioactive isotopes (SLRs) with half-lives between 0.1 to 100 Myr can be used to probe the origin of the Solar System. In this work, we examine the core-collapse supernovae production of the 15 SLRs produced: $^{26}$Al, $^{36}$Cl, $^{41}$Ca, $^{53}$Mn, $^{60}$Fe, $^{92}$Nb, $^{97}$Tc, $^{98}$Tc, $^{107}$Pd, $^{126}$Sn, $^{129}$I, $^{135}$Cs, $^{146}$Sm, $^{182}$Hf, and $^{205}$Pb. We probe the impact of the uncertainties of the core-collapse explosion mechanism by examining a collection of 62 core-collapse models with initial masses of 15, 20, and 25M$_{\odot}$, explosion energies between 3.4$\times$10$^{50}$ and 1.8$\times$10$^{52}$ ergs and compact remnant masses between 1.5M$_{\odot}$and 4.89M$_{\odot}$. We identify the impact of both explosion energy and remnant mass on the final yields of the SLRs. Isotopes produced within the innermost regions of the star, such as $^{92}$Nb and $^{97}$Tc, are the most affected by the remnant mass, $^{92}$Nb varying by five orders of magnitude. Isotopes synthesised primarily in explosive C-burning and explosive He-burning, such as $^{60}$Fe, are most affected by explosion energies. $^{60}$Fe increases by two orders of magnitude from the lowest to the highest explosion energy in the 15M$_{\odot}$model. The final yield of each examined SLR is used to compare to literature models.

We examine the composition of barium stars in the context of mass transfer from an asymptotic giant branch (AGB) companion. We accrete between 0.01 and 0.5 M$_\odot$ of AGB ejecta on to low mass companions of [Fe/H] = -0.25 at the ages expected for the end of the lives of AGB stars of 2.5, 3 and 4M$_\odot$. In each case, we form a star of 2.5 M$_\odot$ which is thought to be a typical barium star mass. We discuss the extent of dilution of accreted material as the star evolves, and describe the impact on the surface abundances. For accretion from a 2.5\ms\ primary, if the secondary's initial mass is 2.45 M$_\odot$ or more, accretion takes place when the secondary is undergoing core helium burning. Using data from the sample of De Castro et al., we attempt to fit the observed properties of 74 barium giants using the models we have computed. We find that all but six of these objects are best fit using ejecta from 2.5 M$_\odot$ (32 objects) or 3 M$_\odot$ (36 objects) AGB stars. Higher accretion masses are typically required when accreting from a lower mass companion. We find accretion masses that are broadly consistent with recent hydrodynamical simulations of wind mass transfer, though the accretion efficiency is toward the upper limit found in these simulations. For the 18 stars with reported orbital periods, we find no strong correlations between period and accretion mass., Comment: 13 pages, 7 figures. Accepted for publication in MNRAS 14/6/21

$\textit{Thailand-UK Python+Astronomy Summer School}$ (ThaiPASS), a collaborative project comprising UK and Thai institutions and assess its impact and possible application to schools in the United Kingdom. Since its inception in 2018, the annual ThaiPASS have trained around 60 Thai high-school students in basic data handling skills using Python in the context of various astronomy topics, using current research from the teaching team. Our impact assessment of the 5 day summer schools show an overwhelmingly positive response from students in both years, with over $80\%$ of students scoring the activities above average in all activities but one. We use this data to suggest possible future improvements. We also discuss how ThaiPASS may inspire further outreach and engagement activities within the UK and beyond., Comment: Accepted for publication into Physics Education; 21 pages, 5 figures and 3 tables

The nature of the progenitors of type Ia supernovae (SNe Ia) remains a mystery. Binary systems consisting of a white dwarf (WD) and a main-sequence (MS) donor are potential progenitors of SNe Ia, in which a thermonuclear explosion of the WD may occur when its mass reaches the Chandrasekhar limit during accretion of material from a companion star. In the present work, we address theoretical rates and delay times of a specific MS donor channel to SNe Ia, in which a He star+MS binary produced from a common envelope event subsequently forms a WD+MS system without the He star undergoing mass transfer by Roche lobe overflow. By combining the results of self-consistent binary evolution calculations with population synthesis models, we find that the contribution of SNe Ia in this channel is around 2.0e-4 per year. In addition, we find that delay times of SNe Ia in this channel cover a range of about 1.0-2.6 Gyr, and almost all SNe Ia produced in this way (about 97%) have a delay time of > 1 Gyr. While the rate of SN Ia in this work is about 10% of the overall SN Ia rate, the channel represents a possible contribution to the old population (1-3 Gyr) of observed SNe Ia., Comment: 6 pages, 4 figures, published in Astronomy & Astrophysics (A&A)

Lead (Pb) is predominantly produced by the slow neutron-capture process (s process) in asymptotic giant branch (AGB) stars. In contrast to significantly enhanced Pb abundances predicted by low-mass, low-metallicity AGB-models, observations of Magellanic post-AGB stars show incompatibly low Pb abundances. Observations of carbon-enhanced metal-poor (CEMP) stars whose s-process enrichments are accompanied by heavy elements traditionally associated with the rapid neutron-capture process (r process) have raised the need for a neutron-capture process operating at neutron densities intermediate to the s and r process: the so-called i process. We study i-process nucleosynthesis with single-zone nuclear-network calculations. Our i-process models can explain the heavy-element abundance patterns measured in Magellanic post-AGB stars including their puzzlingly low Pb abundances. Furthermore, the heavy-element enhancements in the post-AGB and CEMP-i stars, particularly their Pb abundance, allow us to characterise the neutron densities and exposures of the i process that produced the observed abundance patterns. We find that the lower-metallicity CEMP-i stars ($\left[ \mathrm{Fe} / \mathrm{H} \right] \approx -2.5$) have heavy-element abundances best matched by models with higher neutron densities and exposures ($\tau > 2.0 \, \mathrm{mbarn}^{-1}$) compared to the higher-metallicity post-AGB stars ($\left[ \mathrm{Fe} / \mathrm{H} \right] \approx -1.3$, $\tau < 1.3 \, \mathrm{mbarn}^{-1}$). This offers new constraints and insights regarding the properties of i-process sites and demonstrates that the responsible process operates on time scales of the order of a few years or less., Comment: accepted for publication in ApJ

The chemical enrichments detected in CEMP-$s$ stars are believed to be the consequence of a mass-transfer episode from a now extinct AGB primary star. This hypothesis is borne out by the fact that most CEMP-$s$ stars exhibit RV variations suggesting that they are in binary systems. We use the period distribution of an observed sample of CEMP-$s$ stars to investigate the constraints it imposes on our models of binary evolution and on the properties of the metal-poor binary population in the Galactic halo. We simulate binary populations using different assumptions about the initial period distribution and the mass-transfer process, and we compare the predicted period distributions of our synthetic CEMP-$s$ stars with the observed one. With a set of default assumptions often made in binary population-synthesis studies, the percentage of observed CEMP-$s$ systems with periods shorter than 2000 days is underestimated by almost a factor of 3, and by a factor of 2 between 3000 and 10000 days, while 40% of the simulated systems have periods longer than 10000 days, which is approximately the longest measured period among CEMP-$s$ stars. To reconcile the results of the models with the orbital properties of observed CEMP-$s$ stars, one or both of the following conditions are necessary: ($i$) the specific angular momentum carried away by the material that escapes the binary system is 2-5 times higher than currently predicted by analytical models and hydrodynamical simulations of wind mass transfer, and ($ii$) the initial period distribution of very metal-poor binary stars is significantly different from that observed in the solar vicinity and weighted towards periods shorter than about 10000 days. Our simulations show that some, perhaps all, of the observed CEMP-$s$ stars with apparently constant RV could be undetected binaries with periods longer than 10000 days, but (...), Comment: 18 pages, 12 figures, recommended for publication in A&A (section 7)

The nature of the progenitors of Type Ia supernovae (SNe Ia) remains a mystery. Comparing theoretical rates and delay-time distributions of SNe Ia with those inferred observationally can constrain their progenitor models. In this work, taking thermohaline mixing into account in the helium-enriched main-sequence (HEMS) donor scenario, we address rates and delay times of SNe Ia in this channel by combining the results of self-consistent binary evolution calculations with population synthesis models. We find that the Galactic SN Ia rate from the HEMS donor scenario is around 0.6-1.2 per century, which is about 30% of the observed rate. Delay times of SNe Ia in this scenario cover a wide range of 0.1-1.0 Gyr. We also present the pre-explosion properties of companion stars in the HEMS donor scenario, which will be helpful for placing constraints on SN Ia progenitors through analyzing their pre-explosion images., Comment: 11 pages, 8 figures, Accepted for publication in MNRAS

A large fraction of stars in binary systems are expected to undergo mass and angular momentum exchange at some point in their evolution, which can drastically alter the chemical and dynamical properties and fates of the systems. Interaction by stellar wind is an important process in wide binaries. However, the details of wind mass transfer are still not well understood. We perform three-dimensional hydrodynamical simulations of wind mass transfer in binary systems to explore mass accretion efficiencies and geometries of mass outflows, for a range of mass ratios from 0.05 to 1.0. In particular, we focus on the case of a free wind, in which some physical mechanism accelerates the expelled wind material balancing the gravity of the mass-losing star with the wind velocity comparable to the orbital velocity of the system. We find that the mass accretion efficiency and accreted specific angular momentum increase with the mass ratio of the system. For an adiabatic wind, we obtain that the accretion efficiency onto the secondary star varies from about 0.1% to 8% for mass ratios between 0.05 and 1.0., Comment: 15 pages, 9 figures and 1 table, Published in the ApJ. Published version is available here: http://iopscience.iop.org/article/10.3847/1538-4357/aa8622/meta

Hydrodynamical simulations predict that a large amount of hydrogen (>0.1 solar masses) is removed from a hydrogen-rich companion star by the SN explosion in the single-degenerate scenario of Type Ia supernovae (SNe Ia). However, non-detection of hydrogen-rich material in the late-time spectra of SNe Ia suggests that the hydrogen mass stripped from the progenitor system is <0.001-0.058 solar masses. In this letter we include thermohaline mixing into self-consistent binary evolution calculations for the helium-enriched main-sequence (HEMS) donor channel of SNe Ia for the first time. We find that the swept-up hydrogen masses expected in this channel are around 0.10-0.17 solar masses, which is higher than the observational limits, although the companion star is strongly helium-enriched when the SN explodes. This presents a serious challenge to the HEMS donor channel., Comment: 5 pages,2 figures, 1 table, Monthly Notices of the Royal Astronomical Society: Letters, Volume 470, Issue 1, p.L72-L76

Carbon-enhanced metal-poor (CEMP) stars in the Galactic Halo display enrichments in heavy elements associated with either the s (slow) or the r (rapid) neutron-capture process (e.g., barium and europium respectively), and in some cases they display evidence of both. The abundance patterns of these CEMP-s/r stars, which show both Ba and Eu enrichment, are particularly puzzling since the s and the r processes require neutron densities that are more than ten orders of magnitude apart, and hence are thought to occur in very different stellar sites with very different physical conditions. We investigate whether the abundance patterns of CEMP-s/r stars can arise from the nucleosynthesis of the intermediate neutron-capture process (the i process), which is characterised by neutron densities between those of the s and the r processes. Using nuclear network calculations, we study neutron capture nucleosynthesis at different constant neutron densities n ranging from $10^7$ to $10^{15}$ cm$^{-3}$. With respect to the classical s process resulting from neutron densities on the lowest side of this range, neutron densities on the highest side result in abundance patterns that show an increased production of heavy s-process and r-process elements but similar abundances of the light s-process elements. Such high values of n may occur in the thermal pulses of asymptotic giant branch (AGB) stars due to proton ingestion episodes. Comparison to the surface abundances of 20 CEMP-s/r stars show that our modelled i-process abundances successfully reproduce observed abundance patterns that could not be previously explained by s-process nucleosynthesis. Because the i-process models fit the abundances of CEMP-s/r stars so well, we propose that this class should be renamed as CEMP-i., Comment: 16 pages, 28 figures. Accepted for publication in ApJ

A large number of extremely low-mass helium white dwarfs (ELM WDs) have been discovered in recent years. The majority of them are found in close binary systems suggesting they are formed either through a common-envelope phase or via stable mass transfer in a low-mass X-ray binary (LMXB) or a cataclysmic variable (CV) system. Here, we investigate the formation of these objects through the LMXB channel with emphasis on the proto-WD evolution in environments with different metallicities. We study, for the first time, the combined effects of rotational mixing and element diffusion (e.g. gravitational settling, thermal and chemical diffusion) on the evolution of proto-WDs and on the cooling properties of the resulting WDs. We present state-of-the-art binary stellar evolution models computed with MESA for metallicities between Z=0.0002 and Z=0.02, producing WDs with masses between 0.16-0.45 M$_{\odot}$. Our results confirm that element diffusion plays a significant role in the evolution of proto-WDs that experience hydrogen shell flashes. The occurrence of these flashes produces a clear dichotomy in the cooling timescales of ELM WDs, which has important consequences e.g. for the age determination of binary millisecond pulsars. Rotational mixing is found to counteract the effect of gravitational settling in the surface layers of young, bloated ELM proto-WDs and therefore plays a key role in determining their surface chemical abundances. We predict that these proto-WDs have helium-rich envelopes through a significant part of their lifetime, a crucial ingredient for understanding the newly observed ELM proto-WD pulsators. The hydrogen envelope at detachment, although small compared to the total mass of the WD, contains enough angular momentum such that the spin frequency of the resulting WD on the cooling track is well above the orbital frequency., Comment: 28 pages, 10 pages tables, 19 figures, accepted for publication in A&A

The nature of the progenitors of SNe Ia is not yet fully understood. In the single-degenerate (SD) scenario, the collision of the SN ejecta with its companion star is expected to produce detectable UV emission in the first few days after the SN explosion within certain viewing angles. It was recently found that the $B-V$ colour of the nearby SN Ia SN 2012cg at about sixteen days before the maximum $B$-band brightness was about 0.2 mag bluer than those of other normal SNe Ia, which was reported as the first evidence for excess blue light from the interaction of normal SN Ia ejecta with its companion star. In this work, we compare current observations for SN 2012cg from its pre-explosion phase to the late-time nebular phase with theoretical predictions from binary evolution and population synthesis calculations for a variety of popular progenitor scenarios. We find that a main-sequence donor or a carbon-oxygen white dwarf donor binary system is more likely to be the progenitor of SN 2012cg. However, both scenarios also predict properties which are in contradiction to the observed features of this system. Taking both theoretical and observational uncertainties into account, we suggest that it might be too early to conclude that SN 2012cg was produced from an explosion of a Chandrasekhar-mass white dwarf in the SD scenario. Future observations and improved detailed theoretical modelling are still required to place a more stringent constraint on the progenitor of SN 2012cg., Comment: 10 pages, 6 figures, 1 table, accepted for publication in Monthly Notices of the Royal Astronomical Society (MNRAS)

CEMP-$r/s$ stars are metal-poor stars with enhanced abundances of carbon and heavy elements associated with the slow ($s$-) and rapid ($r$-) neutron-capture process. It is believed that carbon and $s$-elements were accreted from the wind of an AGB primary star, a scenario that is generally accepted to explain the formation of CEMP stars that are only enhanced in $s$-elements (CEMP-$s$ stars). The origin of $r$-element-enrichment in CEMP-$r/s$ stars is debated and many formation scenarios have been put forward. We aim to determine the likelihood of the scenarios proposed to explain the formation of CEMP-$r/s$ stars. We calculate the frequency of CEMP-$r/s$ stars among CEMP-$s$ stars for a variety of scenarios, and we compare it with that determined from an observed sample of CEMP-$r/s$ stars collected from the literature. The theoretical frequency of CEMP-$r/s$ stars predicted in most scenarios underestimates the observed ratio by at least a factor of 5. If the enrichments in $s$- and $r$-elements are independent, the model ratio of CEMP-$r/s$ to CEMP-$s$ stars is about 22%, that is approximately consistent with the lowest estimate of the observed ratio. However, this model predicts that about one third of all carbon-normal stars have [Ba/Fe] and [Eu/Fe] higher than 1, and that 40% of all CEMP stars have [Ba/Eu]$\le0$. Stars with these properties are at least ten times rarer in our observed sample. The $intermediate$ or $i$-process, which is supposedly active in some circumstances during the AGB phase, could provide an explanation of the origin of CEMP-$r/s$ stars, similar to that of CEMP-$s$ stars, in the context of wind mass accretion in binary systems. Further calculations of the nucleosynthesis of the $i$-process and of the detailed evolution of late AGB stars are needed to investigate if this scenario predicts a CEMP-$r/s$ star frequency consistent with the observations., Comment: 9 pages, 2 figures, accepted for publication on Astronomy & Astrophysics

The progenitors of many CCSNe are expected to be in binary systems. After the SN explosion, the companion may suffer from mass stripping and be shock heated as a result of the impact of the SN ejecta. If the binary system is disrupted, the companion is ejected as a runaway and hypervelocity star. By performing a series of 3D hydrodynamical simulations of the collision of SN ejecta with the companion star, we investigate how CCSN explosions affect their companions. We use the BEC code to construct the detailed companion structure at the time of SN explosion. The impact of the SN blast wave on the companion is followed by means of 3D SPH simulations using the Stellar GADGET code. For main-sequence (MS) companions, we find that the amount of removed mass, impact velocity, and chemical contamination of the companion that results from the impact of the SN ejecta, strongly increases with decreasing binary separation and increasing explosion energy. Their relationship can be approximately fitted by power laws, which is consistent with the results obtained from impact simulations of SNe~Ia. However, we find that the impact velocity is sensitive to the momentum profile of the outer SN ejecta and, in fact, may decrease with increasing ejecta mass, depending on the modeling of the ejecta. Because most companions to Ib/c CCSNe are in their MS phase at the moment of the explosion, combined with the strongly decaying impact effects with increasing binary separation, we argue that the majority of these SNe lead to inefficient mass stripping and shock heating of the companion star following the impact of the ejecta. Our simulations show that the impact effects of Ib/c SN ejecta on the structure of MS companions, and thus their long-term post-explosion evolution, is in general not dramatic. We find that at most 10% of their mass is lost, and their resulting impact velocities are less than 100 km/s., Comment: Accepted for publication in Astronomy and Astrophysics, some minor typographical errors are fixed, the affiliation of second author is corrected

The progenitors of SNe Ia are not yet fully understood. The two leading progenitor scenarios are the single-degenerate (SD) scenario and the double-degenerate scenario. In the SD scenario, the collision of the SN Ia ejecta with its companion star is expected to produce detectable ultraviolet (UV) emission in the first few days after the SN explosion within certain viewing angles. A strong UV flash has recently been detected in an SN 2002es-like peculiar SN Ia iPTF14atg by Cao et al., which is interpreted as evidence of an early-time UV signature due to SN ejecta interacting with its companion star, supporting the SD scenario. In this paper, we present the expected luminosity distributions of early-time UV emission arising from SN Ia ejecta-companion interaction by performing binary population synthesis calculations for different progenitor systems in the SD scenario. Our theoretical predictions will be helpful for future early-time observations of SNe Ia to constrain their possible progenitors. Assuming the observed strong UV pulse of iPTF14atg was indeed produced by the SN ejecta-companion interaction, our population synthesis model suggests that the progenitor system of iPTF14atg is most likely a red-giant donor binary system, and it is unlikely to have been a main-sequence or helium-star donor system., Comment: 11 pages, 7 figures, accepted for publication in Monthly Notices of the Royal Astronomical Society (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

Type Iax supernovae (SNe Iax) are proposed as one new sub-class of SNe Ia since they present observational properties that are sufficiently distinct from the bulk of SNe Ia. SNe Iax are the most common of all types of peculiar SNe by both number and rate, with an estimated rate of occurrence of about 5-30% of the total SN Ia rate. However, the progenitor systems of SNe Iax are still uncertain. Analyzing pre-explosion images at SN Iax positions provides a direct way to place strong constraints on the nature of progenitor systems of SNe Iax. In this work, we predict pre-explosion properties of binary companion stars in a variety of potential progenitor systems by performing detailed binary evolution calculations with the one-dimensional stellar evolution code STARS. This will be helpful for constraining progenitor systems of SNe Iax from their pre-explosion observations. With our binary evolution calculations, it is found that the non-degenerate helium (He) companion star to both a massive C/O WD (> 1.1 solar mass) and a hybrid C/O/Ne WD can provide an explanation for the observations of SN~2012Z-S1, but the hybrid WD+He star scenario is more favorable., Comment: Journal reference: The Astrophysical Journal, 2015, 808, 138

As part of a larger program aimed at better quantifying the uncertainties in stellar computations, we attempt to calibrate the extent of convective overshooting in low to intermediate mass stars by means of eclipsing binary systems. We model 12 such systems, with component masses between 1.3 and 6.2 solar masses, using the detailed binary stellar evolution code STARS, producing grids of models in both metallicity and overshooting parameter. From these, we determine the best fit parameters for each of our systems. For three systems, none of our models produce a satisfactory fit. For the remaining systems, no single value for the convective overshooting parameter fits all the systems, but most of our systems can be well described with an overshooting parameter between 0.09 and 0.15, corresponding to an extension of the mixed region above the core of about 0.1-0.3 pressure scale heights. Of the nine systems where we are able to obtain a good fit, seven can be reasonably well fit with a single parameter of 0.15. We find no evidence for a trend of the extent of overshooting with either mass or metallicity, though the data set is of limited size. We repeat our calculations with a second evolution code, MESA, and we find general agreement between the two codes. For the extension of the mixed region above the convective core required by the MESA models is about 0.15-0.4 pressure scale heights. For the system EI Cep, we find that MESA gives an overshooting region that is larger than the STARS one by about 0.1 pressure scale heights for the primary, while for the secondary the difference is only 0.05 pressure scale heights., Comment: 9 pages, 15 figures. Accepted for publication in A&A (20/01/2015)

SNe Iax are proposed as one new sub-class of SNe Ia. SNe Iax have been estimated to account for ~5-30% of the total SN Ia rate, and most SNe Iax have been discovered in late-type galaxies. In addition, observations constrain the progenitor systems of some SN Iax progenitors have ages of <80Myr. Recently, the weak deflagration explosions of Ch-mass CO WDs seem to provide a viable physical scenario for SNe Iax. Here, comparing theoretical predictions from BPS calculations with SN Iax observations, we put constraints on the SD Ch-mass model as a possible Iax progenitor. A long delay times of >3Gyr and low SN rates of ~3e-5/yr are found in our RG donor channel, indicating that this channel is unlikely to produce SNe Iax. We predict that the Galactic rate from the MS (He) donor channel is ~1.5e-3/yr (3e-4/yr), which is consistent with the observed Iax rate. The short delay times in the He channel support the young host environments of SNe Iax. However, the relative long delay times in the MS donor channel are less favourable for the observational constraints on the Iax progenitor ages. Finally, we set an upper limit on the pre-SN mass-loss rate at < 10e-4 M_sun/yr (for V_wind=100 km/s). The delay times in the SD Ch-mass model do not account for that most SNe Iax are located in late-type galaxies. However, at least one Iax event (SN 2008ge) is hosted by a S0 galaxy with no signs of star formation. Moreover, current X-ray observations cannot rule out the SD Ch-mass model. Taking all these into account and considering the uncertainty of the observed rate for SNe Iax, we suggest that some SNe Iax may be produced from weak deflagrations of Ch-mass CO WDs in SD progenitors, especially in the He-star channel. This is consistent with recent analysis of HST observations for the SN Iax SN2012Z. However, this SD deflagration model is still unlikely to be the most common progenitor scenario for SNe Iax., Comment: Some minor corrections are made

We revisit the observed frequencies of Carbon-Enhanced Metal-Poor (CEMP) stars as a function of the metallicity in the Galaxy, using data from the literature with available high-resolution spectroscopy. Our analysis excludes stars exhibiting clear over-abundances of neutron-capture elements, and takes into account the expected depletion of surface carbon abundance that occurs due to CN processing on the upper red-giant branch. This allows for the recovery of the initial carbon abundance of these stars, and thus for an accurate assessment of the frequencies of carbon-enhanced stars. The correction procedure we develope is based on stellar-evolution models, and depends on the surface gravity, log g, of a given star. Our analysis indicates that, for stars with [Fe/H]<=-2.0, 20% exhibit [C/Fe]>=+0.7. This fraction increases to 43% for [Fe/H]<=-3.0 and 81% for [Fe/H]<=-4.0, which is higher than have been previously inferred without taking the carbon-abundance correction into account. These CEMP-star frequencies provide important inputs for Galactic and stellar chemical-evolution models, as they constrain the evolution of carbon at early times and the possible formation channels for the CEMP-no stars. We also have developed a public online tool with which carbon corrections using our procedure can be easily obtained., Comment: 27 pages, 16 figures; accepted for publication in ApJ

We report estimated carbon-abundance ratios, [C/Fe], for seven newly-discovered carbon-enhanced metal-poor (CEMP) RR Lyrae stars. These are well-studied RRab stars that had previously been selected as CEMP candidates based on low-resolution spectra. For this pilot study, we observed eight of these CEMP RR Lyrae candidates with the Wide Field Spectrograph (WiFeS) on the ANU 2.3m telescope. Prior to this study, only two CEMP RR Lyrae stars had been discovered: TY Gru and SDSS J1707+58. We compare our abundances to new theoretical models of the evolution of low-mass stars in binary systems. These simulations evolve the secondary stars, post accretion from an AGB donor, all the way to the RR Lyrae stage. The abundances of CEMP RR Lyrae stars can be used as direct probes of the nature of the donor star, such as its mass, and the amount of material accreted onto the secondary. We find that the majority of the sample of CEMP RR Lyrae stars is consistent with AGB donor masses of around 1.5 - 2.0 solar masses and accretion masses of a few hundredths of a solar mass. Future high-resolution studies of these newly-discovered CEMP RR Lyrae stars will help disentangle the effects of the proposed mixing processes that occur in such objects., Comment: Accepted for publication in ApJ

Certain carbon-enhanced metal-poor stars likely obtained their composition via pollution from some of the earliest generations of asymptotic giant branch stars and as such provide important clues to early Universe nucleosynthesis. Recently, Kinman et al. discovered that the highly carbon- and barium-enriched metal-poor star SDSS J1707+58 is in fact an RR Lyrae pulsator. This gives us an object in a definite evolutionary state where the effects of dilution of material during the Main Sequence are minimised owing to the object having passed through first dredge-up. We perform detailed stellar modelling of putative progenitor systems in which we accreted material from asymptotic giant branch stars in the mass range 1-2 solar masses. We investigate how the surface abundances are affected by the inclusion of mechanisms like thermohaline mixing and gravitational settling. While we are able to find a reasonable fit to the carbon and sodium abundances of SDSS J1707+58, suggesting accretion of around 0.1 solar masses from a 2 solar mass companion, the strontium and barium abundances remain problematic and this object may have experienced something other than a standard $s$ process. We have more success in fitting the abundances of the mildly carbon-enriched, metal-poor RR Lyrae pulsator TY Gru (CS 22881-071), which we suggest received 0.1 solar masses of material from a companion of around 1 solar mass., Comment: 10 pages, 5 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

We present a comprehensive abundance analysis of two newly-discovered carbon-enhanced metal-poor (CEMP) stars. HE2138-3336 is a s-process-rich star with [Fe/H] = -2.79, and has the highest [Pb/Fe] abundance ratio measured thus far, if NLTE corrections are included ([Pb/Fe] = +3.84). HE2258-6358, with [Fe/H] = -2.67, exhibits enrichments in both s- and r-process elements. These stars were selected from a sample of candidate metal-poor stars from the Hamburg/ESO objective-prism survey, and followed up with medium-resolution (R ~ 2,000) spectroscopy with GEMINI/GMOS. We report here on derived abundances (or limits) for a total of 34 elements in each star, based on high-resolution (R ~ 30,000) spectroscopy obtained with Magellan-Clay/MIKE. Our results are compared to predictions from new theoretical AGB nucleosynthesis models of 1.3 Mo with [Fe/H] = -2.5 and -2.8, as well as to a set of AGB models of 1.0 to 6.0 Mo at [Fe/H] = -2.3. The agreement with the model predictions suggests that the neutron-capture material in HE2138-3336 originated from mass transfer from a binary companion star that previously went through the AGB phase, whereas for HE2258-6358, an additional process has to be taken into account to explain its abundance pattern. We find that a narrow range of progenitor masses (1.0 < M(Mo) < 1.3) and metallicities (-2.8 < [Fe/H] < -2.5) yield the best agreement with our observed elemental abundance patterns., Comment: Accepted for publication in ApJ

Most carbon-enhanced metal-poor (CEMP) stars are thought to result from past mass transfer of He-burning material from an asymptotic giant branch (AGB) star to a low-mass companion star, which we now observe as a CEMP star. Because AGB stars of intermediate mass efficiently cycle carbon into nitrogen in their envelopes, the same evolution scenario predicts the existence of a population of nitrogen-enhanced metal-poor (NEMP) stars, with [N/Fe] > 1 and [N/C] > 0.5. Such NEMP stars are rare, although their occurrence depends on metallicity: they appear to be more common at [Fe/H] < -2.8 by about a factor of 10 compared to less metal-poor stars. We analyse the observed sample of metal-poor stars with measurements of both carbon and nitrogen to derive firm constraints on the occurrence of NEMP stars as a function of metallicity. We compare these constraints to binary population synthesis calculations in which we vary the initial distributions of mass, mass ratio and binary orbital periods. We show that the observed paucity of NEMP stars at [Fe/H] > -2.8 does not allow for large modifications in the initial mass function, as have been suggested in the literature to account for the high frequency of CEMP stars. The situation at lower metallicity is less clear, and we do not currently have stellar models to perform this comparison for [Fe/H] < -2.8. However, unless intermediate-mass AGB stars behave very differently at such low metallicity, the observed NEMP frequency at [Fe/H] < -2.8 appears incompatible with the top-heavy forms of the initial mass function suggested in the literature., Comment: Accepted for publication in Astronomy & Astrophysics. 11 pages, 5 figures, 3 tables and 1 online appendix

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 present models for the slow neutron-capture process (s process) in asymptotic giant branch (AGB) stars of metallicity [Fe/H]=-2.3 and masses 0.9 Msun to 6 Msun. We encountered different regimes of neutron-capture nucleosynthesis increasing in importance as the stellar mass decreases: the 22Ne(alpha,n)25Mg reaction activated during the thermal pulses, the 13C(alpha,n)16O reaction activated in radiative conditions during the interpulse periods, and the 13C(alpha,n)16O reaction activated during the thermal pulses, also as a result of mild proton ingestion episodes. The models where the 13C burns radiatively (masses ~ 2 Msun) produce an overall good match to carbon-enhanced metal-poor (CEMP) stars showing s-process enhancements (CEMP-s), except they produce too much Na and F. On the other hand, none of our models can provide a match to the composition of CEMP stars also showing rapid-process enhancements (CEMP-s/r). The models fail to reproduce the observed Eu abundances, and they also fail to reproduce the correlation between the Eu and Ba abundances. They also cannot match the ratio of heavy to light s-process elements observed in many CEMP-s/r stars, which can be more than ten times higher than in the solar system. To explain the composition of CEMP-s/r stars we need to invoke the existence of a different "s/r" neutron-capture process either with features in-between the s and the r processes, or generated by superpositions of different neutron-capture processes in the same astrophysical site or in sites linked to each other - for example, in multiple stellar systems., Comment: 51 pages, 10 figures, 6 tables. Accepted for publication on The Astrophysical Journal

Chandra LETG+HRC-S and XMM-Newton RGS spectra of H-like C and N formed in the corona of the lambda And primary star, a mildly metal-poor G8 III-IV first ascent giant that completed dredge-up ~50 Myr ago, have been used to make a precise measurement of its surface C/N ratio. We obtain the formal result [C/N]=0.03+/-0.07, which is typical of old disk giants and in agreement with standard dredge-up theory for stars of about 1 M_sun or lower. In contrast, these stars as a group, including lambda And, have 12C/13C ~< 20, which is much lower than standard model predictions. We show that the abundances of the old disk giants are consistent with models including thermohaline mixing that begins at the red giant branch luminosity function "bump". Instead, lambda And indicates that the 12C/13C anomaly can be present immediately following dredge-up, contrary to current models of extra mixing on the red giant branch., Comment: To appear in the Astronomical 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

(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

GD 356 is unique among magnetic white dwarfs because it shows Zeeman-split Balmer lines in pure emission. The lines originate from a region of nearly uniform field strength (delta B/B is approximately 0.1) that covers 10 per cent of the stellar surface in which there is a temperature inversion. The energy source that heats the photosphere remains a mystery but it is likely to be associated with the presence of a companion. Based on current models we use archival Spitzer IRAC observations to place a new and stringent upper limit of 12 Jupiter masses for the mass of such a companion. In the light of this result and the recent discovery of a 115 min photometric period for GD 356, we exclude previous models that invoke accretion and revisit the unipolar inductor model that has been proposed for this system. In this model a highly conducting planet with a metallic core orbits the magnetic white dwarf and, as it cuts through field lines, a current is set flowing between the two bodies. This current dissipates in the photosphere of the white dwarf and causes a temperature inversion. Such a planet is unlikely to have survived the RGB/AGB phases of evolution so we argue that it may have formed from the circumstellar disc of a disrupted He or CO core during a rare merger of two white dwarfs. GD 356 would then be a white dwarf counterpart of the millisecond binary pulsar PSR 1257+12 which is known to host a planetary system., Comment: 9 pages, 4 figures, accepted by MNRAS

We examine the effects of thermohaline mixing on the composition of the envelopes of low-metallicity asymptotic giant branch (AGB) stars. We have evolved models of 1, 1.5 and 2 solar masses from the pre-main sequence to the end of the thermally pulsing asymptotic giant branch with thermohaline mixing applied throughout the simulations. In agreement with other authors, we find that thermohaline mixing substantially reduces the abundance of helium-3 on the upper part of the red giant branch in our lowest mass model. However, the small amount of helium-3 that remains is enough to drive thermohaline mixing on the AGB. We find that thermohaline mixing is most efficient in the early thermal pulses and its efficiency drops from pulse to pulse. Nitrogen is not substantially affected by the process, but we do see substantial changes in carbon-13. The carbon-12 to carbon-13 ratio is substantially lowered during the early thermal pulses but the efficacy of the process is seen to diminish rapidly. As the process stops after a few pulses, the carbon-12 to carbon-13 ratio is still able to reach values of 10^3-10^4, which is inconsistent with the values measured in carbon-enhanced metal-poor stars. We also note a surprising increase in the lithium-7 abundance, with log epsilon(Li-7) reaching values of over 2.5 in the 1.5 solar mass model. It is thus possible to get stars which are both C- and Li-rich at the same time. We compare our models to measurements of carbon and lithium in carbon-enhanced metal-poor stars which have not yet reached the giant branch. These models can simultaneously reproduced the observed C and Li abundances of carbon-enhanced metal-poor turn-off stars that are Li-rich, but the observed nitrogen abundances still cannot be matched., Comment: Accepted for publication in MNRAS. 12 pages, 7 figures

The carbon-enhanced metal-poor (CEMP) stars constitute approximately one fifth of the metal-poor ([Fe/H] ~< -2) population but their origin is not well understood. The most widely accepted formation scenario, invokes mass-transfer of carbon-rich material from a thermally-pulsing asymptotic giant branch (TPAGB) primary star to a less massive main-sequence companion which is seen today. Recent studies explore the possibility that an initial mass function biased toward intermediate-mass stars is required to reproduce the observed CEMP fraction in stars with metallicity [Fe/H] < -2.5. These models also implicitly predict a large number of nitrogen-enhanced metal-poor (NEMP) stars which is not seen. We investigate whether the observed CEMP and NEMP to extremely metal-poor (EMP) ratios can be explained without invoking a change in the initial mass function. We confirm earlier findings that with current detailed TPAGB models the large observed CEMP fraction cannot be accounted for. We find that efficient third dredge up in low-mass (less than 1.25Msun), low-metallicity stars may offer at least a partial explanation to the large observed CEMP fraction while remaining consistent with the small observed NEMP fraction., Comment: 20 pages, 23 figures, accepted for publication in A&A

We construct binary-star population nucleosynthesis models of carbon-enhanced metal poor (CEMP) stars. We compare the CEMP to EMP (extremely metal poor) ratio of our models to the observed ratio and find it is an order of magnitude too small. Through an increase in the efficiency of third dredge-up in low-mass, low-metallicity thermally-pulsing asymptotic-giant branch (TPAGB) stars our models better match the observations., Comment: Accepted for publication in PASA. Three pages, one figure

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

We have developed a detailed stellar evolution code capable of following the simultaneous evolution of both stars in a binary system, together with their orbital properties. To demonstrate the capabilities of the code we investigate potential progenitors for the Type IIb supernova 1993J, which is believed to have been an interacting binary system prior to its primary exploding. We use our detailed binary stellar evolution code to model this system to determine the possible range of primary and secondary masses that could have produced the observed characteristics of this system, with particular reference to the secondary. Using the luminosities and temperatures for both stars (as determined by Maund et al. 2004) and the remaining mass of the hydrogen envelope of the primary at the time of explosion, we find that if mass transfer is 100 per cent efficient the observations can be reproduced by a system consisting of a 15 solar mass primary and a 14 solar mass secondary in an orbit with an initial period of 2100 days. With a mass transfer efficiency of 50 per cent, a more massive system consisting of a 17 solar mass primary and a 16 solar mass secondary in an initial orbit of 2360 days is needed. We also investigate some of the uncertainties in the evolution, including the effects of tidal interaction, convective overshooting and thermohaline mixing., Comment: Accepted for publication in MNRAS. 11 pages, 5 figures

We investigate the behaviour of asymptotic giant branch (AGB) stars between metallicities Z = 10-4 and Z = 10-8 . We determine which stars undergo an episode of flash-driven mixing, where protons are ingested into the intershell convection zone, as they enter the thermally pulsing AGB phase and which undergo third dredge-up. We find that flash-driven mixing does not occur above a metallicity of Z = 10-5 for any mass of star and that stars above 2 M do not experience this phenomenon at any metallicity. We find carbon ingestion (CI), the mixing of carbon into the tail of hydrogen burning region, occurs in the mass range 2 M to around 4 M . We suggest that CI may be a weak version of the flash-driven mechanism. We also investigate the effects of convective overshooting on the behaviour of these objects. Our models struggle to explain the frequency of CEMP stars that have both significant carbon and nitrogen enhancement. Carbon can be enhanced through flash-driven mixing, CI or just third dredge up. Nitrogen can be enhanced through hot bottom burning and the occurrence of hot dredge-up also converts carbon into nitrogen. The C/N ratio may be a good indicator of the mass of the primary AGB stars., Comment: 15 pages, 13 figures, 1 table, accepted by MNRAS

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 model the evolution of the abundances of light elements in carbon-enhanced metal-poor (CEMP) stars, under the assumption that such stars are formed by mass transfer in a binary system. We have modelled the accretion of material ejected by an asymptotic giant branch star on to the surface of a companion star. We then examine three different scenarios: one in which the material is mixed only by convective processes, one in which thermohaline mixing is present and a third in which both thermohaline mixing and gravitational settling are taken in to account. The results of these runs are compared to light element abundance measurements in CEMP stars (primarily CEMP-s stars, which are rich in $s$-processes elements and likely to have formed by mass transfer from an AGB star), focusing on the elements Li, F, Na and Mg. None of the elements is able to provide a conclusive picture of the extent of mixing of accreted material. We confirm that lithium can only be preserved if little mixing takes place. The bulk of the sodium observations suggest that accreted material is effectively mixed but there are also several highly Na and Mg-rich objects that can only be explained if the accreted material is unmixed. We suggest that the available sodium data may hint that extra mixing is taking place on the giant branch, though we caution that the data is sparse., Comment: Accepted for publication in MNRAS. 9 figures, 1 table

We investigate the formation of carbon-enhanced metal-poor (CEMP) stars via the scenario of mass transfer from a carbon-rich asymptotic giant branch (AGB) primary to a low-mass companion in a binary system. We explore the extent to which material accreted from a companion star becomes mixed with that of the recipient, focusing on the effects of thermohaline mixing and gravitational settling. We have created a new set of asymptotic giant branch models in order to determine what the composition of material being accreted in these systems will be. We then model a range of CEMP systems by evolving a grid of models of low-mass stars, varying the amount of material accreted by the star (to mimic systems with different separations) and also the composition of the accreted material (to mimic accretion from primaries of different mass). We find that with thermohaline mixing alone, the accreted material can become mixed with between 16 and 88 per cent of the pristine stellar material of the accretor, depending on the mass accreted and the composition of the material. If we include the effects of gravitational settling, we find that thermohaline mixing can be inhibited and, in the case that only a small quantity of material is accreted, can be suppressed almost completely., Comment: Accepted for publication in MNRAS. 13 pages, 10 figures. The definitive version is available at www.blackwell-synergy.com

One possible scenario for the formation of carbon-enhanced metal-poor stars is the accretion of carbon-rich material from a binary companion which may no longer visible. It is generally assumed that the accreted material remains on the surface of the star and does not mix with the interior until first dredge-up. However, thermohaline mixing should mix the accreted material with the original stellar material as it has a higher mean molecular weight. We investigate the effect that this has on the surface abundances by modelling a binary system of metallicity Z=0.0001 with a 2 solar mass primary star and a 0.74 solar mass secondary star in an initial orbit of 4000 days. The accretion of material from the wind of the primary leads to the formation of a carbon-rich secondary. We find that the accreted material mixes fairly rapidly throughout 90% of the star, with important consequences for the surface composition. Models with thermohaline mixing predict very different surface abundances after first dredge-up compared to canonical models of stellar evolution., Comment: 4 pages, 5 figures. Accepted for publication in A&A letters

Galactic chemical evolution models which include sodium from type II supernovae alone underestimate the abundance of sodium in the interstellar medium by a factor of two to three over about 3 dex in metallicity and predict a flat behaviour in the evolution of [Na/Fe] at super-solar metallicities. Conversely, recent observations of stars with [Fe/H]~+0.4 suggest that [Na/Fe] increases at high metallicity. We have combined stellar evolution models of asymptotic giant branch and Wolf-Rayet stars with the latest supernova yields in an attempt to resolve these problems ... and have created many more., Comment: Poster given at the conference "Why Galaxies Care About AGB Stars", August 7-11 2006, University of Vienna, Austria; to appear in the proceedings (ASP Conference Series). Winner of best poster prize!

We investigate the uncertainty in surface abundances and yields of asymptotic giant branch (AGB) stars. We apply three different mass loss laws to a 1.5 solar mass star of metallicity Z=0.008 at the beginning of the thermally pulsing asymptotic giant branch (TP-AGB) phase. Efficient third dredge-up is found even at very low envelope mass, contrary to previous simulations with other evolution codes. We find that the yield of carbon is uncertain by about 15% and for most other light elements the yield is uncertain at the level of 20-80%. For iron group elements the uncertainty varies from around 30% for the more abundant species to over a factor of two for the less abundant radioactive species, like iron-60. The post-AGB surface abundances for this mass and metallicity are much more uncertain due to the dilution of dredged-up material in differing envelope masses in the later stages of the models. Our results are compared to known planetary nebula (PN) and post-AGB abundances. We find that the models are mostly consistent with observations but we are unable to reproduce observations of some of the isotopes., Comment: 12 pages, 5 figures. Accepted for publication in MNRAS

A version of the STARS stellar evolution code has been developed that uses a non-simultaneous solution of the equations of stellar structure and evolution. In all other respects it is identical to the normal, fully simultaneous version. It is therefore possible to test the dependence of the solution on how the equations are solved. Two cases are investigated: a 5Msol and a 3Msol star, both of metallicity Z=0.02. Prior to the asymptotic giant branch, the models are almost identical. However once thermal pulses start, the two methods of solution yield diverging results with the non-simultaneous technique predicting longer interpulse periods. This is traced to difficulties associated with hydrogen burning caused by the use of a moving mesh. It is shown that, with careful control of the temporal resolution, the results of the simultaneous technique can be recovered., Comment: 7 pages, 7 figures. Accepted for publication in MNRAS

A long standing problem with asymptotic giant branch (AGB) star models has been their inability to produce the low-luminosity carbon stars in the Large and Small Magellanic Clouds. Dredge-up must begin earlier and extend deeper. We find this for the first time in our models of LMC metallicity. Such features are not found in our models of SMC metallicity. The fully implicit and simultaneous stellar evolution code STARS has been used to calculate the evolution of AGB stars with metallicities of Z=0.008 and Z=0.004, corresponding to the observed metallicities of the Large and Small Magellanic Clouds, respecitively. Third dredge-up occurs in stars of 1Msol and above and carbon stars were found for models between 1Msol and 3Msol. We use the detailed models as input physics for a population synthesis code and generate carbon star luminosity functions. We now find that we are able to reproduce the carbon star luminosity function of the LMC without any manipulation of our models. The SMC carbon star luminosity function still cannot be produced from our detailed models unless the minimum core mass for third dredge-up is reduced by 0.06Msol., Comment: 6 pages, 5 figures. Accepted for publication in MNRAS