A&A 487, 1081-1093 (2008)
DOI: 10.1051/0004-6361:200809670

Four-colour photometry of eclipsing binaries. XLI.

uvby light curves for AD Bootis, HW Canis Majoris, SW Canis Majoris, V636 Centauri, VZ Hydrae, and WZ Ophiuchi[*],[*],[*]

J. V. Clausen1 - L. P. R. Vaz2 - J. M. García3 - A. Giménez4,5 - B. E. Helt1 - E. H. Olsen1 - J. Andersen1,6


1 - Niels Bohr Institute, Copenhagen University, Juliane Maries Vej 30, 2100 Copenhagen Ø, Denmark
2 - Depto. de Física, ICEx, UFMG, CP 702 - 30.123-970 - Belo Horizonte, MG, Brazil
3 - Departamento de Física, E.U.I.T. Industriales, UPM, Ronda de Valencia 3, 28012 Madrid, Spain
4 - Laboratorio de Astrofísica Espacial y Física Fundamental, INTA, Apartado 50727, 28080 Madrid, Spain
5 - Research and Scientific Support Department, ESA, ESTEC, 2200 AG Nordwijk, The Netherlands
6 - Nordic Optical Telescope Scientific Association, Apartado 474, 38 700 Santa Cruz de La Palma, Spain

Received 28 February 2008 / Accepted 25 May 2008

Abstract
Context. Accurate mass, radius, and abundance determinations from binaries provide important information on stellar evolution, fundamental to central fields in modern astrophysics and cosmology.
Aims. Within the long-term Copenhagen Binary Project, we aim to obtain high-quality light curves and standard photometry for double-lined detached eclipsing binaries with late A, F, and G type main-sequence components, needed for the determination of accurate absolute dimensions and abundances, and for detailed comparisons with results from recent stellar evolutionary models.
Methods. Between March 1985 and July 2007, we carried out photometric observations of AD Boo, HW CMa, SW CMa, 636, VZ Hya, and WZ Oph at the Strömgren Automatic Telescope at ESO, La Silla.
Results. We obtained complete ${\it uvby}$ light curves, ephemerides, and standard ${\it uvby}\beta $ indices for all six systems. For V636 Cen and HW CMa, we present the first modern light curves, whereas for AD Boo, SW CMa, VZ Hya, and WZ Oph, they are both more accurate and more complete than earlier data. Due to a high orbital eccentricity (e = 0.50), combined with a low orbital inclination ( $i = 84\hbox{$.\!\!^\circ$ }7$), only one eclipse, close to periastron, occurs for HW CMa. For the two other eccentric systems, V636 Cen (e = 0.134) and SW CMa (e = 0.316), apsidal motion has been detected with periods of $5270 \pm 335$ and $14~900 \pm 3600$ years, respectively.

Key words: stars: evolution - stars: fundamental parameters - stars: binaries: close - stars: binaries: eclipsing - techniques: photometric

   
1 Introduction

Detached, double-lined eclipsing binaries are well-known as the main source of accurate (1-2% or even better) data on stellar masses and radii. During the last two-three decades, such results have been published for about 50 systems within the main-sequence band (e.g. Andersen 1991; Ribas et al. 2000a; Southworth, private communication[*]), but some mass intervals, especially around and below 1 $M_{\odot}$, still need better coverage, as do parts of the main-sequence band at higher masses. A rich sample of accurate binary results is mandatory for detailed tests of stellar models in general, as well as for many specific astrophysical investigations (Ribas 2004, 2006b). Remarkable recent examples are the determination of the helium-to-metal enrichment ratio (Ribas et al. 2000a) and investigations of the mass dependence of convective core overshooting (Ribas et al. 2000b; Claret 2007).

In this paper, we present new accurate and complete ${\it uvby}$ light curves for six detached double-lined eclipsing binaries having late A, F, and G-type main-sequence components. In addition, new spectroscopic observations for radial velocities and abundance analyses have been obtained, and detailed analyses for all the systems will be published separately (Clausen et al. 2008; Clausen et al., in prep.; Torres et al. in prep.).

When we decided to observe AD Boo (F6+G0) and SW CMa (A8), they were lacking modern photometric and spectroscopic observations. This was independently noticed by Lacy (1997a,b) who studied both systems, and by Popper (1998a) who obtained spectroscopic orbits for AD Boo, leading to masses for the components which differ somewhat from those by Lacy. For both systems, a new analysis based on new light curve and radial velocity observations, supplemented with a spectroscopic abundance analysis, is desirable.

HW CMa (A8) is only about $2 \hbox{$^\prime$ }$ distant from SW CMa and was discovered to be a double-lined spectroscopic binary during the radial velocity observations of SW CMa. It was later found to be eclipsing (Liu et al. 1992), and we decided to obtain complete light curves and spectroscopic orbits.

V636 Cen (F8/G0V) was selected as part of a large-scale study of new systems with solar-type components (Clausen et al. 2001). We want e.g. to investigate a serious dilemma that appears to be present in the comparison of predictions from current stellar models with fundamental properties of known 0.7-1.1 $M_{\odot}$ eclipsing binaries (Popper 1997; Clausen et al. 1999b). Our results for V636 Cen support the presence of a real discrepancy - current models, adopting a fixed mixing length parameter calibrated on the Sun for envelope convection, do not predict identical ages for the two components (Larsen 1998; Clausen et al., in prep.). A similar case, V1061 Cyg, was presented by Torres et al. (2006). See also Ribas (2006a) for eclipsing binaries with K and M type components and Morales et al. (2007) for single active late-K and M stars.

VZ Hya (F5) and WZ Oph (F8) are well known systems included by Popper (1980) in his critical review of stellar masses, but excluded by Andersen (1991), since the available results did not fulfill his 2% criteria for masses and radii. We therefore decided to provide the new observations needed for the determination of more accurate dimensions.

In the following, we first give a general description of the observations, the basic photometric reduction, and the light curve formation (Sect. 2), and we then present the results for the individual systems (Sects. 3-8). Throughout the paper, the component eclipsed at the deeper eclipse at phase 0.0 is referred to as the primary, and the other as the secondary component.

   
2 Observations and photometric reduction

The differential ${\it uvby}$ light curve observations were obtained with the Strömgren Automatic Telescope (SAT) at ESO, La Silla and its 6-channel ${\it uvby}\beta $ photometer during several campaigns between March 1985 and March 2002. For SW CMa and V 636 Cen, more eclipse observations were carried out between January 2001 and July 2007 to follow the apsidal motion of their orbits. In addition, on many nights we observed the binaries and the comparison stars together with a large sample of ${\it uvby}$ and $\beta $ standard stars in order to obtain homogeneous standard ${\it uvby}\beta $ indices, primarily for the determination of effective temperatures and interstellar reddening. Details on the photometer and the (semi)automatic mode of the telescope are given by Olsen (1993, 1994b).

   
2.1 Light curve observations

Table 1: rms errors of the magnitude differences (instrumental system) between the comparison stars in units of 0.1 mmag. N is the total number of magnitude differences.

The differential ${\it uvby}$ observations were, depending on the orbital phase of the binary, carried out either continuously over several hours or as shorter series. Automatic centering within a circular diaphragm of $13\hbox{$^{\prime\prime}$ }$ (HW CMa, part of light curves; SW CMa, light curves; V636 Cen, light curves) or $17\hbox{$^{\prime\prime}$ }$ (AD Boo, light curves; HW CMa, part of light curves; SW CMa, additional eclipse observations; V636 Cen, additional eclipse observations; VZ Hya, light curves; WZ Oph, light curves) diameter was used throughout. Nearly all observations were performed at an airmass less than 2.1, and the few at higher airmass passed additional quality checks.

Table 2: New photometric data for the eclipsing binaries and the comparison stars. For the eclipsing binaries, the ${\it uvby}\beta $ information is the mean values outside eclipses, see also Fig 1. N is the total number of observations used to form the mean values, and $\sigma $ is the rms error (per observation) in mmag. See also Table 3.

Table 3: Published photometric data for the eclipsing binaries and the comparison stars. N is the total number of observations used to form the mean values, and $\sigma $ is the rms error (per observation) in mmag.

Two or three comparison (C) stars, preferably matching the spectral type of the binary and positioned on the sky within a few degrees from it, were adopted for each candidate. They were in general selected from systematic searches in the extensive ${\it uvby}\beta $ catalogues published by Olsen (1983, 1988, 1993, 1994a,b), and unpublished $\beta $ photometry (cf. Olsen 1994b), avoiding known and suspected variables, as well as most visual binaries.

The observations were carried out in C1-binary-C2-binary-C3-binary-C1-etc. or C1-binary-C2-binary-C1-etc. sequences, allowing accurate magnitude differences to be formed, even if one of the comparisons were later found to be variable. In the case of HW CMa and SW CMa, which are close on the sky and of similar spectral types, we often used C1-SW CMa-C2-HW CMa-C1-etc. An observation consisted in general of three individual integrations, each with an integration time of 10-60 s, and sky measurements were taken at least once per sequence, normally at a fixed position close to the binary. The resulting rms contributions from photon statistics, including the sky contributions, were in general kept at 5 mmag or lower. The Heliocentric Julian Date (HJD) given for an observation refers to the midpoint of the time interval covered by the integrations.

   
2.2 Photometric reduction and formation of light curves

Linear extinction coefficients were determined individually for each night from the observations of comparison stars and other constant stars. Whenever appropriate, linear or quadratic corrections for drift during the night, caused by changes in the sky transparency and/or the influence of temperature variations on the uncooled photomultipliers, were also applied.

The ${\it uvby}$ instrumental system of the SAT, where the same spectrometer, including its filters, most of the photomultipliers, and the photon counting system have been used since 1985, has proved to be very stable (e.g. Olsen et al., in prep.). Therefore, transformation of the light curve observations, e.g. to the standard ${\it uvby}$ system, is not needed, and the light curves are generated in the instrumental system.

Differential magnitudes were formed for each binary observation from the two comparison star observations closest in time. All comparison star observations were used with C2 and C3, if observed, first shifted to the level of C1. A careful check of the constancy of all two/three stars was performed. For each binary, typical rms errors per light curve point are close to those listed in Table 1 for the magnitude differences between the corresponding comparison stars.

Ephemerides were calculated from published and new times of minima, listed in Tables A.1-A.8. For most of the systems, our literature searches were significantly supplemented by the compilations by Kreiner (private communication); see Kreiner et al. (2001)[*] for further details. The new times of minima have been calculated using the method of Kwee & van Woerden (1956), except for a few cases which could only be determined by 2nd order polynomial fits to the (few) observations.

The individual light curves are shown in Figs. 2-46, 8, and 9, and are presented in Tables 6-11. The additional eclipse observations of SW CMa and 636 are presented in Tables 12 and 13. These eight tables will only be made available in electronic form.

  \begin{figure}
\par\includegraphics[width=16cm,clip]{9670f1.ps}
\end{figure} Figure 1: Phase distribution of the $\beta $ observations, see Table 2. The horizontal lines represent the mean value of $\beta $ outside eclipses (solid) and the rms error of one observation (dotted). The vertical lines show the phase intervals of the eclipses (dashed, central eclipses) and their duration (dotted).
Open with DEXTER

   
2.3 Standard uvby$\beta $ photometry

New standard ${\it uvby}\beta $ indices for the six eclipsing binaries outside eclipses and for all comparison stars are listed in Table 2, and the phase distributions of the $\beta $ observations are shown in Fig. 1.

Part of the standard photometry was obtained at SAT on several nights in 1987-1993, where the binaries and the comparison stars were observed together with a large sample of ${\it uvby}$ and $\beta $ standard stars. The basic reduction and transformation to the standard system was done as described by Olsen (e.g. 1993, 1994b). Furthermore, some of the binaries and their comparisons were also included in the long term search for new solar-type binaries, carried out at SAT since 1994, and additional standard photometry from many nights has been obtained from this project; see Olsen et al. (in prep.) for details.

For comparison, we have listed published ${\it uvby}\beta $ indices for some of the binaries and their comparison stars in Table 3. Individual differences larger than the quoted errors occur; we have used the new homogeneous results in Table 2 for the binary analyses.

   
3 AD Boo

AD  Boo = BD +25$^\circ$2800 = BV 135 was discovered by Strohmeier et al. (1956) to be an eclipsing binary, and times of minima, a linear ephemeris, and a photographic light curve were published by Strohmeier et al. (1963; orbital period half the correct value). A large number of times of minima of AD  Boo has later been observed (see Tables A.1 and A.2), and an improved ephemeris and a visual light curve was published by van Buren (1974; orbital period half the correct value). Koch (1974) presented blue CN-absorption measurements, and Zhai et al. (1982, 1983) obtained photoelectric B,V light curves and derived the true orbital period of about $2\hbox{$.\!\!^{\rm d}$ }07$.

Lacy (1985) discovered AD Boo to be a double-lined spectroscopic binary, and later (Lacy 1997a) he presented absolute dimensions based on spectroscopic elements from 23 KPNO CCD spectra and photometric elements derived from EBOP (Popper & Etzel 1981) re-analyses of the light curves (normal points) by Zhai et al. (1982). AD Boo was included by Popper in his extensive study of F-K binaries (e.g. Popper 1993; Popper & Jeong 1994; Popper 1996), and he obtained a very accurate spectroscopic orbit based on 31 Lick CCD echelle spectra (Popper 1998a). He also pointed out that the components of AD Boo have the largest differences in mass and radius among well-studied systems in its mass range (1.4 and 1.2 $M_{\odot}$, 1.6 and 1.2 $R_{\odot}$). Hence, it may provide stringent tests of stellar models.

   
3.1 uvby light curves

Complete ${\it uvby}$ light curves containing 652 points in each colour were observed on 42 nights during five periods between March 1988 and March 1992 (JD 2 447 234-JD 2 448 695). HD 128185, HD 128369, and HD 129430 = HR 5483 were selected as comparison stars; see Tables 1-3 for further information. The last two stars were found to be constant within the observational accuracy of a few mmag during our observations, whereas signs of a slight long term variability at the level of about $\pm$5 mmag were noticed for HD 128185, and it was consequently decided to reject it for the calculations of the final light curves. We have not studied the possible variations of HD 128185 further; our photometric data for it are available on request.

The light curves of AD  Boo are shown in Fig. 2. The eclipses have been covered several times and most out-of-eclipse phases at least twice. The accuracy per point is about 0.005 mag ($\it vby$) and 0.007 mag (u), but comparison of the data from the five periods reveals that at some phases, systematic differences in the light level exist, of the order of 0.01-0.03 mag, and increasing from y to u. Especially the observations during eclipses from 1989 are systematically fainter. The scatter is probably caused by slight activity of the cooler secondary component, since a weak emission feature in the Ca II H and K lines is seen at the position of the secondary component on high-resolution spectra of AD Boo (Clausen et al. 2008). Zhai et al. (1982) also noticed scatter in their B and V light curves, which were obtained from April to June 1981, but at a much higher level of about 0.075 mag, compared to their observing error of about 0.02 mag.

The new ${\it uvby}$ light curves confirm that AD Boo consists of two well-detached and rather different stars in a circular orbit, and our photometric analysis reveals that the secondary eclipse is in fact total.

  \begin{figure}
\par\includegraphics[width=16cm,clip]{9670f2.ps}
\end{figure} Figure 2: y light curve and b-y and u-b colour curves (instrumental system) for AD Boo.
Open with DEXTER


  \begin{figure}
\par\includegraphics[width=16.5cm,clip]{9670f3.ps}
\end{figure} Figure 3: y light curve and b-y and u-b colour curves (instrumental system) for HW CMa.
Open with DEXTER

   
3.2 Ephemeris

Four times of primary minimum and two of secondary, derived from the ${\it uvby}$ observations, are listed in Tables A.1 and A.2 together with 130 published times.

The linear ephemeris given in Eq. (1) was computed by least squares based on all available visual, photographic, and photoelectric/CCD times of eclipses, assuming a circular orbit.

 \begin{displaymath}
\begin{tabular}{r r c r r}
{\rm Min ~ I} = & 2~449~311.43169...
...ace \\
& $\pm $$22$ & & $\pm $$ 17$\space & \\
\end{tabular}\end{displaymath} (1)

Weights were assigned to the observations based on the individual uncertainties, when available. Typically the visual and photographic minima have no published uncertainties, so those were determined iteratively based on the residuals from the fit, in such a way as to achieve a reduced $\chi^2$ of unity for each type of observation. Similarly, scale factors were determined for the photoelectric/CCD measurements so as to yield a reduced $\chi^2$ of unity, separately for the primary and secondary minima, since they can have different precisions due to the different depth of the eclipses.

Separate least squares fits to all primary and secondary times yield identical periods within errors, and again within errors the same ephemeris is obtained from the subset of photoelectric and CCD data. The method of Lafler & Kinman (1965) applied to the ${\it uvby}$eclipse observations also confirms the period.

Zhai et al. (1983) derived a slightly longer period of $2 \hbox{$.\!\!^{\rm d}$ }0688112 \pm 0.0000002$ from their times of minima and those listed by van Buren (1974), and their ephemeris was used by Lacy (1997a) and Popper (1998a). We see, however, no clear signs of period changes from the data listed in Tables A.1 and A.2.

   
4 HW CMa

HW CMa = HD 54549 was discovered by Liu et al. (1992) to be a double-lined eclipsing binary in a highly eccentric orbit ( $P = 21\hbox{$.\!\!^{\rm d}$ }1$, e = 0.5). Minimum masses of 1.74 and 1.80 $M_{\odot}$ were determined, and a shallow eclipse (about 0.13 mag deep) lasting for about 5 h (0.010 in phase) was detected at conjunction, which occurs close to periastron. The star eclipsed is the less massive one. At apastron, the separation between the components is about three times larger, and no eclipse is seen around the predicted phase (0.517). This is due to the high orbital eccentricity, combined with a low orbital inclination ( $i = 84\hbox{$.\!\!^\circ$ }7$). Houk & Smith-Moore (1988) classify HW CMa as A1 III, but CfA spectra and ${\it uvby}$ photometry reveal that the components are almost identical unevolved main-sequence stars near A8.

   
4.1 uvby light curves

The ${\it uvby}$ light curves, which contain 414 points in each colour, were observed on 53 nights between February 1989 and March 2002 (JD 2 447 576-JD 2 452 335). The observations by Liu et al. (1992), with the same system, have been included. HW CMa is only about $2 \hbox{$^\prime$ }$ distant from SW CMa and of similar spectral type. Therefore, we have used the same comparison stars, HD 56341 = HR 2755 and HD 53123, and as mentioned in Sect. 2.1 HW CMa was often observed in series which also included SW CMa. For the October 2001-March 2002 period, a third star HD55344 was added as a further check of the comparison stars, which were found to be constant within the observational accuracy; see Tables 1-3 for further information.

The light curves are shown in Fig. 3. The eclipse at phase 0.0 was observed on three nights, and the expected interval near phase 0.517 for the other eclipse, which should last significantly longer, was covered closely on one night with a few additional observations from other nights. As mentioned above, no sign of an eclipse is, however, seen. As expected, the light curves are essentially constant outside eclipse, and no significant colour changes are seen during the eclipse, where $\sim$20% of the light of the primary component is blocked.


  \begin{figure}
\par\includegraphics[width=16.5cm,clip]{9670f4.ps}
\end{figure} Figure 4: y light curve and b-y and u-b colour curves (instrumental system) for SW CMa.
Open with DEXTER

   
4.2 Ephemeris

Two minimum times, separated by 170 epochs, have been determined from the ${\it uvby}$ observations; see Table A.3. This limited material yields the following linear ephemeris for HW CMa:

 \begin{displaymath}
\begin{tabular}{r r c r r}
{\rm Min ~ I} = & 2~452~279.6787 ...
... E$\space \\
& $\pm 4$ & & $\pm 33$\space & \\
\end{tabular}\end{displaymath} (2)

For comparison, Liu et al. (1992) derived a period of $21\hbox{$.\!\!^{\rm d}$ }1178 \pm 0\hbox{$.\!\!^{\rm d}$ }0005$, and the spectroscopic orbit based on the CfA velocities, which cover 178 cycles, gives $21\hbox{$.\!\!^{\rm d}$ }11790 \pm 0\hbox{$.\!\!^{\rm d}$ }00025$. Since the orbit is highly eccentric, apsidal motion is expected to be present, although with a long period due to the large separation between the components and the correspondingly small relative radii. The orbital parameters and estimated stellar dimensions yield a theoretically predicted period of about 90 000 yr.

   
5 SW CMa

SW CMa = HD 54520 = HIP 34431 was discovered by Hoffmeister (1932) to be an eclipsing binary, and Florja (1937) observed it visually and determined its ephemeris (Min = $2~426~709.324 + 10 \hbox{$.\!\!^{\rm d}$ }092 \times E$). The period is close to the correct value, but today we know that the epoch corresponds to secondary eclipse. Floja's visual light curve shows only one eclipse, and he therefore speculates whether one of the only two times where Hoffmeister found a minimum, which does not fit the ephemeris, is wrong. Again, today we know that this observation is correct and corresponds to the primary minimum.

Struve (1945) found SW CMa to be a double-lined spectroscopic binary with $\approx$ $2~M_{\odot}$ components in a very eccentric orbit (e = 0.5; significantly higher than the correct value of about 0.3). Later, a photographic light curve and times of minima, based on the old Sonneberg plates used by Hoffmeister, were presented by Wenzel (1952), and several additional times of minima are available in the literature (see Table A.4).

Popper (1965) mentioned SW CMa as one of several southern late B to early A systems that should be studied, and, using modern detectors, Lacy (1984) found narrow, deep lines of nearly equal strength for the components, and encouraged photoelectric observers to obtain accurate light curves. ${\it uvby}\beta $ indices have been published by Wolf & Kern (1983), and $\it UBV$ indices by Lacy (1992). SW CMa was included by Giménez (1994) in his list of eccentric binaries that should be monitored photometrically in order to study apsidal motion, and finally Lacy (1997b) determined absolute dimensions from analyses of $\it UBV$ light curves observed at CTIO 1993-95 and coudé spectra observed at McDonald and KPNO 1982-89. His analysis reveals that SW CMa consists of two somewhat different components (2.2 and 2.0 $M_{\odot}$, 3.0 and 2.5 $R_{\odot}$) which have evolved to the upper part of the main-sequence band.

   
5.1 uvby light curves

Complete ${\it uvby}$ light curves containing 820 points in each colour were observed on 89 nights during five periods between February 1987 and March 1992 (JD 2 446 825-JD 2 448 689). As mentioned in Sect. 2.1, SW CMa was often observed in series which also included HW CMa. For the 1987 observations, we selected HD 55271 and HD 56341 = HR 2755 as comparison stars, and in addition HD 53123 was observed frequently on each night. None of the stars shows any sign of variability within the observational accuracy of a few mmag; see Tables 1 and 2. However, HD 55271 = ADS 5863A = WDS 07113-2148A (A-BC: $\rho = 13 \hbox{$.\!\!^{\prime\prime}$ }6$, $\Delta m = 1.3$, Worley & Douglass 1997) is a member of a double star, and automatic centering of it proved to be less precise than desirable due to the companion. HD 55271 was therefore rejected as comparison for the 1988-92 observations and replaced by HD 53123. We have used the HD 56341 and HD53123 observations for calculation of the light curves.

The new ${\it uvby}$ light curves, which are shown in Fig. 4, contain about twice the number of points per band as the $\it UBV$ light curves by Lacy (1997b), and they are of at least the same photometric quality. Furthermore, they cover the orbit better, especially with respect to the shoulders of the minima, and this allows more accurate photometric elements to be determined. The eclipses have been covered several times and most out-of-eclipse phases at least twice. In addition, eclipse observations were made January-February 2002, October 2005, and December 2006 in order to establish more times of minima and investigate whether apsidal motion could be detected; see below.

   
5.2 Ephemeris and apsidal motion

Five times of primary minimum and six of secondary minimum, derived from the ${\it uvby}$ observations, are given in Table A.4 together with 26 published times. Weighted (weight proportional to the square of the inverse rms) least squares fits to all times of primary and secondary minima yield nearly identical periods of $10\hbox{$.\!\!^{\rm d}$ }091985 \pm 0\hbox{$.\!\!^{\rm d}$ }000002$ and $10\hbox{$.\!\!^{\rm d}$ }091981 \pm 0\hbox{$.\!\!^{\rm d}$ }000003$, respectively, whereas adopting equal weights leads to slightly larger values of $10\hbox{$.\!\!^{\rm d}$ }092003 \pm 0\hbox{$.\!\!^{\rm d}$ }000008$ and $10\hbox{$.\!\!^{\rm d}$ }091999 \pm 0\hbox{$.\!\!^{\rm d}$ }000005$. However, the derived linear ephemerides do not represent the observed times very well. Ziegler (1965) derived $10 \hbox{$.\!\!^{\rm d}$ }091948$, which was adopted by Lacy (1997b). If only the photoelectric times of minima are fitted, periods of $10\hbox{$.\!\!^{\rm d}$ }091983 \pm 0\hbox{$.\!\!^{\rm d}$ }000001$ and $10\hbox{$.\!\!^{\rm d}$ }091977 \pm 0\hbox{$.\!\!^{\rm d}$ }000002$are obtained from primary and secondary eclipses, respectively.

Table 4: Apsidal motion parameters for SW CMa derived from the photoelectric times of minima.

Since the orbit of SW CMA is eccentric ( $e \approx 0.3$), apsidal motion should be present, as also discussed by Lacy (1997b), who found a poorly determined apsidal motion period of several thousand years. We now have more precise times of minima available and derive the apsidal motion parameters presented in Table 4 from a weighted least squares method, following the formalism by Giménez & Garcia-Pelayo (1983) and Giménez & Bastero (1995). The orbital inclination i and eccentricity e were fixed to the values derived from the photometric analysis (Torres et al., in prep.). As seen, also in Fig. 5, a slow but significant motion has been detected, but the apsidal motion period is still very uncertain, although the value obtained is in good agreement with predictions from stellar models. For more precise apsidal motion information, SW CMa should be monitored regularly in the coming decades.

The detected apsidal motion leads to an insignificant phase shift of the secondary eclipse compared to the primary eclipse of only -0.00012 during the 1987-92 ${\it uvby}$ light curve observations. Therefore, a linear ephemeris can safely be used for the light curve analysis, and we adopt


 \begin{displaymath}
\begin{tabular}{r r c r r}
{\rm Min ~ I} = & 2~446~829.6482 ...
...; E$\space \\
& $\pm 1$ & & $\pm 5$\space & \\
\end{tabular}\end{displaymath} (3)

which has been derived from the 1982-87 times of minima. Secondary eclipses occur at phase 0.309. For comparison, the method of Lafler & Kinman (1965) applied to the 1987-92 ${\it uvby}$ eclipse observations gives a period of $10 \hbox{$.\!\!^{\rm d}$ }091982$.

The 2002-2006 ${\it uvby}$ eclipse observations have not been included in the light curve analyses, since their phase shifts are larger.


  \begin{figure}
\par\includegraphics[width=7cm,clip]{9670f5.ps}
\end{figure} Figure 5: Apsidal motion for SW CMa. O-C are the residuals (days) from the linear part of the apsidal motion ephemeris defined by the parameters given in Table 4. The full and dashed curves represent predictions from the apsidal motion parameters for primary (triangles) and secondary (squares) eclipses, respectively. Only the photoelectric observations are shown.
Open with DEXTER

   
6 V636 Cen


  \begin{figure}
\par\includegraphics[width=16.5cm,clip]{9670f6.ps}
\end{figure} Figure 6: y light curve and b-y and u-b colour curves (instrumental system) for V636 Cen.
Open with DEXTER

From 276 visual (1952-53) and 282 photographic (1935-38, 1952-53) observations, Hoffmeister (1958) reported V636 Cen = HD 124784 = HIP 69781 to be a $4 \hbox{$.\!\!^{\rm d}$ }3$ period G0 eclipsing binary, and he gave a linear ephemeris and seven times of minima. Later, Popper (1966) obtained two 22 Å/mm spectra and found sharp but single lines. Today V636 Cen is known to be double-lined, but the lines of the secondary component are much fainter than those from the primary. Dworak (1973) pointed out that V636 Cen is within 100 pc of the Sun, and it was included in the list by Dworak & Oblak (1989) of 95 such eclipsing binaries selected for Hipparcos observations. The parallax determined by Hipparcos (ESA 1997) is $15.36 \pm 1.12$ mas, corresponding to a distance of only 65 pc. With accurate absolute dimensions becoming available, V636 Cen can therefore be added to the (short) list of eclipsing binaries useful for the definition of the radiative flux scale (Popper 1998b).

A preliminary analysis of V636 Cen, based on the ${\it uvby}$ light curves presented here and complete spectroscopic orbits from CORAVEL radial velocity measurements has been done by Larsen (1998). She finds that V636 Cen is one of the few well studied systems with components in the $0.7{-}1.1~M_{\odot}$ range; preliminary masses are 1.04 and 0.85 $M_{\odot}$, respectively. Furthermore, V636 Cen is well detached, and such systems are important for detailed tests of stellar models around and especially below 1 $M_{\odot}$.

   
6.1 uvby light curves

${\it uvby}$ light curves, containing in total 853 points in each colour, were observed on 76 nights during six periods between March 1985 and April 1991 (JD 2 446 144-JD 2 448 377). HD 124829 and HD 125444 were selected as comparison stars, see Tables 1-3 for further information, and they were both found to be constant within the observational accuracy of a few mmag during our observations.

The light curves of V636 Cen are shown in Fig. 6, and variations due to surface activity (from spots) are clearly seen from the ${\it uvby}$ observations. However, the large number of observations obtained during some of the seasons allows a detailed study of the activity and its influence on the accuracy of the photometric elements (Larsen 1998). As seen, V636 Cen consists of two well detached components of very different surface fluxes in an eccentric orbit ( $e \approx 0.13$).

Further eclipse observation were made January 2002- July 2007. They have been used in the apsidal motion analysis presented below, but not in the light curve analyses.

   
6.2 Ephemeris and apsidal motion

Nine times of primary and nine of secondary minimum, derived from the ${\it uvby}$ observations, are given in Table A.5 together with eight published times of primary minimum (none available for the secondary minimum). Weighted (weight proportional to square of inverse rms) least squares fits to all times of primary and secondary minima yield significantly different periods of $4\hbox{$.\!\!^{\rm d}$ }28394311 \pm 0\hbox{$.\!\!^{\rm d}$ }00000013$ and $4\hbox{$.\!\!^{\rm d}$ }28394845 \pm 0\hbox{$.\!\!^{\rm d}$ }00000019$, respectively. Fitting only the ${\it uvby}$ times of primary minima leads to the same period.

As mentioned in Sect. 6.1, the orbit of V636 Cen is eccentric (e = 0.13), and the secondary eclipse occurs near phase 0.52. Apsidal motion must therefore exist, although at a slow rate, since the relative radii of the components are small. From theoretical density concentration coefficients, the absolute dimensions, the eccentricity and inclination of the orbit, and assuming the rotation of the components to be pseudosynchronized, Larsen (1998) derived an expected apsidal motion period of $U \approx 6000$ yr. The relativistic contribution and the classical contribution from tidal effects are of almost the same size ( $P/U \approx 10^{-6}$), making V636 Cen particularly interesting among apsidal motion systems.

With additional ${\it uvby}$ times of minima available, we have repeated the apsidal motion analysis and obtain the parameters listed in Table 5; see also Fig. 7. Well defined apsidal motion parameters have now been determined for meaningful comparison with theoretical density concentrations of stellar models.

The apsidal motion does not lead to significant phase shifts of secondary eclipse compared to primary eclipse during the 1985-1991 ${\it uvby}$ light curve observations. We have adopted the following linear ephemeris, derived from the corresponding times of primary minima for the light curve analyses:

 \begin{displaymath}
\begin{tabular}{r r c r r}
{\rm Min ~ I} = & 2~446~873.80124...
...ce &$\times\; E$\\
& $\pm 14$ & & $\pm 7$ & \\
\end{tabular}\end{displaymath} (4)

Table 5: Apsidal motion parameters for V636 Cen derived from the photoelectric times of minima.


  \begin{figure}
\par\includegraphics[width=7.5cm,clip]{9670f7.ps}
\end{figure} Figure 7: Apsidal motion for V636 Cen. O-C are the residuals (days) from the linear part of the apsidal motion ephemeris defined by the parameters given in Table 5. The full and dashed curves represent predictions from the apsidal motion parameters for primary (triangles) and secondary (squares) eclipses, respectively. Only the ${\it uvby}$ observations are shown.
Open with DEXTER

   
7 VZ Hya


  \begin{figure}
\par\includegraphics[width=16.5cm,clip]{9670f8.ps}
\end{figure} Figure 8: y light curve and b-y and u-b colour curves (instrumental system) for VZ Hya.
Open with DEXTER

VZ Hya = HD 72257 = HIP 41834 was discovered by O'Connell (1932) to be an eclipsing binary, and spectroscopic elements were obtained from 75 Å/mm spectra by Struve (1945), who showed that the period is twice the value given by O'Connell. Early light curves and absolute dimensions were published by Wood (1946, visual data) as part of his thesis, and by Gaposchkin (1953, photographic data). From 20 and 40 Å/mm spectra, Popper (1965) later determined improved elements, leading to masses accurate to about 2.5%.

The first photoelectric ($\it UBV$) light curves were presented by Walker (1970) together with a preliminary photometric solution, whereas Wood (1971) performed a complete analysis as part of the introduction of his then new binary model and the corresponding WINK code. A photometric analysis, also based on WINK, was published by Cester et al. (1978).

Another set of photoelectric $\it UBV$ light curves were published by Padalia & Srivastava (1975). In contrast to Walker, Wood, and Cester et al. they found primary eclipse to be an occultation (larger star in front) rather than a transit (smaller star in front). This picture of VZ Hya was soon shown by Popper (1976) to be wrong, but it was nevertheless defended once more by Padalia (1986).

Because its dimensions were relatively well known, VZ Hya has often been included in binary samples used in various astrophysical investigations, e.g. He abundance in Population I stars (Popper et al. 1970), stellar radii (Lacy 1977; Shallis & Blackwell 1980; Pastori et al. 1985), and mass-radius relations (Zhai & Zhang 1989). The currently most reliable absolute dimensions are those given by Popper (1980) in his critical review of stellar masses. There is, however, still ample space for improvement, provided more accurate light curves and radial velocities are obtained. Since VZ Hya is an important F-type system, we decided to obtain the observations needed for this purpose.

   
7.1 uvby light curves

Complete ${\it uvby}$ light curves containing 1180 points in each colour were observed on 44 nights during four periods between February 1989 and April 1992 (JD 2 447 575-JD 2 448 684). HD 72528, HD 71615, and HD 72782 were used as comparison stars and were found to be constant within the observational accuracy of a few mmag during our observations; see Tables 1-3 for further information.

The light curves of VZ Hya are shown in Fig. 8. As seen, VZ Hya consists of two well-detached components of different surface fluxes in a circular orbit. Both eclipses have been covered several times and most out-of-eclipse phases at least twice. The accuracy per point is about 0.004 mag (vby) and 0.006 mag (u), meaning that the new light curves are of significantly higher quality than those by Walker (1970) and Padalia & Srivastava (1975). We see no signs of intrinsic variability. Photometric analyses reveal that the secondary eclipse is very close to being total (Clausen et al. 2008).


  \begin{figure}
\par\includegraphics[width=16.5cm,clip]{9670f9.ps}
\end{figure} Figure 9: y light curve and b-y and u-b colour curves (instrumental system) for WZ Oph.
Open with DEXTER

   
7.2 Ephemeris

Three times of primary minimum and one of secondary minimum, derived from the ${\it uvby}$ observations, are given in Table A.6 together with 26 published times and times redetermined from the photometry by Walker (1970).

We adopt the following linear ephemeris, which is based on all available times assuming a circular orbit and derived as described in Sect. 3.2:


 \begin{displaymath}
\begin{tabular}{r r c r r}
{\rm Min ~ I} = & 2~448~273.63450...
... E$\space \\
& $\pm 9$ & & $\pm 11$\space & \\
\end{tabular}\end{displaymath} (5)

Within the errors, separate least squares fits to all primary and secondary times yield identical periods. The period is close to the most recent published determination (Srivastava 1987); see also Kreiner[*]. Srivastava found evidence for a slowly increasing period and sinusoidal variations of the $\rm O{-}C$ residuals, indicating a possible third body in the system. The data in Table A.6 provide no clear evidence for such variations, however.

   
8 WZ Oph

WZ Oph = HD 154676 = HIP 83719 was discovered by Metcalf to be an eclipsing binary (Pickering 1917). Visual light curves and times of minima were published by Leiner (1921, 1926), who derived an orbital period of about $2 \hbox{$.\!\!^{\rm d}$ }09$, whereas the correct value of about $4 \hbox{$.\!\!^{\rm d}$ }18$ was proposed by McLaughlin (1929). Lause (1936) presented further visual minima. Sanford (1937) obtained the first radial velocities, and Gaposchkin (1938) established visual (from Leiner's observations) and photographic (from Harvard patrol plates) light curves, and he derived the first absolute dimensions of the two nearly identical components of WZ Oph.

The first modern analysis of WZ Oph, leading to masses and radii accurate to 3-4%, was done by Popper (1965) from photoelectric light curves and 20 Å/mm spectra. WZ Oph is included in his comprehensive review on stellar masses (Popper 1980), but not in the later review on accurate masses of radii of normal stars by Andersen (1991), where a more strict limit of 2% was applied. Very accurate light curves and radial velocities are needed in order to reach such accuracies (e.g. Andersen et al. 1980), and besides standard ${\it uvby}$ indices (Hilditch & Hill 1975) and additional times of minima (see Tables A.7, A.8), practically no new observational data for WZ Oph have been published since those by Popper. We now remedy this lack.

   
8.1 uvby light curves

Complete ${\it uvby}$ light curves containing 697 points in each colour were observed on 41 nights during six periods between May 1991 and March 1997 (JD 2 448 370-JD 2 450 535). HD 155193 and HD 154931 were selected as comparison stars, and they were found to be constant within the observational accuracy of a few mmag during our observations; see Tables 1-3 for further information. The eclipses have been covered on several nights and out-of-eclipse phases, except for a few gaps, at least twice.

As seen in Fig. 9, WZ Oph consists of two well-detached components of almost identical surface fluxes in a circular orbit, and from the analysis by Popper (1965) they are also known to be very similar with respect to masses and radii. The accuracy per point is about 0.004 mag ($\it vby$) and 0.006 mag (u), meaning that they are of significantly higher quality than those by Popper (1965). However, our photometric analyses (Clausen et al. 2008) reveal that throughout all phases the points scatter by 1-4 mmag more than this, highest in u. At a given phase, the observations from different seasons do not differ systematically. The extra scatter may be due to low-level surface activity, since weak emission features in the Ca II H and K lines are seen at the position of both components on high-resolution spectra of WZ Oph.

For the ephemeris adopted by Popper, the component eclipsed at phase 0.0 (his component 1) was found to be (formally) slightly more massive, slightly smaller, and appeared to have slightly weaker lines and to be less luminous than the other. The new accurate ${\it uvby}$ light curves reveal that, adopting Popper's ephemeris, the eclipse at phase 0.5 is in fact slightly deeper in all four colours, meaning that the component eclipsed here has a slightly higher surface flux than the other. Since analyses of new CfA radial velocities (Clausen et al. 2008) show that this component is also marginally the more massive, we have decided to interchange the minima compared to Popper's notation, so for the linear ephemeris presented below, the deeper (primary) eclipse occurs at phase 0.0.

   
8.2 Ephemeris

From the ${\it uvby}$ observations, two times of primary and two times of secondary minimum are available, see Tables A.7 and A.8, where the 104 times of minima we have been able to locate in the literature are also included. We note that those by Gaposchkin (1938) are based on Harvard patrol plates with usual exposure times of one hour.

We adopt the following linear ephemeris, which is based on all available times, assuming a circular orbit, and derived as described in Sect. 3.2:


 \begin{displaymath}
\begin{tabular}{r r c r r}
{\rm Min ~ I} = & 2~450~535.78331...
...E$\space \\
& $\pm 24$ & & $\pm 24$\space & \\
\end{tabular}\end{displaymath} (6)

It represents a significant improvement compared to previous determinations. Excluding the rather inaccurate times by Gaposchkin (1938), which also show large $\rm O{-}C$ residuals in Tables A.7 and A.8, does not change the result significantly. The Lafler & Kinman (1965) method applied to the ${\it uvby}$ observations yields $4 \hbox{$.\!\!^{\rm d}$ }1835070$. Finally, separate least squares fits to all primary and secondary times yield identical periods within errors.

Note that, as mentioned above, the new epoch is shifted by half a period compared to (most) previous ephemerides given for WZ Oph, including that by Kreiner[*].

   
9 Summary and conclusions

Complete ${\it uvby}$ light cuves and ${\it uvby}\beta $ standard indices have been obtained for the six main-sequence eclipsing binaries AD Boo, HW CMa, SW CMa, 636, VZ Hya, and WZ Oph. For 636 and HW CMa, the light curves are the first modern ones obtained, and for the other systems they are both more complete and more accurate than earlier data. Ephemerides have been determined from all available times of minima, and apsidal motion has been detected and quantified for SW CMa and 636.

For AD Boo, the new photometry indicates low-level surface activity, and spectroscopy supports that the cooler secondary component is the most likely candidate (Clausen et al. 2008). The new photometry and spectroscopy also suggests that WZ Oph is mildly active. The ${\it uvby}$ light curves of 636 clearly show variations due to activity at a moderate level, also confirmed by spectroscopy (Clausen et al., in prep.).

In parallel, new spectroscopic observations for accurate radial velocities and abundance determinations were obtained. The combined data yield accurate absolute dimensions and heavy-element abundances for all six systems and allow detailed tests of current stellar evolutionary models. Our analyses will be published in three separate papers (Clausen et al. 2008; the F-type systems AD Boo, VZ Hya, and WZ Oph), (Clausen et al., in prep.; the G-type system 636); (Torres et al., in prep.; the late A-type systems HW CMa and SW CMa).

Acknowledgements

R. Clement, O. Liu, and M. Viskum kindly observed some of the binaries on a few nights as part of their own programs at SAT. Excellent technical support was received from the staff of Copenhagen University and from the staff at ESO, La Silla. We thank J. M. Kreiner for providing complete lists of times of eclipse for the binaries, and G. Torres for refining the ephemeris calculations for AD Boo, VZ Hya, and WZ Oph. The CORAVEL and CfA groups kindly made unpublished radial velocity information and preliminary spectroscopic orbits available for the preparation of this paper.

The projects ``Structure and evolution of stars - new insight from eclipsing binaries and pulsating stars'' and ``Stars: Central engines of the evolution of the Universe'', carried out at Copenhagen University and Aarhus University, are supported by the Danish National Science Research Council. LPRV greatfully acknowledges partial support from the Brazilian agencies CNPq, CAPES and FAPEMIG. Furthermore, this investigation was supported by the Danish Board for Astronomical Research, the Carlsberg Foundation (to JA), and the Spanish Science Research Council.

The following internet-based resources were used in research for this paper: the NASA Astrophysics Data System; the SIMBAD database and the VizieR service operated by CDS, Strasbourg, France; the ar$\chi$iv scientific paper preprint service operated by Cornell University; the BBSAG and Orion bulletins (http://www.bav-astro.de); the VSOLJ Variable Star Bulletin (http://vsolj.cetus-net.org); the AASVO archive (http://www.aavso.org).

References

 

  
Online Material

Appendix A: Times of minima tables

Table A.1: Times of primary minima for AD Boo. O-C values are calculated for the ephemeris given in Eq. (1) assuming a circular orbit. Observing methods are: PG photographic; V visual; PE photoelectric or CCD. Times marked with an asterisk were published after we completed the ephemeris calculations; including them makes no significant changes.

Table A.2: Times of secondary minima for AD Boo. See Table A.1 for details.

Table A.3: Times of primary (P) minima of HW CMa. O-C values are calculated for the ephemeris given in Eq. (2).

Table A.4: Times of primary (P) and secondary (S) minima of SW CMa. O-C values are calculated for the apsidal motion ephemeris given in Table 4. Times given in brackets have not been used for the ephemeris determination.

Table A.5: Times of primary (P) and secondary (S) minima for V636 Cen. O-C values are calculated from the apsidal motion parameters given in Table 5.

Table A.6: Times of primary (P) and secondary (S) minima for VZ Hya. O-C values are calculated for the ephemeris given in Eq. (5) adopting a circular orbit. Observing methods are: PG photographic; V visual; PE photoelectric or CCD. Times marked with an asterisk were published after we completed the ephemeris calculations; including them makes no significant changes.

Table A.7: Times of primary minima for WZ Oph. O-C values are calculated for the ephemeris given in Eq. (6) assuming a circular orbit. Observing methods are: PG photographic; V visual; PE photoelectric or CCD. Times marked with an asterisk were published after we completed the ephemeris calculations; including them makes no significant changes.

Table A.8: Times of secondary minima for WZ Oph. See Table A.7 for details.



Copyright ESO 2008