Analysis of the Activity of the Eclipsing Binary WZ Cephei
Total Page:16
File Type:pdf, Size:1020Kb
ASTRONOMY & ASTROPHYSICS JULY 1998,PAGE17 SUPPLEMENT SERIES Astron. Astrophys. Suppl. Ser. 131, 17–23 (1998) Analysis of the activity of the eclipsing binary WZ Cephei G. Djuraˇsevi´c1,M.Zakirov2, A. Hojaev2, and G. Arzumanyants2 1 Astronomical Observatory, Volgina 7, 11050 Belgrade, Yugoslavia e-mail: [email protected] 2 Astronomical Institute of Uzbek Academy of Sciences, Astronomicheskaya 33, 700052, Tashkent, Uzbekistan Received September 18, 1997; accepted January 28, 1998 Abstract. The present study deals with the problem of 1. Introduction the orbital and physical parameters’ estimation of the ac- tive eclipsing binary WZ Cep based on the interpretation The eclipsing binary WZ Cep is a short period of its photometric observations. The B and V light curves (P =0d.417447)system showing an EB-type light curve. (Hoffmann 1984), as well as the new ones obtained at the AccordingtotheGeneral Catalogue of Variable Stars - Maidanak Observatory in the B, V and R passbands dur- GCVS (Kholopov et al. 1985) the spectral type of the ing 1995, are analysed and discussed. These two groups system is F5. Photoelectric observations (Hoffmann 1984) of light curves show large differences in their shape and show a high degree of asymmetry of the light curve. in the depth of the primary minimum. Hoffmann’s light Namely, the primary maximum (phase 0.25) is signifi- curves show big asymmetry, especially conspicuous in the cantly higher than the secondary one. In Balazs’ (1937) different height of the successive maxima. The new obser- and Detre’s (1940) observations we can find the changes in vations give light curves that are almost symmetric having the light curve shape and the depth of successive minima. significantly shallower primary minimum. New photoelectric observations, obtained at the Maidanak They are analysed in the framework of the Roche Observatory during 1995, show that the light curve of computer model (Djuraˇsevi´c 1992a) by applying the WZ Cep dramatically changed its shape. Compared with inverse problem method (Djuraˇsevi´c 1992b) based on Hoffmann’s light curve, the depth of the primary minimum Marquardt’s (1963) algorithm. The model is generalised is significantly lower and the light curve is quite symmet- for the case of an overcontact configuration also. To ex- ric. This variability in the light curve form may be at- plain these light curves asymmetries and variations, the tributed to a variable star spot activity and to a variable model involved spotted regions on the components. The level of the energy transfer between the components. analysis shows that WZ Cep system is in a slight overcon- By analysing Hoffman’s observations with the Wilson- tact configuration and that best fitting is achieved with Devinney method Kalu˙zny (1986) obtained the solutions spotted regions on the primary star (SpF5), with the tem- that WZ Cep is in the overcontact configuration with the perature contrast between the spotted area and the sur- components having significantly different effective temper- rounding photosphere being As = Ts/T1 ∼ 0.70. The mass atures. He concluded that observed light curves (with large ratio was estimated at the value q = m2/m1 ∼ 0.33, based degree of asymmetry) can be reproduced theoretically if on the light-curve analysis. one assumes the existence of a relatively small hot spot re- The basic parameters of the system and of the spot- gion placed on the less massive component, near the neck ted areas are estimated for all foregoing mentioned light between stars. curves. The obtained solutions show that there are no sig- The reanalysis of this observational material in the nificant variations of the system’s basic parameters es- frame of Roche computer model (Djuraˇsevi´c 1992a) gener- timated by analysing these two groups of very different alised to the case of the overcontact configuration showed light curves. Consequently, the main variations in the light that the hot spot hypothesis cannot offer a satisfactory curves can be explained by the change of the position and solution. Moreover, the analysis of Hoffman’s and new ob- the size of the spotted areas on the primary. servations presented in this paper proved that the changes Key words: stars — activity — W UMa type — eclipsing in the light curves of WZ Cep are caused by the presence binaries-stars — individual — WZ Cep and the development of the dark spotted areas on the system’s components. This type of activity is relatively well established within great number of stars with convec- Send offprint requests to:G.Djuraˇsevi´c tive envelopes. Based on the spectral type (F5) and the 18 G. Djuraˇsevi´c et al.: Analysis of the activity of the eclipsing binary WZ Cephei Table 1. Basic points on the averaged light curve of WZ Cep the inverse problem method (Djuraˇsevi´c 1992b) based on obtained during 1995 Marquardt’s (1963) algorithm. Phase V U − B B − V V − R According to this method, the stellar size in the model Max 12.45 0.00 0.51 0.58 is described by the filling coefficients for the critical Roche Min I 13.01 0.27 0.62 0.57 lobes F1,2 of the primary and secondary component, re- Min II 12.91 0.12 0.61 0.55 spectively, which indicate to what degree the stars of the system fill their corresponding critical lobes. For syn- chronous rotation of the components these coefficients are solutions obtained by analysing the light curves we think expressed via the ratio of the stellar polar radii, R1,2,and that the WZ Cep system exhibits this type of activity also. the corresponding polar radii of the critical Roche lobes, i.e., F1,2 = R1,2/RRoche1,2 . The model is generalised to the case of an overcontact configuration, which is characteris- 2. New observations tic for the W UMa-type systems. In this case the poten- tial Ω1,2 characterising the common photosphere is derived New photoelectric observations of WZ Cep are obtained with the filling coefficient of the critical Roche lobe F1 > 1 by using the telescope with the aperture of 0.6 m at of the primary, while the coefficient F2 may be excluded Maydanak Observatory (Uzbekistan) during 1995. SAO from further consideration. The degree of overcontact, is 10752 (V =9m.045, U − B = −0m. 042, B − V =0m.458 defined in the classical way (Lucy & Wilson 1979) as: and V − R =0m.422) is chosen as a comparison star. 428 (U), 438 (B), 452 (V ) and 436 (R) observations are ob- fover[%] = 100 · (Ω1,2 − Ωi)/(Ω0 − Ωi), tained during 29 nights. The data are available on request where Ω1,2,Ωi,andΩ0 are the potentials of the com- from the authors. Some basic features of these light curves mon photosphere and of the inner and outer contact sur- are presented in Table 1. Due to the large error in the faces, respectively. Since the potentials Ωi and Ω0 include photon registration with the small telescope aperture in Lagrangean equilibrium points L1 and L2 respectively, the U passband the corresponding light curve is poorly they can be easily derived. Namely, in the mentioned defined. Therefore, we excluded the U passband measure- points the net force acting in the system is equal to zero. ments from further analysis. In a spherical coordinate system the surfaces of the The use of the ephemeris given in the GCVS shows components are divided into a large number of elemen- that the secondary minimum of the light curve is deeper tary cells, whose intensity and angular distribution of ra- than the primary one. Doubting about this result, we re- diation are determined by the stellar effective tempera- analysed all published observations of the light curve’s ture, limb-darkening, gravity-darkening and by the effect minima for WZ Cep. We established that Balazs (1937) of reflection in the system. The presence of spotted ar- had correctly found the epoch of the primary minimum. eas (dark or bright) enables one to explain the asymme- Later, in the ephemeris given by Detre (1940), the epochs tries, the light curve anomalies and the O’Connell effect of the primary and the secondary minima were switched. on the light curves of active CBs. In our programme these Hoffmann (1984) also noticed this error in the ephemeris active regions are approximated by circular spots, char- by Detre. Since Hoffman’s paper has not drawn much at- acterised by the temperature contrast of the spot with tention, the incorrect ephemeris by Detre were adopted respect to the surrounding photosphere (As = Ts/T∗), by in GCVS. Based on all available observations of the pri- the angular dimension (radius) of the spot (θ)andbythe mary minimum our analysis give the corrected value for longitude (λ) and latitude (ϕ) of the spot centre. The the ephemeris of WZ Cep: longitude (λ) is measured clockwise (as viewed from the direction of the +Z-axis) from the +X-axis (line connect- Min I = Hel. JD 2449890.3552 + 0d. 41744677 × E ing star centres) in the range 0o − 360o. The latitude (ϕ) ±0.0014 ± 0d. 00000009. is measured from 0o at the stellar equator (orbital plane) to +90o towards the “north” (+Z)and−90o towards the “south” (−Z) pole. The generalisation of the computer 3. Roche model and procedure of light curve analysis model (Djuraˇsevi´c 1992a) to the case of an overcontact configuration reveals the high importance of the “visibil- In the present paper, the WZ Cep photoelectric obser- ity” of elementary cells around the neck region. In order vations were analysed to study the activity of the sys- that, at a chosen orbital phase, the radiation from elemen- tem. For analysing these asymmetric light curves, proba- tary cells, in which the stellar surface is divided, reaches bly deformed by the presence of spotted areas on the com- the observer, a classical procedure is to check the fulfilling ponents, we used the programme of Djuraˇsevi´c (1992a), of the condition: cos γ>0; no eclipsing by other sys- which is based on the Roche model and the princi- tem’s component.