A&A 376, 158-164 (2001)
DOI: 10.1051/0004-6361:20010914
B. Albayrak1 - O. Demircan2 - G. Djurasevic3 - S. Erkapic3 - H. Ak 1
1 - Ankara University, Science Faculty, Dept. of Astronomy
and Space Sciences, 06100 Tandogan, Ankara, Turkey
2 -
Çanakkale Onsekiz Mart University, Faculty of
Arts and Sciences, Dept. of Physics, 17100 Çanakkale, Turkey
3 -
Astronomical Observatory, Volgina 7, 11050 Belgrade, Yugoslavia
Received 17 May 2001 / Accepted 22 June 2001
Abstract
New BVR light curves and times of minimum light for the
short period RS CVn system SV Cam were analysed to derive the physical
parameters of the system and the parameters of the third body orbit. The
light curves obtained at the TÜBITAK National Observatory
during two nights in 2000 show considerable asymmetry and night-to-night
variations. The analysis of the light curves is made using Djurasevic's inverse
problem method. The Roche model with spotted areas on the hotter primary
component yields a good fit to observations. The extensive series of
published photoelectric minima times indicate that the eclipsing pair orbits
around the common mass center of the triple system with a period of 41.32
yr.
Key words: stars: binaries: eclipsing - stars: individual: SV Cam
The short-period (
)
RS CVn type totally eclipsing binary
system SV Cam (
,
HD 44982, BD+82
174,
HIC 32015) has been extensively observed since its discovery by Guthnick
(1929). Light variations and asymmetry in the overall shape of the
light curves, as much as 0.1 mag in time scales as short as one month, had
been noticed by all observers (e.g. van Woerden 1957; Hilditch et al.
1979; Patkos 1982a,b; Cellino et al.
1985; Milano et al. 1986). Hilditch et al.
reported that the cool secondary component is the source of distortion in
the light curves. On the contrary, Cellino et al. claimed that the active spotted
component should be the hotter primary. This hypothesis was supported by the
brightening they observed at phase 0.5, at the same time as the increased
reddening at the primary eclipse due to the passage of the spotted area through
phase 0.0. On the other hand, Özeren et al. (2001) found
an excess emission in both H
and H
lines connected
with the secondary, cooler component of the system.
The light curve solutions indicate that the primary component is very close to its Roche limit and provides at least 90% of the total light of the system. It is therefore not possible for the fainter cool secondary (which contributes less than 10% to the system light) to be the source of large distortion in the light curves.
The spectral type of the component stars were estimated to be G2-3 V
and K4 V by Hilditch et al. (1979), giving a mass
ratio of about 0.71. The later G8 V classification of the primary
component by Popper (1996) is based on the Na D lines.
From a radial velocity study, Pojmanski (1998) estimated a new
mass ratio of
and found an earlier spectral type for
the primary. The indices (h, c) refer to the hotter (more-massive) primary
and cooler (less-massive) secondary component respectively. From
the stellar masses, he derived spectral types F5 V and KO V
for the components of SV Cam system.
The light-time effect possibility of orbital motion about a third body was first suggested by Friebos-Conde & Herczeg (1973) to explain the period variation of SV Cam. The existence of the third body in the system was later confirmed by Hilditch et al. (1979), Cellino et al. (1985), Sarma et al. (1985), and Albayrak et al. (1999) with extended data. However, the parameters of the third body orbit were found to be very different in each study.
This paper presents three new light curves in B, V, and R,
obtained in 2000. These were analysed with a fixed revised mass
ratio of
.
We used two different hypotheses
concerning the temperature of the hotter primary component:
K
(Hyp. I.) and
K (Hip. II.) corresponding to
F5 V and F8 V stars, and assumed
cool spots on the hotter component. The period variation of the system was
also studied again by using the newly extended data in the hope of finding
more reliable parameters of the third body orbit.
New observations of SV Cam in BVR were obtained at the TÜBITAK National Observatory (Turkey) on the nights of 22 January
and 2 February 2000, by using a SSP-5A photometer attached to a 0.4 m
Cassegrain telescope. The same comparison star BD+82
0176 was chosen
as Milano et al. (1986).
A total of 242 observations were secured in each filter. The nightly extinction
coefficients for each colour were determined by using the observations of the
comparison star. The instrumental differential B, V, and R
magnitudes in the sense variable minus comparison corrected for the
atmospheric extinction, are given in Table 1 (accessible in electronic form).
The differential B, V, and R light and B-V and
V-R colour curves are shown on the upper right-hand of Fig.
2. The photometric phases of the light and colour curves
are calculated with the following ephemeris given by Pojmanski
(1998):
![]() |
(1) |
During the observations three primary and two secondary minima were obtained. Their timings were calculated by using the well-known method of Kwee & van Woerden (1956). The photoelectric minima published before mid-1997 were collected by Albayrak et al. (1999). We added minima obtained later by Agerer & Hübscher (1998), Pribulla et al. (1999, 2001) and by the present study. Altogether 152 primary and 18 secondary photoelectric minima are given in Table 2. Since the secondary minima are very shallow, the precise determination of the time of each minimum is strongly affected by the light-curve asymmetry. Thus, four secondary minima times were omitted as their scatter suggests they are not as well determined as the rest of the values. They are marked in italics in Table 2. The linear light elements given by Eq. (1) were used in estimating C values. The (O-C) diagrame in Fig. 1 can be represented by a continuous oscillatory variation covering one minimum and almost two maxima by the present observational data. This oscillation is produced by the light time effect due to a third body in the system.
To derive the light-time orbit and the parameters of the invisible component,
the analysis was performed by using Mayer's (1990) equations. We
used the scientific graphic, curve fitting and statistic program PRISIM
(Cass 2000) to obtain the parameters, listed in Table 3 with their
standard deviations. Table 3 also compares the results of this analysis with
those of different authors for SV Cam. The mass of the third body may be
evaluated from the mass function as the (O-C) analysis allows
us to obtain f(M3). But, it depends on
the inclination of the long-period orbit. The inclination i of the third
body orbit was assumed to be equal to the inclination of the eclipsing
binary orbit found as
in this study (see Table 4). Since masses of the
components of the eclipsing pair were derived by us are 1.42 and 0.80 M
(see Table 4), the mass of the third component is
M
.
If we take its mass and the period as 0.26
M
and 41.32 yr, respectively, Kepler's third law gives the
semi-major axis of the orbit to be 16.18 AU. Assuming the Hipparcos distance
of SV Cam d = 84.96 pc, we get the angular separation of the third
body from the eclipsing pair as
.
HJD* | Ref. | HJD* | Ref. | HJD* | Ref. | HJD* | Ref. | HJD* | Ref. | HJD* | Ref. |
primary | |||||||||||
32253.8190 |
1 | 33791.6582 | 1 | 41835.4800 | 9 | 42634.3500 | 9 | 43880.3930 | 9 | 45043.4053 | 18 |
32265.6819 | 1 | 33803.5200 | 1 | 41905.4670 | 9 | 42771.3500 | 11 | 43892.2540 | 15 | 45627.5824 | 19 |
32268.6483 | 1 | 33844.4420 | 1 | 41930.3770 | 9 | 42777.2786 | 11 | 43927.2450 | 9 | 45658.4215 | 19 |
32281.6963 | 1 | 33895.4451 | 1 | 41931.5635 | 9 | 42777.2792 | 11 | 43928.4300 | 9 | 45696.3787 | 19 |
32287.6263 | 1 | 33911.4575 | 1 | 41933.3437 | 9 | 42829.4700 | 9 | 44049.4184 | 9 | 45741.4540 | 19 |
32878.9214 | 1 | 33921.5415 | 1 | 41934.5278 | 9 | 42830.6550 | 9 | 44081.4440 | 9 | 46334.5271 | 20 |
32883.6648 | 1 | 33923.3209 | 1 | 41959.4370 | 9 | 42836.5850 | 9 | 44103.3850 | 9 | 46827.3718 | 21 |
32911.5381 | 1 | 33927.4731 | 1 | 41960.6230 | 9 | 42836.5905 | 12 | 44167.4387 | 16 | 50096.9835 | 22 |
32949.4949 | 1 | 33928.6591 | 1 | 41962.4024 | 9 | 42852.6010 | 12 | 44167.4388 | 16 | 50259.4818 | 23 |
33179.6094 | 1 | 33931.6243 | 1 | 41978.4160 | 9 | 42855.5658 | 12 | 44190.5665 | 16 | 50268.3809 | 23 |
33180.7937 | 1 | 33943.4847 | 1 | 41981.3800 | 9 | 42871.5780 | 9 | 44190.5676 | 16 | 50673.4498 | 24 |
33183.7583 | 1 | 38652.4790 | 2 | 41982.5665 | 9 | 42947.4910 | 9 | 44278.3442 | 16 | 50839.5105 | 25 |
33314.8284 | 1 | 38667.3120 | 2 | 41984.3449 | 9 | 43061.3603 | 9 | 44278.3443 | 16 | 50849.5929 | 25 |
33741.8398 | 1 | 38671.4600 | 2 | 42019.3347 | 9 | 43077.3740 | 9 | 44285.4614 | 9 | 51160.3641 | 25 |
33741.8405 | 1 | 39681.4567 | 3 | 42106.5154 | 9 | 43078.5608 | 9 | 44291.3915 | 17 | 51166.2944 | 25 |
33761.4114 | 1 | 39776.3460 | 3 | 42304.6040 | 9 | 43115.3310 | 6 | 44345.3617 | 9 | 51179.3434 | 25 |
33762.5979 | 1 | 39945.3718 | 4 | 42366.2810 | 10 | 43135.4945 | 9 | 44371.4559 | 9 | 51435.5490 | 26 |
33768.5288 | 1 | 39977.3980 | 4 | 42405.4243 | 9 | 43138.4610 | 13 | 44454.4869 | 9 | 51566.6192 | 27 |
33769.7144 | 1 | 40092.4561 | 4 | 42460.5805 | 9 | 43192.4293 | 9 | 44477.6150 | 9 | 51572.5499 | 27 |
33769.7159 | 1 | 40127.4466 | 4 | 42465.3249 | 9 | 43198.3599 | 9 | 44582.5892 | 9 | 51577.2936 | 27 |
33769.7174 | 1 | 40528.3615 | 5 | 42517.5142 | 11 | 43218.5236 | 9 | 44642.4879 | 16 | 51878.5755 | 26 |
33775.6465 | 1 | 40593.5983 | 6 | 42523.4453 | 9 | 43263.5948 | 14 | 44642.4888 | 16 | 51921.2793 | 26 |
33775.6470 | 1 | 40857.5139 | 7 | 42545.3880 | 9 | 43288.5050 | 9 | 44642.4893 | 16 | ||
33777.4241 | 1 | 41681.2879 | 8 | 42545.3895 | 11 | 43393.4808 | 9 | 44731.4486 | 16 | ||
33784.5405 | 1 | 41695.5194 | 9 | 42545.3909 | 11 | 43849.5518 | 9 | 44731.4491 | 16 | ||
33791.6579 | 1 | 41697.2980 | 9 | 42603.5100 | 9 | 43878.6130 | 9 | 44731.4494 | 16 | ||
secondary | |||||||||||
44134.5281 | 16 | 44278.6418 | 16 | 44644.5686 | 16 | 50273.4290 | 23 | 51166.5939 | 25 | 51577.5847 | 27 |
44134.5298 | 16 | 44278.6424 | 16 | 44644.5687 | 16 | 50849.2992 | 25 | 51550.3106 | 26 | 51924.5385 | 26 |
44134.5327 | 16 | 44644.5685 | 16 | 50257.4228 | 23 | 51158.2894 | 25 | 51566.3233 | 27 | 51997.4875 | 26 |
References: 1 - van Woerden (1957), 2 - Pohl & Kizilirmak (1966), 3 - Kizilirmak & Pohl (1969), 4 - Pohl & Kizilirmak (1970), 5 - Muthsam (1972), 6 - Hilditch et al. (1979), 7 - Kizilirmak & Pohl (1971), 8 - Kizilirmak & Pohl (1974), 9 - Patkos (1982a), 10 - Pohl & Kizilirmak (1975), 11 - Pohl & Kizilirmak (1976), 12 - Mallama et al. (1977), 13 - Pohl & Kizilirmak (1977), 14 - Mallama (1979), 15 - Pohl & Gülmen (1981), 16 - Milano et al. (1986), 17 - Pohl et al. (1982), 18 - Braune & Mundry (1982), 19 - Cellino et al. (1985), 20 - Pohl et al. (1987), 21 - Kundera (2001), 22 - Heckert (1996), 23 - Albayrak et al. (1996), 24 - Agerer & Hübscher (1998), 25 - Pribulla et al. (1999), 26 - Pribulla et al. (2001), 27 - this study.
Parameters | Sommer | Friebos-Conde | Hilditch | Cellino | Sarma | Albayrak | present |
(1956) | & Herczeg (1973) | et al. (1979) | et al. (1985) | et al. (1985) | et al. (1999) | study | |
P
![]() |
57.50 | 72.75 | 64.07 | 74.70 | 54.23 | 43.81 |
![]() |
A (days) | 0.0079 | 0.0085 |
![]() |
||||
e | 0 | 0 | 0.60 | 0.40 | 0.59 | 0.30 |
![]() |
w (![]() |
90 | 90 | 174 | 227 |
![]() |
||
f( M3) ( M![]() |
0.00114 | 0.00130 | 0.00163 | 0.00170 |
![]() |
||
M3 ( M![]() |
0.16a | 0.19a | 0.18a | 0.18a |
![]() |
To analyse the asymmetric light curves formed by the new observations, we used Djurasevic's (1992a) programme, which is based on the Roche model and the principles from the paper by Wilson & Devinney (1971). The light-curve analysis was made by applying the inverse-problem method (Djurasevic 1992b) based on Marquardt's (1963) algorithm.
The stellar size in the model is described
by the filling factors for the critical Roche lobes F1,2 of
the primary and secondary components, respectively, which tell us to what
degree the stars in the system fill their corresponding critical
lobes. For synchronous rotation of the components these factors
are expressed as the ratio of the stellar polar radii, R1,2,
and the corresponding polar radii of the critical Roche lobes, i.e.,
.
To achieve more reliable estimates of the model parameters in the
light-curve analysis programme, we applied a quite dense coordinate grid, having
individual elementary cells per star. The intensity and angular
distribution of the radiation of individual cells are determined by the stellar
effective temperature, limb-darkening, gravity-darkening and by the
reflection effect are in the system.
The asymmetries in the light curves are assumed to be formed by cool spots on
the hotter more massive component. In the code the spotted regions are
approximated by circular spots, characterised by the temperature contrast of
the spot with respect to the surrounding photosphere (
),
by the angular dimension (radius) of the spot (
)
and by the
longitude (
)
and latitude (
)
of the spot
centre. The longitude (
)
is measured clockwise (as viewed
from the direction of the +Z-axis) from the +X-axis (line connecting the stars'
centres) in the range
-
.
The latitude (
)
is measured from
at the stellar equator (orbital plane) to
towards the "north'' (+Z) and
towards the "south'' (-Z) pole.
For a successful application of this model in the analysis of the observed
light curves, the method of Djurasevic (1992b)
was used. Optimum model parameters are obtained through the minimization of
,
where
is the residual between the
observed (LCO) and synthetic (LCC) light curves for a given orbital
phase. The minimization of S is performed in an iterative cycle of corrections of
the model parameters. In this way the inverse-problem method gives us the
estimates of system parameters and their standard errors.
We considered that Patkos's (1982a) symmetrical reference light curve is very probably clean with respect to spot effects. The maximum light level of this reference light curve coincides with the light level of the secondary maximum in our light curves. Thus we decided to normalise our light curves with respect to the secondary maximum.
The mass-ratio of the components was fixed to a revised value of
estimated by Pojmanski
(1998) from the radial velocity solution. From his radial
velocity study, for the mass of the primary we have
.
The appropriate spectral type for this mass is
F5 V, with an effective temperature
K (Lang
1992).
This simple calibration may not work for spotted systems (in which
large cool spots probably affect the spectral type determination).
The light-curve analysis is made also under the assumption that
primary's spectral type is close to F8 V, as suggested by
Pojmanski (1998). Thus we used two different values for
the temperature of the hotter primary
component:
K (Hyp. I.) and
K (Hip. II.) corresponding to spectral types of F5 V
and F8 V, respectively.
The values of the limb-darkening coefficients were derived from the
stellar effective temperature and surface gravity, according to the
given spectral type, by using the polynomial proposed by
Díaz-Cordovés et al. (1995). During the process of
optimisation, with the temperature changes, we have an
automatic recomputation of the limb-darkening. For the R-filter the
limb-darkening was taken from the tables of Al Naimy
(1977). Following Lucy (1967), Rucinski
(1969) and Rafert & Twigg (1980), the
gravity-darkening coefficients of the stars,
,
and
their albedos,
,
were set at the values of 0.08 and 0.5,
respectively, appropriate for stars with convective envelopes.
Previous versions of our programme chose between two treatment of the radiation
law: simple black-body theory, or stellar atmosphere models by Carbon
& Gingerich (1969) (CG). The current version of the programme
(e.g. Djurasevic et al. 2001) employs the Basel Stellar
Library (BaSeL). We have explored the "corrected'' BaSeL model flux distributions,
consistent with extant empirical calibrations (Lejeune et al. 1997,
1998), and with a large range of effective temperatures
2000
K, surface gravities,
and metallicity,
,
where [Fe/H] is the logarithmic
metal abundance. The surface gravities can be derived very accurately from the
masses and radii of close binary (CB) stars by solving the inverse problem of
the light-curve analysis, but the temperature determination is related to the
assumed metallicity and strongly depends on the photometric calibration.
In the inverse problem the fluxes are calculated in each iteration for the
current values of temperatures and
by interpolation in both
of these quantities in the atmosphere tables, as an input, for a given metallicity
of the CB components. If there is a good physical reason, the
metallicity of the CB components could be different. But, normally, one uses
the same value for both components. The two-dimensional flux interpolation
in
and
is based on the application of the
bicubic spline interpolation (Press et al. 1992). This proved
to be a good choice.
By selecting the particular input switch, the programme for the
light-curve analysis can be simply redirected to the Planck or CG approximation,
or to the more realistic BaSeL model atmospheres. The disagreements between individual
B, V, and R solutions decrease if we introduce the "corrected''
BaSeL model flux distributions. A change in the assumed metallicity causes a
noticeable change in the predicted stellar effective temperature. The mean
value of the chemical abundances of the components was obtained
by experimenting with values close to solar metallicity. For SV Cam, the best
fit of the B, V, and R light curves was obtained with
for each components. With this
value, the individual B, V, and R solutions are in good
agreement. The results presented here use this stellar atmosphere
approximation.
In the analysis, the inclination of the orbit was estimated to be
.
The present light-curve analysis shows that during the
deeper (primary) minimum the cooler (less-massive and smaller) component
eclipses the hotter (more-massive and larger) one.
Since the stars in the system have external convective envelopes, which can exhibit magnetic activity, we started the "spotted solution'' by assuming that the components of SV Cam have cool spots, of the same nature as solar magnetic spots. Moreover, since the results of the light-curve analysis depend on the choice of the adopted working hypothesis, the analysis was carried out within the framework of several hypotheses with spotted areas on the components.
We rejected those hypotheses which produced significantly different values of
the parameters for the system and the active spotted areas, estimated by
analysing the individual light curves in the B, V, and R
passbands. Finally we chose the Roche model with dark spotted areas on the
more massive (hotter) component as the optimum solution. Within this
hypothesis the analysis of the light curves yields mutually consistent
parameters of the system and active region in the B, V, and
R passbands. In this case we obtained a very good fit to the
observations.
![]() |
- accepted metallicity of the components,
- temperature of the more-massive (hotter) star,
- nonsynchronous rotation coefficients of the components,
- mass ratio of the components,
- gravity-darkening coefficients of the components,
- albedo coefficients of the components,
- spots' temperature coefficient,
BaSeL approximation of the stellar atmospheres (
- accepted
metallicity of the components).
Note: n - number of observations,
-
final sum of squares of residuals between observed (LCO) and synthetic (LCC)
light curves,
- standard deviation of the observations,
,
and
- spots' angular dimensions, longitudes and latitudes (in arc degrees),
- filling factors for the critical Roche lobe of the
hotter (more-massive) and cooler (less-massive) star,
- temperature
of the less-massive cooler star, i - orbit inclination (in arc degrees),
- limb-darkening coefficients of the components,
- dimensionless surface potentials of the primary and
secondary,
- polar radii of the components in
units of the distance between the component centres,
- luminosity of the hotter star (including spots on the cooler one),
1 - absolute bolometric magnitudes of SV Cam components,
- stellar masses in solar units,
- mean radii of stars in solar units,
- logarithm (base 10) of the mean surface
acceleration (effective gravity) for system stars and
- orbital semi-major axis in units of
solar radius.
The light-time effect due to the orbit of the eclipsing pair around their
center of mass with a third body has been computed with our derived
parameters given in Table 3. The eclipsing pair completes a revolution on
this orbit in 41.32 yr. The projectional angular separation between the
third body and the eclipsing pair is
which is within
the observing limits of modern technology. We are certain that the presence of a
third body in the system of SV Cam is now well established. But further
observations of minima are needed for a more precise determination of
the second maxima of the (O-C) variation. It is desirable to try to confirm
detection by speckle interferometry which can support our position.
The parameters derived from the light curve analysis are listed in Table
4. The errors of the parameter estimates arise from the nonlinear
least squares method, on which the inverse problem method is based. The fixed
parameters and the meaning of the symbols in the first column are listed in a
footnote to Table 4. The spot characteristics (spot temperature
factor,
;
angular radius,
;
longitude,
;
latitude,
)
are also given in
Table 4. The determination of these parameters is based on a
simultaneous fitting of the available light curves in the B, V
and R photometric bands.
Finally, in Table 4 we provide some important absolute parameters.
They were derived from the revised mass ratio of the components
(
), orbital period (
),
and the semimajor axis (
)
found by
Pojmanski (1998) on the basis of the solution of the radial
velocity curve. Our estimate of the accuracy in the determination of these
parameters is based on the influence of formal errors arising from the nonlinear
method of the light-curve analysis, and also on the differences between
individual B, V and R solutions. The errors of
the input parameters of the model, are treated as fixed in the
inverse-problem method. They are taken from Pojmanski (1998).
Using the inverse-problem solutions for individual light curves, Fig. 2 (left) presents the optimum fit of the observed light curves (LCO) by the synthetic ones (LCC). The reference light level at phase 0.75, for the assumed unspotted configuration of the system is denoted by a dashed line. The final residuals (O-C) between the observed (LCO) and optimum synthetic (LCC) light curves are also given. A view of the Roche model of SV Cam obtained with the parameters estimated by analysing the new light curves (see Table 4) is provided in the lower right-hand panel of Fig. 2. Two cool spots were also shown on the surface of the hotter more massive primary component. We give only the graphical presentation of the results obtained within (Hyp. I.) since the differences in the fitting quality and estimated geometrical parameters of the system for both hypotheses are insignificant.
It is evident from Table 4 and on the left-hand panel of Fig.
2 that the Roche model with cool spot areas on the hotter
primary component gives a satisfactory fit to the analysed new light
curves. Thus the complex nature of the light curve
variations of SV Cam can be explained by cool spot areas on the hotter
primary component. The presence of two large cool spots leads to a
conspicuous asymmetry in the light curves. Since the system's period is
short (P
0
59), the presence of spots at high
latitudes (near the polar regions) can be explained by the dynamo
mechanism for rapid rotators (Schüssler & Solanski
1992). It is possible that large spot areas at high latitudes
may correspond to an enhanced activity state of the system.
Acknowledgements
This work was supported partly by the Scientific and Thechnical Research Council of Turkey (TÜBITAK) and the Research Foundation of Çanakkale University. We appreciate the anonymous referee for her/his constructive and generous comments, Saul J. Adelman for improvements of English. We also thank the TÜBITAK National Observatory for observing time and to K. Yüce for her assistance during the observations.