Free Access
Issue
A&A
Volume 642, October 2020
Article Number A183
Number of page(s) 7
Section Stellar structure and evolution
DOI https://doi.org/10.1051/0004-6361/202038026
Published online 20 October 2020

© ESO 2020

1. Introduction

Neglecting tides and magnetic braking, Packet (1981) showed that mass gainers continue to rotate with the critical velocity after the accretion of only 5%–10% of their own mass. This critical velocity is in most cases much faster than observed.

Van Rensbergen & De Greve (2008) published results of binary evolutionary computations aimed at reproducing the observed orbital parameters, masses, location in the HRD and rotational velocities of a collection of well-studied Algol-type semi-detached binaries. In this binary evolution code, the effects of mass and angular momentum losses from the system by stellar wind were included, assuming that stellar wind carries the orbital angular momentum of both stars into space. The mass losses by stellar wind were from Vink et al. (2001) for stars hotter than 12 500 K, and De Jager et al. (1988) for cooler stars. In the case of liberal evolution of the system we assumed that the mass-losing gainer carries its specific orbital angular momentum into space.

Our previous work (Van Rensbergen & De Greve 2008) used tidal interaction according to Wellstein (2001). In a subsequent paper, Van Rensbergen & De Greve (2016) included a more refined treatment of tidal interactions. Hereinafter, tides for a star with a radiative atmosphere are very different to those of a star in convection (see e.g. Hilditch 2001). The convective mode was always applied during Roche Lobe Over Flow (henceforth: RLOF). Meridional circulation, as proposed by Tassoul (2000), was added as a significant contributor to the tidal action. This could not, however, explain the observed rotational velocities of the mass gainers published by Van Hamme & Wilson (1990), Miller et al. (2007), Glazunova & Yushcenko (2008) and Dervisoglu et al. (2010).

The influence of magnetic braking on the evolution of binaries has been acknowledged repeatedly in the literature (see e.g. Eggleton & Kiseleva-Eggleton 2002 and Song et al. 2018). Dervisoglu et al. (2010) showed how the rotation of gainers in binaries is braked by assumed magnetic fields of a few kG. Vanbeveren et al. (2018) determined the evolution of the spin of the O components in (WR+O) binaries in a similar way.

The effects of magnetic braking were so far not included in our code. In the present paper, we calculated the magnetic field strength, and hence the magnetic braking, using a simple model of differential rotation. The magnetic fields are variable with time and different for different stars, depending on the evolution of physical parameters at the interface between the inner core and the outer shell of the mass-gaining star. We report the results of magnetic braking included, in order to better understand the observed rotational velocities.

2. Generation of the magnetic field of the gainer

A solid body rotator with constant angular velocity Ω from the centre to the edge of the star cannot develop a magnetic field. For this paper, we used magnetic fields that are produced by the dynamo model of Spruit (2002) which is at work when the angular velocity rises with an amount ΔΩ over a distance Δr. The magnetic field is then directly proportional to the quantity q:

q = Δ Ω Ω Δ r r . $$ \begin{aligned} {q}={{\Delta \Omega }\over {\Omega }} \left. {{\Delta {r}}\over {r}}\right.. \end{aligned} $$(1)

In order to calculate q, we used a simplified model to introduce differential rotation into the gainer. We considered the gainer as a star that is composed as a core surrounded by a shell. Before the start of RLOF, both parts rotate synchronously and produce no magnetic field. Magnetic fields can develop from the start of RLOF onwards. The core has a radius Rcore and is a rigid rotator with angular velocity Ωcore. The differential rotation of the shell rises from Ωcore at Rcore, to Ωedge at Rg: the physical radius of the gainer. A reasonable evaluation of Ωedge is hereafter proposed in relation (8). Relation (1) can now be read as:

q = ( Ω edge Ω core ) Ω core ( R g R core ) R core . $$ \begin{aligned} \left.{q}={{\left({\Omega _\mathrm{edge} }-\Omega _\mathrm{core} \right)}\over {\Omega _\mathrm{core} }}~\right.~{{\left({{R}_\mathrm{g} }-{R}_\mathrm{core} \right)}\over {{R}_\mathrm{core} }}. \end{aligned} $$(2)

Apart from the evaluation of the quantity q in relation (2), the following quantities needed to compute the magnetic field with the dynamo model of Spruit (2002) are: the Brunt-Väsälä frequency N, the thermal conductivity κ, the magnetic diffusivity η and the density ρ.

These quantities are calculated in the gainer at every stage of the evolution of the binary.

Spruit (2002) calculated the radial component Br from the azimuthal component Bϕ of the magnetic field with relations (5) and (6).

The minimum shear rate required for a dynamo to operate is given by the quotient of the magnetic field strength (B) needed to produce the Alfvén frequency divided by the minimum field strength (Bc) for the Tayler instability. These quotients are shown in relations (3) and (4).

If thermal conductivity can be ignored (subscript “0”) one obtains:

B 0 B c 0 = q ( Ω core N ) 7 4 ( R core 2 N η ) 1 4 < 1 . $$ \begin{aligned} {{{B}_{0}}\over {{B}_{c0}}}~=~{q}~{\left({{\Omega _\mathrm{core} }\over {N}}\right)^{{7}\over {4}}}~\left({{{R}_\mathrm{core} ^{2}~{N}}\over {\eta }}\right)^{{1}\over {4}}~ < ~1. \end{aligned} $$(3)

If this is the case, one includes the thermal conductivity κ in the evaluation and calculates the quotient (subscript “1”):

B 1 B c 1 = q 1 2 ( Ω core N ) 7 8 ( R core 2 N κ ) 1 8 ( κ η ) 1 2 . $$ \begin{aligned} {{{B}_{1}}\over {{B}_{c1}}}~=~{q}^{{1}\over {2}}~{\left({{\Omega _\mathrm{core} }\over {N}}\right)^{{7}\over {8}}}~\left({{{R}_\mathrm{core} ^{2}~{N}}\over {\kappa }}\right)^{{1}\over {8}}~{\left({{\kappa }\over {\eta }}\right)}^{{1}\over {2}}. \end{aligned} $$(4)

If B 1 B c 1 $ {B}_{1}\over{{B}_{c1}} $ > 1, the magnetic field is calculated with relation (6).

We followed this logic. However, when the quotients in relations (3) and (4) are both larger than 1, both modes are able to trigger the working of the dynamo. In that case, the mean value of Br0 in relation (5) and Br1 in relation (6) were taken for the magnetic field that governs the magnetic braking:

B ϕ 0 = R core ( 4 π ρ 2 ) 1 2 q Ω core 2 N ; B r 0 B ϕ 0 = q ( Ω core N ) 2 $$ \begin{aligned} {B}_{{\phi }0}={{R}_\mathrm{core} ~{{(4~\pi }{\rho ^{2}})^{{1}\over {2}}}~{q}~{{\Omega _\mathrm{core} ^{2}}}\over {{N}}};~{{{B}_{r0}}\over {{B}_{{\phi }0}}}={q}~\left({{\Omega _\mathrm{core} }\over {{N}}}\right)^{2} \end{aligned} $$(5)

B ϕ 1 = R core ( 4 π ρ 2 ) 1 2 q 1 2 Ω core ( Ω core N ) 1 8 ( κ N R core 2 ) 1 8 ; B r 1 B ϕ 1 = ( Ω core N ) 1 4 ( κ N R core 2 ) 1 4 . $$ \begin{aligned} \begin{aligned} {B}_{{\phi }1}={R}_\mathrm{core} ~{(4~\pi }{\rho ^{2}})^{{1}\over {2}}~{q}^{{1}\over {2}}~ \Omega _\mathrm{core} ~{\left({{\Omega _\mathrm{core} }\over {{N}}}\right)}^{{1}\over {8}}~{\left({{\kappa }\over {{N}~{R}_{\rm core}^{2}}}\right)}^{{1}\over {8}} ; \\ {{{B}_{r1}}\over {{B}_{{\phi }1}}}={\left({{\Omega _\mathrm{core} }\over {{N}}}\right)}^{{1}\over {4}}~{\left({{\kappa }\over {{N}~{R}_\mathrm{core} ^{2}}}\right)}^{{1}\over {4}}. \end{aligned} \end{aligned} $$(6)

The values of ρ, κ, N, and η in relations (3)–(6) are evaluated at the interface between the core and the shell.

From the start of RLOF onwards, the core of the gainer is not influenced by the matter coming from the donor and impinging on the surface of the gainer. The amount of angular momentum that is added only to the shell is given by Packet (1981), reduced by the quotient of the impact-parameter d over the radius of the gainer Rg, rendering the spin-up inefficient at close approach:

Δ J spin , shell + = 6.04534 × 10 54 R g ( M g + Δ M g 2 ) 1 2 ( d R g ) . $$ \begin{aligned} \Delta {J}_\mathrm{spin,shell} ^{+}=6.04534 \times 10^{54}~{R}_\mathrm{g} ~ \left({M}_\mathrm{g} + {\Delta {M}_\mathrm{g} \over 2}\right)^{1\over 2}~\left({{d}\over {{R}_\mathrm{g} }} \right). \end{aligned} $$(7)

Δ J spin,shell + $ \Delta {J}_\mathrm{spin,shell}^{+} $ is expressed in cgs units, whereas masses and radii are, respectively, in M and R.

With every value of Jspin, shell corresponds a value of Ωshell = Jspin, shell/Ishell, which is a characteristic value of the angular velocity in the shell. A magnetic field is created when Ωshell  >   Ωcore. The value of Ω cannot, however, change discontinuously from Ωcore to Ωshell at the interface between core and shell.

In order to accomplish the continuity of Ω(r) and to avoid friction between the outer edge of the core and the inner surface of the shell, we chose a radius-dependent angular velocity Ω(r) raising continuously from Ωcore at the interface between core and shell to Ωedge at the edge of the gainer. Although any other choice would have been possible, we assumed that this growth follows the shape of an ellipse, somewhat similar to the shape of the tacholine in the sun as used by Spruit (2002) to explain the observed general magnetic field of the sun (Br about 1 Gauss).

Requiring that the amount of angular momentum in the shell remains the same as for a rigidly rotating shell with uniform angular velocity Ωshell, one obtains:

Ω edge = 4 π ( Ω shell Ω core ) + Ω core . $$ \begin{aligned} \Omega _\mathrm{edge} = {4\over \pi } \left(\Omega _\mathrm{shell} -\Omega _\mathrm{core} \right) + \Omega _\mathrm{core} . \end{aligned} $$(8)

3. Evolution of a binary including magnetic braking

3.1. Conservation of angular momentum

The total angular momentum JΣ of a binary is the sum of the orbital angular momentum Jorb and the spin angular momenta of gainer Jg and donor Jd. Tides continuously exchange amounts of ΔJorb, ΔJg and ΔJd. In the conservative case, one obtains:

Δ J Σ = Δ J orb + Δ J g + Δ J d = 0 . $$ \begin{aligned} \Delta {J}_{\Sigma }=\Delta {J}_\mathrm{orb} + \Delta {J}_\mathrm{g} + \Delta {J}_{d} = 0. \end{aligned} $$(9)

Angular momentum loss due to stellar wind (SW), mass loss from the system during liberal eras of RLOF (OUT) and angular momentum loss due to magnetic braking (MAG) do not violate the law of conservation of angular momentum, but change the angular momentum balance into:

Δ J orb + Δ J g + Δ J d Δ J d + g SW + OUT Δ J g MAG = 0 . $$ \begin{aligned} \Delta {J}_\mathrm{orb} + \Delta {J}_\mathrm{g} + \Delta {J}_{d} - \Delta {J}_{d+\mathrm{g}}^\mathrm{SW+OUT} - \Delta {J}_\mathrm{g} ^\mathrm{MAG} = 0. \end{aligned} $$(10)

The values of Δ J d+g SW+OUT $ \Delta {J}_{d+{\rm g}}^\mathrm{SW+OUT} $ and Δ J g MAG $ \Delta {J}_\mathrm{g}^\mathrm{MAG} $ are negative. The amount of angular momentum Δ J g MAG $ \Delta {J}_ \mathrm{g}^\mathrm{MAG} $ has only been lost by the gainer’s shell.

Since the core rotates as a solid body, no magnetic field (and hence no magnetic braking) is created in the core.

In this paper, we emphasise the effect of magnetic braking, since this was not included in our previous papers. The amount of angular momentum lost by the gainer’s shell due to magnetic braking is given by Dervisoglu et al. (2010) in cgs units:

Δ J g MAG = ( 4.807043 × 10 40 ( d M d t ) 3 7 ( R ) 24 7 ( B ) 8 7 ( M ) 2 7 Ω edge Δ t ) . $$ \begin{aligned} \Delta {J}_\mathrm{g} ^\mathrm{MAG} =-\left(4.807043 \times 10^{40}\left(\mathrm{d} M\over \mathrm{d} t\right)^{3\over 7}\left({R}\right)^{24\over 7}\left({B}\right)^{8\over 7}\left({M}\right)^{-2\over 7}\Omega _\mathrm{edge} \Delta {t}\right). \end{aligned} $$(11)

Here, the quantity B is the radial component of the magnetic field produced by the dynamo model of Spruit (2002) expressed in Gauss. The quantity d M d t $ \mathrm{d}M\over \mathrm{d}t $ is the stellar wind mass loss of the gainer in M year $ {M}_{\odot}\over\mathrm{year} $. This wind carries matter and angular momentum into space. Furthermore, R is the radius of the gainer in R, M is the mass of the gainer expressed in M and Ωedge is the angular velocity at the edge of the gainer, defined in relation (8). So, veq, g = ΩedgeReq, g.

Relation (11) is the expression derived by Dervisoglu et al. (2010) for a magnetic dipole, which is valid in this case.

3.2. Calculating the orbital period

The orbital angular momentum was calculated from relation (10). From this value the orbital period is calculated for circular orbits in cgs units with:

J orb = 1.045064 × 10 51 P 1 3 M d M g ( M d + M g ) 1 3 . $$ \begin{aligned} {J}_\mathrm{orb} =1.045064 \times 10^{51}~{P}^{1\over 3}~\frac{{M}_{d}~{M}_\mathrm{g} }{\left({M}_{d}+{M}_\mathrm{g} \right)^{1\over 3}}. \end{aligned} $$(12)

The orbital period P is in days and the masses in M.

3.3. Extent of the shell applied to the evolution of a progenitor of λ Tau

The gainer is now composed of an inner core in solid rotation surrounded by a shell that rotates differentially. The core is not spun up nor slowed down magnetically. Its rotation is only modulated by tides. The shell is spun up following relation (7). Its rotation is modulated by tides and undergoes magnetic braking following relation (11).

A core that takes a core fraction CF of the total mass of the gainer leaves (1-CF) of the mass for the shell.

Figure 1 shows the evolution with time of the equatorial velocity of the gainer starting from a (6.36 + 2,7) M binary with an initial orbital period of 1.89745 d. This is a plausible progenitor for λ Tau. The calculation was performed for the system undergoing tides and magnetic braking with initial value of CF = 0.95. This evolution with time is compared with the evolution of the gainer in the same system with a rigidly rotating and hence never magnetic gainer for which the rotation is affected by tides only. The evolution towards the point when neither tides nor magnetic braking are active is also included in the figure, which shows then a gainer rotating continuously with critical velocity once this velocity is achieved.

thumbnail Fig. 1.

Evolution with age (in years) of the equatorial velocity (in km s−1) of the gainer of a (6.36 M + 2.7 M) binary with an initial period of 1.89745 days. A likely progenitor of λ Tau. With the black line, tides and magnetic braking are not at work. With the red line tides act alone. With the blue line tides and magnetic braking act together.

The plausible progenitor for λ Tau starts RLOF during core hydrogen burning of the donor. Critical rotation of the gainer is achieved during a phase of rapid RLOF at around 48 million years after ZAMS. After that, RLOF occurs at a lower speed. Tidal interaction and magnetic braking are then strong enough to synchronize the rotation of the gainer. In Figure 1 one sees clearly that synchronization settles more rapidly when magnetic braking helps tidal interaction. After around 65 million years, during the hydrogen shell burning of the donor, RLOF restarts. Critical rotation is again attained. But the rush to critical rotation is slowed down by the combined action of tides and magnetic braking. The present observed state of λ Tau occurs during this stage of rising rotational velocity. The gainer of λ Tau now has an equatorial velocity of 178 km s−1, far below the critical value of 552 km s−1.

When critical rotation is achieved during RLOF B the orbital period is about 30 days. Figure 1 shows that tides are then too weak to prevent critical rotation. Magnetic braking will, however, decrease the velocity below the critical value. This work is done at the very end of, and after, RLOF B.

Figure 2 shows the evolution of the magnetic field of the gainer in the binary shown in Fig. 1. The magnetic field is generated by the difference of angular velocity Ω between the shell and core. Before RLOF A, the shell and core of the gainer have the same synchronous rotational velocity. The core always continues to rotate almost synchronously. When RLOF A starts, only the shell is spun up. When the shell rotates at critical velocity, the magnetic field is at a maximum of approximately 3000 Gauss. When tides and magnetic braking have synchronised the rotation of the shell, the magnetic field disappears. From the beginning of RLOF B, the shell is again spun-up, and the magnetic field starts to be built up again. The value of the magnetic field is lower than during rapid RLOF A. A value of about 750 Gauss is reached. Since magnetic braking is then nevertheless very active, the shell is gradually synchronised. After that, the magnetic field disappears again.

thumbnail Fig. 2.

Evolution of the magnetic field strength (in Gauss) produced by the Spruit mechanism for the binary mentioned in Fig. 1.

In summary, the magnetic field is large at high equatorial velocities (large ΔΩ) and small at low equatorial velocities (small ΔΩ). Figure 2 shows the eras where the gainer can be considered as a magnetic star.

A progenitor is plausible when the calculated positions of donor and gainer in the HRD are found to be close to the observations. Figure 3 shows the evolutionary path of the plausible progenitor of λ Tau for which the evolution of the equatorial velocity and of the magnetic field strength of the gainer are shown in Figs. 1 and 2.

thumbnail Fig. 3.

Evolution through the HRD of the binary taken as an example in Figs. 1 and 2. The donor follows the blue path and its current observed position is a triangle, whereas the observed rejuvenated gainer follows the red path and its current position is a filled circle.

Calculations show that choices different from an initial value of CF = 0.95 (5% of the mass of the gainer in the shell), would hardly change the results shown in Figs. 1 and 2. Only an initial value of CF = 0 is of course not justified, because this option produces a solid rotator in which no magnetic field is produced.

Figure 4 shows the evolution of the quantity CF for the gainer of the binary shown in Figs. 13. The gainer has an initial mass of 2.7 M. With CF = 0.95 one has a shell with an initial mass of 0.135 M. The solidly rotating core keeps its mass of 2.565 M during the entire evolution. At the end of RLOF A, the mass of the gainer has grown up to 6.936 M. Since the core maintains its initial mass, the mass of the shell raises up to 4.371 M, containing 63% of the gainer’s mass (CF = 0.37). At the end of RLOF B, the mass of the gainer is again somewhat larger and the value of CF shrinks to 0.31.

thumbnail Fig. 4.

Mass fraction of the core (CF) decreases with time during the evolution of the system shown in Figs. 13.

This is not so different from the evolution of a shell with, for example, an initial value of CF = 0.5 that ultimately evolves into a massive shell with CF = 0.2 (80% of the mass of the gainer being in the shell), giving rise to magnetic fields and magnetic braking that are very similar to those obtained with an initial value of CF = 0.95. This value was used in all our calculations.

3.4. Calculated cases

We found 44 binaries for which the equatorial velocity of the gainer was published. These systems are shown in Tables 1 and 2.

Table 1.

Calculated equatorial velocities of gainers that best fit observations.

Table 2.

Calculated equatorial velocities of gainers fitting observations less good than in Table 1.

For every binary in these tables, we calculated the evolution of a number of possible progenitors. In the case of conservative evolution, the initial orbital period P0 of every progenitor was determined from the presently observed period P1 with relation (13), which obeys to the rule of conservation of angular momentum:

P 0 P 1 = ( M d , 1 × M g , 1 M d , 0 × M g , 0 ) 3 . $$ \begin{aligned} \frac{{P}_{0}}{{P}_{1}}=\left(\frac{{M}_{d,1}\times {{M}_\mathrm{g,1} }}{{M}_{d,0}\times {{M}_\mathrm{g,0} }}\right)^{3}. \end{aligned} $$(13)

The subscript d is for the donor and g for the gainer. Although we calculated the orbital period with relations (10)–(12) using our code, relation (13) yields correct initial periods in the case of conservative evolution. This demonstrates that neglecting the spin angular momenta of the individual stars, as done in relation (13), is of minor importance.

In the case of liberal evolution initial periods are determined with relation (14). This relation is valid when all matter that leaves the system takes only the orbital angular momentum of the involved stars.

P 0 P 1 = ( M d , 1 × M g , 1 M d , 0 × M g , 0 ) 3 × ( M 1 M 0 ) 2 . $$ \begin{aligned} \frac{{P}_{0}}{{P}_{1}}=\left(\frac{{M}_{d,1}\times {{M}_\mathrm{g,1} }}{{M}_{d,0}\times {{M}_\mathrm{g,0} }}\right)^{3}\times \left(\frac{{M}_{1}}{{M}_{0}}\right)^{2}. \end{aligned} $$(14)

Here, M is the total mass of the system, so that:

M 1 = M 0 | Δ M SW | | Δ M OUT | . $$ \begin{aligned} {M}_{1}~={M}_{0}~-|\Delta {M}^\mathrm{SW} |-|\Delta {M}^\mathrm{OUT} |. \end{aligned} $$(15)

The quantity ΔMSW, the amount of mass that leaves the system as stellar wind, is always small compared to the stellar masses. The quantity ΔMOUT is the amount of RLOF mass that is not accreted by the gainer. This amount can be large and is difficult to predict. In our sample of 44 binaries, only two of them (V356 Sgr and TU Mon) experienced a liberal era during their life. The calculated masses of their components correspond less well to the observed values when compared to the cases with conservative evolution. For the liberal evolutions, relation (14) turns out to give less reliable initial periods, except in the case that the assumed value of ΔMOUT proves to have exactly the same value as calculated by our code.

When the initial period is sufficiently large, RLOF A, as shown in Sect. 3.3 for a plausible progenitor of λ Tau, is skipped from the evolution. The first RLOF starts then during hydrogen shell burning of the donor. In this case, calculations show that it is more difficult to avoid critical rotation.

4. Results and future work

Van Hamme & Wilson (1990) defined the quantities F = v eq v sync $ F=\frac{{{v}_{\mathrm{eq}}}}{{{v}_{\mathrm{sync}}}} $ and R = F gainer 1 F crit , gainer 1 $ R=\frac{{F}_{\mathrm{gainer}}-1}{{F}_{\mathrm{crit,gainer}}-1} $. In this case, R ∈ [0–1] is a measure of rotation. Specifically, R is 0 at synchronous rotation and 1 at critical rotation. The use of the quantity R as a characteristic measure of rotation has the advantage that minor differences in the evaluation of the synchronous and critical velocities by different authors are ironed out.

Tables 1 and 2 compare the masses of the components and the equatorial velocity characteristics of the gainers as obtained by our code with the observations. For every system a consistent number of possible progenitors was calculated from birth to the moment that the gainer leaves the main sequence. The initial periods were calculated with relation (13) or (14). From the many possible progenitors only the most plausible were included in the tables. For a few systems two plausible progenitors are mentioned. The calculated numbers are taken at the presently observed orbital period. The ranking runs from perfect determination to bad determination of veq from the top of Table 1 to the bottom of Table 2. A determination of veq is considered to be good when ΔR = |Rmodel − Robs| is small (near to 0) and is unacceptable when the same quantity is large (near to 1).

The results and reservations about the tables are:

– Results are shown in Tables 1 and 2 for 44 binaries.

– The observed equatorial velocities of 16 systems are very well confirmed by theory (ΔR <  0.1).

– 17 more systems can be considered as well explained by the theory (ΔR∈ [0.11–0.35]).

– Eight more systems are weakly reproduced by theory (ΔR∈ [0.50–0.90]).

– Only three systems are poorly reproduced with the theory (ΔR >  0.90). The observed equatorial velocities are far below critical in these cases. However, model calculations from all possible progenitors result in critical rotators. What they have in common is that the calculated models appear to be in the phase of rapid RLOF at the current observed period, when up-spinning cannot be slowed down by the combined action of tides and magnetic braking. In that case, our model fails to explain the observed spin characteristics of the binary.

– Poorly explained cases are shown at the bottom of Table 2.

– Braking of rapidly rotating gainers can be strengthened by the action of a magnetic torque. We invite the reader to consult, for example Tout & Pringle (1992) who spun-down rapidly rotating convective single stars with a magnetic torque due to the interaction with their accretion disc. Martin et al. (2011) used decretion discs for the down-spinning of Be stars. Wickramasinghe et al. (2014) dissipated the differential rotation into rigid rotation of a recently formed star resulting from the merging of two protostars into a star that owes its magnetic field to its differential rotation that has disappeared in the mean time. These examples cannot be simply applied to the gainer in a binary that evolves in a very different way. Additionally, RLOF will maintain differential rotation in the shell of the gainer. Nevertheless, the exchange of angular momentum between the (in our simple model) fast differentially rotating magnetic shell and the slower rigidly rotating non-magnetic core can reduce the value of the equatorial velocity of the gainer, and this will be included in our future research.

– Current orbital periods are well known. This is, however, less so the case for current stellar masses. It is therefore not impossible that some observed stellar masses need to be revised. In that case, they have to be approached with a different set of possible progenitors, (perhaps) yielding “better” results.

– Observed values of equatorial velocities below synchronous velocity are usually not reproduced by our calculations since tides always tend to synchronize the rotation.

5. Conclusions

For systems with small initial orbital periods, our models predict two stages with strong magnetic fields. The first one occurs during phase A of RLOF at a high mass-transfer rate and lasts rather a short period of time, In the example discussed in Sect. 3.3 the mass transfer rate goes up to (5 × 10−5) ( M yr ) $ \left(\frac{{M}_{\odot}}{\mathrm{yr}}\right) $, and the duration of a large magnetic field (above 1 kG with a maximum of 3 kG) lasts for 3 × 105 years. At that stage, the spectral lines are heavily rotationally broadened by the critical rotation.

The second stage, during phase B of RLOF, lasts longer during a period of lower mass-transfer. However, the mass transfer rate goes up to 10−6 M yr $ {M}_{\odot}\over\mathrm{yr} $ and the duration of a large magnetic field lasts for 107 years. This magnetic field strength is less high compared to the one that is active during RLOF A and shrinks from a maximum of 750–100 Gauss. Then, the spectral lines then are again rotationally broadened due to the critical rotation of the gainer.

The detection of Zeeman effect at work in the spectra of gainers rotating at large equatorial velocities might be problematic. But it would offer a test for the validity of our model.

For systems with larger initial orbital periods RLOF A will be skipped from the evolution. Straight RLOF B will spin the gainer up, more or less quickly, to critical velocity. This critical velocity remains until the end of RLOF B. At that moment, the orbital period is large and critical velocities cannot be dragged below this limit by weak tidal interaction. Magnetic braking will be the major mechanism that will ultimately synchronize the gainer.

Despite the fact that we used a simple model for differential rotation of the gainer in order to generate magnetic braking, we found that 75% of the observed rotational velocities of gainers in interacting binaries were reproduced reasonably well by our model. In less than 10% of the examined cases, our model failed completely to reproduce the observations. This might show that our model is not applicable in these cases, but also that observers could evaluate new values for equatorial velocities of gainers and masses of the components in interacting binaries. Different masses indeed yield different sets of possible progenitors that reproduce the characteristics of the current binary. On the other hand the observed orbital periods are correctly listed.

Researchers who want to obtain calculations with our current binary evolutionary code are welcome to request this from us.

Acknowledgments

We thank Henk Spruit for explaining patiently how his dynamo works. We thank the referee and Edward Van den Heuvel for their constructive and helpful comments.

References

  1. Brancewicz, H., & Dworak, T. 1980, Acta Astron., 30, 501 [Google Scholar]
  2. Budding, E., Erdem, A., Cicek, C., et al. 2004, A&A, 417, 263 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  3. De Jager, C., Nieuwenhuyzen, H., & Van der Hucht, K. 1988, A&AS, 72, 259 [Google Scholar]
  4. Dervisoglu, A., Tout, C., & Ibanglu, C. 2010, MNRAS, 406, 1071 [Google Scholar]
  5. Dominis, D., Mimica, P., Pavlovski, K., & Tamajo, E. 2005, ApSS, 296, 189 [Google Scholar]
  6. Eggleton, P., & Kiseleva-Eggleton, L. 2002, ApJ, 575, 461 [NASA ADS] [CrossRef] [Google Scholar]
  7. Glazunova, L., & Yushcenko, A. 2008, AJ, 135, 1736 [NASA ADS] [CrossRef] [Google Scholar]
  8. Hilditch, R. 2001, An Introduction to Close Binary Stars (Cambridge University Press) [CrossRef] [Google Scholar]
  9. Martin, R., Pringle, J., Tout, C., & Lubow, H. 2011, MNRAS, 416, 2827 [NASA ADS] [CrossRef] [Google Scholar]
  10. Miller, B., Budaj, J., & Richards, M. 2007, ApJ, 656, 1075 [NASA ADS] [CrossRef] [Google Scholar]
  11. Packet, P. 1981, A&A, 102, 17 [Google Scholar]
  12. Song, H., Meynet, G., Maeder, A., Ekström, S., et al. 2018, A&AS, 609, A3 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  13. Spruit, H. R. 2002, A&A, 381, 923 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  14. Tassoul, J. L., 2000, Stellar Rotation (Cambridge University Press) [CrossRef] [Google Scholar]
  15. Tout, C., & Pringle, J. 1992, MNRAS, 256, 269 [CrossRef] [Google Scholar]
  16. Vanbeveren, D., Mennekens, N., Shara, M., & Moffat, A. 2018, A&A, 615, A65 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  17. Van Hamme, W., & Wilson, R. 1990, AJ, 100, 1982 [NASA ADS] [CrossRef] [Google Scholar]
  18. Van Rensbergen, W., & De Greve, J. P. 2016, A&A, 592, A151 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  19. Van Rensbergen, W., De Greve, J. P., et al. 2008, A&A, 487, 1129 Vizier Online Data Catalog: J/A&A/487/1129 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  20. Van Rensbergen, W., De Greve, J. P., Mennekens, N., et al. 2010, A&A, 510, A13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  21. Van Rensbergen, W., De Greve, J. P., Mennekens, N., et al. 2011, A&A, 528, A16 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  22. Vink, J., De Koter, A., & Lamers, H. 2001, A&A, 369, 574 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  23. Wellstein, S. 2001, PhD Thesis, Potsdam University [Google Scholar]
  24. Wickramasinghe, D., Tout, C., & Ferrario, L. 2014, MNRAS, 437, 675 [NASA ADS] [CrossRef] [Google Scholar]

All Tables

Table 1.

Calculated equatorial velocities of gainers that best fit observations.

Table 2.

Calculated equatorial velocities of gainers fitting observations less good than in Table 1.

All Figures

thumbnail Fig. 1.

Evolution with age (in years) of the equatorial velocity (in km s−1) of the gainer of a (6.36 M + 2.7 M) binary with an initial period of 1.89745 days. A likely progenitor of λ Tau. With the black line, tides and magnetic braking are not at work. With the red line tides act alone. With the blue line tides and magnetic braking act together.

In the text
thumbnail Fig. 2.

Evolution of the magnetic field strength (in Gauss) produced by the Spruit mechanism for the binary mentioned in Fig. 1.

In the text
thumbnail Fig. 3.

Evolution through the HRD of the binary taken as an example in Figs. 1 and 2. The donor follows the blue path and its current observed position is a triangle, whereas the observed rejuvenated gainer follows the red path and its current position is a filled circle.

In the text
thumbnail Fig. 4.

Mass fraction of the core (CF) decreases with time during the evolution of the system shown in Figs. 13.

In the text

Current usage metrics show cumulative count of Article Views (full-text 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 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.