Issue 
A&A
Volume 626, June 2019



Article Number  A68  
Number of page(s)  19  
Section  Stellar structure and evolution  
DOI  https://doi.org/10.1051/00046361/201834598  
Published online  14 June 2019 
Slowly, slowly in the wind
3D hydrodynamical simulations of wind mass transfer and angularmomentum loss in AGB binary systems
^{1}
Department of Astrophysics/IMAPP, Radboud University, PO Box 9010, 6500 GL Nijmegen, The Netherlands
email: m.saladino@astro.ru.nl, o.pols@astro.ru.nl
^{2}
Leiden Observatory, Leiden University, PO Box 9513, 2300 RA Leiden, The Netherlands
^{3}
ArgelanderInstitute für Astronomie (AIfA), University of Bonn, Auf dem Hügel 71, 53121 Bonn, Germany
email: cabate@unibonn.de
Received:
7
November
2018
Accepted:
11
March
2019
Wind mass transfer in binary systems with asymptotic giant branch (AGB) donor stars plays a fundamental role in the formation of a variety of objects, including barium stars and carbonenhanced metalpoor (CEMP) stars. In an attempt to better understand the properties of these systems, we carry out a comprehensive set of smoothedparticle hydrodynamics (SPH) simulations of windlosing 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 the ratio of wind velocity to orbital velocity (v_{∞}/v_{orb}) and on 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 angularmomentum 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 and the angularmomentum 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, spinorbit 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 Rochelobe overflow. In this case the massaccretion efficiency and angularmomentum loss differ from those found for a nonrotating donor. Finally, we provide relations for the massaccretion efficiency and angularmomentum loss as a function of v_{∞}/v_{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, such as cataclysmic binaries and type Ia supernovae.
Key words: binaries: general / stars: AGB and postAGB / hydrodynamics / methods: numerical / stars: winds, outflows / stars: lowmass
© ESO 2019
1. Introduction
In recent years, major efforts have been made to characterise the properties of the outflow of nearby asymptotic giant branch (AGB) stars (e.g. Decin et al. 2010, 2015; Blind et al. 2011; Jeffers et al. 2014; Kervella et al. 2016). These studies shed unprecedented light on the structure of the winds, their chemical composition and dusttogas ratios, wind velocity profiles, and massloss rates. In many cases they reveal the presence of a binary companion that interacts with the outflow and modifies its structure (Blind et al. 2011; Jeffers et al. 2014; Decin et al. 2015; Bujarrabal et al. 2018). As many AGB stars are known to be in binary systems, several authors have attempted to model the complex interactions between the outflow and the two stellar components, with the aim of reproducing their observed morphologies (e.g. Mohamed & Podsiadlowski 2007; Kim et al. 2017; Ramstedt et al. 2017), and to understand how the evolution of the system is affected. To investigate this evolution, several authors have studied what fraction of wind material can be transferred onto the lowmass companion (Theuns et al. 1996; Nagae et al. 2004; de ValBorro et al. 2009, 2017; Mohamed & Podsiadlowski 2007, 2012; Liu et al. 2017; Chen et al. 2017) and how much angular momentum is lost by wind ejection (Brookshaw & Tavani 1993; Jahanara et al. 2005; Chen et al. 2018; Saladino et al. 2018).
A thorough understanding of wind mass transfer in binary systems is of fundamental importance to determine the formation mechanism of a variety of objects, such as symbiotic binaries (Merrill 1942, 1944, 1948, 1950; Kenyon 1992), barium (Ba) stars (Bidelman & Keenan 1951), carbon and sprocess enhanced metalpoor (CEMPs) stars (Beers & Christlieb 2005), and blue stragglers in old open clusters (Mathieu & Geller 2015). Binary populationsynthesis studies typically assume that the AGB star loses mass in the form of a fast spherically symmetric wind and that the mass accretion process is described by the BondiHoyleLyttlelton (BHL) formalism (Hoyle & Lyttleton 1939; Bondi & Hoyle 1944; Edgar 2004). In this simplified model, the orbit of the binary expands in response to mass loss. As a consequence of these assumptions, theoretical studies of the progeny of AGB binary systems, such as CEMPs and Ba stars, predict a dearth of systems with orbital periods between a few hundred and a few thousand days (Pols et al. 2003; Izzard et al. 2009, 2010; Nie et al. 2012; Abate et al. 2018) which is at odds with the observations (e.g. Jorissen et al. 1998, 2016; Hansen et al. 2016). In addition, in order to reproduce the chemical abundances of some of the progeny of AGB binary systems, massaccretion efficiencies much higher than those predicted by the BHL formalism are needed (see e.g. Abate et al. 2015a,b, for CEMP stars).
In Saladino et al. (2018, hereafter Paper I), we performed smoothedparticle hydrodynamic (SPH) simulations including cooling of the gas to model interactions of nonrotating AGB binary systems via their winds. We used a predefined constant wind velocity and adopted the same initial stellar parameters of the binary components in all models. We studied the effect of the terminal velocity of the wind on the massaccretion efficiency and on the angularmomentum loss. We found that these quantities show a strong dependence on the ratio of the wind velocity to the relative orbital velocity, v_{∞}/v_{orb}. Despite differences in the assumptions, by comparing our models to the hydrodynamical models by Jahanara et al. (2005) and Chen et al. (2017), we found a similar trend for the angularmomentum loss as a function of v_{∞}/v_{orb}. In Paper I, we also presented a preliminary analytical expression that describes the angularmomentum loss in terms of v_{∞}/v_{orb}.
In the context of CEMPs stars, Abate et al. (2018) performed binary population synthesis simulations in an attempt to reproduce the observed orbitalperiod distribution of the CEMP stars in the sample obtained by Hansen et al. (2016). In the simulations of Abate et al. (2018) massaccretion efficiency and angularmomentum are treated independently because hydrodynamical simulations computing simultaneously the mass transfer rates and the angularmomentum loss via winds were not available at the time. In their model sets M6 and M7, Abate et al. (2018) implement the fit we provide in Paper I for angularmomentum loss, and adopt the wind Rochelobe overflow (WRLOF) model for the massaccretion efficiency (Mohamed & Podsiadlowski 2007; Abate et al. 2013). With these assumptions they find a slight increase in the proportion of evolved binaries at short orbital periods (<2500 day) compared to their default model set assuming isotropic winds; however, the period distribution of their synthetic population does not reproduce the observations. A selfconsistent model in which angularmomentum loss and massaccretion efficiencies are calculated from simulations where the same set of conditions are adopted may help to improve the results of binary population synthesis codes on the populations of CEMPs stars. In addition, most hydrodynamical studies have focused on investigating the consequences of modifying one orbital parameter at a time on the massaccretion efficiency or angularmomentum loss. No hydrodynamical study has yet been performed with the same set of assumptions for a large grid of binary and wind parameters. The aim of this work is to perform these hydrodynamical simulations for a much wider range of binary and wind parameters in order to have a better understanding of the important quantities driving the interaction of AGB winds with the binary systems.
Another aspect that may be of importance and that is sometimes overlooked is the effect of the spin of the donor star on the interacting binary. During the AGB phase the rotational velocity of a single star is greatly reduced due to the expansion of its envelope and because angular momentum is carried away by the mass the star is losing. For this reason, the angular velocity of a single AGB star is negligibly low. However, when the AGB star is in a binary system tidal effects can trigger spinorbit coupling. In the equilibriumtide model the synchronisation timescale depends on the ratio (R_{d}/a)^{−6}, where R_{d} is the radius of the donor star and a the orbital separation of the binary (Zahn 1977). This implies that for AGB binaries where the donor star is close to filling its Roche lobe, tidal effects can be important since the synchronisation timescale is usually much shorter than the typical evolution timescale of an AGB star. Most of the previous hydrodynamical studies of AGB binaries are performed in the corotating frame, thus corotation of the AGB star is implicit in the simulations. Only a few works (e.g. Theuns et al. 1996; Saladino et al. 2018) have performed hydrodynamical simulations in the inertial frame of the binary system, assuming the donor star is not rotating. The morphologies found in some of the models where corotation of the donor is assumed correspond to the WRLOF geometry, i.e. the dense wind is confined to the Roche lobe of the donor star and flows through the inner Lagrange point towards the companion star (e.g. de ValBorro et al. 2009; Mohamed & Podsiadlowski 2012; Chen et al. 2017). In contrast, the morphologies found when considering a nonrotating donor star are somewhat different, mainly showing spiral patterns formed by the accretion wake of the companion star. The models by Mohamed & Podsiadlowski (2012) and Chen et al. (2017) show that when interaction of the binary stars occurs via WRLOF the massaccretion efficiency onto the companion star is enhanced. Considering the differences in the results found in numerical models, it is necessary to better understand how rotation modifies the morphology of the AGB outflow in binary systems, and to constrain the range of orbital parameters within which it is realistic to assume corotation of the AGB donor.
In this work we present a comprehensive study of how angularmomentum loss and accretion efficiency in lowmass binary stars vary as a function of mass ratio, orbital separation, and windtoorbital velocity ratios. We provide two analytical expressions for angularmomentum loss and massaccretion efficiency that can be easily implemented in binary population synthesis codes. In addition, we explore how the results are changed if we assume that the AGB donor is corotating with the binary system instead of nonrotating, and under what circumstances corotation can be expected.
2. Method
2.1. Numerical method
The numerical method used in this paper for the SPH simulations is similar to that described in Paper I. Here we only give a brief summary. We use the Astrophysical Multipurpose Software Environment AMUSE (Portegies Zwart et al. 2009, 2013; Pelupessy et al. 2013) to couple the smoothedparticle hydrodynamics (SPH) code FI (Hernquist & Katz 1989; Gerritsen & Icke 1997; Pelupessy et al. 2004) with the Nbody code HUAYNO (Pelupessy et al. 2012) using the BRIDGE module (Fujii et al. 2007). We use the Nbody code to evolve the orbits of the stars. We choose an SPH code because these codes are arguably better at conserving angular momentum than gridbased codes (Price 2012). Unlike in Paper I, here we allow the stars to feel the gravitational field of the gas and the selfgravity of the gas is taken into account.
The AGB wind is modelled using the “accelerating wind” mode of the STELLAR_WIND.PY routine (van der Helm et al. 2019, also available in AMUSE). Wind particles are created in a shell around the star losing mass and are injected with initial velocity v_{init}. The initial temperature of the wind particles is equal to the effective temperature T_{eff} of the donor star. The acceleration given to the SPH particles balances the gravity of the donor star:
Due to gas pressure, the wind will experience an outward acceleration close to the donor star until it reaches its terminal velocity v_{∞} > v_{init} at some distance from the surface of the AGB star. This acceleration mechanism differs from that presented in Paper I where the acceleration was set to compensate for each term in the equation of motion of the gas, such that the velocity of the wind was constant as a function of distance form the star (dv/dr = 0).
Here, as in Paper I, the gas is assumed to be monoatomic, with an adiabatic index γ = 5/3. The equation governing the cooling or heating of the gas is given by
where k_{B} is the Boltzmann constant, T is the gas temperature, ρ is the gas density, μ is the mean molecular weight, and C a constant parameter which value is 10^{−5} g s cm^{−3}. The first term assumes that cooling comes from gas radiating away or absorbing thermal energy trying to reach the equilibrium temperature at radius r given by the Eddington approximation for a grey spherical atmosphere (Chandrasekhar 1934)
where κ_{g} is the gas opacity and κ_{d} is the dust opacity. Unlike in Paper I, where we assumed a constant value of the total opacity (κ_{g} + κ_{d}), in this work we only allow for the dust opacity to come into play at distances larger than three times the radius of the donor star (3R_{d}), where dust is expected to form according to AGB wind models (Höfner & Olofsson 2018). The values of the opacities are constant in the calculations and are equal to κ_{g} = 2 × 10^{−4} cm^{2} g^{−1} and κ_{d} = 5 cm^{2} g^{−1} (Bowen 1988). The second term in Eq. (2) corresponds to the cooling rate prescription for high temperatures (log T ≥ 3.8 K) of Schure et al. (2009). Because the abundances in this paper are different from those in Paper I (see Sect. 2.2), we compute a new cooling table using Eq. (3) of Schure et al. (2009) and the abundances n_{i} from Table A.1.
The resolution in SPH is defined by the smoothing length h_{i} ∝ (m_{g}/ρ_{i})^{1/3}, where m_{g} is the SPH gas particle mass and ρ_{i} the density of the gas at the position r_{i} of the particle i. In Paper I, we found that the angularmomentum loss is independent of the resolution used, while the accretion efficiency is only weakly dependent on the resolution (see Fig. 2 of that paper). Because of this, and in order to minimise the computational time, the resolution we use in this work corresponds to that used in model R3 of Paper I. For an orbital separation of 5 AU, and scaling the density to the value of the massloss rate Ṁ_{d} used (see Sect. 2.2), the corresponding mass of an SPH particle is m_{g} = 1.2 × 10^{−9} M_{⊙}. Furthermore, in order to optimise the numerical computation at large orbital separations, we choose the typical smoothing length to be proportional to the semimajor axis of the binary a. Since the average gas density is expected to decrease with the inverse square of the distance to the donor star, i.e. to scale with a^{−2}, we can achieve this by choosing the SPH particle mass m_{g} to be proportional to a.
As seen in Paper I (see Fig. 1 of that paper) the transfer of angular momentum from the binary orbit to the gas occurs within a few times the orbital separation. For this reason, and in order to minimise the computational time, we remove the SPH particles once they have crossed a boundary at 5a from the centre of mass of the system. Finally, as in Paper I, the values of the artificial viscosity parameters are α_{SPH} = 0.5 and β_{SPH} = 1. As shown in Paper I, adopting different values for these parameters does not modify the results in a significant way.
2.2. Binary setup
The motivation of this paper is to study the effect of angularmomentum loss and massaccretion efficiency for a range of parameters relevant for lowmass binary stars. In this work we explore the lowmetallicity case of the donor star, which applies to the progenitors of CEMPs stars. The low metallicity impacts several properties of the donor star during the AGB phase, such as stellar radius, effective temperature, and consequently the massloss rate. In order to obtain realistic parameters for the AGB donor, we compute the evolution of a single star of initial mass M_{d, init} = 1.5 M_{⊙} and initial metallicity Z_{init} = 10^{−4} using the population nucleosynthesis code binary_c^{1} (Izzard et al. 2004, 2006, 2009, 2018; Abate et al. 2013). The initial metal abundances of the star are equal to solar (Asplund et al. 2009) scaled down to Z_{init}. The massloss rate of the AGB star is calculated using the method described in Vassiliadis & Wood (1993). As starting conditions for the simulations, we adopt the stellar mass M_{d}, radius R_{d}, effective temperature T_{eff}, massloss rate Ṁ_{d}, and chemical abundances at the moment the star has reached the superwind phase on the AGB, and the radius is close to its maximum value. This is the phase of evolution when the strongest interaction between the wind and the companion occurs. At this moment the mass is reduced to 1.2 M_{⊙}, its age is 1796 Myr, and its surface is strongly enriched in carbon (see Tables 1 and A.1).
Constant parameters of the donor star in the hydrodynamical simulations.
Similar to Paper I, the companion star is modelled as a sink particle with radius equal to 0.1R_{L, 2}, where R_{L, 2} is the Rochelobe radius of the accretor star. For the companion star mass M_{a}, we take M_{a} = 1.2 M_{⊙}, 0.9 M_{⊙}, 0.6 M_{⊙}, 0.4 M_{⊙}, and 0.3 M_{⊙}. The mass ratio q = M_{d}/M_{a} varies accordingly between 1 and 4. These mass ratios encompass the mass ratios of CEMPs progenitors from Abate et al. (2015a,c). We should note that these are not the initial mass ratios, but the mass ratios at the time when AGB mass transfer takes place.
Using the setup described above, we performed 31 simulations (Table 2) in which the initial orbit of the binary system is circular and its initial orbital separation ranges between 4 and 20 AU. These orbital separations are chosen so that they correspond to the typical orbital periods of CEMPs progenitor stars found by Abate et al. (2015a),c). The stars are placed in the inertial frame and the centre of mass of the binary system is positioned at the origin. We simulate models with different initial velocities of the wind v_{init} = 1, 5 and 12 km s^{−1}. We note that in the models with the closest orbits (4−5 AU) the donor star fills a substantial fraction of its Roche lobe and its surface layers and atmosphere will be deformed by tides from the companion to some extent. We do not take this tidal deformation into account in the simulations since the wind particles are always injected from a spherical shell with inner radius R_{d} (see Sect. 2.1).
Parameters of the simulations.
2.3. Rotation and spinorbit coupling
In this paper we explore the possibility of rotation of the AGB donor. In order to check for which orbital separations corotation of the AGB donor with the orbit is expected, we perform a set of simulations within binary_c of binary systems interacting via isotropic winds, and we compute their tidal evolution. The physics of tidal evolution in binary_c is implemented following the work of Hurley et al. (2002, Sect. 2.3 and references therein) and the equation for the tidal circularisation and synchronisation timescale is based on Rasio et al. (1996). The equations of tidal evolution are originally reported by Hut (1981). We evolve systems with initial orbital separations between 3 and 30 AU, initial masses as in Sect. 2.2 and with the initial angular velocity of the donor star Ω_{spin} set equal to zero, i.e. the star is initially nonrotating. We check the angular velocity of the donor star at the time the mass of the AGB star and the mass ratios are equal to those assumed in our hydrodynamical simulations. We find that for orbital separations up to ≈15−20 AU, depending on the mass ratio, tidal interactions have spun up the donor star to a nonzero angular velocity. For initial separations smaller than 7−10 AU a state of nearcorotation, with Ω_{spin} comparable to the angular velocity of the binary Ω_{bin} is reached at some point during the AGB phase. However, except for very close orbital separations (a ≲ 4−6 AU), corotation is lost again by the time we start our hydrodynamical simulations as a result of the high massloss rate of the AGB star. The ratio Ω_{spin}/Ω_{bin} at this moment is a decreasing function of separation and mass ratio (see Appendix D for details).
To account for spinorbit coupling, and in order to investigate the effect of rotation on the angularmomentum loss and massaccretion efficiency, in a subset of our simulations we assume both nonrotation and corotation with the orbit for models with the same mass ratio and the same orbital separation (see Table 2; the corotating models are designated by the symbol Ω). In the models where corotation of the binary is assumed, we add a tangential component to the velocity of the wind particles as they leave the donor star, v_{T} = Ω_{spin} × r, with Ω_{spin} = Ω_{bin} and r the position of the gas particles with respect to the centre of mass of the donor star. For the other simulations we assume the donor star to be nonrotating even if Ω_{spin} may be nonzero.
3. Results
3.1. Terminal velocity of the wind
As mentioned in Sect. 2.1, in this work the acceleration of the wind only balances the gravitational deceleration by the AGB star. For this reason, the terminal velocity of the wind is different from that at which the wind is injected. In Paper I, where we imposed a constant velocity on the wind particles with a predefined terminal velocity, we only explored velocities as low as 10 km s^{−1} for numerical reasons^{2}. Here we explore lower values of the initial wind velocities. We determine the terminal velocities that can be achieved by performing simulations of single AGB stars. In these models the star is allowed to eject matter for 40 yr in order to reach stable conditions. The average terminal velocities are measured at a distance of 25R_{d}, and averaged over the last 1.5 yr. Figure 1 shows the velocity and density profiles of the wind as a function of distance for a single star (see Table 1 for stellar parameters) and an initial velocity of the wind of v_{init} = 12 km s^{−1}. Because of gas pressure, the wind accelerates to an average terminal velocity v_{∞} = 15.1 km s^{−1}. For v_{init} = 5 km s^{−1}, we find that the average v_{∞} = 10.8 km s^{−1} and for v_{init} = 1 km s^{−1}, v_{∞} = 6.0 km s^{−1}. The injection of particles with such low initial wind velocities allows the wind to achieve terminal velocities that lie within the observed range given by Vassiliadis & Wood (1993) for the winds of AGB stars. We describe the results of our binary simulations in terms of the terminal velocity achieved by the wind of a single star with the same initial velocity as assumed for the donor star.
Fig. 1.
Velocity and density profiles of the SPH gas particles leaving a single star with v_{init} = 12 km s^{−1}. The top dashed line corresponds to a constant velocity of 15.1 km s^{−1}. The lower dashed line is the corresponding analytical density in a steady state. The red rectangle delimits the inner boundary of the SPH particles around the donor star. 
3.2. Morphology of the outflow
In the first part of this section we present the results of the outflow morphologies for the nonrotating models. In the second part, we discuss the differences between the models in which corotation of the AGB donor is assumed and the corresponding models in which rotation is neglected.
3.2.1. Nonrotating models
The outflow morphologies in the simulations with a nonrotating donor star are very similar to those found in Paper I, and are in agreement with the results of similar works (e.g. Mastrodemos & Morris 1998; Liu et al. 2017). In general, for all the models we observe two spiral arms wrapped around the binary system. These spiral arms delimit the accretion wake behind the companion star. In the following we describe the observed morphologies as a function of the orbital separation (Fig. 2), the mass ratio (Fig. 3), and the wind velocity (Fig. 4).
Fig. 2.
Gas density in the orbital plane for models Q1a5, Q1a10, and Q1a15 (same mass ratio, but different orbital separations). The companion star is located on the left and the AGB donor on the right. Each figure is scaled to display the full simulation, i.e. a sphere with radius equal to 5a, and shows the situation after 9.5 orbital periods. 
Fig. 3.
Same as Fig. 2, but for models with the same orbital separation (a = 5 AU) and different mass ratios. 
Fig. 4.
Gas density on the orbital plane for for models with the same orbital separation and mass ratio (a = 5 AU, q = 1), but different initial wind velocities (1 km s^{−1} and 12 km s^{−1}). For clarity, the density scale here is different from the previous figures. For model Q1a5 (right) this figure is a closeup version of Fig. 3a. 
Figure 2 shows the gas density profile in the orbital plane for models Q1a5, Q1a10, and Q1a20, which have the same mass ratio (q = 1) but different orbital separations (a = 5, 10, and 20 AU), at the same orbital phase. We observe that the complexity of the morphology decreases with increasing orbital separations, which suggests that as the orbital separation increases, the interaction between the companion star and the wind becomes less strong. In relatively close systems, the gas density in the spiral arms is higher than in wider systems. Also, the spiral arms that delimit the accretion wake behind the companion star are less tightly wound around the binary for wider systems. One more feature that varies in these models is the opening angle θ of the accretion wake. We find that for close orbital separations θ is large, and decreases when the orbital separation increases. This occurs because θ is a function of the Mach number ℳ, i.e. θ = sin^{−1}ℳ^{−1} (see Appendix B). We find that for the model with a = 5 AU, ℳ ≈ 2 in the vicinity of the accretor, while the Mach number increases to ℳ ≈ 4.5 and ℳ ≈ 5 for the systems with a = 10 AU and a = 20, respectively. In addition, we find that for model Q1a5 an accretion disk is formed that is not observed in the wider models. However, as is extensively discussed in Paper I, since the radius of the sink particle is large, the disk may not be able to form in these cases.
The same trend in the morphology of the outflow with orbital separation is found for systems with higher mass ratios. For close orbital separations a stronger interaction between the gas and the stars occurs, as shown by the presence of an accretion disk around the companion star and dense tightly wound spiral patterns around the binary.
In order to illustrate the dependence of the morphology on the mass ratio, Fig. 3 shows the density profile in the orbital plane for models Q13a5, Q2a5, Q3a5, and Q4a5 (same orbital separation, but different mass ratio). This can be compared with model Q1a5 (Fig. 2a), which has the same orbital separation. However, these models do not show an accretion disk around the companion star, for the reasons explained above. Overall, the gas densities are similar and all the models have two arms tightly wound around the binary. The main differences occur in the accretion wake behind the companion star, which has lower density and a smaller opening angle as the mass ratio increases. Similarly to Liu et al. (2017) ^{3}, we find that for higher mass ratios the spiral arm is less tightly wound. These differences in the geometry of the outflow can be explained by the lower gravity of the companion star. When the companion star is more massive, its gravitational attraction is stronger and so is the spiral shock formed in the vicinity of the accretor. At the same time, the opening angle of the spiral shock becomes wider. This is explained by the dependence on the Mach number, as discussed above. As the mass ratio increases from q = 1 to q = 4, the Mach number near the companion increases from ℳ ≈ 2.5 to ℳ ≈ 5. On the other hand, as the mass ratio increases, the density along the inner spiral arm decreases. As is shown in Sect. 3.3, a lower density in the accretion wake implies a weaker torque on the binary system, which affects the amount of angular momentum exchanged between the orbit and the gas.
In general we observe that as the orbital separation increases, the structures in the morphology of the outflow become less prominent and less complex. This suggests that the interaction between the gas and the stars decreases as a function of distance. A similar trend was found in Paper I, where the free parameter in the simulations was the velocity of the wind. In the models discussed above the terminal velocity of the wind is fixed and we vary the orbital separation, which is equivalent to changing the orbital velocity of the system. As a result, the ratio v_{∞}/v_{orb} of the models changes. We find very similar flow morphologies to those in Paper I when comparing models with the same v_{∞}/v_{orb}. Furthermore, we confirm that the relationship found in Paper I on the morphology as a function of v_{∞}/v_{orb} also holds for different mass ratios.
Finally, in order to study how the wind velocity influences the morphology of the outflow, in Fig. 4 we compare two models with a mass ratio q = 1, a small orbital separation (a = 5 AU), and initial wind velocities of 1 km s^{−1} and 12 km s^{−1}. Since the sound speed close to the donor star is about 6 km s^{−1}, in the first model (Q1a5v1) the wind leaves the donor star in the subsonic regime, whereas for the higher velocity model (Q1a5) the wind is supersonic when it leaves the donor star. In contrast to Figs. 2 and 3, only the inner region of the simulation is shown in Fig. 4. In model Q1a5 (Fig. 4b) the wind leaving the donor remains spherically symmetric until it reaches the accretion shock and the inner spiral arm. However, when the initial wind velocity is subsonic (model Q1a5v1, Fig. 4a), the outflow in the vicinity of the donor star’s surface is clearly asymmetric and partly focused into a broad stream that is deflected behind the companion’s accretion wake. The deflection is related to the nonrotation of the donor star (see discussion in Sect. 3.2.2). We note that due to the small orbital separation in these models, the surface of the donor star and the wind launch region will be somewhat deformed by tides from the companion. This is not taken into account in our simulations, but is likely to enhance the asymmetry in the inner part of the outflow seen in Fig. 4b. On larger scales than shown in Fig. 4, both models show a spiral pattern similar to that observed in Fig. 2a. However, as expected because of the low wind velocity, the density in the spiral arms of model Q1a5v1 is higher than in model Q1a5.
3.2.2. Rotating models
In the models where the AGB donor is in corotation with the binary, we find two different regimes in the geometry of the outflow. This difference appears to be related to the initial velocity of the wind. For models Q2a5Ω and Q4a5Ω, where v_{init} = 12 km s^{−1}, the outflows are very similar to their nonrotating counterparts (Q2a5 and Q4a5; see Fig. 3). For model Q1a5Ω (not shown), we observe a similar morphology to that of Q1a5 (Fig. 2a). Both systems have an accretion disk and two spiral arms tightly wound around the binary. However, the inner spiral arm in model Q1a5 is slightly denser than the one in model Q1a5Ω.
On the other hand, when we compare the corotating models with small orbital separations (4 AU and 5 AU) and low initial velocities of the wind (v_{init} = 1 km s^{−1} and v_{init} = 5 km s^{−1}) with their nonrotating counterparts, we note clear differences in the morphologies of the outflow. At these orbital separations, the companion star is located at a position where the gas has not yet reached its terminal velocity and the donor star is very close to filling its Roche lobe (R_{d}/R_{L, 1} = 0.86 for Q2a4v1Ω and R_{d}/R_{L, 1} = 0.69 for Q2a5v5Ω). Interestingly, the geometries in these corotating models resemble the WRLOF morphology found in simulations performed in the corotating frame (e.g. de ValBorro et al. 2009; Mohamed & Podsiadlowski 2012, 2007; Chen et al. 2017). In Sect. 4.1 we compare the results of our corotating models Q2a4v1Ω and Q2a5v5Ω to the results of simulations performed in the corotating frame.
Figure 5 shows the gas density on the planes z = 0 (left) and y = 0 (right) for models Q2a4v1 (top) and Q2a4v1Ω (bottom). The velocity field in the inertial frame is plotted in the same figure. In spite of the accretion disk around the companion star which is a common feature in both models, there are many dissimilarities between the models. A clear difference that can be observed in the orbital plane, z = 0 (Figs. 5a and c), is the high gas density along the inner spiral arm for model Q2a4v1 (Fig. 5a; x ≈ −2.5 AU, 0 AU ≲ y ≲ 0.1 AU), which is substantially lower in the corotating model Q2a4v1Ω (Fig. 5c). Also, the opening angle of the inner arm is larger in the nonrotating model and wraps around the donor star at closer distances than in the corotating model. Another clear difference is the stream of gas flowing from the donor star to the companion star in model Q2a4v1Ω, which can be observed in both planes (z = 0 and y = 0). In the orbital plane, this stream is seen to flow from the inner Lagrangian point L_{1} at about x ≈ −0.49 AU to x ≈ −2.5 AU, where it encounters the accretion wake of the companion star. The presence of the stream between the two stars in Fig. 5c, which resembles RLOF, creates the impression that the effect of the rotating donor is to focus the wind in the direction of the accretor. Conversely, model Q2a4v1 does not show this stream of gas moving from L_{1} towards the companion star, but instead the gas moving from the AGB star towards the companion is deflected behind this star, where it is captured by the accretion wake and then falls onto the accretion disk. However, unlike model Q2a4v1, in model Q2a4v1Ω the gas appears to pass through the accretion wake once it reaches it. This might be due to the fact that for Q2a4v1 the density in the accretion wake is higher than for the rotating model, and is able to trap the gas more easily when it reaches it. In addition, from the flow in the y = 0 plane (Figs. 5b and d), we observe a region in model Q2a4v1Ω where material is moving away from the binary (−3 AU ≲ x ≲ −5 AU and −1 AU ≲ z ≲ 1 AU). Finally, in model Q1a4v1, the material moving away is confined to a narrower band around the z = 0 plane, whereas the gas at z ≳ 0.5 AU and z ≲ −0.5 AU is observed to be moving towards the companion star.
Fig. 5.
Density of the gas for models Q2a4v1 and Q2a4v1Ω after 9.5 orbital periods. The images on the left show the density and velocity field in the inertial frame in the orbital plane (z = 0). The images on the right show the same quantities, but on the y = 0 plane. The dashed lines correspond to the Rochelobes of the stars. 
The morphology of model Q2a5v5Ω (not shown) is somewhat similar to that of model Q2a4v1Ω with gas being focused in the direction of the companion star. We find a stream of gas flowing through the inner Lagrangian point and moving towards the inner spiral arm. As in Q1a4v1Ω, the inner spiral arm is less dense than its nonrotating analogue. The main difference we find with model Q1a4v1Ω is that the region in the y = 0 plane where the gas moves away from the binary is narrower in the z = 0 plane giving the impression that less material escapes the system.
3.3. Angularmomentum loss
The loss of angular momentum for nonconservative mass transfer in binaries can be parameterised as^{4}
The first term corresponds to the change in the orbital angular momentum J̇_{orb}, where β = −Ṁ_{a}/Ṁ_{d} is the fraction of ejected mass that is accreted by the companion star, Ṁ_{d} < 0 is the massloss rate, and η is the specific angular momentum of the material lost in units of the orbital angular momentum of the system per reduced mass J/μ = a^{2}Ω_{bin}. The second term in Eq. (4) corresponds to the angular momentum lost due to rotation of the donor star. If the donor star is assumed to be nonrotating, this term is zero.
As we did in Paper I, here we take the values of the mass and angular momentum lost from the system as those taken away by gas particles which cross a boundary of 3a from the centre of mass of the system. In order to reduce statistical fluctuations, we average the angularmomentum loss over intervals of one orbital period from the time the simulation reaches the steady state, which occurs after about four orbits (cf. Saladino et al. 2018, Sect. 3.5).
In the simulations where rotation is included, the numerical value measured for J̇ is the combination of the angular momentum removed from the orbit and the angular momentum taken from the spin of the donor star, J̇_{spin}, which we cannot disentangle. In order to apply Eq. (4) and derive the value of η, J̇_{spin} is assumed to be given by the analytical equation for isotropic mass loss from a spherical surface
where Ω_{spin} = Ω_{bin} in the case of corotation. In order to verify that the analytical value for J̇_{spin} given by Eq. (5) accurately represents the spin angularmomentum loss, we perform a numerical simulation of a single star with the properties of Table 1, rotating at Ω_{spin} = 2.76 × 10^{−8} s^{−1}. We find that the numerical value obtained for J̇_{spin} agrees with Eq. (5) within 2%. For this reason, for corotating systems in which the flow morphology is not strongly modified by rotation (see Sect. 3.3.2), we expect Eq. (4) to yield η values similar to those of their nonrotating analogues.
In the following we present the results for the angularmomentum loss of our models in terms of the parameter η. First we discuss the models in which the donor star is assumed to be nonrotating (J̇_{spin} = 0) and in the second part we present the results for the angularmomentum loss of the corotating models.
3.3.1. Nonrotating models
Figure 6 shows the specific angular momentum lost η in units of J/μ as a function of the velocity ratio v_{∞}/v_{orb} for the different nonrotating models studied in this paper. Colours represent equivalent mass ratios. This figure only shows those systems with nonrotating donor stars. When compared to the results of Paper I (open stars in Fig. 6), where angularmomentum loss was studied only for q = 2, we observe a similar trend of the η values as a function of v_{∞}/v_{orb}. We note that in Paper I the parameter we varied was v_{∞}, whereas here a is varied. For high velocity ratios, the values of η approach the isotropic wind^{5} values given by
which are shown as dashed lines in Fig. 6. On the other hand, when v_{∞}/v_{orb} decreases, the angularmomentum loss of the models in enhanced. A notable difference with Paper I is that the η values found in this work for low v_{init} are lower than those found in Paper I at the same v_{∞}/v_{orb}. However, it should be kept in mind that in this work there are a number of differences in the numerical method compared to Paper I. The stellar radius is larger and the effective temperature of the donor star is higher in this work. This leads to a different temperature profile in the outflow, which may influence the interaction between the gas and the stars. In addition, the mechanism of injection of the wind particles is different from the one used in Paper I. This may affect the angular momentum transfer, in particular when the companion star is located at a position where the wind is still being accelerated.
Fig. 6.
Orbital angularmomentum loss in units of J/μ as a function of v_{∞}/v_{orb} for models assumed to be nonrotating. Colours represent equivalent mass ratios. Dashed lines show the value of η_{iso} for different mass ratios. Filled circles correspond to models in which v_{∞} = 15.1 km s^{−1}, triangles to models where v_{∞} = 10.9 km s^{−1}, and squares to models where v_{∞} = 6 km s^{−1}. The solid lines correspond to Eq. (7) for different mass ratios. 
From Fig. 6, we also observe for similar values of v_{∞}/v_{orb} that as the mass ratio increases, the angularmomentum loss decreases. This behaviour is expected according to the results discussed in Sect. 3.2.1. Given the lower densities behind the accretion wake for higher mass ratios, the torque exerted by the wake on the binary is weaker, resulting in a smaller exchange of angular momentum. An interesting feature we find is that for low v_{∞}/v_{orb} and q ≲ 2, the value of η appears to level off at a maximum of about 0.6. However, when q ≳ 3, η never reaches such high values. Finally, although for most of our simulations η seems to grow monotonically with decreasing v_{∞}/v_{orb} until reaching its maximum for the lowest v_{∞}/v_{orb}, this does not seem to be the case for the q = 4/3 models. For these models the maximum value of η occurs for v_{∞}/v_{orb} ≈ 0.6 and it decreases for the lowest v_{∞}/v_{orb}.
In order to apply our results in binary population synthesis simulations, the angularmomentum loss needs to be expressed as a function of the physical parameters of the systems. In Paper I, we provide an analytic relation for the angularmomentum loss as a function of v_{∞}/v_{orb}, but independent of q. However, we find that this relation only matches our current results for high v_{∞}/v_{orb} and mass ratios between 1 and 2. For this reason, we provide a new function that describes our results as a function of the mass ratio and the ratio between the terminal velocity of the wind and the system’s orbital velocity
where
and η_{iso} is given by Eq. (6).
This fit is overplotted in Fig. 6 for several values of the mass ratio, corresponding to our simulations. For comparison, we also plot the results for the angularmomentum loss from Jahanara et al. (2005, blue crosses), where the wind is accelerated in a manner similar to that used in this work, and the results of Chen et al. (2018, orange plus signs), who used a different method to model the stellar wind. For both works, we apply Eq. (4) to the published results to correct for the spin angular momentum due to corotation of the system.
3.3.2. Rotating models
Figure 7 shows the orbital angularmomentum loss, as measured by the parameter η, as a function of the ratio v_{∞}/v_{orb} for the corotating models (open shapes) with their nonrotating counterparts (filled shapes). We observe that for v_{∞}/v_{orb} ≳ 0.8 the orbital angularmomentum loss is equal for the rotating and nonrotating models. This means that the extra angularmomentum loss measured in the corotating models indeed comes from the spin of the donor star and is correctly described by Eq. (5). However, for low v_{∞}/v_{orb} the orbital angularmomentum loss in the corotating models is lower than that measured for the nonrotating models. The difference is small for model Q1a5Ω, and somewhat larger for Q2a4v1Ω (≈8%) and Q2a5v5Ω (≈13%). We note that η measures the exchange of angular momentum between the escaping gas and the binary orbit. The fact that η is independent of rotation for higher v_{∞}/v_{orb} is consistent with the fact that we find no differences in the morphology between rotating and nonrotating models at these velocities (see Sect. 3.2.2). The lower angular momentum exchange in models Q2a4v1Ω and Q2a5v5Ω compared to Q2a4v1 and Q2a5v5, respectively, could be due to a lower torque exerted by the accretion wake on the binary in the rotating models, given the lower density in the inner spiral arms for these models compared to their nonrotating analogues.
Fig. 7.
Orbital angularmomentum loss in terms of η for models in corotation with the binary and their nonrotating analogues. Colours represent equivalent mass ratios. Circles correspond to models where the initial velocity was set as v_{init} = 12 km s^{−1} (v_{∞} = 15.1 km s^{−1}) and squares to models where v_{init} = 1 km s^{−1} (v_{∞} = 6 km s^{−1}). Open shapes are the corresponding corotating models. 
3.4. Massaccretion efficiency
The BHL approximation provides an estimate of the massaccretion efficiency onto a body moving at given speed in a highvelocity wind assuming supersonic flow. Under this approximation, in the case of a binary system in which one of the stars is losing mass via winds, the companion star will accrete a fraction of the mass lost by the donor given by
where α_{BHL} is the efficiency parameter, which has a value between 0.5 and 1, and is the relative wind velocity seen by the accretor (Theuns et al. 1996).
Considering the assumptions made in the canonical BHL approximation and that AGB winds are very slow, a more suitable way to estimate the massaccretion efficiencies is to measure them directly from detailed hydrodynamical simulations. However, we note that since we cannot resolve the surface of the companion star computationally, our numerical results for the massaccretion efficiency should be taken as upper limits. As discussed in Paper I, rather than measuring the mass flux that crosses the boundary of the sink particle, a more reliable measurement for the accretion rate onto the companion star is derived from the sum of the mass accretion rate on the sink particle and the rate at which mass is added to the accretion disk. For this reason, and following the method used in Paper I, we adopt as the massaccretion efficiency the net flux crossing a shell centred on the companion star with radius equal to 0.4 R_{L, 2}. As we did for the angularmomentum loss, to avoid statistical fluctuations, we measure the average massaccretion efficiency over several orbital periods once the system has reached the steady state.
3.4.1. Nonrotating models
Figure 8 shows the fraction of mass accreted β = Ṁ_{0.4RL,2}/Ṁ_{d} as a function of v_{∞}/v_{orb}. The dotted lines correspond to the predicted BHL value (Eq. (8)) with α_{BHL} = 0.75 (as Abate et al. 2013). As in Paper I, for individual mass ratios we find that the lower the v_{∞}/v_{orb} ratio, the higher the measured massaccretion efficiency. In general, we find that for different mass ratios and v_{∞}/v_{orb} > 0.5, the accretion efficiency measured in our models is always higher than the BHL prediction (see Table 3). The maximum difference between the BHL approximation and our results occurs for models Q1a5 and Q13a5v5, where we find that the value of β from the hydrodynamical models is about 1.6 times higher than the BHL prediction. Only when v_{∞}/v_{orb} is very high does the accretion efficiency approach the BHL approximation. In the cases where this occurs (see Sect. 3.2.1) the morphologies of the outflow are less complex, suggesting less interaction between the wind and the companion star. On the other hand, for models with v_{init} = 1 km s^{−1} there is no clear trend in the accretion efficiencies we measure, which can be either higher or lower than the β values predicted by the BHL approximation. The maximum accretion efficiencies we obtain never exceed about 30%, in contrast with the accretion efficiencies in Paper I, where for low velocity ratios, we found values up to about 40%. We should note that the stellar parameters used in this study differ from Paper I, and that the current numerical resolution is much lower than the resolution used in Paper I, which could be influencing our results (see Sect. 4.3).
Fig. 8.
Massaccretion efficiency for the nonrotating models. Colours represent equivalent mass ratios. Circles correspond to models where v_{∞} = 15.1 km s^{−1}, triangles to models where v_{∞} = 10.9 km s^{−1}, and squares to models where v_{∞} = 6 km s^{−1}. The dotted lines correspond to the values predicted by BHL as given by Eq. (10) assuming α_{BHL} = 0.75. Continuous lines correspond to the fit given by Eq. (11). 
Results from the simulations.
In order to implement our results in a binary population synthesis code, we construct a function that describes the behaviour found from our hydrodynamical simulations for the massaccretion efficiency. As we did for the angularmomentum loss, we fit our results as a function of mass ratio and the ratio v_{∞}/v_{orb}. A relation that describes our numerical results is given by
where
and

–
k_{1} = 1.7 + 0.3q,

–
k_{2} = 0.5 + 0.2q,

–
β_{max} = min(0.3, 1.4q^{−2}),
and β_{BHL} is given by Eq. (10) with α_{BHL} = 1. This relation is shown in Fig. 8 for different mass ratios.
3.4.2. Rotating models
Figure 9 shows the mass accretion efficiencies for the models assumed to be in corotation. We find that the massaccretion efficiency for models Q2a5Ω and Q4a5Ω is slightly higher (by a factor of about 1.1) than for their nonrotating counterparts, while for Q1a5Ω and Q2a5v5Ω it is slightly lower. However, for model Q2a4v1Ω, where the injection velocity of the wind is v_{init} = 1 km s^{−1} and the geometry of the outflow resembles WRLOF, we find that the accretion efficiency is less than half the value of its nonrotating analogue (see Table 3). Although we may attribute this low massaccretion efficiency to material crossing the lowdensity accretion wake of the companion star (see Sect. 3.2.2), this appears to be in contrast to the SPH simulations performed by Mohamed & Podsiadlowski (2007) in which they find high accretion efficiencies when mass transfer occurs via WRLOF. For comparison, in Fig. 9 we plot the mass accretion efficiencies from Abate et al. (2013) based on the hydrodynamical simulations of Mohamed (2010) as a function of v_{∞}/v_{orb} (brown Y shapes). We compute v_{orb} from the orbital periods provided in Abate et al. (2013) and we take v_{∞} = 4 km s^{−1} as in Mohamed & Podsiadlowski (2012). In the same figure, we also show the massaccretion efficiencies from the hydrodynamical models of Chen et al. (2017; orange plus signs). We observe that in both works the trend of the massaccretion efficiency as a function of v_{∞}/v_{orb} is similar to the one we find, i.e. β decreases with increasing v_{∞}/v_{orb}. However, for v_{∞}/v_{orb} in the range between 0.5 and 1, both Mohamed & Podsiadlowski (2012) and Chen et al. (2017) find much higher values than predicted by BHL. This differs from the values we find in the same velocity range, where our β values are only slightly higher than BHL (see Sect. 4 for a discussion). On the other hand, for v_{∞}/v_{orb} < 0.5, the β values found by Mohamed & Podsiadlowski (2012) decrease again and, as in our results, for the lowest velocity ratio β is much lower than the BHL expected value, even though WRLOF occurs.
3.5. Change in orbital separation
Due to the loss of angular momentum via stellar winds, the orbital separation of the binary will be affected. The rate of change of a can be derived from angular momentum conservation (see Paper I):
where the parameters η and β are computed from detailed hydrodynamical simulations. We should note that when there is efficient spinorbit coupling by tidal interactions, another term should be added to Eq. (13) to take this into account.
Another method that can be used to compute the change in a is to measure it dynamically from our numerical simulations. To compute it in this way, we use the Newtonian description of the motion of two masses under their mutual gravitational field (see Appendix C). From our numerical simulations we know the position r and the relative velocity of the stars v at any time. From these we compute the orbital energy and angular momentum per reduced mass, and by approximating the orbits as Keplerian, we compute instantaneous values of a and e. We then calculate ȧ as the average rate of change of a over the last orbits. We find that e is not exactly zero, but its value shows a sinusoidal variation for all models. However, this variation is so small (∼10^{−5}) that it can be neglected.
Figure 10 shows the change in orbital separation calculated in the two different ways described above for the nonrotating models. On the xaxis we plot the values of ȧ/a measured directly from the simulations (Eqs. (C.1) and (C.3)) and on the yaxis we show the values computed by substituting the measured values of η and β (Table 3) in Eq. (13). For most systems in which the orbit expands the two values agree quite well. In these binaries, the wind leaves the system in an almost spherically symmetric way. However, in systems with shrinking orbits we find a substantial difference between the two ȧ values. The discrepancy appears to be correlated with the speed at which the orbit shrinks, and in most cases the systems shrink more slowly if measured directly from the simulations. We note that the cases where the discrepancy is large correspond to models where the interaction between the gas and the stars is strong (squares and triangles in Fig. 10), i.e. where more angular momentum is transferred from the binary orbit to the gas. This means that if angular momentum is not conserved precisely in these cases, an error in this quantity may build up.
Fig. 10.
Evolution of the orbital separation ȧ/a, measured with two different methods. The values on the xaxis are measured directly from the evolution of the orbits of the stars. The values on the yaxis are obtained by substituting the parameters η and β measured from the gas dynamics in Eq. (13). 
For a TreeSPH code like FI angular momentum is not conserved exactly (Wijnen et al. 2016). For this reason, we compute how well angular momentum is conserved for our models. In order to do so, we compute the total initial angular momentum, J_{tot, init} = J_{stars, init} + J_{gas, init}, where J_{stars, init} is the initial orbital angular momentum of the binary and J_{gas, init} is the angular momentum of the gas which is initially in the simulation. We also compute the total angular momentum at the moment we stop the simulations, J_{tot, end} = J_{stars, end} + J_{gas, end} + J_{acc} + J_{esc}, where J_{stars, end} is the orbital angular momentum at the last timestep of our simulation, J_{gas, end} is the angular momentum of the gas at the end of the simulation, J_{acc} is the cumulative angular momentum of the gas accreted by the sink particle, and J_{esc} the cumulative angular momentum of the gas we removed during the simulation. We express the error δJ as the difference J_{tot, end} − J_{tot, init} in units of the total initial angular momentum per unit of time. We find that in most of our models angular momentum is conserved fairly well. For model Q3a20, where the difference between the two methods used to compute ȧ/a is only 4%, we find that δJ = 4.46 × 10^{−7}J_{tot, init} yr^{−1}. However, for model Q2a4v1 where the discrepancy in ȧ/a is ≈37%, angular momentum is conserved up to 2.47 × 10^{−5}J_{tot, init} per year. Since the error in the angular momentum of the stars and the gas is a fraction of the total error computed and we cannot separate the errors, it is not possible to precisely correct the ȧ/a values computed either dynamically or from the parameters β and η. However, we have verified by error propagation that both (ȧ/a)_{bin} and (ȧ/a)_{dyn} are consistent within the total error in the angular momentum budget for each system.
4. Discussion
4.1. Morphology of the outflow: rotating vs nonrotating systems
Observations of lowmass binary systems interacting via winds show spiral patterns wrapped around the stars (e.g. Kim et al. 2017; Ramstedt et al. 2017). Numerical simulations also produce these features in the outflow (e.g. Mastrodemos & Morris 1998; Liu et al. 2017; Saladino et al. 2018). In addition, some planetary nebulae, which are thought to be the final evolutionary stage of low and intermediatemass stars, show morphologies that cannot be explained by singlestar evolution. For these reasons, as discussed in Paper I, understanding how the morphologies are influenced by the interacting binary can provide an insight into the mechanisms driving mass transfer between the stars.
In this work we find that the spin of the star plays a significant role in shaping the morphology of the outflow. For cases in which the angular velocity of the donor star is zero, we find similar morphologies to those in Paper I: two spiral arms surrounding the binary system and in some cases an accretion disk. As in Paper I, a strong interaction of the wind with the stars can be observed as an increasingly complex morphology in the outflow, i.e. the morphology of the wind differs significantly from the spherically symmetric case. In addition, we find that the strength of the interaction depends on the orbital separation and mass ratio of the binary system. With increasing orbital separation the morphology of the outflow becomes less complex, and for large orbital separations the wind escapes the binary in a nearly isotropic way. Varying the orbital separation of the system results in geometries as a function of v_{∞}/v_{orb} that are similar to those when v_{∞} is varied for a fixed separation (as was done in Paper I). In terms of the mass ratio, we find spiral patterns that become less prominent as the mass ratio is further from unity. Similar structures whose complexity decreases as a function of the mass ratio have also been found by Liu et al. (2017), who use a different equation of state (EoS). These results on the morphology can help to constrain the orbital separations and mass ratios of binary observations.
As in the zerospin models, in the rotating models we find that the geometry of the outflow becomes less complex with increasing v_{∞}/v_{orb}, and more closely resembles the spherically symmetric wind case. A similar result is observed in the hydrodynamical models in the corotating frame of Chen et al. (2017), where a spiral outflow morphology is found for wide binaries, but it becomes less prominent with increasing orbital separation (i.e. with increasing v_{∞}/v_{orb}). We find that the effect of rotation on the morphology of the wind is negligible when v_{∞}/v_{orb} ≳ 0.7 and that for these models the geometries of the outflow are similar to those observed when the star is not rotating. However, when the wind velocity is low (v_{init} = 1 km s^{−1} and v_{init} = 5 km s^{−1}) and the system is relatively close, we find structures in the outflow that resemble WRLOF. The morphologies we find are similar to those observed by de ValBorro et al. (2009), Mohamed & Podsiadlowski (2012), and Chen et al. (2017) in the corotating frame. By comparing our results with these works, we note that these geometries are generally found when the velocity of the wind is much lower than the orbital velocity of the binary. For instance, in their more recent hydrodynamical models of binary stars interacting via winds, de ValBorro et al. (2017) do not find the prominent stream flowing between the stars discussed in de ValBorro et al. (2009), where the adopted wind velocities were lower.
In the physical interpretation of Mohamed & Podsiadlowski (2007), WRLOF is expected to occur when the dust formation region is larger than the Roche lobe of the star, confining material to the gravitational potential of the donor star and filling the Roche lobe. Material is then transferred to the companion star via the inner Lagrangian point. Our models do not include dust formation and the acceleration of the wind simply balances the gravity of the star. However, given that some acceleration via gas pressure still occurs, the velocity of the wind within the Roche lobe is very low (see Table 3) allowing the wind material to fill the Roche lobe and part of the material to be transferred in a way that resembles RLOF. As can be seen in the models by de ValBorro et al. (2009) and Mohamed & Podsiadlowski (2012), we observe that the stream of gas is not exactly focused towards the companion star, but slightly deflected towards the accretion wake. Morphologies like those described here are also observed. For instance, recent observations by Bujarrabal et al. (2018) of the symbiotic binary R Aqr show an arc of gas joining the two stars, which is attributed to an episode of mass transfer taking place in the system.
When we compare our models on spatial scales larger than about 5 AU (as shown in Fig. 5) to Fig. 4 in Chen et al. (2017), we observe that generally the outflow geometries are very similar, with two differences. In the inner region of the binary system our models show the flow to be very steady compared to theirs. A likely explanation for this difference is that in their models pulsations of the AGB star are included, which make this region more dynamic. The second difference we note is that in some of their models, especially where the morphologies are similar to WRLOF, a circumbinary disk is found. In our simulations we do not see such a structure. However, we should keep in mind that in our models we remove the escaping material at a relatively short distance from the binary and that we do not include dust formation or radiative transfer, which may be relevant for the formation of this disk. Longer numerical simulations that include this physics would be needed to check whether a circumbinary disk is formed.
Finally, we note that the peculiarities in the geometries of the outflow discussed here are likely to have an impact on the evolution of the binary via differences in the massaccretion efficiency and angularmomentum loss. For instance, a denser accretion wake exerts a higher torque on the binary resulting in a larger exchange of angular momentum between the orbit of the stars and the gas.
4.2. Angularmomentum loss and the impact on the orbit
Recent hydrodynamical studies of angularmomentum loss during wind mass transfer have explored its dependence on only a few parameters. Jahanara et al. (2005) performed gridbased hydrodynamical models of binary systems undergoing mass loss via radiatively driven winds, i.e. as in this work the acceleration of the wind only balances the gravitational force of the donor star. In their models, the EoS used is adiabatic, the mass ratio is q = 1, and they varied the injection velocity of the wind. Chen et al. (2018) also performed gridbased hydrodynamical simulations that include pulsations of the AGB donor star, dust formation, cooling of the gas, and radiative transfer. Except for one model, the mass ratio (q = 2) was fixed and they studied the dependence of the angularmomentum loss as a function of the orbital separation of the binary. In Paper I we performed a few exploratory simulations which did not include acceleration of the wind, but included cooling of the gas. The works of Jahanara et al. (2005), Chen et al. (2018), and Paper I show a very similar dependence of the angularmomentum loss on v_{∞}/v_{orb}. In this work, we explore a larger grid of binary parameters and wind velocities, and we include the possibility of corotation of the donor star with the binary. In the following we extend the comparison to our present work.
4.2.1. Nonrotating models
Figure 6 shows the orbital angularmomentum loss, expressed in terms of η, as a function of v_{∞}/v_{orb} for the nonrotating models. When we compare this to the results of Jahanara et al. (2005, their radiatively driven wind case), Chen et al. (2018), and Paper I, we see that all models show a good agreement for v_{∞}/v_{orb} > 1, where η tends towards the isotropicwind value (see Eq. (6)). For the cases where v_{∞}/v_{orb} < 1, there is a clear trend of η increasing with decreasing v_{∞}/v_{orb}. This behaviour is expected because the velocity of the wind is much lower than the relative orbital velocity of the stars, thus the wind has a stronger interaction with the companion star, which allows it to draw more angular momentum from the orbit. However, we find that even though it follows the same trend, there is a slight discrepancy between Paper I and this work in the η values found as a function of v_{∞}/v_{orb}. We associate this difference with the different assumptions made for the acceleration of the wind, and with the different stellar parameters of the donor star. The different acceleration mechanism can lead to close binaries interacting in a region where the wind has not reached the terminal velocity. The difference in stellar masses and metallicity can influence the temperature profile of the gas, since in this work the effective temperature of the donor star is higher. These two effects influence the interaction between the wind and the binary. Furthermore, unlike what we anticipated in Paper I, our current models show that the values of η converge to a maximum of η ≈ 0.6 for low v_{∞}/v_{orb} and q ≤ 2.
4.2.2. Rotating models
For models in which rotation is included, we find two regimes for the measured η values. For models with v_{∞}/v_{orb} ≳ 0.7, where the spin of the star does not affect the geometry of the gas (see Sect. 3.2.2), we find that the orbital angular momentum carried away by the escaping gas is independent of the spin of the donor star, i.e. the values we obtain for η as derived from Eq. (4) are the same as in their nonrotating counterparts. This implies that the terms J̇_{orb} and J̇_{spin} act independently and that the orbital evolution of the system can be predicted by considering separately the change in angular momentum due to mass lost from the system and the change in angular momentum due to tides, as usually assumed in binary population synthesis codes. However, for v_{∞}/v_{orb} ≲ 0.7, the spin of the star modifies the outflow which also changes the interaction between the gas and the stars compared to the nonrotating models. Consequently, we find different η values for the corotating and nonrotating models. This implies that for cases in which the outflow is strongly modified by the spin of the star, detailed simulations including rotation of the donor are needed to correctly predict the amount of orbital angularmomentum loss and therefore the evolution of the orbit.
In the numerical models performed by Chen et al. (2018) in the corotating frame the parameter that is varied is the orbital separation of the system. They explain why synchronisation of the donor star with the binary is expected only for some of the orbital separations they explore. In order to investigate the cases in which the donor star is not rotating, they provide a prescription where they simply subtract the contribution of the spin of the star from the total angularmomentum loss they measure (similar to the method we use in this work). However, because of the results discussed above, the prescription they use to study the zerospin cases for the AGB star can only be applied when the rotation of the donor does not modify the outflow of the gas. In addition, due to the nature of their numerical models, they find that gas gains angular momentum as it moves radially away from the masslosing star, which they attribute to numerical errors due to the intrinsic viscosity of Eulerian codes. This is not observed in our numerical models since SPH codes are better at conserving angular momentum than gridbased codes.
As mentioned in Sect. 3.5, from Eq. (4) we can predict the rate of change of the orbital separation of the binary system due to orbital angularmomentum loss. However, if tidal interactions are effective in the binary, additional terms need to be taken into account in that equation. As shown in Appendix D, for orbital separations up to ≈10 AU, tidal friction is so strong that spinorbit coupling is expected at some point during the AGB phase of the donor star. Furthermore, even when tidal interactions are not strong enough for synchronisation they can effectively spin up the donor star (see Fig. D.1b) for orbital separations between 10−20 AU. In order to keep the donor star in corotation with the binary via tidal effects, angular momentum needs to be taken from the orbit of the binary. This means that for some of the orbital separations considered in this work tidal evolution plays an important role, implying that additional shrinking of the orbit will take place.
4.3. Massaccretion efficiency
Observations of the progeny of AGB binary systems show enhanced chemical abundances in AGB nucleosynthesis elements, which sometimes cannot be explained by the standard BHL formalism (e.g. the chemical abundances of observed CEMPs stars; Abate et al. 2015a,c). Since the assumptions made for the flow in the BHL formalism usually do not hold for the winds of AGB stars, hydrodynamical models are needed to obtain better estimates of the massaccretion efficiency of the companion star. Several previous hydrodynamical studies have focused on computing this quantity (e.g. Theuns et al. 1996; de ValBorro et al. 2009, 2017; Mohamed & Podsiadlowski 2007; Liu et al. 2017). However, from the numerical point of view this is not an easy task. Since the radius of the companion star cannot be resolved numerically, these works have used a variety of methods to model the accretion process on the surface of the star. In addition, some works have shown that the massaccretion efficiency depends on the EoS used and on the numerical resolution of the simulation (Theuns et al. 1996; Saladino et al. 2018). For these reasons comparing results from different studies is not straightforward. However, given some similarities in the trends of massaccretion efficiency as a function of mass ratio and v_{∞}/v_{orb}, we can attempt to do so.
Figure 8 shows that in agreement with Paper I for the nonrotating models, the mass accretion efficiency can be described as a function of v_{∞}/v_{orb}, where for high v_{∞}/v_{orb} values β approaches the BHL approximation. As do Liu et al. (2017), we find that β depends on the mass ratio, with higher accretion efficiencies for binary stars with more comparable masses. Liu et al. (2017) provide an analytical expression that describes their numerical results for the massaccretion efficiency as a function of the mass ratio. Since in their models a constant orbital separation of 3 AU is adopted and the stellar parameters used were different than in our models, we cannot compare our β results directly with their analytical expression. However, even for our models with the closest separation (a = 5 AU), we find β values that are much higher than predicted by their analytical expression. In their models they use an adiabatic EoS to model the thermodynamics of the gas. As shown by Theuns et al. (1996) and Saladino et al. (2018) for their models with q = 2 and a = 3 AU, an adiabatic EoS usually leads to lower accretion efficiencies than predicted by BHL. For these reasons, the fit we provide in Eq. (11) is more appropriate for binaries interacting via winds.
We can compare the β values of our corotating models with Mohamed & Podsiadlowski (2012) and Chen et al. (2017), who made similar assumptions for the spin of the donor star. Similar to the nonrotating models, these models show a dependence on the ratio v_{∞}/v_{orb}. However, in our models β is only a factor of ∼1.6 higher than in the standard BHL model, whereas in both Mohamed & Podsiadlowski (2012) and Chen et al. (2017) the difference is much larger (a factor of 2 or more). Although each of these works uses a different numerical method to measure the accretion rate onto the companion star, this is unlikely to explain the differences between the results. For instance, the accretion region in Chen et al. (2017) is smaller than the sink particle assumed in this work, which should result in a lower accretion rate. Mohamed (2010) finds that the smooth method used in Mohamed & Podsiadlowski (2012) to model the accretion process produces lower accretion efficiencies than immediately removing particles from the simulation after they enter the sink region, as we do. On the other hand, we note that the SPH resolution chosen in our models is quite low in order to make the simulations computationally efficient. As shown in Paper I, a resolution like the one adopted in our models can underestimate the accretion efficiency by a small factor compared to a higher resolution. It seems more likely that the different accretion efficiency found in our work compared to Mohamed & Podsiadlowski (2012) and Chen et al. (2017) originates from the differences in the assumed input physics. Both Mohamed & Podsiadlowski (2012) and Chen et al. (2017) include pulsations of the AGB star which make the outflow very dynamic in the vicinity of the stars compared to the steady outflow in our simulations. Moreover, in their works the acceleration of the wind is driven by the pulsations of the AGB donor star and radiation pressure on dust grains, which leads to a different wind velocity profile than in Fig. 1. These physical mechanisms may influence the amount of material that is able to reach the accretion region in their models.
In our two models in which the outflow of the gas resembles WRLOF we find that the mass accretion efficiencies are relatively low, either close to the BHL value (model Q2a5v5Ω) or below it (model Q2a4v1Ω; see Table 3). From Fig. 9, we observe that for v_{∞}/v_{orb} < 0.3 the results of the WRLOF simulations of Mohamed & Podsiadlowski (2012), as reported by Abate et al. (2013), also show a much lower value than the one expected from the BHL approximation. This may have important consequences for the results based on this model, which is used to study the population of CEMPs stars (Abate et al. 2015b, 2018). From our understanding, not only are the mass ratio and the orbital separation important parameters that influence binary interaction, but as pointed out in Paper I and confirmed in the present study, the velocity of the wind plays a major role as well. Therefore, the fit for the massaccretion efficiency by Abate et al. (2013) may only be valid for the low wind velocities studied by Mohamed & Podsiadlowski (2012), v_{∞} ≈ 4 km s^{−1}, rather than the 15 km s^{−1} they used in their population synthesis study. This implies that Abate et al. (2013) overestimate the massaccretion efficiency, which may change the number of CEMP stars they find, as well as their final orbitalperiod distribution.
4.4. Possible implications for lowmass binary evolution
In order to quantify the impact of our results for the massaccretion efficiency and angularmomentum loss on the evolution of low and intermediatemass binary stars, the relations for these quantities given in Sects. 3.3.1 and 3.4.1 need to be implemented in binary population synthesis codes. In the following we discuss how our results may impact the evolution of these binary systems.
Abate et al. (2018) apply the fit for angularmomentum loss as a function of v_{∞}/v_{orb} that we provide in Paper I in combination with the analytical formula of Abate et al. (2013) for the WRLOF massaccretion efficiency to a population of progenitors of CEMP stars. They find that this model is not able to reproduce the observed orbital periods of CEMP binary stars. A selfconsistent model where both the angularmomentum loss and massaccretion efficiency are derived under the same conditions may help to improve their results. However, we can already foresee some problems that may be encountered if we were to apply our current results. With the canonical BHL model, population synthesis codes usually underestimate the observed frequency of CEMP stars among very metal poor stars. Abate et al. (2013) showed that when mass transfer occurs via WRLOF over a wide range of separations, the high associated accretion efficiencies increase the frequency of CEMP systems. However, as discussed in Sect. 3.4, our results yield massaccretion efficiencies that are at most a factor of ≈1.6 higher than the BHL estimates and much lower than in the WRLOF model of Abate et al. (2013). As discussed by Abate et al. (2013, 2018), the main effect of a lower accretion efficiency is that CEMP stars are produced in a narrower range of orbital periods and masses. Consequently, the resulting CEMP fraction will be lower with our model than in the WRLOF prescription and will likely underestimate the observations. Another aspect that may influence the results of current binary population synthesis studies is that these typically do not include a velocity profile of the wind, but assume a constant wind velocity. As is shown in Paper I and in this work, the velocity of the wind plays a major role in driving the interaction of the binary system.
We find that the angularmomentum loss is relatively enhanced above the isotropicwind value for systems with lower mass ratios (i.e. more equalmass components). With the angularmomentum loss derived from our hydrodynamical models, systems with q = 1 and v_{∞}/v_{orb} ≲ 1.4 will shrink instead of widening, as would be predicted if the wind were spherically symmetric. However, as q increases, the maximum value of v_{∞}/v_{orb} for which the orbit shrinks decreases. For instance, for q = 4, only systems with v_{∞}/v_{orb} ≲ 0.82 will shrink. We find that 15 out of our 26 models will shrink, whereas the isotropic wind mode usually adopted in binary population synthesis codes would have predicted a widening of their orbit. Compared to the isotropicwind mode, these systems will evolve towards shorter orbital periods, implying that equalmass systems are likely to undergo a common envelope (CE) phase (see e.g. Ivanova et al. 2013, for a review on CE) for a wider range of initial orbital separations. In a test of our fits for the angularmomentum loss and massaccretion efficiencies in the binary population synthesis code binary_c we find that the maximum initial orbital separation at which systems will enter a common envelope increases by a factor between 1.25 and 1.5 AU compared to the isotropic wind case, depending on the mass ratio of the binary. The smallest increase from 4.38 AU to 5.46 AU occurs for q = 4, and the largest shift, from 4.32 AU to 6.58 AU for q = 4/3. This effect may be important in determining the formation rate of the progeny of low and intermediatemass binary systems, such as type Ia supernovae, cataclysmic binaries, and double white dwarfs.
5. Conclusions
We find that the two main parameters that determine how mass transfer proceeds in lowmass binary stars interacting via AGB winds are the mass ratio of the binary system and the ratio v_{∞}/v_{orb}. The morphology of the outflow depends on these parameters and determines the amount of mass accreted by the companion of the AGB star and the specific angular momentum that is carried away by the ejected material. Furthermore, because v_{orb} depends on the orbital separation, we find that modifying the orbital separation of the system produces results similar to those in Paper I where the orbital separation was fixed and the wind velocity was varied.
For any mass ratio, a low value of v_{∞}/v_{orb} produces high massaccretion efficiencies and large angularmomentum loss, whereas for high v_{∞}/v_{orb} the angularmomentum loss approaches the isotropicwind value and the massaccretion efficiency approaches the BHL approximation. When the masses of the stars are similar, the interaction between the wind and the companion star is stronger resulting in a larger angularmomentum loss and high accretion efficiency.
Using binary_c we model the evolution of binary systems for a large grid of orbital separations that span our hydrodynamical models. We find that for binary stars with initial orbital separations up to 7−10 AU, depending on the mass ratio, tidal forces are efficient and spinorbit coupling is expected at some point during the AGB phase. In particular, for systems with stellar parameters similar to those in our hydrodynamical models and orbital separations smaller than ≈4−6 AU the donor star is nearly in corotation with the orbit. Our hydrodynamical simulations show that for systems with v_{∞}/v_{orb} ≳ 0.7 rotation of the donor star has little effect on the morphology of the outflow, the massaccretion efficiency, and the orbital angularmomentum loss. On the other hand, for lower velocity ratios, corotation of the AGB donor star plays an important role in the way the stars interact. The outflow morphologies for the two corotating models in which v_{∞}/v_{orb} ≲ 0.7 resemble the WRLOF geometry and their orbital angularmomentum loss as well as the massaccretion efficiencies are different from their nonrotating counterparts. This implies that for v_{∞}/v_{orb} ≲ 0.7, the orbital evolution of the system cannot be predicted by treating independently the tidal interactions and the angularmomentum lost due to material escaping the binary system, as is usually done in binary population synthesis codes.
Finally, we fit an analytical relation to our results for the angularmomentum loss and massaccretion efficiency as a function of the mass ratio and v_{∞}/v_{orb}. In order to test whether our results help to reproduce the orbitalperiod distribution of the progeny of low and intermediatemass binaries interacting via AGB winds, our fits for angularmomentum loss and massaccretion efficiency need to be applied in binary population synthesis studies.
Equation (4) ignores the rotational velocity of the secondary star, which could accrete angular momentum from the wind and be spun up. However, since we model the accretor star as a sink we cannot account for the real angular momentum that is transferred (see Paper I for a discussion).
Acknowledgments
The authors thank the anonymous referee for the useful comments that improved this paper. MIS thanks Frank Verbunt and Richard Stancliffe for the very useful comments during the writing of this manuscript. MIS also thanks GlennMichael Oomen for the very interesting discussions during the development of this project. MIS is very grateful to the Stellar Dynamics and Computational Astrophysics group at Leiden, especially to Simon PortegiesZwart and Inti Pelupessy for the the very helpful input during the progress of this project. This work was supported by the Netherlands School for Research in Astronomy (NOVA). CA is the recipient of an Alexander von Humboldt Fellowship.
References
 Abate, C., Pols, O. R., Izzard, R. G., Mohamed, S. S., & de Mink, S. E. 2013, A&A, 552, A26 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Abate, C., Pols, O. R., Karakas, A. I., & Izzard, R. G. 2015a, A&A, 576, A118 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Abate, C., Pols, O. R., Stancliffe, R. J., et al. 2015b, A&A, 581, A62 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Abate, C., Pols, O. R., Izzard, R. G., & Karakas, A. I. 2015c, A&A, 581, A22 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Abate, C., Pols, O. R., & Stancliffe, R. J. 2018, A&A, 620, A63 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Asplund, M., Grevesse, N., Sauval, A. J., & Scott, P. 2009, ARA&A, 47, 481 [NASA ADS] [CrossRef] [Google Scholar]
 Beers, T. C., & Christlieb, N. 2005, ARA&A, 43, 531 [NASA ADS] [CrossRef] [Google Scholar]
 Bidelman, W. P., & Keenan, P. C. 1951, ApJ, 114, 473 [NASA ADS] [CrossRef] [Google Scholar]
 Blind, N., Boffin, H. M. J., Berger, J.P., et al. 2011, A&A, 536, A55 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Bondi, H., & Hoyle, F. 1944, MNRAS, 104, 273 [NASA ADS] [CrossRef] [Google Scholar]
 Bowen, G. H. 1988, ApJ, 329, 299 [NASA ADS] [CrossRef] [Google Scholar]
 Brookshaw, L., & Tavani, M. 1993, ApJ, 410, 719 [NASA ADS] [CrossRef] [Google Scholar]
 Bujarrabal, V., Alcolea, J., Mikołajewska, J., CastroCarrizo, A., & Ramstedt, S. 2018, A&A, 616, L3 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Chandrasekhar, S. 1934, MNRAS, 94, 444 [NASA ADS] [CrossRef] [Google Scholar]
 Chen, Z., Frank, A., Blackman, E. G., Nordhaus, J., & CarrollNellenback, J. 2017, MNRAS, 468, 4465 [NASA ADS] [CrossRef] [Google Scholar]
 Chen, Z., Blackman, E. G., Nordhaus, J., Frank, A., & CarrollNellenback, J. 2018, MNRAS, 473, 747 [NASA ADS] [CrossRef] [Google Scholar]
 de ValBorro, M., Karovska, M., & Sasselov, D. 2009, ApJ, 700, 1148 [NASA ADS] [CrossRef] [Google Scholar]
 de ValBorro, M., Karovska, M., Sasselov, D. D., & Stone, J. M. 2017, MNRAS, 468, 3408 [NASA ADS] [CrossRef] [Google Scholar]
 Decin, L., Justtanont, K., De Beck, E., et al. 2010, A&A, 521, L4 [Google Scholar]
 Decin, L., Richards, A. M. S., Neufeld, D., et al. 2015, A&A, 574, A5 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Edgar, R. 2004, New Astron. Rev., 48, 843 [NASA ADS] [CrossRef] [Google Scholar]
 Fujii, M., Iwasawa, M., Funato, Y., & Makino, J. 2007, PASJ, 59, 1095 [NASA ADS] [Google Scholar]
 Gerritsen, J. P. E., & Icke, V. 1997, A&A, 325, 972 [NASA ADS] [Google Scholar]
 Hansen, T. T., Andersen, J., Nordström, B., et al. 2016, A&A, 588, A3 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Hernquist, L., & Katz, N. 1989, ApJS, 70, 419 [NASA ADS] [CrossRef] [Google Scholar]
 Höfner, S., & Olofsson, H. 2018, A&ARv, 26, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Hoyle, F., & Lyttleton, R. A. 1939, Proc. Cambridge Philos. Soc., 35, 405 [Google Scholar]
 Hurley, J. R., Tout, C. A., & Pols, O. R. 2002, MNRAS, 329, 897 [NASA ADS] [CrossRef] [Google Scholar]
 Hut, P. 1981, A&A, 99, 126 [NASA ADS] [Google Scholar]
 Ivanova, N., Justham, S., Chen, X., et al. 2013, A&ARv, 21, 59 [Google Scholar]
 Izzard, R. G., Tout, C. A., Karakas, A. I., & Pols, O. R. 2004, MNRAS, 350, 407 [NASA ADS] [CrossRef] [Google Scholar]
 Izzard, R. G., Dray, L. M., Karakas, A. I., Lugaro, M., & Tout, C. A. 2006, A&A, 460, 565 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Izzard, R. G., Glebbeek, E., Stancliffe, R. J., & Pols, O. R. 2009, A&A, 508, 1359 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Izzard, R. G., Dermine, T., & Church, R. P. 2010, A&A, 523, A10 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Izzard, R. G., Preece, H., Jofre, P., et al. 2018, MNRAS, 473, 2984 [NASA ADS] [CrossRef] [Google Scholar]
 Jahanara, B., Mitsumoto, M., Oka, K., et al. 2005, A&A, 441, 589 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Jeffers, S. V., Min, M., Waters, L. B. F. M., et al. 2014, A&A, 572, A3 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Jorissen, A., Van Eck, S., Mayor, M., & Udry, S. 1998, A&A, 332, 877 [NASA ADS] [Google Scholar]
 Jorissen, A., Van Eck, S., Van Winckel, H., et al. 2016, A&A, 586, A158 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Kenyon, S. J. 1992, IAU Symp., 151, 137 [NASA ADS] [Google Scholar]
 Kervella, P., Homan, W., Richards, A. M. S., et al. 2016, A&A, 596, A92 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Kim, H., Trejo, A., Liu, S.Y., et al. 2017, Nat. Astron., 1, 0060 [NASA ADS] [CrossRef] [Google Scholar]
 Liu, Z.W., Stancliffe, R. J., Abate, C., & Matrozis, E. 2017, ApJ, 846, 117 [NASA ADS] [CrossRef] [Google Scholar]
 Mastrodemos, N., & Morris, M. 1998, ApJ, 497, 303 [NASA ADS] [CrossRef] [Google Scholar]
 Mathieu, R. D., & Geller, A. M. 2015, Astrophys. Space Sci. Lib., 413, 29 [NASA ADS] [CrossRef] [Google Scholar]
 Merrill, P. W. 1942, ApJ, 95, 386 [NASA ADS] [CrossRef] [Google Scholar]
 Merrill, P. W. 1944, ApJ, 99, 15 [NASA ADS] [CrossRef] [Google Scholar]
 Merrill, P. W. 1948, ApJ, 107, 317 [NASA ADS] [CrossRef] [Google Scholar]
 Merrill, P. W. 1950, ApJ, 111, 484 [NASA ADS] [CrossRef] [Google Scholar]
 Mohamed, S. S. 2010, PhD Thesis, St Edmund Hall, University of Oxford [Google Scholar]
 Mohamed, S., & Podsiadlowski, P. 2007, ASP Conf. Ser., 372, 397 [Google Scholar]
 Mohamed, S., & Podsiadlowski, P. 2012, Balt. Astron., 21, 88 [Google Scholar]
 Nagae, T., Oka, K., Matsuda, T., et al. 2004, A&A, 419, 335 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Nie, J. D., Wood, P. R., & Nicholls, C. P. 2012, MNRAS, 423, 2764 [NASA ADS] [CrossRef] [Google Scholar]
 Pelupessy, F. I., van der Werf, P. P., & Icke, V. 2004, A&A, 422, 55 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Pelupessy, F. I., Jänes, J., & Portegies Zwart, S. 2012, New Astron., 17, 711 [NASA ADS] [CrossRef] [Google Scholar]
 Pelupessy, F. I., van Elteren, A., de Vries, N., et al. 2013, A&A, 557, A84 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Pols, O. R., Karakas, A. I., Lattanzio, J. C., & Tout, C. A. 2003, ASP Conf. Ser., 303, 290 [Google Scholar]
 Portegies Zwart, S., McMillan, S., Harfst, S., et al. 2009, New Astron., 14, 369 [NASA ADS] [CrossRef] [Google Scholar]
 Portegies Zwart, S., McMillan, S. L. W., van Elteren, E., Pelupessy, I., & de Vries, N. 2013, Comput. Phys. Commun., 183, 456 [Google Scholar]
 Price, D. J. 2012, ASP Conf. Ser., 453, 249 [NASA ADS] [Google Scholar]
 Ramstedt, S., Mohamed, S., Vlemmings, W. H. T., et al. 2017, A&A, 605, A126 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Rasio, F. A., Tout, C. A., Lubow, S. H., & Livio, M. 1996, ApJ, 470, 1187 [NASA ADS] [CrossRef] [Google Scholar]
 Saladino, M. I., Pols, O. R., van der Helm, E., Pelupessy, I., & Portegies Zwart, S. 2018, A&A, 618, A50 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Schure, K. M., Kosenko, D., Kaastra, J. S., Keppens, R., & Vink, J. 2009, A&A, 508, 751 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Theuns, T., Boffin, H. M. J., & Jorissen, A. 1996, MNRAS, 280, 1264 [NASA ADS] [CrossRef] [Google Scholar]
 van der Helm, E., Saladino, M. I., Portegies Zwart, S. F., & Pols, O. R. 2019, A&A, 625, A85 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Vassiliadis, E., & Wood, P. R. 1993, ApJ, 413, 641 [NASA ADS] [CrossRef] [Google Scholar]
 Wijnen, T. P. G., Pols, O. R., Pelupessy, F. I., & Portegies Zwart, S. 2016, A&A, 594, A30 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Zahn, J.P. 1977, A&A, 57, 383 [NASA ADS] [Google Scholar]
Appendix A: CEMP abundances
Table A.1 shows the surface abundances for an AGB star with an initial metallicity of Z = 10^{−4}, obtained with the method described in Sect. 2. These abundances are used in the implementation of the cooling of the gas.
Surface abundances of the AGB donor star, for initial metallicity 10^{−4}.
Appendix B: Opening angle
An object moving in a fluid causes small disturbances that propagate as sound waves with speed c. In the supersonic case (v > c) the sound waves propagate within a cone delimited by tangents to the sound wave spheres. The angle between the line delimiting the cone and the flow moving at velocity v forms what is known as the opening angle θ. In Fig. B.1 we see that from trigonometry this angle is
The Mach number ℳ is defined as the ratio v/c.
Fig. B.1.
Schematic representation of the opening angle. The green lines represent the cone formed by the sound waves (light blue) propagating at different times (Adapted from Prof. Jerry M. Seitzman lecture notes http://www.seitzman.gatech.edu/classes/ae2010/machanglenumber.pdf). 
Appendix C: Dynamical change in a
In the Newtonian approximation, the total energy of the relative orbit of two point masses moving under their mutual gravity is given by the sum of the kinetic and potential energies. The energy per reduced mass is defined as
where r is the distance between the bodies and v their relative velocity. The angular momentum of the binary per reduced mass is
Since the change in the orbit is small over one orbital period, we approximate the orbits as Keplerian. In this approximation, the relative orbit of the two masses is always a conic section. In order to check whether the orbits of our models remain circular, we assume that the orbit of the system is an ellipse for which the orbital energy per reduced mass is
and the angular momentum per reduced mass
At any time in our models, we estimate e from Eqs. (C.1), (C.4), and (C.2). The value of a can be determined at any time from the orbital energy using Eqs. (C.1) and (C.3), or if the orbit remains circular, a can be computed from the angular momentum using Eqs. (C.4), (C.3), and (C.2). We estimate ȧ as the average rate of change of a over the last five completed orbital periods.
Appendix D: Spin of the AGB star
In this section we present the results obtained for spinorbit coupling for the binary models evolved with binary_c as described in Sect. 2.3. We assume the stars to be initially nonrotating and we follow the evolution of the spin and orbital angular velocity (Ω_{spin} and Ω_{bin}) under the influence of tidal interactions.
Figure D.1a shows the ratio Ω_{spin}/Ω_{bin} as a function of the orbital separation at the moment the AGB donor star has the stellar parameters used as input for our hydrodynamical simulations (Table 1). The different coloured lines correspond to different mass ratios. We see that Ω_{spin}/Ω_{bin} never reaches a value equal to one. However, for small orbital separations (a ≈ 4−6, depending on mass ratio) Ω_{spin}/Ω_{bin} ≈ 0.95, which indicates that the AGB donor is nearly in corotation with the orbit of the binary. We also observe that for lower mass ratios, the tidal effects are strong at relatively large orbital separations. For instance, for a ≈ 10 AU and q = 1, the donor star rotates at ≈60% of the angular velocity of the binary at the time when the mass of the donor is M_{d} = 1.2 M_{⊙} (Fig. D.1a). For q = 4 at the same orbital separation, the donor star rotates only at ≈30% of the orbital angular velocity. It is important to note that even for large orbital separations, tidal interaction can trigger rotation of the AGB star. This effect drops considerably with orbital separation and becomes negligible for a ≳ 20 AU.
Fig. D.1.
Panel a: Ω_{spin}/Ω_{bin} as a function of the semimajor axis at the time the AGB star has stellar parameters as in Table 1. The coloured lines represent different mass ratios. The region where the vertical lines start corresponds to orbital separations for which the binary system had gone into a common envelope at the time we computed the binary parameters. Panel b: maximum Ω_{spin}/Ω_{bin} as a function of the initial orbital separation of the system. Colours represent equivalent mass ratios. 
Because at the time we measure Ω_{spin}/Ω_{bin}, the donor star is at the superwind phase and has already gone through previous episodes of mass loss and interaction with its less evolved companion before and during the AGB phase, we also compute the maximum angular velocity that the donor star reaches during the AGB phase. Figure D.1b shows the maximum Ω_{spin}/Ω_{bin} as a function of the initial orbital separation of the binary for different mass ratios. We see that spinorbit coupling can lead to nearcorotation at some point during the AGB phase even for relatively wide binary stars (up to ≈10 AU), and that the orbital separation at which it occurs decreases with increasing mass ratio. By the time the superwind stage is reached (Fig. D.1a) the star has spun down as a result of strong mass loss and corotation is lost again, except in the closest orbits.
Finally, we also explore the case in which the primary star is initially rotating. We arbitrarily set the initial rotational velocity of the donor star to be 200 km s^{−1}. We find that at the moment the donor star has the stellar parameters used as initial setup of our SPH simulations (see Table 1), the angular velocity of the AGB star has a similar trend to that shown in Fig. D.1a as a function of the orbital separation for a given mass ratio, up to about 20 AU. However, for orbital separations larger than 20 AU the spin of the AGB star does not tend to zero, but the star is rotating slowly at about ≈0.1Ω_{bin}. This shows that for a ≲ 20 AU the spin of the star is dominated during the AGB phase by tidal evolution and is independent of the initial spin of the donor star.
All Tables
All Figures
Fig. 1.
Velocity and density profiles of the SPH gas particles leaving a single star with v_{init} = 12 km s^{−1}. The top dashed line corresponds to a constant velocity of 15.1 km s^{−1}. The lower dashed line is the corresponding analytical density in a steady state. The red rectangle delimits the inner boundary of the SPH particles around the donor star. 

In the text 
Fig. 2.
Gas density in the orbital plane for models Q1a5, Q1a10, and Q1a15 (same mass ratio, but different orbital separations). The companion star is located on the left and the AGB donor on the right. Each figure is scaled to display the full simulation, i.e. a sphere with radius equal to 5a, and shows the situation after 9.5 orbital periods. 

In the text 
Fig. 3.
Same as Fig. 2, but for models with the same orbital separation (a = 5 AU) and different mass ratios. 

In the text 
Fig. 4.
Gas density on the orbital plane for for models with the same orbital separation and mass ratio (a = 5 AU, q = 1), but different initial wind velocities (1 km s^{−1} and 12 km s^{−1}). For clarity, the density scale here is different from the previous figures. For model Q1a5 (right) this figure is a closeup version of Fig. 3a. 

In the text 
Fig. 5.
Density of the gas for models Q2a4v1 and Q2a4v1Ω after 9.5 orbital periods. The images on the left show the density and velocity field in the inertial frame in the orbital plane (z = 0). The images on the right show the same quantities, but on the y = 0 plane. The dashed lines correspond to the Rochelobes of the stars. 

In the text 
Fig. 6.
Orbital angularmomentum loss in units of J/μ as a function of v_{∞}/v_{orb} for models assumed to be nonrotating. Colours represent equivalent mass ratios. Dashed lines show the value of η_{iso} for different mass ratios. Filled circles correspond to models in which v_{∞} = 15.1 km s^{−1}, triangles to models where v_{∞} = 10.9 km s^{−1}, and squares to models where v_{∞} = 6 km s^{−1}. The solid lines correspond to Eq. (7) for different mass ratios. 

In the text 
Fig. 7.
Orbital angularmomentum loss in terms of η for models in corotation with the binary and their nonrotating analogues. Colours represent equivalent mass ratios. Circles correspond to models where the initial velocity was set as v_{init} = 12 km s^{−1} (v_{∞} = 15.1 km s^{−1}) and squares to models where v_{init} = 1 km s^{−1} (v_{∞} = 6 km s^{−1}). Open shapes are the corresponding corotating models. 

In the text 
Fig. 8.
Massaccretion efficiency for the nonrotating models. Colours represent equivalent mass ratios. Circles correspond to models where v_{∞} = 15.1 km s^{−1}, triangles to models where v_{∞} = 10.9 km s^{−1}, and squares to models where v_{∞} = 6 km s^{−1}. The dotted lines correspond to the values predicted by BHL as given by Eq. (10) assuming α_{BHL} = 0.75. Continuous lines correspond to the fit given by Eq. (11). 

In the text 
Fig. 9.
Similar to Fig. 8, but for the rotating models. 

In the text 
Fig. 10.
Evolution of the orbital separation ȧ/a, measured with two different methods. The values on the xaxis are measured directly from the evolution of the orbits of the stars. The values on the yaxis are obtained by substituting the parameters η and β measured from the gas dynamics in Eq. (13). 

In the text 
Fig. B.1.
Schematic representation of the opening angle. The green lines represent the cone formed by the sound waves (light blue) propagating at different times (Adapted from Prof. Jerry M. Seitzman lecture notes http://www.seitzman.gatech.edu/classes/ae2010/machanglenumber.pdf). 

In the text 
Fig. D.1.
Panel a: Ω_{spin}/Ω_{bin} as a function of the semimajor axis at the time the AGB star has stellar parameters as in Table 1. The coloured lines represent different mass ratios. The region where the vertical lines start corresponds to orbital separations for which the binary system had gone into a common envelope at the time we computed the binary parameters. Panel b: maximum Ω_{spin}/Ω_{bin} as a function of the initial orbital separation of the system. Colours represent equivalent mass ratios. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.