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.

In 1981, the idea of a superwind that ends the life of cool giant stars was proposed. Extreme OH/IR-stars develop superwinds with the highest mass-loss rates known so far, up to a few 10^(-4) Msun/yr, informing our understanding of the maximum mass-loss rate achieved during the Asymptotic Giant Branch (AGB) phase. A condundrum arises whereby the observationally determined duration of the superwind phase is too short for these stars to become white dwarfs. Here, we report on the detection of spiral structures around two cornerstone extreme OH/IR-stars, OH26.5+0.6 and OH30.1-0.7, identifying them as wide binary systems. Hydrodynamical simulations show that the companion's gravitational attraction creates an equatorial density enhancement mimicking a short extreme superwind phase, thereby solving the decades-old conundrum. This discovery restricts the maximum mass-loss rate of AGB stars around the single-scattering radiation-pressure limit of a few 10^(-5) Msun/yr. This brings about crucial implications for nucleosynthetic yields, planet survival, and the wind-driving mechanism.

Metal-poor stars in the Galactic halo often show strong enhancements in carbon and/or neutron-capture elements. However, the Galactic bulge is notable for its paucity of carbon-enhanced metal-poor (CEMP) and/or CH-stars, with only two such objects known to date. This begs the question whether the processes that produced their abundance distribution were governed by a comparable nucleosynthesis in similar stellar sites as for their more numerous counterparts in the halo. Recently, two contenders of such stars were discovered in the bulge, at [Fe/H] = $-1.5$ and $-$2.5 dex, both of which show enhancements in [C/Fe] of 0.4 and 1.4 dex, [Ba/Fe] in excess of 1.3 dex, and also elevated nitrogen. The more metal-poor of the stars is matched by standard $s$-process nucleosynthesis in low-mass Asymptotic Giant Branch (AGB) polluters. The other star shows an abnormally high [Rb/Fe] ratio. Here, we investigate the origin of the abundance peculiarities in the Rb-rich star by new, detailed measurements of heavy element abundances and by comparing the chemical element ratios of 36 species to models of neutron-capture nucleosynthesis. The $i$-process with intermediate neutron densities between those of the $s$- and $r$-processes has been previously found to provide good matches of CEMP stars with enhancements in both $r$- and $s$-process elements, rather than invoking a superposition of yields from the respective individual processes. However, the peculiar bulge star is incompatible with a pure $i$-process from a single ingestion event. Instead, it can, statistically, be better reproduced by models accounting for two proton ingestion events, or by an $i$-process component in combination with $s$-process nucleosynthesis in low-to-intermediate mass AGB stars, indicating multiple polluters. [abridged]

The progenitors of many core-collapse supernovae (CCSNe) are expected to be in binary systems. By performing a series of three-dimensional hydrodynamical simulations, we investigate how CCSN explosions affect their binary companion. We find that the amount of removed stellar mass, the resulting impact velocity, and the chemical contamination of the companion that results from the impact of the SN ejecta, strongly increases with decreasing binary separation and increasing explosion energy. Also, it is foud that the impact effects of CCSN ejecta on the structure of main-sequence (MS) companions, and thus their long term post-explosion evolution, is in general not be dramatic.

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

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.

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.

Carbon-enhanced metal-poor stars with s-process enrichment (CEMP-s) are believed to be the products of mass transfer from an AGB companion, which has long since become a white dwarf. The surface abundances of CEMP-s stars are thus commonly assumed to reflect the nucleosynthesis output of the first AGB stars. We have previously shown that, for this to be the case, some physical mechanism must counter atomic diffusion in these nearly fully radiative stars, which otherwise leads to surface abundance anomalies clearly inconsistent with observations. Here we take into account angular momentum accretion by these stars. We compute in detail the evolution of typical CEMP-s stars from the ZAMS, through the mass accretion, and up the RGB for a wide range of specific angular momentum of the accreted material, corresponding to rotation velocities between about 0.3 and 300 km/s. We find that only for specific angular momentum above 1e+17 cm2/s (rotation velocities above 20 km/s) angular momentum accretion directly causes chemical dilution of the accreted material. This could nevertheless be relevant to CEMP-s stars, which are observed to rotate more slowly, if they undergo continuous angular momentum loss akin to solar-like stars. In models with rotation velocities characteristic of CEMP-s stars, rotational mixing primarily serves to inhibit atomic diffusion, such that the maximal surface abundance variations (with respect to the composition of the accreted material) prior to first dredge-up remain within about 0.4 dex without thermohaline mixing or about 0.5-1.5 dex with thermohaline mixing. Even in models with the lowest rotation velocities (under a km/s), rotational mixing is able to severely inhibit atomic diffusion, compared to non-rotating models. We thus conclude that it offers a natural solution to the problem posed by atomic diffusion and cannot be neglected in models of CEMP-s stars.

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.

The chemically peculiar barium stars, CH stars, and most CEMP stars are all believed to be the products of mass transfer in binary systems from a now extinct AGB primary star. The mass of the AGB star and the orbital parameters of the system are the key factors usually considered when determining how much mass is transferred onto the lower-mass main-sequence companion. What is usually neglected, however, is the angular momentum of the accreted material, which should spin up the accreting star. If the star reaches critical rotation, further accretion should cease until the excess angular momentum is somehow dealt with. If the star cannot redistribute or lose the angular momentum while the primary is on the AGB, the amount of mass accreted could be much lower than otherwise expected. Here we present calculations, based on detailed stellar evolution models, of the mass that can be accreted by putative progenitors of Ba and CEMP stars before they reach critical rotation under the assumption that no angular momentum loss occurs during the mass transfer. We consider different accretion rates and values of specific angular momentum. The most stringent limits on the accreted masses result from considering accretion from a Keplerian accretion disk, which is likely present during the formation of most extrinsically-polluted carbon-enriched stars. Our calculations indicate that in this scenario only about 0.05 solar masses of material can be added to the accreting star before it reaches critical rotation, which is much too low to explain the chemical enrichment of many Ba and CEMP stars. Either the specific angular momentum of the accreted material has to effectively be lower by about a factor of ten than the Keplerian value, or significant angular momentum losses must occur for substantial accretion to take place.

The production of the elements heavier than iron via slow neutron captures (the s process) is a main feature of the contribution of asymptotic giant branch (AGB) stars of low mass (< 5 Msun) to the chemistry of the cosmos. However, our understanding of the main neutron source, the 13C(alpha,n)16O reaction, is still incomplete. It is commonly assumed that in AGB stars mixing beyond convective borders drives the formation of 13C pockets. However, there is no agreement on the nature of such mixing and free parameters are present. By means of a parametric model we investigate the impact of different mixing functions on the final s-process abundances in low-mass AGB models. Typically, changing the shape of the mixing function or the mass extent of the region affected by the mixing produce the same results. Variations in the relative abundance distribution of the three s-process peaks (Sr, Ba, and Pb) are generally within +/-0.2 dex, similar to the observational error bars. We conclude that other stellar uncertainties - the effect of rotation and of overshoot into the C-O core - play a more important role than the details of the mixing function. The exception is at low metallicity, where the Pb abundance is significantly affected. In relation to the composition observed in stardust SiC grains from AGB stars, the models are relatively close to the data only when assuming the most extreme variation in the mixing profile.

The aim of this paper is to investigate the $^{17}$O/$^{18}$O ratio for a sample of AGB stars, containing M-, S- and C-type stars. These ratios are evaluated in relation to fundamental stellar evolution parameters: the stellar initial mass and pulsation period. Circumstellar $^{13}$C$^{16}$O, $^{12}$C$^{17}$O and $^{12}$C$^{18}$O line observations were obtained for a sample of nine stars with various single-dish long-wavelength facilities. Line intensity ratios are shown to relate directly to the surface $^{17}$O/$^{18}$O abundance ratio. Stellar evolution models predict the $^{17}$O/$^{18}$O ratio to be a sensitive function of initial mass and to remain constant throughout the entire TP-AGB phase for stars initially less massive than 5\,$M_{\odot}$. This makes the measured ratio a probe of the initial stellar mass. Observed $^{17}$O/$^{18}$O ratios are found to be well in the range predicted by stellar evolution models that do not consider convective overshooting. From this, accurate initial mass estimates are calculated for seven sources. For the remaining two sources two mass solutions result, though with a larger probability that the low-mass solution is the correct one. Finally, hints at a possible separation between M/S- and C-type stars when comparing the $^{17}$O/$^{18}$O ratio to the stellar pulsation period are presented.

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.

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.

A significant fraction of all metal-poor stars are carbon-rich. Most of these carbon-enhanced metal-poor (CEMP) stars also show enhancement in elements produced mainly by the s-process (CEMP-s stars) and evidence suggests that the origin of these non-standard abundances can be traced to mass transfer from a binary asymptotic giant branch (AGB) companion. Thus, observations of CEMP-s stars are commonly used to infer the nucleosynthesis output of low-metallicity AGB stars. A crucial step in this exercise is understanding what happens to the accreted material after mass transfer ceases. Here we present models of the post-mass-transfer evolution of CEMP-s stars considering the physics of thermohaline mixing and atomic diffusion, including radiative levitation. We find that stars with typical CEMP-s star masses (M ~ 0.85 Msun) have very shallow convective envelopes (Menv < 1e-7 Msun). Hence, the surface abundance variations arising from the competition between gravitational settling and radiative levitation should be orders of magnitude larger than observed (e.g. [C/Fe]<-1 or [C/Fe]>+4). We are therefore unable to reproduce the spread in the observed abundances with these models and conclude that some other physical process must largely suppress atomic diffusion in the outer layers of CEMP-s stars. We demonstrate that this could be achieved by some additional (turbulent) mixing process operating at the base of the convective envelope, as found by other authors. Alternatively, mass-loss rates around 1e-13 Msun/yr could also negate most of the abundance variations by eroding the surface layers and forcing the base of the convective envelope to move inwards in mass. Since atomic diffusion cannot have a substantial effect on the surface abundances of CEMP-s stars, the dilution of the accreted material, while variable in degree from one star to the next, is most likely the same for all elements.

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.

We assess the systematic uncertainties in stellar evolutionary calculations for low- to intermediate-mass, main-sequence stars. We compare published stellar tracks from several different evolution codes with our own tracks computed using the stellar codes STARS and MESA. In particular, we focus on tracks of 1 and 3 solar masses at solar metallicity. We find that the spread in the available 1 solar mass tracks (computed before the recent solar composition revision by Asplund et al.) can be covered by tracks between 0.97-1.01 solar masses computed with the STARS code. We assess some possible causes of the origin of this uncertainty, including how the choice of input physics and the solar constraints used to perform the solar calibration affect the tracks. We find that for a 1 solar mass track, uncertainties of around 10% in the initial hydrogen abundance and initial metallicity produce around a 2% error in mass. For the 3 solar mass tracks, there is very little difference between the tracks from the various different stellar codes. The main difference comes in the extent of the main sequence, which we believe results from the different choices of the implementation of convective overshooting in the core. Uncertainties in the initial abundances lead to a 1-2% error in the mass determination. These uncertainties cover only part of the total error budget, which should also include uncertainties in the input physics (e.g., reaction rates, opacities, convective models) and any missing physics (e.g., radiative levitation, rotation, magnetic fields). Uncertainties in stellar surface properties such as luminosity and effective temperature will further reduce the accuracy of any potential mass determinations.

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.

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.

We present late-time Hubble Space Telescope (HST) ultraviolet (UV) and optical observations of the site of SN 2011dh in the galaxy M51, ~1164 days post-explosion. At the SN location, we observe a point source that is visible at all wavelengths, that is significantly fainter than the spectral energy distribution (SED) of the Yellow Supergiant progenitor observed prior to explosion. The previously reported photometry of the progenitor is, therefore, completely unaffected by any sources that may persist at the SN location after explosion. In comparison with the previously reported late-time photometric evolution of SN 2011dh, we find that the light curve has plateaued at all wavelengths. The SED of the late-time source is clearly inconsistent with a SED of stellar origin. Although the SED is bright at UV wavelengths, there is no strong evidence that the late-time luminosity originates solely from a stellar source corresponding to the binary companion, although a partial contribution to the observed UV flux from a companion star can not be ruled out.

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.

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.

The stellar population in the Galactic halo is characterised by a large fraction of CEMP stars. Most CEMP stars are enriched in $s$-elements (CEMP-$s$ stars), and some of these are also enriched in $r$-elements (CEMP-$s/r$ stars). One formation scenario proposed for CEMP stars invokes wind mass transfer in the past from a TP-AGB primary star to a less massive companion star which is presently observed. We generate low-metallicity populations of binary stars to reproduce the observed CEMP-star fraction. In addition, we aim to constrain our wind mass-transfer model and investigate under which conditions our synthetic populations reproduce observed abundance distributions. We compare the CEMP fractions and the abundance distributions determined from our synthetic populations with observations. Several physical parameters of the binary stellar population of the halo are uncertain, e.g. the initial mass function, the mass-ratio and orbital-period distributions, and the binary fraction. We vary the assumptions in our model about these parameters, as well as the wind mass-transfer process, and study the consequent variations of our synthetic CEMP population. The CEMP fractions calculated in our synthetic populations vary between 7% and 17%, a range consistent with the CEMP fractions among very metal-poor stars recently derived from the SDSS/SEGUE data sample. The results of our comparison between the modelled and observed abundance distributions are different for CEMP-$s/r$ stars and for CEMP-$s$ stars. For the latter, our simulations qualitatively reproduce the observed distributions of C, Na, Sr, Ba, Eu, and Pb. Contrarily, for CEMP-$s/r$ stars our model cannot reproduce the large abundances of neutron-rich elements such as Ba, Eu, and Pb. This result is consistent with previous studies, and suggests that CEMP-$s/r$ stars experienced a different nucleosynthesis history to CEMP-$s$ stars.

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.

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

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.

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.

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

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.

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.

Calculations from stellar evolutionary models of low- and intermediate-mass asymptotic giant branch (AGB) stars provide predictions of elemental abundances and yields for comparison to observations. However, there are many uncertainties that reduce the accuracy of these predictions. One such uncertainty involves the treatment of low-temperature molecular opacities that account for the surface abundance variations of C, N, and O. A number of prior calculations of intermediate-mass AGB stellar models that incorporate both efficient third dredge-up and hot bottom burning include a molecular opacity treatment which does not consider the depletion of C and O due to hot bottom burning. Here we update the molecular opacity treatment and investigate the effect of this improvement on calculations of intermediate-mass AGB stellar models. We perform tests on two masses, 5 M$_{\odot}$ and 6 M$_{\odot}$, and two metallicities, $Z~=~0.001$ and $Z~=~0.02$, to quantify the variations between two opacity treatments. We find that several evolutionary properties (e.g. radius, $T_{\rm eff}$ and $T_{\rm bce}$) are dependent on the opacity treatment. Larger structural differences occur for the $Z~=~0.001$ models compared to the $Z~=~0.02$ models indicating that the opacity treatment has a more significant effect at lower metallicity. As a consequence of the structural changes, the predictions of isotopic yields are slightly affected with most isotopes experiencing changes up to 60 per cent for the $Z~=~0.001$ models and 20 per cent for the $Z~=~0.02$ models. Despite this moderate effect, we conclude that it is more fitting to use variable molecular opacities for models undergoing hot bottom burning.

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.

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.

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.

To understand the chemical evolution of the Galaxy, we need to understand the contribution of PNe in the 3He abundance. 3He abundances have been detected only in a couple of PNe, their abundances are consistent with the standard models, in which 3He is produced in significant quantities by stars of 1-2M_sun. However, for all the other PNe observed to date there have been no detections, therefore only upper limits in their abundance can be calculated. Observations of the 3He+ emission line using the VLA towards the PNe IC 418, NGC 6572 and NGC 7009 were used to obtain upper limits for their 3He abundance. Because the abundance of 3He in HII regions, the ISM and the proto-solar system is much lower than what is predicted, new chemical models were proposed. The resulting evolution of 3He, based on stellar evolution models, can be consistent with the values determined in pre-solar material and the ISM if 96 percent of the population of stars with mass below 2.5M_sun has undergone enhanced 3He depletion. This implies that unless the combined sample of PNe that has been observed so far is very atypical, the solution to the "The 3He Problem" lies elsewhere. However, the results presented here suggest that more observations are needed in order to make a strong conclusion about the stellar evolution models.

The handful of available observations of AGB stars in Galactic Globular Clusters suggest that the GC AGB populations are dominated by cyanogen-weak stars. This contrasts strongly with the distributions in the RGB (and other) populations, which generally show a 50:50 bimodality in CN band strength. If it is true that the AGB populations show very different distributions then it presents a serious problem for low mass stellar evolution theory, since such a surface abundance change going from the RGB to AGB is not predicted by stellar models. However this is only a tentative conclusion, since it is based on very small AGB sample sizes. To test whether this problem really exists we have carried out an observational campaign specifically targeting AGB stars in GCs. We have obtained medium resolution spectra for about 250 AGB stars across 9 Galactic GCs using the multi-object spectrograph on the AAT (2df/AAOmega). We present some of the preliminary findings of the study for the second parameter trio of GCs: NGC 288, NGC 362 and NGC 1851. The results indeed show that there is a deficiency of stars with strong CN bands on the AGB. To confirm that this phenomenon is robust and not just confined to CN band strengths and their vagaries, we have made observations using FLAMES/VLT to measure elemental abundances for NGC 6752.We present some initial results from this study also. Our sodium abundance results show conclusively that only a subset of stars in GCs experience the AGB phase of evolution. This is the first direct, concrete confirmation of the phenomenon.

Our understanding of magnetic fields in late-type stars is strongly driven by what we know of the solar magnetic field. For this reason, it is crucial to understand how typical the solar dynamo is. To do this we need to compare the solar magnetic field with that of other stars as similar to the Sun as possible, both in stellar parameters and age, hence activity. We present here the detection of a magnetic field in three planet-hosting solar-like stars having a mass, age, and activity level comparable to that of the Sun. We used the HARPSpol spectropolarimeter to obtain high-resolution high-quality circularly polarised spectra of HD 70642, HD 117207, and HD 154088, using the Least-Squares Deconvolution technique to detect the magnetic field. From the Stokes I spectra, we calculated the logR activity index for each star. We compared the position of the stars in the Hertzsprung-Russell diagram to evolutionary tracks, to estimate their mass and age. We used the lithium abundance, derived from the Stokes I spectra, to further constrain the ages. We obtained a definite magnetic field detection for both HD 70642 and HD 154088, while for HD 117207 we obtained a marginal detection. Due to the lower signal-to-noise ratio of the observations, we were unable to detect the magnetic field in the second set of observations available for HD 117207 and HD 154088. On the basis of effective temperature, mass, age, and activity level the three stars can be considered solar analogs. HD 70642, HD 117207, and HD 154088 are ideal targets for a comparative study between the solar magnetic field and that of solar analogs.

The Galactic globular cluster NGC 1851 has raised much interest since HST photometry revealed that it hosts a double subgiant branch. Here we report on our homogeneous study into the cyanogen (CN) bandstrengths in the RGB population (17 stars) and AGB population (21 stars) using AAOmega/2dF spectra with R $\sim 3000$. We discover that NGC 1851 hosts a quadrimodal distribution of CN bandstrengths in its RGB and AGB populations. This result supports the merger formation scenario proposed for this cluster, such that the CN quadrimodality could be explained by the superposition of two `normal' bimodal populations. A small sample overlap with an abundance catalogue allowed us to tentatively explore the relationship between our CN populations and a range of elemental abundances. We found a striking correlation between CN and [O/Na]. We also found that the four CN peaks may be paired -- the two CN-weaker populations being associated with low Ba and the two CN-stronger populations with high Ba. If true then s-process abundances would be a good diagnostic for disentangling the two original clusters in the merger scenario. More observations are needed to confirm the quadrimodality, and also the relationship between the subpopulations. We also report CN results for NGC 288 as a comparison. Our relatively large samples of AGB stars show that both clusters have a bias towards CN-weak AGB populations.

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.

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.

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.

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.

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.

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.

We provide a brief review of thermohaline physics and why it is a candidate extra mixing mechanism during the red giant branch (RGB). We discuss how thermohaline mixing (also called $\delta$ $\mu$ mixing) during the RGB due to helium-3 burning, is more complicated than the operation of thermohaline mixing in other stellar contexts (such as following accretion from a binary companion). We try to use observations of carbon depletion in globular clusters to help constrain the formalism and the diffusion coefficient or mixing velocity that should be used in stellar models. We are able to match the spread of carbon depletion for metal poor field giants but are unable to do so for cluster giants, which may show evidence of mixing prior to even the first dredge-up event.

A previously reported literature search suggested that the AGB stars in Galactic globular clusters may be showing different distributions of CN-strong and CN-weak stars as compared to their RGB stars. In most cases the second giant branches of GCs appeared to be deficient in stars with strong CN bands. However the sample sizes of AGB stars at that time were too small to give a definitive picture. Thus an observing campaign targeting GC AGB stars was proposed. We now have medium resolution spectral observations of about 250 GC AGB stars across 9 globular clusters, obtained with the 2dF/AAOmega instrument on the Anglo-Australian Telescope. In this paper we report some preliminary results regarding the distributions of CN-strong and CN-weak stars on the two giant branches of a selection of globular clusters. We find that some GCs show a total lack of CN-strong stars on the AGB, whilst some show a reduction in CN-strong stars as compared to the RGB. Standard stellar evolution does not predict this change in surface abundance between the two giant branches. We discuss some possible causes of this unexpected phenomenon.

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

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.

Diffusion of atoms can be important during quiescent phases of stellar evolution. Particularly in the very thin inert envelopes of subdwarf B stars, diffusive movements will considerably change the envelope structure and the surface abundances on a short timescale. Also, the subdwarfs will inherit the effects of diffusion in their direct progenitors, namely giants near the tip of the red giant branch. This will influence the global evolution and the pulsational properties of subdwarf B stars. We investigate the impact of gravitational settling, thermal diffusion and concentration diffusion on the evolution and pulsations of subdwarf B stars. Our diffusive stellar models are compared with models evolved without diffusion. We constructed subdwarf B models with a mass of 0.465 Msun from a 1 and 3 Msun ZAMS progenitor. The low mass star ignited helium in an energetic flash, while the intermediate mass star started helium fusion gently. For each progenitor type we computed series with and without atomic diffusion. Atomic diffusion in red giants causes the helium core mass at the onset of helium ignition to be larger. We find an increase of 0.0015 Msun for the 1 Msun model and 0.0036 Msun for the 3 Msun model. The effects on the red giant surface abundances are small after the first dredge up. The evolutionary tracks of the diffusive subdwarf B models are shifted to lower surface gravities and effective temperatures due to outward diffusion of hydrogen. This affects both the frequencies of the excited modes and the overall frequency spectrum. Especially the structure and pulsations of the post-non-degenerate sdB star are drastically altered, proving that atomic diffusion cannot be ignored in these stars.

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.

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.

We present the results of binary population simulations of carbon- and nitrogen-enhanced metal-poor (CEMP and NEMP) stars. We show that the observed paucity of very nitrogen-rich stars puts strong constraints on possible modifications of the initial mass function at low metallicity.

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.

We present MONTAGE, a post-processing nucleosynthesis code that combines a traditional network for isotopes lighter than calcium with a rapid algorithm for calculating the s-process nucleosynthesis of the heavier isotopes. The separation of those parts of the network where only neutron-capture and beta-decay reactions are significant provides a substantial advantage in computational efficiency. We present the yields for a complete set of s-process isotopes for a 3 Mo, Z = 0.02 stellar model, as a demonstration of the utility of the approach. Future work will include a large grid of models suitable for use in calculations of Galactic chemical evolution.

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.

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.

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.

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.

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.

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.

We use the Cambridge stellar evolution code STARS to model the evolution of 5-7 solar mass zero-metallicity stars. With enhanced resolution at the hydrogen and helium burning shell in the AGB phases, we are able to model the entire thermally pulsing asymptotic giant branch (TP-AGB) phase. The helium luminosities of the thermal pulses are significantly lower than in higher metallicity stars so there is no third dredge-up. The envelope is enriched in nitrogen by hot-bottom burning of carbon that was previously mixed in during second dredge-up. There is no s-process enrichment owing to the lack of third dredge up. The thermal pulses grow weaker as the core mass increases and they eventually cease. From then on the star enters a quiescent burning phase which lasts until carbon ignites at the centre of the star when the CO core mass is 1.36 solar mass. With such a high degeneracy and a core mass so close to the Chandrasekhar mass, we expect these stars to explode as type 1.5 supernovae, very similar to Type Ia supernovae but inside a hydrogen rich envelope.

We use the Cambridge stellar evolution code STARS to model the evolution and nucleosynthesis of zero-metallicity intermediate-mass stars. We investigate the effect of duplicity on the nucleosynthesis output of these systems and the potential abundances of the secondaries. The surfaces of zero-metallicity stars are enriched in CNO elements after second dredge up. During binary interaction, such as Roche lobe overflow or wind accretion, metals can be released from these stars and the secondaries enriched in CNO isotopes. We investigate the formation of the two most metal poor stars known, HE 0107-5240 and HE 1327-2326. The observed carbon and nitrogen abundances of HE 0107-5240 can be reproduced by accretion of material from the companion-enhanced wind of a seven solar star after second dredge-up, though oxygen and sodium are underproduced. We speculate that HE 1327-2326, which is richer in nitrogen and strontium, may similarly be formed by wind accretion in a later AGB phase after third dredge-up.

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.

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.

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.

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.

In the light of recent recalculations of the \el{19}{F}(alpha,p)\el{22}{Ne} reaction rate we present results of the expected yield of \el{19}{F} from Wolf-Rayet (WR) stars. We have computed models using the upper and lower limits for the rate in addition to the recommended rate and hence we constrain the uncertainty in the yield with respect to this reaction. We find a yield of 3.1x10^{-4}M_sol of \el{19}{F} with our recommended rate and a difference of a factor of two between the yields computed with the upper and lower limits. In comparison with previous work we find a difference in the yield of approximately a factor of 4, connected to a different choice of mass loss. Model uncertainties must be carefully evaluated in order to obtain a reliable estimate of the yield of fluorine from WR stars together with its uncertainties.

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.

We present results of the evolution of AGB stars of 3M_sol and 5M_sol with solar metallicity calculated with the Eggleton stellar evolution code (STARS), which has a fully implicit and simultaneous method for solving for the stellar structure, convective mixing and nuclear burning. We introduce the concept of a viscous mesh in order to improve the numerical stability of the calculations. For the 5M_sol star, we evolve through 25 thermal pulses and their associated third dredge-up events. We obtain a maximum helium luminosity of 1.7x10^9 L_sol and significantly deep dredge-up after the second pulse. Strong hot-bottom burning is observed after the 5th pulse. The 3M_sol model is evolved through 20 thermal pulse events and we find third dredge-up after the 7th pulse. During the 14th pulse sufficient carbon has been brought to the surface to produce a carbon star. We find that dredge-up and the transformation into a carbon star occur at significantly smaller core masses (0.584M_sol and 0.608M_sol, respectively) than in previous calculations for 3M_sol.