WZ Cephei: A Dynamically Active W UMa-Type Binary Star

Journal of Astronomy and Space Sciences.
2011.
Sep,
28(3):
163-172

This is an Open Access article distributed under the terms of theCreative Commons Attribution Non-Commercial License(http://creativecommons.org/licenses/by-nc/3.0/)which permits unrestrictednon-commercial use, distribution, and reproduction in any medium,provided the original work is properly cited.

- Received : June 27, 2011
- Accepted : June 30, 2011
- Published : September 15, 2011

Download

PDF

e-PUB

PubReader

PPT

Export by style

Share

Article

Metrics

Cited by

TagCloud

An intensive analysis of 185 timings of WZ Cep, including our new three timings, was made to understand the dynamical picture of this active W UMa-type binary. It was found that the orbital period of the system has complexly varied in two cyclical components superposed on a secularly downward parabola over about 80
^{y}
. The downward parabola, corresponding to a secular period decrease of -9.
^{d}
97 × 10
^{-8}
y
^{-1}
, is most probably produced by the action of both angular momentum loss (AML) due to magnetic braking and mass-transfer from the massive primary component to the secondary. The period decrease rate of -6.
^{d}
72 × 10
^{-8}
y
^{-1}
due to AML contributes about 67% to the observed period decrease. The mass flow of about 5.16 × 10
^{-8}
M
_{⊙}
y
^{-1}
from the primary to the secondary results the remaining 33% period decrease. Two cyclical components have an 11.
^{y}
8 period with amplitude of 0.
^{d}
0054 and a 41.
^{y}
3 period with amplitude of 0.
^{d}
0178. It is very interesting that there seems to be exactly in a commensurable 7:2 relation between their mean motions. As the possible causes, two rival interpretations (i.e., light-time effects (LTE) by additional bodies and the Applegate model) were considered.In the LTE interpretation, the minimum masses of 0.30
M
_{⊙}
for the shorter period and 0.49
M
_{⊙}
for the longer one were calculated. Their contributions to the total light were at most within 2%, if they were assumed to be main-sequence stars. If the LTE explanation is true for the WZ Cep system, the 7:2 relation found between their mean motions would be interpreted as a stable 7:2 orbit resonance produced by a long-term gravitational interaction between two tertiary bodies. In the Applegate model interpretation, the deduced model parameters indicate that the mechanism could work only in the primary star for both of the two period modulations, but could not in the secondary. However, we couldn't find any meaningful relation between the light variation and the period variability from the historical light curve data. At present, we prefer the interpretation of the mechanical perturbation from the third and fourth stars as the possible cause of two cycling period changes.
^{d}
260984. Nine years later, a photographic light curve was obtained by Balázs (1937), who correctly classified WZ Cep as a W UMa binary star with an orbital period of 0.
^{d}
41745. Detre (1940) also made photographic observations of the system, confirmed Balázs'result, and gave basic system parameters from the analysis of his light curve using Russell's classical model.The first photoelectric
BV
light curves were secured by Hoffmann (1984) who detected a strong O'Connell effect in his light curves, as the brightness in the 0.25 phase (Max I) was significantly higher than that in the 0.75 phase(Max II). The spectral type was assigned as F5 in the General Catalogue of Variable Stars (Kholopov et al. 1987). Hoffmann's light curves were reanalyzed by Kałużny(1986) with the Wilson Devinney binary model (WD;Wilson & Devinney 1971). His solutions showed that WZ Cep is a contact binary with the components of unequal surface temperatures (
ΔT
≈ 1,000 K). The O'Connell effect in the light curve was ascribed to a small hot spot on the less massive component near the neck between stars. These interpretations were questioned by Djurašević et al. (1998), who obtained new
BVR
light curves that were almost symmetrical. Their analysis of both their own light curves and Hoffmann's using their own Roche computer model (Djurašević 1992a, b) showed that the WZ Cep system is a shallow contact binary (
f
= 9~14%) with the components of slightly different temperatures (
ΔT
≈ 140 K) and with two large dark spots on the surface of the massive, hotter primary star. Recently, WZ Cep was revisited by Lee et al. (2008) and Zhu & Qian (2009). Lee et al. (2008) presented symmetrical
BVR
light curves in the observing season in 2005, while Zhu & Qian (2009) obtained an asymmetrical
V
light curve with a strong O'Connell effect. Both authors analyzed their light curves using the WD method and confirmed the photometric results of Djurašević et al. (1998). At the same time, Zhu & Qian (2009) made the first intensive period study of WZ Cep with all timings available to them, and found that the period variation consists of both a periodic and secularly decreasing terms. The periodic term has a period of 34.
^{y}
2 with amplitude of 0.
^{d}
013, while the secular period decrease rate is about -8.8 × 10
^{-8}
day/year.
Numerous timings just before and after the period study of Zhu & Qian (2009) have been published in the literature. In this paper, all published and newly observed timings have been reanalyzed. It is noticed that the ephemeris suggested by Zhu & Qian (2009) is not compatible with recent timing residuals. New ephemeris is derived with a detection of a secondary cyclic modulation of the O-C residuals by a period of about 11.8 year. A discussion of the implications on the existing working hypotheses follows.
^{m}
02. From our observations, three primary times of minimum light were determined using the conventional Kwee & van Woerden (1956) method, and are listed in
Table 1
.
O-C
) residuals of all timings were calculated with the linear ephemeris of Kreiner et al. (2001), as follows:
The (
O-C
_{1}
) diagram is shown in
Fig. 1
, where the times of minimum are marked by assorted symbols that differ in size and shape according to observational method and type of eclipse. The (
O-C
_{1}
) residuals were listed in the sixth column of
Table 1
. As can be seen in the figure, it
The (O-C_{1}) diagram of WZ Cep drawn together with the non-linear terms (solid and dashed lines) of Eq. (2). The recent residuals from Eq. (2) in the bottom are remarkably deviating from the non-linear terms of Eq. (2).
is clear that the variation patterns of the residuals before and after 1995
^{y}
are quite different from each other. The orbital period before 1995
^{y}
seemed to be nearly constant,
Photoeletric and CCD times of minima of WZ Cep not listed in Zhu & Qian (2009) or published after their study. ^{a}Me: method, CCD and CC: charge coupled device, PE: photoelectric ^{b}Ty: type, I: primary eclipse, II: secondary eclipse. ^{c}The timing was listed at the web (http://var.astro.cz/ocgate/ocgate.php?star=wz+cep&submit=Submit&lang=en). submit=Submit&lang=en).
although there is a large gap of about 25 years between 1939
^{y}
and 1964
^{y}
. However, the period from 1995
^{y}
to the present has dramatically decreased in a continuous way.
This fact was noticed by Zhu & Qian (2009), who adopted a quadratic plus sine ephemeris to explain the (
O-C
_{1}
) residuals constructed with the timings available to their time, as follows:
In Eq. (2) the phase angle among the arguments in sine function should be read as +103.
^{o}
5 rather than -76.
^{o}
5, which must be a typo. The value of +103.
^{o}
5 was calculated by us with the (
O-C
)
_{1}
and (
O-C
)
_{2}
values listed in the fifth and six columns of
Table 2
in the paper of Zhu & Qian (2009). The dashed and solid lines in
Fig. 1
denote the linear plus quadratic and full terms in Eq. (2), respectively. In the bottom of
Fig. 1
, the residuals from Eq. (2) where the recent timings since 1995y have been clearly deviating from the ephemeris are plotted, implying that Eq. (2) should be revised and/or the period change of WZ Cep may be more complicated. The standard deviation (σ) for the photoelectric and CCD residuals was calculated to be σ = ±0.
^{d}
0029. Such a value cannot be neglected because of the accuracy of modern photoelectric and CCD timings (Kim et al. 2005).
As a first step in understanding the general behavior of period change of the WZ Cep system, we tried to fit all (
O-C
_{1}
) residuals to a quadratic plus sine ephemeris by using the Levenberg-Marquardt method (Press et al. 1992). The final solution was quickly converged as follows:
In this and subsequent calculations, we assigned different weights to each of the data according to the observational method: 10 to photoelectric or CCD observations,3 to photographic, and 1 to visual data, which were the same weights as those used by Zhu & Qian (2009). The linear plus quadratic and full terms in Eq. (3) were drawn with the (
O-C
_{1}
) residuals as the dashed and solid lines in
Fig. 2
, respectively. The residuals from full terms of Eq. (3) were plotted in the bottom of
Fig. 2
. As shown in
Fig. 2
, Eq. (3) gives a more satisfactory fit to the observed timings than Eq. (2). From the argument of 0.o0089 in the sine term of Eq. (3), the cyclic period is calculated as 46.y4 (±0.
^{y}
2), which is 12.
^{y}
2 longer than that of Eq. (2) given by Zhu & Qian (2009). It is interesting to note that the residuals from Eq. (3) show an oscillatory pattern, especially in the photoelectric and CCD timings, with small amplitude of about 0.
^{d}
005 and with a short period of about 10
^{y}
.
The (O-C_{1}) diagram of WZ Cep drawn together with the non-linear terms (solid and dashed lines) of Eq. (3). The recent residuals from Eq. (3) in the bottom clearly show an oscillatory pattern with a short period of about 10^{y} and a small amplitude of about 0.^{d}004.
Final solution of Eq. (4) and related parameters.
To find the variation characteristics of the small oscillation in
Fig. 2
, several non-linear terms in addition to a linear plus quadratic ephemeris were considered, such as: 1) a light-time effect (LTE), 2) two sine curves, 3) a sine curve plus an LTE, 4) two LTEs. Next, attempts were made to fit all timings to each of the four ephemeris models. The result showed that the fits with the ephemeris models of both 3) and 4) were superior to those with 1) and 2). In addition, the solution with the ephemeris 4) gave zero eccentricity for the LTE orbit, with a shorter period and smaller amplitude in two LTE orbits. This means that the ephemeris 3) is best for our purpose. The ephemeris 3) is expressed as:
where τ
_{3}
is the light-time due to the assumed third body, and includes five parameters (
a
_{12}
sin
i
_{12}
,
e
_{12}
,
ω
_{12}
,
P
_{12}
,
T
_{12}
) which are the orbital parameters of the eclipsing pair around the mass center of the triple system. The parametric and differential forms of the orbital elements for the light-time orbit were taken from Irwin (1952, 1959). The Levenberg-Marquardt method was used again to solve Eq. (4) simultaneously. In this case, there are 11 parameters to be adjusted in Eq. (4). The calculations converged quickly to yield the solution listed in
Table 2
. The quality of the final solution represented by Eq. (4) is very satisfactory, and may be judged from
Fig. 3
, in which the (
O-C
_{1}
) residuals were plotted with Eq. (1), and the dashed
The (O-C_{1}) diagram of WZ Cep drawn together with the non-linear terms (solid and dashed lines) in Table 2 .
The (O-C) diagram of WZ Cep phased with the shorter period of 11.^{y}8 in Table 2 .
and solid lines represent the linear plus quadratic and full terms of Eq. (4), respectively. The residuals from all terms of Eq. (4) were listed in the seventh column of
Table 1
, and were plotted in the bottom of
Fig. 3
. The standard deviation for the photoelectric and CCD residuals in the bottom of
Fig. 3
is σ = ±0.
^{d}
0002, which is excellently compatible with that of the modern photoelectric and CCD timing accuracy.
Figs. 4
and
5
show two cyclic (
O-C
) diagrams phased with each of their periods, in which the continuous lines were drawn with the solution parameters in
Table 2
. Summing up our analysis of all timings of WZ Cep, the period of the system has varied in two cyclical components superposed on a global downward parabola over about 80
^{y}
.
The (O-C) diagram of WZ Cep phased with the longer period of 41.^{y}3 in Table 2 .
^{d}
97 (±0.14) × 10
^{-8}
y
^{-1}
. In principle, there are two working mechanisms for the secular period decrease observed in a close binary star such as the WZ Cep system: 1) secular mass transfer from the more massive(“primary” hereinafter) to the less massive (“secondary”hereinafter) stars, 2) angular momentum loss (AML) by magnetic braking via a magnetically controlled stellar wind from a magnetically active star in a binary. WZ Cep has been known as a shallow contact system with a fill-out factor of about 9-24% (Djurašević et al. 1998, Lee et al. 2008, Zhu & Qian 2009) and with the massive primary hotter than the secondary. Moreover, the primary star has a large and variable dark spot(s) (Djurašević et al. 1998, Lee et al. 2008, Zhu & Qian 2009), indicating that it is a magnetically active star. This would mean that both a secular transfer of a gaseous mass stream from the primary to the secondary star and the AML from the magnetically active primary star are occurring concurrently in the WZ Cep system. On this basis, the observed secular period decrease is expressed as follows:
where the subscripts mt and mb denote mass transfer and magnetic braking, respectively. The period change by conservative mass transfer is expressed with the following well-known formula,
where
M
_{1}
and
M
_{2}
denote the masses of the primary and secondary stars, respectively. The period change due to the AML by magnetic braking is approximately expressed by the following formula (Bradstreet & Guinan 1994):
where
R
_{1}
and
R
_{2}
are the radii of the components in solar units,
q
is the mass ratio (
q
=
M
_{2}
/
M
_{1}
),
k
^{2}
is the gy-
Absolute parameters of WZ Cep.
ration constant typically ranging from 0.07 to 0.15 for solar-type stars, and
P
is the orbital period in days. Before solving Eqs. (5)-(7) above, we have to estimate the absolute parameters of WZ Cep, because no reliably published absolute parameters currently exist. Thus, the absolute parameters of WZ Cep were estimated with the photometric solution of Djurašević et al. (1998) and with Harmanec's (1988) empirical relation between mass and temperature.
Table 3
lists the resultant absolute parameters.By using the parameters in
Table. 2
and
3
and k2 = 0.1 (Webbink 1976), Eq. (7) was solved to determine the rate of decrease of the period due to magnetic braking as -6.
^{d}
72 × 10
^{-8}
y
^{-1}
, which contributes about 67% to the observed period decrease. Therefore, we note that magnetic braking may be a major mechanism contributing to the observed secular period change. Inserting these values into Eqs. (5) and (6), the mass transfer rate from the primary to the secondary star was calculated as 1.69 × 10
^{-8}
M
_{⊙}
y
^{-1}
. If the observed decreasing rate of period is assumed to occur purely by a continuous mass transfer, the mass flow rate is calculated with Eq. (6) to yield about 5.16 × 10
^{-8}
M
_{⊙}
y
^{-1}
, which is about three times larger than that obtained when the magnetic braking and the mass transfer are occurring simultaneously. From our discussions above, one may conclude that 1) the secular period decrease of WZ Cep has resulted from the combination of AML by magnetic braking and mass transfer from the primary to the secondary components, and 2) AML is about two times more effective than mass transfer in terms of its contribution to the observed secular period decrease.
O-C
_{1}
) residuals, shown in
Fig. 3
and listed in
Table 2
: One is a sine term with a semi-amplitude of 0.
^{d}
0054 and a period (
P
_{sine}
) of 11.
^{y}
8, while the other an LTE term with a semi-amplitude of 0.
^{d}
0178, a period (
P
_{LTE}
) of 41.
^{y}
3, and an eccentricity of 0.54. It is very interesting to note that there exists a commensurable 2:7 relation between
P
_{sine}
and
P
_{LTE}
; namely, 7
P
_{sine}
= 2
P
_{LTE}
.
There have been two rival theories to explain these cyclic period changes: One explanation is the LTEs due to additional third- and fourth-bodies in the WZ Cep system,and the other is the Applegate (1992) mechanism. Firstly, if the LTEs are assumed to be the mechanism to produce two cyclic period changes, then mass functions, masses and bolometric luminosities of two tertiary bodies were calculated for different inclinations by using the absolute parameters in
Table 3
. The results were listed in
Table 4
. In the calculation, the bolometric luminosities were obtained by using Harmanec's (1988) empirical formula for main-sequence stars. As listed in
Table 4
, the minimum masses (
i
' = 90
^{o}
) of 0.30
M
_{⊙}
for the third body and of 0.49
M
_{⊙}
for the fourth body would contribute very small fractions of about 0.1% and 2% lights, respectively, to the total luminosity (
l
_{t}
=
l
_{1}
+
l
_{2}
+
l
_{3}
+
l
_{4}
). Therefore, l3 and l4 would be hardly detectable in the light curve syntheses by the previous investigators (Djurašević et al. 1998, Lee et al. 2008, Zhu & Qian 2009).
Secondly, as an alternative explanation for the cyclical components of the period variability of WZ Cep, Applegate's (1992) model deserves attention. According to this theory, the orbital period can be modulated from changes in the distribution of angular momentum (and corresponding changes in stellar oblateness) as the magnetically active star undergoes changes in magnetic activity.
Basic parameters for the hypothetical third- and fourth-bodies in the WZ Cep system.
Applegate model parameters for the shorter and longer periods.
Using the modulation periods and amplitudes listed in
Table 2
and adopting the absolute dimensions provided in
Table 3
, the variations of parameters needed to change the orbital period of WZ Cep by specified amounts can be obtained from the formulae given by Applegate (1992). The calculations were made with the assumption that two cyclical changes of period are produced by both components. The resultant parameters were listed in
Table 5
, and show that the Applegate mechanism could operate only in the primary star for both two cyclic changes of period, but cannot in the secondary because the calculated luminosity change (
ΔL
_{rms}
) is about 11 and 3 times larger for the shorter and longer periods, respectively, than the values required by the model. Our determination that the primary may be the magnetically active star could be supported by Mullan's (1975) conclusion that it is the more massive components of W UMa systems that should preferentially manifest star spots, and thus, magnetic activity (Eaton 1986). Previous investigators (Djurašević et al. 1998, Lee et al. 2008, Zhu & Qian 2009) also invoked a cool spot (s) on the surface of the more massive primary component to explain the light curve asymmetry of WZ Cep.
Since the Applegate model predicts that the rms luminosity of the active primary star would be variable with the same period as that of the orbital period modulation,the relative lights (e.g., Max II [0.
^{p}
75] - Max I [0.
^{p}
25]) from four
V
bandpass light curves published so far were deduced, and their time variability was carefully investigated.But no meaningful relation between the light variation and the period variability was found, although the light asymmetry at a given epoch is remarkable and variable with different epochs. Two major reasons for this finding might be stressed: 1) The number of the light curves is only four, which is too few to test the Applegate model in detail, and 2) strong light variations produced by the localized spot activity (activities) may reside concurrently with the rms luminosity variations, or may overwhelm the rms luminosity itself. This negative conclusion somewhat strengthens the interpretation of the mechanical perturbation from the third- and fourth stars that was developed above.
^{d}
97 (±0.14) × 10
^{-8}
y
^{-1}
was calculated and interpreted as a combination of both AML due to magnetic braking and mass-transfer from the massive primary component to the secondary. The decreasing rate of period due to AML was obtained as -6.
^{d}
72 × 10
^{-8}
y
^{-1}
, which contributes about 67% to the observed period decrease. The remaining 33% contribution results from the mass flow of about 5.16 × 10
^{-8}
M
_{⊙}
y
^{-1}
from the massive primary component to the secondary. In this picture, the major driven contributor to the observed period decrease is AML by magnetic braking, which is about two times larger than mass-transfer. At the moment, WZ Cep may be considered as on an evolutionary track from the present shallow contact state of about 15% to a more deep contact state during a thermal time scale, which is expected by the magnetic braking theories (van't Veer 1979, Rucinski 1982, Maceroni & van't Veer 1991, Bradstreet& Guinan 1994, Stȩpień 1995, 2006, Demircan 1999) and the thermal relaxation oscillation theories (Flannery 1976, Lucy 1976, Robertson & Eggleton 1977).
In addition to the secular period decrease, WZ Cep has suffered from two cyclical period changes consisting of a sine curve with a shorter period of 11.
^{y}
8 and a semi-amplitude of 0.
^{d}
0054, and an LTE orbit with a longer period of 41.
^{y}
3, a semi-amplitude of 0.
^{d}
0178, and an eccentricity of 0.54. Very interestingly, there seem to exist in a commensurable 7:2 relation between their mean motion and a 33:10 relation between their amplitudes. As the possible causes for the cyclical period changes, two rival interpretations(LTE and Applegate model) were considered, and the parameters related with these were derived in
Table. 4
and 5, respectively. In the LTE interpretation, the obtained minimum masses of 0.30
M
_{⊙}
and 0.49
M
_{⊙}
for the third and fourth bodies, respectively, would contribute very small fractions of about 0.1% and 2% lights, respectively,to the total luminosity (
l
_{t}
=
l
_{1}
+
l
_{2}
+
l
_{3}
+
l
_{4}
), if they are assumed to be main-sequence stars. This explains why previous investigators (Djurašević et al. 1998, Lee et al. 2008, Zhu & Qian 2009) could not detect the third light (
l
_{3}
) in their light curve syntheses. In the LTEs explanation, the 7:2 relation found between their mean motions would be interpreted as a stable 7:2 orbit resonance produced by a long-term gravitational interaction between two tertiary bodies (Peale 1976, Kley et al. 2004)
In the interpretation with the Applegate model, the model parameters in
Table 5
indicate that the Applegate mechanism could work only in the primary component for both two cyclic changes of period, and cannot in the secondary. This means that the primary may be the magnetically active star, which is consistent with the studies of Mullan (1975) and previous investigators (Djurašević et al. 1998, Lee et al. 2008, Zhu & Qian 2009). However, we couldn't find any meaningful relation between the light variation and the period variability from the historical light curve data observed at four different epochs, although the light asymmetry at a given epoch is remarkable and variable with different epochs. We wait for future observations to resolve these unsolved and veiled matters. At the present, we somewhat prefer the interpretation of mechanical perturbation from the third- and fourth stars.

1. INTRODUCTION

The W UMa-type binary star, WZ Cep (AN 244.1928,2MASS J23222421+7254566, 1RXS J232216.6+725505)was discovered by Schneller (1928), as an RR Lyr type variable star with a period of 0.
2. OBSERVATIONS

Charge-coupled device (CCD) photometric observations of WZ Cep for the purpose of timing determination were made on the three nights of May 28 and October 9, 2005 and September 23, 2006 with the 35-cm reflector of the campus station of Chungbuk National University Observatory in Korea. The telescope was equipped with an SBIG ST-8 CCD imaging system, electrically cooled, with a 19' × 12' field of view. No filters were used in our observations. GSC 4486-1402 and GSC 4486-1352 were chosen as the comparison and check stars, respectively. Our comparison is the same one used by Lee et al. (2008). Camera exposure time ranged between 40 seconds and 60 seconds according to the quality of the night.
The instrumentation used and the reduction method for the raw CCD frames have been described in detail by Kim et al. (2006). The resultant standard errors of our observations in terms of comparison minus check star were about ±0.
3. PERIOD STUDY

To investigate the period variation of WZ Cep, a total of 185 (91 visual, 4 photographic, and 90 photoelectric and CCD) times of minimum light, including ours, have been collected from a modern database (Kreiner et al. 2001) and from the recent literature.
Table 1
lists only photoelectric and CCD minima, which were not compiled by Zhu & Qian (2009) or published after their study.
The (
Lager Image

Lager Image

Photoeletric and CCD times of minima of WZ Cep not listed in Zhu & Qian (2009) or published after their study.aMe: method, CCD and CC: charge coupled device, PE: photoelectricbTy: type, I: primary eclipse, II: secondary eclipse.cThe timing was listed at the web (http://var.astro.cz/ocgate/ocgate.php?star=wz+cep&submit=Submit&lang=en).submit=Submit&lang=en).

Lager Image

Lager Image

Lager Image

Lager Image

Final solution of Eq. (4) and related parameters.

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

4. THE SECULAR PERIOD DECREASE

From the coefficient (A) of the quadratic term in
Table 2
, the secularly decreasing rate of the period of WZ Cep is calculated to be -9.
Lager Image

Lager Image

Lager Image

Absolute parameters of WZ Cep.

Lager Image

5. TWO CYCLIC PERIOD CHANGES

In Section 3 we have decoupled two cyclic terms via Eq. (4) from the observed (
Basic parameters for the hypothetical third- and fourth-bodies in the WZ Cep system.

Lager Image

Applegate model parameters for the shorter and longer periods.

Lager Image

6. DISCUSSION AND CONCLUSION

We analyzed a total of 185 timings of WZ Cep, including our new three timings, to understand the dynamical picture of this very active W UMa binary system in detail. It was found that the orbital period of the system has varied in two cyclical components superposed on a global downward parabola over about 80y. Three decoupled components of the period variability were intensively investigated. For further investigations of these components,the absolute dimensions of WZ Cep were estimated with the photometric solutions of Djurašević et al. (1998) and with Harmanec's (1988) empirical relation between mass and temperature.
The secular period decrease of -9.
Acknowledgements

This research has made use of the Simbad database maintained at CDS, Strasbourg, France. This work was supported by a research grant from Chungbuk National University in 2009.

Balázs J
,
Cephei WZ
1937
Beob Zirk d A N
19
7 -

Baldwin ME
,
Samolyk G
2004
Observed minima timings of eclipsing binaries number 9
American Association of Variable Star Observers
Cambridge

Baldwin ME
,
Samolyk G
2007
Observed minima timings of eclipsing binaries number 12
American Association of Variable Star Observers
Cambridge

Bradstreet DH
,
Guinan EF
1994
Stellar mergers and acquisitions: the formation and evolution of W Ursae Majoris binaries
ASPC
56
228 -
243

Demircan O
1999
Dynamical evolution of the RS CVn-type binaries
TJPh
23
425 -
432

Detre L
1940
Das system WZ Cephei
Budapest Mitt
10
3 -

Diethelm R
2010
Timings of minima of eclipsing binaries
IBVS
5920
1 -

Diethelm R
2011
Timings of minima of eclipsing binaries
IBVS
5960
1 -

Eaton JA
1986
SW Lacertae W Urase Majoris YY Eridani and the prevalence of starspots in cool contact binaries
AcA
36
79 -
103

Harmanec P
1988
Stellar masses and radii based on modern binary data
BAICz
39
329 -
345

Hoffmann M
1984
Photoelectric observations of contact binary stars
VeBon
96
1 -
9

Kałużny J
1986
Contact binaries with components in poor thermal contact. II. WZ Cephei
AcA
36
105 -
111

Kholopov PN
1987
General catalogue of variable stars 4th ed. vol. III Constellations pavo-vulpecula
4th ed.
Nauka Publishing House
Moscow

Kim C-H
,
et al.
,
Lee C-U
,
et al.
,
Yoon Y-N
,
et al.
,
Park S-S
,
et al.
,
Kim D-H
,
et al.
2006
New CCD times of minima of eclipsing binary systems
et al.
5694
1 -

Kreiner JM
,
Kim C-H
,
Nha I-S
2001
An atlas of O-C diagrams of eclipsing binary stars
Wydawnictwo Naukowe Akademii Pedagogicznej
Kraków

Kwee KK
,
van Woerden H
1956
A method for computing accurately the epoch of minimum of an eclipsing variable
BAN
12
327 -

Maceroni C
,
van’t Veer F
1991
The galactic cluster NGC 188: W Ursae Majoris contact binaries as a clue to two separate bursts of star formation
A&A
248
430 -
434

Nagai K
2010
Visual and CCD minima of eclipsing binaries during 2009
Var Star Bull
50
1 -
10

Parimucha Š
,
et al.
,
Dubovský P
,
et al.
,
Baluďanský D
,
et al.
,
Pribulla T
,
et al.
,
Ham-bálek Ľ
,
et al.
2009
Minima times of selected eclipsing binaries
IBVS
et al.
5898
1 -

Parimucha Š
,
et al.
,
Dubovský P
,
et al.
,
Vaňko M
,
et al.
,
Pribulla T
,
et al.
,
Kudzej I
,
et al.
2011
Minima times of selected eclipsing binaries
IBVS
et al.
5980
1 -

Press WH
,
Teukolsky SA
,
Vetterling WT
,
Flannery BP
1992
Numerical recipes in C
2nd ed.
Cambridge University Press
Cambridge
Chap.15

Robertson JA
,
Eggleton PP
1977
The evolution of W Ursae Majoris systems
MNRAS
179
359 -
375

Rucinski SM
1982
Contact binaries: angular momentum loss in and out of contact
A&A
112
273 -
276

Samolyk G
2008a
Recent minima of 155 eclipsing binary stars
JAAVSO
36
171 -
185

Samolyk G
2008b
Recent minima of 184 eclipsing binary stars
JAAVSO
36
186 -
206

Samolyk G
51
Recent minima of 154 eclipsing binary stars
JAAVSO
37
44 -
51

Samolyk G
2010a
Recent minima of 185 eclipsing binary Stars
JAAVSO
38
1 -
10

Samolyk G
2010b
Recent minima of 161 eclipsing binary stars
JAAVSO
38
85 -
92

Stȩpień K
1995
Loss of angular momentum of cool close binaries and formation of contact systems
MNRAS
274
1019 -
1028

Stȩpień K
2006
The low-mass limit for total mass of W UMa-type binaries
AcA
5
347 -
364

van’t Veer F
1979
The angular momentum controlled evolution of solar type contact binaries
A&A
80
287 -
295

Citing 'WZ Cephei: A Dynamically Active W UMa-Type Binary Star
'

@article{ OJOOBS_2011_v28n3_163}
,title={WZ Cephei: A Dynamically Active W UMa-Type Binary Star}
,volume={3}
, url={http://dx.doi.org/10.5140/JASS.2011.28.3.163}, DOI={10.5140/JASS.2011.28.3.163}
, number= {3}
, journal={Journal of Astronomy and Space Sciences}
, publisher={The Korean Space Science Society}
, author={Jeong, Jang Hae
and
Kim, Chun-Hwey}
, year={2011}
, month={Sep}