Advanced
The First Comprehensive Photometric Study of the Neglected Binary System V345 Cassiopeiae
The First Comprehensive Photometric Study of the Neglected Binary System V345 Cassiopeiae
Journal of Astronomy and Space Sciences. 2013. Dec, 30(4): 213-221
Copyright ©2013, The Korean Space Science Society
This is an open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/3.0/) which premits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.
  • Received : November 11, 2013
  • Accepted : November 11, 2013
  • Published : December 15, 2013
Download
PDF
e-PUB
PubReader
PPT
Export by style
Share
Article
Author
Metrics
Cited by
TagCloud
About the Authors
Min-Ji Jeong
Chun-Hwey Kim
kimch@chungbuk.ac.kr
Abstract
We present the first BVRI CCD photometric observations of V345 Cas made on 22 nights during the observing seasons in 2007-2008. Our light curves, resembling those of EB-type binary stars, are analyzed with the 2003 version of the Wilson-Devinney binary code. The photometric study shows that V345 Cas is a typical semi-detached binary system with the primary star being about Δ T = 2400 K hotter, two times more massive, but only 17% larger than the secondary star being filled with it's Roche-lobe. The orbit has a relatively large inclination of about 88°. A cool spot on the secondary component is modeled to explain a small light curve asymmetry. Absolute dimensions and related radiometric parameters of the eclipsing pair are calculated and their evolutionary states are discussed with the HR diagrams of mass-radius and temperature-luminosity. A period analysis of all available times of minima, including our measurements, indicates that the orbital period may vary in a cyclical way, unfortunately the secondary period for the variation can not be uniquely determined because of lack of present timing data.
Keywords
1. INTRODUCTION
The light variability of V345 Cas (AN 190.1943, 2MASS J23083986+5406545, USNO B1.0 1441-0442463) was discovered by Hoffmeister (1943). Huth (1956) classified the system as a β Lyr-type eclipsing binary showing the light-variation of 13. m 1 to 14. m 2. He also determined the first light elements as Min I = HJD 2428835.307 + 0. d 6887589 E . Since his study numerous times of minimum light were published by many investigators (see Table 3 for the references). By using all the times of minimum light available to them, Kreiner et al. (2001) gave an ephemeris ( Min I = HJD 2428835.3107 + 0. d 68875832E). Neither any light curves nor radial velocity curves have been published so far.
As reviewed above, V345 Cas is a neglected binary system whose photometric properties were less known. Thus we have included V345 Cas in our target stars to study photometric properties of close binary stars. In this paper we present, for the first time, two sets of BVRI light curves of V345 Cas secured in 2007-2008, and analyze them to obtain the global photometric properties as well as the dynamical behavior of the orbital period. Also we estimate the absolute dimensions of V345 Cas, and finally we discuss the evolutionary status of two components in the system.
2. OBSERVATIONS AND LIGHT CURVES
- 2.1 Observations
Our photometric observations were made on 22 nights in the 2007-2008 observing seasons with a FLI 2K CCD camera and a Johnson's standard BVRI filter set attached to a 1-m reflecting telescope at Mt. Lemon Astronomical Observatory (LOAO) in Arizona, USA. The telescope was controlled remotely at the Korea Astronomy and Space Science Institute (KASI). The CCD camera, which is thermoelectrically cooled during the observations, has the field
The coordinates and magnitudes of variable, comparison and check stars.††All data were extracted from the USNO-B ) catalog (see Monet et al. 2003)
Lager Image
The coordinates and magnitudes of variable, comparison and check stars. All data were extracted from the USNO-B ) catalog (see Monet et al. 2003)
of view of about 20.'5 × 20.'5 at the f/7.5 Cassegrain focus of the telescope. Each of the frames was 2 × 2 binned and differently exposed from 30 s to 140 s according to the filters used and weather conditions of each nights. All measured frames were corrected with bias, dark, and flat-field images and preprocessed in the IRAF/DIPHO S/W package. The IRAF Apphot was used to determine instrumental magnitudes. The USNO B1.0 1440-0439477 and USNO B1.0 1440-0439661 near V345 Cas were chosen as a comparison and check stars, respectively. Table 1 lists the coordinates and photometric characteristics of V345 Cas, comparison and check stars and Fig. 1 shows one of the observed CCD frames where the program stars are indicated. Differential magnitudes of the check star relative to the comparison star show that our mean photometric accuracy during the observations is about 0.01 mag in all BVRI bandpasses. A total of 862 measurements ( B : 230, V : 230, R : 200, and I : 202) in 2007 and 1020 ( B : 253, V : 255, R : 261, and I : 251) in 2008 were secured. All individual observations are available on request from the second author. Two sets of complete BVRI
Lager Image
The finding chart of V345 Cas. The field of view is about 22.5' × 22.5'. 'V', 'C' and 'K' denote V345 Cas, the comparison (USNO B1.0 1440-0439477) and check (USNO B1.0 1440-0439661) stars, respectively. East is left and north is up.
light and color curves were shown in Fig. 2 where the phase was calculated with a new light elements as follows;
Lager Image
Eq. (1) was determined with a least-squares method by using recent 21 times of minima of which the (O-C) residuals show a linear pattern to fit well to our observations.
- 2.2 Light curves
As shown in Fig. 2 , two sets of BVRI light curves of V345 Cas are consistent each other and have deep eclipses in the primary minimum and shallow in the secondary, similar to other short-period Algol-type eclipsing binaries (e.g., CL Aur (Lee et al. 2010), SW Lyn (Kim et al. 2010) and RU UMi (Lee et al. 2008)). The depth of the primary eclipse decreases with wavelength, being equal to about 1. m 42, 1. m 30, 1. m 23 and 1. m 11 for B, V, R , and I , respectively. However, the depth of the secondary eclipse increases with wavelength, being equal to about 0. m 23, 0. m 25, 0. m 30 and 0. m 33 for B, V, R and I , respectively. Another interesting feature is that all color curves of B-V, V-R and R-I show redder during the primary eclipse and bluer during the secondary eclipse. These two phenomena (wavelength-dependence of eclipse depth and opposite color variations at the primary and secondary eclipses) would be strong indications that 1) primary eclipse is an occultation (Sipahi 2005) and 2) temperature difference between two components is relatively large. In addition it is seen that there is a small asymmetry between the descending and the ascending branches during secondary eclipse. The descending branch is more sharp than the ascending one and the asymmetry is more pronounced as the wavelength is shorter. Furthermore the ascending branches in the 2007 light curves are lower than those in the 2008 light curves, indicating that the asymmetry is time-variable.
3. LIGHT CURVE SYNTHESIS
As explained in the previous section, there is a small difference due to light curve asymmetry and its variation between the 2007 and 2008 light curves. Therefore we solved two sets of filtered light curves separately with the 2003 version of the Wilson-Devinney binary code (hereafter WD, Wilson&Devinney 1971, van Hamme&Wilson 2003). Unfortunately, we have no information on the temperature of the primary star ( T 1 ) as well as the mass ratio ( q = M 2 / M 1 ),
Lager Image
The observed two sets of BVRI light and color curves of V345 Cas. Left: the 2007 light curves, right: the 2008ones. The light curves of comparison minus check stars are plotted in the top of both panels. The solid curves denote thetheoretical model light and color curves.
although these parameters are very essential to any light curve synthesis based on a Roche geometry (such as the WD code) which is sensitive to them. In order to estimate the temperature of the primary star which is eclipsed at primary eclipse, we collected information on magnitudes measured at various bandpasses from some catalogues. From the magnitude data we calculated color indices and estimated temperatures corresponding to each of the color indices by using the Allen's tables (Cox 2000). These are listed in Table 2 . The resultant temperatures of the V345 Cas system are distributed in the wide range from 6548 K to 17273 K. So the temperatures were averaged to be 7366 K when the unreasonably large temperatures of 11042 K from GSC2.3 catalog and 17273 K from NOMAD-1 were excluded. It may be hasty to use the averaged temperature of 7366 K as the temperature of the primary star because it was obtained from too much widely distributed temperatures. Thus, by using our 2008 light curves, a temperature-searching (hereafter T -search) of the primary star with the WD code was intensively made for both of cases that the primary star
Color indices and temperatures of V345 Cas extracted from various catalogs.†http://vizier.u-strasbg.fr/
Lager Image
Color indices and temperatures of V345 Cas extracted from various catalogs. http://vizier.u-strasbg.fr/
has a radiative envelope or a convective envelope, because the averaged temperature is near the The result of our T -search for a series of models with varying temperature was shown in the upper panel of Fig. 3 where open squares with lines represent the case with a radiative envelope and open circles a convective envelope. For the radiative models the weighted sum of the squared residuals, ΣW (O-C) 2 (hereafter Σ ) has a global minimum at T 1 = 7400 K while no convergence of the Σ was achieved for the convective cases, indicating that the primary star has a radiative envelope. It is very impressive to note that the temperature obtained with the T -search is practically the same as that averaged from the catalogues. As a next step, after fixing the temperature of the primary star as T 1 = 7400 K, we made an intensive q -search which shows to have a global convergence at q = 0.5, as shown in the lower panel of Fig. 3 . In these procedures, the differential correction (DC) program of the WD code started from mode 2 (detached system) but always converged to mode 5 (semi-detached system with the secondary component being filled with it's inner Roche lobe) during the computation. Therefore subsequent light curve syntheses were made with mode 5, employing the method of multiple subsets (Wilson&Biermann 1976) to ensure the stability of our solutions. Gravity darkening and albedo coefficients were adopted as g 1 = 1.00 and g 2 = 0.32 (Lucy 1967) and A 1 = 1.0 and A 2 = 0.5 (Rucinski 1969), respectively. Here, the subscript"1"denotes the hotter and massive primary star. The limb-darkening (LD) coefficients from the logarithmic law were
Lager Image
(a) The T-search diagram of V345 Cas. The dash line indicates the global minimum of ΣW(O-C)2 at T = 7400 K for the case with a radiative envelope. (b) The q-search diagram of V345 Cas. An arrow indicates the global minimum at q = 0.5.
The photometric solutions of V345 Cas.
Lager Image
The photometric solutions of V345 Cas.
Absolute dimensions of V345 Cas.
Lager Image
Absolute dimensions of V345 Cas.
computed at each iteration from the LD tables given by van Hamme (1993). The options of model atmosphere (Kurucz 1993) and synchronous rotation in the WD program were also adopted. The temperature was fixed as T 1 = 7400K and the mass ratio of q = 0.5 was used as an initial value. After the photometric solutions of the 2007 and 2008 light curves were separately obtained, a spot model (cool or hot spots on each of two components) was considered to model a small asymmetry appeared in our light curves. Only a cool spot on the secondary star has satisfactory fits to both 2007 and 2008 light curves.
At this moment a third-light was also included in the WD code, but the resultant values for the third-light in the BVRI bandpasses were negligibly small within their assigned errors. Our final photometric solutions of V345 Cas for the 2007 and 2008 light curves were listed in the second and third columns of Table 3 , respectively, where two solutions are practically consistent each other. The solid lines in Fig. 2 represent the model light and color curves drawn with the solution parameters in Table 3 . The weighted mean values of parameters of two solutions were listed in the fourth column of Table 3 and adopted as our final photometric solution of V345 Cas.
Our solution shows that V345 Cas is a typical semi-detached binary system with the cool secondary component filling it's inner Roche lobe, the mass ratio of q = 0.49, the inclination of i = 88. o 2 and the temperature difference of ΔT = 2440 K between two components. The secondary star has a cool spot with a large radius of about 21 o ~ 24 o near the neck connecting two stars, which makes the light curves slightly asymmetrical and variable near 0.5 phase.
Using the mean solution parameters in Table 3 , the Harmanec's (1988) main-sequence relationship between temperature and mass and some well-established relations (Flower 1996, Hilditch 2001, Torres 2010), the absolute dimensions and related radiative parameters of V345 Cas were calculated and listed in Table 4 . The luminosities and the bolometric magnitudes were computed by adopting T =5780 K and M bol⊙ = +4. m 69 given by Popper (1980). Fig. 4 shows the Roche geometry of V345 Cas where the spectral types of two components are deduced from the table of Cox (2000).
Lager Image
The Roche geometry of V345 Cas.
Lager Image
The eclipse timing diagram of V345 Cas. The solid and dashed curves represent the light-time orbits with two possible periods of 29.5 and 40.0 years, respectively.
4. PERIOD STUDY
From our LOAO observations and Super Wide Angle Search for Planets (Butters et al. 2010, Street et al. 2003) measurements, new 14 times of minimum light were determined with the method of Kwee&van Woerden (1956). A total of 69 times of minima, including ours, were collected from the recent database of Kreiner et al. (2001) and from recent literatures. These minima were listed in Table 5 and analyzed to figure out the dynamical behavior of V345 Cas for the first time. An eclipse timing diagram constructed with the Kreiner et al.'s ephemeris was shown in Fig. 5 which contains a timing history of V345 Cas over 66 years. The period change prior to 1990 is not impressive due to large scatter of the residuals of earlier non-CCD minima. However, it is clearly seen that the period, after 1990, is changing in a cyclical pattern, owing to recent high precision CCD minima. Under the assumption that the period change is due to the light-time effect (LITE) by a third-body, all the timings were fitted to an LITE ephemeris as follows:
Lager Image
Observed times of minimum light for V345 Cas.*P: plate, PG: photographic, VI: visual : CCD: charge-coupled device
Lager Image
Observed times of minimum light for V345 Cas. *P: plate, PG: photographic, VI: visual : CCD: charge-coupled device
where τ 3 is the LITE term of which the details were given by Irwin (1952). The seven unknown parameters in Eq. (3)
The solution of Eq. (2) for V345 Cas.
Lager Image
The solution of Eq. (2) for V345 Cas.
(Irwin 1959) were determined with a Levenberg-Marquardt method (Press et al. 1992) which was frequently used by the authors (e.g., Kim et al. 2010, Jeong&Kim 2011, Kim&Jeong 2012). The result shows that two solutions rather than an unique solution are possible in parallel mainly because the time-span of the CCD timing data with high precision is short while non-CCD minima with longer time-span have large scatters. As listed in Table 6 , the one LITE solution has P 3 = 29.5 yr, e = 0.61, and K = 0.0055 d and the other P 3 = 40.0 yr, e = 0.51, and K = 0.0050 d. The 1- σ value of the solution for the shorter period is slightly smaller than that for the longer one, however, the difference is so small that we can not decern, at the moment, which one between two solutions is more realistic. More timings would be needed to resolve the ambiguity in the secondary periods for the cyclical period change. Solid and dashed curves in Fig. 5 denote the shorter and longer LITE solutions, respectively. When two cyclical period changes suggested here are assumed to be caused by the Applegate (1992) mechanism in the convective region of the late-type secondary star, the resultant Applegate parameters are quite small (e.g., 2 × 10 -2 Ls / Ls <4 × 10 -2 ), indicating that the Applegate model could not work well.
5. DISCUSSION AND CONCLUSION
In order to investigate the evolutionary status of each of components of V345 Cas, mass-radius and mass-luminosity diagrams for the components of V345 Cas together with six other short-period Algol binaries were
Absolute dimensions of some short-period Algol binary stars.*1: Lee et al. (2010), 2: Lee et al. (2009), 3: Kim et al. (2010), 4: Sürgit et al. (2010), 5: İbanoğlu et al. (2000), 6: Lee et al. (2008)
Lager Image
Absolute dimensions of some short-period Algol binary stars. *1: Lee et al. (2010), 2: Lee et al. (2009), 3: Kim et al. (2010), 4: Sürgit et al. (2010), 5: İbanoğlu et al. (2000), 6: Lee et al. (2008)
Lager Image
The HR diagrams of some short-period Algol binary stars. (a) The mass-radius diagram. (b) The temperature-luminosity diagram. The open and filled symbols represent the primary and secondary components. Two components of V345 Cas were designated as star symbols. The thick solid and dotted lines represent the beginning (ZAMS) and ending (TAMS) phases of core hydrogen burning, respectively, which were calculated with the Giradi et al.'s (2000) evolutionary stellar models for Population I stars (X = 0.708, Y = 0.273, Z = 0.019). Two evolutionary tracks corresponding to 1.6 and 0.8 M were also drawn as short-dashed lines in the temperature-luminosity diagram.
drawn in the upper and lower parts of Fig. 6 , respectively. The information on stellar parameters of other binaries was listed in Table 7 . In the figure, open and filled star symbols represent the primary and secondary components of V345 Cas,ely, while the primary and secondary components of other binary stars are designated as open and filled circles,ely. The thick solid and dotted lines represent the beginning (ZAMS) and ending (TAMS) phases of core hydrogen burning,ely, which were calculated with the Girardi et al.'s (2000) evolutionary stellar models for Population I stars ( X = 0.708, Y = 0.273, Z = 0.019). Two evolutionary tracks corresponding to 1.6 and 0.8 M were also drawn in the mass-luminosity diagram. As shown in the figures, the components of V345 Cas are located on the positions that the primary is in slightly upper region from the ZAMS while the secondary is located in the sub-giant domain, which is similar to other short-period Algols, indicating that the primary star is a A7 main-sequence star but the secondary is an evolved K3 sub-giant star. The ratio of the primary star's radius to it's Roche volume radius is about 0.85. Therefore, V345 Cas is not a near-contact binary and would be expected to evolve into the near-contact stage.
Although the photometric properties of V345 Cas were comprehensively understood by our study, however the details are far from the reality. The lack of the present timing data can not constrain any secular period change which would be expected from the semi-detached configuration of the system and resulted in two possible secondary periods for a cyclical period change. More photometric and spectroscopic observations are needed to resolve what we have suggested as some problems and further to improve the understanding of the system.
Acknowledgements
We would like to thank the operating staff of LOAO 1m telescope of KASI. We are also thankful for the reviewer’s constructive comments which improve our original manuscript. This research has been made using the SIMBAD and VizieR databases operated by the Centre de Donnees Astronomiques (Strasbourg). This publication makes use of data products from the Two Micron All Sky Survey, which is a joint project of the University of Massachusetts and the Infrared Processing and Analysis Center/California Institute of Technology, funded by the National Aeronautics and Space Administration and the National Science Foundation. We acknowledge the WASP consortium which comprises of the University of Cambridge, Keele University, University of Leicester, The Open University, the Queen's University Belfast, St. Andrews University and the Isaac Newton Group. Funding for WASP comes from the consortium universities and from the UK's Science and Technology Facilities Council. This work was supported by the research grant of the Chungbuk National University in 2011 and by the Korea Research Foundation (KRF) grant funded by the Korea government (NRF-2012R1A1A4A01012467).
References
Applegate JH (1992) A mechanism for orbital period modulation in close binaries. ApJ 385 621 - 629
Brát L , Zejda M , Svoboda P (2007) B.R.N.O. Contributions 34 OEJV 74 1 - 75
Butters OW , West RG , Anderson DR (2010) Collier Cameron A, Clarkson WI, et al., The first WASP public data release A&A 1-4 520 L10 -
Cox AN Allen's Astrophysical Quantities fourth edition Springer New York ed., A. N. Cox 151, 388 -
Diethelm R (1999) 152. List of Minima of Eclipsing Binaries BBSAG Bull 119 1 - 10
Diethelm R (2002) 160. List of Minima of Eclipsing Binaries BBSAG Bull. 127 1 - 11
Diethelm R (2003) 162-nd list of minima timings of eclipsing binaries by BBSAG observers IBVS 5438 1 - 11
Diethelm R (2010) Timings of Minima of Eclipsing Binaries IBVS 5920 1 - 7
Diethelm R (2011) Timings of Minima of Eclipsing Binaries. IBVS 5960 1 - 11
Diethelm R (2012) Timings of Minima of Eclipsing Binaries IBVS 6011 1 - 13
Flower PJ (1996) Transformations from Theoretical Hertzsprung-Russell Diagrams to Color-Magnitude Diagrams: Effective Temperatures, B-V Colors, and Bolometric Corrections ApJ 469 355 - 365
Götz W , Huth H , Hoffmeister C (1957) Die Veraenderlichen Sterne der nordlichen Milchstrasse Teil IX. VeSon 4 123 - 238
Girardi L , Bressan A , Bertelli G , Chiosi C (2000) Evolutionary tracks and isochrones for low- and intermediate-mass stars: From 0.15 to 7 M⊙, and from Z = 0.0004 to 0.03 A&AS 141 371 - 383
Harmanec P (1988) Stellar masses and radii based on modern binary data. Bull. Astron. Inst. Czech. 39 329 - 345
Häussler K (1990) Veranderliche im Feld 23h09m52o5 (supplement). VeSon 10 374 - 388
Hilditch RW 2001 An Introduction to Close Binary Stars Cambridge University Press Cambridge Chap. 6.
Hoffmeister C (1943) 213 neue Veranderliche. AN 274 36 - 40
Huth H (1956) Herausgegeben von der sternwarte Sonneberg der Deutschen Akademie der Wissenschaften MVS 229
Hübscher J (2011) BAV-results of observations-photoelectric minima of selected eclipsing binaries and maxima of pulsating stars. IBVS 5984 1 - 16
Hübscher J , Lehmann PB (2012) BAV-Results of Observations-Photoelectric Minima of Selected Eclipsing Binaries and Maxima of Pulsationg Stars IBVS 6026 1 - 23
Hübscher J , Lehmann PB , Monninger G , Steinbach H-M , Walter F (2010) BAV-results of observations-photoelectric minima of selected eclipsing binaries and maxima of pulsating stars IBVS 5941 1 - 16
Hübscher J , Lehmann PB , Walter F (2012) BAV-Results of Observations-Photoelectric Minima of Selected Eclipsing Binaries and Maxima of Pulsationg Stars. IBVS 6010 1 - 22
HübscherHubscher J , Steinbach H-M , Walter F (2008) BAV-Results of Observations-Photoelectric Minima of Selected Eclipsing Binaries and Maxima of Pulsationg Stars. IBVS 5830 1 - 8
Hübscher J , Steinbach H-M , Walter F (2009a) BAV-Results of Observations-Photoelectric Minima of Selected Eclipsing Binaries and Maxima of Pulsationg Stars. IBVS 5874 1 - 13
Hübscher J , Steinbach H-M , Walter F (2009b) BAV-results of observations-photoelectric minima of selected eclipsing binaries and maxima of pulsating stars. IBVS 5889 1 - 13
Irwin JB (1952) The determination of a light-time orbit. ApJ 116 211 - 217
Irwin JB (1959) Standard light-time curves. AJ 64 149 - 155
İbanoğlu C , Çakirli Ö , Değirmenci Ö , Saygan S , Ulaş B (2000) The BV light and it O-C curves analyses of the triple system V505 Sagittarii. A&A 354 188 - 192
Jeong JH , Kim C-H (2011) WZ Cephei: a dynamically active W UMa-type binary star. JASS 28 163 - 172
Kim C-H , Jeong JH (2012) V700 Cygni: A Dynamically Active W UMa-type Binary Star II. JASS 29 151 - 161
Kim C-H , Kim H-I , Yoon TS , Han W , Lee JW (2010) SW Lyncis-advances and questions. JASS 27 263 - 278
Kreiner JM , Kim C-H , Nha IS 2001 An atlas of O-C diagrams of eclipsing binary stars Wydawnictwo Naukowe Akademii Pedagogicznej Krakow
Kurucz RL 1993 in Light Curve Modeling of Eclipsing Binary Stars, ed., E. F. Milone Springer-Verlag New York 93 -
Kwee KK , van Woerden H (1956) method for computing accurately the epoch of minimum of an eclipsing variable BAN 12 327 - 330
Lee JW , Kim C-H , Kim DH , Kim S-L , Lee CU (2010) The First Comprehensive Photometric Study of the Algol-type System CL Aurigae. AJ 139 2669 - 2678
Lee JW , Kim C-H , Kim S-L , Lee C-U , Han W (2008) A Long-term Photometric Study of the Near-contact Binary RU Ursae Minoris. PASP 120 720 - 729
Lee JW , Kim SL , Lee CU , Kim HI , Park JH (2009) Photometric Properties of the Near-Contact Binary GW Geminorum. PASP 121 114 - 110
Locher K (1981a) 89th List of Minima of Eclipsing Binaries. BBSAG Bull. 56 1 - 7
Locher K (1981b) 90th List of Minima of Eclipsing Binaries. BBSAG Bull. 57 1 - 7
Locher K (1982a) 91th List of Minima of Eclipsing Binaries. BBSAG Bull. 58 1 - 7
Locher K (1982b) 95th List of Minima of Eclipsing Binaries. BBSAG Bull. 62 1 - 7
Lucy LB (1967) Gravity-Darkening for Stars with Convective Envelopes ZA 65 89 - 92
Mikulášek Z , Šilhán J , Zejda M (1992) Contributions Nicholas Copernicus Observatory. Brno 30 4 -
Monet DG , Levine SE , Canzian B , Ables HD , Bird AR (2003) The USNO-B catalog. AJ 125 984 - 993
Popper DM (1980) Stellar masses ARA&A 18 115 - 164
Press WH , Teukolsky SA , Vetterling WT , Flannery BP 1992 Numerical recipes in C 2nd ed. Cambridge University Press Cambridge Chap. 15.
Rucinski SM (1969) The Proximity Effects in Close Binary Systems. II. The Bolometric Reflection Effect for Stars with Deep Convective Envelopes AcA 19 245 - 255
Šafář J , Zejda M (2002) CCD times of minima of faint eclipsing binaries III IBVS 5263 1 - 9
Sipahi E (2005) 2004 UBVR Photometry of the Eclipsing Binary KR Cyg IBVS 5635 1 - 4
Street RA , Pollaco DL , Fitzsimmons A , Keenan FP , Horne K (2003) SuperWASP: Wide Angle Search for Planets. ASPC 294 405 - 408
Sürgit D , Erdem A , Budding E 2010 The Algol binary system KZ Pavonis revisited MNRAS 407 497 - 506
Torres G (2010) On the Use of Empirical Bolometric Corrections for Stars. AJ 140 1158 - 1162
van Hamme W (1993) New limb-darkening coefficients for modeling binary star light curves. AJ 106 2096 - 2117
van Hamme W , Wilson RE (2003) Stellar atmospheres in eclipsing binary models. ASPC 48 323 - 328
Wilson RE , Biermann P (1976) TX CANCRI - Which component is hotter A&A 48 349 - 357
Wilson RE , Devinney EJ (1971) Realization of accurate close-binary light curves: application to MR Cygni. ApJ 166 605 - 619
Zejda M (1995) Contributions Nicholas Copernicus Observatory. Brno 31 4 -