Free Access
Issue
A&A
Volume 506, Number 3, November II 2009
Page(s) 1269 - 1275
Section Stellar structure and evolution
DOI https://doi.org/10.1051/0004-6361/200810112
Published online 03 September 2009

A&A 506, 1269-1275 (2009)

The massive Wolf-Rayet binary LSS 1964 (=WR 29)[*]

II. The V light curve[*]

R. C. Gamen1,2,[*] - E. Fernández-Lajús1,2,[*] - V. S. Niemela[*] - R. H. Barbá3,4,[*]

1 - Instituto de Astrofísica de La Plata, CONICET, Paseo del bosque s/n, B1900FWA, La Plata, Argentina
2 - Facultad de Ciencias Astronómicas y Geofísicas, Universidad Nacional de La Plata, Paseo del bosque s/n, B1900FWA, La Plata, Argentina
3 - Instituto de Ciencias Astronómicas de la Tierra y del Espacio, CONICET, Avda. España 1512 Sur, J5402DSP, San Juan, Argentina
4 - Departamento de Física, Universidad de La Serena, Benavente 980, La Serena, Chile

Received 2 May 2008 / Accepted 1 August 2009

Abstract
Context. WR 29 is a known WN7h+O double-lined binary system with a rather short period (3.164 days).
Aims. We search for light variations to determine the inclination of the system and thus the absolute masses of both components.
Methods. We observed photometrically the field of WR 29 between December, 2002, and February, 2006.
Results. We find that the V light of WR 29 varies in phase with the spectroscopic period of 3.16412 days, presenting two minima corresponding to the conjunctions of the binary components. Numerical models fitted to the light curve indicate an orbital inclination of about 44$^\circ$, and masses of 53 $M_\odot$ and 42 $M_\odot$ for the O- and WN-type components, respectively.

Key words: stars: binaries: close - stars: binaries: eclipsing - stars: binaries: spectroscopic - stars: Wolf-Rayet

1 Introduction

Stellar mass is one of the most important parameters of massive stars. Together with the mass-loss rate, chemical abundance, and rotation (Meynet & Maeder 2005), it defines how stars live and die. In spite of its importance, there are few determinations of masses in massive stars. A direct way to estimate stellar masses is by analysis of eclipsing and spectroscopic binary systems.

Among the 299 Galactic Wolf-Rayet stars listed in the VIIth Catalogue (van der Hucht 2001) and its annex (van der Hucht 2006), only 48 are known as spectroscopic binaries, and 26 are double-lined (SB2) binaries (Rauw et al. 2004; Gamen & Niemela 2002; Niemela et al. 2008; van der Hucht 2001; Gamen et al. 2008). Only 17 of them present photospheric or nonphotospheric eclipses, and in many cases the light curves appear poorly defined (see e.g. Lamontagne et al. 1996).

The binary nature of LSS 1964 (WR 29, van der Hucht et al. 1981) was revealed by Niemela & Gamen (2000, afterwards Paper I). They presented the first orbital analysis of both components of WR 29 based on low-resolution digital spectra, and obtained very high minimum masses (40 and 33 $M_\odot$ for the WN and O components, respectively). This led us to carry out photometric observations of WR 29 looking for light variations that could help to determine the inclination of the binary system and thus the absolute masses.

In this paper, we present for the first time the V light curve of the WN+O binary system WR 29. We performed a Wilson-Devinney (WD) model analysis of the V-band light curve and the radial velocity orbits, which include new radial velocities measured in low- and high- resolution spectra.

2 Observations and data reduction

2.1 Photometry

Digital CCD images were acquired with the 2.15-m J. Sahade telescope at Complejo Astronómico El Leoncito[*] and the 0.9-m telescope at Cerro Tololo Interamerican Observatory (CTIO), Chile, between 2002 and 2006 (See Table 1 for a summary of the observations).

Table 1:   Observing logs of photometric runs.

In the CASLEO observations, we used a focal reducer to image a wider field. The data were obtained through a V filter with two different CCD cameras. A camera with a Tektronics back-illuminated chip of $1024 \times 1024$ pixels (24 $\mu$m square pixel) was used between 2002, December, and 2004, February. Another camera with a Roper 1300B EEV CCD36-40 chip of $1340 \times 1300$ pixels (20 $\mu$m square pixels) was used during 2004-2005. The V filter used is that recommended by Bessell (1990). Our instrumental configurations with the two cameras result in samplings of 0 $.\!\!^{\prime\prime}$813 per pixel (Tek), and of 0 $.\!\!^{\prime\prime}$678 per pixel (Roper), respectively. In the majority of the images, the FWHM of the stellar images was between 1 $.\!\!^{\prime\prime}$9 and 3 $.\!\!^{\prime\prime}$1.

V-band images from the 0.9-m Cassegrain telescope[*] at CTIO were taken during 7 consecutive nights in February 2006. This telescope is equipped with a dedicated $2048\times2048$Tektronix CCD camera, with a $\sim$ $0\hbox{$.\!\!^{\prime\prime}$ }4$ per pixel scale and a whole field size of about $13\hbox{$.\mkern-4mu^\prime$ }5 \times 13\hbox{$.\mkern-4mu^\prime$ }5$. Typical seeing during these nights at CTIO was about $1\hbox{$.\!\!^{\prime\prime}$ }6$.

All of the CCD images were corrected for bias, flat-fielded with lamp-flats or twilight-flats, and reduced using aperture photometry with IRAF[*]. The radius of the aperture for all the stars was 7 $.\!\!^{\prime\prime}$32. The magnitudes of WR 29 in this aperture include the faint visual neighbor star to the South-East (See Fig. 1). We performed a PSF fit of both stars individually in one of our best frames (where the two stars are well separated) and obtained a difference of $3.75 \pm 0.009$ mag, which means that the visual companion contributes about 3% of the light intensity within the radius used for aperture photometry.

Differential photometry of WR 29 was determined using the star S1113013266, in the GSC2.2 catalog (STScI, 2001), as the comparison because its brightness is similar to that of WR 29 and is quite constant during the epochs of our observations. We checked this star with a fainter star, S111301326256 in GSC2.2, obtaining rms errors in differential magnitudes of about 0.009 mag, and the differences between the maximum and minimum values are less than 0.056 mag. We determine that the rms errors in the differential V-magnitudes of WR 29 are smaller than 0.02 mag. These data are presented in Table 2 (the complete dataset is available as an electronic file at the CDS) where in successive columns we quote the Heliocentric Julian Date of each observation, the differential magnitude of WR 29 (WR 29-C1), the orbital phase (from Eq. (1)), and the configuration used (Tek+2.15-m, CASLEO, Roper+2.15-m, CASLEO, or Tek+0.9-m, CTIO). A typical CASLEO frame of the field of WR 29 is shown in Fig. 2, where the comparison and check stars are labeled.

Table 2:   Differential photometry of WR 29.

\begin{figure}
\par\resizebox{5.0 cm}{!}{\includegraphics{10112fig1.eps}}
\end{figure} Figure 1:

A zoom of the vicinity of WR 29, indicating the radius of the photometric aperture including the close visual neighbor to the SE.

Open with DEXTER

\begin{figure}
\par\includegraphics[width=12cm,clip]{10112fig2.eps}
\end{figure} Figure 2:

A typical frame in the direction of WR 29 obtained at CASLEO in 2004, February. The comparison star S1113013266 is indicated as ``C1'' and the check star S111301326256 is labeled as ``C2''.

Open with DEXTER

2.2 Spectroscopy

We have obtained new spectral images of WR 29 between 2000 and 2006, aiming to refine the orbital solution published in Paper I.

Fifteen low-resolution spectra were collected at CASLEO, between 2000 and 2003, using the REOSC echelle spectrograph in its single dispersion mode, attached to the 2.15-m reflector. A Tek  $1024 \times 1024$ pixel CCD, with pixel size of 24 $\mu$m, was used as the detector. The reciprocal dispersion was $\sim$1.8 Å pixel-1, and the observed wavelength region was 3850-5450 Å.

We have also acquired 8 spectra with the echelle spectrograph attached to the 2.5-m du Pont telescope at LCO, during April, 2005. The detector used was a $2048\times2048$ format Tek#5 CCD with 24 $\mu$m pixels, providing simultaneous wavelength coverage from 3700 to 7000 Å with a resolving power of about 26000 (the FWHM of comparison-arc emission lines, with the binned CCD, is about 0.19 Å pixel-1 at 5000 Å).

Eleven low-resolution spectra were obtained with the R-C grating spectrograph attached to the 4-m telescope at CTIO, during 6 nights of March, 2007. We used the KPGL1 grating which provides a reciprocal dispersion of 1.0 Å per pixel, and a wavelength coverage of about 3650-6700 Å. The detector used is a Loral 3k CCD with 3k $\times$ 1k 15 $\mu$ pixels.

We measured the same lines with radial velocities (RV) published in Paper I. New measured RV along with their corresponding Heliocentric Julian date are shown in Table 3 for completeness.

Table 3:   Observed Heliocentric radial velocities (in km s-1) of WR 29.

3 Results and analysis

3.1 The V light curve

The differential V magnitude of WR 29 clearly presents variations with time. We phased the observation times with the already known spectroscopic period (Paper I) and revealed that these variations are orbit related. Later, we tried to determine the period using only these data, finding a period P=3.16410 days, as will be described in Sect. 3.3. The light curve is presented in Fig. 3, as a function of the orbital ephemeris given in Eq. (1). There, it is seen that the system presents two eclipses, i.e. when the WR is in front (phase $\phi =0.0$) and also at the other conjunction ($\phi=0.50$), when the disc of the O-type star eclipses the WN component.

The shape of the light curve is such that the beginning and end of each eclipse are not well-defined. The broad wings of the minima indicate that both stars have quite large radii and since this is a short period binary, its components could be in near contact and the system might be at least semi-detached. The eclipse minimum depths are clearly different, i.e. $\Delta_1 \sim 0.\!{\!^{\rm mag}}$16 and $\Delta_2 \sim 0.\!{\!^{\rm mag}}$12 for the primary and secondary eclipses respectively. Also, both light maxima are rather different; we estimate that the maximum at $\phi=0.25$ is $\sim$0.02 mag brighter than the other at around $\phi=0.75$. Although both minimum times are not completely sampled in our observations, the eclipses seem to be equidistantly separated, in concordance with the shape of the radial-velocity curves. The wings of the eclipses seem to be symmetrical. Some short time scale variability is present around the mean light curve at different phases, but especially during eclipses. We also found that the shape of the light curve shows variations for different epochs. We believe that this variability is real, and could be associated with phenomena due to circumstellar material and/or to random wind variability.

3.2 The radial velocity orbit

In order to improve the orbital solution already published in Paper I, we measured the same lines in our new spectra.

We ran a program to search for orbital elements, which is based on Bertiau & Grobben (1968), with the complete RV dataset, i.e. our new measurements plus those of Paper I. We took into account the different instrumental configurations and the spectral S/N involved in our dataset, assigning to the LCO echelle spectra a weighting of 1 and to the CASLEO lower-resolution spectra a weighting of 0.6. We inspected the individual $\rm O-C$ values derived by the program, but no systematic RV-shift was found. Thus, we did not apply any correction to the data. Leaving the period and eccentricity free, we obtained periods of 3.16404 $\pm$ 0.00002 days, 3.16414 $\pm$ 0.00006 days, 3.16410 $\pm$ 0.00002 days, and 3.16407 $\pm$ 0.00003 days for the N  IV $\lambda $4058 emission, N  IV $\lambda $5203 absorption, He  II $\lambda $4686 emission, and N  V $\lambda $4606 absorption line, respectively. We obtained rather low eccentricities in each dataset (between $0.03\pm 0.01$ and $0.10\pm 0.03$), but in view of the short period we prefer to be cautious and to consider a circular orbit, as the higher values were obtained in the solutions with more dispersion, i.e. obtained with the N  IV $\lambda $5203 and N  V $\lambda $4606 absorptions. Very similar values were obtained when we searched for periodicities imposing e=0. Looking at the similar periods found from the different lines, we adopted a value of 3.16404 days, because it is the period found using the best defined line in the spectra of WR 29, i.e. N  IV $\lambda $4058. We performed new determinations of orbital parameters fixing the eccentricity. These were similar to those published in Paper I, except for:

  • the semi-amplitude of the radial-velocity variations of the N  IV $\lambda $5203 absorption line. The new solution gives $K = 236\pm 8$ km s-1 (38 km s-1 smaller than in Paper I). We think this difference is due to the great scatter that exists at quadrature when this absorption line reaches its maximum negative velocity and blends with a solar line at $\lambda $ 5199, which is present in our oldest spectra obtained during nearly full moon observing runs (used in Paper I);
  • the new systemic velocity found for the N  V $\lambda $4606 absorption line is 27 km s-1 less negative than the published value;
  • the new systemic velocity found for the He  II $\lambda $4686 emission line is different by tens of km s-1. We think, as in the previous item, that these differences could be instrumental in origin but are within the errors, without consequences on the physics of the system;
  • the orbital variations of the absorption-line RVs from the O component are quite different from those of Paper I. We think the absorption lines measured in the lower-resolution spectra (obtained at CASLEO) may be contaminated by the P-Cygni profiles of the H$\delta $ (+N III) emission (see Fig. 4).
To definitively be sure of the O-component feature's existence, we performed an estimate on a dynamic spectrum constructed with the CTIO spectra (a specific wavelength range is depicted in Fig. 6). The most conspicuous absorption lines belonging to the O-type companion are (from strongest to faintest): H$\gamma $, He  I  $\lambda4471$, He  II $\lambda4200$, and H$\delta $. The last two are very poorly defined and, H$\delta $ could be contaminated by other profiles like Si  IV  $\lambda4089$ and $\lambda4116$ or N  III  $\lambda4097$ and $\lambda4107$ (from the WN-type component). Therefore, we only take into account the RVs of the H$\gamma $ absorption line. We also tried to extract the mean spectra of each of the two stars, by means of the technique of Demers et al. (2002). We note that individual spectra of the ``disentangled'' O-type component vary with orbital phase, i.e. the O features become stronger when the WN component is in front of the system (See Fig. 5), which could be related to absorption of the O-star's continuum at line frequencies by ions in the WR wind. This behaviour invalidates a reliable spectral classification using this technique.

\begin{figure}
\par\resizebox{8.5 cm}{!}{\includegraphics{10112fig3.eps}}
\end{figure} Figure 3:

The V light curve of WR 29 (in differential magnitudes). The primary minimum occurs when the WN star is at inferior conjuction.

Open with DEXTER
\begin{figure}
\resizebox{14.4 cm}{!}{\includegraphics[angle=270]{10112fig4.eps}}
\end{figure} Figure 4:

Comparison between low-resolution (in red or grey) and echelle spectra (black), in the wavelength domain of H$\delta $ and H$\gamma $, of WR 29 at two representative orbital phases, 0.15 and 0.75. Note how the blue-shifted absorption lines indicated by lines can be deblended in the echelle spectrum.

Open with DEXTER

As a first step to obtain values of the physical parameters of the components of the system WR 29, we determined a double-lined RV orbit for the binary. We used the RVs of the N  IV 4058 Å emission line (as representing the WN motion), and the H$\gamma $ absorption line (O motion). Although it is more appropiate to use the intrinsic absorption line N  IV $\lambda $5203 (see Paper I) as representing the WN motion, as its orbital solution (semi-amplitude) is very similar to that obtained using the N  IV 4058 Å emission line, we decided to use the latter because its solution presents smaller errors. We performed independent orbital fits to each dataset of RVs, i.e. N  IV emission and mean absorption, to characterize their respective systemic velocities. The results of both datasets were rather different, so a small correction of 17 km s-1 was made to the RVs of the H$\gamma $ absorption line (not introduced in Table 3). All the new orbital parameters adjusted for the circular orbits of the components of WR 29 are presented in Table 4, and they were used as constraints in the subsequent fit to Wilson-Deviney models.

As can be seen from Table 4, this new orbital solution implies that the O-type component is slightly more massive than the WN star.

\begin{figure}
\par\resizebox{9.5 cm}{!}{\includegraphics[angle=0]{10112fig5.eps}}
\end{figure} Figure 5:

Residual spectra (individual minus the mean WN7h) of WR 29 at both conjunctions, with the O star behind ($\phi =0.0$) and in front ($\phi =0.5$) of the system. The spectrum of the O-type star seems to vary with the orbital phase, i.e. the absorption features are stronger when the light of the O star passes through the expanded envelope of the WN star.

Open with DEXTER

Table 4:   Circular orbital solution for both component of WR 29.

3.3 The WD analysis of WR 29

In order to obtain fairly reliable absolute masses of each binary component, we have adjusted a numerical eclipsing binary model to the observations, using the Wilson-Devinney (WD) code (Wilson 1990; Wilson & Devinney 1971; Wilson & van Hamme 2004). We made use of the PHOEBE package tool (Prsa & Zwitter 2005) for the light curve (LC) fitting and differential correction (DC) fits. Our photometric data and RVs of both components of the binary system were used for this purpose.

The system was assumed to be an over-contact binary not in thermal contact. This mode fixes the $\Omega_{\rm WR}$ potential to make its value that of the lobe filling. The luminosity of the WR component is coupled to the temperatures. Some simple considerations were applied for the emission parameters of the stars in the model, i.e. the stars were considered as black bodies and the approximate reflection model (MREF=1) was adopted. A third light element was applied according to that mentioned in Sect. 2.1. No spots were included in this work, but they should be considered in future and refined works about WR 29, to take into account the gas that probably exists in between both stars, and could be producing the many asymmetries seen in the light curve of WR 29, like different maximum brightness or irregularities in the eclipse wings. Gravity darkening exponents g1 = g2 = 1 and bolometric albedos $\rm Alb_{1} = Alb_{2} = 1$ were set for radiative envelopes. Limb-darkening was obtained with a square-root law. Bolometric limb darkening coefficients for the V band wavelengths were extrapolated from van Hamme (1993). The temperature for the O-type component T1 was fixed from the Spectral Type- $T_{\rm eff}$ calibration tables published by Martins et al. (2005), because we see some He II absorption lines belonging to the O component. Since we are not sure of the subtype (or luminosity class) of the O-type component, we adopted a conservative mean $T_{\rm O}=45~000$ K and adjusted the temperature of the WN component $T_{\rm WR}$ but we are conscious that only the temperature ratio has astrophysical relevance. We considered that the system has a circular orbit and that both components rotate synchronously ( F1 = F2 = 1). The RV of the center of mass of the binary system was adopted from the RV orbit. Weights were applied to photometric points taking into account the internal dispersion within defined phase-bins. Spectroscopic points were weighted according to the resolving power and S/N of spectra.

\begin{figure}
\par\resizebox{16 cm}{!}{\includegraphics[angle=-90]{10112fig6.eps}}
\end{figure} Figure 6:

Dynamic grey-scaled spectra of WR 29 (from black to white) constructed with the 11 spectra obtained at CTIO. The most important features are indicated. Note how the ``S-shape'' of the He I $\lambda $4471 absorption line is mounted on a wider emission and moves in anti-phase with respect to other emission lines (the S-shapes of H$\gamma $, H$\delta $, and He II $\lambda $4200 are more difficult to see).

Open with DEXTER

We adopted the spectroscopic period, and obtained a reasonable but simple model fit for the light curve (LC). Such a model does not account for the systematic deviations from the mean LC. This would require adding discrete mechanisms such as star spots. We do not do this, however, since this would be far from our principal aim, i.e. to determine absolute stellar masses.

We used the DC code of WD complemented by the Levenberg-Marquardt algorithm for differential corrections within PHOEBE, to get an integrated solution for the whole dataset. Then, we fitted P, T0, a, q, $V_{\gamma}$, i, $T_{\rm WR}$, and $R_{\rm O}$. The ephemeris of the adopted solution is:

\begin{displaymath}{\rm Min} \mbox{ }I = \rm JD~ 2453039.87 + 3.\!{\!^d}16412 \cdot {\it E} .\\
\end{displaymath} (1)

The inclination of the orbital plane of the system relative to the sky plane was found to be $44.\!{\!^\circ}0 \pm 1.\!{\!^\circ}5$. The obtained orbital size is $41.4\pm1$ $R_{\odot}$ and the effective radii, $R_{\rm O}=18.4\pm 0.5$ and $R_{\rm WR}=14.9\pm0.5$ $R_{\odot}$. The Wolf-Rayet star is 1.2 times cooler than the O-type component. For the adopted $T_{\rm O}$ a rather lower $T_{\rm WR}=37~000\pm1000^\circ$K is derived. However, we empirically verified that the temperatures do not influence the stellar radius determinations. The absolute stellar masses are $M_{\rm O}=53.3\pm3.5$ and $M_{\rm WR}=42.1\pm3.5$ solar masses for the O-type and the WN-type, respectively. The result of our best fitting model is presented in Fig. 7 together with the photometric points of the LC of WR 29, and the measured RV of both components. Although preliminary, in Table 5 we show some of the relevant parameters of the WD solution.

Table 5:   Preliminary absolute parameters of WR 29.

\begin{figure}
\center
\resizebox{9 cm}{!}{\includegraphics[angle=0]{10112fig7.eps}}
\end{figure} Figure 7:

Simultaneous WD fittings to the V light curve (in relative intensities) and RV orbits of WR 29. The lower panel shows the RV curves obtained from the model, and the RV variations of the N IV 4058 Å emission line from the WN star and the Balmer H$\gamma $ line from the O star. Full squares depict the RVs of the absorption line measured in the few echelle spectra. Note the apparent phase-shift that seems to exist between phases 0.75 and 1.25 in the orbit of the N  IV emission line (see text for explanation).

Open with DEXTER

4 Concluding remarks

We have revealed that the V-band light of the binary system WR 29 varies in a periodic way, showing partial eclipses at both conjunctions of the system. We performed a simple WD analysis constraining some fundamental parameters of the spectroscopic orbital elements (Table 4). We should emphasize here that this is a very simplified model which is unable to consider stellar winds. However, we think it can be used as a rough initial estimation of the masses and radii of the stellar components of WR 29, as applied to the analysis of other WR+O systems such as WR 20a (Bonanos et al. 2004), and WR 155 (Leung et al. 1983). From our best fits we found WR 29 to be a contact binary, with a moderate percentage of over-contact, at least, of the O-type component. Indeed, both stars lie very close to each other, but retain quite different temperatures ( $T_{\rm O}/T_{\rm WR} \sim 1.2$).

Within the uncertainties, the first values of physical parameters of the O-type component of WR 29, which we have determined, are in fair agreement with the tabulations of a Galactic O5 I star[*] based on models of stellar atmospheres (cf. Martins et al. 2005). We intended to fit many models with smaller radius, but they did not reproduce the observed wider primary minimum. Thus, the O-type component appears to be in a supergiant stage.

The RVs of the N  IV $\lambda4058$ emission line (and also the other measured lines) show a phase-shift with respect to the model between orbital phases $\phi=0.75{-}1.25$. If true, this can be due to asymmetries in the line forming region caused by temperature differences in the WN winds (on its side facing the O star, etc.), or by the presence of the shock cone (the O star wind replaces portions of the WR wind). Similar phase-shifts are also seen in other WR+O binary systems, namely WR 21, WR 31, WR 47, and WR 97 (Gamen 2004).

We note the similarity of the light curve of WR 29 with that of CQ Cep ($\equiv$WR 155), a WN6+O9II-Ib (Smith et al. 1996) binary with a period $P\sim1.64$d (See Figs. 6 and 7 in Leung et al. 1983). Moreover, both these binary systems share even more features in common, e.g. RVs measured in emission lines present phase-shifts, and the O-type companion seems to be filling its Roche lobe. The derived radius of the WN component in CQ Cep is $R \sim 11$ $R_{\odot}$ (that in WR 29 is 4 $R_{\odot}$ larger). This could indicate that the radius is limited by the smaller size of the binary system. The most intriguing difference is that the WN component in CQ Cep is more massive than the O-type secondary (Harries & Hilditch 1997). Then, WR 29 could be the future evolutionary stage of a system like CQ Cep.

We derived absolute masses of 53 $M_\odot$ and 42 $M_\odot$ for the O-type and the WN7h-type components, respectively. Although these masses are preliminary, that for the WN star agrees with that derived for the other WN7h-type star in a SB2 system, WR 22, whose WN component is $55\pm7 ~M_\odot$ (Schweickhardt et al. 1999).

The hypothesis that two stars of $\sim$ $50 ~M_\odot$ should be related to some known cluster was explored. The closest (in angle) young clusters to WR 29 are Bochum 11 and Sher 1 (as found in the WEBDA database, Mermilliod 1995). The former is located at 35 arcmin from WR 29 and its earliest population is a compoundof an O4 III star and at least 7 OB-type stars, its derived age is about 5 Myr but its distance is 2.4 kpc (Smith 2006), which seems be too low for WR 29. The other cluster, Sher 1, whose center is 1.32 degrees from WR 29, is also a young cluster (age = 107) and its derived distance is 10 kpc (Moffat et al. 1991); then the projected linear distance to WR 29 is about 230 pc which seems be too large. A tailored study to explore the presence of a hitherto unknown cluster or association in the field of WR 29 is needed. In a circle of 1 degree around WR 29 another four WR stars exist (i.e. numbers 30a, 31a, 31b, 31c) with no relationship to any cluster, too. WR 31a (Hen 3-519) and WR 31b (AG Car) are known luminous blue variables located 23$^\prime$ and 40$^\prime$ from WR 29, respectively.

Acknowledgements
We thank the referee, Anthony Moffat, for his suggestions and comments which contributed to considerably improve this paper. We thank the Director and staff of CASLEO, LCO, and CTIO for the use of their facilities, and acknowledge the use, at CASLEO, of the CCD and data acquisition system supported under US NSF grant AST-90-15827 to R.M. Rich. R. B. acknowledges partial support from Universidad de La Serena, Project DIULS CD08102. This research has made use of the WEBDA database, operated at the Institute for Astronomy of the University of Vienna.

References

  • Bertiau, F., & Grobben, J. 1968, Ric. Astr. Spec. Vat., 8, 1
  • Bessell, M. S. 1990, PASP, 102, 1181 [NASA ADS] [CrossRef]
  • Bonanos, A. Z., Stanek, K. Z., Udalski, A., et al. 2004, ApJ, 611, L33 [NASA ADS] [CrossRef]
  • Demers, H., Moffat, A. F. J., Marchenko, S. V., Gayley, K. G., & Morel, T. 2002, ApJ, 577, 409 [NASA ADS] [CrossRef]
  • Gamen, R. 2004, Ph.D. Thesis, La Plata University
  • Gamen, R., Gosset, E., Morrell, N., et al. 2008, Rev. Mex. Astron. Astrofis. Conf. Ser., in press
  • Gamen, R. C., & Niemela, V. S. 2002, New Astron., 7, 511 [NASA ADS] [CrossRef]
  • Harries, T. J., & Hilditch, R. W. 1997, MNRAS, 291, 544 [NASA ADS]
  • Lamontagne, R., Moffat, A. F. J., Drissen, L., Robert, C., & Matthews, J. M. 1996, AJ, 112, 2227 [NASA ADS] [CrossRef]
  • Leung, K.-C., Seggewiss, W., & Moffat, A. F. J. 1983, ApJ, 265, 961 [NASA ADS] [CrossRef]
  • Martins, F., Schaerer, D., & Hillier, D. J. 2005, A&A, 436, 1049 [NASA ADS] [CrossRef] [EDP Sciences]
  • Mermilliod, J.-C. 1995, in Information & On-Line Data in Astronomy, ed. D. Egret, & M. A. Albrecht, Ap&SS Library, 203, 127
  • Meynet, G., & Maeder, A. 2005, in The Nature and Evolution of Disks Around Hot Stars, ed. R. Ignace, & K. G. Gayley, ASP Conf. Ser., 337, 15
  • Moffat, A. F. J., Shara, M. M., & Potter, M. 1991, AJ, 102, 642 [NASA ADS] [CrossRef]
  • Niemela, V., Gamen, R., Barbá, R., et al. 2008, MNRAS, submitted
  • Niemela, V. S., & Gamen, R. 2000, A&A, 362, 973 [NASA ADS]
  • Prsa, A., & Zwitter, T. 2005, ApJ, 628, 426 [NASA ADS] [CrossRef]
  • Rauw, G., De Becker, M., Nazé, Y., et al. 2004, A&A, 420, L9 [NASA ADS] [CrossRef] [EDP Sciences]
  • Schweickhardt, J., Schmutz, W., Stahl, O., Szeifert, T., & Wolf, B. 1999, A&A, 347, 127 [NASA ADS]
  • Smith, L. F., Shara, M. M., & Moffat, A. F. J. 1996, MNRAS, 281, 163 [NASA ADS]
  • Smith, N. 2006, MNRAS, 367, 763 [NASA ADS] [CrossRef]
  • van der Hucht, K., Conti, P., Lundstrom, I., & Stenholm, B. 1981, Space Sci. Rev., 28, 227 [NASA ADS]
  • van der Hucht, K. A. 2001, New Astron. Rev., 45, 135 [NASA ADS] [CrossRef]
  • van der Hucht, K. A. 2006, A&A, 458, 453 [NASA ADS] [CrossRef] [EDP Sciences]
  • van Hamme, W. 1993, AJ, 106, 2096 [NASA ADS] [CrossRef]
  • Wilson, R. E. 1990, ApJ, 356, 613 [NASA ADS] [CrossRef]
  • Wilson, R. E., & Devinney, E. J. 1971, ApJ, 166, 605 [NASA ADS] [CrossRef]
  • Wilson, R. E., & van Hamme, W. 2004, in Computing Binary Star Observables, available at ftp://ftp.astro.ufl.edu/pub/wilson/lcdc2003/
 

Footnotes

... (=WR 29)[*]
Full Table 2 is only available in electronic form at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/506/1269
... curve[*]
Based on data collected at CASLEO, CTIO, and LCO.
...[*]
Visiting Astronomer, CASLEO, Argentina, and CTIO and LCO, Chile.
...[*]
Visiting Astronomer, CASLEO, Argentina.
... Niemela[*]
In memoriam (1936-2006).
...[*]
Visiting Astronomer, LCO and CTIO, Chile.
... Leoncito[*]
CASLEO is operated under agreement between CONICET, SeCyT, and the Universities of La Plata, Córdoba and San Juan, Argentina.
... telescope[*]
Operated by the SMARTS consortium.
... IRAF[*]
IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Association of Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science Foundation.
... star[*]
Martins et al's values for O5 I: R=19.4 $R_{\odot}$; $M_{\rm spec}=50.7$ $M_\odot$.
Copyright ESO 2009

All Tables

Table 1:   Observing logs of photometric runs.

Table 2:   Differential photometry of WR 29.

Table 3:   Observed Heliocentric radial velocities (in km s-1) of WR 29.

Table 4:   Circular orbital solution for both component of WR 29.

Table 5:   Preliminary absolute parameters of WR 29.

All Figures

  \begin{figure}
\par\resizebox{5.0 cm}{!}{\includegraphics{10112fig1.eps}}
\end{figure} Figure 1:

A zoom of the vicinity of WR 29, indicating the radius of the photometric aperture including the close visual neighbor to the SE.

Open with DEXTER
In the text

  \begin{figure}
\par\includegraphics[width=12cm,clip]{10112fig2.eps}
\end{figure} Figure 2:

A typical frame in the direction of WR 29 obtained at CASLEO in 2004, February. The comparison star S1113013266 is indicated as ``C1'' and the check star S111301326256 is labeled as ``C2''.

Open with DEXTER
In the text

  \begin{figure}
\par\resizebox{8.5 cm}{!}{\includegraphics{10112fig3.eps}}
\end{figure} Figure 3:

The V light curve of WR 29 (in differential magnitudes). The primary minimum occurs when the WN star is at inferior conjuction.

Open with DEXTER
In the text

  \begin{figure}
\resizebox{14.4 cm}{!}{\includegraphics[angle=270]{10112fig4.eps}}
\end{figure} Figure 4:

Comparison between low-resolution (in red or grey) and echelle spectra (black), in the wavelength domain of H$\delta $ and H$\gamma $, of WR 29 at two representative orbital phases, 0.15 and 0.75. Note how the blue-shifted absorption lines indicated by lines can be deblended in the echelle spectrum.

Open with DEXTER
In the text

  \begin{figure}
\par\resizebox{9.5 cm}{!}{\includegraphics[angle=0]{10112fig5.eps}}
\end{figure} Figure 5:

Residual spectra (individual minus the mean WN7h) of WR 29 at both conjunctions, with the O star behind ($\phi =0.0$) and in front ($\phi =0.5$) of the system. The spectrum of the O-type star seems to vary with the orbital phase, i.e. the absorption features are stronger when the light of the O star passes through the expanded envelope of the WN star.

Open with DEXTER
In the text

  \begin{figure}
\par\resizebox{16 cm}{!}{\includegraphics[angle=-90]{10112fig6.eps}}
\end{figure} Figure 6:

Dynamic grey-scaled spectra of WR 29 (from black to white) constructed with the 11 spectra obtained at CTIO. The most important features are indicated. Note how the ``S-shape'' of the He I $\lambda $4471 absorption line is mounted on a wider emission and moves in anti-phase with respect to other emission lines (the S-shapes of H$\gamma $, H$\delta $, and He II $\lambda $4200 are more difficult to see).

Open with DEXTER
In the text

  \begin{figure}
\center
\resizebox{9 cm}{!}{\includegraphics[angle=0]{10112fig7.eps}}
\end{figure} Figure 7:

Simultaneous WD fittings to the V light curve (in relative intensities) and RV orbits of WR 29. The lower panel shows the RV curves obtained from the model, and the RV variations of the N IV 4058 Å emission line from the WN star and the Balmer H$\gamma $ line from the O star. Full squares depict the RVs of the absorption line measured in the few echelle spectra. Note the apparent phase-shift that seems to exist between phases 0.75 and 1.25 in the orbit of the N  IV emission line (see text for explanation).

Open with DEXTER
In the text


Copyright ESO 2009

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.