Research Article  Open Access
M. M. Elkhateeb, M. I. Nouh, E. Elkholy, B. Korany, "An Extensive Photometric Investigation of the W UMa System DK Cyg", Journal of Astrophysics, vol. 2015, Article ID 590673, 8 pages, 2015. https://doi.org/10.1155/2015/590673
An Extensive Photometric Investigation of the W UMa System DK Cyg
Abstract
DK Cyg is a contact binary system that undergoes complete eclipses. All the published photoelectric data have been collected and utilized to reexamine and update the period behavior of the system. A significant period increase with rate of 12.590 × 10^{−11} days/cycle was calculated. New period and ephemeris have been calculated for the system. A long term photometric solution study was performed and a light curve elements were calculated. We investigated the evolutionary status of the system using theoretical evolutionary models.
1. Introduction
The eclipsing binary DK Cyg (BD +330 4304, 10.37–10.93 m_{v}) is a well known contact binary system with a period of about 0.4707 days. It was discovered as variable star earlier by Guthnick and Prager [1], so their epoch of intensive observations is very long. The earliest photographic light curve classified the system as a W Ursae Majoris type. Rucinski and Lu [2] carried out the first spectroscopic observations and estimated the mass ratio as and classified the system as an Asubtype contact system with spectral type of A8V. Visual light curves were published by Piotrowski [3] and Tsesevitch [4] from Klepikova’s observations.
First photoelectric observations for the system were carried out by Hinderer [5], while Binnendijk [6] observed the system photoelectrically in B and Vbands and derived least squares orbital solution. The system DK Cyg was classified in the General Catalogue of Variable Stars as A7V [7], while Binnendijk [6] adopted it as A2. Mochnacki and Doughty [8] showed that the color index of the system judged its spectral type and found that the spectral type of the system is more likely to be about F0 to F2. Because the system DK Cyg is a summer object in the Northern hemisphere with 11.5hour period and short durations of night, it is bound to remain illobserved [9]. Only four complete light curves by Binnendijk [6], Paparo et al. [10], Awadalla [9], and Baran et al. [11] are published. Photoelectric observations and new times of minima have been carried out by many authors: Borkovits et al. [12], Sarounová and Wolf [13], Drozdz and Ogloza [14], Hübscher et al. [15], Dogru et al. [16], Hubscher et al. [17], Hubscher et al. [18], Erkan et al. [19], Diethelm [20], Dogru et al. [21], Simmons [22], Diethelm [23], and Diethelm [24].
In the present paper we are going to perform comprehensive photometric study for the system DK Cyg. The structures of the paper are as follows: Section 2 deals with the period change, Section 3 is devoted to the light curve modeling, and Section 4 presents the discussion and conclusion reached.
2. Period Change
Although the period variation of contact binary systems of the W UMatype is a controversial issue of binary star astrophysics, the cause of the variations (long as well as short term) is still a mystery for a discussion of possible physical mechanisms [25]. Magnetic activity cycle is one of the main mechanisms that caused a period variation together with the mass exchange between the components of each system. Kaszas et al. [26] stated that the long term period variation may be interpreted by a perturbation of the third companion or surface activity of the system components.
Observations by Binnendijk [6] showed a change of the secondary minimum depth and a new linear light element was derived. Period study by Paparo et al. [10] showed that the orbital period of the system DK Cyg increases and the first parabolic light elements were calculated, which confirmed the light curve variability. Kiss et al. [25] updated the linear ephemeris of DK Cyg, while Awadalla [9] recalculated a new quadratic element for the system and confirmed the light curve variability suggested by Paparo et al. [10]. Wolf et al. [27] used a set of 101 published times of minimum covering the interval between 1926 and 2000 in order to update the quadratic element calculated by Awadalla [9]. They showed that the period increases by the rate days/cycle. Borkovits et al. [28] follow the period behavior of the system using set of published minima from HJD 2424760 to HJD 2453302.
In this paper we studied the orbital period behavior of the system DK Cyg using the diagram based on more complete data set collected from the literatures and databases of BAV, AAVSO, and BBSAG observers. Part of our collected data set was given by Kreiner et al. [29]; unpublished Hipparcos observations and main part were downloaded from website (http://astro.sci.muni.cz/variables/ocgate/); Table 1 listed only those minima not listed on the mentioned websites. A total of 195 minima times were incorporated in our analysis covering about 86 years (66689 orbital revolutions) from 1927 to 2013. It is clear that our set of data added about 94 of new minima and increases the interval limit of the orbital period study about 13 years more than the data of Wolf et al. [27], which may give more accurate insight on the period behavior of the system. The different types of the collected minima (i.e., photographic, visual, photoelectric, and CCD) were weighted according to their type. The residual ’s were computed using Binnendijk [6] ephemeris (1) and represented in Figure 1. No distinction has been made between primary and secondary minima: It can be seen from the figure that the behavior of the orbital period of the system DK Cyg shows a parabolic distribution which is generally interpreted by the transfer of mass from one component to the other of binary. Reasonable linear least squares fit of the data available improved the light elements given in (1) to The linear element yields a new period of days which is longer by 0.13 seconds with respect to the value given by Binnendijk [6]. Quadratic least squares fit gives The rate of period increasing resulting from the quadratic elements (3) is days/cycle or days/year or 0.84 seconds/century. More future systematic and continuous photometric observations are needed to follow a continuous change in the orbital period of the system DK Cyg which may show a periodic behavior. The fourth column of Table 1 represents the quadratic residuals calculated using the new element of (2) and represented in Figure 2. All published linear and quadratic elements together with that resulting from our calculations are listed in Table 2. It is noted from the table that the quadratic term resulting from our calculations has slightly higher value than that calculated by Awadalla [9], Wolf et al. [27], and Borkovits et al. [28]. This can be interpreted by the increasing of the set of minima in our study compared to the one they used (nearly double); also we covered an interval larger than the one they used.
 
(1) Szafraniec [30]; (2) Hubscher et al. [31]; (3) Wolf et al. [27]; (4) Hipparcos observations (unpublished); (5) Baldwin and Samolyk [32]; (6) Kiss et al. [25]; (7) Borkovits et al. [33]; (8) Borkovits et al. [28]; (9) Gerner [34]; (10) Menzies [35]; (11) Samolyk [36]; (12) Simmons [22]. 

3. Light Curve Modeling
Light curve modeling for the system DK Cyg by Mochnacki and Doughty [8] using Binnendijk [6] observations in Vband showed nonmatching between the theoretical curve and the observations. The photometric mass ratio calculated from their accepted model was , while the spectroscopic value estimated using radial velocity study by Rucinski and Lu [2] is . Baran et al. [11] estimated an alternating model of spectroscopic and photometric data based on iterative solutions. Their model shows a better fit by introducing a cool spot on the surface of the more luminous component and adopted the third light as free parameter in the computations. On the other hand Rucinski and Lu [2] stated that they did not find any evidence for the existence of a third component in the system during their spectroscopic study. They refer the probability for the presence of a third star in the system to the diagram [29], which shows sinusoidal variation. This evidence is weak because only one cycle is covered up to date.
In the present work we used complete light curves published by Binnendijk [6], Paparo et al. [10], Awadalla [9], and Baran et al. [11] in Vband through a long term photometric solution study in order to estimate the physical parameters of the system and to follow its evolutionary status. The collected light curves showed a flatbottom minima and O’Connell effect. Observations by Paparo et al. [10] and Awadalla [9] displayed some scattering specially at the two maxima of their light curves. Also the observed light curve by Awadalla [9] shows sudden increase in the light level at secondary minimum with respect to the other collected curves.
Photometric analysis for the studied light curves of the system DK Cyg was carried out using Mode 3 (overcontact) of WDint56a Package [43] based on the 2009 version of Wilson and Devinney (WD) code with Kurucz model atmospheres [44–46]. The observed light curves were analyzed using all individual observations. Appropriate gravity darkening and bolometric albedo exponents were assumed for the convective envelope. We adopted [47] and [48]. Bolometric limb darkening values are adopted using the table of van Hamme [49]. Temperature of the primary star was adopted according to Baran et al. [11] model (°K).
The adjustable parameters are the mean temperature of the secondary component , orbital inclination , and the potential of the two components , while the spectroscopic mass ratio () by Baran et al. [11] was fixed for all calculated models together with the primary star’s temperature ().
We started modeling using as initial values the parameters of Baran et al. [11] solution based on cool spot on the luminous components and a third light as a free parameter. The used parameters show disagreement between the theoretical and observed light curves, except for Baran et al. observations. Regarding the conclusion of Rucinski and Lu [2], which stated a weak evidence of presence of a third component, we tried to construct a spotted model without the third light. We constructed a model including two spots; the first one is a cool spot located on surface of the more massive component, while the other is a hot spot located on the surface of the other component. The accepted model reveals good agreement between theoretical and observed light curves for all collected data. Table 3 lists the calculated parameters for the four light curves, while Figure 3 represented the theoretical light curves according to the accepted solution together with the reflected points in Vband. The values in Table 3 are indicative of comparisons in future studies, since the number of observations and the accuracy are not the same in the four light curves. Absolute physical parameters for each component of the system DK Cyg were calculated based on the results of the radial velocity data of Baran et al. [11] and our new photometric solution for each light curve. The calculated parameters are listed in Table 3. The results show that the primary component is more massive and hotter than the secondary component. A threedimensional geometrical structure for the system DK Cyg is displayed in Figure 4 using the software Package Binary Maker 3.03 [50] based on the calculated parameters resulting from our models.
 
Not adjusted. 
4. Discussion and Conclusion
Studying of the period behavior of the system DK Cyg based on all available published times of minima, covering 86 yr of observations including 195 times of light minima, shows a continuous period increase with the rate days/cycle or days/year or 0.84 seconds/century. New linear and quadratic elements were calculated using all available published data and yield a new period of days. A long term photometric study was performed using published observations by Binnendijk [6], Paparo et al. [10], Awadalla [9], and Baran et al. [11]. More systematic and continuous photometric observations for the system DK Cyg are needed to confirm a continuous change in the period and follow its light curve variation.
One of the difficulties for W UMa binaries is to use stellar models of single stars to investigate the evolutionary status of these systems. However, using these theoretical models may give approximate view about the evolutionary status of the system.
We used the physical parameters listed in Table 4 to investigate the current evolutionary status of DK Cyg. In Figures 5 and 6, we plotted the components of DK Cygon on the massluminosity () and massradius () relations along with the evolutionary tracks computed by Girardi et al. [51] for both zero age main sequence stars (ZAMS) and terminal age main sequence stars (TAMS) with metallicity . As it is clear from the figures, the primary component of the system is located nearly on the ZAMS for both the  and  relations. The secondary component is above the TAMS track for  and the  relations. For the sake of comparison, we plotted sample of Atype contact binaries listed in Table 5. The components of DK Cyg have the same behavior of the selected Atype systems.
 
Note: subscripts 1 and 2 mean primary and secondary component, respectively. 
 
(1) Essam et al. [38], (2) Elkhateeb and Nouh [39], (3) Elkhateeb and Nouh [40], (4) Yang et al. [41], and (5) Degirmenci [42]. 
The masseffective temperature relation () for intermediate and low mass stars [37] is displayed in Figure 7. The location of our mass and radius on the diagram revealed a good fit for the primary and poor fit for the secondary components. This gave the same behavior of the system on the massluminosity and massradius relations.
Conflict of Interests
The authors declare that they have no conflict of interests regarding the publication of this paper.
Acknowledgments
This research has made use of NASA’s ADS and the available online material of the IBVS. The authors sincerely thank Dr. Bob Nelson, who allowed them to use his windows interface code WDwint56a, for his helpful discussions and advice.
References
 P. Guthnick and R. Prager, “Benennung von veränderlichen Sternen,” Astronomische Nachrichten, vol. 1011, pp. 161–206, 1936. View at: Publisher Site  Google Scholar
 S. M. Rucinski and W. Lu, “Radial velocity studies of close binary stars. II,” The Astronomical Journal, vol. 118, no. 5, pp. 2451–2459, 1999. View at: Publisher Site  Google Scholar
 S. Piotrowski, “Observations photometriques de DK cyg et de BB peg,” Bulletin de L'Observatoire Astronomique de Belgrade, vol. 1, no. 2, p. 9, 1936. View at: Google Scholar
 L. Binnendijk, “The light variation and orbital elements of DK Cygni,” The Astronomical Journal, vol. 69, pp. 157–164, 1964. View at: Publisher Site  Google Scholar
 F. Hinderer, “Lichtelektrische untersuchungen an W ursae majorissternen,” Journal des Observateurs, vol. 43, p. 161, 1960. View at: Google Scholar
 L. Binnendijk, “The light variation and orbital elements of DK Cygni,” The Astronomical Journal, vol. 69, pp. 157–164, 1964. View at: Google Scholar
 N. Samus, O. Durlevich, V. Goranskij et al., 2014, http://www.sai.msu.su/gcvs/cgibin/search.htm.
 S. W. Mochnacki and N. A. Doughty, “Models for five W ursae majoris systems,” Monthly Notices of the Royal Astronomical Society, vol. 156, no. 2, pp. 243–252, 1972. View at: Publisher Site  Google Scholar
 N. Awadalla, “Photoelectric light curve study of the DK Cygni system,” Astronomy & Astrophysics, vol. 289, pp. 137–140, 1994. View at: Google Scholar
 M. Paparo, M. Hamdy, and I. Jankovics, “Photoelectric observation of DK Cyg,” Information Bulletin on Variable Stars, no. 2838, 1985. View at: Google Scholar
 A. Baran, S. Zola, S. Rucinski, J. Kreiner, and M. Drozdz, “Physical parameters of components in close binary systems: II,” Acta Astronautica, vol. 54, pp. 195–206, 2004. View at: Google Scholar
 T. Borkovits, I. B. Bíró, Sz. Csizmadia et al., “New times of minima of eclipsing binary systems,” Information Bulletin on Variable Stars, vol. 5579, 2004. View at: Google Scholar
 L. Sarounová and M. Wolf, “Precise CCD times of minima of selected eclipsing binaries,” Information Bulletin on Variable Stars, no. 5594, 2005. View at: Google Scholar
 M. Drozdz and W. Ogloza, “Photoelectric minima of eclipsing binaries,” Information Bulletin on Variable Stars, vol. 5623, p. 1, 2005. View at: Google Scholar
 J. Hübscher, A. Paschke, and F. Walter, “Photoelectric minima of selected eclipsing binaries and maxima of pulsating stars,” Information Bulletin on Variable Stars, vol. 5731, pp. 1–31, 2006. View at: Google Scholar
 S. Dogru, A. Dönmez, M. Tüysüz et al., “New times of minima of some eclipsing binary stars,” Information Bulletin on Variable Stars, vol. 5746, 2007. View at: Google Scholar
 J. Hubscher, H.M. Steinbach, and F. Walter, “Minima and maxima of 292 variables (Hubscher+, 2008),” Information Bulletin on Variable Stars, vol. 5830, 2008. View at: Google Scholar
 J. Hubscher, P. Lehmann, G. Monninger, H. Steinbach, and F. Walter, “BAVresults of observations—photoelectric minima of selected eclipsing binaries and maxima of pulsating stars,” Information Bulletin on Variable Stars, vol. 5941, 2010. View at: Google Scholar
 N. Erkan, A. Erdem, T. Akin, F. Aliçavus, and F. Soydugan, “New times of minima of some eclipsing binary stars,” Infornation Bulletin on Variable Stars, vol. 5924, 2010. View at: Google Scholar
 R. Diethelm, “Timings of minima of eclipsing binaries,” Information Bulletin on Variable Stars, vol. 5920, 2010. View at: Google Scholar
 S. Dogru, A. Erdem, F. Aliçavus, T. Akin, and Ç. Kanvermez, “CCD times of Minima of some eclipsing variables,” Information Bulletin on Variable Stars, vol. 5988, 2011. View at: Google Scholar
 N. Simmons, 2011, http://www.aavso.org/datadownload.
 R. Diethelm, “Timings of minima of eclipsing binaries,” Information Bulletin on Variable Stars, vol. 6011, p. 1, 2012. View at: Google Scholar
 R. Diethelm, “Timings of minima of eclipsing binaries,” Information Bulletin on Variable Stars, vol. 6042, 2013. View at: Google Scholar
 L. Kiss, G. Kaszás, G. Fürész, and J. Vinkó, “New times of minima and updated ephemerides of selected contact binaries,” Information Bulletin on Variable Stars, vol. 4681, p. 1, 1999. View at: Google Scholar
 G. Kaszas, J. Vinko, K. Szatmary et al., “Period variation and surface activity of the contact binary VW Cephei,” Astronomy & Astrophysics, vol. 331, pp. 231–243, 1998. View at: Google Scholar
 M. Wolf, P. Molik, K. Hornoch, and L. Sarounova, “Period changes in W UMatype eclipsing binaries: DK Cygni, V401 Cygni, AD Phoenicis and Y Sextantis,” Astronomy and Astrophysics Supplement Series, vol. 147, pp. 243–249, 2000. View at: Google Scholar
 T. Borkovits, M. M. Elkhateeb, S. Csizmadia et al., “Indirect evidence for short period magnetic cycles in W UMa stars Period analysis of five overcontact systems,” Astronomy & Astrophysics, vol. 441, no. 3, pp. 1087–1097, 2005. View at: Publisher Site  Google Scholar
 J. Kreiner, C. Kim, and I. Nha, An Atlas of O–C Diagrams of Eclipsing Binary Stars, Wydawnictwo Naukowe Akademii Pedagogicznej, Krakow, Poland, 2001.
 R. Szafraniec, “Minima of eclipsing variables observed in 1952,” Acta Astronomica. Series C, vol. 5, p. 51, 1953. View at: Google Scholar
 J. Hubscher, F. Agerer, and E. Wunder, “Beobachtungsergebnisse der Berliner Arbeitsgemeinschaft für Veränderliche Sterne e.V,” Bundesdeutsche Arbeitsgemeinschaft für Veränderliche Sterne Mitteilungen, vol. 59, 1990. View at: Google Scholar
 M. Baldwin and G. Samolyk, AAVSO, No. 8, 2003.
 T. Borkovits, I. B. Bíró, T. Hegedüs et al., “New times of Minima of eclipsing binary systems,” Information Bulletin on Variable Stars, vol. 5313, 2002. View at: Google Scholar
 H. Gerner, “Recent minima of 155 eclipsing binary stars,” The Journal of the American Association of Variable Star Observers, vol. 36, 2008. View at: Google Scholar
 G. Samolyk, “Recent minima of 154 eclipsing binary star,” Journal of the American Association of Variable Star Observers, vol. 37, no. 1, p. 44, 2009. View at: Google Scholar
 G. Samolyk, “Recent minima of 154 eclipsing binary stars,” Journal of the American Association of Variable Star Observers, vol. 37, pp. 44–51, 2009. View at: Google Scholar
 O. Y. Malkov, “Massluminosity relation of intermediatemass stars,” Monthly Notices of the Royal Astronomical Society, vol. 382, no. 3, pp. 1073–1086, 2007. View at: Publisher Site  Google Scholar
 A. Essam, S. M. Saad, M. I. Nouh, A. Dumitrescu, M. M. ElKhateeb, and A. Haroon, “Photometric and spectroscopic analysis of YY CrB,” New Astronomy, vol. 15, no. 2, pp. 227–233, 2010. View at: Publisher Site  Google Scholar
 M. M. Elkhateeb and M. I. Nouh, “Comprehensive photometric study of the eclipsing binary AW UMa,” Astrophysics and Space Science, vol. 352, no. 2, pp. 673–689, 2014. View at: Publisher Site  Google Scholar
 M. M. Elkhateeb and M. I. Nouh, “A holistic study of the eclipsing binary EQ Tau,” Journal of Physics and Astronomy Research, vol. 1, no. 3, p. 15, 2014. View at: Google Scholar
 Y. Yang, S. Qian, L. Zhu, J. He, and J. Yuan, “Photometric investigations of three shortperiod binary systems: GSC 07630572, RR centauri, and ϵ Coronae australis,” Publications of the Astronomical Society of Japan, vol. 57, no. 6, pp. 983–993, 2005. View at: Publisher Site  Google Scholar
 O. L. Degirmenci, “Photometric analysis of the W UMa type binary V566 ophiuchi,” Information Bulletin on Variable Stars, vol. 5726, pp. 1–4, 2006. View at: Google Scholar
 R. Nelson, 2009, http://members.shaw.ca/bob.nelson/software1.htm.
 R. Wilson and E. Devinney, “Realization of accurate closebinary light curves: application to MR cygni,” The Astrophysical Journal, vol. 166, p. 605, 1971. View at: Publisher Site  Google Scholar
 R. E. Wilson, “Accuracy and efficiency in the binary star reflection effect,” The Astrophysical Journal, vol. 356, no. 2, pp. 613–622, 1990. View at: Publisher Site  Google Scholar
 J. Kallarath and E. Milone, “A study of the O'Connell effect in the light curves of eclipsing binaries,” The Astrophysical Journal Supplement Series, vol. 55, pp. 571–584, 1984. View at: Publisher Site  Google Scholar
 L. Lucy, “Gravitydarkening for stars with convective envelopes,” Zeitschrift für Astrophysik, vol. 65, p. 89, 1967. View at: Google Scholar
 R. Ruciński, “The proximity effects in close binary systems. II. The bolometric reflection effect for stars with deep convective envelopes,” Acta Astronomica, vol. 19, p. 245, 1969. View at: Google Scholar
 W. van Hamme, “New limbdarkening coefficients for modeling binary star light curves,” The Astronomical Journal, vol. 106, no. 5, pp. 2096–2117, 1993. View at: Publisher Site  Google Scholar
 D. Bradstreet and D. Steelman, “Binary maker 3.0—an interactive graphicsbased light curve synthesis program written in java,” Bulletin of the American Astronomical Society, vol. 34, p. 1224, 2002. View at: Google Scholar
 L. Girardi, A. Bressan, G. Bertelli, and C. Chiosi, “Evolutionary tracks and isochrones for low and intermediatemass stars: from 0.15 to 7 ${M}_{\odot}$, and from $Z=0.0004$ to 0.03,” Astronomy and Astrophysics Supplement Series, vol. 141, pp. 371–383, 2000. View at: Google Scholar
Copyright
Copyright © 2015 M. M. Elkhateeb et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.