Metallicity is known to significantly affect the radial expansion of a massive star: the lower the metallicity, the more compact the star, especially during its post-MS evolution. We study this effect in the context of binary evolution. Using the stellar-evolution code MESA, we computed evolutionary tracks of stars at different metallicities, exploring variations of factors known to affect the radial expansion (e.g. semiconvection, overshooting, rotation). We find observational support for an evolution in which already at metallicity $0.2Z_{\odot}$ massive stars remain relatively compact during the Hertzprung-Gap (HG) phase and most of their expansion occurs during core-helium burning (CHeB). Consequently, we show that metallicity has a strong influence on the type of mass transfer evolution in binary systems. At solar metallicity, a case-B mass transfer is initiated shortly after the end of MS, and a giant donor is almost always a rapidly expanding HG star. At lower metallicity, the parameter space for mass transfer from a more evolved CHeB star increases dramatically. This means that envelope stripping and formation of helium stars in low-metallicity environments occurs later in the evolution of the donor, implying a much shorter duration of the Wolf-Rayet phase (even by an order of magnitude) and higher final core masses. This metallicity effect is independent of the impact of metallicity-dependent stellar winds. At very low metallicities, a significant fraction of massive stars in binaries engages in the first episode of mass transfer very late into their evolution, when they already have a well-developed CO core. The remaining lifetime ($< 10^4$ yr) is unlikely to be enough to strip the entire H-rich envelope. We also briefly discuss the extremely small parameter space for mass transfer from massive convective-envelope donors in the context of binary black hole merger formation.

Aluminium-26 is a short-lived radionuclide with a half-life of 0.72Myr, which is observed today in the Galaxy via gamma-ray spectroscopy and is inferred to have been present in the early Solar System via analysis of meteorites. Massive stars are considered the main contributors of Al26. Although most massive stars are found in binary systems, the effect, however, of binary interactions on the Al26 yields have not been investigated since Braun & Langer (1995). Here we aim to fill this gap. We have used the MESA stellar evolution code to compute massive (10Msun<=M<=80Msun), non-rotating, single and binary stars of solar metallicity (Z=0.014). We computed the wind yields for the single stars and for the binary systems where mass transfer plays a major role. Depending on the initial mass of the primary star and orbital period, the Al26 yield can either increase or decrease in a binary system. For binary systems with primary masses up to ~35-40Msun, the yield can increase significantly, especially at the lower mass-end, while above ~45Msun the yield becomes similar to the single star yield or even decreases. Our preliminary results show that compared to supernova explosions, the contribution of mass-loss in binary systems to the total Al26 abundance produced by a stellar population is minor. On the other hand, if massive star mass-loss is the origin of Al26 in the early Solar System, our results will have significant implications for the identification of the potential stellar, or stellar population, source.

Many disc-type post-asymptotic giant branch (post-AGB) stars are chemically peculiar, showing underabundances of refractory elements in their photospheres that correlate with condensation temperature. The aim of this paper is to investigate how accretion from a circumbinary disc can cause this phenomenon of depletion and how this impacts the evolution of post-AGB stars. We used the \texttt{MESA} code to evolve stars in the post-AGB phase, while including accretion of metal-poor gas. We compared the models to a sample of 58 observed disc-type post-AGB stars with chemical abundance data. For each of these stars, we estimated the luminosity and the mass using the Gaia distance. We modelled the accretion rate onto the binary from a viscously evolving disc for a range of initial accretion rates and disc masses. We find that large initial accretion rates ($\gtrsim 3\times10^{-7}$ $M_\odot$/yr) and large initial disc masses ($\sim10^{-2}$ $M_\odot$) are needed to reproduce the observed depleted post-AGB stars. Based on these high accretion rates, the evolution timescale of post-AGB stars can be significantly extended by a factor between two and five. We distinguish depletion patterns that are unsaturated (plateau profile) from those that are saturated, and we expect that post-red giant branch (post-RGB) stars are much more likely to show an unsaturated abundance pattern compared to post-AGB stars. Finally, because of the slower evolution of the low-mass post-RGB stars, we find that these systems can become depleted at lower effective temperatures ($< 5000$ K). We conclude that accretion from a circumbinary disc successfully accounts for the chemical peculiarity of post-AGB stars.

Carbon-enhanced metal-poor stars, CH stars and barium stars, among other classes of chemically peculiar stars, are thought to be products of the interaction of low- and intermediate-mass binaries which occurred when the most evolved star was in the asymptotic giant branch (AGB) phase. Binary evolution models predict that if the initial orbital periods of such systems are shorter than a few thousand days, their orbits should have circularised due to tidal effects. However, observations of the progeny of AGB binaries show that many of these objects have substantial eccentricities, up to about 0.9. In this work we explore the impact of wind mass transfer on the orbits of AGB binaries by performing numerical simulations in which the AGB wind is modelled using a hydrodynamical code and the stellar dynamics is evolved using an N-body code. We find that in most models wind mass transfer contributes to the circularisation of the orbit, but on longer timescales than tidal circularisation if the eccentricity is less than about 0.4. For low initial wind velocities and pseudo-synchronisation of the donor star, we find a structure resembling wind Roche-lobe overflow near periastron. In this case, the interaction between the gas and the stars is stronger than for high initial wind velocities and the orbit shrinks while the eccentricity decreases. In one of our models wind interaction is found to pump the eccentricity on a similar timescale as tidal circularisation. Although our study is based on a small sample of models, it offers some insight into the orbital evolution of eccentric binaries interacting via winds. A larger grid of numerical models for different binary parameters is needed to test if a regime exists where hydrodynamical eccentricity pumping can effectively counteract tidal circularisation, and if this can explain the puzzling eccentricities of the descendants of AGB binaries.

We present stellar_wind.py, a module that provides multiple methods of simulating stellar winds using smoothed particle hydrodynamics codes (SPH) within the astrophysical multipurpose software environment (AMUSE) framework. With the simple wind mode, we create SPH wind particles in a spherically symmetric shell. We inject the wind particles with a velocity equal to their terminal velocity. The accelerating wind mode is similar, but with this method particles can be injected with a lower initial velocity than the terminal velocity and they are accelerated away from the star according to an acceleration function. With the heating wind mode, SPH particles are created with zero initial velocity with respect to the star, but instead wind particles are given an internal energy based on the integrated mechanical luminosity of the star. This mode is designed to be used on longer timescales and larger spatial scales compared to the other two modes and assumes that the star is embedded in a gas cloud. For fast winds, we find that both the simple and accelerating mode can reproduce the desired velocity, density and temperature profiles. For slow winds, the simple wind mode is insufficient due to dominant hydrodynamical effects that change the wind velocities. The accelerating mode, with additional options to account for these hydrodynamical effects, can still reproduce the desired wind profiles. We test the heating mode by simulating both a normal wind and a supernova explosion of a single star in a uniform density medium. The stellar wind simulation results matches the analytical solution for an expanding wind bubble. The supernova simulation gives qualitatively correct results, but the simulated bubble expands faster than the analytical solution predicts. We conclude with an example of a triple star system which includes the colliding winds of all three stars.

Wind mass transfer in binary systems with AGB donor stars plays a fundamental role in the formation of a variety of objects, including barium stars and CEMP stars. We carry out a comprehensive set of SPH simulations of wind-losing AGB stars in binaries, for a variety of binary mass ratios, orbital separations, initial wind velocities and rotation rates of the donor star. The initial parameters of the simulated systems are chosen to match the expected progenitors of CEMP stars. We find that the strength of interaction between the wind and the stars depends on both the wind-velocity-to-orbital-velocity ratio ($v_\infty/v_\mathrm{orb}$) and the binary mass ratio. Strong interaction occurs for close systems and comparable mass ratios, and gives rise to a complex morphology of the outflow and substantial angular-momentum loss, which leads to a shrinking of the orbit. As the orbital separation increases and the mass of the companion star decreases, the morphology of the outflow, as well as the angular-momentum loss, become more similar to the spherically symmetric wind case. We also explore the effects of tidal interaction and find that for orbital separations up to 7-10 AU, depending on mass ratio, spin-orbit coupling of the donor star occurs at some point during the AGB phase. If the initial wind velocity is relatively low, we find that corotation of the donor star results in a modified outflow morphology that resembles wind Roche-lobe overflow. In this case the mass-accretion efficiency and angular-momentum loss differ from those found for a non-rotating donor. Finally, we provide a relation for both the mass-accretion efficiency and angular-momentum loss as a function of $v_\infty/v_\mathrm{orb}$ and the binary mass ratio that can be easily implemented in a population synthesis code to study populations of barium stars, CEMP stars and other products of interaction in AGB binaries.

Neutron star mergers (NSM) are likely to be the main production sites for the rapid (r-) neutron capture process elements. We study the r-process enrichment of the stellar halo of the Milky Way through NSM, by tracing the typical r-process element Eu in the Munich-Groningen semi-analytic galaxy formation model, applied to three high resolution Aquarius dark matter simulations. In particular, we investigate the effect of the kick velocities that neutron star binaries receive upon their formation, in the building block galaxies (BBs) that partly formed the stellar halo by merging with our Galaxy. When this kick is large enough to overcome the escape velocity of the BB, the NSM takes place outside the BB with the consequence that there is no r-process enrichment. We find that a standard distribution of NS kick velocities decreases [Eu/Mg] abundances of halo stars by $\sim 0.5$~dex compared to models where NS do not receive a kick. With low NS kick velocities, our simulations match observed [Eu/Mg] abundances of halo stars reasonably well, for stars with metallicities [Mg/H]$\geq -1.5$. Only in Aquarius halo B-2 also the lower metallicity stars have [Eu/Mg] values similar to observations. We conclude that our assumption of instantaneous mixing is most likely inaccurate for modelling the r-process enrichment of the Galactic halo, or an additional production site for r-process elements is necessary to explain the presence of low-metallicity halo stars with high Eu abundances.

Binary post-asymptotic giant branch (post-AGB) stars are thought to be the products of a strong but poorly-understood interaction during the AGB phase. The aim of this contribution is to update the orbital elements of a sample of galactic post-AGB binaries observed in a long-term radial-velocity monitoring campaign. Radial velocities are computed from high signal-to-noise spectra by use of a cross-correlation method. The radial-velocity curves are fitted by using both a least-squares algorithm and a Nelder-Mead simplex algorithm. We use a Monte Carlo method to compute uncertainties on the orbital elements. The resulting mass functions are used to derive a companion mass distribution by optimising the predicted to the observed cumulative mass-function distributions, after correcting for observational bias. As a result, we derive and update orbital elements for 33 galactic post-AGB binaries, among which 3 are new orbits. The orbital periods of the systems range from 100 to about 3000 days. Over 70 percent (23 out of 33) of our binaries have significant non-zero eccentricities ranging over all periods. Their orbits are non-circular despite the fact that the Roche-lobe radii are smaller than the maximum size of a typical AGB star and tidal circularisation should have been strong when the objects were on the AGB. We derive a distribution of companion masses that is peaked around 1.09 $M_\odot$ with a standard deviation of 0.62 $M_\odot$. The large spread in companion masses highlights the diversity of post-AGB binary systems. Furthermore, we find that only post-AGB stars with high effective temperatures (> 5500 K) in wide orbits are depleted in refractory elements, suggesting that re-accretion of material from a circumbinary disc is an ongoing process. It appears, however, that chemical depletion is inefficient for the closest orbits irrespective of the actual surface temperature.

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

In low-mass binary systems, mass transfer is likely to occur via a slow and dense stellar wind when one of the stars is in the AGB phase. Observations show that many binaries that have undergone AGB mass transfer have orbital periods of 1-10 yr, at odds with the predictions of binary population synthesis models. We investigate the mass-accretion efficiency and angular-momentum loss via wind mass transfer in AGB binary systems. We use these quantities to predict the evolution of the orbit. We perform 3D hydrodynamical simulations of the stellar wind lost by an AGB star using the AMUSE framework. We approximate the thermal evolution of the gas by imposing a simple effective cooling balance and we vary the orbital separation and the velocity of the stellar wind. We find that for wind velocities $v_{\infty}$ larger than the relative orbital velocity of the system $v_\mathrm{orb}$ the flow is described by the Bondi-Hoyle-Lyttleton approximation and the angular-momentum loss is modest, leading to an expansion of the orbit. For low wind velocities an accretion disk is formed around the companion and the accretion efficiency as well as the angular-momentum loss are enhanced, implying that the orbit will shrink. We find that the transfer of angular momentum from the orbit to the outflowing gas occurs within a few orbital separations from the center of mass of the binary. Our results suggest that the orbital evolution of AGB binaries can be predicted as a function of the ratio $v_{\infty}/v_\mathrm{orb}$. Our results can provide insight into the puzzling orbital periods of post-AGB binaries and suggest that the number of stars entering into the common-envelope phase will increase. The latter can have significant implications for the expected formation rates of the end products of low-mass binary evolution, such as cataclysmic binaries, type Ia supernova and double white-dwarf mergers. [ABRIDGED]

The chemical yields of supernovae and the metal enrichment of the hot intra-cluster medium (ICM) are not well understood. This paper introduces the CHEmical Enrichment RGS Sample (CHEERS), which is a sample of 44 bright local giant ellipticals, groups and clusters of galaxies observed with XMM-Newton. This paper focuses on the abundance measurements of O and Fe using the reflection grating spectrometer (RGS). The deep exposures and the size of the sample allow us to quantify the intrinsic scatter and the systematic uncertainties in the abundances using spectral modeling techniques. We report the oxygen and iron abundances as measured with RGS in the core regions of all objects in the sample. We do not find a significant trend of O/Fe as a function of cluster temperature, but we do find an intrinsic scatter in the O and Fe abundances from cluster to cluster. The level of systematic uncertainties in the O/Fe ratio is estimated to be around 20-30%, while the systematic uncertainties in the absolute O and Fe abundances can be as high as 50% in extreme cases. We were able to identify and correct a systematic bias in the oxygen abundance determination, which was due to an inaccuracy in the spectral model. The lack of dependence of O/Fe on temperature suggests that the enrichment of the ICM does not depend on cluster mass and that most of the enrichment likely took place before the ICM was formed. We find that the observed scatter in the O/Fe ratio is due to a combination of intrinsic scatter in the source and systematic uncertainties in the spectral fitting, which we are unable to disentangle. The astrophysical source of intrinsic scatter could be due to differences in AGN activity and ongoing star formation in the BCG. The systematic scatter is due to uncertainties in the spatial line broadening, absorption column, multi-temperature structure and the thermal plasma models. (Abbreviated).

Misalignments between the orbital planes of planets and the equatorial planes of their host stars have been observed in our solar system, in transiting exoplanets, and in the orbital planes of debris discs. We present a mechanism that causes such a spin-orbit misalignment for a protoplanetary disc due to its movement through an ambient medium. Our physical explanation of the mechanism is based on the theoretical solutions to the Stark problem. We test this idea by performing self-consistent hydrodynamical simulations and simplified gravitational $N$-body simulations. The $N$-body model reduces the mechanism to the relevant physical processes. The hydrodynamical simulations show the mechanism in its full extent, including gas-dynamical and viscous processes in the disc which are not included in the theoretical framework. We find that a protoplanetary disc embedded in a flow changes its orientation as its angular momentum vector tends to align parallel to the relative velocity vector. Due to the force exerted by the flow, orbits in the disc become eccentric, which produces a net torque and consequentially changes the orbital inclination. The tilting of the disc causes it to contract. Apart from becoming lopsided, the gaseous disc also forms a spiral arm even if the inclination does not change substantially. The process is most effective at high velocities and observational signatures are therefore mostly expected in massive star-forming regions and around winds or supernova ejecta. Our $N$-body model indicates that the interaction with supernova ejecta is a viable explanation for the observed spin-orbit misalignment in our solar system.

Observations indicate that the dispersal of protoplanetary discs in star clusters occurs on time scales of about 5 Myr. Several processes are thought to be responsible for this disc dispersal. Here we compare two of these processes: dynamical encounters and interaction with the interstellar medium, which includes face-on accretion and ram pressure stripping. We perform simulations of embedded star clusters with parameterisations for both processes to determine the environment in which either of these processes is dominant. We find that face-on accretion, including ram pressure stripping, is the dominant disc truncation process if the fraction of the total cluster mass in stars is $\lesssim 30\,\%$ regardless of the cluster mass and radius. Dynamical encounters require stellar densities $\gtrsim 10^4$ pc$^{-3}$ combined with a mass fraction in stars of $\approx 90\,\%$ to become the dominant process. Our results show that during the embedded phase of the cluster, the truncation of the discs is dominated by face-on accretion and dynamical encounters become dominant when the intra-cluster gas has been expelled. As a result of face-on accretion the protoplanetary discs become compact and their surface density increases. In contrast, dynamical encounters lead to discs that are less massive and remain larger.

Observations indicate that stars generally lose their protoplanetary discs on a timescale of about 5 Myr. Which mechanisms are responsible for the disc dissipation is still debated. Here we investigate the movement through an ambient medium as a possible cause of disc dispersal. The ram pressure exerted by the flow can truncate the disc and the accretion of material with no azimuthal angular momentum leads to further disc contraction. We derive a theoretical model from accretion disc theory that describes the evolution of the disc radius, mass, and surface density profile as a function of the density and velocity of the ambient medium. We test our model by performing hydrodynamical simulations of a protoplanetary disc embedded in a flow with different velocities and densities. We find that our model gives an adequate description of the evolution of the disc radius and accretion rate onto the disc. The total disc mass in the simulations follows the theoretically expected trend, except at the lowest density where our simulated discs lose mass owing to continuous stripping. This stripping may be a numerical rather than a physical effect. Some quantitative differences exist between the model predictions and the simulations. These are at least partly caused by numerical viscous effects in the disc and depend on the resolution of the simulation. Our model can be used as a conservative estimate for the process of face-on accretion onto protoplanetary discs, as long as viscous processes in the disc can be neglected. The model predicts that in dense gaseous environments, discs can shrink substantially in size and can, in theory, sweep up an amount of gas of the order of their initial mass. This process could be relevant for planet formation in dense environments.

The hot intra-cluster medium (ICM) is rich in metals, which are synthesised by supernovae (SNe) and accumulate over time into the deep gravitational potential well of clusters of galaxies. Since most of the elements visible in X-rays are formed by type Ia (SNIa) and/or core-collapse (SNcc) supernovae, measuring their abundances gives us direct information on the nucleosynthesis products of billions of SNe since the epoch of the star formation peak (z~2-3). In this study, we compare the most accurate average X/Fe abundance ratios (compiled in a previous work from XMM-Newton EPIC and RGS observations of 44 galaxy clusters, groups, and ellipticals), representative of the chemical enrichment in the nearby ICM, to various SNIa and SNcc nucleosynthesis models found in the literature. The use of a SNcc model combined to any favoured standard SNIa model (deflagration or delayed-detonation) fails to reproduce our abundance pattern. In particular, the Ca/Fe and Ni/Fe ratios are significantly underestimated by the models. We show that the Ca/Fe ratio can be reproduced better, either by taking a SNIa delayed-detonation model that matches the observations of the Tycho supernova remnant, or by adding a contribution from the Ca-rich gap transient SNe, whose material should easily mix into the hot ICM. On the other hand, the Ni/Fe ratio can be reproduced better by assuming that both deflagration and delayed-detonation SNIa contribute in similar proportions to the ICM enrichment. In either case, the fraction of SNIa over the total number of SNe (SNIa+SNcc) contributing to the ICM enrichment ranges within 29-45%. This fraction is found to be systematically higher than the corresponding SNIa/SNe fraction contributing to the enrichment of the proto-solar environnement (15-25%). We also discuss and quantify two useful constraints on both SNIa and SNcc that can be inferred from the ICM abundance ratios.

Globular clusters (GCs) are known to harbor multiple stellar populations. To explain these observations Bastian et al. suggested a scenario in which a second population is formed by the accretion of enriched material onto the low-mass stars in the initial GC population. The idea is that the low-mass, pre-main sequence stars sweep up gas expelled by the massive stars of the same generation into their protoplanetary disc as they move through the GC core. We perform simulations with 2 different smoothed particle hydrodynamics codes to investigate if a low-mass star surrounded by a protoplanetary disc can accrete the amount of enriched material required in this scenario. We focus on the gas loading rate onto the disc and star as well as on the lifetime of the disc. We find that the gas loading rate is a factor of 2 smaller than the geometric rate, because the effective cross section of the disc is smaller than its surface area. The loading rate is consistent for both codes, irrespective of resolution. The disc gains mass in the high resolution runs, but loses angular momentum on a time scale of 10^4 yrs. Two effects determine the loss of (specific) angular momentum in our simulations: 1) continuous ram pressure stripping and 2) accretion of material with no azimuthal angular momentum. Our study and previous work suggest that the former, dominant process is mainly caused by numerical rather than physical effects, while the latter is not. The latter process causes the disc to become more compact, increasing the surface density profile at smaller radii. The disc size is determined in the first place by the ram pressure when the flow first hits the disc. Further evolution is governed by the decrease in the specific angular momentum of the disc. We conclude that the size and lifetime of the disc are probably not sufficient to accrete the amount of mass required in Bastian et al.'s scenario.

The X-ray source HLX-1 near the spiral galaxy ESO 243-49 is currently the best intermediate-mass black hole candidate. It has a peak bolometric luminosity of $10^{42}$ erg s$^{-1}$, which implies a mass inflow rate of $\sim10^{-4}$ MSun yr$^{-1}$, but the origin of this mass is unknown. It has been proposed that there is a star on an eccentric orbit around the black hole which transfers mass at pericentre. To investigate the orbital evolution of this system, we perform stellar evolution simulations using mesa and SPH simulations of a stellar orbit around an intermediate-mass black hole using fi. We run and couple these simulations using the amuse framework. We find that mass is lost through both the first and second Lagrange points and that there is a delay of up to 10 days between the pericentre passage and the peak mass loss event. The orbital evolution timescales we find in our simulations are larger than what is predicted by analytical models, but these models fall within the errors of our results. Despite the fast orbital evolution, we are unable to reproduce the observed change in outburst period. We conclude that the change in the stellar orbit with the system parameters investigated here is unable to account for all observed features of HLX-1.

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.

Many observed CEMP stars are found in binary systems and show enhanced abundances of $s$-elements. The origin of the chemical abundances of these CEMP-$s$ stars is believed to be accretion in the past of enriched material from a primary star in the AGB phase. We investigate the mechanism of mass transfer and the process of nucleosynthesis in low-metallicity AGB stars by modelling the binary systems in which the observed CEMP-$s$ stars were formed. For this purpose we compare a sample of $67$ CEMP-$s$ stars with a grid of binary stars generated by our binary evolution and nucleosynthesis model. We classify our sample CEMP-$s$ stars in three groups based on the observed abundance of europium. In CEMP$-s/r$ stars the europium-to-iron ratio is more than ten times higher than in the Sun, whereas it is lower than this threshold in CEMP$-s/nr$ stars. No measurement of europium is currently available for CEMP-$s/ur$ stars. On average our models reproduce well the abundances observed in CEMP-$s/nr$ stars, whereas in CEMP-$s/r$ stars and CEMP-$s/ur$ stars the abundances of the light-$s$ elements are systematically overpredicted by our models and in CEMP-$s/r$ stars the abundances of the heavy-$s$ elements are underestimated. In all stars our modelled abundances of sodium overestimate the observations. This discrepancy is reduced only in models that underestimate the abundances of most of the $s$-elements. Furthermore, the abundance of lead is underpredicted in most of our model stars. These results point to the limitations of our AGB nucleosynthesis model, particularly in the predictions of the element-to-element ratios. Finally, in our models CEMP-$s$ stars are typically formed in wide systems with periods above 10000 days, while most of the observed CEMP-$s$ stars are found in relatively close orbits with periods below 5000 days.

AGB stars are responsible for producing a variety of elements, including carbon, nitrogen, and the heavy elements produced in the slow neutron-capture process ($s$-elements). There are many uncertainties involved in modelling the evolution and nucleosynthesis of AGB stars, and this is especially the case at low metallicity, where most of the stars with high enough masses to enter the AGB have evolved to become white dwarfs and can no longer be observed. The stellar population in the Galactic halo is of low mass ($\lesssim 0.85M_{\odot}$) and only a few observed stars have evolved beyond the first giant branch. However, we have evidence that low-metallicity AGB stars in binary systems have interacted with their low-mass secondary companions in the past. The aim of this work is to investigate AGB nucleosynthesis at low metallicity by studying the surface abundances of chemically peculiar very metal-poor stars of the halo observed in binary systems. To this end we select a sample of 15 carbon- and $s$-element-enhanced metal-poor (CEMP-$s$) halo stars that are found in binary systems with measured orbital periods. With our model of binary evolution and AGB nucleosynthesis, we determine the binary configuration that best reproduces, at the same time, the observed orbital period and surface abundances of each star of the sample. The observed periods provide tight constraints on our model of wind mass transfer in binary stars, while the comparison with the observed abundances tests our model of AGB nucleosynthesis.

Understanding the formation of planetary nebulae remains elusive because in the preceding asymtotic giant branch (AGB) phase these stars are heavily enshrouded in an optically thick dusty envelope. To further understand the morphology of the circumstellar environments of AGB stars we observe the closest carbon-rich AGB star IRC+10216 in scattered light. When imaged in scattered light at optical wavelengths, IRC+10216 surprisingly shows a narrow equatorial density enhancement, in contrast to the large-scale spherical rings that have been imaged much further out. We use radiative transfer models to interpret this structure in terms of two models: firstly, an equatorial density enhancement, commonly observed in the more evolved post-AGB stars, and secondly, in terms of a dust rings model, where a local enhancement of mass-loss creates a spiral ring as the star rotates. We conclude that both models can be used to reproduce the dark lane in the scattered light images, which is caused by an equatorially density enhancement formed by dense dust rather than a bipolar outflow as previously thought. We are unable to place constraints on the formation of the equatorial density enhancement by a binary system.

Aims: We study single and binary white dwarfs in the inner halo of the Milky Way in order to learn more about the conditions under which the population of halo stars was born, such as the initial mass function (IMF), the star formation history, or the binary fraction. Methods: We simulate the evolution of low-metallicity halo stars at distances up to ~ 3 kpc using the binary population synthesis code SeBa. We use two different white dwarf cooling models to predict the present-day luminosities of halo white dwarfs. We determine the white dwarf luminosity functions (WDLFs) for eight different halo models and compare these with the observed halo WDLF of white dwarfs in the SuperCOSMOS Sky Survey. Furthermore, we predict the properties of binary white dwarfs in the halo and determine the number of halo white dwarfs that is expected to be observed with the Gaia satellite. Results: By comparing the WDLFs, we find that a standard IMF matches the observations more accurately than a top-heavy one, but the difference with a bottom-heavy IMF is small. A burst of star formation 13 Gyr ago fits slightly better than a star formation burst 10 Gyr ago and also slightly better than continuous star formation $10-13$ Gyrs ago. Gaia will be the first instrument to constrain the bright end of the field halo WDLF, where contributions from binary WDs are considerable. Many of these will have He cores, of which a handful have atypical surface gravities ($\log g < 6$) and reach luminosities $\log(L/L_\odot) > 0$ in our standard model for WD cooling. These so called pre-WDs, if observed, can help us to constrain white dwarf cooling models and might teach us something about the fraction of halo stars that reside in binaries.

It is thought that type Ia supernovae (SNe Ia) are explosions of carbon-oxygen white dwarfs (CO WDs). Two main evolutionary channels are proposed for the WD to reach the critical density required for a thermonuclear explosion: the single degenerate scenario (SD), in which a CO WD accretes from a non-degenerate companion, and the double degenerate scenario (DD), in which two CO WDs merge. However, it remains difficult to reproduce the observed SN Ia rate with these two scenarios. With a binary population synthesis code we study the main evolutionary channels that lead to SNe Ia and we calculate the SN Ia rates and the associated delay time distributions. We find that the DD channel is the dominant formation channel for the longest delay times. The SD channel with helium-rich donors is the dominant channel at the shortest delay times. Our standard model rate is a factor five lower than the observed rate in galaxy clusters. We investigate the influence of ill-constrained aspects of single- and binary-star evolution and uncertain initial binary distributions on the rate of type Ia SNe. These distributions, as well as uncertainties in both helium star evolution and common envelope evolution, have the greatest influence on our calculated rates. Inefficient common envelope evolution increases the relative number of SD explosions such that for $\alpha_{\rm ce} = 0.2$ they dominate the SN Ia rate. Our highest rate is a factor three less than the galaxy-cluster SN Ia rate, but compatible with the rate determined in a field-galaxy dominated sample. If we assume unlimited accretion onto WDs, to maximize the number of SD explosions, our rate is compatible with the observed galaxy-cluster rate.

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.

In young dense clusters repeated collisions between massive stars may lead to the formation of a very massive star (above 100 Msun). In the past the study of the long-term evolution of merger remnants has mostly focussed on collisions between low-mass stars (up to about 2 Msun) in the context of blue-straggler formation. The evolution of collision products of more massive stars has not been as thoroughly investigated. In this paper we study the long-term evolution of a number of stellar mergers formed by the head-on collision of a primary star with a mass of 5-40 Msun with a lower mass star at three points in its evolution in order to better understand their evolution. We use smooth particle hydrodynamics (SPH) calculations to model the collision between the stars. The outcome of this calculation is reduced to one dimension and imported into a stellar evolution code. We follow the subsequent evolution of the collision product through the main sequence at least until the onset of helium burning. We find that little hydrogen is mixed into the core of the collision products, in agreement with previous studies of collisions between low-mass stars. For collisions involving evolved stars we find that during the merger the surface nitrogen abundance can be strongly enhanced. The evolution of most of the collision products proceeds analogously to that of normal stars with the same mass, but with a larger radius and luminosity. However, the evolution of collision products that form with a hydrogen depleted core is markedly different from that of normal stars with the same mass. They undergo a long-lived period of hydrogen shell burning close to the main-sequence band in the Hertzsprung-Russell diagram and spend the initial part of core helium burning as compact blue supergiants.

Carbon-enhanced metal-poor (CEMP) stars count for 9-25% of all the very metal-poor stars of the halo. In at least some CEMP stars the chemical enrichment is believed to be due to wind mass transfer in the past from an AGB donor star on to a low-mass companion. However, binary population synthesis models predict much lower CEMP fractions. As an alternative to the canonical Bondi-Hoyle-Lyttleton (BHL) wind accretion model, recent hydrodynamical simulations suggest an efficient mode of wind mass transfer, called wind Roche-lobe overflow (WRLOF), can reproduce observations of AGB winds in binary systems. We use our population synthesis model to test the consequences of WRLOF on a population of CEMP stars. Compared to previous predictions based on the BHL model we find a modest increase of the fraction of CEMP stars and substantial differences in the distributions of carbon and periods in the population of CEMP stars.

Carbon-enhanced metal-poor stars (CEMP) are observed as a substantial fraction of the very metal-poor stars in the Galactic halo. Most CEMP stars are also enriched in s-process elements and these are often found in binary systems. This suggests that the carbon enrichment is due to mass transfer in the past from an asymptotic giant branch (AGB) star on to a low-mass companion. Models of binary population synthesis are not able to reproduce the observed fraction of CEMP stars without invoking non-standard nucleosynthesis or a substantial change in the initial mass function. This is interpreted as evidence of missing physical ingredients in the models. Recent hydrodynamical simulations show that efficient wind mass transfer is possible in the case of the slow and dense winds typical of AGB stars through a mechanism called wind Roche-lobe overflow (WRLOF), which lies in between the canonical Bondi-Hoyle-Lyttleton (BHL) accretion and Roche-lobe overflow. WRLOF has an effect on the accretion efficiency of mass transfer and on the angular momentum lost by the binary system. The aim of this work is to understand the overall effect of WRLOF on the population of CEMP stars. To simulate populations of low-metallicity binaries we combined a synthetic nucleosynthesis model with a binary population synthesis code. In this code we implemented the WRLOF mechanism. We used the results of hydrodynamical simulations to model the effect of WRLOF on the accretion efficiency and we took the effect on the angular momentum loss into account by assuming a simple prescription. As a result the number of CEMP stars predicted by our model increases by a factor 1.2-1.8 compared to earlier results that consider the BHL prescription. Moreover, higher enrichments of carbon are produced and the final orbital period distribution is shifted towards shorter periods.

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.

Hierarchical triple systems are common among field stars yet their long-term evolution is poorly understood theoretically. In such systems Kozai cycles can be induced in the inner binary system during which the inner orbit eccentricity and the inclination between both binary orbits vary periodically. These cycles, combined with tidal friction and gravitational wave emission, can significantly affect the inner binary evolution. To investigate these effects quantitatively we perform a population synthesis study of triple systems and focus on evolutionary paths that lead to mergers of carbon-oxygen (CO) white dwarfs (WDs), which constitute an important candidate progenitor channel for type Ia supernovae (SNe Ia). We approach this problem by Monte Carlo sampling from observation-based distributions of systems within the primary mass range 1.0 - 6.5 M_Sun and inner orbit semi-major axes a_1 and eccentricities e_1 satisfying a_1 (1-e_1^2) > 12 AU, i.e. non-interacting in the absence of a tertiary component. We evolve these systems by means of a newly developed algorithm that couples secular triple dynamics with an existing binary population synthesis code. We find that the tertiary significantly alters the inner binary evolution in about 24% of all sampled systems. In particular, we find several channels leading to CO WD mergers. Amongst these is a novel channel in which a collision occurs in wide inner binary systems as a result of extremely high eccentricities induced by Kozai cycles. With assumptions on which CO WD mergers lead to a SN Ia explosion we estimate the SNe Ia delay time distribution resulting from triples and compare to a binary population synthesis study and to observations. Although we find that our triple rate is low, we have determined a lower limit of the triple-induced SNe Ia rate and further study is needed that includes triples with initially tighter inner orbits.

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.

Massive stars that lose their hydrogen-rich envelope down to a few tenths of a solar mass explode as extended type IIb supernovae, an intriguing subtype that links the hydrogen-rich type II supernovae with the hydrogen-poor type Ib and Ic. The progenitors may be very massive single stars that lose their envelope due to their stellar wind, but mass stripping due to interaction with a companion star in a binary system is currently considered to be the dominant formation channel. We computed an extensive grid of binary models with the Eggleton binary evolution code. The predicted rate from our standard models, which assume conservative mass transfer, is about 6 times smaller than the current rate indicated by observations. It is larger but still comparable to the rate expected from single stars. To recover the observed rate we must generously allow for uncertainties and low accretion efficiencies in combination with limited angular momentum loss from the system. Motivated by the claims of detection and non-detection of companions for a few IIb supernovae, we investigate the properties of the secondary star at the moment of explosion. We identify three cases: (1) the companion appears as a hot O star in about 90% of the cases, (2) the companion becomes an over-luminous B star in about 3% of the cases, or (3) the companion will appear as a K supergiant. The second case, which applies to the well-studied case of SN 1993J and possibly to SN 2001ig, is the least common case and requires that the companion very efficiently accretes the transferred material -in contrast to what is required to recover the overall IIb rate. The fast increasing number of type IIb detections may lead to more solid constraints on model uncertainties, like stellar wind-mass loss and internal mixing.

We show the preliminary results of our search for the progenitor systems of type Ia supernovae (SNe Ia). We model binary populations our aim being to compare these models with the observations of detailed element abundances of the hot Intra-Cluster Medium.

Rotation can have severe consequences for the evolution of massive stars. It is now considered as one of the main parameters, alongside mass and metallicity that determine the final fate of single stars. In massive, fast rotating stars mixing processes induced by rotation may be so efficient that helium produced in the center is mixed throughout the envelope. Such stars evolve almost chemically homogeneously. At low metallicity they remain blue and compact, while they gradually evolve into Wolf-Rayet stars and possibly into progenitors of long gamma-ray bursts. In binaries this type of evolution may occur because of (I) tides in very close binaries, as a result of (II) spin up by mass transfer, as result of (III) a merger of the two stars and (IV) when one of the components in the binary was born with a very high initial rotation rate. As these stars stay compact, the evolutionary channels are very different from what classical binary evolutionary models predict. In this contribution we discuss examples of nearly chemically homogeneous evolution in very close tidally-locked binaries. Even in such very close massive binaries, the stars may remain compact and avoid mass transfer, while Roche lobe overflow and a merger would be inevitable in the classical picture. This type of evolution may provide an alternative path to form tight Wolf-Rayet binaries and massive black hole binaries.

Abundance anomalies observed in globular cluster stars indicate pollution with material processed by hydrogen burning. Two main sources have been suggested: asymptotic giant branch stars and massive stars rotating near the break-up limit. We discuss the potential of massive binaries as an interesting alternative source of processed material. We discuss observational evidence for mass shedding from interacting binaries. In contrast to the fast, radiatively driven winds of massive stars, this material is typically ejected with low velocity. We expect that it remains inside the potential well of a globular cluster and becomes available for the formation or pollution of a second generation of stars. We estimate that the amount of processed low-velocity material that can be ejected by massive binaries is larger than the contribution of two previously suggested sources combined.

The formation of close massive black-hole binaries is a challenge for binary evolutionary models, especially the intriguing system M33 X-7 which harbours one of the most massive stellar-mass black holes (16 solar masses) orbiting a 70 solar mass O-star every 3.5 days. In standard binary evolution theory an episode of mass transfer or common envelope is inevitable in a binary with such a small orbital period, which complicates the formation of a black hole with such a high mass. To explain this system, we discuss a new binary evolution channel (De Mink et al. 2009), in which rotational mixing plays an important role. In very massive close binaries, tides force the rotation rate of the stars to be so high that rotationally induced mixing becomes very efficient. Helium produced in the center is mixed throughout the envelope. Instead of expanding during their mainsequence evolution (with the inevitable consequence of mass transfer), these stars stay compact, and avoid filling their Roche lobe. They gradually evolve into massive helium stars. This scenario naturally leads to the formation of very massive black holes in a very close orbit with a less evolved massive companion such as M33 X-7.

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.

Abundance anomalies observed in globular cluster stars indicate pollution with material processed by hydrogen burning. Two main sources have been suggested: asymptotic giant branch (AGB) stars and massive stars rotating near the break-up limit (spin stars). We propose massive binaries as an alternative source. We compute the evolution of a 20 Msun star in a close binary and find that it sheds about 10 Msun. The ejecta are enriched in He, N, Na, and Al and depleted in C and O, similar to the abundance patterns observed in gobular cluster stars. In contrast to the fast, radiatively driven winds of massive stars, this material is typically ejected with low velocity. We expect that it remains inside the potential well of a globular cluster and becomes available for the formation or pollution of a second generation of stars. We estimate that the amount of processed low-velocity material ejected by massive binaries is greater than the contribution of AGB stars and spin stars combined, assuming that the majority of massive stars in a proto-globular cluster interact. If we take the possible contribution of intermediate mass stars in binaries into account and assume that the ejecta are diluted with an equal amount of unprocessed material, we find that this scenario can provide enough material to form a second generation of low-mass stars, which is as numerous as the first generation of low-mass stars. In this scenario there is no need to make commonly adopted assumptions, such as preferential loss of the first generation of stars, external pollution of the cluster, or an anomalous initial mass function. [Abridged]

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.

In the cores of young dense star clusters repeated stellar collisions involving the same object can occur, which has been suggested to lead to the formation of an intermediate-mass black hole. In order to verify this scenario we compute the detailed evolution of the merger remnant of three sequences. We follow the evolution until the onset of carbon burning and estimate the final remnant mass to determine the ultimate fate of a runaway merger sequence. We use a detailed stellar evolution code to follow the evolution of the collision product. At each collision, we mix the two colliding stars, taking account of mass loss during the collision. During the stellar evolution we apply mass loss rates from the literature, as appropriate for the evolutionary stage of the merger remnant. We compute models for high ($Z=0.02$) and low ($Z=0.001$) metallicity to quantify metallicity effects. We find that the merger remnant becomes a Wolf-Rayet star before the end of core hydrogen burning. Mass loss from stellar winds dominates over the mass increase due to repeated mergers for all three merger sequences that we consider. In none of our high metallicity models an intermediate-mass black hole is formed, instead our models have a mass of 10--14 \Msun{} at the onset of carbon burning. For low metallicity we expect the final remnant of the merger sequence to explode as a pair creation supernova. We find that our metal-rich models become inflated as a result of developing an extended low-density envelope. This may increase the probability of further collisions, but self-consistent $N$-body calculations with detailed evolution of runaway mergers are required to verify this.

Models of rotating single stars can successfully account for a wide variety of observed stellar phenomena, such as the surface enhancements of N and He. However, recent observations have questioned the idea that rotational mixing is the main process responsible for the surface enhancements, emphasizing the need for a strong and conclusive test. We investigate the consequences of rotational mixing for massive main-sequence stars in short-period binaries. In these systems the tides spin up the stars to rapid rotation. We use a state-of-the-art stellar evolution code including the effect of rotational mixing, tides, and magnetic fields. We discuss the surface abundances expected in massive close binaries (M1~20 solar masses) and we propose using such systems to test the concept of rotational mixing. As these short-period binaries often show eclipses, their parameters can be determined with high accuracy, allowing for a direct comparison with binary evolution models. In more massive close systems (M1~50 solar masses, Porb<~2 days) we find that helium is efficiently mixed throughout the envelope. The star remains blue and compact during the main-sequence phase. It stays within its Roche lobe while it gradually becomes a helium star. It is the less massive star, in which the effects of rotational mixing are less pronounced, which fills its Roche lobe first. We propose that this evolution path provides an alternative channel for the formation of tight Wolf-Rayet binaries with a main-sequence companion and might explain massive black hole binaries such as the intriguing system M33 X-7.

One of the main uncertainties in evolutionary calculations of massive stars is the efficiency of internal mixing. It changes the chemical profile inside the star and can therefore affect the structure and further evolution. We demonstrate that eclipsing binaries, in which the tides synchronize the rotation period of the stars and the orbital period, constitute a potentially strong test for the efficiency of rotational mixing. We present detailed stellar evolutionary models of massive binaries assuming the composition of the Small Magellanic Cloud. In these models we find enhancements in the surface nitrogen abundance of up to 0.6 dex.

We present the results of binary population simulations of carbon-enhanced metal-poor (CEMP) stars. We show that nitrogen and fluorine are useful tracers of the origin of CEMP stars, and conclude that the observed paucity of very nitrogen-rich stars puts strong constraints on possible modifications of the initial mass function at low metallicity. The large number fraction of CEMP stars may instead require much more efficient dredge-up from low-metallicity asymptotic giant branch stars.

Stellar collisions are an important formation channel for blue straggler stars in globular and old open clusters. Hydrodynamical simulations have shown that the remnants of such collisions are out of thermal equilibrium, are not strongly mixed and can rotate very rapidly. Detailed evolution models of collision products are needed to interpret observed blue straggler populations and to use them to probe the dynamical history of a star cluster. We expand on previous studies by presenting an efficient procedure to import the results of detailed collision simulations into a fully implicit stellar evolution code. Our code is able to evolve stellar collision products in a fairly robust manner and allows for a systematic study of their evolution. Using our code we have constructed detailed models of the collisional blue stragglers produced in the $N$-body simulation of M67 performed by Hurley \emph{et al.} in 2005. We assume the collisions are head-on and thus ignore the effects of rotation in this paper. Our detailed models are more luminous than normal stars of the same mass and in the same stage of evolution, but cooler than homogeneously mixed versions of the collision products. The increased luminosity and inefficient mixing decrease the remaining main-sequence lifetimes of the collision products, which are much shorter than predicted by the simple prescription commonly used in $N$-body simulations.

In a companion paper we studied the detailed evolution of stellar collision products that occurred in an $N$-body simulation of the old open cluster M67 and compared our detailed models to simple prescriptions. In this paper we extend this work by studying the evolution of the collision products in open clusters as a function of mass and age of the progenitor stars. We calculated a grid of head-on collisions covering the section of parameter space relevant for collisions in open clusters. We create detailed models of the merger remnants using an entropy-sorting algorithm and follow their subsequent evolution during the initial contraction phase, through the main sequence and up to the giant branch with our detailed stellar evolution code. We compare the location of our models in a colour-magnitude diagram to the observed blue straggler population of the old open clusters M67 and NGC 188 and find that they cover the observed blue straggler region of both clusters. For M67, collisions need to have taken place recently. Differences between the evolution tracks of the collision products and normal main sequence stars can be understood quantitatively using a simple analytic model. We present an analytic recipe that can be used in an $N$-body code to transform a precomputed evolution track for a normal star into an evolution track for a collision product.

Rotational mixing is a very important but uncertain process in the evolution of massive stars. We propose to use close binaries to test its efficiency. Based on rotating single stellar models we predict nitrogen surface enhancements for tidally locked binaries. Furthermore we demonstrate the possibility of a new evolutionary scenario for very massive (M > 40 solar mass) close (P < 3 days) binaries: Case M, in which mixing is so efficient that the stars evolve quasi-chemically homogeneously, stay compact and avoid any Roche-lobe overflow, leading to very close (double) WR binaries.

A super-solar fluorine abundance was observed in the carbon-enhanced metal-poor (CEMP) star HE 1305+0132 ([F/Fe] = +2.90, [Fe/H] = -2.5). We propose that this observation can be explained using a binary model that involve mass transfer from an asymptotic giant branch (AGB) star companion and, based on this model, we predict F abundances in CEMP stars in general. We discuss wether F can be used to discriminate between the formation histories of most CEMP stars: via binary mass transfer or from the ejecta of fast-rotating massive stars. We compute AGB yields using different stellar evolution and nucleosynthesis codes to evaluate stellar model uncertainties. We use a simple dilution model to determine the factor by which the AGB yields should be diluted to match the abundances observed in HE 1305+0132. We further employ a binary population synthesis tool to estimate the probability of F-rich CEMP stars. The abundances observed in HE 1305+0132 can be explained if this star accreted 3-11% of the mass lost by its former AGB companion. The primary AGB star should have dredged-up at least 0.2 Msun of material from its He-rich region into the convective envelope via third dredge-up, which corresponds to AGB models of Z ~ 0.0001 and mass ~ 2 Msun. Many AGB model uncertainties, such as the treatment of convective borders and mass loss, require further investigation. We find that in the binary scenario most CEMP stars should also be FEMP stars, that is, have [F/Fe] > +1, while fast-rotating massive stars do not appear to produce fluorine. We conclude that fluorine is a signature of low-mass AGB pollution in CEMP stars, together with elements associated with the slow neutron-capture process.

Many binary stellar systems in which the primary star is beyond the asymptotic giant branch (AGB) evolutionary phase show significant orbital eccentricities whereas current binary interaction models predict their orbits to be circularised. We analyse how the orbital parameters in a system are modified under mass loss and mass exchange among its binary components and propose a model for enhanced mass-loss from the AGB star due to tidal interaction with its companion, which allows a smooth transition between the wind and Roche-lobe overflow mass-loss regimes. We explicitly follow its effect along the orbit on the change of eccentricity and orbital semi-major axis, as well as the effect of accretion by the companion. We calculate timescales for the variation of these orbital parameters and compare them to the tidal circularisation timescale. We find that in many cases, due to the enhanced mass loss of the AGB component at orbital phases closer to the periastron, the net eccentricity growth rate in one orbit is comparable to the rate of tidal circularisation. We show that with this eccentricity enhancing mechanism it is possible to reproduce the orbital period and eccentricity of the Sirius system, which under the standard assumptions of binary interaction is expected to be circularised. We also show that this mechanism may provide an explanation for the eccentricities of most barium star systems, which are expected to be circularised due to tidal dissipation. By proposing a tidally enhanced model of mass loss from AGB stars we find a mechanism which efficiently works against the tidal circularisation of the orbit, which explains the significant eccentricities observed in binary systems containing a white dwarf and a less evolved companion, such as Sirius and systems with barium stars.

Rapid mass transfer in a binary system can drive the accreting star out of thermal equilibrium, causing it to expand. This can lead to a contact system, strong mass loss from the system and possibly merging of the two stars. In low metallicity stars the timescale for heat transport is shorter due to the lower opacity. The accreting star can therefore restore thermal equilibrium more quickly and possibly avoid contact. We investigate the effect of accretion onto main sequence stars with radiative envelopes with different metallicities. We find that a low metallicity (Z<0.001), 4 solar mass star can endure a 10 to 30 times higher accretion rate before it reaches a certain radius than a star at solar metallicity. This could imply that up to two times fewer systems come into contact during rapid mass transfer when we compare low metallicity. This factor is uncertain due to the unknown distribution of binary parameters and the dependence of the mass transfer timescale on metallicity. In a forthcoming paper we will present analytic fits to models of accreting stars at various metallicities intended for the use in population synthesis models.

The evolution of stellar collision products in cluster simulations has usually been modelled using simplified prescriptions. Such prescriptions either replace the collision product with an (evolved) main sequence star, or assume that the collision product was completely mixed during the collision. It is known from hydrodynamical simulations of stellar collisions that collision products are not completely mixed, however. We have calculated the evolution of stellar collision products and find that they are brighter than normal main sequence stars of the same mass, but not as blue as models that assume that the collision product was fully mixed during the collision.

When two stars collide and merge they form a new star that can stand out against the background population in a starcluster as a blue straggler. In so called collision runaways many stars can merge and may form a very massive star that eventually forms an intermediate mass blackhole. We have performed detailed evolution calculations of merger remnants from collisions between main sequence stars, both for lower mass stars and higher mass stars. These stars can be significantly brighter than ordinary stars of the same mass due to their increased helium abundance. Simplified treatments ignoring this effect give incorrect predictions for the collision product lifetime and evolution in the Hertzsprung-Russell diagram.

The evolution of single stars at low metallicity has attracted a large interest, while the effect of metallicity on binary evolution remains still relatively unexplored. We study the effect of metallicity on the number of binary systems that undergo different cases of mass transfer. We find that binaries at low metallicity are more likely to start transferring mass after the onset of central helium burning, often referred to as case C mass transfer. In other words, the donor star in a metal poor binary is more likely to have formed a massive CO core before the onset of mass transfer. At solar metallicity the range of initial binary separations that result in case C evolution is very small for massive stars, because they do not expand much after the ignition of helium and because mass loss from the system by stellar winds causes the orbit to widen, preventing the primary star to fill its Roche lobe. This effect is likely to have important consequences for the metallicity dependence of the formation rate of various objects through binary evolution channels, such as long GRBs, double neutron stars and double white dwarfs.

In close binaries mass and angular momentum can be transferred from one star to the other during Roche-lobe overflow. The efficiency of this process is not well understood and constitutes one of the largest uncertainties in binary evolution. One of the problems lies in the transfer of angular momentum, which will spin up the accreting star. In very tight systems tidal friction can prevent reaching critical rotation, by locking the spin period to the orbital period. Accreting stars in systems with orbital periods larger than a few days reach critical rotation after accreting only a fraction of their mass, unless there is an effective mechanism to get rid of angular momentum. In low mass stars magnetic field might help. In more massive stars angular momentum loss will be accompanied by strong mass loss. This would imply that most interacting binaries with initial orbital periods larger than a few days evolve very non-conservatively. In this contribution we wish to draw attention to the unsolved problems related to mass and angular momentum transfer in binary systems. We do this by presenting the first results of an implementation of spin up by accretion into the TWIN version of the Eggleton stellar evolution code.

One of the major uncertainties in close binary evolution is the efficiency of mass transfer beta: the fraction of transferred mass that is accreted by a secondary star. We attempt to constrain the mass-transfer efficiency for short-period massive binaries undergoing case A mass transfer. We present a grid of about 20,000 detailed binary evolution tracks with primary masses 3.5-35 Msun, orbital periods 1-5 days at a metallicity Z=0.004, assuming both conservative and non-conservative mass transfer. We perform a systematic comparison, using least-squares fitting, of the computed models with a sample of 50 double-lined eclipsing binaries in the Small Magellanic Cloud, for which fundamental stellar parameters have been determined. About 60% of the systems are currently undergoing slow mass transfer. In general we find good agreement between our models and the observed detached systems. However, for many of the semi-detached systems the observed temperature ratio is more extreme than our models predict. For the 17 semi-detached systems that we are able to match, we find a large spread in the best fitting mass-transfer efficiency; no single value of beta can explain all systems. We find a hint that initially wider systems tend to fit better to less conservative models. We show the need for more accurate temperature determinations and we find that determinations of surface abundances of nitrogen and carbon can potentially constrain the mass-transfer efficiency further.

Thermally pulsating asymptotic giant branch (AGB) stars are the main producers of slow neutron capture (s-) process elements, but there are still large uncertainties associated with the formation of the main neutron source, 13C, and with the physics of these stars in general. Observations of s-process element enhancements in stars can be used as constraints on theoretical models. For the first time we apply stellar population synthesis to the problem of s-process nucleosynthesis in AGB stars, in order to derive constraints on free parameters describing the physics behind the third dredge-up and the properties of the neutron source. We utilize a rapid evolution and nucleosynthesis code to synthesize different populations of s-enhanced stars, and compare them to their observational counterparts to find out for which values of the free parameters in the code the synthetic populations fit best to the observed populations. These free parameters are the amount of third dredge-up, the minimum core mass for third dredge-up, the effectiveness of 13C as a source of neutrons and the size in mass of the 13C pocket. We find that galactic disk objects are reproduced by a spread of a factor of two in the effectiveness of the 13C neutron source. Lower metallicity objects can be reproduced only by lowering by at least a factor of 3 the average value of the effectiveness of the 13C neutron source needed for the galactic disk objects. Using observations of s-process elements in post-AGB stars as constraints we find that dredge-up has to start at a lower core mass than predicted by current theoretical models, that it has to be substantial ($\lambda$ >~ 0.2) in stars with mass M <~ 1.5 M_sun and that the mass of the 13C pocket must be about 1/40 that of the intershell region.

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.

We investigate the formation of the ten double-lined double white dwarfs that have been observed so far. A detailed stellar evolution code is used to calculate grids of single-star and binary models and we use these to reconstruct possible evolutionary scenarios. We apply various criteria to select the acceptable solutions from these scenarios. We confirm the conclusion of Nelemans et al. (2000) that formation via conservative mass transfer and a common envelope with spiral-in based on energy balance or via two such spiralins cannot explain the formation of all observed systems. We investigate three different prescriptions of envelope ejection due to dynamical mass loss with angular-momentum balance and show that they can explain the observed masses and orbital periods well. Next, we demand that the age difference of our model is comparable to the observed cooling-age difference and show that this puts a strong constraint on the model solutions. However, the scenario in which the primary loses its envelope in an isotropic wind and the secondary transfers its envelope, which is then re-emitted isotropically, can explain the observed age differences as well. One of these solutions explains the DB-nature of the oldest white dwarf in PG1115+116 along the evolutionary scenario proposed by Maxted et al. (2002a), in which the helium core of the primary becomes exposed due to envelope ejection, evolves into a giant phase and loses its hydrogen-rich outer layers.

We study s-process element abundance ratios in stars by carrying out stellar population synthesis, using a rapid synthetic stellar evolution code which includes an up-to-date treatment of AGB nucleosynthesis and evolution. In contrast to other studies, we find that a large spread in the 13C efficiency parameter (13C_eff) is not needed to explain the observed spread in the ratios of heavy s-process to light s-process elements ([hs/ls]), but this comes naturally from the range of different initial stellar masses and their time evolution. As a result, the 13C efficiency needed for fitting most stars in the galactic disk is constrained to 1 < 13C_eff < 2.5. In the same fashion we also study the [Pb/Ce] ratios of lead stars and find out that for low metallicities 13C_eff is approximately 0.5.

Double neutron stars (DNSs), binary systems consisting of a radio pulsar and a generally undetected second neutron star (NS), have proven to be excellent laboratories for testing the theory of general relativity. The seven systems discovered in our Galaxy exhibit a remarkably well-defined relation between the pulsar spin period and the orbital eccentricity. Here we show, using a simple model where the pulsar is spun up by mass transfer from a helium-star companion, that this relation can only be produced if the second neutron star received a kick that is substantially smaller (with a velocity dispersion of less than 50 km/s) than the standard kick received by a single radio pulsar. This demonstrates that the kick mechanism depends on the evolutionary history of the NS progenitor and that the orbital parameters of DNSs are completely determined by the evolution in the preceding helium star - neutron star phase. This has important implications for estimating the rates of NS-NS mergers, one of the major potential astrophysical sources for the direct detection of gravitational waves, and for short-period gamma-ray bursts.

The old open cluster M67 is an ideal testbed for current cluster evolution models because of its dynamically evolved structure and rich stellar populations that show clear signs of interaction between stellar, binary and cluster evolution. Here we present the first truly direct N-body model for M67, evolved from zero age to 4 Gyr taking full account of cluster dynamics as well as stellar and binary evolution. Our preferred model starts with 12000 single stars and 12000 binaries placed in a Galactic tidal field at 8.0 kpc from the Galactic Centre. Our choices for the initial conditions and for the primordial binary population are explained in detail. At 4 Gyr, the age of M67, the total mass has reduced by 90% as a result of mass loss and stellar escapes. The mass and half-mass radius of luminous stars in the cluster are a good match to observations although the model is more centrally concentrated than observations indicate. The stellar mass and luminosity functions are significantly flattened by preferential escape of low-mass stars. We find that M67 is dynamically old enough that information about the initial mass function is lost, both from the current luminosity function and from the current mass fraction in white dwarfs. The model contains 20 blue stragglers at 4 Gyr which is slightly less than the 28 observed in M67. Nine are in binaries. The blue stragglers were formed by a variety of means and we find formation paths for the whole variety observed in M67. Both the primordial binary population and the dynamical cluster environment play an essential role in shaping the population. A substantial population of short-period primordial binaries (with periods less than a few days) is needed to explain the observed number of blue stragglers in M67.

2S 0918-549 is a low-mass X-ray binary (LMXB) with a low optical to X-ray flux ratio. Probably it is an ultracompact binary with an orbital period shorter than 60 min. Such binaries cannot harbor hydrogen rich donor stars. As with other (sometimes confirmed) ultracompact LMXBs, 2S 0918-549 is observed to have a high neon-to-oxygen abundance ratio (Juett et al. 2001) which has been used to argue that the companion star is a CO or ONe white dwarf. However, type-I X-ray bursts have been observed from several of these systems implying the presence of hydrogen or helium on the neutron star surface. In this paper, we argue that the companion star in 2S 0918-549 is a helium white dwarf We first present a Type I X-ray burst from 2S 0918-549 with a long duration of 40 minutes. We show that this burst is naturally explained by accretion of pure helium at the inferred accretion rate of ~0.01 times the Eddington accretion rate. At higher accretion rates of ~0.1 Eddington, hydrogen is required to explain long duration bursts. However, at low rates the long duration is due to the large amount of helium that accumulates prior to the burst. We show that it is possible to form a helium white dwarf donor in an ultracompact binary if accretion starts during the first ascent of the giant branch, when the core is made of predominantly helium. Furthermore, this scenario naturally explains the high neon-to-oxgen ratio, without requiring a CO or ONe white dwarf companion. The only observational aspect of 2S 0918-549 that we cannot explain is the absence of helium lines in the optical spectrum. Model calculations of optical accretion disk spectra need to be carried out in order to obtain limits on the helium abundance.

A binary in which a slightly evolved star starts mass transfer to a neutron star can evolve towards ultra-short orbital periods under the influence of magnetic braking. This is called magnetic capture. In a previous paper we showed that ultra-short periods are only reached for an extremely small range of initial binary parameters, in particular orbital period and donor mass. Our conclusion was based on one specific choice for the law of magnetic braking, and for the loss of mass and angular momentum during mass transfer. In this paper we show that for less efficient magnetic braking it is impossible to evolve to ultra-short periods, independent of the amount of mass and associated angular momentum lost from the binary.

A binary in which a slightly evolved star starts mass transfer to a neutron star can evolve towards ultra-short orbital periods under the influence of magnetic braking. This is called magnetic capture. We investigate in detail for which initial orbital periods and initial donor masses binaries evolve to periods less than 30-40 minutes within the Hubble time. We show that only small ranges of initial periods and masses lead to ultra-short periods, and that for those only a small time interval is spent at ultra-short periods. Consequently, only a very small fraction of any population of X-ray binaries is expected to be observed at ultra-short period at any time. If 2 to 6 of the 13 bright X-ray sources in globular clusters have an ultra-short period, as suggested by recent observations, their formation cannot be explained by the magnetic capture model.

A binary in which a slightly evolved star starts mass transfer to a neutron star can evolve towards ultra-short orbital periods under the influence of magnetic braking. This is called magnetic capture. We investigate in detail for which initial orbital periods and initial donor masses binaries evolve to periods less than 30-40 minutes within the Hubble time. We show that only small ranges of initial periods and masses lead to ultra-short periods, and that for those only a small time interval is spent at ultra-short periods. Consequently, only a very small fraction of any population of X-ray binaries is expected to be observed at ultra-short period at any time. If 2 to 6 of the 13 bright X-ray sources in globular clusters have an ultra-short period, as suggested by recent observations, their formation cannot be explained by the magnetic capture model.

The post-Asymptotic Giant Branch (AGB) star HR4049 is in an eccentric binary system with a relatively short period probably surrounded by a dusty circumbinary disk. Extremely anomalous oxygen isotopic ratios, O16/O17 ~ O16/O18 ~ 7, have been measured from CO_2 molecules likely residing in the disk. Such a composition cannot be explained in the framework of AGB and post-AGB evolution while it can be qualitatively associated with the nucleosynthesis occurring during nova outbursts. We discuss nova models, the presence of a white dwarf companion to HR4049 and possible scenarios for the dynamical evolution of this binary system. Circumbinary disks in which mixing occurs between red-giant and nova material may also be invoked as the site of formation of some rare types of meteoritic presolar grains.

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.

With a view to understanding the formation of double neutron-stars (DNS), we investigate the late stages of evolution of helium stars with masses of 2.8 - 6.4 Msun in binary systems with a 1.4 Msun neutron-star companion. We found that mass transfer from 2.8 - 3.3 Msun helium stars and from 3.3 - 3.8 Msun in very close orbits (P_orb > 0.25d) will end up in a common-envelope (CE) and spiral-in phase due to the development of a convective helium envelope. If the neutron star has sufficient time to complete the spiraling-in process before the core collapses, the system will produce very tight DNSs (P_orb ~ 0.01d) with a merger timescale of the order of 1 Myr or less. These systems would have important consequences for the detection rate of GWR and for the understanding of GRB progenitors. On the other hand, if the time left until the explosion is shorter than the orbital-decay timescale, the system will undergo a SN explosion during the CE phase. Helium stars with masses 3.3 - 3.8 Msun in wider orbits (P_orb > 0.25d) and those more massive than 3.8 Msun do not go through CE evolution. The remnants of these massive helium stars are DNSs with periods in the range of 0.1 - 1 d. This suggests that this range of mass includes the progenitors of the galactic DNSs with close orbits (B1913+16 and B1534+12). A minimum kick velocity of 70 km/s and 0 km/s (for B1913+16 and B1534+12, respectively) must have been imparted at the birth of the pulsar's companion. The DNSs with wider orbits (J1518+4904 and probably J1811-1736) are produced from helium star-neutron star binaries which avoid RLOF, with the helium star more massive than 2.5 Msun. For these systems the minimum kick velocities are 50 km/s and 10 km/s (for J1518+4904 and J1811-1736, respectively).

We summarize the main results from MODEST-1, the first workshop on MOdeling DEnse STellar systems. Our goal is to go beyond traditional population synthesis models, by introducing dynamical interactions between single stars, binaries, and multiple systems. The challenge is to define and develop a software framework to enable us to combine in one simulation existing computer codes in stellar evolution, stellar dynamics, and stellar hydrodynamics. With this objective, the workshop brought together experts in these three fields, as well as other interested astrophysicists and computer scientists. We report here our main conclusions, questions and suggestions for further steps toward integrating stellar evolution and stellar (hydro)dynamics.

The progenitors of Type Ia supernovae (SNe Ia) have not been identified. Though they are no longer fashionable we investigate the consequences if a significant number of SNe Ia were edge-lit detonations (ELDs) of carbon/oxygen white dwarfs that have accreted a critical mass of helium. Our best understanding of the Phillips relation between light curve speed and peak luminosity assigns both these phenomena to the amount of $^{56}$Ni produced. In ELDs there are two sites of $^{56}$Ni synthesis. If the peak luminosity is determined primarily by the C/O ratio in the core it is primarily a function of its progenitor's initial mass. If the light curve decay speed is determined by the total mass of iron group elements ejected this is a function of the total mass of the ELD at the time of explosion. In general these two masses are correlated and an empirical relation between peak luminosity and light curve shape can be expected. However when we perform population synthesis for progenitors of different metallicities we find a systematic shift in this relation that would make distant SNe Ia fainter than those nearby. The abundances of alpha-rich isotopes, such as $^{44}$Ca, in the solar system indicate that only about 40 per cent of SNe Ia are edge-lit so any systematic effect that could be present would be correspondingly diluted. If we examine only the small subset of ELDs that accrete from a naked helium star, rather than a He white dwarf, the systematic effect disappears.

We present new evolutionary sequences for low and intermediate mass stars for three different metallicities, Z = 0.02,0.008 and 0.004. We evolve the models from the pre-main sequence to the thermally-pulsing asymptotic giant branch phase. We have two sequences of models for each mass, one which includes mass loss and one without mass loss. Typically 20 or more pulses have been followed for each model, allowing us to calculate the third dredge-up parameter for each case. Using the results from this large and homogeneous set of models, we present an approximate fit for the core mass at the first thermal pulse, as well as for the third dredge-up efficiency parameter, and the core mass at the first third dredge-up episode as a function of metallicity and total mass. We also examine the effect of a reduced envelope mass on the value of the third dredge-up efficiency parameter.

Recently the observationally derived stellar-wind mass-loss rates for Wolf-Rayet stars, or massive naked helium stars, have been revised downwards by a substantial amount. We present evolutionary calculations of helium stars incorporating such revised mass-loss rates, as well as mass transfer to a close compact binary companion. Our models reach final masses well in excess of 10 Msun, consistent with the observed masses of black holes in X-ray binaries. This resolves the discrepancy found with previously assumed high mass-loss rates between the final masses of stars which spend most of their helium-burning lifetime as Wolf-Rayet stars (about 3 Msun) and the minimum observed black hole masses (6 Msun). Our calculations also suggest that there are two distinct classes of progenitors for Type Ic supernovae: one with very large initial masses (> 35 Msun), which are still massive when they explode and leave black hole remnants, and one with moderate initial masses (about 12 - 20 Msun) undergoing binary interaction, which end up with small pre-explosion masses and leave neutron star remnants.

The evolution of helium stars with masses of 1.5 - 6.7 M_sun in binary systems with a 1.4 M_sun neutron-star companion is presented. Such systems are assumed to be the remnants of Be/X-ray binaries with B-star masses in the range of 8 - 20 M_sun which underwent a case B or case C mass transfer and survived the common-envelope and spiral-in process. The orbital period is chosen such that the helium star fills its Roche lobe before the ignition of carbon in the centre. We distinguish case BA (in which mass transfer is initiated during helium core burning) from case BB (onset of Roche-lobe overflow occurs after helium core burning is terminated, but before the ignition of carbon). We found that the remnants of case BA mass transfer from 1.5 - 2.9 M_sun helium stars are heavy CO white dwarfs. This implies that a star initially as massive as 12 M_sun is able to become a white dwarf. CO white dwarfs are also produced from case BB mass transfer from 1.5 - 1.8 M_sun helium stars, while ONe white dwarfs are formed from 2.1 - 2.5 M_sun helium stars. Case BB mass transfer from more-massive helium stars with a neutron-star companion will produce a double neutron-star binary. We are able to distinguish the progenitors of type Ib supernovae (as the high-mass helium stars or systems in wide orbits) from those of type Ic supernovae (as the lower-mass helium stars or systems in close orbits). Finally, we derive a "zone of avoidance" in the helium star mass vs. initial orbital period diagram for producing neutron stars from helium stars.

We present a rapid binary evolution algorithm that enables modelling of even the most complex binary systems. In addition to all aspects of single star evolution, features such as mass transfer, mass accretion, common-envelope evolution, collisions, supernova kicks and angular momentum loss mechanisms are included. In particular, circularization and synchronization of the orbit by tidal interactions are calculated for convective, radiative and degenerate damping mechanisms. We use this algorithm to study the formation and evolution of various binary systems. We also investigate the effect that tidal friction has on the outcome of binary evolution. Using the rapid binary code, we generate a series of large binary populations and evaluate the formation rate of interesting individual species and events. By comparing the results for populations with and without tidal friction we quantify the hitherto ignored systematic effect of tides and show that modelling of tidal evolution in binary systems is necessary in order to draw accurate conclusions from population synthesis work. Tidal synchronism is important but because orbits generally circularize before Roche-lobe overflow the outcome of the interactions of systems with the same semi-latus rectum is almost independent of eccentricity. It is not necessary to include a distribution of eccentricities in population synthesis of interacting binaries, however, the initial separations should be distributed according to the observed distribution of semi-latera recta rather than periods or semi-major axes.

We present a state-of-the-art N-body code which includes a detailed treatment of stellar and binary evolution as well as the cluster dynamics. This code is ideal for investigating all aspects relating to the evolution of star clusters and their stellar populations. It is applicable to open and globular clusters of any age. We use the N-body code to model the blue straggler population of the old open cluster M67. Preliminary calculations with our binary population synthesis code show that binary evolution alone cannot explain the observed numbers or properties of the blue stragglers. On the other hand, our N-body model of M67 generates the required number of blue stragglers and provides formation paths for all the various types found in M67. This demonstrates the effectiveness of the cluster environment in modifying the nature of the stars it contains and highlights the importance of combining dynamics with stellar evolution. We also perform a series of N = 10000 simulations in order to quantify the rate of escape of stars from a cluster subject to the Galactic tidal field.

We present analytic formulae that approximate the evolution of stars for a wide range of mass and metallicity. Stellar luminosity, radius and core mass are given as a function of age, M and Z, for all phases from the zero-age main-sequence up to, and including, the remnant stages. For the most part we find continuous formulae accurate to within 5% of detailed models. These formulae are useful for purposes such as population synthesis that require very rapid but accurate evaluation of stellar properties, and in particular for use in combination with N-body codes. We describe a mass loss prescription that can be used with these formulae and investigate the resulting stellar remnant distribution.

We find the distribution of coalescence times, birthrates, spatial velocities, and subsequent radial offsets of coalescing neutron stars (NSs) in various galactic potentials accounting for large asymmetric kicks introduced during a supernovae. The birthrates of bound NS-NS binaries are quite sensitive to the magnitude of the kick velocities but are, nevertheless, similar (~10 per Galaxy per Myr) to previous population synthesis studies. The distribution of merger times since zero-age main sequence is, however, relatively insensitive to the choice of kick velocities. With a median merger time of ~100 Myr, we find that compact binaries should closely trace the star formation rate in the Universe. In a range of plausible galactic potentials (M_galaxy > 3 x 10^10 M_solar) the median radial offset of a NS-NS mergers is less than 10 kpc. At a redshift of z=1 (with H_0 = 65 km/s/Mpc and Omega = 0.2), this means that half the coalescences should occur within ~1.3 arcsec from the host galaxy. In all but the most shallow potentials, ninety percent of NS-NS binaries merge within 30 kpc of the host. We find that although the spatial distribution of coalescing neutron star binaries is consistent with the close spatial association of known optical afterglows of gamma-ray bursts (GRBs) with faint galaxies, a non-negligible fraction (~15 percent) of GRBs should occur well outside (>30 kpc) dwarf galaxy hosts. Extinction due to dust in the host, projection of offsets, and a range in interstellar medium densities confound the true distribution of NS-NS mergers around galaxies with an observable set of optical transients/galaxy offsets.

We present a simple and efficient, yet reasonably accurate, equation of state, which at the moderately low temperatures and high densities found in the interiors of stars less massive than the Sun is substantially more accurate than its predecessor by Eggleton, Faulkner & Flannery. Along with the most recently available values in tabular form of opacities, neutrino loss rates, and nuclear reaction rates for a selection of the most important reactions, this provides a convenient package of input physics for stellar modelling. We briefly discuss a few results obtained with the updated stellar evolution code.