Search

Precise masses of red-giant stars enable a robust inference of their ages, but there are cases where these age estimates are highly precise yet very inaccurate. Examples are core-helium-burning (CHeB) stars that have lost more mass than predicted by standard single-star evolutionary models. Members of star clusters in the ${\it Kepler}$ database represent a unique opportunity to identify such stars, because they combine exquisite asteroseismic constraints with independent age information. In this work we focus on the single, metal-rich, Li-rich, low-mass, CHeB star KIC4937011, which is a member of the open cluster NGC 6819 (turn-off mass of $\approx 1.6 \, M_\odot$, i.e. age of $\approx 2.4$ Gyr). This star has $\approx 1 \, M_\odot$ less mass than expected for its age and metallicity, which could be explained by binary interactions or mass-loss along the red-giant branch (RGB). To infer formation scenarios for this object, we perform a Bayesian analysis by combining the binary stellar evolutionary framework $\texttt{binary_c v2.2.3}$ with the dynamic nested sampling approach contained in the $\texttt{dynesty v2.1.1}$ package. We find that this star is likely the result of a common-envelope evolution (CEE) phase during the RGB stage of the primary star in which the low-mass ($<0.71 \, M_\odot$) main sequence companion does not survive. During the CEE phase $\approx 1 \, M_\odot$ of material is ejected from the system, and the final star reaches the CHeB stage after helium flashes as if it were a single star of mass $\approx 0.7 \, M_\odot$, which is what we observe today. Although the proposed scenario is consistent with photometric and spectroscopic observations, a quantitative comparison with detailed stellar evolution calculations is needed to quantify the systematic skewness of radius, luminosity, and effective temperature distributions towards higher values than observations., Comment: Accepted for publication by A&A

Context. The explosive burning that drives nova eruptions results in unique nucleosynthesis that heavily over-produces certain isotopes relative to the solar abundance. However, novae are often ignored when considering the chemical evolution of our Galaxy due to their low ejecta masses. Aims. In this work, we use previously computed synthetic nova populations and the galactic chemical evolution code OMEGA+ to assess the impact that novae have on the evolution of stable elemental and isotopic abundances. Methods. We combine populations of novae computed using the binary population synthesis code binary_c with the galactic chemical evolution code OMEGA+ and detailed, white dwarf mass-dependent nova yields to model the nucleosynthetic contributions of novae to the evolution of the Milky Way. We consider three different nova yield profiles, each corresponding to a different set of nova yield calculations. Results. Despite novae from low-mass white dwarfs (WDs) dominating nova ejecta contributions, we find that novae occurring on massive WDs are still able to contribute significantly to many isotopes, particularly those with high mass numbers. We find that novae can produce up to 35% of the Galactic 13C and 15N mass by the time the model Galaxy reaches [Fe/H] = 0, and earlier in the evolution of the Galaxy (between [Fe/H] = -2 and -1) novae may have been the dominant source of 15N. Predictions for [13C/Fe], [15N/Fe], 12C/13C, and 14N/15N abundances ratios vary by up to 0.2 dex at [Fe/H] = 0 and by up to 0.7 dex in [15N/Fe] and 14N/15N between [Fe/H] = -2 and -1 (corresponding approximately to Galactic ages of 170 Myr and 1 Gyr in our model). The Galactic evolution of other stable isotopes (excluding Li) is not noticeably affected by including novae.

We develop a rapid algorithm for the evolution of stable, circular, circumbinary discs suitable for parameter estimation and population synthesis modelling. Our model includes disc mass and angular momentum changes, accretion on to the binary stars, and binary orbital eccentricity pumping. We fit our model to the post-asymptotic giant branch (post-AGB) circumbinary disc around IRAS 08544-4431, finding reasonable agreement despite the simplicity of our model. Our best-fitting disc has a mass of about $0.01\, \mathrm{M}_{\odot }$ and angular momentum $2.7\times 10^{52}\, \mathrm{g}\, \mathrm{cm}^{2}\, \mathrm{s}^{-1}\simeq 9 \,\mathrm{M}_{\odot }\, \mathrm{km}\, \mathrm{s}^{-1}\, \mathrm{au}$, corresponding to 0.0079 and 0.16 of the common-envelope mass and angular momentum, respectively. The best-fitting disc viscosity is $\alpha _\mathrm{disc} = 5 \times 10^{-3}$ and our tidal torque algorithm can be constrained such that the inner edge of the disc $R_{\mathrm{in}}\sim 2a$. The inner binary eccentricity reaches about 0.13 in our best-fitting model of IRAS 08544-4431, short of the observed 0.22. The circumbinary disc evaporates quickly when the post-AGB star reaches a temperature of $\sim \! 6\times 10^4\, \mathrm{K}$, suggesting that planetismals must form in the disc in about $10^{4}\, \mathrm{yr}$ if secondary planet formation is to occur, while accretion from the disc on to the stars at about 10 times the inner-edge viscous rate can double the disc lifetime., Comment: Accepted for publication in Monthly Notices of the Royal Astronomical Society

A common requirement in science is to store and share large sets of simulation data in an efficient, nested, flexible and human-readable way. Such datasets contain number counts and distributions, i.e. histograms and maps, of arbitrary dimension and variable type, e.g. floating-point number, integer or character string. Modern high-level programming languages like Perl and Python have associated arrays, knowns as dictionaries or hashes, respectively, to fulfil this storage need. Low-level languages used more commonly for fast computational simulations, such as C and Fortran, lack this functionality. We present libcdict, a C dictionary library, to solve this problem. Libcdict provides C and Fortran application programming interfaces (APIs) to native dictionaries, called cdicts, and functions for cdicts to load and save these as JSON and hence for easy interpretation in other software and languages like Perl, Python and R., Comment: Accepted for publication in JOSS (The Journal of Open-Source Software)

We present detailed implementations of (a) binary stellar evolution (using binary_c) and (b) dust production and destruction into the cosmological semi-analytic galaxy evolution simulation, L-Galaxies. This new version of L-Galaxies is compared to a version assuming only single stars and to global and spatially-resolved observational data across a range of redshifts ($z$). We find that binaries have a negligible impact on the stellar masses, gas masses, and star formation rates of galaxies only if the total mass ejected by massive stars is unchanged. This is because massive stars determine the strength of supernova (SN) feedback, which in turn regulates galaxy growth. Binary effects, such as common envelope ejection and novae, affect carbon and nitrogen enrichment in galaxies, however heavier alpha elements are more affected by the choice of SN and wind yields. Unlike many other simulations, the new L-Galaxies reproduces observed dust-to-metal (DTM) and dust-to-gas (DTG) ratios at $z\sim{}0-4$. This is mainly due to shorter dust accretion timescales in dust-rich environments. However, dust masses are under-predicted at $z>4$, highlighting the need for enhanced dust production at early times in simulations, possibly accompanied by increased star formation. On sub-galactic scales, there is very good agreement between L-Galaxies and observed dust and metal radial profiles at $z=0$. A drop in DTM ratio is also found in diffuse, low-metallicity regions, contradicting the assumption of a universal value. We hope that this work serves as a useful template for binary stellar evolution implementations in other cosmological simulations in future., Comment: 19 pages, 10 figures. Submitted to MNRAS. Comments welcome!

Aluminium-26 is a radioactive isotope which can be synthesized within asymptotic giant branch (AGB) stars, primarily through hot bottom burning. Studies exploring $^{26}$Al production within AGB stars typically focus on single-stars; however, observations show that low- and intermediate-mass stars commonly exist in binaries. We use the binary population synthesis code binary_c to explore the impact of binary evolution on $^{26}$Al yields at solar metallicity both within individual AGB stars and a low/intermediate-mass stellar population. We find the key stellar structural condition achieving most $^{26}$Al overproduction is for stars to enter the thermally-pulsing AGB (TP-AGB) phase with small cores relative to their total masses, allowing those stars to spend abnormally long times on the TP-AGB compared to single-stars of identical mass. Our population with a binary fraction of 0.75 has an $^{26}$Al weighted population yield increase of $25\%$ compared to our population of only single-stars. Stellar-models calculated from the Mt Stromlo/Monash Stellar Structure Program, which we use to test our results from binary_c and closely examine the interior structure of the overproducing stars, support our binary_c results only when the stellar envelope gains mass after core-He depletion. Stars which gain mass before core-He depletion still overproduce $^{26}$Al, but to a lesser extent. This introduces some physical uncertainty into our conclusions as $55\%$ of our $^{26}$Al overproducing stars gain envelope mass through stellar wind accretion onto pre-AGB objects. Our work highlights the need to consider binary influence on the production of $^{26}$Al., Comment: 20 pages, 17 figures, and 6 tables. This article has been accepted for publication in MNRAS Published by Oxford University Press on behalf of the Royal Astronomical Society

X-ray binaries (XRBs) are thought to regulate cosmic thermal and ionization histories during the Epoch of Reionization and Cosmic Dawn ($z\sim 5-30$).Theoretical predictions of the X-ray emission from XRBs are important for modelling such early cosmic evolution. Nevertheless, the contribution from Be-XRBs, powered by accretion of compact objects from decretion disks around rapidly rotating O/B stars, has not been investigated systematically. Be-XRBs are the largest class of high-mass XRBs (HMXBs) identified in local observations and are expected to play even more important roles in metal-poor environments at high redshifts. In light of this, we build a physically motivated model for Be-XRBs based on recent hydrodynamic simulations and observations of decretion disks. Our model is able to reproduce the observed population of Be-XRBs in the Small Magellanic Cloud with appropriate initial conditions and binary stellar evolution parameters. We derive the X-ray output from Be-XRBs as a function of metallicity in the (absolute) metallicity range $Z\in [10^{-4},0.03]$ with a large suite of binary population synthesis (BPS) simulations. The simulated Be-XRBs can explain a non-negligible fraction ($\gtrsim 30\%$) of the total X-ray output from HMXBs observed in nearby galaxies for $Z\sim 0.0003-0.02$. The X-ray luminosity per unit star formation rate from Be-XRBs in our fiducial model increases by a factor of $\sim 8$ from $Z=0.02$ to $Z=0.0003$, which is similar to the trend seen in observations of all types of HMXBs. We conclude that Be-XRBs are potentially important X-ray sources that deserve greater attention in BPS of XRBs., Comment: 26 pages, 16+7 figures, accepted for publication in MNRAS, added discussions on supernova kicks, stellar winds and treatments of super-Eddington accretion and angular momentum loss during mass transfer

Binary stars evolve into chemically-peculiar objects and are a major driver of the Galactic enrichment of heavy elements. During their evolution they undergo interactions, including tides, that circularize orbits and synchronize stellar spins, impacting both individual systems and stellar populations. Using Zahn's tidal theory and MESA main-sequence model grids, we derive the governing parameters $\lambda_{lm}$ and $E_2$, and implement them in the new MINT library of the stellar population code BINARY_C. Our MINT equilibrium tides are 2 to 5 times more efficient than the ubiquitous BSE prescriptions while the radiative-tide efficiency drops sharply with increasing age. We also implement precise initial distributions based on bias-corrected observations. We assess the impact of tides and initial orbital-parameter distributions on circularization and synchronization in eight open clusters, comparing synthetic populations and observations through a bootstrapping method. We find that changing the tidal prescription yields no statistically-significant improvement as both calculations typically lie within 0.5$\sigma$. The initial distribution, especially the primordial concentration of systems at $\log_{10}(P/{\rm d}) \approx 0.8, e\approx 0.05$ dominates the statistics even when artificially increasing tidal strength. This confirms the inefficiency of tides on the main sequence and shows that constraining tidal-efficiency parameters using the $e-\log_{10}(P/{\rm d})$ distribution alone is difficult or impossible. Orbital synchronization carries a more striking age-dependent signature of tidal interactions. In M35 we find twice as many synchronized rotators in our MINT calculation as with BSE. This measure of tidal efficiency is verifiable with combined measurements of orbital parameters and stellar spins., Comment: 24 pages, 29 figures includings appendices. Accepted for publication in MNRAS

We use the rapid binary stellar evolution code BINARY_C to estimate the rate of merging neutron stars with numerous combinations of envelope ejection efficiency and natal kick dispersion. We find a peak in the local rate of merging neutron stars around $\alpha \approx 0.3$$-$$0.4$, depending on the metallicity, where $\alpha$ is the efficiency of utilising orbital energy to unbind the envelope. The peak height decreases with increasing electron-capture supernova kick dispersion $\sigma_\mathrm{ECSN}$. We explain the peak as a competition between the total number of systems that survive the common-envelope phase increasing with $\alpha$ and their separation, which increases with $\alpha$ as well. Increasing $\alpha$ reduces the fraction of systems that merge within a time shorter than the age of the Universe and results in different mass distributions for merging and non-merging double neutron stars. This offers a possible explanation for the discrepancy between the Galactic double neutron star mass distribution and the observed massive merging neutron star event GW190425. Within the $\alpha$$-$$\sigma_\mathrm{ECSN}$ parameter space that we investigate, the rate of merging neutron stars spans several orders of magnitude up to more than $1\times 10^{3} \, \mathrm{Gpc}^{-3}\,\mathrm{yr}^{-1}$ and can be higher than the observed upper limit or lower than the observed lower limit inferred thus far from merging neutron stars detected by gravitational waves. Our results stress the importance of common-envelope physics for the quantitative prediction and interpretation of merging binary neutron star events in this new age of gravitational wave astronomy., Comment: 16 pages, 10 figures, this is an author-produced version of an article accepted for publication in MNRAS following peer review. The version of record is now available online at https://doi.org/10.1093/mnras/stad971 . (N.b. the tables are better formatted in the arXiv article.) Our data are available online at https://zenodo.org/record/7811486#.ZDJt6I5BzJU

Of all the light elements, the evolution of lithium (Li) in the Milky Way is perhaps the most difficult to explain. Li is difficult to synthesize and is easily destroyed, making most stellar sites unsuitable for producing Li in sufficient quantities to account for the proto-solar abundance. For decades, novae have been proposed as a potential explanation to this 'Galactic Li problem', and the recent detection of 7Be in the ejecta of multiple nova eruptions has breathed new life into this theory. In this work, we assess the viability of novae as dominant producers of Li in the Milky Way. We present the most comprehensive treatment of novae in a galactic chemical evolution code to date, testing theoretical- and observationally-derived nova Li yields by integrating metallicity-dependent nova ejecta profiles computed using the binary population synthesis code binary c with the galactic chemical evolution code OMEGA+. We find that our galactic chemical evolution models which use observationally-derived Li yields account for the proto-solar Li abundance very well, while models relying on theoretical nova yields cannot reproduce the proto-solar observation. A brief exploration of physical uncertainties including single-stellar yields, the metallicity resolution of our nova treatment, common-envelope physics, and nova accretion efficiencies indicates that this result is robust to physical assumptions. Scatter within the observationally-derived Li yields in novae is identified as the primary source of uncertainty, motivating further observations of 7Be in nova ejecta., Comment: Accepted for publication in ApJL 28/7/2022

The metallicity of a star affects its evolution in a variety of ways, changing stellar radii, luminosities, lifetimes, and remnant properties. In this work, we use the population synthesis code binary_c to study how metallicity affects novae in the context of binary stellar evolution. We compute a 16-point grid of metallicities ranging from $Z=10^{-4}$ to 0.03, presenting distributions of nova white dwarf masses, accretion rates, delay-times, and initial system properties at the two extremes of our 16-point metallicity grid. We find a clear anti-correlation between metallicity and the number of novae produced, with the number of novae at $Z=0.03$ roughly half that at $Z=10^{-4}$. The white dwarf mass distribution has a strong systematic variation with metallicity, while the shape of the accretion rate distribution is relatively insensitive. We compute a current nova rate of approximately 33 novae per year for the Milky Way, a result consistent with observational estimates relying on extra-Galactic novae but an under-prediction relative to observational estimates relying on Galactic novae. However, the shape of our predicted Galactic white dwarf mass distribution differs significantly to existing observationally derived distributions, likely due to our underlying physical assumptions. In M31, we compute a current nova rate of approximately 36 novae per year, under-predicting the most recent observational estimate of $65^{+15}_{-16}$. Finally, we conclude that when making predictions about currently observable nova rates in spiral galaxies, or stellar environments where star formation has ceased in the distant past, metallicity can likely be considered of secondary importance compared to uncertainties in binary stellar evolution.

Novae are some of the most commonly detected optical transients and have the potential to provide valuable information about binary evolution. Binary population synthesis codes have emerged as the most effective tool for modelling populations of binary systems, but such codes have traditionally employed greatly simplified nova physics, precluding detailed study. In this work, we implement a model treating H and He novae as individual events into the binary population synthesis code \binaryc. This treatment of novae represents a significant improvement on the `averaging' treatment currently employed in modern population synthesis codes. We discuss the evolutionary pathways leading to these phenomena and present nova event rates and distributions of several important physical parameters. Most novae are produced on massive white dwarfs, with approximately 70 and 55 per cent of nova events occurring on O/Ne white dwarfs for H and He novae respectively. Only 15 per cent of H-nova systems undergo a common-envelope phase, but these systems are responsible for the majority of H nova events. All He-accreting He-nova systems are considered post-common-envelope systems, and almost all will merge with their donor star in a gravitational-wave driven inspiral. We estimate the current annual rate of novae in M31 (Andromeda) to be approximately $41 \pm 4$ for H novae, underpredicting the current observational estimate of $65^{+15}_{-16}$, and $0.14\pm0.015$ for He novae. When varying common-envelope parameters, the H nova rate varies between 20 and 80 events per year., Comment: Accepted, MNRAS. 7 Jun 2020: Minor correction regarding AM CVn masses at period bounce, courtesy of P. Neuteufel

Post-asymptotic giant branch (post-AGB) stars with discs are all binaries. Many of these binaries have orbital periods between 100 and 1000 days so cannot have avoided mass transfer between the AGB star and its companion, likely through a common-envelope type interaction. We report on preliminary results of our project to model circumbinary discs around post-AGB stars using our binary population synthesis code binary_c. We combine a simple analytic thin-disc model with binary stellar evolution to estimate the impact of the disc on the binary, and vice versa, fast enough that we can model stellar populations and hence explore the rather uncertain parameter space involved with disc formation. We find that, provided the discs form with sufficient mass and angular momentum, and have an inner edge that is relatively close to the binary, they can both prolong the life of their parent post-AGB star and pump the eccentricity of orbits of their inner binaries., Comment: 8 pages, refereed conference proceedings of the excellent Asymmetric Planetary Nebulae VII conference in Hong Kong, 2018 https://apn7.com/ published in Galaxies

Many aspects of the evolution of stars, and in particular the evolution of binary stars, remain beyond our ability to model them in detail. Instead, we rely on observations to guide our often phenomenological models and pin down uncertain model parameters. To do this statistically requires population synthesis. Populations of stars modelled on computers are compared to populations of stars observed with our best telescopes. The closest match between observations and models provides insight into unknown model parameters and hence the underlying astrophysics. In this brief review, we describe the impact that modern big-data surveys will have on population synthesis, the large parameter space problem that is rife for the application of modern data science algorithms, and some examples of how population synthesis is relevant to modern astrophysics., Comment: Accepted for publication by Cambridge University Press as book chapter in 'The Impact of Binary Stars on Stellar Evolution' by Giacomo Beccari and Henri Boffin. 19 pages, 2 figures

We evolve stellar models to study the common envelope (CE) interaction of an early asymptotic giant branch star of initial mass $5\,\rm M_{\odot}$ with a companion star of mass ranging from $0.1$ to $2\,\rm M_{\odot}$. We model the CE as a fast stripping phase in which the primary experiences rapid mass loss and loses about 80 per cent of its mass. The post-CE remnant is then allowed to thermally readjust during a Roche-lobe overflow (RLOF) phase and the final binary system and its orbital period are investigated. We find that the post-CE RLOF phase is long enough to allow nuclear burning to proceed in the helium shell. By the end of this phase, the donor is stripped of both its hydrogen and helium and ends up as carbon-oxygen white dwarf of mass about $0.8\,\rm M_{\odot}$. We study the sensitivity of our results to initial conditions of different companion masses and orbital separations at which the stripping phase begins. We find that the companion mass affects the final binary separation and that helium-shell burning causes the star to refill its Roche lobe leading to post-CE RLOF. Our results show that double mass transfer in such a binary interaction is able to strip the helium and hydrogen layers from the donor star without the need for any special conditions or fine tuning of the binary parameters., Comment: Accepted to MNRAS 10 pages, 8 figures

The first gravitational wave detections of mergers between black holes and neutron stars represent a remarkable new regime of high-energy transient astrophysics. The signals observed with LIGO-Virgo detectors come from mergers of extreme physical objects which are the end products of stellar evolution in close binary systems. To better understand their origin and merger rates, we have performed binary population syntheses at different metallicities using the new grid-based binary population synthesis code ComBinE. Starting from newborn pairs of stars, we follow their evolution including mass loss, mass transfer and accretion, common envelopes and supernova explosions. We apply the binding energies of common envelopes based on dense grids of detailed stellar structure models, make use of improved investigations of the subsequent Case BB Roche-lobe overflow and scale supernova kicks according to the stripping of the exploding stars. We demonstrate that all the double black hole mergers, GW150914, LVT151012, GW151226, GW170104, GW170608 and GW170814, as well as the double neutron star merger GW170817, are accounted for in our models in the appropriate metallicity regime. Our binary interaction parameters are calibrated to match the accurately determined properties of Galactic double neutron star systems, and we discuss their masses and types of supernova origin. Using our default values for the input physics parameters, we find a double neutron star merger rate of about 3.0 Myr^-1 for Milky-Way equivalent galaxies. Our upper limit to the merger-rate density of double neutron stars is R=400 yr^-1 Gpc^-3 in the local Universe (z=0)., Comment: 36 pages, 26 figures, 8 tables, plus 9 pages appendix. Accepted 2018 August 6 by MNRAS after revision according to referee report (in particular, including further discussions on the progenitor binary of GW170817)

The combination of asteroseismologically-measured masses with abundances from detailed analyses of stellar atmospheres challenges our fundamental knowledge of stars and our ability to model them. Ancient red-giant stars in the Galactic thick disc are proving to be most troublesome in this regard. They are older than 5 Gyr, a lifetime corresponding to an initial stellar mass of about $1.2{\mathrm{M}_{\odot}}$. So why do the masses of a sizeable fraction of thick-disc stars exceed $1.3{\mathrm{M}_{\odot}}$, with some as massive as $2.3{\mathrm{M}_{\odot}}$ ? We answer this question by considering duplicity in the thick-disc stellar population using a binary population-nucleosynthesis model. We examine how mass transfer and merging affect the stellar mass distribution and surface abundances of carbon and nitrogen. We show that a few per cent of thick-disc stars can interact in binary star systems and become more massive than $1.3{\mathrm{M}_{\odot}}$. Of these stars, most are single because they are merged binaries. Some stars more massive than $1.3{\mathrm{M}_{\odot}}$ form in binaries by wind mass transfer. We compare our results to a sample of the APOKASC data set and find reasonable agreement except in the number of these thick-disc stars more massive than $1.3{\mathrm{M}_{\odot}}$. This problem is resolved by the use of a logarithmically-flat orbital-period distribution and a large binary fraction., Comment: This is a pre-copyedited, author-produced version of an article accepted for publication in Monthly Notices of the Royal Astronomical Society following peer review. The version of record [Robert Izzard, 'Binary stars in the Galactic thick disc', Monthly Notices of the Royal Astronomical Society, accepted 7th Sep. 2017] is at https://academic.oup.com/mnras/article-lookup/doi/10.1093/mnras/stx2355

Stellar astrophysicists are increasingly taking into account the effects of orbiting companions on stellar evolution. New discoveries, many thanks to systematic time-domain surveys, have underlined the role of binary star interactions in a range of astrophysical events, including some that were previously interpreted as due uniquely to single stellar evolution. Here, we review classical binary phenomena such as type Ia supernovae, and discuss new phenomena such as intermediate luminosity transients, gravitational wave-producing double black holes, or the interaction between stars and their planets. Finally, we examine the reassessment of well-known phenomena in light of interpretations that include both single and binary stars, for example supernovae of type Ib and Ic or luminous blue variables. At the same time we contextualise the new discoveries within the framework and nomenclature of the corpus of knowledge on binary stellar evolution. The last decade has heralded an era of revival in stellar astrophysics as the complexity of stellar observations is increasingly interpreted with an interplay of single and binary scenarios. The next decade, with the advent of massive projects such as the Large Synoptic Survey Telescope, the Square Kilometre Array, the James Webb Space Telescope and increasingly sophisticated computational methods, will see the birth of an expanded framework of stellar evolution that will have repercussions in many other areas of astrophysics such as galactic evolution and nucleosynthesis., Comment: 40 pages, 18 figures accepted for publication as a Dawes Review of the Publications of the Astronomical Society of Australia

The candidate Thorne-\.Zytkow object (T\.ZO), HV2112, is becoming a well-studied if enigmatic object. A key point of its candidacy as a T\.ZO is whether or not it resides in the Small Magellanic Cloud (SMC). HV2112 has detections in a series of photometric catalogues which have resulted in contradictory estimates of its proper motion and, therefore, its membership within the SMC. This letter seeks to resolve the issue of the SMC membership of HV2112 through a reanalysis of extant photometric data. We also demonstrate the difficulties and downfalls inherent in considering a range of catalogue proper motions. We conclude that the proper motion, and associated ancillary radial velocity, positional and photometric properties, are fully consistent with HV2112 being within the SMC and thus it remains a candidate T\.ZO., Comment: 5 pages, 5 figure, 1 table

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

Classical Cepheids are key probes of both stellar astrophysics and cosmology as standard candles and pulsating variable stars. It is important to understand Cepheids in unprecedented detail in preparation for upcoming GAIA, JWST and extremely-large telescope observations. Cepheid eclipsing binary stars are ideal tools for achieving this goal, however there are currently only three known systems. One of those systems, OGLE-LMC-CEP1812, raises new questions about the evolution of classical Cepheids because of an apparent age discrepancy between the Cepheid and its red giant companion. We show that the Cepheid component is actually the product of a stellar merger of two main sequence stars that has since evolved across the Hertzsprung gap of the HR diagram. This post-merger product appears younger than the companion, hence the apparent age discrepancy is resolved. We discuss this idea and consequences for understanding Cepheid evolution., Comment: 5 pages, 3 figures, accepted to A&A

Many core-collapse supernova progenitors are presumed to be in binary systems. If a star explodes in a binary system, the early supernova light curve can be brightened by the collision of the supernova ejecta with the companion star. The early brightening can be observed when the observer is in the direction of the hole created by the collision. Based on a population synthesis model, we estimate the fractions of core-collapse supernovae in which the light-curve brightening by the collision can be observed. We find that 0.19% of core-collapse supernova light curves can be observed with the collisional brightening. Type Ibc supernova light curves are more likely to be brightened by the collision (0.53%) because of the high fraction of the progenitors being in binary systems and their proximity to the companion stars. Type II and IIb supernova light curves are less affected (~1e-3% and ~1e-2%, respectively). Although the early, slow light-curve declines of some Type IIb and Ibc supernovae are argued to be caused by the collision with the companion star (e.g. SN 2008D), the small expected fraction, as well as the unrealistically small separation required, disfavour the argument. The future transient survey by the Large Synoptic Survey Telescope is expected to detect ~10 Type Ibc supernovae with the early collisional brightening per year, and they will be able to provide information on supernova progenitors in binary systems., Comment: 6 pages, 3 figures, 1 table, accepted by Monthly Notices of the Royal Astronomical Society

Classical Cepheids, like binary stars, are laboratories for stellar evolution and Cepheids in binary systems are especially powerful ones. About one-third of Galactic Cepheids are known to have companions and Cepheids in eclipsing binary systems have recently been discovered in the Large Magellanic Cloud. However, there are no known Galactic binary Cepheids with orbital periods less than one year. We compute population synthesis models of binary Cepheids to compare to the observed period and eccentricity distributions of Galactic Cepheids as well as to the number of observed eclipsing binary Cepheids in the LMC. We find that our population synthesis models are consistent with observed binary properties of Cepheids. Furthermore, we show that binary interaction on the red giant branch prevents some red giant stars from becoming classical Cepheids. Such interactions suggest that the binary fraction of Cepheids should be significantly less than that of their main-sequence progenitors, and that almost all binary Cepheids have orbital periods longer than one year. If the Galactic Cepheid spectroscopic binary fraction is about 35%, then the spectroscopic binary fraction of their intermediate mass main sequence progenitors is about 40-45%., Comment: 7 pages, 3 figures, resubmitted to A&A

Half or more of stars more massive than our Sun are orbited by a companion star in a binary system. Many binaries have short enough orbits that the evolution of both stars is greatly altered by an exchange of mass and angular momentum between the stars. Such mass transfer is highly likely on the asymptotic giant branch (AGB) because this is when a star is both very large and has strong wind mass loss. Direct mass transfer truncates the AGB, and its associated nucleosynthesis, prematurely compared to in a single star. In wide binaries we can probe nucleosynthesis in the long-dead AGB primary star by today observing its initially lower-mass companion. The star we see now may be polluted by ejecta from the primary either through a wind or Roche-lobe overflow. We highlight recent quantitative work on nucleosynthesis in (ex-)AGB mass-transfer systems, such as carbon and barium stars, the link between binary stars and planetary nebulae, and suggest AGB stars as a possible source of the enigmatic element, lithium., Comment: To appear in the proceedings of "Why Galaxies Care About AGB Stars III" http://www.univie.ac.at/galagb/

We use our new population synthesis code BONNFIRES to test how surface abundances predicted by rotating stellar models depend on the numerical treatment of rotational mixing, such as spatial resolution, temporal resolution and computation of mean molecular weight gradients. We find that even with identical numerical prescriptions for calculating the rotational mixing coefficients in the diffusion equation, different timesteps lead to a deviation of the coefficients and hence surface abundances. We find the surface abundances vary by 10-100% between the model sequences with short timestep of 0.001Myr to model sequences with longer timesteps. Model sequences with stronger surface nitrogen enrichment also have longer main-sequence lifetimes because more hydrogen is mixed to the burning cores. The deviations in main-sequence lifetimes can be as large as 20%. Mathematically speaking, no numerical scheme can give a perfect solution unless infinitesimally small timesteps are used. However, we find that the surface abundances eventually converge within 10% between modelling sequences with sufficiently small timesteps below 0.1Myr. The efficiency of rotational mixing depends on the implemented numerical scheme and critically on the computation of the mean molecular weight gradient. A smoothing function for the mean molecular weight gradient results in stronger rotational mixing. If the discretization scheme or the computational recipe for calculating the mean molecular weight gradient is altered, re-calibration of mixing parameters may be required to fit observations. If we are to properly understand the fundamental physics of rotation in stars, it is crucial that we minimize the uncertainty introduced into stellar evolution models when numerically approximating rotational mixing processes., Comment: 8 pages, 6 figures, accepted by A&A

Powerful telescopes equipped with multi-fibre or integral field spectrographs combined with detailed models of stellar atmospheres and automated fitting techniques allow for the analysis of large number of stars. These datasets contain a wealth of information that require new analysis techniques to bridge the gap between observations and stellar evolution models. To that end, we develop BONNSAI (BONN Stellar Astrophysics Interface), a Bayesian statistical method, that is capable of comparing all available observables simultaneously to stellar models while taking observed uncertainties and prior knowledge such as initial mass functions and distributions of stellar rotational velocities into account. BONNSAI can be used to (1) determine probability distributions of fundamental stellar parameters such as initial masses and stellar ages from complex datasets, (2) predict stellar parameters that were not yet observationally determined and (3) test stellar models to further advance our understanding of stellar evolution. An important aspect of BONNSAI is that it singles out stars that cannot be reproduced by stellar models through $\chi^{2}$ hypothesis tests and posterior predictive checks. BONNSAI can be used with any set of stellar models and currently supports massive main-sequence single star models of Milky Way and Large and Small Magellanic Cloud composition. We apply our new method to mock stars to demonstrate its functionality and capabilities. In a first application, we use BONNSAI to test the stellar models of Brott et al. (2011a) by comparing the stellar ages inferred for the primary and secondary stars of eclipsing Milky Way binaries. Ages are determined from dynamical masses and radii that are known to better than 3%. We find that the stellar models reproduce the Milky Way binaries well. BONNSAI is available through a web-interface at http://www.astro.uni-bonn.de/stars/bonnsai., Comment: Accepted for publication in A&A; 15 pages, 10 figures, 4 tables; BONNSAI is available through a web-interface at http://www.astro.uni-bonn.de/stars/bonnsai

Massive stars rapidly change their masses through strong stellar winds and mass transfer in binary systems. We show that such mass changes leave characteristic signatures in stellar mass functions of young star clusters which can be used to infer their ages and to identify products of binary evolution. We model the observed present day mass functions of the young Galactic Arches and Quintuplet star clusters using our rapid binary evolution code. We find that shaping of the mass function by stellar wind mass loss allows us to determine the cluster ages to 3.5$\pm$0.7 Myr and 4.8$\pm$1.1 Myr, respectively. Exploiting the effects of binary mass exchange on the cluster mass function, we find that the most massive stars in both clusters are rejuvenated products of binary mass transfer, i.e. the massive counterpart of classical blue straggler stars. This resolves the problem of an apparent age spread among the most luminous stars exceeding the expected duration of star formation in these clusters. We perform Monte Carlo simulations to probe stochastic sampling, which support the idea of the most massive stars being rejuvenated binary products. We find that the most massive star is expected to be a binary product after 1.0$\pm$0.7 Myr in Arches and after 1.7$\pm$1.0 Myr in Quintuplet. Today, the most massive 9$\pm$3 stars in Arches and 8$\pm$3 in Quintuplet are expected to be such objects. Our findings have strong implications for the stellar upper mass limit and solve the discrepancy between the claimed 150 $\mathrm{M}_\odot$ limit and observations of fours stars with initial masses of 165-320 $\mathrm{M}_\odot$ in R136 and of SN 2007bi, which is thought to be a pair-instability supernova from an initial 250 $\mathrm{M}_\odot$ star. Using the stellar population of R136, we revise the upper mass limit to values in the range 200-500 $\mathrm{M}_\odot$., Comment: 13 pages (+4 pages appendix), 8 figures, 2 tables; accepted for publication in ApJ

BONNFIRES, a new generation of population synthesis code, can calculate nuclear reaction, various mixing processes and binary interaction in a timely fashion. We use this new population synthesis code to study the interplay between binary mass transfer and rotation. We aim to compare theoretical models with observations, in particular the surface nitrogen abundance and rotational velocity. Preliminary results show binary interactions may explain the formation of nitrogen-rich slow rotators and nitrogen-poor fast rotators, but more work needs to be done to estimate whether the observed frequencies of those stars can be matched., Comment: 4 pages, 2 figures, conference proceeding of "Setting a new standard in the analysis of binary stars"

CEMP-s stars are very metal-poor stars with enhanced abundances of carbon and s-process elements. They form a significant proportion of the very metal-poor stars in the Galactic halo and are mostly observed in binary systems. This suggests that the observed chemical anomalies are due to mass accretion in the past from an asymptotic giant branch (AGB) star. Because CEMP-s stars have hardly evolved since their formation, the study of their observed abundances provides a way to probe our models of AGB nucleosynthesis at low metallicity. To this end we included in our binary evolution model the results of the latest models of AGB nucleosynthesis and we simulated a grid of 100,000 binary stars at metallicity Z=0.0001 in a wide range of initial masses and separations. We compared our modelled stars with a sample of 60 CEMP-s stars from the SAGA database of metal-poor stars. For each observed CEMP-s star of the sample we found the modelled star that reproduces best the observed abundances. The result of this comparison is that we are able to reproduce simultaneously the observed abundance of the elements affected by AGB nucleosynthesis (e.g. C, Mg, s-elements) for about 60% of the stars in the sample., Comment: 4 pages, proceedings of the Conference "Setting a new standard in the analysis of binary stars" held 16-19 September 2013 in Leuven. The original publication will be available at http://www.edpsciences.org/eas

It is likely that at least some planetary nebulae are composed of matter which was ejected from a binary star system during common-envelope (CE) evolution. For these planetary nebulae the ionizing component is the hot and luminous remnant of a giant which had its envelope ejected by a companion in the process of spiralling-in to its current short-period orbit. A large fraction of CE phases which end with ejection of the envelope are thought to be initiated by low-mass red giants, giants with inert, degenerate helium cores. We discuss the possible end-of-CE structures of such stars and their subsequent evolution to investigate for which structures planetary nebulae are formed. We assume that a planetary nebula forms if the remnant reaches an effective temperature greater than 30 kK within 10^4 yr of ejecting its envelope. We assume that the composition profile is unchanged during the CE phase so that possible remnant structures are parametrized by the end-of-CE core mass, envelope mass and entropy profile. We find that planetary nebulae are expected in post-CE systems with core masses greater than about 0.3 solar masses if remnants end the CE phase in thermal equilibrium. We show that whether the remnant undergoes a pre-white dwarf plateau phase depends on the prescribed end-of-CE envelope mass. Thus, observing a young post-CE system would constrain the end-of CE envelope mass and post-CE evolution., Comment: Published in MNRAS. 12 pages, 12 figures. Minor changes to match published version

Globular clusters contain many stars with surface abundance patterns indicating contributions from hydrogen burning products, as seen in the anti-correlated elemental abundances of e.g. sodium and oxygen, and magnesium and aluminium. Multiple generations of stars can explain this phenomenon, with the second generation forming from a mixture of pristine gas and ejecta from the first generation. We show that massive binary stars may be a source of much of the material that makes this second generation of stars. Mass transfer in binaries is often non-conservative and the ejected matter moves slowly enough that it can remain inside a globular cluster and remain available for subsequent star formation. Recent studies show that there are more short-period massive binaries than previously thought, hence also more stars that interact and eject nuclear-processed material., Comment: Contribution to the proceedings of "Reading the book of globular clusters with the lens of stellar evolution", Rome, 26-28 November 2012

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

Circumbinary discs are commonly observed around post-asymptotic giant branch (AGB) systems and are known to play an important role in their evolution. Several studies have pointed out that a circumbinary disc interacts through resonances with the central binary and leads to angular momentum transfer from the central binary orbit to the disc. This interaction may be responsible for a substantial increase in the binary eccentricity. We investigate whether this disc eccentricity-pumping mechanism can be responsible for the high eccentricities commonly found in post-AGB binary systems., Comment: 7 pages, 3 figures

The formation mechanism of the barium stars is thought to be well understood. Barium-rich material, lost in a stellar wind from a thermally-pulsing asymptotic-giant branch star in a binary system, is accreted by its companion main-sequence star. Now, many millions of years later, the primary is an unseen white dwarf and the secondary has itself evolved into a giant which displays absorption lines of barium in its spectrum and is what we call a barium star. A similar wind-accretion mechanism is also thought to form the low-metallicity CH and carbon-enhanced metal-poor stars. Qualitatively the picture seems clear but quantitatively it is decidedly murky: several key outstanding problems remain which challenge our basic understanding of binary-star physics. Barium stars with orbital periods less than about 4,000 days should -- according to theory -- be in circular orbits because of tidal dissipation, yet they are often observed to be eccentric. Only one barium-star period longer than 10^4 days has been published although such stars are predicted to exist in large numbers. In this paper we attempt to shed light on these problems. First, we consider the impact of kicking the white dwarf at its birth, a notion which is supported by independent evidence from studies of globular clusters. Second, we increase the amount of orbital angular momentum loss during wind mass transfer, which shrinks barium-star binaries to the required period range. We conclude with a discussion of possible physical mechanisms and implications of a kick, such as the break up of wide barium-star binaries and the limits imposed on our models by observations., Comment: 7 pages, 6 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

We compare our latest single and binary stellar model results from the Cambridge STARS code to several sets of observations. We examine four stellar population ratios, the number of blue to red supergiants, the number of Wolf-Rayet stars to O supergiants, the number of red supergiants to Wolf-Rayet stars and the relative number of Wolf-Rayet subtypes, WC to WN stars. These four ratios provide a quantitative measure of nuclear burning lifetimes and the importance of mass loss during various stages of the stars' lifetimes. In addition we compare our models to the relative rate of type Ib/c to type II supernovae to measure the amount of mass lost over the entire lives of all stars. We find reasonable agreement between the observationally inferred values and our predicted values by mixing single and binary star populations. However there is evidence that extra mass loss is required to improve the agreement further, to reduce the number of red supergiants and increase the number of Wolf-Rayet stars., Comment: Accepted for publication in MNRAS, 11 pages, 10 figures

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

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 present Window To The Stars, a graphical user interface to the popular TWIN single/binary stellar evolution code, for novices, students and professional astrophysicists. It removes the drudgery associated with the traditional approach to running the code, while maintaining the power, output quality and flexibility a modern stellar evolutionist requires. It is currently being used for cutting edge research and interactive teaching., Comment: Accepted for publication in New Astronomy

We describe our first attempt at modelling nucleosynthesis in massive AGB stars which have undergone core carbon burning, the super-AGB stars. We fit a synthetic model to detailed stellar evolution models in the mass range 9<=M/Msun<=11.5 (Z=0.02), and extrapolate these fits to the end of the AGB. We determine the number of thermal pulses and AGB lifetime as a function of mass and mass-loss prescription. Our preliminary nucleosynthesis calculations show that, for a reasonable mass-loss rate, the effect of hot-bottom burning in super-AGB stars on the integrated yield of a stellar population is not large. There are many uncertainties, such as mass-loss and convective overshooting, which prevent accurate yield calculations. However, as potential progenitors of electron-capture supernovae, these stars may contribute 7% of non-type-Ia supernovae., Comment: To appear in the proceedings of the Eighth Torino Workshop on Nucleosynthesis in AGB Stars (6 pages, 3 figures)

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

We study the evolution of stars that may be the progenitors of long-soft gamma-ray burst (GRBs) -- rotating naked helium stars presumed to have lost their envelopes to winds or companions. Our aim is to investigate the formation and development of single and binary systems and from this population evaluate the rates of interesting individual species. Using a rapid binary evolution algorithm, that enables us to model the most complex binary systems and to explore the effect of metallicity on GRB production, we draw the following conclusions. First we find that, if we include an approximate treatment of angular momentum transport by mass loss, the resulting spin rates for single stars become too low to form a centrifugally supported disc that can drive a GRB engine. Second massive stars in binaries result in enough angular momentum -- due to spin-orbit tidal interactions -- to form a centrifugally supported disc and are thus capable of supplying a sufficient number of progenitors. This holds true even if only a small fraction of bursts are visible to a given observer and the GRB rate is several hundred times larger than the observed rate. Third low-metallicity stars aid the formation of a rapidly rotating, massive helium cores at collapse and so their evolution is likely to be affected by the local properties of the ISM. This effect could increase the GRB formation rate by a factor of 5-7 at Z=Z_solar/200. Finally we quantify the effects of mass loss, common-envelope evolution and black-hole formation and show that more stringent constraints to many of these evolution parameters are needed in order to draw quantitative conclusions from population synthesis work., Comment: 37 pages, 14 figures, MNRAS in press