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.

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]