286 Acta Polytechnica CTU Proceedings 2(1): 286–290, 2015 286 doi: 10.14311/APP.2015.02.0286 Rotation of the Mass Donors in High-mass X-ray Binaries and Symbiotic Stars K. Stoyanov1, R. Zamanov1 1Institute of Astronomy and National Astronomical Observatory, Bulgarian Academy of Sciences, 72 Tsarigradsko Shousse Blvd., 1784 Sofia, Bulgaria Corresponding author: kstoyanov@astro.bas.bg Abstract Our aim is to investigate the tidal interaction in High-mass X-ray Binaries and Symbiotic stars in order to determine in which objects the rotation of the mass donors is synchronized or pseudosynchronized with the orbital motion of the compact companion. We find that the Be/X-ray binaries are not synchronized and the orbital periods of the systems are greater than the rotational periods of the mass donors. The giant and supergiant High-mass X-ray binaries and symbiotic stars are close to synchronization. We compare the rotation of mass donors in symbiotics with the projected rotational velocities of field giants and find that the M giants in S-type symbiotics rotate on average 1.5 times faster than the field M giants. We find that the projected rotational velocity of the red giant in symbiotic star MWC 560 is v sin i= 8.2 ± 1.5 km s−1, and estimate its rotational period to be Prot = 144 - 306 days. Using the theoretical predictions of tidal interaction and pseudosynchronization, we estimate the orbital eccentricity e = 0.68 − 0.82. Keywords: stars: binaries: symbiotic - stars: rotation - stars: late type - stars: binaries: close - X-rays: binaries. 1 Introduction A High-mass X-ray binary system consists of a com- pact object (a neutron star or a black hole) accreting material from an O or B companion star. They are di- vided into Be/X-ray binaries (main-sequence star as a companion) and Giant/Supergiant X-ray binaries (gi- ant or supergiant star as a companion). Accretion of matter is different for both types of X-ray binaries. In the Be/X-ray binaries, the compact object crosses the circumstellar disc and accretes matter from that disk. In the giant/supergiant X-ray binaries, the mass donor ejects a slow and dense wind radially outflowing from the equator and the compact object directly accretes the stellar wind through Bondi-Hoyle-Lyttleton accre- tion. Symbiotic stars are interacting binaries, consisting of an evolved giant (either a normal red giant in S-types symbiotics or a Mira-type variable in D-types symbi- otics) transferring mass to a hot and luminous white dwarf or neutron star. The symbiotic stars are sur- rounded by a rich and luminous nebula resulting from the presence of both an evolved giant with a heavy mass-loss and of a hot companion abundant in ioniz- ing photons and often emanating its own wind. 2 Synchronization and Pseudosynchronization In a binary with a circular orbit the rotational period of the primary, Prot, reaches an equilibrium value at the orbital period, Porb = Prot. In other words the synchronous rotation (synchronization) means that the rotational period is equal to the orbital period. In a binary with an eccentric orbit, the corresponding equi- librium is reached at a value of Prot which is less than Porb, the amount less being a function solely of the or- bital eccentricity e. In practice, in a binary with an eccentric orbit the tidal force acts to synchronize the rotation of the mass donor with the motion of the com- pact object at the periastron - pseudosynchronous rota- tion. To calculate the period of pseudosynchronization, Pps, we use (Hut 1981): Pps = (1 + 3e2 + 3 8 e4)(1 − e2) 3 2 1 + 15 2 e2 + 45 8 e4 + 5 16 e6 Porb. (1) 2.1 Stars with radiative envelopes Following Hurley, Tout & Pols (2002) the circulariza- tion timescale for stars with radiative envelopes is: 1 τcirc = 21 2 ( GM1 R31 )1 2 q2 (1 + q2) 11 6 E2 ( R1 a )21 2 , (2) 286 http://dx.doi.org/10.14311/APP.2015.02.0286 Rotation of the Mass Donors in High-mass X-ray Binaries and Symbiotic Stars where M1 and R1 are the mass and the radius of the pri- mary respectively, q2 is the mass ratio M2/M1, and a is the semi-major axis. The second-order tidal coefficient E2 = 1.592 ×10−9M2.841 . The synchronization time scale is given as, τsync = Kτcirc, (3) where K is: K ≈ 0.015 rg 1 + q2 q2 ( R1 a )2 . (4) For the gyration radius of the primary rg we adopt rg ≈ 0.16 for giants, and rg ≈ 0.25 for main sequence stars (Claret & Gimenez, 1989). 2.2 Stars with convective envelopes Following Hurley, Tout & Pols (2002) the synchroniza- tion timescale for stars with convective envelopes is: τsyn ≈ 800 ( M1R1 L1 )1/3 M21 ( M1M2 + 1)2 R61 P4orb yr, (5) where L1 is the luminosity of the giant. The circularization time scale is: 1 τcirc = 21 2 ( k T ) q2(1 + q2) ( R1 a )8 . (6) In Eq. 6, (k/T) is:( k T ) = 2 21 fconv τconv Menv M1 yr−1, (7) where Renv is the depth of the convective envelope, Menv is the envelope’s mass, and τconv = 0.4311 ( MenvRenv(R1 − 12Renv) 3L1 )1 3 yr (8) is the eddy turnover time scale (the time scale on which the largest convective cells turnover). The numerical factor fconv is fconv = min [ 1, ( Ptid 2τconv )2] , (9) where Ptid is the tidal pumping time scale given by 1 Ptid = ∣∣∣∣ 1Porb − 1Prot ∣∣∣∣ . (10) The pseudosynchronization timescale is τps = (7/3(α - 3)) τcirc, where α is a dimensionless quantity, representing the ratio of the orbital and rotational an- gular momentum: α = q2 1 + q2 1 r2g ( a R1 )2 . (11) For a red giant we adopt rg ≈ 0.3 (Claret, 2007). In all equations, the masses, the radii and the lu- munosities are in solar units. 3 High-mass X-ray Binaries The orbital and stellar parameters of 13 High-mass X- ray binaries are given in Table 1 and Table 2 in Stoy- anov & Zamanov (2009). We add 2 more objects - 4U 2206+54 and MWC 148. The orbital and stellar parameters are taken from Ribó et. al. (2006) and Casares et al. (2012) respectively. Using Eq.2 and Eq.3 we estimate the circularization and synchronization timescales. The results are given in Table 1. The lifetime of a star on the main sequence can be estimated as τMS = 10 10(M�/M) 2.5 years (Hansen & Kawaler, 1994). Comparing these lifetimes with τsync from Table 1, we see that among the Be/X-ray binaries only for LSI+610303 is τsync ∼ τMS. This is the only Be/X-ray binary for which we can expect considerable changes of the rotation of the primary during the life- time of the Be star. The lifetime of the giant is comparable or longer then τcirc and τsync for the giant/supergiant sys- tems with short orbital periods. The exceptions are V725 Tau and BP Cru, for which τsync and τcirc are longer than the lifetime of the giant/supergiant stage. On Fig.1 in Stoyanov & Zamanov (2009) is plotted Prot versus Pps. The giant/supergiant systems are lo- cated close to the line Pps = Prot, while those with mass donors from spectral class V are far away from the equilibrium. In the Be/X-ray systems BQ Cam, V635 Cas, V725 Tau and 4U 2206+54, the tidal force spinning down the donor star. For the system LSI+610303, the rotation of the mass donor is close to pseudosynchronization. This is the only Be/X-ray binary in which τsync is compa- rable with the life-time of the binary. In the binaries X Per and MWC 148, the neutron star is far away from the Be star and the tidal force is weak. Giant and supergiant systems are close to (pseudo)synchronization. In these binaries the rota- tion of the mass donors is influenced by the presence of the compact object. In LMC X-4 and Cen X-3, the mass donors are synchronized and the orbits are circu- larized. With respect to the rotation of the mass donor, V725 Tau is similar to the Be/X-ray binaries. Cyg X-1 is synchronized and almost circularized. V830 Cen is pseudosynchronized but not circularized yet. The sys- tems LSI+650010 and Vela X-1 are close to pseudosyn- chronization and the tidal force accelerates the rotation of the mass donors. In the case of SMC X-1, the tidal force acts as a decelerator of the rotation of the mass donor. In BP Cru, a gas stream from the mass donor 287 K. Stoyanov, R. Zamanov exists, probably resulting from the strong tidal force and spin-up of the mass donor (Leahy & Kostka, 2008). Table 1: Calculated time scales. Given here are the name of the object, synchronization time scale, circu- larization time scale, and the lifetime. object τsync [yr] τcirc [yr] lifetime [yr] Be/X-ray binaries LSI+610303* 2.8×107 2.4×108 5.6×107 X Per 6.2×1017 1.8×1021 1.1×107 BQ Cam 3.5×1011 7.6×1013 3.9×106 V635 Cas 1.4×1011 9.5×1012 7.3×106 4U 2206+54 4.9×109 3.7×1011 7.3×107 MWC 148 1.2×1017 5.1×1020 9.8×107 Giant systems V725 Tau 2.8×1012 8.0×1014 4×105 LMC X-4 4.5×102 7.7×102 1×106 Cen X-3 2.3×103 4.2×103 5×105 Supergiant systems V830 Cen 7.5×103 1.4×104 1×106 LSI+650010 1.3×104 3.9×104 1×106 Vela X-1 1.0×104 2.8×104 3.9×105 SMC X-1 3.3×104 8.2×104 8.8×105 BP Cru 1.8×106 8.8×106 8×104 Cyg X-1 < 1 < 1 1×105 *assuming neutron star as a secondary component. 4 S-type Symbiotic Stars 43 symbiotic stars have been observed with FEROS spectrograph at the 2.2m ESO telescope of the La Silla Observatory (Zamanov et al. 2007). The data for the rotation of 55 field red giants are taken from the lit- erature. M giants in S-type symbiotics rotate faster than the field M giants. Histograms of the available v sin i data for the red giants are plotted in Fig.2 in Za- manov & Stoyanov (2012). For the field M0III-M6III giants we calculate a mean v sin i=5.0 km s−1, median v sin i=4.3 km s−1, and standard deviation of the mean σ = 4.0 km s−1. For the M0III-M6III giants in sym- biotics, we get a mean v sin i= 7.8 km s−1, median v sin i=8.0 km s−1, and standard deviation of the mean σ =2.1 km s−1. There are 5 objects in our sample that deviate from the synchronization. These objects are RS Oph, MWC 560, CH Cyg, CD-43◦14304 and Z And. In three of them collimated jets are detected: Z And (Skopal et al. 2009); CH Cyg (Crocker et al. 2002), MWC 560 (Tomov et al. 1990). Additionally to the jets, ejec- tion of blobs are detected from RS Oph and CH Cyg (Iijima et al. 1994). This confirms the suggestions that in the jet-ejecting symbiotics the mass donors rotate faster than the orbital periods. Probably there is a link between the jets and the mass donor rotation. On Fig.1 are plotted together the High-mass X-ray binaries and the S-type symbiotic stars. It shows that none of the objects in our sample is above the line of synchronization. Figure 1: Prot versus Porb on a logarithmic scale. The blue circles indicates the High-mass X-ray Binaries - filled symbols for giant/supergiant systems and open symbols for Be/X-ray binaries. The red triangles indi- cates the S-type symbiotic stars. 5 Orbital Eccentricity of MWC 560 MWC 560 is a symbiotic star, which consists of a red giant and a white dwarf (Tomov et al. 1990). The most spectacular features of this object are the col- limated ejections of matter with velocities of up to ∼ 6000 km s−1 (Tomov et al. 1992) and the resemblance of its emission line spectrum to that of the low-redshift quasars (Zamanov & Marziani, 2002). The jet ejections are along the line of sight and the system is seen almost pole-on (i < 16◦). This makes it difficult to obtain the orbital eccentricity of the system in a conventional way. For the system we adopt Rg = 140 ± 7 R�, Lg ∼ 2400 L�, Mg = 1.7 M�, Mwd = 0.65 M�, Porb = 1931 ± 162 day (Gromadzki et al. 2007), Renv = 0.9 Rg 288 Rotation of the Mass Donors in High-mass X-ray Binaries and Symbiotic Stars and Menv = 1.0 M� (Herwig 2005). With the above values of the parameters assumed, we derive the semi- major axis of the orbit to be a ≈ 860 R�. Using these parameters, we calculate from Eq. 5 and Eq. 6 the syn- chronization and circularization time scales: τsync = 2.6 × 104 yr and τcirc = 3.1 × 106 yr. The typical life- time of a symbiotic star is τss ∼ 105 yr (Yungelson et al., 1995). From the rate of accretion on the white dwarf, Ṁacc ≈ 5 × 10−7 M� (Schmid et al. 2001), we can esti- mate, that it will take 106 yr to accrete ∼ 0.5 M� from the envelope of the red giant companion. Because the giant also losses mass via stellar wind, we find that the lifetime of the symbiotic phase of MWC 560 should be τss ≤ 106 yr. For MWC 560 we have therefore the situ- ation in which τps < τsyn < τss < τcirc. This means that the symbiotic phase is long enough that the tidal forces can (pseudo)synchronize the rotation of the red giant. On the other hand, the value of τcirc demonstrates that the symbiotic lifetime of MWC 560 is shorter than the circularization time, and therefore the orbit can be ec- centric. This is in agreement with the observational evidences found by Fekel et al. (2007) that the sym- biotic stars with Porb > 800 days tend to have eccen- tric orbits. The above implies that in MWC 560, the red giant is probably synchronized, but the orbit is not circularized. To determine the orbital eccentricity of MWC 560, we need to calculate Prot for the mass donor. We analyzed 21 high resolution spectra of MWC 560 and obtained value for v sin i= 8.2 ± 1.5 km s−1. Us- ing Rg = 140 ± 7 R� and i = 12◦ − 16◦, we calcu- late Prot = 144 − 306 days. This value is less than the orbital period. MWC 560 should be close to syn- chronization or pseudosynchronization, and Prot = Pps. Using Eq. 1 we can estimate the orbital eccentricity to be e = 0.68 − 0.82. 6 Conclusions Using rotational velocity measurements and the theory of synchronization/pseudosynchronization we: (1) find that the Be/X-ray binaries are far away from (pseudo)synchronization. The tidal force in the Be/X-ray binaries acts as a decelerator of the rotation of the mass donors. The only Be/X-ray binary which is close to pseudosynchronization is the LSI+610303. The objects containing mass donors of spectral class I and III typically have Prot ∼ Pps and are close to (pseudo)synchronization; (2) demonstrate that the M giants in symbiotic stars rotate faster than the field giants. Most symbiotics with orbital period less than 1000 d are synchronized; (3) show that the High-mass X-ray binaries and the S-type symbiotic stars are either on the line of synchro- nization or they are under the line. None of the objects in our sample is above the line of synchronization. (4) calculate that the orbit of the symbiotic star MWC 560 should be highly eccentric, with e ∼ 0.7. Acknowledgement We thank the anonymous referee for constructive com- ments. This work was supported by the OP “HRD“, ESF and Bulgarian Ministry of Education, Youth and Science under the contract BG051PO001-3.3.06-0047. References [1] Casares, J., Ribó, M., Ribas, I., Paredes, J. M., Vi- lardell, F., & Negueruela, I.: 2012, MNRAS, 421, 1103 doi:10.1111/j.1365-2966.2011.20368.x [2] Claret, A.: 2007, A&A 467, 1389 [3] Claret, A., & Gimenez, A.: 1989, A&AS, 81, 37 [4] Crocker, M. M., Davis, R. J., Spencer, R. E., et al.: 2002, MNRAS, 335, 1100 doi:10.1046/j.1365-8711.2002.05705.x [5] Fekel, F. C., Hinkle, K. H., Joyce, R. R., Wood, P. R., Lebzelter, T.: 2007, AJ 133, 17 [6] Gromadzki, M., Miko lajewska, J., Whitelock, P. A., Marang, F.: 2007, A&A 463, 703 [7] Hansen, C. J., & Kawaler, S. D.: 1994, Stellar In- teriors. Physical Principles, Structure, and Evo- lution, Springer-Verlag [8] Herwig, F.: 2005, ARA&A 43, 435 doi:10.1146/annurev.astro.43.072103.150600 [9] Hurley, J. R., Tout, C. A., & Pols, O. R.: 2002, MNRAS, 329, 897 doi:10.1046/j.1365-8711.2002.05038.x [10] Hut, P.: 1981, A&A, 99, 126 [11] Iijima, T., Strafella, F., Sabbadin, F., & Bianchini, A.: 1994, A&A, 283, 919 [12] Leahy, D. A., & Kostka, M.: 2008, MNRAS, 384, 747 doi:10.1111/j.1365-2966.2007.12754.x [13] Ribó, M., Negueruela, I., Blay, P., Torrejón, J. M., & Reig, P.: 2006, A&A, 449, 687 [14] Schmid, H. M., Kaufer, A., Camenzind, M., Riv- inius, T., Stahl, O., Szeifert, T., Tubbesing, S., Wolf, B.: 2001, A&A 377, 206 [15] Skopal, A., Pribulla, T., Budaj, J., et al.: 2009, ApJ, 690, 1222 doi:10.1088/0004-637X/690/2/1222 [16] Stoyanov, K. A., & Zamanov, R. K.: 2009, As- tronomische Nachrichten, 330, 727 289 http://dx.doi.org/10.1111/j.1365-2966.2011.20368.x http://dx.doi.org/10.1046/j.1365-8711.2002.05705.x http://dx.doi.org/10.1146/annurev.astro.43.072103.150600 http://dx.doi.org/10.1046/j.1365-8711.2002.05038.x http://dx.doi.org/10.1111/j.1365-2966.2007.12754.x http://dx.doi.org/10.1088/0004-637X/690/2/1222 K. Stoyanov, R. Zamanov [17] Tomov, T., Zamanov, R., Kolev, D., Georgiev, L., Antov, A., Mikolajewski, M., Esipov, V.: 1992, MNRAS 258, 23 doi:10.1093/mnras/258.1.23 [18] Tomov, T., Kolev, D., Zamanov, R., Georgiev, L., Antov, A.: 1990, Nature 346, 637 doi:10.1038/346637a0 [19] Yungelson, L., Livio, M., Tutukov, A., & Kenyon, S. J.: 1995, ApJ 447, 656 doi:10.1086/175908 [20] Zamanov, R. K., Bode, M. F., Melo, C. H. F., et al.: 2007, MNRAS, 380, 1053 doi:10.1111/j.1365-2966.2007.12150.x [21] Zamanov, R., Marziani, P.: 2002, ApJL 571, L77 doi:10.1086/341367 [22] Zamanov, R. K., & Stoyanov, K. A.: 2012, Bulgar- ian Astronomical Journal, 18, 41 290 http://dx.doi.org/10.1093/mnras/258.1.23 http://dx.doi.org/10.1038/346637a0 http://dx.doi.org/10.1086/175908 http://dx.doi.org/10.1111/j.1365-2966.2007.12150.x http://dx.doi.org/10.1086/341367 Introduction Synchronization and Pseudosynchronization Stars with radiative envelopes Stars with convective envelopes High-mass X-ray Binaries S-type Symbiotic Stars Orbital Eccentricity of MWC 560 Conclusions