A&A 405, 1087-1093 (2003)
K. Yakut 1 - A. E. Tarasov 2,3 - C. Ibanoglu1 - P. Harmanec4,5 - B. Kalomeni1 - D. E. Holmgren 6,7 - H. Bozic8 - P. Eenens9
1 - University of Ege, Faculty of Science, Department of Astronomy and Space Sciences, 35100 Bornova, Izmir, Turkey
2 - Crimean Astrophysical Observatory, 98409 Nauchny, Crimea, Ukraine
3 - Isaac Newton Institute of Chile, Crimean Branch, 98409 Nauchny, Crimea, Ukraine
4 - Astronomical Institute of the Charles University, V Holesovickách 2, 180 00 Praha 8, Czech Republic
5 - Astronomical Institute of the Academy of Sciences, 251 65 Ondrejov, Czech Republic
6 - SMART Technologies, Inc., Bay 8, 1460 - 28 Street NE, Calgary, Alberta, T2A 7W6, Canada
7 - Dominion Astrophysical Observatory, Canada
8 - Hvar Observatory, Faculty of Geodesy, Zagreb University, Kaciceva 26, 10000 Zagreb, Croatia
9 - Departamento de Astronomia, Universidad de Guanajuato, Apartado 144, 36000 Guanajuato, GTO, Mexico
Received 28 March 2003 / Accepted 6 May 2003
New light and radial-velocity curves of V497 Cep, a binary in the open cluster NGC 7160, were obtained and the linear ephemeris of the system was refined to . The first light and radial-velocity curve solutions allowed us to derive the basic physical properties of this astrophysically important binary. It was found that the observed light variation of V497 Cep consists of a strong ellipticity effect and a small contribution from grazing eclipses. A comparison of masses and radii of V497 Cep with theoretical evolutionary tracks indicates that both binary components are very close to the zero-age main sequence. A comparison of disentangled line profiles of the He I 6678 line with synthetic, rotationally broadened line profiles indicates that the rotation of both stars is synchronized with the orbital revolution as expected. This finding increases the credibility of our solutions. We find . The distance to the cluster NGC 7160 was found to be about 760 (100) pc which agrees well with other available estimates.
Key words: stars: binaries: eclipsing - stars: binaries: spectroscopic - stars: individual: V497 Cep - Galaxy: open clusters and associations: individual: NGC 7160
BD+612213 was first reported to be a light variable by Hill (1967a,b). He tentatively classified it as an eclipsing variable, not on the basis of his photometry but on the basis of an earlier report of its radial-velocity (RV hereafter) variability by Hayford (1932). Hill et al. (1976) published individual photometric observations of the star, again indicating mild light changes.
Overlooking these reports, Harmanec et al. (1994) used the star occasionally as one of the UBV standards in the course of their systematic observations of Be and Ap stars and some binary systems at Hvar and Skalnaté Pleso observatories. Alerted by larger systematic deviations of this "standard" during the 1998 UBV observations, Harmanec et al. (1999) homogenized and re-analyzed all available photometric observations, including a rich set of 107 good Hipparcos observations, and concluded that the star varies periodically, with a period of 1 2028251. Harmanec et al. suggested that the star is more probably an ellipsoidal variable than an eclipsing binary. They appealed to spectroscopists and photometrists to collect new accurate observations which would allow to improve the knowledge of the basic physical elements of the binary and, consequently, also of the distance to cluster NGC 7160. Subsequently, Kazarovets et al. (2000) assigned the object a variable-star name V497 Cep.
This prompted us to include V497 Cep into our spectroscopic and photometric observing programs at several observatories. Preliminary results of an analysis of spectral observations secured at Crimean Astrophysical Observatory and Dominion Astrophysical Observatory (DAO hereafter) were published in a poster paper by Holmgren & Tarasov (2000). They arrived at RV semiamplitudes of 165 and 228 km s-1 , for the primary and secondary, respectively.
The aim of this paper is to present and analyze new light and RV curves of the system and to derive for the first time also the absolute dimensions of the system.
New systematic BV observations of V497 Cep were carried out at Ege University Observatory. The 0.48-m Cassegrain reflector equipped with a SSP-5 photometer and Johnson B and V filters was used. The comparison star already used by Harmanec et al. (1999), HD 208440 = NGC 7160*3, was used. Its all-sky UBV magnitudes derived by Harmanec et al. were adopted. V497 Cep was observed on 10 nights between July 31, 1999 and October 17, 2000. All differential observations were corrected for atmospheric extinction, using the extinction coefficient derived for each night of observation from the observations of the comparison star. A total of 270 observations in B, and 290 in V band were obtained.
Additionally, all UBV observations homogenized and used by Harmanec et al. (1999), together with a few new observations secured at Hvar and San Pedro Mártir observatories since then, were also used. Altogether, we accumulated 492 V, 365 B and 95 U observations, covering the period HJD 2441499.8-2451843.8. For convenience of future investigators, we publish all these observations in Table 1.
Initially, the phase of each observation was calculated using
the light elements given by Harmanec et al. (1999):
|Figure 1: The B and V light curves of V497 Cep, based on Ege University observations.|
|Open with DEXTER|
|Figure 2: Selected He I 6678 profiles of V497 Cep. The profiles are identified by phase according to ephemeris (1). Longer tick marks on the flux scale are separated for 0.2 of relative flux units.|
|Open with DEXTER|
|Figure 3: Selected H 6563 profiles of V497 Cep. The profiles are identified by phase according to ephemeris (1). Longer tick marks on the flux scale are separated for 0.2 of relative flux units.|
|Open with DEXTER|
Crimean observations were secured at two spectral regions: near H and near He I 6678 Å. Only three H spectra were obtained, in June and July 1997. In the He I 6678 Å region, 13 spectra were accumulated: 9 of them in October-November 1997 and 4 in July 1998. The detector was a CDS 9000 CCD array with pixels. All the observations were obtained in the second order of a diffraction grating with a reciprocal dispersion of 3 'A mm-1. The spectral resolution was about 30 000 and each spectrum covers about 60 Å. The exposure time ranged from 30 to 60 min, depending on the weather conditions. All spectra were reduced using the standard techniques: dark-current subtraction and flat-field division. The subsequent reduction included the following procedures: removal of cosmic-ray spikes; sky background subtraction; normalization to the continuum approximated by a third-degree polynomial; and wavelength calibration using a ThAr comparison lamp, with an accuracy better than 1.5 km s-1. Wavelength scale and calculated Julian dates (JDs hereafter) were corrected to the barycenter of solar system.
First 17 DAO spectra were obtained in July 1999 using the 1.2-m telescope, 3261H spectrograph, and the UBC 4096 CCD with pixels. The wavelength region was centred on H, and each spectrum covered approximately 1200 Å. The reciprocal dispersion was approximately 20 'A mm-1. An exposure time of about 1800 s was used. Four additional spectra were secured in the year 2000. The data were reduced in the usual fashion using the MIDAS package. The bias and flat-field frames were averaged and used to remove the bias level and to remove the detector response from all stellar exposures. Each stellar frame was also cleaned of cosmic-ray spikes. One-dimensional spectra were then extracted and wavelength-calibrated. Finally, all stellar spectra were rectified using a spline interpolation. All wavelength scales were adjusted for heliocentric radial velocity.
In all analyses, the Crimean spectra were given 5 times higher weight than the DAO spectra to take their much better linear dispersion into account. All spectra were imported into reduction program SPEFO (Horn et al. 1996) and the radial velocities were determined from the comparison of direct and mirror image of each line profile. Heavily blended profiles near conjunctions were not measured. Following Horn et al. (1996) we also measured a selection of reliable telluric lines for all DAO and Crimean H spectra and used the differences of the mean telluric velocity and the calculated heliocentric correction to bring the data onto a joint RV zero point. Applying these corrections, we created input data for KOREL disentangling program (Hadrava 1995,1997) and derived orbital solution and disentangled line profiles for the He 6678 line which is more suitable for the mass determination than the intrinsically broader H line. For comparison, we also used program FOTEL (Hadrava 1990) for the orbital solution based on classical RV measurements.
Individual radial velocities derived by direct measurements and by
KOREL are given in detail in Table 2. Orbital elements based on a
FOTEL solution using classical RV measurements and those derived
via KOREL are summarized in Table 3 and the
radial-velocity curves of the He I 6678 line derived via direct
measurements and FOTEL solution, and via KOREL disentangling are
compared in Fig. 4. Note that KOREL can use even blended
line profiles which accounts for a higher number of data points in
the bottom panel.
|Figure 4: A comparison of radial-velocity curves of the primary and secondary measured directly in SPEFO (upper panel) and based on KOREL disentangling (bottom panel). Filled symbols refer to the primary and the Crimean spectra are shown by circles, and the DAO spectra by triangles. The final ephemeris (2) was used.|
|Open with DEXTER|
He I 6678 line profiles of both binary components, disentangled by KOREL, are shown in detail in Fig. 5.
We give two FOTEL solutions in Table 3, allowing for a difference in the zero point of RV scale in the first case (the region near the He I 6678 line is free of telluric lines and no check on the zero point is therefore possible there). One can see that the classical and KOREL solutions agree within the limits of quoted errors. However, the KOREL solution gives a much smaller rms error and should be preferred for the determination of the masses of both binary components.
Radial velocities of the H line give a similar solution but with slightly reduced amplitudes which is a known effect. There are too few spectra as yet to obtain a complete removal of the telluric spectrum near H with KOREL. We therefore postpone the study of the H profiles for a future study.
We arrived at the following improved ephemeris
Table 3: Radial-velocity solutions for the He I 6678 line, based on classical RV measurements (solutions 1 and 2) and on KOREL disentangling (solution 3). Period was kept fixed at 1 2028287 and zero eccentricity was assumed. All epochs are in HJD-2451022 and the quoted errors of the elements are rms errors. The rms error of the solution is the rms error of one point of unit weight.
|Figure 5: He I 6678 line profiles of primary and secondary disentangled via KOREL (solid lines) are compared with synthetic LTE line profiles shifted for the systemic velocity of -19 km s-1 and rotationally broadened to projected rotational velocities of 132 and 104 km s-1, for the primary and secondary, respectively, which correspond to spin-orbit synchronization. Wavelengths are in Å and the disentangled profiles from KOREL were normalized using the luminosity ratio from the V light-curve solution - see the text for details.|
|Open with DEXTER|
The observed character of new light curves does not correspond to what would be expected for an ellipsoidal variable, advocated by Harmanec et al. (1999).
The solution of the new light curves were obtained using the latest version of the Wilson-Devinney code (Wilson & Devinney 1971; Wilson 1994). For the main-sequence stars of spectral types B0 to A0, one can estimate the interstellar reddening, using the reddening-free quantity Q=(U-B) - 0.72 (B-V). According to the standard observed UBV colours of V497 Cep obtained by Harmanec et al. (1999) we get Q=-0.604 which corresponds to a B3V star (Drilling & Landolt 2000). According to the mass of the primary component its effective temperature was assumed to be 19 500 K, following Drilling & Landolt (2000), and it was kept fixed during the light-curve solution. The limb darkening coefficients of x1=0.441 and x2=0.388(Claret 2000), albedos A1 = A2 = 1.0 (Rucinski 1969) and gravity exponents of g1 = g2 = 1.00 (von Zeipel 1924) were assumed and kept fixed. The brightness of V497 Cep varies continuously in all orbital phases. A model light curve of a pure reflection effect has always an almost sinusoidal shape. Different choices of gravity- and limb-darkening coefficients led to only minute modifications. The deviations from a pure sinusoidal shape increase slightly with increasing orbital inclination in the sense that maximum become narrower. However, the observed light curve has a maximum which is almost equally broad as the minimum. Therefore, a satisfactory description of the light curve can only be achieved if one allows for some ellipticity, acting together with the reflection.
In order to match the observed light curve, the model requires a trade off between several parameters. First of all, the inclination of the orbit must be such to allow grazing eclipses to occur. The secondary component should have a similar radius and effective temperature as the primary, given the observed similarity in the depths of both eclipses. Having these constraints in mind, we used mode 2 of the Wilson-Devinney code (suitable for detached binaries with no constraints on surface potentials) to find the elements leading to the optimal agreement between model light curve and observations. A circular orbit and a black-body radiation were assumed for both stars. The following elements were adjusted during the solution: orbital inclination i, temperature of the secondary T2, dimensionless surface potentials and , and the relative luminosities of the components in respective passbands, L1and L2. The mass ratio of 0.783 is well determined by our RV solution 3. It was, therefore, kept fixed in the light-curve solution. Since the V light curve has unequal maxima, only the B light-curve solution was used to obtain the geometrical and physical elements of the system. The functional dependencies of the effective temperature on mass and spectral type differ significantly from one calibration to another, especially for the massive stars. The results of the light curve analysis will of course depend on the adopted temperature of the primary component. Various effective temperatures for the primary component from 17 000 to 23 000 K were taken for the analysis of the light curve. The best fit, minimum squared residuals, was obtained at 19 500 K, which corresponds to a B3V star appropriate to the Q value. The best-fit values obtained after several dozen of iterations with the DC code are listed in Table 4. The uncertainties of the parameters are formal errors which are probably unrealistically low. The WD code tends to give small formal errors of the derived parameters even for data with a larger scatter. The theoretical and observed B light curves are compared in Fig. 6. We carried out, however, solutions of both, Band V light curves and obtained the fractional luminosities of the primary component as 0.654 and 0.651 in B and V, respectively. In spite of all effort, we were unable to obtain a reliable solution for the limited number of more noisier U-passband observations and postpone that task until new, more accurate photoelectric observations will become available.
Our solution indicates that the orbital inclination is about and only grazing partial eclipses occur. Our model explains the observed light curve as being due to the combined effect of partial eclipses and ellipsoidal variations in a system consisting of two early B type components with similar temperatures.
Table 4: The B light-curve solution carried out with WD code.
|Figure 6: A comparison of the observed (dots) and computed (solid line) B light curves of V497 Cep.|
|Open with DEXTER|
Table 5: Absolute parameters for V497 Cep. Errors of the last two digits are given in parentheses.It should be stated here that the uncertainties of the absolute parameters are underestimated and mostly depended on the fractional parameters derived light curve analysis. In Fig. 7, the positions of binary components are compared with theoretical evolutionary tracks in the HR diagram and radius vs. effective temperature diagram based on Schaller's et al. (1992) evolutionary models for the main sequence stars. One can see that the masses of both binary components agree very satisfactorily with their positions in both diagrams. It is also clear that both components of V497 Cep are little evolved main-sequence stars, with the primary being more evolved from the zero-age main sequence than the secondary as expected. The inherent remaining uncertainties of our solution does not allow us to use this comparison to an accurate determination of the age of V497 Cep. Yet, it is clear that our result does not contradict the evolutionary age of NGC 7160 cluster, yr, derived by Harris & Harris (2000).
|In the text||Figure 7: The positions of the stellar components of the eclipsing, double-lined binary V497 Cep in the a) log - log and b) log - log diagrams. Masses of corresponding evolutionary models are indicated near the corresponding points with both, zero-age (ZAMS) and terminal-age (TAMS) main sequence locations.|
|Open with DEXTER|
The line of He I 6678 is known to be prone to NLTE effects, therefore the less perfect agreement in the line core of the primary between the observed and synthetic profile is not very disturbing. Note also that we had to use a synthetic profile corresponding to somewhat higher effective temperature than what we believe is the true effective temperature of the primary. For the secondary, the agreement is nearly perfect.
In any case, one can see that the widths of the observed and model line profiles agree perfectly for both stars. This implies that both binary components indeed rotate with angular velocities synchronized with orbital revolution as expected for such a close binary. In turn, it represents another consistency check of our orbital and light-curve solutions.
Using the fractional luminosities of the components from
the B and V light-curve solutions, quoted above,
and adopting the observed V and B magnitudes
of the system at maximum light,
standard dereddening, which leads to
the following dereddened values for the primary and secondary:
PH gratefully acknowledges the warm hospitality during his 2002 visit to Department of Astronomy and Space Sciences of Ege University where some aspects of this paper were discussed. DH would like to thank the director of the Dominion Astrophysical Observatory for generous allocation of observing time. This study has been supported by Ege University Research Fund (Project No. 99/FEN/016 and 2002/FEN/002). It is also partly based on the observational data secured by PH and PE in the course of the Czech-Mexican KONTAKT/CONACyT collaborative project E130.734. In the early stages, it was carried out as a part of the research supported by the grants 205/96/0162 and 205/02/0445 of the Grant Agency of the Czech Republic. In later stages, PH was also supported via research plans J13/98: 113200004 of MSMT, and AV 0Z1 003909 and via project K2043105 of AV CR. Finally, we express our thanks to the referee, Dr. C. Lazaro, for his helpful, constructive criticism which helped to improve the paper.