New observations for the times of minimum lights of a well-known apsidal motion star CW Cephei were made using a 0.6 m wide field telescope at Jincheon station of Chungbuk National University Observatory, Korea during the 2015 observational season. We determined new times of minimum lights from these observations and analyzed
O
-
C
diagrams together with collected times of minima to study both the apsidal motion and the Light Time Effect (LTE) suggested in the system. The new periods of the apsidal motion and the LTE were calculated as 46.6 and 39.3 years, respectively, which were similar but improved accuracy than earlier ones investigated by
Han et al. (2002)
,
Erdem et al. (2004)
and
Wolf et al. (2006)
.
1. INTRODUCTION
CW Cep is a well-known apsidal motion binary system observed by many investigators
(Abrami & Cester 1960
;
Nha 1975
;
Han et al. 2002
;
Wolf et al. 2006)
since its discovery by
Petrie (1947)
as a double lined spectroscopic eclipsing binary. The spectral types of two components of CW Cep were determined as B0.5 + B0.5, and masses were calculated as 11.7 M
⊙
and 11.0 M
⊙
by
Popper (1974)
. The light curves were analyzed by many authors using well observed full-light curves with different models for the parameters of the binary system, related with its apsidal motion
(Cester et al. 1978)
with the WINK model,
Clausen & Giménez (1991)
by the EBOP model, and
Terrell (1991)
, and
Han et al. (2002)
by the WD model. While
Popper & Hill (1991)
discussed the radial velocity curve based on the spectroscopic observations by
Popper (1974)
, high resolution IUE observations for the radial velocity curve were made by
Stickland et al. (1992)
. Most of these studies suggested that CW Cep is a detached main sequence binary system, composed of two stars of slightly different mass and radius.
With regard to the apsidal motion study of the close eclipsing binary systems, CW Cep showed one of the shortest apsidal period, among similar binary systems, which leads particular interests for this binary system. Based on intensive photoelectric observations with
UBV
filters,
Nha (1975)
calculated the apsidal motion period of CW Cep as 39 years, while
Han (1984)
derived the apsidal motion period as 43.1 years with more times of minimum lights based on the
UBV
photoelectric observations. After their studies, more times of minimum lights were available by many investigators recently, according to the new observations. Therefore, the apsidal motion period of CW Cep was improved in time order as 45.5 years by
Gimenez et al. (1987)
, 45.6 years by
Clausen & Giménez (1991)
, and 45.7 years by
Han et al. (2002)
. The last authors, for the first time, noticed a new type of complication confusing the apsidal motion behavior, which shows about 34 years periodic variation with an amplitude of 0.003 days.
Erdem et al. (2004)
calculated the apsidal motion period of CW Cep as 46.1 years through the variation of the longitude of periastron (
ω
) of the system's orbit, with new observations using
BVR
filters. More recently, with new observations by
Wolf et al. (2006)
, they suggested the apsidal motion period of CW Cep as 46 years, which is similar to those of other previous authors
(Clausen & Giménez 1991
;
Han et al. 2002)
. They confirmed an extra periodic change suggested by
Han et al. (2002)
and proposed a Light Time Effect (LTE) with the period of 38.5 years and the small semi-amplitude of 0.002 days caused by an unseen third body.
It is nearly 10 years since the observations by
Wolf et al. (2006)
, and this indicates the necessity of new observations for this binary system to examine more accurate apsidal motion period. In line with this, we made new observations for CW Cep in the 2015 season and obtained three new times of minimum lights to confirm the apsidal motion period of this system, together with available minimum times from the published literatures.
2. NEW OBSERVATIONS
New CCD photometric observations presented in this work were made in the 2015 observational season at the Jincheon station of Chungbuk National University Observatory (CBNUOJ), Korea. The altitude of the Jincheon station is 87 m sea level and the 60-cm reflecting telescope is installed by Korea Astronomy and Space Science Institute (KASI), and operated by CBNUOJ. An electronically cooled SBIG STX- 16803 4K CCD camera is equipped with a Johnson standard
BVRI
filter set, and is attached at the reflector f/2.9 prime focus of the 60 cm wide field telescope. Since our observations are focused on the acquisition of accurate eclipse timings, only B filter is used for fast time-resolved photometry for the first two nights. The observations of the last night were made with
BVR
filters to derive the weighted average of each filters for higher accuracy. The field of view of the CCD camera is 72’ × 72’ with a pixel scale of 1.05 arcsec/pixel, which is wide enough to cover the many field stars around CW Cep including the comparison star GSC 04282-00330, and differential photometry was applied for this observation. All of the observed CCD frames were exposed 15 ~ 80 seconds depending on the observational circumstances, and corrected with bias, dark and flat-field images for the differential photometry. The IRAF/DIPHO software package was applied to process the observed CCD frames. The details of the data reduction were extensively described by
Kim et al. (2014)
based on the same observational circumstances. We obtained three times of minimum lights, two for primary and one secondary minima, and determined the times of minimum lights by conventional
Kwee & van Woerden (1956)
method as listed in
Table 1
together with other times of minima determined by previous investigators.
The times of minimum lights of CW Cep
1) PE: Photoelectric photometry, CCD photometry 2) CBNUOJ : Chungbuk National University Observatory at Jincheon
3. MINIMUM LIGHTS ANDO-CDIAGRAM
With our new times of minimum lights of CW Cep, we collected all published times of minimum lights from the literature to analyze the apsidal motion of this binary system. Among the available times of minimum lights, the visual and photographic data observed at the early stage of apsidal motion were excluded in this study because those data could increase uncertainty of period variations. Only photoelectric and CCD times of minimum lights were used in our analysis and listed in
Table 1
, including our timings.
These times of minima were simultaneously analyzed with the apsidal motion plus the light-time effect ephemeris represented by:
where
Ps
is the sidereal period,
τap
and
τ
3
denote the apsidal motion and the LTE terms, respectively.
τap
is a function of four parameters (anomalous period
Pa
, eccentricity
e
, longitude of the periastron
ω
0
, angular velocity of the line of apsides
) The formalism for
τ
3
was given up to sixth power of the eccentricity
Giménez & Bastero (1995)
.
Irwin (1952)
gave an exact form for
τ
3
with five unknown parameters which are the semi-amplitude
K
, the eccentricity
e
3
, the period
P
3
, the periastron passage time
T
3
, the longitude of periastron
ω
3
for the light orbit of mass center of the eclipsing pair around the mass center of the triple system.
All the timings were fitted to Eq. (1) by using the code (available at
http://sirrah.troja.mff.cuni.cz/~zasche/Programs.html
) given by
Zasche et al. (2009)
and
Zasche & Wolf (2013)
. We made small corrections of the code by considering the inclination of the eclipsing pair in it involving up to sixth power of the eccentricity. With the initial values of the apsidal motion and LTE parameter given by
Wolf et al. (2006)
, we improved iteratively the apsidal motion and LTE parameters. the calculation we fixed the inclination of the eclipsing pair as 82.5° derived by
Clausen & Giménez (1991)
and
Han et al. (2002)
. The final solution was listed in the sixth column of
Table 2
, together with those et al. (2002),
Wolf et al. (2006)
and
Bulut et al. (2011)
for comparison. The resultant O-C diagram was drawn with the linear term of our result in the top of
Fig. 1
where the solid continuous curves denote the combination of
τap
and
τ
3
and the dashed curve represents only
τ
3
term. In the middle the residuals from the linear plus the apsidal motion ephemeris were plotted with the theoretical LTE
τ
3
term. The residuals from the full contribution of Eq. (1) were plotted in the bottom panel. The standard deviation of the residuals with a relatively large scatter was calculated to be ±0.0025 days which is slightly larger than the small semi-amplitude of 0.0021 days for the LTE orbit as listed in
Table 2
.
Eclipse timing variation diagram with theoretical lines calculated with the parameters in Table 2. Open and cross symbols refer to primary and secondary minimum timings, respectively, where circle and square represent Photoelectric (PE) and CCD timings, respectively. The black and grey solid lines on the top panel represent the theoretical apsidal motion plus LTE terms for primary and secondary eclipses, respectively. The dash lines on the top and middle panels also represent the theoretical LTE (τ3). The bottom panel indicates the residual from Eq (1).
Parameters of the apsidal motion and LTE for CW Cep
Parameters of the apsidal motion and LTE for CW Cep
Because the Zasche code we used was based on a differential least-square method
Zasche et al. (2009)
, our solution of Eq. (1) may be a solution with a local minimum rather than a global solution. It would be very necessary to check whether our solution is globally valid for all parameter space. In line with this validity, we subsequently attempted a Monte Carlo simulation to obtain optimized values of the parameters both of the apsidal motion and the LTE orbit as performed by
Lee et al. (2015)
. After 800 simulations all the apsidal motion parameters showed a quick convergence to nearly the same values listed in
Table 2
, while the LTE parameters show relatively wide variation as shown in the histograms of
Fig. 2
. Particularly, the distribution of
e
3
showed remarkably wide variation, indicating that the proposed LTE orbit is not well defined with reliable accuracy by present observational data.
Histograms for distributions of the third body orbit parameters. The black lines represent the calculated elements of the third body orbit.
4. DISCUSSION OF APSIDAL MOTION PERIOD
In this study, the new times of minimum lights were presented by the observations of CBNUOJ, together with the collected times of minimum lights as shown in
Table 1
. The
O-C
diagram of a differential least-square fit to the minimum lights were shown in
Fig. 1
. We calculated the new apsidal motion period as 46.27 ± 0.04 years. The apsidal parameters are listed in
Table 2
, which shows similar apsidal motion period with those of earlier investigators
(Han at al 2002
;
Wolf et al. 2006)
, however with much improved accuracy. The existence of a third body in the CW Cep system was also suggested by many previous investigators (Han et al. 2000;
Wolf et al. 2006
;
Bulut et al. 2011)
. Based on the analysis of currently available minimum lights in
Table 1
, we calculated the period of the LTE caused by the third body as 39.3 ± 4.4 years as shown dash line in
Fig. 1
. However, the LTE caused by the third body may not be convinced due to a wide variation of orbital eccentricity (
e
3
) in
Fig. 2
and a small semi-amplitude compatible to a large scatter band of the residuals in
Fig. 1
.
The masses and radii of CW Cep by the photometric and spectroscopic parameters suggested by
Clausen & Giménez (1991)
are
M
1
=11.82 ± 0.14
M
⊙
,
M
2
=11.09 ± 0.14
M
⊙
,
R
1
=5.48 ± 0.12
R
⊙
, and
R
2
=4.99 ± 0.12
R
⊙
. With these parameters, we calculated the internal structure constant of the system, as log
k
2
,obs
=-2.106 that is slightly small value rather than the theoretical internal structure constant -2.08 suggested by
Wolf et al. (2006)
. If the third body component exists around CW Cep, the minimum mass of the third can be calculated by the mass function (
f(M
3
)) of 4×10
-5
. The results show that the third body component should be 0.279 M
⊙
when
i
3
is 90°, and in the case of
i
3
is 30°, the mass is 0.562 M
⊙
. These are approximately 2% range of mass of CW Cep binary system, which is too small compared to CW Cep and the masses are close to red dwarfs range. Therefore detecting the third body by photometry and spectroscopy would be very difficult. Apsidal motion period is very important because it can provide the important constant for the stellar internal structure. In this regards, many investigators have paid attention to the apsidal motion of CW Cep, and more accurate observations will be necessary to analyze the apsidal motion of the binary system and periodic variations by the third body.
Acknowledgements
We thank the staff of CBNUOJ for their assistance with our observations. WH acknowledges partial support by grant number NRF-2014M1A3A3A02034746 from the National Research Foundation of Korea, and C-H K was partially supported by the National Research Foundation of Korea (NRF) grant funded by the Korean Government (NRF- 2015R1D1A1A01058924).
Abrami A
,
Cester B
(1960)
Studio photoelectric spectrophotometric CW Cephei binary system, in Osservatorio astronomico di Trieste, vol. 300
Pavia
Trieste
24 -
Agerer F
,
Hübscher J
(1996)
Photoelectric Minima and Maxima of Selected Eclipsing and Pulsating Variables
Inform. Bull. Var. Stars
4382
1 -
4
Agerer F
,
Hübscher J
(1997)
Photoelectric Minima of Selected Eclipsing Binaries and Maxima of Pulsating Stars
Inform. Bull. Var. Stars
4472
1 -
5
Agerer F
,
Hübscher J
(2002)
Photoelectric Minima of Selected Eclipsing Binaries and Maxima of Pulsating Stars
Inform. Bull. Var. Stars
5296
1 -
16
Agerer F
,
Hübscher J
(2003)
Photoelectric Minima of Selected Eclipsing Binaries
Inform. Bull. Var. Stars
5484
1 -
13
Bíró IB
,
Borkovits T
(2000)
CCD Times of Minima of Eclipsing Binary Systems
Inform. Bull. Var. Stars
4967
1 -
4
Brancewicz H
,
Kreiner JM
(1976)
Photoelectric Minima of Eclipsing Binaries
Inform. Bull. Var. Stars
1119
1 -
Brát L
,
Šmelcer L
,
Kučáková H
,
Ehrenberger R
,
Kocián R
(2008)
B.R.N.O Times of minima
Open Eur. J. Var. Stars
94
1 -
28
Brelstaff T
(1994)
Photoelectric Minima of Eclipsing Binaries 1992- 1993
British Astron. Assoc. Var. Star Sect. Circ.
81
4 -
5
Brelstaff T
(1995)
Photoelectric Minima of Eclipsing Binaries 1994
British Astron. Assoc. Var. Star Sect. Circ.
84
(2)
2 -
Bulut A
,
Çiçek C
,
Bulut İ
,
Erdem A
(2011)
LIGHT-TIME EFFECT and APSIDAL MOTION of TWO ECLIPTING BINARIES: AO Vel and CW Cep
Balkan. Phys. Lett.
19
283 -
288
Caton DB
,
Burns WC
(1993)
Times of Minimum Light for 35 Eclipses of 21 Apsidal Motion Binaries
IBVS
3900
1 -
4
Cester B
,
Fedel B
,
Giuicin G
,
Mardirossian F
,
Mezzetti M
(1978)
Revised photometric elements of 14 detached system
Astron. Astrophys. Suppl.
33
91 -
102
Clausen JV
,
Giménez A
(1991)
Absolute dimensions of eclipsing binaries XVIII. The Cepheus OB 3 member CW Cephei
Astron. Astrophys.
241
98 -
106
Diethelm R
(1988)
122thList of Minima of Eclipsing Binaries, Bull. Bedeckungsveränderliche-Beobachter Schweiz.
Astron. Ges.
89
1 -
12
Diethelm R
(1992)
135thList of Minima of Eclipsing Binaries, Bull. Bedeckungsveränderliche-Beobachter Schweiz.
Astron. Ges.
102
1 -
10
Diethelm R
(2012)
Timings of Minima of Eclipsing Binaries
Inform. Bull. Var. Stars
6011
1 -
13
Doĝru SS
,
Erdem A
,
Dönmez A
,
Bulut A
,
Akin T
(2009)
New Times of Minima of Some Eclipsing Binary Stars
Inform. Bull. Var. Stars
5893
1 -
3
Erdem A
,
Soydugan E
,
Soydugan F
,
Özdemir S
,
Çiçek C
(2004)
A new photometric study of CW Cephei and its apsidal motion
Astron. Note
http://dx.doi.org/10.1002/asna.200310204
325
336 -
342
Giménez A
,
Bastero M
(1995)
A revision of the ephemeris-curve equations for eclipsing binaries with apsidal motion
Astrophys. Space Sci.
226
99 -
107
DOI : 10.1007/BF00626903
Giménez A
,
Kim CH
,
Nha IS
(1987)
Apsidal motion in the early-type eclipsing binaries CW Cephei, Y Cygni and AG Persei
Mon. Not. Roy. Astron. Soc.
224
543 -
555
DOI : 10.1093/mnras/224.3.543
Han W
(1984)
Photoelectric Observations of the Close Eclipsing Binary System CW Cephei
J. Astron. Space Sci.
1
Han W
,
Kim CH
,
Lee WB
,
Koch RH
(2002)
Photometric studies of CW Cephei
Astron. J.
123
2724 -
2732
DOI : 10.1086/340072
Hübscher J
(2005a)
Photoelectric Minima of Selected Eclipsing Binaries and Maxima of Pulsating Stars
Inform. Bull. Var. Stars
5643
1 -
20
Hübscher J
(2011a)
BAV Result of Observations – Photoelectric Minima of Selected Eclipsing Binaries and Maxima of Pulsating Stars
Inform. Bull. Var. Stars
5984
1 -
16
Hübscher J
(2014)
BAV Result of Observations – Photoelectric Minima of Selected Eclipsing Binaries and Maxima of Pulsating Stars
Inform. Bull. Var. Stars
6118
1 -
17
Hübscher J
,
Walter F
(2007)
Photoelectric Minima of Selected Eclipsing Binaries and Maxima of Pulsating Stars
Inform. Bull. Var. Stars
5761
1 -
12
Hübscher J
,
Monninger G
(2011b)
BAV Result of Observations – Photoelectric Minima of Selected Eclipsing Binaries and Maxima of Pulsating Stars
Inform. Bull. Var. Stars
5959
1 -
16
Hübscher J
,
Lehmann PB
(2012a)
BAV Result of Observations – Photoelectric Minima of Selected Eclipsing Binaries and Maxima of Pulsating Stars
Inform. Bull. Var. Stars
6026
1 -
24
Hübscher J
,
Lehmann PB
(2013)
BAV Result of Observations – Photoelectric Minima of Selected Eclipsing Binaries and Maxima of Pulsating Stars
Inform. Bull. Var. Stars
3070
1 -
15
Hübscher J
,
Agerer F
,
Wunder E
(1992)
Beobachtungsergebnisse der Berliner Arbeitsgemeinschaft für Veränderliche Sterne e.V. (BAV) (BAV-Mitteilungen Nr. 60)
Self-published
1 -
16
Hübscher J
,
Agerer F
,
Frank P
,
Wunder E
(1994)
Beobachtungsergebnisse der Berliner Arbeitsgemeinschaft für Veränderliche Sterne e.V. (BAV) (BAV-Mitteilungen Nr. 68)
Self-published
1 -
24
Hübscher J
,
Paschke A
,
Walter F
(2005b)
Photoelectric Minima of Selected Eclipsing Binaries and Maxima of Pulsating Stars
Inform. Bull. Var. Stars
5657
1 -
24
Hübscher J
,
Steinbach H-M
,
Walter F
(2008)
BAV Result of Observations – Photoelectric Minima of Selected Eclipsing Binaries and Maxima of Pulsating Stars
Inform. Bull. Var. Stars
5830
1 -
8
Hübscher J
,
Steinbach H-M
,
Walter F
(2009)
BAV Result of Observations – Photoelectric Minima of Selected Eclipsing Binaries and Maxima of Pulsating Stars
Inform. Bull. Var. Stars
5889
1 -
13
Hübscher J
,
Lehmann PB
,
Walter F
(2012b)
BAV Result of Observations – Photoelectric Minima of Selected Eclipsing Binaries and Maxima of Pulsating Stars
Inform. Bull. Var. Stars
6010
1 -
22
Irwin JB
(1952)
The determination of a light-time orbit
Astrophys. J.
116
211 -
217
DOI : 10.1086/145604
Jordi C
,
Ribas I
,
Gracía JM
(1996)
Times of Minima of Five Eclipsing Binaries
Inform. Bull. Var. Stars
4300
1 -
3
Kim CH
,
Song MH
,
Yoon JN
,
Han W
,
Jeong MJ
(2014)
BD Andromedae: A New Short-period RS CVn Eclipsing Binary Star with a Distant Tertiary Body in a Highly Eccentric Orbit
Astrophys. J.
788
134 -
155
DOI : 10.1088/0004-637X/788/2/134
Kwee KK
(1956)
van Woerden H, A method for computing accurately the epoch of minimum of an eclipsing variable
Bull. Astron. Inst. Neth.
12
327 -
330
Marino G
,
Arena C
,
Bellia I
,
Benintende G
,
Cremaschini C
(2010)
CCD Minim.
Eclips. Bin. Stars
5917
1 -
4
Nha IS
(1975)
CW Cephei: An Important Close Binary Member of the III Cephei Association
Astron. J.
80
232 -
238
DOI : 10.1086/111736
Ogloza W
(1995)
Photoelectric Minima of Eclipsing Binaries
Inform. Bull. Var. Stars
4236
1 -
2
Petrie, RM
(1947)
The spectrographic orbits and dimensions of H.D. 218066
Publ. Dom. Astrophys. Obs. Ott.
7
305 -
309
Popper DM
(1974)
Radial velocity measurements and spectrographic orbits
Astron. J.
79
1307 -
1313
DOI : 10.1086/111679
Popper DM
,
Hill G
(1991)
Rediscussion of eclipsing binaries. XVII – Spectroscopic orbits of OB systems with a cross-correlation procedure
Astron. J.
101
600 -
615
DOI : 10.1086/115709
Söderhjelm S
(1976)
A Re-discussion of CW Cep
Astron. Astrophys. Suppl.
25
151 -
158
Soydugan E
,
Demircan O
,
Keskin V
,
Erdem A
,
Çiçek C
(2001)
BVR Photometry of CW Cephei
Inform. Bull. Var. Stars
5154
1 -
2
Stickland DJ
,
Koch RH
,
Pfeiffer RJ
(1992)
Spectroscopic binary orbits from ultraviolet radial velocities. X - CW Cephei (HD 218066)
Obs.
112
277 -
281
Terrell D
(1991)
A re-discussion of four early-type eclipsing binary systems
Mon. Not. Roy. Astron. Soc.
250
209 -
214
DOI : 10.1093/mnras/250.1.209
Wolf M
,
Kucáková H
,
Kolasa M
,
Stastný P
,
Bozkurt Z
(2006)
Apsidal Motion in Eccentric Eclipsing Binaries: CW Cephei, V478 Cygni, AG Persei, and IQ Persei
Astron. Astrophys.
456
1077 -
1083
DOI : 10.1051/0004-6361:20065327
Zasche P
,
Wolf M
2013)
Apsidal motion and absolute parameters for five LMC eccentric eclipsing binaries
Astron. Astronphys.
558
A51 -
DOI : 10.1051/0004-6361/201322054
Zasche P
,
Liakos A
,
Niarchos P
,
Wolf M
,
Manimanis V
(2009)
Period changes in six contact binaries: WZ And, V803 Aql, DF Hya, PY Lyr, FZ Ori, and AH Tau
New Astron.
14
121 -
128
DOI : 10.1016/j.newast.2008.06.002
Zasche P
,
Uhlař R
,
Kučáková H
,
Svoboda P
,
Mašek M
(2014)
Collection of Minima of Eclipsing Binaries
Inform. Bull. Var. Stars
6114
1 -
19