Properties of six short-period massive binaries: a study of the effects of binarity on surface chemical abundances††thanks: Based on observations made with the SOPHIE spectrograph on the 1.93-m telescope at Observatoire de Haute-Provence (OHP, CNRS/AMU), France
Key Words.:Stars: early-type – Stars: atmospheres – Stars: fundamental parameters – Stars: abundances – Stars: binaries: general
Context:A significant fraction of massive stars are found in multiple systems. The effect of binarity on stellar evolution is poorly constrained. In particular, the role of tides and mass transfer on surface chemical abundances is not constrained observationally.
Aims:The aim of this study is to investigate the effect of binarity on the stellar properties and surface abundances of massive binaries.
Methods:We perform a spectroscopic analysis of six Galactic massive binaries. The spectra of individual components are obtained from a spectral disentangling method and are subsequently analyzed by means of atmosphere models. The stellar parameters and CNO surface abundances are determined.
Results:Most systems are made of main-sequence stars. Three systems are detached, two are in contact and no information is available for the sixth system. For eleven out of the twelve stars studied the surface abundances are only mildly affected by stellar evolution and mixing. They are not different from those of single stars, within the uncertainties. The secondary of XZ Cep is strongly chemically enriched. Considering previous determinations of surface abundance in massive binary systems suggests that the effect of tides on chemical mixing is limited, whereas mass transfer and removal of outer layers of the mass donor leads to the appearance of chemically processed material at the surface, although this is not systematic. The evolutionary masses of the components of our six systems are on average 16.5% higher than the dynamical masses. Some systems seem to have reached synchronization, while others may still be in a transitory phase.
The evolution of massive stars depends on several physical processes, the main ones being mass loss and rotation. Mass loss removes material from the outer layers, reducing the mass of the star and thus affecting its internal structure and luminosity (Chiosi & Maeder 1986). Main-sequence massive stars rotate on average faster than any other main-sequence star. Their average projected rotational velocity is found to be between 100 and 150 km s depending on metallicity (Penny & Gies 2009; Huang et al. 2010; Ramírez-Agudelo et al. 2013). Rotation flattens the star, triggering the transport of angular momentum and chemical species from the core to the surface (Maeder & Meynet 2000). Rotation also modifies mass loss rates. Other processes affect the evolution of massive stars. The presence of a strong dipolar magnetic field (present in about 7% of OB stars, Grunhut et al. 2017) can alter both mass loss and rotation (Ud-Doula et al. 2009; Meynet et al. 2011). The presence of a companion can also modify the evolution of a star through tides and mass transfer (e.g. Mandel & de Mink 2016).
A number of massive OB stars are found in binary systems. However the fraction of multiple systems is not fully established and varies depending on samples and environment. Kobulnicky et al. (2014) reported that 35% (55% after corrections for biases) of the 128 OB stars they studied in Cyg OB2 definitely have a companion ; Sana et al. (2012) found a binary fraction of 56% (69% after corrections) for six Galactic open clusters and a total of 71 systems (40 definite binaries); Sana et al. (2014) detected a companion in 53% of a sample of 96 Galactic O stars and claimed a binary fraction of 91% accounting for known systems not detected in their survey (many systems having large separations though); in the Large Magellanic Cloud, Sana et al. (2013) determined an observed fraction of 35% (corrected to 51%) among 360 O stars (126 definite binaries); a minimum fraction of 21% was determined by Mahy et al. (2013) for four Cygnus OB associations (19 systems, 4 definite binaries). In addition, Kobulnicky et al. (2014) and Sana et al. (2012) found a rather uniform distribution of mass ratios and a weak decrease of the period distribution among definite binaries. The latter distribution extends up to several thousands of days.
The presence of a companion around a significant number of OB stars raises the question of its effects on stellar evolution. Interactions are expected in close systems. For wide separations, both components evolve as single stars. There are two main categories of effects caused by the presence of a close companion: those due to tides and those due to mass transfer. The former trigger energy exchange between the system and its components. The consequence is a change of the internal structure under the influence of energy dissipation (and thus heating) through radiative damping, viscous friction and (gravito-)inertial waves (Zahn 1989; Mathis & Remus 2013). The geometry of the star and its internal rotation profile can also be modified. Hence the transport of angular momentum and chemical species are affected. In practice, these effects lead to synchronization of the rotational and orbital periods (Zahn 1977). This can spin-up or spin-down individual components and thus amplify or reduce the effects of rotation on their evolution. Mass transfer also drastically affects the evolution of massive stars. Due to the mass-luminosity relation, a decrease of mass implies a reduction of luminosity, and thus a change of the path in the Hertzsprung-Russell diagram (Wellstein et al. 2001). The mass donor loses its external layers, revealing at the surface hotter internal regions with different chemical properties. The mass gainer is polluted by material from the donor and its envelope is mixed with new material (Langer et al. 2008). Together with mass, angular momentum is exchanged. This affects the rotational properties of both stars, in a way that depends on the efficiency of mass transfer. The outcome of mass transfer depends on the orbital properties of the system (mass ratio, separation) and on the intrinsic evolution of the components.
Stellar evolution in (massive) binaries is thus complex. Constraints from analysis of dedicated systems are therefore highly requested. One of the expectations of binary evolution is that surface abundances are modified compared to single stars. Indeed, mixing processes due to rotation depend on the star’s internal structure which is affected by tides (de Mink et al. 2009). Contamination by accretion in mass transfer episodes obviously affects surface chemical patterns (Langer et al. 2008). In this case, the mass donor loses its external layers so that its surface is physically moved to deeper layers where chemical composition corresponds to more processed material. These effects of binarity on surface abundances are one of the possibilities quoted to explain peculiar abundance patterns. For instance, the nitrogen-rich stars with low projected rotational velocity (as well as their counterparts, nitrogen-poor stars with high sin) reported by Hunter et al. (2008) – see also Grin et al. (2017) – may be the result of such interactions according to the authors. However the difference between the chemical patterns of binary and single stars has not been established observationally. The attempt performed by Garland et al. (2017) on a sample of B-type binaries did not show any clear difference, but only the primary component of their systems was analyzed.
Here we present a study of six Galactic massive binaries. Our goal is to determine the stellar parameters and surface CNO abundances of their components in order to investigate whether they are different from presumably single stars. Section 2 describes our sample and the spectroscopic observations. Data analysis and the associated results are gathered in Sect. 3. Results are discussed in Sect. 4. Finally we summarize our conclusions in Sect. 5.
2 Sample and observations
In order to test the effects of binarity on the properties and evolution of massive stars, we selected six short period systems ( days) from the compilation of Gies (2003). All systems except DH Cep are eclipsing binaries. The targets presented in this study were observed at Observatoire de Haute-Provence with the SOPHIE instrument (Bouchy et al. 2013) during the nights of 22 to 26 August 2013. SOPHIE delivers high spectral resolution échelle spectra in the 3900-6900Å wavelength range. The main properties of the observed binaries are listed in Table 1. The High Efficiency mode of SOPHIE, corresponding to a spectral resolution , was used. The spectra were obtained in 15 to 40 minutes depending on the target’s brightness and have a signal-to-noise ratio close to 200 in the continuum. The observing strategy (number of observations per system) was adjusted each night on order to ensure an optimal coverage of the phases of maximum separation between the systems’ components as requested for a good spectral disentangling. Data reduction was performed automatically by the SOPHIE pipeline (Bouchy et al. 2009) adapted from the ESO/HARPS software.
3 Modelling and spectroscopic analysis
3.1 Orbital solutions and spectral disentangling
To compute the orbital solution of each system, we first measured the radial velocities of each component. For this purpose, we fitted the spectral lines of helium with the highest ionization stage. The higher the ionization stage, the closer to the photosphere the lines are formed. The computations of the radial velocities (RVs) were done by adopting the rest wavelength from Conti et al. (1977). We then used the Liège Orbital Solution Package (LOSP111LOSP is developed and maintained by H. Sana. The algorithm is based on the generalization of the SB1 method of Wolfe et al. (1967) to the SB2 case along the lines described in Rauw et al. (2000) and Sana et al. (2006). It is available at http://www.stsci.edu/~hsana/losp.html.) to determine the SB2 orbital solution of the system. Since the orbital periods of these systems are well known, we decided to fix them in the computations. Except for Y Cyg, the other systems have all a circular orbit. The orbital solutions are displayed in Fig. 1. The corresponding parameters are gathered in Table 2.
We then used the orbital parameters as input of our spectral disentangling code. We applied the Fourier approach of Hadrava (1995) to separate the spectral contributions of each component in the different systems. This method uses the Nelder & Meadâs Downhill Simplex method on the multidimensional parameter space to reach the best fit between the recombined component spectra and the observed data. Discussions about the spectral disentangling method used in this paper can be found in Pavlovski & Hensberge (2010) and in Mahy et al. (2017). Because the eclipses are not total in the different systems, the output composite spectra must be renormalized and corrected for the brightness ratios. A first approximation of the brightness ratio was estimated by computing the equivalent widths of several spectral lines (Si iv 4089, He i 4143, He i 4471, He ii 4542, O iii 5592) and by comparing these values to those calculated from synthetic spectra corresponding to stars with the same spectral types as the components of our sample. Then, the brightness ratios have been iteratively adjusted by computing them from the ratios between the luminosities of the two components (luminosities that have been obtained from the effective temperature and the radius of each star).
|AH Cep||DH Cep||XZ Cep|
|[day]||1.774727 (fixed)||2.11095 (fixed)||5.097253 (fixed)|
|[HJD 2 450 000]||6526.783 0.008||6525.564 0.006||6522.656 0.014|
|1.132 0.056||1.149 0.020||2.013 0.092|
|[km s]||236.46 7.33||267.57 8.29||234.81 3.97||269.70 4.56||122.86 3.49||247.29 7.03|
|||8.29 0.26||9.38 0.29||9.79 0.17||11.24 0.19||12.37 0.35||24.89 0.71|
|[M]||12.49 0.78||11.04 0.66||15.01 0.66||13.07 0.57||17.88 1.14||8.88 0.46|
|rms [km s]|
|V382 Cyg||V478 Cyg||Y Cyg|
|[day]||1.885545 (fixed)||2.880867 (fixed)||2.9963316 (fixed)|
|[HJD 2 450 000]||6527.020 0.003||6525.763 0.006||6526.454 0.041|
|1.376 0.009||1.046 0.040||1.053 0.017|
|[km s]||257.32 1.46||354.11 2.01||222.11 5.43||232.28 5.67||232.41 3.00||244.81 3.16|
|||9.58 0.05||13.19 0.08||12.64 0.31||13.22 0.32||13.64 0.18||14.37 0.19|
|[M]||25.85 0.38||18.78 0.27||14.31 0.72||13.68 0.68||16.87 0.53||16.02 0.50|
|rms [km s]|
3.2 Stellar parameters
We used the Fourier transform method (Gray 1976; Simón-Díaz & Herrero 2007) to determine the projected rotational velocity sin. When the combined spectra showed sufficiently wide separations between both components of the system so that individual lines could be clearly isolated, we used directly the combined spectrum. Otherwise, we relied on the disentangled spectra. The uncertainty of our determinations is 10 km s. We found that all stars have sin larger than 100 km s.
To estimate the amount of extra-broadening (macroturbulence) we used synthetic spectra (see below) with and consistent with the spectral type and the luminosity class of each star (Martins et al. 2005; Nieva & Przybilla 2014). We convolved these spectra with a rotational profile adopting the sin derived from the Fourier transform method. We added an extra-broadening by means of a radial-tangential profile parameterized by a velocity . The comparison with the observed, disentangled spectrum of each component revealed that for most systems there was no need for an extra-broadening. We thus adopted = 0 km s. Only for DH Cep did we choose = 100 km s in order to correctly reproduce the line wings.
and were determined from the classical diagnostics: the helium ionization balance for and the Balmer lines broadening for the surface gravity (e.g. Martins 2011). For B stars we replaced helium by silicon for the effective temperature determination. We relied on a grid of synthetic spectra computed with the code CMFGEN (Hillier & Miller 1998). Non-LTE atmosphere models were computed for stars with 21000 45000 K and 3.0 4.3. Luminosities were assigned according to the calibration of Martins et al. (2005) and mass loss rates were computed from Vink et al. (2001). Wind terminal velocities were estimated from = 3.0 (Garcia et al. 2014). A solar composition (Grevesse et al. 2010) was assumed and the following elements were included: H, He, C, N, O, Ne, Mg, Si, S, Ar, Ca, Fe, Ni. Once atmosphere models were converged a formal solution of the radiative transfer equation lead to the synthetic emergent spectrum, which was computed assuming a microturbulent velocity ranging from 10 km s at the photosphere to 0.1 in the outermost parts of the atmosphere. For each star, a subset of the model grid (encompassing the values of and estimated from the spectral type) was convolved with the appropriate sin and as described above. The resulting synthetic spectra, correctly shifted in radial velocities, were subsequently compared to the disentangled spectrum of each target. The quality of the fit was quantified by means of a analysis renormalized to 1.0 for the best-fit combination of and . The output of such a process is illustrated in Fig. 2 in the case of the primary component of Y Cyg. The uncertainties on and were estimated from such figures. In practice, and are correlated and giving a formal error on and is not correct. However, we decided to take the width of the contour as a representative error. These values are reported in Table 3.
Surface abundances of C, N and O were determined as in Martins et al. (2015): for a given set of and models with different C, N and O content were computed and subsequently compared to disentangled spectra. The goodness of fit was quantified by means of a analysis. In practice, we relied on the following lines:
for carbon: C iii 4068-70, C iii 4153, C iii 4156, C iii 4163, C iii 4187, C ii 4267, C ii 6578, and C ii 6583.
for nitrogen: N ii 3995, N ii 4004, N ii 4035, N ii 4041, N iii 4044, N ii 4236, N ii 4447, N iii 4511, N iii 4515, N iii 4518, N iii 4524, N ii 4530, N ii 4601, N ii 4607, N ii 4614, N ii 4621, N ii 4630, N ii 4788, N ii 4803, N ii 5001, N ii 5005 and N ii 5680.
for oxygen: O ii 3913, O ii 3955, O ii 4120, O ii 4284, O ii 4305, O ii 4318, O ii 4321, O ii 4368, O ii 4592, O ii 4597, O ii 4603, O ii 4611, O ii 4663, O ii 4678, O ii 4700, O ii 4707, and O iii 5592.
The number of lines used for each star depends on the effective temperature and on the quality of the spectra. We adopted He/H=0.1 in our computations. No significant deviation from this value was observed in our fits (see Fig. 3).
The derived parameters are given in Table 3. The final best-fit models are shown in Fig. 3. For each system, the best fit of the individual components were co-added taking into account the luminosity ratio and the radial velocities determined at the date at which the comparison spectrum was obtained. In general, the quality of the fits is good. Notable problems are encountered around Balmer lines, especially H, illustrating the difficulty to correctly disentangle broad lines that are never completely separated observationally. This is even worse when their wings contain lines from other elements.
|[kK]||[km s]||[km s]|||||||
Except for DH Cep, the five other systems display eclipses. We thus used the PHOEBE (PHysics Of Eclipsing BinariEs, Prša & Zwitter 2005, v0.31a) software to model the Hipparcos lightcurves (Perryman et al. 1997). For their analysis, we kept the parameters of the orbit fixed at their values obtained through the spectroscopic analysis (Section 3.1). PHOEBE allows to model the lightcurve and the RV curves of an object at the same time. This software is based on Wilson & Devinneyâs code (Wilson & Devinney 1971) and uses Nelder & Meadâs Simplex fitting method to adjust all the input parameters to find the best fit to the lightcurve. In the search for the best-fit model, we include in the computation of the lightcurve the reflection effect. The latter is reported as the brightening of one star due to irradiation by its companion and has a non-negligible impact on very luminous systems when the two component orbit very close to each other. The reflection effect can affect the curvature of the lightcurve in zones outside the eclipse. Reflection is treated according to the theory of Wilson (1990) in the code PHOEBE. We also fixed the effective temperatures to the values obtained in the previous section. The photometric parameters derived by PHOEBE are given in Table 4. The error bars are computed by fixing one parameter at a time and we allow the others to vary to reach the minimum of the corresponding to a 68.3% confidence level (1-).
|AH Cep||XZ Cep||V382 Cyg||V478 Cyg||Y Cyg|
4.1 Comparison to previous analysis
|V478 Cyg||15.30.8||14.60.8||7.60.1||8.50.1||this work|
|16.60.9||16.30.9||7.430.12||7.430.12||Popper & Hill (1991)|
|Y Cyg||16.90.5||16.00.5||5.90.1||5.10.1||this work|
|16.80.5||17.70.4||–||–||Burkholder et al. (1997)|
|17.720.35||17.730.30||5.7850.091||5.8160.063||Harmanec et al. (2014)|
|17.570.27||17.040.26||5.930.07||5.780.07||Simon et al. (1994)|
|V382 Cyg||26.10.4||19.00.3||9.40.2||8.70.2||this work|
|32.61.8||22.91.3||8.80.6||7.40.6||Popper & Hill (1991)|
|29.71.1||20.30.9||–||–||Burkholder et al. (1997)|
|26.00.7||19.30.4||9.60.1||8.40.1||Harries et al. (1997)|
|XZ Cep||18.71.3||9.30.5||14.20.1||14.20.1||this work|
|15.80.4||6.40.3||7.00.2||140.50.2||Harries et al. (1997)|
|AH Cep||14.31.0||12.60.9||5.60.1||4.70.1||this work|
|16.2||13.3||–||–||Burkholder et al. (1997)|
|18.10.9||15.90.8||6.70.2||6.20.2||Bell et al. (1986)|
|DH Cep||38.42.5||33.42.2||–||–||this work|
|34.4||29.8||–||–||Burkholder et al. (1997)|
|32.71.7||29.61.6||–||–||Hilditch et al. (1996)|
In this section we compare our results with those available in the literature. Table 5 gathers the published values of masses and radii for the systems we have analyzed.
V478 Cyg was studied by Popper & Hill (1991). Within the error bars, their results are in excellent agreement with our findings with the exception of the radius of the secondary which is slightly larger in our study.
Popper & Hill (1991) also analyzed V382 Cyg. Our radii are marginally consistent with their results, while our mass estimates are lower. The masses of Burkholder et al. (1997) are in better agreement with our findings. Burkholder et al. (1997) concluded that V382 Cyg was an interacting system since both stars fill their Roche lobe. Similar conclusions were reached by Harries et al. (1997) who obtained masses and radii in excellent agreement with our results.
The mass and radius of the primary component of Y Cyg obtained by Burkholder et al. (1997), Harmanec et al. (2014) and Simon et al. (1994) are consistent with our values. On the contrary, we find a less massive and more compact secondary. We note that this system is the only one for which the phase coverage of our observations does not sample the maximum separation (Fig. 1). This probably affects our results and explains the differences with previous studies. The effective temperature and surface gravity derived by Simon et al. (1994) – = 34500 (34200) K and = 4.16 (4.18) for the primary (secondary) – are similar to our determinations, within the error bars. The same conclusion applies to the determination of the inclination of the system (=8637012 for Simon & Sturm, =8802.0 for us).
The inclination of DH Cep was estimated by Sturm & Simon (1994): they found =471. Adopting this value, we obtain masses of 38.42.5 and 33.42.2 M for the primary and secondary. Within the error bars these values are in agreement with those of Burkholder et al. (1997) and Hilditch et al. (1996).
Our effective temperatures are consistent with those of Sturm & Simon (1994) – = 44000 (43000) K for the primary (secondary) – but we find surface gravities larger by 0.3 dex (Simon & Sturm report = 3.95 (4.03) for the primary (secondary)). This likely reflects the difficulty to correctly disentangle the individual profiles of Balmer lines that are the main gravity indicators.
Finally, XZ Cep was studied by Harries et al. (1997) who obtained masses lower by 3 M compared to our determinations. On the contrary their radii are slightly larger. These differences are probably due to the contact nature of the system (see below).
4.2 Evolutionary status
Fig. 4 shows the classical Hertzsprung-Russell (HR) diagram and the - diagram. All systems but XZ Cep have both of their components on the main sequence as defined by single-star evolutionary tracks. XZ Cep is the most evolved system of our sample. For almost all stars the same track can explain the position of stars in both the HR and - diagrams. The secondary component of XZ Cep may be the only exception since its luminosity places it close to the 15 M track in the HR diagram while a mass of 20 M better reproduces its position in the - diagram. But we also note that within the error bars the 15 and 20 M tracks are relevant in both diagrams.
Table 6 provides an estimate of the Roche lobe size for each component of the five eclipsing systems. This estimate was made using equation 2 of Eggleton (1983). In the third column of Table 6 we compare the size of the Roche lobe to the stellar radius in the direction of the companion () when possible, or to the mean radius when is not available (see Table 4). V478 Cyg, Y Cyg and AH Cep are detached systems. The latter two have the largest ratios R/ which is consistent with their position very close to the ZAMS. V478 Cyg, which is more evolved, is on the verge of filling its Roche lobe (secondary star). V382 Cyg and XZ Cep are contact systems, in agreement with the findings of Harries et al. (1997). Given the above remarks regarding the position on the HR and - diagrams, XZ Cep has most likely experienced mass transfer and certainly follows an evolution different from that of single stars (see also Sect. 4.3). The properties of V382 Cyg in both evolutionary diagrams do not indicate deviation from single star evolution as strong as for XZ Cep: V382 Cyg is thus probably only in an early contact phase so that deviation from single-star evolution has not appeared yet.
It is interesting to compare Y Cyg and V478 Cyg in the light of Table 6: both stars have a mass ratio close to 1.0, dynamical masses around 15 M, and similar orbital periods (just below 3 days). According to Fig. 4 Y Cyg is located closer to the ZAMS and is less evolved than V478 Cyg which lies in the middle of the main sequence. This is reflected in Table 6 where the former system has its components filling a smaller fraction of its Roche lobe than the latter. We conclude that for an orbital period of about 3 days and stars with similar initial masses of about 15 M contact does not happen before the final portion of the main sequence.
4.3 Surface abundances
Fig. 5 shows that most of the sample stars are barely evolved in terms of surface CNO abundances. Almost all stars have C, N and O abundance patterns consistent with no processing. The only exception is the secondary of XZ Cep which is significantly enriched. The ratios N/C and N/O are consistent with CNO burning.
Fig. 6 shows a comparison of the sample stars with presumably single stars analyzed by Martins et al. (2015). We chose this study as a reference for comparison since it is based on data of similar quality and the method employed to determine the stellar parameters and surface abundances is the same as the one we used in the present study. The left panel of Fig. 6 shows the - diagram. The comparison stars are separated in two sub-samples depending on their initial mass (above or below 28 M) according to the tracks of Ekström et al. (2012). In the middle panel, the surface abundances of V478 Cyg, Y Cyg, XZ Cep and AH Cep can be compared to the single stars with masses below 28 M. V382 Cyg and DH Cep are compared to more massive single stars in the right panel. We see that in both mass ranges there is no difference between the binary and single star samples: they cover the same area of the log(N/C) – diagram. As above, the only possible difference is the secondary star of XZ Cep. Although there is no direct comparison single stars close to it, we see that its N/C ratio corresponds to predictions of single star models with initial masses above 25 M (and close to 40 M), while its estimated mass is in the range 8-20 M (see Sect. 4.4). It is thus likely that in this evolved system binarity has affected the surface chemical composition. For the other five systems, there is no evidence for such an effect, even for V382 Cyg in which both components (barely) fill their Roche lobe. In Fig. 6 we have also separated comparison stars according to whether they have sin higher than 120 km s or not. The high- sin sample encompass the range of projected rotational velocities of the binary components. Here again, there is no clear distinction between binaries and single stars.
We thus conclude that in the binary sample we have studied, detached systems do not show clear differences in their surface abundances compared to single stars. For interacting systems, large N/C may be encountered (XZ Cep) or not (V382 Cyg). A certain time may be required after the beginning of the mass transfer in order to see its effects on surface abundances.
Garland et al. (2017) studied a sample of 33 B-type binary systems in the Large Magellanic Cloud, all of them located on the main sequence. They did not perform a disentangling of the individual spectra but determined the stellar parameters of the primaries with atmosphere models. They used various assumptions regarding contamination of their spectra by the secondaries to infer the impact on the derived parameters. They estimated nitrogen surface abundances and projected rotational velocities and concluded that primary stars of their binary systems do not show differences in terms of surface enrichment compared to B-type single stars. They argue that they may have evolved effectively as single stars with low rotational velocities.
Mahy et al. (2017) determined the stellar parameters and surface abundances of the O7.5If+O9I(f) system HD 166734. Their analysis indicates that this is a detached system and that the degree of chemical processing is not different from what is observed in single O supergiants.
Mahy et al. (2011) analyzed the properties of the two components of the massive binary LZ Cep. They found that the secondary is extremely nitrogen- and helium-rich, while being carbon- and oxygen-poor. Quantitatively, (N/C)=1.6 which, together with a of 3.1, places the star in the upper right corner of Fig. 7 (the present-day dynamical mass of the secondary is 6 M). On the other hand the primary is barely chemically evolved. Mahy et al. (2011) concluded that the surface abundances could be explained by an inefficient mass transfer from the secondary (which was initially the primary) towards the secondary. The dynamical masses obtained by Mahy et al. for LZ Cep ( 16 and 6 M for the primary and secondary) are significantly lower than the evolutionary masses (25.3 and 18.0 M) determined using single star evolutionary tracks. This is another indication that mass transfer occurred in that system according to Mahy et al. (2011). In any case, the degree of chemical enrichment is higher than expected for single stars of masses 5 to 25 M– it would be consistent with the enrichment of a 50 M star.
The surface abundances of the components of the system HD149404 are similar to that of LZ Cep according to Raucq et al. (2016) – see also Fig. 7. Since HD 149404 is not an eclipsing binary, there is no determination of its dynamical mass. The evolutionary masses reported by Raucq et al. (2016) are of the order 30 M assuming a standard initial rotation rate. As for LZ Cep, the secondary is thus highly enriched, as confirmed by its position in Fig. 7: the observed N/C ratio is larger by 0.8 dex compared to the prediction of a 32 M model. Given these abundances and the asynchronous rotation of the system, Raucq et al. (2016) suggested that HD149404 is in a post-Roche lobe overflow state. Such an interaction may have revealed internal layers of the secondary star, explaining the rather strong enrichment.
The massive evolved system LSS 3074 was studied by Raucq et al. (2017). The dynamical masses of the components are 14.8 and 17.2 M for the primary and secondary respectively. These values (together with the estimated radii) are low for stars with spectral type O4f and O6-7(f). Raucq et al. determined lower limits of 0.94 (1.26) for the primary (secondary) (N/C) ratios. The primary is also helium-rich (He/H=0.25). The surface gravities of both components being close to 3.8, the stars are located in the upper part of Fig. 7. Raucq et al. (2017) indicate that LSS 3074 is likely an overcontact system which experienced Roche lobe overflow. Their simulations favour a relatively massive primary (M=35 M). In that case, the N/C ratio determined for that star is high compared to predictions of evolutionary models (see Fig. 7). Compared to the two systems described above, it is worth noting that in the case of LSS 3074 both components show a high N/C ratio.
Both components of Plaskett’s star (HD 47129) have initial masses above 40 M according to Linder et al. (2008). The primary has a high N/C ratio, while the secondary is surprisingly chemically unevolved. The secondary has a high projected rotational velocity, which prompted Linder et al. to argue that this system experienced a recent mass transfer in which a small fraction of the mass removed from the primary was accreted by the secondary. This could explain the high rotation and absence of nitrogen enrichment / carbon depletion of the secondary. In Fig. 7 the primary is located in the upper part at a position consistent with high mass single stars. We note that a strong magnetic field is present at the surface of the secondary Grunhut et al. (2013), which may have affected the interaction history and surface abundances of that system (e.g. Meynet et al. 2011).
From the above discussion, and keeping in mind that samples of binary stars for which surface abundances are available remain small, the following trend tends to emerge:
systems for which no evidence of mass transfer exist do not show surface enrichment significantly different from single stars, at least within the uncertainties on the derived surface abundances. Said differently, the additional role of tides on mixing in detached systems has a limited effect on surface abundances compared to mixing in single stars.
on the contrary in three systems which experienced mass transfer (HD 47129, LZ Cep and HD 149404), the mass donor has a high N/C ratio, even higher than what is observed in single stars of the same initial mass and rotation rate for two systems (LZ Cep and HD 149404). Our results for XZ Cep follow the same trend: the secondary fills its Roche lobe and has a larger N/C than expected for its mass. In these systems mass transfers probably removed the external layers of the donor. This process exposes at the surface internal regions more chemically mixed (an effect similar to stellar winds for Wolf-Rayet stars). The mass gainer do not show peculiar surface abundances compared to single stars. This latter result is consistent with predictions of Langer et al. (2008) who show that the accreted material is mixed and diluted in the pre-existing envelope. For HD 149404, both the primary and secondary show relatively high N/C ratios, possibly indicating that the mass gainer was already chemically enriched (because of rotational mixing) before interaction. Finally, V382 Cyg just fills its Roche lobe but do not show signs of peculiar N/C.
Peculiar N/C ratios would thus be observed in the mass donor of post-mass transfer systems, possibly because of the removal of external layers and hence the appearance of internal layers deeply chemically processed at the surface. In absence of mass transfer, surface abundances of binary components would be indistinguishable from those of single stars. The case of the contact system V382 Cyg, for which no peculiar abundances are determined, suggests that the above picture is only partial or that the system was caught at the very beginning of the mass transfer episode, when significant envelope removal did not occur yet. Obviously, this tentative picture needs to be confirmed by analysis of additional evolved systems. Such a task is difficult because the estimates of initial masses in interacting systems are based on assumptions on the efficiency of mass transfer, on the initial orbital parameters, on the initial mass ratio (e.g. Wellstein et al. 2001; de Mink et al. 2009; Raucq et al. 2017). In addition, evolved systems are not spherically symmetric and the use of 1D atmosphere models may be questioned (Palate et al. 2013).
4.4 Stellar masses
In Table 7 we present three mass estimates: the dynamical mass (M) results from the orbital solution; the evolutionary mass (M) is estimated from the position in the HR diagram (Fig. 4, right panel); the spectroscopic mass (M) is obtained from the surface gravity and mean radius.
Fig. 8 shows a comparison between these three masses. In the left panel, we see that M and M are well correlated. However, there is a systematic upward shift indicating that evolutionary masses tend to be slightly larger (by a few solar masses) than dynamical masses. Quantitatively, the former are 16.510.6% higher than the latter, with differences ranging between 8 and 38%.
Weidner & Vink (2010) studied a large sample of binaries and concluded that there was no difference between dynamical and evolutionary masses. The only system in common with our study is V478 Cyg for which Weidner & Vink (2010) give M=16.69.0 M (16.39.0 M for the secondary) and M=18 M. The dynamical masses are in good agreement with our measurements. The evolutionary masses of Weidner & Vink (2010) are slightly lower than ours, but consistent within the error bars. Vink & Weidner relied on the tracks of Meynet & Maeder (2003) which are less luminous, for a given initial mass, than those of Ekström et al. (2012) that we use. With the 2003 Geneva tracks, one finds M=20.0 M for the primary of V478 Cyg, in reasonable agreement with Weidner & Vink (2010). The different conclusions reached by us and Vink & Weidner regarding the discrepancy between evolutionary and dynamical masses is thus not clear. A dedicated analysis of their sample with our method is required to shed more light on this issue.
The right panel of Fig. 8 shows that the uncertainty on the spectroscopic masses are too large to allow any meaningful conclusion regarding the comparison between evolutionary and spectroscopic masses. This is partly due to the difficulty of the disentangling technique to correctly separate the Balmer lines of the systems’ components, and thus to the difficulty of determining surface gravities.
4.5 Periods and synchronization
Table 8 summarizes the rotational and orbital periods of each star/system. All systems but DH Cep being eclipsing, and assuming that the rotational axis is parallel to the orbital axis, the rotational period can be obtained from sin and the radius obtained from the fit of the light curve. In two systems (Y Cyg and AH Cep) the rotational periods are shorter than the orbital one. Since tides tend to synchronize the rotation of stars with the orbital period (e.g. Zahn 1977), one can conclude that in these two systems the stars have not finished yet to spin-down and that they were initially rotating with equatorial velocities larger then 150 km s. Interestingly, these systems are also the closest to the ZAMS (together with DH Cep for which no rotational period can be determined).
For V382 Cyg both periods are the same within the error bars, indicating that synchronization has been reached. This is also partly true for V478 Cyg.
XZ Cep is a peculiar case since the primary has P P while the opposite is true for the secondary. Given that this is the most evolved system of our sample and that the secondary fills its Roche lobe, this could indicate that mass and angular momentum transfer has lead to spin-up of the mass gainer and spin-down of the mass donor.
A better investigation of the effect of stellar evolution and mass transfer on rotational and orbital periods requires dedicated models. Such a study will be presented in a subsequent paper.
In this paper we have studied six Galactic massive binaries. We have obtained high-resolution spectroscopy at different orbital phases. A spectral disentangling technique was applied to yield the individual spectra of each component. These spectra were subsequently analyzed with atmosphere models and synthetic spectra to provide the stellar parameters and CNO surface abundances. The results are:
Five of the six systems have their components located on the main sequence (as defined by single-star evolutionary tracks). For the sixth system, one component is evolved.
Evolutionary masses determined using the single star evolutionary tracks of Ekström et al. (2012) are on average 16.5% higher than dynamical masses.
Three systems are detached (Y Cyg, V478 Cyg and AH Cep), two have at least one star filling its Roche lobe (V382 Cyg and XZ Cep) and no information is available for DH Cep since it is not eclipsing.
The CNO surface abundances of the three detached systems are little affected by stellar evolution and binarity. They are consistent with abundances of single stars within the uncertainties. Similar conclusions apply to the components of the contact system V382 Cyg and to the primary of XZ Cep. The N/C ratio of the secondary of XZ Cep is higher than what is expected for single stars with similar initial masses (according to the evolutionary tracks of Ekström et al. 2012).
Comparison of the derived surface abundances to those of systems known to have experienced mass transfer suggests the following trend: tides do not significantly affect the surface abundances of detached systems, while mass transfer, through the removal of external layers, can make chemically processed material appear at the surface. It may take some time after the onset of mass transfer before such effects are visible.
Some systems have reached synchronization, while other may still be in a spin-down phase. In the contact binary XZ Cep mass transfer has probably affected the rotational periods.
The trends emerging from our results, especially regarding the effect of binarity on surface abundances, need to be confirmed by analysis of additional detached and post-mass transfer systems.
We thank an anonymous referee for a prompt and positive report. We thank the Observatoire de Haute Provence for assistance with the T-193 telescope and the SOPHIE instrument. We warmly thank John Hillier for making the code CMFGEN available to the community and for constant help with it. L.M. acknowledges support from the Fonds National de la Recherche Scientifique (F.R.S.-F.N.R.S.) and through the ARC grant for Concerted Research Actions, financed by the French Community of Belgium (Wallonia-Brussels Federation).
- Bell et al. (1986) Bell, S. A., Hilditch, R. W., & Adamson, A. J. 1986, MNRAS, 223, 513
- Bouchy et al. (2013) Bouchy, F., Díaz, R. F., Hébrard, G., et al. 2013, A&A, 549, A49
- Bouchy et al. (2009) Bouchy, F., Hébrard, G., Udry, S., et al. 2009, A&A, 505, 853
- Burkholder et al. (1997) Burkholder, V., Massey, P., & Morrell, N. 1997, ApJ, 490, 328
- Chiosi & Maeder (1986) Chiosi, C. & Maeder, A. 1986, ARA&A, 24, 329
- Conti et al. (1977) Conti, P. S., Leep, E. M., & Lorre, J. J. 1977, ApJ, 214, 759
- de Mink et al. (2009) de Mink, S. E., Cantiello, M., Langer, N., et al. 2009, A&A, 497, 243
- Eggleton (1983) Eggleton, P. P. 1983, ApJ, 268, 368
- Ekström et al. (2012) Ekström, S., Georgy, C., Eggenberger, P., et al. 2012, A&A, 537, A146
- Garcia et al. (2014) Garcia, M., Herrero, A., Najarro, F., Lennon, D. J., & Alejandro Urbaneja, M. 2014, ApJ, 788, 64
- Garland et al. (2017) Garland, R., Dufton, P. L., Evans, C. J., et al. 2017, ArXiv e-prints [\eprint[arXiv]1704.07131]
- Gies (2003) Gies, D. R. 2003, in IAU Symposium, Vol. 212, A Massive Star Odyssey: From Main Sequence to Supernova, ed. K. van der Hucht, A. Herrero, & C. Esteban, 91
- Gray (1976) Gray, D. F. 1976, The observation and analysis of stellar photospheres
- Grevesse et al. (2010) Grevesse, N., Asplund, M., Sauval, A. J., & Scott, P. 2010, Ap&SS, 328, 179
- Grin et al. (2017) Grin, N. J., Ramírez-Agudelo, O. H., de Koter, A., et al. 2017, A&A, 600, A82
- Grunhut et al. (2013) Grunhut, J. H., Wade, G. A., Leutenegger, M., et al. 2013, MNRAS, 428, 1686
- Grunhut et al. (2017) Grunhut, J. H., Wade, G. A., Neiner, C., et al. 2017, MNRAS, 465, 2432
- Hadrava (1995) Hadrava, P. 1995, A&AS, 114, 393
- Harmanec et al. (2014) Harmanec, P., Holmgren, D. E., Wolf, M., et al. 2014, A&A, 563, A120
- Harries et al. (1997) Harries, T. J., Hilditch, R. W., & Hill, G. 1997, MNRAS, 285, 277
- Hilditch et al. (1996) Hilditch, R. W., Harries, T. J., & Bell, S. A. 1996, A&A, 314, 165
- Hillier & Miller (1998) Hillier, D. J. & Miller, D. L. 1998, ApJ, 496, 407
- Huang et al. (2010) Huang, W., Gies, D. R., & McSwain, M. V. 2010, ApJ, 722, 605
- Hunter et al. (2008) Hunter, I., Brott, I., Lennon, D. J., et al. 2008, ApJ, 676, L29
- Kobulnicky et al. (2014) Kobulnicky, H. A., Kiminki, D. C., Lundquist, M. J., et al. 2014, ApJS, 213, 34
- Langer et al. (2008) Langer, N., Cantiello, M., Yoon, S.-C., et al. 2008, in IAU Symposium, Vol. 250, Massive Stars as Cosmic Engines, ed. F. Bresolin, P. A. Crowther, & J. Puls, 167–178
- Linder et al. (2008) Linder, N., Rauw, G., Martins, F., et al. 2008, A&A, 489, 713
- Maeder & Meynet (2000) Maeder, A. & Meynet, G. 2000, ARA&A, 38, 143
- Mahy et al. (2017) Mahy, L., Damerdji, Y., Gosset, E., et al. 2017, ArXiv e-prints [\eprint[arXiv]1707.02060]
- Mahy et al. (2011) Mahy, L., Martins, F., Machado, C., Donati, J.-F., & Bouret, J.-C. 2011, A&A, 533, A9
- Mahy et al. (2013) Mahy, L., Rauw, G., De Becker, M., Eenens, P., & Flores, C. A. 2013, A&A, 550, A27
- Mandel & de Mink (2016) Mandel, I. & de Mink, S. E. 2016, MNRAS, 458, 2634
- Martins (2011) Martins, F. 2011, Bulletin de la Societe Royale des Sciences de Liege, 80, 29
- Martins et al. (2015) Martins, F., Hervé, A., Bouret, J.-C., et al. 2015, A&A, 575, A34
- Martins et al. (2005) Martins, F., Schaerer, D., & Hillier, D. J. 2005, A&A, 436, 1049
- Mathis & Remus (2013) Mathis, S. & Remus, F. 2013, in Lecture Notes in Physics, Berlin Springer Verlag, Vol. 857, Lecture Notes in Physics, Berlin Springer Verlag, ed. J.-P. Rozelot & C. . Neiner, 111–147
- Meynet et al. (2011) Meynet, G., Eggenberger, P., & Maeder, A. 2011, A&A, 525, L11
- Meynet & Maeder (2003) Meynet, G. & Maeder, A. 2003, A&A, 404, 975
- Nieva & Przybilla (2014) Nieva, M.-F. & Przybilla, N. 2014, A&A, 566, A7
- Palate et al. (2013) Palate, M., Rauw, G., Koenigsberger, G., & Moreno, E. 2013, A&A, 552, A39
- Pavlovski & Hensberge (2010) Pavlovski, K. & Hensberge, H. 2010, in Astronomical Society of the Pacific Conference Series, Vol. 435, Binaries - Key to Comprehension of the Universe, ed. A. Prša & M. Zejda, 207
- Penny & Gies (2009) Penny, L. R. & Gies, D. R. 2009, ApJ, 700, 844
- Perryman et al. (1997) Perryman, M. A. C., Lindegren, L., Kovalevsky, J., et al. 1997, A&A, 323, L49
- Popper & Hill (1991) Popper, D. M. & Hill, G. 1991, AJ, 101, 600
- Prša & Zwitter (2005) Prša, A. & Zwitter, T. 2005, ApJ, 628, 426
- Ramírez-Agudelo et al. (2013) Ramírez-Agudelo, O. H., Simón-Díaz, S., Sana, H., et al. 2013, A&A, 560, A29
- Raucq et al. (2017) Raucq, F., Gosset, E., Rauw, G., et al. 2017, A&A, 601, A133
- Raucq et al. (2016) Raucq, F., Rauw, G., Gosset, E., et al. 2016, A&A, 588, A10
- Rauw et al. (2000) Rauw, G., Sana, H., Gosset, E., et al. 2000, A&A, 360, 1003
- Sana et al. (2013) Sana, H., de Koter, A., de Mink, S. E., et al. 2013, A&A, 550, A107
- Sana et al. (2012) Sana, H., de Mink, S. E., de Koter, A., et al. 2012, Science, 337, 444
- Sana et al. (2006) Sana, H., Gosset, E., & Rauw, G. 2006, MNRAS, 371, 67
- Sana et al. (2014) Sana, H., Le Bouquin, J.-B., Lacour, S., et al. 2014, ApJS, 215, 15
- Simon et al. (1994) Simon, K. P., Sturm, E., & Fiedler, A. 1994, A&A, 292, 507
- Simón-Díaz & Herrero (2007) Simón-Díaz, S. & Herrero, A. 2007, A&A, 468, 1063
- Sturm & Simon (1994) Sturm, E. & Simon, K. P. 1994, A&A, 282, 93
- Ud-Doula et al. (2009) Ud-Doula, A., Owocki, S. P., & Townsend, R. H. D. 2009, MNRAS, 392, 1022
- Vink et al. (2001) Vink, J. S., de Koter, A., & Lamers, H. J. G. L. M. 2001, A&A, 369, 574
- Weidner & Vink (2010) Weidner, C. & Vink, J. S. 2010, A&A, 524, A98
- Wellstein et al. (2001) Wellstein, S., Langer, N., & Braun, H. 2001, A&A, 369, 939
- Wilson (1990) Wilson, R. E. 1990, ApJ, 356, 613
- Wilson & Devinney (1971) Wilson, R. E. & Devinney, E. J. 1971, ApJ, 166, 605
- Wolfe et al. (1967) Wolfe, Jr., R. H., Horak, H. G., & Storer, N. W. 1967, The machine computation of spectroscopic binary elements, ed. O. Struve & M. Hack, 251
- Zahn (1977) Zahn, J.-P. 1977, A&A, 57, 383
- Zahn (1989) Zahn, J.-P. 1989, A&A, 220, 112