Biperiodic orbital period modulation of the W UMa binary system CK Bootis

Abstract

A study of the orbital period variation of the W UMa system CK Bootis is made using an extended observational time base. The biperiodicity of the orbital period modulation is emphasized. Both detected periodicities (24.14 yr and 10.62 yr) cannot be explained through the light-time effect unless the companion would be a white dwarf as suggested by other authors, too.

Moreover, we also argue that, nowadays at least, it seems that there is no causal relation between the orbital period variation and the recently discovered visual companion.

Consequently, we infer that at least one of the two periodicities may be related to the magnetic activity cycles in the component stars of CK Boo, while the other periodicity could be related to the presence of a fourth companion in the system

Authors

C. Vamoș
Tiberiu Popoviciu Institute of Numerical Analysis, Romanian Academy

A. Pop
Astronomical Institute, Astronomical Observatory Cluj-Napoca

Keywords

Stars; binaries; eclipsing Stars; binaries; close Stars; individual; CK Boo Methods; data analysis

Cite this paper as:

A. Pop, C. Vamoş, Biperiodic orbital period modulation of the W UMa binary system CK Bootis, New Astronomy, 17 (2012), 667-672, doi: 10.1016/j.newast.2012.04.003

PDF

About this paper

Print ISSN

1384-1076

Online ISSN

MR

?

ZBL

?

[1] Abhyankar, K.D., Panchatsaram, T., 1984. Light time effect in TW Draconis. MNRAS 211, 75.

[2] Applegate, J.H., 1992. A mechanism for orbital period modulation in close binaries. ApJ 385, 621.

[3] Aslan, Z., Derman, E., 1986. Photoelectric observations of CK bootis. A&AS 66, 281.

[4] Barning, F.J.M., 1963. The numerical analysis of the light-curve of 12 Lacertae. BAN 17, 22.

[5] Borkovits, T., Hegedüs, T., 1996. On the invisible components of some eclipsing binaries. A&AS 12, 63.

[6] Borkovits, T., Csizmadia, Sz., Hegedüs, T., Bíró, I.B., Sándor, Zs., Opitz, A., 2002. Complex period variations in the binary system IM Aurigae. A&A 392, 895.

[7] Borkovits, T., Elkhateeb, M.M., Csizmadia, Sz., Nuspl, J., Biró, I.B., Hegedüs, T., Csorvási, R., 2005. Indirect evidence for short period magnetic cycles in W UMa stars. Period analysis of five overcontact systems. A&A 441, 1087.

[8] Breger, M., Huang, L., Jiang, S.-Y., Guo, Z.-H., Antonello, E., Mantegazza, L., 1987. Multiple close frequencies of the Delta Scuti star h2 Tau. A&A, 175, 117.

[9] D’Angelo, C., van Kerkwijk, M.H., Rucinski, S.M., 2006. Contact binaries with additional components. II. A spectroscopic search for faint tertiaries. AJ 132, 650.

[10] Demircan, O., 1987. BV observations of W UMa-type binaries: CK Bootis, BI CVn, and AH Vir. Ap&SS 135, 169.

[11] Demircan, O., Selam, S.O., 1993. Long-term behaviour of the orbital period of Algoltype binary ST Persei. A&AS 98, 513.

[12] Derman, E., Demircan, O., 1992. A period study of GK Cephei. AJ 103, 599.

[13] Diethelm, R., 2011.Timings of minima of eclipsing binaries. IBVS, 5992.

[14] Gazeas, K.D., Niarchos, P.G., Zola, S., Kreiner, J.M., Rucinski, S.M., 2006. Physical parameters of components in close binary systems: VI. AcA 56, 127.

[15] Hall, D.S., 1989. The relation between RS CVn and Algol. SSR 50, 219.

[16] Hall, D.S., 1990. Period changes and magnetic cycles. In:  Ibanoglu, C. (Ed.), Active Close Binaries. Kluwer Academic Publishers, Dordrecht, p. 95.

[17] Huang, S.-S., 1963. Modes of mass ejection by binary stars and effect on their orbital periods. ApJ 138, 471.

[18] Jia, G., Liu, X., Huang, H., 1992. 1991 Photometry of the W UMa type binary CK Boo. IBVS, 3727.

[19] Kalci, R., Derman, E., 2005. CK Bootis: a W UMa system with a small mass ratio. AN 326, 342.

[20] Kalimeris, A., Rovithis-Livaniou, H., Rovithis, P., 1994a. On the orbital changes in contact binaries. A&A 282, 775.

[21] Kalimeris, A., Rovithis-Livaniou, H., Rovithis, P., Oprescu, G., Dumitrescu, A., Suran, M.D., 1994b. An orbital period study of the contact system AB Andromedae. A&A 291, 765.

[22] Kopal, Z., 1978. Dynamics of Close Binary Systems. D. Reidel Publishing Company, Dordrecht. Kreiner, J.M., 1971. Investigation of changes in periods of eclipsing variables. AcA 21, 365.

[23] Krzesin´ ski, J., Mikolajewski, M., Pajdosz, G., Zola, S., 1991. An observational study of a spotted W UMa-type binary star CK Bootis. Ap&SS 184, 37.

[24] Lanza, A.F., Rodonò, M., 1999. Orbital period modulation and quadrupole moment changes in magnetically active close binaries. A&A 349, 887.

[25] Lanza, A.F., Rodonò, M., Rosner, R., 1998. Orbital period modulation and magnetic cycles in close binaries. MNRAS 296, 893.

[26] Li, Y., Qian, S.B., 2005. GMG 2.4-m telescope and light-time effect research for W UMa-type binaries. In: Sterken, C. (Ed.), ASP Conference Series, vol. 335, p. 245.

[27] Liao, W.-P., Qian, S.-B., 2010. The most plausible explanation of the orbital period changes in close binaries: the case of the RS CVn-type binary WW Dra. MNRAS 405, 1930.

[28] López de Coca, P., Garrido, R., Rolland, A., 1984. Mode identification for the Delta Scuti variable HR 7222. A&AS 58, 441.

[29] Panchatsaram, T., 1981. Light-time effect in RT Persei. Ap&SS 77, 179.

[30] Panchatsaram, T., Abhyankar, K.D., 1981. SW Lacertae – a quadruple system. Bull. Astron. Soc. India 9, 31.

[31] Pop, A., 1996. The use of characteristic time to the study of period variability phenomena. Romanian Astron. J. 6, 147.

[32] Pop, A., 1999. A multiperiodic ephemeris for RZ Cephei. In: IBVS, vol. 4801.

[33] Pop, A., 2005. On the orbital period modulation of the eclipsing binary system Y Leonis. In: Sterken, C. (Ed.), ASP Conference Series, vol. 335, p. 263.

[34] Pop, A., 2007. On the detection and diagnosis of low-level stellar variability. Romanian Astron. J. 17, 35.

[35] Pop, A., Roman, R., 2005. RT and: multiperiodic orbital period variability, resonance, and the hypothesis of the light-time effect. Romanian Astron. J. 15, 163, Supplement.

[36] Pop, A., Vamos, C., 2007. The detection and diagnosis of low-level stellar variability through resampling methods. Romanian Astron. J. 17, 101.

[37] Pop, A., Todoran, I., Agerer, F., 1996. Multiperiodic behaviour of RT Persei. Romanian Astron. J. 6, 141.

[38] Pop, A., Liteanu, V., Pavel, C., 2000. On the multiperiodicity of the orbital period variation of ST Persei. Romanian Astron. J. 10, 25.

[39] Pop, A., Liteanu, V., Moldovan, D., 2003. On the period variability of GP Andromedae. Ap&SS 284, 1207.

[40] Pop, A., Vamos, C., Turcu, V., 2010. Deterministic components in the light curve amplitude of Y Oph. AJ 139, 425.

[41] Pop, A., Turcu, V., Marcu, A., 2011. The intriguing orbital period variability of Y Leonis. Ap&SS 333, 17.

[42] Pribulla, T., Rucinski, S.M., 2006. Contact binaries with additional components. I. The extant data. AJ 131, 2986.

[43] Pringle, J.E., 1975. Period changes in eruptive binaries. MNRAS 170, 633.

[44] Qian, S., 2001. Possible mass and angular momentum loss in Algol-type binaries. V RT Persei and TX Ursae Majoris. AJ 122, 2686.

[45] Qian, S., Liu, Q., 2000. A possible connection between the variability of light curve and the change of the orbital period in the contact binary CK Bootis. Ap&SS 271, 331.

[46] Rafert, J.B., 1982. Periodic ephemerides for 49 eclipsing binary-star systems. PASP 94, 485.

[47] Rucinski, S., Lu, W., 1999. Radial velocity studies of close binary stars. II. AJ 118, 2451.

[48] Rucinski, S.M., Pribulla, T., van Kerkwijk, M.H., 2007. Contact binaries with additional components. III. A search using adaptive optics. AJ 134, 2353.

[49] Samolyk, G., 2011. Recent minima of 146 eclipsing binary stars. J. Amer. Assoc. Var. Star Obs. 39, 177.

[50] Sârbu, C., 1997. A new median – based robust regression method. Croat. Chem. Acta 70, 795.

[51] Szebehely, V., Zare, K., 1977. Stability of classical triplets and their hierarchy. A&A 58, 145.

[52] Wehlau, W., Leung, K.-C., 1964. The multiple periodicity of Delta Delphini. ApJ 139, 843.

[53] Wood, D.B., Forbes, J.E., 1963. Ephemerides of eclipsing stars. AJ 68, 257.

[54] Yüce, K., Selam, S.O., Albayrak, B., Ak, T., 2006.    . Ap&SS 304, 67.

Biperiodic orbital period modulation of the W UMa binary system CK Bootis Alexandru Pop a, , Ca˘lin Vamos ß b a Astronomical Institute of the Romanian Academy, Astronomical Observatory Cluj-Napoca, Str. Cires ßilor 19, RO-400487 Cluj-Napoca, Romania b ‘‘T. Popoviciu’’ Institute of Numerical Analysis, Romanian Academy, P.O. Box 68, 400110 Cluj-Napoca, Romania article info Article history: Received 17 December 2011 Received in revised form 25 March 2012 Accepted 11 April 2012 Available online 20 April 2012 Communicated by E.P.J. van den Heuvel Keywords: Stars: binaries: eclipsing Stars: binaries: close Stars: individual: CK Boo Methods: data analysis abstract A study of the orbital period variation of the W UMa system CK Bootis is made using an extended obser- vational time base. The biperiodicity of the orbital period modulation is emphasized. Both detected peri- odicities (24.14 yr and 10.62 yr) cannot be explained through the light-time effect unless the companion would be a white dwarf as suggested by other authors, too. Moreover, we also argue that, nowadays at least, it seems that there is no causal relation between the orbital period variation and the recently dis- covered visual companion. Consequently, we infer that at least one of the two periodicities may be related to the magnetic activity cycles in the component stars of CK Boo, while the other periodicity could be related to the presence of a fourth companion in the system. Ó 2012 Elsevier B.V. All rights reserved. 1. Introduction Usually, a first approach of the analysis of astronomical time series obtained through stellar variability studies (light and radial velocity curves), implicitly assumes that we are dealing with the superposition of a deterministic component (secular trend or (multi) periodic signal) and a random component (observational noise or possible correlated variability physically related to the observed object). Very often, the secular trend can be modelled by low order polynomials. The residuals obtained after removing the polynomial trend have to be further analysed in order to investigate the pres- ence of possible (multi)periodic variability phenomena. The meth- odology suitable for such studies has been developed alongside the techniques for analysis of the pulsating stars light curves (e.g. Bar- ning, 1963; Wehlau and Leung, 1964; López de Coca et al., 1984; Breger et al., 1987). This methodology relies on the successive pre- whitening in the time domain. In the case of period variability phenomena, we are dealing with OC curves relying on timing data of light curves extrema of the observed variable star. It is well-known that concurrently with the gradual increase of the observing time base, the shape of many variable stars OC curves became more and more intricate. Corre- spondingly, more complex ephemerides were needed to describe the distribution of the observed timing data. Thus, from low order polynomial ephemerides (e.g. Wood and Forbes, 1963), one evolved either to superposition of (multi)periodic term(s) on low order polynomial models (e.g. Kreiner, 1971; Panchatsaram and Abhyankar, 1981; Rafert, 1982; Abhyankar and Panchatsaram, 1984; Demircan and Selam, 1993; Pop, 1996; Pop et al., 1996, 2000, 2003; Borkovits et al., 2002, 2005; Li and Qian, 2005), or to high order polynomial ephemerides (e.g. Derman and Demircan, 1992; Kalimeris et al., 1994a,b). CK Bootis (BD + 9°2916, HD 128141, HIP 71319) is an A type W UMa binary system with a spectral type of F7/F8 V (Rucinski and Lu, 1999), and an orbital period of about 0.355 d. Examining the run of the OC curve of CK Boo, Aslan and Derman (1986) hypoth- esized that a period variation could occur in this binary system, but without assuming a given mathematical law. Demircan (1987) also noted some change in the run of the OC residuals, but he ex- pressed his doubt (like Aslan and Derman, 1986) concerning the reality of the corresponding period change. Jia et al. (1992) empha- sized for the first time a linear increase of the orbital period of CK Boo. Later, Qian and Liu (2000) confirmed the presence of a para- bolic trend in the OC residuals. Moreover, they found evidence for a small amplitude period modulation phenomenon superposed on this trend. They estimated the involved periodicity at about 14 yr. On the other hand, they also found that the luminosity of the secondary component varies in phase, approximately with the same period. They concluded that the observed orbital period modulation in CK Boo is related to cyclic magnetic activity of the secondary component through Applegate’s (1992) mechanism. Kalci and Derman (2005) confirmed the previous obtained results. The period of the orbital modulation was found to be about 15.75 yr. They also concluded that the amplitude of this period var- iation might increase with time and, for this reason, they rejected the hypothesis of a third body in the system CK Boo. Consequently, 1384-1076/$ - see front matter Ó 2012 Elsevier B.V. All rights reserved. http://dx.doi.org/10.1016/j.newast.2012.04.003 Corresponding author. Tel.: +40 264594592; fax: +40 264596889. E-mail addresses: andi_pop@yahoo.com (A. Pop), cvamos@ictp.acad.ro (C. Vamosß). New Astronomy 17 (2012) 667–672 Contents lists available at SciVerse ScienceDirect New Astronomy journal homepage: www.elsevier.com/locate/newast
taking into account the spectral type of this system, they adopted the mechanism of magnetic activity cycles occurring in both com- ponent stars as a cause for the observed orbital period variation. A year later, Yüce et al. (2006) resumed the study of orbital period variation of CK Boo. They supplied a new confirmation of the par- abolic trend and of the periodic modulation in the OC curve. They assumed that the last is caused by the presence of a third body in the system, with a mass of 0.48 M Sun , revolving on an eccentric or- bit (e = 0.55) with a period of 17.41 yr. Concerning the nature of this companion, more massive than the secondary component, Yüce et al. (2006) considered that it could be either a compact ob- ject, or a close binary. The light-time effect in case of CK Boo was also approached by Pribulla and Rucinski (2006). The involved periodicity – considered by the authors to be very uncertain – is 21.3 yr, while the corresponding orbit was found to be highly eccentric e = 0.96. The minimum mass of the hypothetical compan- ion was estimated to be extremely high with respect to that of the secondary star of the system, i.e., 1.50 M Sun . Taking into account the small UV excess detected by Krzesin´ ski et al. (1991), Pribulla and Rucinski (2006) concluded that it could be a white dwarf. Valuable complementary information has been supplied by Pri- bulla and Rucinski (2006), D’Angelo et al. (2006), and Rucinski et al. (2007) through both spectroscopic and adaptive optics observa- tions. They found that CK Boo has a close visual companion at a separation of 0.12 arcsec, with a late spectral type M0/M1 V. According to Rucinski et al. (2007) its mass might be 0.45 M Sun . They estimated the projected separation between the close binary and companion at 19 AU, and the orbital period of the wide system at 54 yr. From methodological viewpoint, we have to note that both Qian and Liu’s (2000) and Kalci and Derman’s (2005) approaches, which emphasized the presence of a wavelike component in the run of the OC residuals, relies on the method of Kalimeris et al. (1994a). This makes this binary system an interesting target for a study using a complementary approach of (multi)periodicity anal- ysis namely that described in the papers of Pop (1996, 1999, 2005) and Pop et al. (2000, 2003, 2011) (see also Borkovits et al., 2005). We reanalysed the period variability of the W UMa binary system CK Boo, relying on an extended time base: from DT 29 yr (Kalci and Derman, 2005), to DT 36 yr. This enabled us to emphasize the biperiodic character of the orbital period variation, superposed on its linear increase, previously detected by Qian and Liu (2000), and Kalci and Derman (2005). The interpretation of these results and their astrophysical implications are investigated. 2. Observational data For the present study of the contact binary CK Boo, we collected 70 times of primary minimum light and 52 times of secondary minimum light from the Lichtenknecker – Database of the BAV (www.bav-astro.de/LKDB) and ‘‘OC Gateway Database of times of minima (E) and maxima (RR)’’ of Variable Star Section of Czech Astronomical Society (http://var.astro.cz/ocgate/index.php), as well as from the papers of Kalci and Derman (2005), Diethelm (2011), and Samolyk (2011). All timing data have been checked according to their original sources. The observed minima corre- sponding to the same event (eclipse) have been averaged. The final data set containing both primary and secondary minima consists in 101 data points which cover a time base of about 36 yr. 3. Analysis methods The methodology used for the analysis and modelling of the timing data of CK Boo was already applied in the recent paper of Pop et al. (2011) (see also Pop, 1996, 1999 and Pop et al., 2000, 2003). It relies on the classical OC method and takes into account the nonlinear fitting of an ephemeris consisting in the superposi- tion of a K order polynomial term (P[K]) and a multiperiodic term (F[L];[M 1 ],[M 2 ], ... ,[M L ]) t n ¼ t 0 þ X K k¼1 s k n k þ X L l¼1 X M l m¼1 s lm sinðX lm n þ U lm Þ; ð3:1Þ where f 0l = P p /P sl , P p s 1 being the orbital (primary) period, while P sl (l = 1, 2, ... , L) being the secondary periodicities previously detected in the OC curve. The periodicity detection was performed using the Monte Carlo method used by Pop et al. (2010) (see also Pop, 2005, 2007 and Pop and Vamos ß, 2007). 4. Data analysis and results The following linear ephemeris has been established for the analyzed data set of CK Boo t n ¼ 2449135:2816 þ 0:35515426n 0:0013 0:00000013 ð4:1Þ Within the limits of the available data, the OC residuals of this con- tact binary may be well described by the superposition of two mod- ulating signals with periodicities of 24.14 and 10.62 yr, on a parabolic trend corresponding to a linear increase of the orbital per- iod (Fig. 1). The values of the parameters of the established ephem- eris obtained through nonlinear least squares fitting, are given in Table 1, while the steps of prewhitening process, in the amplitude spectra of the OC residuals, are displayed in Fig. 2. The statistical significance of the highest peak appearing in the amplitude spectra of the successive obtained residuals has been estimated using the Monte Carlo method mentioned in Section 3. The involved numer- ical experiments relied on null hypotheses built up using Gaussian noise, bootstrap resampling, and random permutations. In all these cases, the two above periodicities have been detected at confidence levels of 100%. The remaining OC residuals proved to consist of observational noise with Gaussian character only. 4.1. The secular trend in the orbital period The presence of a secular parabolic trend in the OC curve of CK Boo was already emphasized by Qian and Liu (2000), Kalci and Der- man (2005), and Yüce et al. (2006). Our present study confirmed the linear increase of the orbital period (Table 1). Like the above mentioned authors, we interpreted this period variation as an ef- fect of the conservative mass transfer occurring in the system (e.g. Huang, 1963; Pringle, 1975). Accordingly, we found that the primary component receives matter from the secondary one at a rate of _ M 1 ¼ 1:62 10 8 M Sun yr 1 , i.e., a slower rate than estimated by Kalci and Derman (2005) ð _ M 1 ¼ 5:73 10 8 M Sun yr 1 ), and Yüce et al. (2006) ð _ M 1 ¼ 3:92 10 8 M Sun yr 1 ). Our estimation relies on the mass values of the two components given by these last authors. 4.2. The magnetic activity cycles hypothesis As we already mentioned, according to its spectral type (F7/F8), the contact system CK Boo should display an orbital period modu- lation phenomenon induced by the cyclic magnetic activity of the component stars (e.g. Hall, 1989, 1990; Applegate, 1992). Interest- ing, complementary information, which is useful for interpretation of the detected periodicities from the viewpoint of magnetic activ- ity cycles may be obtained from the correlation between the length of the activity cycle (P cyc ) and the angular velocity (X) for active stars found by Lanza and Rodonò (1999), which relies on the dyna- mo theory. They assumed: (i) that the length of the modulation cy- 668 A. Pop, C. Vamos ß / New Astronomy 17 (2012) 667–672
cle in the OC curve (P mod P S in our notations) coincides with P cyc , (ii) the synchronism between the orbital motion and the rotation of the component stars (X =2p/P orb , with P orb P p in our nota- tions). Pop et al. (2011) reconsidered this correlation through lin- ear least squares regression and through a robust linear regression method (the mean median method) proposed by Sârbu (1997). Here, in order to find out the implications of the above correla- tion in case of the orbital period modulation in CK Boo, we further reconsider it by using recently estimated P p and P S values for one Algol system, two RS CVn systems and three W UMa systems taken into account by Lanza and Rodonò (1999). We also included two additional systems: one of Algol type (ST Per), and one of W UMa type (OO Aql) (see Table 2). Both considered Algol systems (see Ta- ble 2), display longer periods previously attributed to the light- time effect related to the presence of an unseen companion (see Panchatsaram, 1981; Qian, 2001 for RT Per, and Demircan and Se- lam, 1993; Borkovits and Hegedüs, 1996; Pop et al., 2000 for ST Per). In case of RT Per, Pop et al., 1996 emphasized the existence of a secondary, shorter periodicity, while in case of ST Per an addi- tional, shorter periodicity has been considered by Demircan and Selam (1993) (see also Borkovits and Hegedüs, 1996) and re-eval- uated by Pop et al. (2000). In case of the W UMa systems investi- gated by Borkovits et al. (2005), these authors were able to prove the coexistence of a light-time effect and an additional short period modulation in case of AB And, OO Aql, and U Peg. Using the actu- alized and extended sample (48 binary systems) of Lanza and Ro- donò (1999), we found the following regression relationships log P mod ½yr¼0:051 0:376 log X½s 1 0:246 0:061; ð4:1Þ in case of ordinary least squares method (OLS), and log P mod ½yr¼ 0:27 0:29 log X½s 1 0:47 0:12; ð4:2Þ in case of robust regression method of Sârbu (1997). Unlike OLS, the robust method is less sensitive to the outlying points existing in the considered sample, and therefore the slopes estimated through these methods may differ from one another. Our changes in the data sample led to a slight increase of the correlation coefficient (in mod- ulus) from r = 0.616 ± 0.093, to r = 0.672 ± 0.080. According to the Eqs. (4.1) and (4.2) we got the following pre- dictions for the case of CK Boo: (i) P mod = 21.68 yr in case of ordin- ary least squares method, and (ii) P mod = 21.86 yr in case of robust regression method. The first periodicity detected in the OC curve of CK Boo is P S1 = 24.14 yr, i.e., higher, but close to the above two predictions. This result could support the hypothesis that magnetic activity cycles are occurring in this system with a length of about 24.14 yr. This hypothesis may also be considered as a cause for the second detected periodicity (P S2 = 10.62 yr). Further information concerning the features related to the mechanism of the cyclic magnetic activity would be obtained through Applegate’s (1992) model. However, we have to take into account that this model (see also Lanza et al., 1998) approaches the case on which just one component star is magnetically active (in Fig. 1. The parabolic+biperiodic model of the OC curve of CK Boo. Table 1 Parameters of the polynomial+multiperiodic (P2+ F2; 1, 1) ephemeris for the OC curve of CK Boo. t 0 ¼ HJD2449135:27712 0:00035 s 1 Pp ¼ 0:355154216dc 1 0:000000042 s 2 ¼ 4:87 10 11 dc 2 0:31 10 11 ð1=PpÞðdPp=dtÞ¼ 2:82 10 7 yr 1 0:27 10 7 f 01 ¼ 0:00004029 0:00000080 P S1 ¼ 24:14 yr 0:48 s 11 ¼ 0:01303 d 0:00035 U 11 ¼ 5:953 rad 0:033 f 02 ¼ 0:0000916 0:0000012 P S2 ¼ 10:62 yr 0:14 s 21 ¼ 0:00449 d 0:00034 U 21 ¼ 0:561 rad 0:069 r res = 0.00206 d A. Pop, C. Vamos ß / New Astronomy 17 (2012) 667–672 669
the case of CK Boo both stars are expected to be magnetically ac- tive). Other complications could arise from the complex geometry of the contact binary systems. Consequently, we restricted our- selves to the estimation of some overall parameters. We shall use the following notations according to the papers of Applegate (1992), Lanza et al. (1998), Lanza and Rodonò (1999): A OC – the semi-amplitude of the OC curve, DP p /P p – the relative semi-ampli- tude of the orbital period modulation, DQ – the variation of the quadrupole moment of the active star which would be related to the observed DP p /P p value. Estimating the semi-major axis of the system’s relative orbit the value a = 2.468 R Sun (see Gazeas et al., 2006), we got: P S1 ¼ 24:14 yr : A OC ¼ 0:02606 d; DP p =P p ¼ 9:29 10 6 ; DQ 1 ¼8:72 10 42 kg m 2 ; DQ 2 ¼9:37 10 41 kg m 2 ; P S2 ¼ 10:62 yr : A OC ¼ 0:00898 d; DP p =P p ¼ 7:27 10 6 ; DQ 1 ¼6:83 10 42 kg m 2 ; DQ 2 ¼7:34 10 41 kg m 2 ; where the subscripts 1 and 2 refers to the primary and secondary components, respectively. The following remarks can be formu- lated: (i) the values of the modulation periods are of the same order of magnitude as other W UMa systems (e.g. Lanza and Rodonò, 1999); (ii) the values of the quadrupole moment estimated for the two component stars and for both modulation periods are of the or- der of 10 41 – 10 42 kg m 2 (see also Borkovits et al., 2005), i.e., of the same order of magnitude as cataclysmic variables (10 42 kg m 2 , see Lanza and Rodonò, 1999). 4.3. The third body hypothesis Other possible mechanism which may be involved in the orbital period modulation is that of the light-time effect due to the pres- ence of one or two unseen companions revolving on circular or- bit(s) around the barycentre of the system. We did not take into account the case of an arbitrary periodic modulation, which corre- sponds to an eccentric orbit (e.g. Kopal, 1978), for the following two reasons: (i) the nonlinear fitting of a P2+ F1;2 model was not convergent, and (ii) in case of the P2+ F2;1,1 model, we found f 02 /f 01 = 2.274 ± 0.054, which means that, within the limits of observational errors, f 02 2f 01 , i.e. f 02 is not the first harmonic of f 01 . For the two detected periodicities we estimated the values of the corresponding parameters: a 12 sini – the projection of the semi- major axis of the close binary absolute orbit on the line of sight, i inclination of the normal to the orbit plane on the line of sight, f(M) – the mass function, M 3 – the mass of the hypothetical companion, and K – the semi-amplitude of the radial velocity curve. Thus, tak- ing into account the values of the masses of the binary component stars estimated by Gazeas et al. (2006, see Table 6 therein), the bin- ary system’s mass is M 12 = 1.597 M Sun , and for zero eccentricity of the system, we obtained: P S1 ¼ 24:14 yr : a 12 sin i ¼ 2:256AU; f ðMÞ¼ 0:0197M Sun ; K ¼ 2:78 km s 1 ; P S2 ¼ 10:62 yr : a 12 sin i ¼ 0:777AU; f ðMÞ¼ 0:0042M Sun ; K ¼ 2:18 km s 1 : We present in Fig. 3 the values for the mass (M 3 ) and the orbit ra- dius (a 3 ) of the companion as functions of the inclination value, computed for the two periodicities. Adopting from the same authors the inclination value estimated for photometric solution with a cool spot on the primary star (see Table 4 of Gazeas et al., 2006) i = 63.7°, we got: P S1 ¼ 24:14 yr : M 3 ¼ 0:492M Sun ; a 3 ¼ 8:16AU; P S2 ¼ 10:62 yr : M 3 ¼ 0:272M Sun ; a 3 ¼ 5:08AU; assuming the coplanarity of orbital planes. In both cases, the respective systems may be considered as hier- archical triple systems: an inner close binary, and a distant com- panion, the semimajor axis of their relative orbit being 2296.645 R Sun (P S1 = 24.14 yr), and 1280.090 R Sun (P S2 = 10.62 yr), respec- tively. Having in view a two-body approximation, we investigated the stability of these hypothetical triplets using the simple stability criterion of Szebehely and Zare (1977). It resulted that both of Table 2 Eclipsing binary systems used to bring up-to-date and to extend the sample of Lanza and Rodonò (1999). Name Class Spectrum X (10 4 s 1 ) P mod (yr) Ref. RT Per Algol F5 + G0 0.86 26.5 Pop et al. (1996) ST Per A3 V + K1–2 IV 0.27 23.8 Pop et al. (2000) RT And RS CVn F8 V + G0 – K1–3 V 1.16 29.4 Pop and Roman (2005) WW Dra G2 + K0 0.16 112.2 Liao and Qian (2010) AB And W UMa G5 V 2.19 18.3 Borkovits et al. (2005) OO Aql G5 V 1.43 20.4 V566 Oph F5 V 1.78 19.8 U Peg G2 V 1.94 17.9 Fig. 2. Amplitude spectra corresponding to different stages of prewhitening of the OC curve. 670 A. Pop, C. Vamos ß / New Astronomy 17 (2012) 667–672
them would be stable either for direct or retrograde motion, for all inclination angles, but ‘‘more stable‘‘, in case of direct motion. It also resulted that the configuration corresponding to the longer periodicity would be ‘‘more stable’’. Our attempt to investigate the orbital period variability through the light-time effect also has to take into account the possible involvement of the recently discovered visual companion. The pro- jected separation of the two visual components of the system (ex- pressed in AU) is d obs = q/p, where q = 0.12 00 is the angular separation (see Pribulla and Rucinski, 2006; Rucinski et al., 2007), and p Hp = 0.00696 00 is the system parallax as determined by Hipparcos. According to the orientation of the relative orbit of the wide system, there are two extreme situations for the value of d obs : (i) maximum apparent separation, when the semimajor axis is perpendicular on the line of sight, i.e., d obs = d max = a w (1 + e), and herefrom a w = d obs /(1 + e); and (ii) minimum apparent separation, when the semimajor axis coincides with the line of sight, i.e., d obs = d min = a w (1 e) cosi, and thus a w = d obs /((1 e) cosi). We de- noted: a w = a 12 + a 3 – the semimajor axis of the relative orbit of the wide system, e – the eccentricity of the relative orbit of the wide system, and i – the wide orbit inclination with respect to the line of sight. Using Kepler’s third law one can estimate the value of the orbital period of the visual binary system: P w ¼ 2pða 3 w =ðGðM 12 þ M 3 ÞÞÞ 1=2 , where M 12 = 1.597M Sun (Gazeas et al., 2006), and M 3 = 0.45 M Sun (Rucinski et al., 2007). Within the hypothesis d obs = d max , for e = 0, we got a w = 17.24 AU, P w = 50.04 yr, and for e = 0.99, a w = 8.66 AU, P w = 17.82 yr. Taking into account the other hypothesis, i.e., d obs = d min , then for i ? 90 , we obviously have a w , P w ? 1. As we previously argued, presently, the light- time effect can be taken into account only for circular orbits. In this case (e = 0), for i =0° we have: a w = 17.24 AU, P w = 50.04 yr, while for i = 63.7° (the orbits of the close pair and that of the wide system are coplanar): a w = 38.91 AU, P w = 169.66 yr. Comparing these esti- mates with the values of P S1 , P S2 , as well as with those of a 3 dis- played in Fig. 3, it becomes obvious that the observed orbital period variability of CK Boo is not an effect of the gravitational cou- pling with its visual companion. 5. Concluding remarks The main result of our study on the orbital period variability of the contact binary CK Boo is that it is featured by a biperiodic mod- ulation superposed on an already known linear increase. The in- volved periodicities are P S1 = 24.14 yr and P S2 = 10.62 yr. We have to remark that the actual time base covers only 1.5 cycles of P S1 . In case of P S2 the coverage is better, about 3.4 cycles. Consequently, P S1 has to be confirmed by future observations. According to the prediction supplied by the correlation found by Lanza and Rodonò (1999) and updated here by us, the longest periodicity detected in the OC curve of CK Boo could be that of magnetic activity cycles occurring in this system. However, as Lan- za and Rodonò (1999) mentioned, this correlation did not take ac- count the situation of multiperiodicity existing in the OC curves of some eclipsing binary systems. They also could not exclude a pos- sible contamination of their sample with some systems displaying an orbital period modulation caused by the light-time effect. Therefore, neither their regression relationship nor those estab- lished by us (Eqs. (4.1) and (4.2)) can be considered more than pre- liminary and approximate ones. Consequently, for the moment, through this approach, it is not possible to decide whether P S1 , or P S2 , or even both of them are involved in the magnetic cycles occur- ring in the system CK Boo. The interpretation of the two periodic modulations of the orbital period as being caused by the light-time effect allowed us to esti- mate the masses and the semimajor axes of the orbits of the hypo- thetical companions. According to Szebehely and Zare (1977) stability criterion, stable triple systems involving low-mass com- panions corresponding either to P S1 or to P S2 , would be allowed (see Section 4.3). However, for i = 63.7°, in the hypothesis of orbits coplanarity, the masses of the hypothetical companions (M 3 = 0.492 M Sun , and M 3 = 0.272 M Sun , respectively) would be larger than the mass of the secondary component of CK Boo (M 2 = 0.155 M Sun ). Moreover, according to Fig. 3, even for i = 90°, the lowest mass of the companion would be about 0.433 M Sun , or 0.242 M Sun , respec- tively. Especially in case of the hypothetical companion corre- sponding to P S1 , the corresponding mass would be very close to that estimated for the visual companion, and consequently they would have quite similar evolutionary stages, and relatively close luminosities. The hypothesis of a white dwarf companion, which would be in agreement with the small ultraviolet excess discovered by Krzesin´ ski et al. (1991), may also be taken into account. Concern- ing the visual companion of CK Boo, it appears not to be involved in the orbital period modulation of this binary system. Taking into account the above considerations, it seems that at least one of the periodicities detected in the OC curve of CK Boo (we recall that P S1 needs future observational confirmation) could be related to the development of magnetic activity cycles in the component stars of this system. However, the possibility that CK Boo is a quadruple system cannot be ruled out. Acknowledgments This research has made use of NASA’s Astrophysics Data System Bibliographic Services, the Simbad database, operated at CDS, Fig. 3. The mass (lower panel) and the semimajor axis of the absolute orbit (upper panel) values of the hypothetical ‘‘third body’’ corresponding to the two detected periodicities, for different inclination values. A. Pop, C. Vamos ß / New Astronomy 17 (2012) 667–672 671
Strasbourg, France, the Lichtenknecker-Database of the BAV oper- ated by the Bundesdeutsche Arbeitsgemeinschaft für Veränderli- che Sterne e.V. (BAV), and ‘‘OC Gateway Database of times of minima (E) and maxima (RR)’’ of Variable Star Section of Czech Astronomical Society. References Abhyankar, K.D., Panchatsaram, T., 1984. Light time effect in TW Draconis. MNRAS 211, 75. Applegate, J.H., 1992. A mechanism for orbital period modulation in close binaries. ApJ 385, 621. Aslan, Z., Derman, E., 1986. Photoelectric observations of CK bootis. A&AS 66, 281. Barning, F.J.M., 1963. The numerical analysis of the light-curve of 12 Lacertae. BAN 17, 22. Borkovits, T., Hegedüs, T., 1996. On the invisible components of some eclipsing binaries. A&AS 12, 63. Borkovits, T., Csizmadia, Sz., Hegedüs, T., Bíró, I.B., Sándor, Zs., Opitz, A., 2002. Complex period variations in the binary system IM Aurigae. A&A 392, 895. Borkovits, T., Elkhateeb, M.M., Csizmadia, Sz., Nuspl, J., Bı ´ ró, I.B., Hegedüs, T., Csorvási, R., 2005. Indirect evidence for short period magnetic cycles in W UMa stars. Period analysis of five overcontact systems. A&A 441, 1087. Breger, M., Huang, L., Jiang, S.-Y., Guo, Z.-H., Antonello, E., Mantegazza, L., 1987. Multiple close frequencies of the Delta Scuti star h 2 Tau. A&A, 175, 117. D’Angelo, C., van Kerkwijk, M.H., Rucinski, S.M., 2006. Contact binaries with additional components. II. A spectroscopic search for faint tertiaries. AJ 132, 650. Demircan, O., 1987. BV observations of W UMa-type binaries: CK Bootis, BI CVn, and AH Vir. Ap&SS 135, 169. Demircan, O., Selam, S.O., 1993. Long-term behaviour of the orbital period of Algol- type binary ST Persei. A&AS 98, 513. Derman, E., Demircan, O., 1992. A period study of GK Cephei. AJ 103, 599. Diethelm, R., 2011.Timings of minima of eclipsing binaries. IBVS, 5992. Gazeas, K.D., Niarchos, P.G., Zola, S., Kreiner, J.M., Rucinski, S.M., 2006. Physical parameters of components in close binary systems: VI. AcA 56, 127. Hall, D.S., 1989. The relation between RS CVn and Algol. SSR 50, 219. Hall, D.S., 1990. Period changes and magnetic cycles. In: _ Ibanog˘lu, C. (Ed.), Active Close Binaries. Kluwer Academic Publishers, Dordrecht, p. 95. Huang, S.-S., 1963. Modes of mass ejection by binary stars and effect on their orbital periods. ApJ 138, 471. Jia, G., Liu, X., Huang, H., 1992. 1991 Photometry of the W UMa type binary CK Boo. IBVS, 3727. Kalci, R., Derman, E., 2005. CK Bootis: a W UMa system with a small mass ratio. AN 326, 342. Kalimeris, A., Rovithis-Livaniou, H., Rovithis, P., 1994a. On the orbital changes in contact binaries. A&A 282, 775. Kalimeris, A., Rovithis-Livaniou, H., Rovithis, P., Oprescu, G., Dumitrescu, A., Suran, M.D., 1994b. An orbital period study of the contact system AB Andromedae. A&A 291, 765. Kopal, Z., 1978. Dynamics of Close Binary Systems. D. Reidel Publishing Company, Dordrecht. Kreiner, J.M., 1971. Investigation of changes in periods of eclipsing variables. AcA 21, 365. Krzesin´ ski, J., Mikolajewski, M., Pajdosz, G., Zola, S., 1991. An observational study of a spotted W UMa-type binary star CK Bootis. Ap&SS 184, 37. Lanza, A.F., Rodonò, M., 1999. Orbital period modulation and quadrupole moment changes in magnetically active close binaries. A&A 349, 887. Lanza, A.F., Rodonò, M., Rosner, R., 1998. Orbital period modulation and magnetic cycles in close binaries. MNRAS 296, 893. Li, Y., Qian, S.B., 2005. GMG 2.4-m telescope and light-time effect research for W UMa-type binaries. In: Sterken, C. (Ed.), ASP Conference Series, vol. 335, p. 245. Liao, W.-P., Qian, S.-B., 2010. The most plausible explanation of the orbital period changes in close binaries: the case of the RS CVn-type binary WW Dra. MNRAS 405, 1930. López de Coca, P., Garrido, R., Rolland, A., 1984. Mode identification for the Delta Scuti variable HR 7222. A&AS 58, 441. Panchatsaram, T., 1981. Light-time effect in RT Persei. Ap&SS 77, 179. Panchatsaram, T., Abhyankar, K.D., 1981. SW Lacertae – a quadruple system. Bull. Astron. Soc. India 9, 31. Pop, A., 1996. The use of characteristic time to the study of period variability phenomena. Romanian Astron. J. 6, 147. Pop, A., 1999. A multiperiodic ephemeris for RZ Cephei. In: IBVS, vol. 4801. Pop, A., 2005. On the orbital period modulation of the eclipsing binary system Y Leonis. In: Sterken, C. (Ed.), ASP Conference Series, vol. 335, p. 263. Pop, A., 2007. On the detection and diagnosis of low-level stellar variability. Romanian Astron. J. 17, 35. Pop, A., Roman, R., 2005. RT and: multiperiodic orbital period variability, resonance, and the hypothesis of the light-time effect. Romanian Astron. J. 15, 163, Supplement. Pop, A., Vamosß, C., 2007. The detection and diagnosis of low-level stellar variability through resampling methods. Romanian Astron. J. 17, 101. Pop, A., Todoran, I., Agerer, F., 1996. Multiperiodic behaviour of RT Persei. Romanian Astron. J. 6, 141. Pop, A., Liteanu, V., Pavel, C., 2000. On the multiperiodicity of the orbital period variation of ST Persei. Romanian Astron. J. 10, 25. Pop, A., Liteanu, V., Moldovan, D., 2003. On the period variability of GP Andromedae. Ap&SS 284, 1207. Pop, A., Vamosß, C., Turcu, V., 2010. Deterministic components in the light curve amplitude of Y Oph. AJ 139, 425. Pop, A., Turcu, V., Marcu, A., 2011. The intriguing orbital period variability of Y Leonis. Ap&SS 333, 17. Pribulla, T., Rucinski, S.M., 2006. Contact binaries with additional components. I. The extant data. AJ 131, 2986. Pringle, J.E., 1975. Period changes in eruptive binaries. MNRAS 170, 633. Qian, S., 2001. Possible mass and angular momentum loss in Algol-type binaries. V RT Persei and TX Ursae Majoris. AJ 122, 2686. Qian, S., Liu, Q., 2000. A possible connection between the variability of light curve and the change of the orbital period in the contact binary CK Bootis. Ap&SS 271, 331. Rafert, J.B., 1982. Periodic ephemerides for 49 eclipsing binary-star systems. PASP 94, 485. Rucinski, S., Lu, W., 1999. Radial velocity studies of close binary stars. II. AJ 118, 2451. Rucinski, S.M., Pribulla, T., van Kerkwijk, M.H., 2007. Contact binaries with additional components. III. A search using adaptive optics. AJ 134, 2353. Samolyk, G., 2011. Recent minima of 146 eclipsing binary stars. J. Amer. Assoc. Var. Star Obs. 39, 177. Sârbu, C., 1997. A new median – based robust regression method. Croat. Chem. Acta 70, 795. Szebehely, V., Zare, K., 1977. Stability of classical triplets and their hierarchy. A&A 58, 145. Wehlau, W., Leung, K.-C., 1964. The multiple periodicity of Delta Delphini. ApJ 139, 843. Wood, D.B., Forbes, J.E., 1963. Ephemerides of eclipsing stars. AJ 68, 257. Yüce, K., Selam, S.O., Albayrak, B., Ak, T., 2006. Monitoring secular orbital period variations of some eclipsing binaries at the Ankara University Observatory. Ap&SS 304, 67. 672 A. Pop, C. Vamos ß / New Astronomy 17 (2012) 667–672
2012

Related Posts