MASS TRANSFER AND LIGHT TIME EFFECT STUDIES FOR AU SERPENTIS

Journal of The Korean Astronomical Society.
2015.
Feb,
48(1):
1-7

- Received : June 02, 2014
- Accepted : October 07, 2014
- Published : February 28, 2015

Download

PDF

e-PUB

PubReader

PPT

Export by style

Share

Article

Metrics

Cited by

TagCloud

The orbital period changes of the W UMa eclipsing binary AU Ser are studied using the (
O−C
) method. We conclude that the period variation is due to mass transfer from the primary star to the secondary one at a very low and decreasing rate
dP
/
dt
= −8.872×10
^{−8}
, superimposed on the sinusoidal variation due to a third body orbiting the binary with period 42.87±3.16 years, orbital eccentricity
e
= 0.52±0.12 and a longitude of periastron passage
ω
= 133°.7±15. On studying the magnetic activity, we have concluded that the Applegate mechanism failed to describe the cycling variation of the (
O−C
) diagram of AU Ser.
α
_{2000}
= 15
^{h}
56
^{m}
49
^{s}
,
δ
_{2000}
= +22°15′ 42′′ .3) was discovered by
Hoffmeister (1935)
and classified as a W-type binary system according to
Binnendijk (1972)
. It is a short period (
P
= 0
^{d}
.386) W UMa system with magnitude
V_{max}
. = 10
^{m}
.9 and spectral type G5V (
Kukarkin, 1970
).
Visual observations were made by
Soloviev (1936
,
1951)
, while photographic observations by
Huth (1964)
. The radial velocity curve was obtained using the cross correlation technique by
Hrivnak (1993)
. He found a mass ratio
q
=
m
_{2}
/
m
_{1}
= 0.71. He also obtained a projected total mass of AU Ser of (
m
_{1}
+
m
_{2}
)
sin
^{3}
i
= 1.51
M
_{⊙}
, a value similar to that obtained by
Pribulla et al. (2009)
from their spectroscopic observations.
The first photoelectric observations were made by
Binnendijk (1972)
and subsequently
Rucinski (1974)
,
Kennedy (1985)
,
Li, et al. (1992
,
1998)
,
Djurasevic (1993)
and
Gürol (2005)
. They obtained photoelectric light curves in different wavelengths. In addition,
Gürol (2005)
studied the period variation of the system suggesting its triplicity with a third body that orbiting the binary in about 94 years. However, many photoelectric and CCD minima times were observed and published during the last decade (39 minima) which motivated us to re-construct and analyze the
O − C
diagram in order to obtain more precise orbital parameters for the third body. The goal of the present study is to discuss the causes of variation in the orbital period; mainly to re-determine the third body orbital parameters of AU Ser.
O − C
residuals of
Table 2
have been calculated using the ephemeris given by
Kreiner et al. (2001)
:
where
E
is the number of integer cycles.
To predict timing of new minima, the last 147 minima times, which cover the interval from June 1976 till April 2012, have been linearly fitted (
Figure 1
) and used to obtain the new light elements:
The (O-C) diagram. The dashed line represents the linear fit to the last 143 minimum data.
with residual mean squares = 0.0002.
We have constructed the (
O − C
) diagram of AU Ser (
Figure 2
). It shows complicated period variations, which may be due to a combination of more than one of the following effects: mass transfer, light time effect, magnetic activity. Thus, we examine each reason affecting orbital period behavior.
The (O − C ) variation with quadratic and sinusoidal fit and their residuals of AU Ser. Filled circles represent primary minima while open circles are for secondary minima.
Q
= −4.69×10
^{−11}
and consequently we calculated the rate of change of the orbital period
dP
/
dt
= −8.872 × 10
^{−8}
day/year. The quadratic fit is represented as the dashed line on
Figure 2
.
In the simplest case of conservative mass transfer, if the more massive component loses mass, the orbital size will decrease and the period of the system must decrease too (
Pringle, 1985
). So, the obtained orbital period decrease may be interpreted in terms of mass transfer from the more massive star to the less one (
Kwee and van Woerden, 1958
). The rate of mass transfer in the conservative case could be estimated by using the formula derived by
Kreiner and Ziolkowski (1978)
:
where the quadratic term coefficient
Q
, and the period
P
are in days. Adopting the values of the absolute parameters of AU Ser from
Gürol (2005)
:
M
_{1}
= 0.895 M
_{⊙}
and
M
_{2}
= 0.635 M
_{⊙}
, the rate of mass transfer:
which is of the same order as the values given in the literature for contact binaries (see
Liu & Yang 2003
).
O − C
) diagram (
Figure 2
). The time of mid eclipse can be computed as follows:
where
e
_{3}
,
ω
_{3}
,
ν
,
a
_{12}
sin i and c are the eccentricity, longitude of the periastron, true anomaly of the binary orbit around the center of mass of the triple system, projected semi-major axis, and the speed of light, respectively.
Gürol was the first who considered the sine-like variation of the
O − C
curve of AU Ser. He calculated the light time effect (LITE) and obtained the orbital parameters of a third body of at least 0.53 solar mass orbiting the binary with an orbital period of 94.15 years. The recent photoelectric and CCD minima times (39 minima) together with the earlier minima could enable us to re-calculate the LITE and obtain a more real and precise solution.
Using the program prepared by
Zasche, et al. (2009)
, based on
Irwin’s method (1959)
, we have calculated the orbital parameters of the third body, and its LITE on the binary system. The program contains three modes 0, 1 and 2. These three modes correspond to computing the LITE due to the third body, the LITE together with the quadratic term and only the quadratic term, respectively.
Three different weights 1, 5 and 10 were applied to the data points for: visual (v), photographic (pg) and photoelectric (pe) or CCD minima times, respectively. Applying code 1, we obtain the orbital parameters of the third body as in
Table 1
, and as represented by the solid sinusoidal (red) line in
Figure 2
.
The light–time effect solution and the corresponding quadratic ephemeris of the binary system is also presented
On applying the well known mass function relation (cf.,
Albayrak et al., 1999
):
where
M
_{12}
and
M
_{3}
are the masses (in solar units) of the eclipsing pair and the third body,
A
is the amplitude in days, and
P
_{3}
is the period of the third body in years. The inclination
i
of the third body orbit was assumed to be equal to the inclination of the eclipsing binary orbit. The minimal mass
M
_{3}
_{min}
is then corresponding to
i
_{3}
= 90°.
From our new analysis and including all the recent minima times, we obtain a third body orbital period of about 43 years instead of the very large period (≃ 94 yr.) given by
Gürol (2005)
. Also, we obtain a significantly different value for the semi-amplitude. The comparison between the newly obtained set of parameters and those obtained earlier is given in
Table 1
.
O − C
) residual diagram for AU Ser contains a cycle of about 43 years. Assuming this long period
P
_{3}
to be the modulation period,
P_{mod}
, of the stellar magnetic activity of the convective secondary star, with semi amplitudes
O − C
= 0.0197 day, and accepting the parameters given by
Gürol (2005)
[
M
_{2}
= 0.635 M
_{⊙}
,
R
_{2}
= 0.94 R
_{⊙}
,
L
_{2}
= 3.8 L
_{⊙}
and the orbital semi-major axis
a
= 1.19 R
_{⊙}
] one can follow the Applegate procedure (see Applegate 1992).
The observed amplitude of the period modulation of the cycle, ∆
P
/
P
= 2
π
(
O − C
)/
P_{mod}
= 7.92 × 10
^{−6}
gives the variation of the orbital period ∆
P
= 0.264 second. The angular momentum transfer is ∆
J
= 3.69 × 10
^{+46}
g cm
^{2}
s
^{−1}
. If the mass of the shell is
M_{shell}
= 0.1 M
_{2}
, the moment of inertia of the shell is
I_{shell}
= 3.604 × 10
^{+53}
g cm
^{2}
, and the variable part of the differential rotation of the active star is ∆Ω/Ω = 0.00036. The energy budget needed to transfer the ∆
J
is ∆
E
= 7.564 × 10
^{+39}
ergs. The luminosity change is ∆
L_{RMS}
= 1.757 × 10
^{+31}
. This luminosity variation is ∆
L_{RMS}
/
L
= 0.0011 ≃ 0.0 of the luminosity of the active star. This value is inconsistent with the values suggested by Applegate (1992) model which should be around 10% to prove the presence of magnetic activity on similar chromospherically active stars.
O − C
) diagrams of eclipsing binaries and showed that migrating star spots can only introduce high frequency, low amplitude disturbances. Because of this
Gürol (2005)
deduced that the main causes of the (
O − C
) variations of AU Ser may be only mass transfer between the two stars and/or light time effect due to the presence of a third body. This shows that our result is in agreement to that obtained by Gürol.
The present analysis of the
O − C
diagram of AU Ser, suggests a decrease in the orbital period due to mass transfer from the more massive primary component to the less massive secondary one by a rate of about
dP
/
dt
= −8.87 × 10
^{−8}
day/year with a mass transfer rate of
dM
/
dt
= −1.67 × 10
^{−7}
M
_{⊙}
/year. This period decrease is superimposed on a sinusoidal variation, as seen in
Figure 2
, due to the presence of a third body orbiting the binary AU Ser in about 42.9 years with an orbital eccentricity
e
_{3}
= 0.52 and a longitude of periastron passage
ω
_{3}
= 133.7°.
Pribulla et al. (2009)
observed the system spectroscopically to study its radial velocity. They did not see evidence for a third component when using the broadening functions technique described by
Rusinski (2002)
. However, this result does not dismiss the third body hypothesis proposed in the present work and by
Gürol (2005)
. Hence, more precise photoelectric and CCD observations are still needed to verify the obtained results.

1. INTRODUCTION

AU Ser (
2. PERIOD VARIATION

To investigate the period changes of the W UMa system AU Ser, we have collected all the available times of minima since HJD 2428318.5 (≡ 30 May 1936) until 2456034.5 (≡ 17 April 2012) which cover about 75.9 years. The times of minima are listed in
Table 2
(Appendix). The successive columns of the table are: HJD, the number of integer cycles, type of the minimum, the O-C residuals and the references.
The
PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

- 2.1. Mass Transfer

A quadratic least square fit concerning the first three terms of equation (5) has been performed. We obtained the coefficient of the quadratic term
PPT Slide

Lager Image

PPT Slide

Lager Image

- 2.2. Light Time Effect

Observational detection of a periodic orbital period variation of a binary star system can be considered as a strong evidence of the existence of a third body around the binary system. This body causes the Light Time Effect (LITE) as displacement of the times of eclipse minimum light in a sinusoidal form with a period equal to the period of the third body, (
Woltjer 1922
;
Irwin 1959
).
At first, it may be noticed that both the primary and the secondary minima have the same trend on the (
PPT Slide

Lager Image

The light–time effect solution and the corresponding quadratic ephemeris of the binary system is also presented

PPT Slide

Lager Image

PPT Slide

Lager Image

- 2.3. Magnetic Activity and Star Spots

Kalužny (1986)
analyzed the light curves that were made by
Binnendijk (1972)
using the W-D code. He reported that the depth of the primary and the secondary minimum differs by about 0.2 mag and such difference is unusually large for a W UMa-type system. He also noticed the degree of asymmetry in maxima of the light curve and suggested that such pronounced asymmetry of the light curves may be due to a hot spot located at the neck between the stars.
Gürol (2005)
has studied the published light curves (LCs) and reported that all the LCs can be modelled by hot or cool spots located on the secondary components. He also studied the O’Connell effect and attributed such complex nature of the light curve to variable starspot activity with 32 to 35 years as a probable period for the system. For a detailed discussion for the light curve analysis concerning the starspot activity one may review
Gürol (2005)
.
Changes in the magnetic field distribution result in changes of angular momentum distribution. Gravitational quadrupole coupling produces changes in the internal structure of the active star which results in a period variation. The Applegate (1992) model involves variations of the subsurface magnetic field. Such subsurface magnetic field may be compared to solar activity cycles. The model can give a plausible explanation of the cyclic period variations of late type active stars.
For AU Ser, the star-spots are expected to be presented on the cooler member i.e., the secondary less massive star (Sp. Type G5V) was considered as the active component when applying the Applegate (1992) mechanism. For more details about the mechanism see
Applegate and Patterson (1987)
, Applegate (1992) and references therein. We applied the Applegate mechanism using of all data.
The present (
3. DISCUSSION AND CONCLUSIONS

In case of AU Ser especially for the observed light curves in 1992 and 1995,
Kalimeris et al. (2002)
studied the effects of star spots on the (
Acknowledgements

I acknowledge the variable star observations from the BBSAG International Database, IBVS, and the NASA Astrophysics Data System Bibliographic Services. Thanks go also to Dr. Petr Zasche for his program for the determination of the LITE due to the presence of the third body.

Albayrak B.
,
Özeren F.
,
Ekmekçi F.
,
Demircan O.
1999
Period Variation of Six RS CVn-Type Binaries with Possible Light-Time Effect
RMxAA
35
3 -

Applegate J. H.
,
Patterson J.
1987
Magnetic Activity, Tides, and Orbital Period Changes in Close Binaries
ApJ
322
L99 -
** DOI : 10.1086/185044**

Binnendijk L.
1972
The Light Variation of AU Serpentis
AJ
77
603 -
** DOI : 10.1086/111325**

Budding E.
,
Demircan O.
2007
Introduction to Astronomical Photometry, Chap. 8
2nd edn.
Cambridge University Press
New York

Djurasevic G.
1993
An Analysis of Close Binaries (CB) based on photometric measurements (7)
Ap&SS
206
207 -
** DOI : 10.1007/BF00658146**

Gürol B.
2005
Long Term Photometric and Period Study of AU Serpentis
NA
10
653 -
** DOI : 10.1016/j.newast.2005.04.004**

Hoffmeister C.
1935
162 Neue Veränderliche
AN
255
403 -

Hrivnak B. J.
,
Leung K. C.
,
Nha I. S.
1993
New Radial Velocity Curves of Six W UMa Binaries, New frontiers in binary star research
ASP Conference Series
38
269 -

Huth H.
1964
Photographische Beobachtungen von regelmässigen und Halbregelmässigen Veränderlichen auf Platten der Sonneberger Himmelsüberwachung
Mitt. Sonneberg
2
(5)
126 -

Irwin J. B.
1959
Standard light-Time curves
AJ
64
149 -
** DOI : 10.1086/107913**

Kalimeris A.
,
Rovithis- Livaniou H.
,
Rovithis P.
2002
Starspots and Photometric Noise on Observed minus Calculated (O-C) Diagrams
A&A
387
969 -

Kalužny J.
1986
Contact Binaries with Components in Poor Thermal Contact - Part Three - Au-Serpentis and Ft-Lupi
AcA
36
113 -

Kennedy H. D.
1985
The Period Behaviour of AU Serpintis
IBVS
2742

Kreiner J. M.
,
Ziolkowski J.
1978
Period Changes and Evolutionary Status of 18 Algol-Type Systems
AcA
28
497 -

Kreiner J. M.
,
Kim C.-H.
,
Nha I.-S.
2001
An Atlas of O-C Diagrams of Eclipsing Binary Star
Wydawnictwo

Kukarkin B. V.
1970
General Catalogue of Variable Stars
third ed.
Academy of Sciences
Moscow
373 -

Kwee K. K.
,
van Woerden H.
1958
Investigation of Variations in the Period of Sixteen Bright Short-Period Eclipsing Binary Stars
BAN
14
131 -

Liu Q.-Y.
,
Yang Y.-L.
2003
A Possible Explanation of the O’Connell Effect in Close Binary Stars
ChJAA
3
(2)
142 -

Li Z.-y.
,
Ding Y.-r.
,
Zhang Z.-s.
1998
The Variations in Light and Color Curves and Possible Short Period Oscillations of the Binary AU Serpentis
A&AS
131
115 -

Li Z.-y.
,
Zhan Z.-s.
,
Li Y.-l.
1992
The Variation of the Light Curves of AU Ser
IBVS
3802

Mayer P.
1990
Eclipsing Binaries with Light-Time Effect
BAICz
41
231 -

Pribulla T.
,
Rucinski S.
,
DeBond H.
,
DeRidder A.
,
Karmo T.
,
Thomson J.
,
Croll B.
,
Ogloza W.
,
Pilecki B.
,
Siwak M.
2009
Radial Velocity Studies of Close Binary Stars. XIV
AJ
137
3646 -
** DOI : 10.1088/0004-6256/137/3/3646**

Pringle J. E.
,
Pringle J. E.
,
Wade R. A.
1985
Interacting Binary Stars
Cambridge University Press
Cambridge
1 -

Rucinski S. M.
1974
Binaries. II. A- and W-type Systems. The W UMa-type Systems as Contact
AcA
24
119 -

Rucinski S. M.
2002
Radial Velocity Studies of Close Binary Stars. VII. Methods and Uncertainties
AJ
124
1746 -
** DOI : 10.1086/342342**

Soloviev A. V.
1936
Observations of Variables of RR Lyrae Type, Tadjik Obs. Circ., No. 21

Soloviev A. V.
1951
W Ursae Majoris
PZ
8
64 -

Woltjer J.
1922
On a Special Case of Orbit Determination in the Theory of Eclipsing Variables
Bull. Astron. Inst. Netherlands
1
93 -

Zasche P.
,
Liakos A.
,
Niarchos P.
,
Wolf M.
,
Manimanis V.
,
Gazeas K.
2009
Period Changes in Six Contact Binaries: WZ And, V803 Aql, DF Hya, PY Lyr, FZ Ori, and AH Tau
New A.
14
121 -
** DOI : 10.1016/j.newast.2008.06.002**

Citing 'MASS TRANSFER AND LIGHT TIME EFFECT STUDIES FOR AU SERPENTIS
'

@article{ CMHHBA_2015_v48n1_1}
,title={MASS TRANSFER AND LIGHT TIME EFFECT STUDIES FOR AU SERPENTIS}
,volume={1}
, url={http://dx.doi.org/10.5303/JKAS.2015.48.1.1}, DOI={10.5303/JKAS.2015.48.1.1}
, number= {1}
, journal={Journal of The Korean Astronomical Society}
, publisher={The Korean Astronomical Society}
, author={AMin, S. M.}
, year={2015}
, month={Feb}