Free Access
Issue
A&A
Volume 503, Number 3, September I 2009
Page(s) 873 - 881
Section Stellar structure and evolution
DOI https://doi.org/10.1051/0004-6361/200811514
Published online 29 April 2009

The triple system HIP 96515: a low-mass eclipsing binary with a DB white dwarf companion[*]

N. Huélamo1 - L. P. R. Vaz2 - C. A. O. Torres3 - P. Bergeron4 - C. H. F. Melo5 - G. R. Quast3 - D. Barrado y Navascués1 - M. F. Sterzik6 - G. Chauvin7 - H. Bouy8,[*] - N. R. Landin2

1 - LAEX-CAB (INTA-CSIC); Postal address: LAEFF, PO Box 78, 28691 Villanueva de la Cañada, Madrid, Spain
2 - Depto. de Física, Universidade Federal de Minas Gerais, C.P.702, 30161-970 - Belo Horizonte, MG, Brazil
3 - Laboratório Nacional de Astrofísica/MCT, Rua Estados Unidos 154, 37504-364 Itajubá, Brazil
4 - Département de Physique, Université de Montréal, C.P. 6128, Succ. Centre-Ville, Montréal, Québec H3C 3J7, Canada
5 - European Southern Observatory, Karl-Schwarzschild-Strasse 2, 85748 Garching bei Muenchen, Germany
6 - European Southern Observatory, Alonso de Cordova 3107, Casilla 19, Santiago, Chile
7 - Laboratoire d'Astrophysique, Observatoire de Grenoble, BP 53, 38041 Grenoble, Cedex 9, France
8 - Instituto de Astrofísica de Canarias, C Vía Láctea, s/n, E38205 - La Laguna (Tenerife), Spain

Received 11 December 2008 / Accepted 27 January 2009

Abstract
Context. HIP 96515 A is a double-lined spectroscopic binary included in the SACY catalog as a potential young star. It has a visual companion (CCDM 19371-5134 B, HIP 96515 B) at 8 $.\!\!^{\prime\prime}$6. If bound to the primary, the optical and infrared colors of this wide companion are consistent with those of a white dwarf.
Aims. We attempt to characterize the system HIP 96515 A&B by studying each of its components.
Methods. We analyzed spectroscopic and photometric observations of HIP 96515 A and its visual companion, HIP 96515 B. To confirm the system as a common proper-motion pair, we analyzed the astrometry of the components using high-angular resolution infrared observations obtained within a time span of two years, and archival astrometry.
Results. The high-resolution optical spectrum of HIP 96515 A was used to derive a mass ratio, M2/M1, close to 0.9. The optical lightcurve of HIP 96515 A shows periodic variations with $P_{\rm orbital}$ = 2.3456 days, revealing that HIP 96515 A is an eclipsing binary with preliminary orbital parameters of i = $89\hbox{$.\!\!^\circ$ }0$ $\pm $ $0\hbox{$.\!\!^\circ$ }2$, and M1 = 0.59 $\pm $ 0.03 $M_{\odot }$ and M2 = 0.54 $\pm $ 0.03 $M_{\odot }$, for the primary and secondary, respectively, at an estimated distance of 42 $\pm $ 3 pc. This is a new eclipsing binary with component masses below 0.6 $M_{\odot }$. Multi-epoch observations of HIP 96515 A&B show that the system is a common proper-motion pair. The optical spectrum of HIP 96515 B is consistent with a pure helium atmosphere (DB) white dwarf. The comparison with evolutionary cooling sequence models provides $T_{{\rm eff,WD}}$ = 19 126 $\pm $ 195 K, $\log g_{{\rm WD}}$ = 8.08, $M_{{\rm WD}}/M_{\odot}$ = 0.6, and a distance of $\sim $46 pc. The estimated WD cooling age is $\sim $100 Myr and the total age of the object (including the main-sequence phase) is $\sim $400 Myr. Finally, if HIP 96515 A&B are coeval, and assuming a common age of $\sim $400 Myr, the comparison of the masses of the eclipsing binary members with evolutionary tracks shows that they are underestimated by $\sim $15% and 10%, for the primary and secondary, respectively.

Key words: stars: white dwarfs - stars: binaries: eclipsing - stars: individual: HIP 96515

1 Introduction

In recent years, several infrared surveys have focused on the detection of substellar companions around young and nearby stars using high-contrast imaging techniques, (e.g., Lafreniere et al. 2007; Lowrance et al. 2005; Chauvin et al. 2003; Masciadri et al. 2005). Young nearby stars are privileged targets for these kind of studies due to their proximity (below 100 pc), which allows to study very small separations, and their youth, which enhances the contrast between the central star and the potential companion.

In 2006, we started an infrared survey to detect substellar objects around young stars from the SACY catalog (Torres et al. 2006,2003). Our sample included HIP 96515, a nearby (d=43.9 $\pm $ 9 pc, van Leeuwen 2007) M1-type star classified as a double-lined spectroscopic binary (SB2) by Torres et al. (2006). The object is catalogued as a visual binary (CR23, Worley & Douglass 1997), with a secondary companion at an angular separation of $\rho$ $\sim $ 8 $.\!\!^{\prime\prime}$6 (CCDM 19371-5134 B or HIP 96515 B Dommanget & Nys 2000). This visual binary has been detected in different infrared and optical surveys, although its properties have never been studied in detail.

As a first step in unveiling the nature of HIP 96515 B, we compared our $K_{\rm s}$-band photometry with evolutionary tracks (Baraffe et al. 2003; Chabrier et al. 2000). The comparison suggested that HIP 96515 B may be of substellar nature if the age of the primary is in the range 0.07-1 Gyr and the two objects are comoving. However, the optical and infrared colors retrieved from archival observations are inconsistent with being a substellar object, but consistent with those of a white dwarf.

Motivated by this preliminary result, we decided to perform a detailed study of HIP 96515. We analyzed new and archival observations of the triple system and derived the main properties of each source. Section 2 describes the main observational properties of HIP 96515 A&B. In Sect. 3, we analyze spectroscopic and photometric data of the primary star HIP 96515 A, while the true nature of the visual component, HIP 96515 B is unveiled in Sect. 4. A discussion of the system is provided in Sect. 5, and the main results are summarized in Sect. 6.

Table 1:   Main properties of HIP 96515 A.

Table 2:   Photometric data of the visual binary HIP 96515 A&B.

2 Observational properties of HIP 96515 A&B

The main stellar properties of HIP 96515 A are displayed in Table 1, while its optical and infrared photometry are shown in Table 2.

The analysis of the optical spectrum of HIP 96515 A obtained during the SACY survey, allowed Torres et al. (2006) to classify the object as an SB2 with almost equal-mass components (see Table 3). From the cross-correlation function (CFF), we estimated a magnitude difference of $\Delta~V$ = 0.55 mag between the binary components.

HIP 96515 was included in the ROSAT All-Sky Bright Star catalog (Voges et al. 1999). The source was detected by ROSAT with a count rate of 0.172 $\pm $ 0.036 counts s-1. We derived the X-ray flux using the conversion factor definition by Schmitt et al. (1995). The X-ray luminosity, assuming a distance of 44 pc (see Table 1), is $L_{{\rm x}}$ = (3.61 $\pm $ 0.75) $\times $ 1029 erg s-1. Using the bolometric correction (BC) by Kenyon & Hartmann (1995) for an M1-type star, and the V-mag listed in Table 2, we derived a $\log(L_{\rm x}/L_{\rm bol})$ ratio of -3.07. If we assume that the two components of the spectroscopic binary (HIP 96515 Aa & Ab) emit in X-rays, such that Lx,Aa=Lx,Ab=Lx/2, and using the V-mag of each of the binary members (see Table 3) and the BC for a M1 and M2 star, we derive $\log(L_{\rm x}/L_{\rm bol})$ ratios of -3.17 and -3.04 for Aa and Ab, respectively. These values are comparable to those found in very active pre-main sequence (PMS) late-type stars (e.g., Stelzer & Neuhäuser 2001; Barrado y Navascués et al. 1999) and in main-sequence (MS) late-type stars in very close binaries (see e.g. López-Morales 2007).

The age of HIP 96515 A is unknown. The lithium content is normally used as a youth indicator in the case of late-type stars. However, it is not always present in young M-type stars since it is depleted too fast (e.g., Basri et al. 1996; Torres et al. 2008a; Neuhäuser 1997). In fact, none of the binary members of HIP 96515 A show the lithium absorption doublet ($\lambda$ = 6707.8 Å) in the optical spectrum. Torres et al. (2008a) studied the lithium content as a function of the V-I color in different young associations with ages in the range 5-70 Myr. By comparing HIP 96515 A with similar M-type stars (with V-I $\sim $ 1.8, Torres et al. 2006) in those associations, we can compute a lower limit to the age of $\ge$70 Myr. By following the temporal evolution of lithium in stellar models (e.g., those by Landin et al. 2006) for stars with the masses of HIP 96515 A (Sect. 3), Li should be exhausted by an age between 30 Myr (for a mixing length parameter of $\alpha$ = 2) and 50 Myr (for $\alpha=1$).

Finally, HIP 96515 B is included in the AC catalog (Urban et al. 1998) as a visual companion to HIP 96515 A. The object was detected in different infrared and optical surveys and multi-wavelength photometry is available in different public catalogs. We included all the available data in Table 2. We note that, in most cases, the photometric measurements show large uncertainties.

3 Characterization of HIP 96515 A: a M1+M2 eclipsing binary

To characterize the SB2 in more detail, we retrieved photometric data of HIP 96515 A from both the Hipparcos archive (Perryman et al. 1997, 96 data points) and the All Sky Automated Survey, ASAS[*], (Pojmanski & Maciejewski 2005, 818 data points). The Hipparcos observations, in its own photometric system (Hip), were corrected with the relation $V_{\rm J} = Hip - 0.08$, suggested for the V-I of this object (V-I = 1.85; Torres et al. 2006).

The observations of the source are displayed in Fig. 1, showing significant variability. The object was classified as an eclipsing binary by Paczynski et al. (2006), using ASAS data alone. We applied the method by Lafler & Kinman (1965) to search for any periodicity in the data. The most likely period, 2.345 days, was used to phase-fold both the ASAS and Hipparcos data, confirming the SB2 as an eclipsing system. We now describe our full analysis of the photometric data and derive the main properties of the eclipsing components, Aa and Ab.

Table 3:   Spectroscopic data of the eclipsing binary HIP 96515 A.

 \begin{figure}
\par\includegraphics[width=9cm,clip]{1514fig1.eps}
\end{figure} Figure 1:

Hipparcos and ASAS V-band light-curves of HIP 96515 A.

Open with DEXTER

3.1 Photometric data: period search

The method by Lafler & Kinman (1965) was used to search for periods in the Hipparcos V-band observations of HIP 96515 A. We searched for periods between 0.5 day and 63 days, by using a relative period step of $\delta P/P$ = 4 $\times $ 10-7. Since M stars may flare, and HIP 96515 is a strong X-ray emitter, we excluded the 4 brightest measurements (from Hipparcos measurements) in this period search, to avoid clearly non-periodic events, if present. Approximately 60 significant periods were found, but the most significant was $P_{\rm orb}$ = (2.344451 $\pm $ 0.000002) d. In the analysis of the photometry with a version of the WD model (see Sect. 3.2), we adjusted a phase shift to the folded observations, and obtained the ephemeris (HJD; the digits in parentheses affect the last digits)

\begin{displaymath}%
{\rm Min~I} = 2~448~493.80515(10) + 2.3444507(15) E,
\end{displaymath} (1)

which we use in our analysis. According to Eq. (1), the secondary minimum occurs close to phase 0.5, indicating an orbit with a small (if any) eccentricity, and wihout further evidence, we assume in the preliminary analysis a circular orbit. We note that circular orbits are consistent with short orbital (and rotational) periods, as found for HIP 96515 A (e.g., Zahn 1989). The V-band lightcurves (Hipparcos and ASAS), folded with these ephemeris, are shown in Fig. 2. The stellar eclipses are clearly seen in both datasets.

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

The phase-folded lightcurves of HIP 96515 A from ASAS (+ symbols) and Hipparcos ($\times $ symbols), consisting of 818 and 96 observations, respectively, with our theoretical solution for q=0.9 obtained in Sect. 3.2. The Hipparcos observations were shifted by 0.6 mag, in order to separate the light curves. Note that ASAS observations were not corrected from the long term variation of the maxima, shown in Fig. 1.

Open with DEXTER

3.2 Preliminary analysis of the HIP 96515 A light-curve

We analyzed the Hipparcos optical lightcurve with a version of the WD model (Wilson & Devinney 1971; Wilson 1993,1979) extensively improved as described in Vaz et al. (2007) and references therein. Our version uses stellar atmospheres and models, the radiated flux of both components using the PHOENIX atmosphere models (Allard & Hauschildt 1995; Hauschildt et al. 1997a; Allard et al. 1997; Hauschildt et al. 1997b). The (linear-law) limb-darkening coefficients for both components were taken from Claret (2000), and interpolated using a bilinear scheme for the current values of $\log g$ and  $T_{\rm eff}$ at each iteration. The reflection albedos for both components were fixed to a value of 0.5, which is appropriate for stars with convective envelopes. The gravity-brightening exponent, $\beta$, was computed using the local value of  $T_{\rm eff}$ for each point on the stellar surface and taking into account mutual illumination following Alencar & Vaz (1997) and Alencar et al. (1999). Our procedure combines the WD modified code with several UNIX scripts and auxiliary programs to guarantee the consistency of the solutions in all steps.

We have a single measurement of the radial velocity for each component at the phase 0.7715, according to Eq. (1), by chance, close to the quadrature (0.75 in a circular orbit), when the line separation is maximum. Since we do not have any estimation of the center of mass velocity of the system, the mass ratio is undefined. Therefore, to investigate the solution space, we decided to compile a grid of solutions for different values of the mass ratio, q(=M2/M1), varying from 1.000 to 0.850 in steps of 0.025. In compiling this grid, we can determine the velocity of the center of mass of the system for each q, and adjust the remaining parameters to reproduce the individual heliocentric velocities measured at phase 0.7715.

We assigned an arbitrary initial value of $89\hbox{$^\circ$ }$ to the system inclination and used the WD model to find initial values of the gravitational pseudo-potentials (which define the size and form of the components) to reproduce the optical lightcurve. Finally, the initial values of the effective temperatures were adopted from the tabulated values listed in Popper (1980), and $T_{{\rm eff,pri}}$ was kept fixed at the adopted value of 3714 $\pm $ 150 K. In the lack of further information, the eclipsing components were considered to be synchronized with the orbital motion, an assumption that may not be true (see below).

Despite the presence of a visual companion to HIP 96515 A at  $8\hbox{$.\!\!^{\prime\prime}$ }6$, we do not expect any contamination in the Hipparcos lightcurve of the primary, since HIP 96515 B is outside the Hipparcos field of view. On the other hand, the resolution of the ASAS setup is 14 $.\!\!^{\prime\prime}$2/pixel and the light from HD 96515 B is definitely included in their measurements. Therefore, we assume a small amount of third light contribution in the analysis of the ASAS data lightcurve. From the appearance of the folded light curve (using our ephemeris) and in the absence of a well determined radial velocity curve, we assumed a circular orbit in our analysis.

Taking into account the afore mentioned initial values and assumptions, we modeled the orbital inclination, the gravitational pseudo-potentials, the secondary effective temperature, a phase shift in the primary minimum, and the primary internal luminosity of each point in the grid. The most interesting results from our modeling are displayed in Table 4.

Table 4:   Grid of solutions for 7 values of the mass ratio, q(=M2/M1), for HIP 96515 A.

After convergence of the solutions in Table 4, we adjusted the amount of third light in the q=0.900 solution, because of evidence that the visual companion, HIP 96515 B, contaminates the ASAS measurements. In units of the eclipsing system light at quadrature, the third light parameter was determined to be ${\cal L}_3$ = 0.005 $\pm $ 0.038, and the corresponding best fit solution was found to be almost indistinguishable from the solution with  ${\cal L}_3 = 0$.

To reproduce the observed radial velocities, we adjusted the semi-major orbital axis which, together with the orbital inclination and period, fixed the theoretical radial velocity amplitudes, K1 and K2 (in km s-1). The apparent V-magnitude retrieved from SIMBAD is 11.10 mag. However, the Hipparcos measurement converted to Johnson's V is equal to 11.73 mag, and our measurement is 11.78 mag (Torres et al. 2006), the value adopted for the distance estimation. By using bolometric corrections for both components from Popper (1980), we can derive the absolute dimensions and, finally, estimate the distance to be (42 $\pm $ 3) pc. The interstellar absorption in such a short distance should be negligible. The parallax for the system measured by Hipparcos is (22.79 $\pm $ 4.87) mas, which infers (44 $\pm $ 9) pc, in good agreement with our determination.

All solutions of Table 4 are almost indistinguishable from the curve-fitting point of view, as can be seen by the rms of the residuals ($\sigma$, 1 obs) in magnitudes. The probable intrinsic variability of one or both components makes it necessary to follow the system photometrically closely in time, preferably by observations in different passbands, and simultaneously with spectroscopic data to obtain radial velocity measurements. Only by obtaining more precise lightcurves in different photometric bands and a precise radial velocity curve will it be possible to eliminate this degeneracy in the solutions for the mass ratio. In this paper, we assume that the components are on the main-sequence, and make use of the mass-luminosity relation to fix q.

Figure 2 shows our solution for q=0.9 and Fig. 3 shows the corresponding residuals O-C. The steady and slow decreasing maxima found in ASAS observations was not taken into account in this analysis. It is common for M dwarfs to exhibit variations (e.g., Rockenfeller et al. 2006), mainly due to spots on their surfaces. However, we are unable to determine for the data that we analyzed, if the variability is caused by one or both eclipsing components. We searched for periodicities in the O-C residuals between 0.08 day and 50 days, using the phase-dispersion method by Lafler & Kinman (1965) and relative period steps of $\delta P/P$ = 3.d-3, and found a significant period of (0.24965 $\pm $ 0.00001) day (about 6 h), as can be seen in Fig. 4. We note that the variation pattern is present in both the Hipparcos and the ASAS data sets. The cadence of these observations (one or two observations a day, long intervals without observations being common) is inadequate for detecting these systematic variations in the O-C, which has the unfortunate result of introducing a multitude of aliases in the analysis. However, it is clear that the period of 0.24965 day, if real, is not the only one: Fig. 1 infers longer period variations, and we estimate in the O-C at least another possible period of (4100 $\pm $ 200) days, but of low significance level due to the relatively short duration of the observations ($\sim $6630 days).

 \begin{figure}
\par\includegraphics[width=9cm,clip]{1514fig3.eps}
\end{figure} Figure 3:

The residuals O-C between the V light curves of HIP 96515 A from Hipparcos ($\times $ symbols) and ASAS (+ symbols) shown in Fig. 2, and our theoretical solution for q = 0.9 obtained in Sect. 3.2. The standard error (1 observation) is shown in the upper left part of each panel.

Open with DEXTER

 \begin{figure}
\par\includegraphics[width=9cm,clip]{1514fig4.eps}
\end{figure} Figure 4:

The residuals O-C between the V light curves of HIP 96515 A from ASAS (+ symbols) and Hipparcos (grey circles) shown in Fig. 3, folded with the period of (0.24965 $\pm $ 0.00001) day. Note that the pattern is present in both Hipparcos and ASAS data sets (two cycles are plotted).

Open with DEXTER

The period of 0.24965 day detected may be either due to pulsation or spots on the surface of one or both of the eclipsing components. Rockenfeller et al. (2006) listed other M dwarfs of 6 h period. However, if the variation is due to spots on HIP 96515 A, it would indicate that the spotted component(s) spins significantly (more than a factor of 9) faster that the orbital motion (see e.g., CV Boo Torres et al. 2008b). This is probably not the case in HIP 96515 A: using an orbital period of 2.345 days, and the radius of the binary components from Table 5, we estimate rotational velocities of $\sim $14 km s-1 and $\sim $12 km s-1 for Aa and Ab, respectively. These values are comparable with those derived from the optical spectrum (Table 3, and assuming an inclination of 89 degrees), indicating that the system is synchronized. Therefore, the origin of the 6h-period remains unclear with the current dataset.

Table 5:   Preliminary physical parameters for HIP 96515 Aa&Ab.

By combining the results of Table 4, and taking into account the maximum variations in the variables along the grid, we can derive the preliminary absolute dimensions for HIP 96515 A, which are listed in Table 5. The masses of the binary components, HIP 96515 Aa and Ab, are below 0.6 $M_{\odot }$. Hence, HIP 96515 A is the fourteenth confirmed EB with component masses below 0.7 $M_{\odot }$ (see e.g., Shkolnik et al. 2008). All these EBs are important when calibrating theoretical evolutionary models.

Finally, we note that the eclipsing system definitely deserves more observations, both multi-band photometry with a cadence that allows the determination of precise light curves, and spectroscopy, to obtain reliable radial velocity curves and spectral types for both components. If possible, HIP 96515 B, should also be observed simultaneously, to ensure consistency in the analysis.

4 Characterization of HIP 96515 B

We now focus on the characterization of HIP 96515 B. First, we investigated if the object forms a common proper-motion pair with HIP 96515 A by analyzing diffraction-limited infrared observations of the pair obtained in the course of our program to detect substellar objects around SACY targets (Huelamo et al. 2008, in prep.). As a second step, we obtained an optical spectrum to derive the main properties of the object.

4.1 NACO/VLT observations of HIP 96515

We observed HIP 96515 twice with NAOS-CONICA (NACO), the Adaptive Optics facility at the Very Large Telescope (VLT), during the nights of the 2006 May 25 and 2008 June 16. The target was observed with the visible wavefront sensor and the $K_{\rm s}$-band filter. We used the S27 objective (nominal plate scale of $\sim $0.027 $^{\prime\prime}$/pixel), which provides a total field-of-view of $\sim $27 $^{\prime\prime}$ $\times $ 27 $^{\prime\prime}$. The total exposure time was $\sim $12 min on source.

The data were obtained using a random jitter between the exposures to compute the sky emission. The images were reduced with the Eclipse reduction package (Devillard 1997) following the standard steps: dark-subtraction, flat-field division, and image shifting and stacking. The final image of HIP 96515 is displayed in Fig. 5. Apart from HIP 96515 A, we detected only the already known visual companion, HIP 96515 B, at a projected separation of 8 $.\!\!^{\prime\prime}$6 ($\sim $378 AU at 44 pc) from the primary. We note that the eclipsing binary is too tight to be resolved in the NACO images.

The difference in magnitude between HIP 96515 A&B is 6.98 $\pm $ 0.07 mag. Taking into account the brightness of the primary ($K_{\rm s}$ = 7.9 mag, see Table 2), we estimate a  $K_{\rm s}$-mag of 14.94 $\pm $ 0.08 for the secondary.

The NACO images were used to derive the separation and position angle of the two objects at two different epochs. To derive accurate astrometry of HIP 96515 B in the two images, we calibrated the plate scale and orientation of the infrared detector, CONICA, with the astrometric calibrator $\theta^1$ Ori C (McCaughrean & Stauffer 1994). The results are displayed in Table 6.

4.2 HIP 96515 A & B: a common proper-motion pair

The astrometric accuracy of the NACO/VLT data is high enough to investigate if HIP 96515 A&B are comoving. For completeness, we compared the NACO/VLT data with older (and less accurate) astrometric measurements from public surveys (e.g., Dommanget & Nys 2000; Cutri et al. 2003). The visual binary is included in the Washington Double Star Catalog (WDSC, Worley & Douglass 1997) with three observations between 1910 and 2000. Table 6 includes the angular separation and position angle of the two components at five different epochs. As can be seen, both remain almost constant during a time interval of 98 years.

 \begin{figure}
\par\includegraphics[width=9cm,clip]{1514fig5.ps}
\end{figure} Figure 5:

VLT/NACO $K_{\rm s}$-band image of HIP 96515 A and its visual companion, HIP 96515 B, obtained in June 2008. The projected separation is of 378 AU at a distance of 44 pc. Note that the eclipsing binary (Aa & Ab) is not spatially resolved by NACO.

Open with DEXTER

 \begin{figure}
\par\includegraphics[width=7.5cm,clip]{1514fig6.eps}
\end{figure} Figure 6:

Differences in right ascention and declination between HIP 96515 A & B at different epochs (black circles). The expected position of a background object (grey circles and labels) at each epoch is also displayed. Top: all the astrometric data from 1910 to 2008. Bottom: only the NACO data, which provides the highest astrometric accuracy. The figure shows that the two objects are comoving, that is, they are a common proper motion pair.

Open with DEXTER

Table 6:   Relative astrometry of HIP 96515 A & B.

Figure 6 (top panel) shows the difference in right ascension (RA) and declination (Dec) of the binary members measured at five epochs between 1910 and 2008. We have also overplotted the expected difference in RA and Dec of a background object taking into account the proper motion of the primary. If we focus on the NACO/VLT data only (Fig. 6, bottom panel), which has the highest astrometric accuracy, we can confirm that HIP 96515 A & B are comoving, that is, they form a common proper-motion pair.

4.3 The optical spectrum of HIP 96515 B: observations and modeling

To ascertain the nature of HIP 96515 B, we obtained an optical spectrum of the target with the ESO Multi-Mode Instrument (EMMI) on the New Technology Telescope (NTT) on October 2007. The wavelength range of the spectrum was between 380 nm and 700 nm, and the spectral resolution was $R \sim 1100$. The data were reduced with standard routines within IRAF, including bias subtraction and flat-fielding, and flux-calibrated using a spectrophotometric standard observed during the night.

 \begin{figure}
\par\includegraphics[width=8.5cm,clip]{1514fig7.ps}
\end{figure} Figure 7:

Top: NTT/EMMI Optical spectrum of HIP 96515 B. We have marked all the identified He I lines. Bottom: fit of the optical spectrum using pure helium model atmospheres. The best fit parameters are displayed at the top right corner.

Open with DEXTER

The optical spectrum of HIP 96515 B is displayed in Fig. 7. The continuum increases towards bluer wavelengths and the most obvious detected spectral features are the He  I absorption lines, which are typical of helium (DB) white dwarfs.

We compared the optical spectrum of HIP 96515 B with atmospheric models of DB white dwarfs, to derive the physical parameters of the object. Our grid of model atmospheres and synthetic spectra for DB stars is described in Beauchamp et al. (1996). These include the improved Stark profiles of neutral helium from Beauchamp et al. (1997). The model atmospheres that we used here assume a pure helium composition and cover a range of effective temperatures of $T_{\rm eff}=10~000~(1000)$, 16 000 (2000), 30 000 K, and surface gravities of $\log g=7.0~(0.5)~9.0$. Our fitting technique relies on the non-linear least squares method of Levenberg-Marquardt (Press et al. 1986), which is based on a steepest descent method. The model spectra (convolved with a Gaussian instrumental profile) and the optical spectrum were first normalized to a continuum defined to be unity. The calculation of $\chi^2$ was then completed for these normalized line profiles only. We adopt a pure helium composition based on the lack of any H$\alpha$ absorption feature; additional models with mixed helium and hydrogen compositions allow us to set a limit of  $N({\rm H})/N({\rm
He}) <10^{-5}$.

As usual, there is a cool and a hot solution on each side of the temperature at which the line strengths reach their maximum, but the hot solution close to $T_{\rm eff}$ = 25 000 K can easily be excluded based on the slope of the observed energy distribution. Our final solution, $T_{\rm eff}$ = 19 130 $\pm $ 190 K and $\log g$ = 8.08 $\pm $ 0.07, is displayed in Fig. 7. These atmospheric parameters can be converted into a mass of M = 0.65 $\pm $ 0.04 $M_{\odot }$ using evolutionary models with C/O cores, $q({\rm He})\equiv \log M_{\rm He}/M_{\star}$ = 10-2, and $q({\rm H})$ = 10-10, which are representative of helium-atmosphere white dwarfs[*].

By using the best-fit parameters of the model, we can derive the absolute magnitudes in different photometric bands (see Table 7). Using $M_{\rm J}$ and the 2MASS J-band brightness of HIP 96515 B, which is the measurement with the smallest photometric error, we estimate a distance to the source of d = 45.8 pc, which is in good agreement with both the Hipparcos parallax and the distance derived from the analysis of HIP 96515 A.

Table 7:   Best-fit parameters for HIP 96515 B.

The evolutionary sequences described above, as well as empirical relations (see e.g., Fontaine et al. 2001), allowed us to estimate the mass of the white dwarf progenitor and its cooling age. These data are provided in Table 7. We obtained two different estimates of the white dwarf progenitor mass. First, we used the empirical initial-final mass relationship (IFMR) from Wood (1992), which is based on the study of the white dwarf luminosity function and mass distribution (see also Fontaine et al. 2001; Leggett et al. 1998). This yielded a value of $M_{\rm MS}=3.8$ $\pm $ $0.5~M_{\odot}$. Alternatively, we used the IFMR determined by Ferrario et al. (2005) based on white dwarfs observed in open clusters (see in particular their Fig. 1). This second relation yielded a significantly smaller final mass of $M_{\rm MS}=2.5$ $\pm $ $0.5~M_{\odot}$. However, as discussed by Ferrario et al., the IMFR of most clusters appears to have an intrinsic spread (i.e., there is a range of initial masses that gives rise to a given final mass). Hence, a value of $M_{\rm MS}$ as high as 3.8 $M_{\odot }$ could not be necessarily ruled out. In the following, we thus adopt a mean value of $M_{\rm MS}=3.1$ $\pm $ $0.8~M_{\odot}$ to allow for the possible mass range of the progenitor.

We derived the main-sequence lifetime of HIP 96515 B using two sets of evolutionary tracks for a 3.1 $\pm $ 0.8 $M_{\odot }$ star. First, we used the ATON 2.4 models, which are a modified version of models by D'Antona & Mazzitelli (D'Antona & Mazzitelli 1994,1997), and are fully described in Landin et al. (2006). According to these evolutionary models, a 3.1 $\pm $ 0.8 $M_{\odot }$ star spends 280+346-125 Myr on the main-sequence, i.e., from the zero-age main-sequence (ZAMS) to the terminal age main-sequence (TAMS). On the other hand, Claret (2004) models (Z=0.02) predict a span of $\sim $317 Myr for the MS time of a 3.1 $M_{\odot }$. Assuming a MS lifetime of $\sim $300 Myr and a cooling time of $\sim $100 Myr (see Table 7), the total age of the system would be $\sim $400 Myr, where the errors are dominated by the uncertainty in the mass.

 \begin{figure}
\par\includegraphics[width=9cm,clip]{1514fig8.eps}
\end{figure} Figure 8:

Stellar luminosity ( top panel) and surface gravity ( bottom panel) versus effective temperature. HIP 96515 Aa and Ab are represented by filled circles. The grey and black solid lines represent evolutionary tracks by Baraffe et al. (1998) and Landin et al. (2006), respectively, for 0.6 and 0.5 $M_{\odot }$ stellar masses. The dotted lines represent the uncertainties in the stellar masses for the ATON models. The comparison of HIP 96515 A with the two sets of evolutionary tracks provides an age of $\sim $60 Myr, that is, places the eclipsing binary members on the pre-main sequence, although with large uncertainties.

Open with DEXTER

5 Discussion

The comparison of HIP 96515 Aa & Ab with evolutionary tracks by Baraffe et al. (1998) provides an age of $\sim $50 Myr for masses of 0.5 $M_{\odot }$ (with alpha = 1), and places the binary on the PMS (see Fig. 8). The PMS nature of HIP 96515 A is also supported by comparison with tracks generated by the ATON 2.4 stellar evolutionary code, where we adopted the solar chemical composition (Z=0.0175 and Y=0.27) and alpha equal to 1. The tracks in Fig. 8 generated with the ATON2.4 code cover the PMS, MS and earlier stages of post main-sequence phases. The most probable age, according to ATON 2.4, is $\sim $60 Myr (Fig. 8), although with large uncertainties.

On the other hand, our astrometric study shows that the white dwarf and the EB are co-moving. If we assume that HIP 96515 A&B are coeval, the age of the triple system should be $\sim $400 Myr, which is significantly older than the age provided by the PMS tracks for HIP 96515 A, although within the observational uncertainties of the primary star (see Fig. 8).

To shed light on the evolutionary status of HIP 96515 A, we have studied the kinematical properties of the object. Using an average distance of 45 pc (from the EB and WD independent studies presented here), the RV derived by Torres et al. (2006), and the proper motions measured by Hipparcos (see Table 1), we have derived the UVW components of the Galactic space velocity vector: UVW = -6.0, -0.5, -21.6 km s-1. We have compared these values with the UVW vectors of young moving groups in the solar neighborhood included in Zuckerman & Song (2004) and Torres et al. (2008a). We conclude that HIP 96515 A does not share the kinematical properties of any of these associations.

Assuming an age of 400 Myr for the triple system, a possible explanation of the age of 60 Myr derived from PMS evolutionary tracks is the different evolution of both EBs and very magnetically active stars, and non-active objects. EBs are normally fast rotators ( $P_{\rm rot}$ < 3 days) and exhibit strong X-ray emission related to the presence of magnetic fields (López-Morales 2007). As explained by e.g., Chabrier et al. (2007), rapid rotation and/or magnetic fields reduce the efficiency of large-scale thermal convection in their interior leading to less efficient heat transport. The reduction in the area of the stellar radiating surface due to the presence of cold spots can cause lower values of $T_{\rm eff}$ and larger values of the radius, thus departing significantly from model predictions (see e.g., Torres & Ribas 2002).

To investigate this problem in more detail, we have studied the radius vs. mass relation for HIP 96515 A and compared this with the predictions of evolutionary models by Baraffe et al. (1998). The adopted models correspond to a solar metallicity of [Fe/H] = 0, and a mixing length of $\alpha$ = 1. The results are displayed in Fig. 9, which includes all the low-mass eclipsing binaries confirmed so far. The data was adopted from Shkolnik et al. (2008) and references therein.

By assuming an age of 400 Myr for the triple system, the radii of HIP 96515 Aa & Ab were found to be $\sim $15% and 10% larger, respectively, than predicted by the evolutionary tracks (see Fig. 9). These discrepancies were previously observed in similar eclipsing binaries and attributed to the effect of rapid rotation, magnetic activity, and/or metallicity (e.g., Berger et al. 2006; López-Morales 2007; Chabrier et al. 2007).

 \begin{figure}
\par\includegraphics[width=8.5cm,clip]{1514fig9.eps}
\end{figure} Figure 9:

Mass-luminosity relation for all known eclipsing binaries with masses M < 0.7 $M_{\odot }$ (see Shkolnik et al. 2008, and references therein). HIP 96515 Aa and Ab are represented by filled circles. The dotted, dashed, and solid lines represent evolutionary tracks by Baraffe et al. (1998) for 50, 100, and 500 Myr, respectively.

Open with DEXTER

6 Results and conclusions

We have studied in detail HIP 96515, a triple system consisting of a double-lined spectroscopic binary (SB2) and a comoving white dwarf with a projected separation of 8 $.\!\!^{\prime\prime}$6. We have analyzed both new and archival observations of the system. Our main results can be summarized as follows:

  • The analysis of multi-epoch optical photometry of HIP 96515 A has revealed that this SB2 is an eclipsing binary with almost equal-mass components: MAa = 0.59 $M_{\odot }$ and MAb = 0.54 $M_{\odot }$. This is the fourteenth confirmed eclipsing binary with component masses below 0.7 $M_{\odot }$, all of them being of extreme importance to calibrate theoretical evolutionary tracks.

  • Archival astrometry and new NACO/VLT data show that HIP 96515 A&B is a common proper-motion pair.

  • The optical spectrum of HIP 96515 B is consistent with that of a hot DB white dwarf. We have modeled the spectrum using evolutionary cooling sequence models, deriving a temperature, surface gravity, and mass of $T_{\rm eff}$ (K) = 19 126, $\log g$ = 8.08, and M = 0.65 $M_{\odot }$, respectively. Empirical relations provide a progenitor mass of 3.1 $\pm $ 0.8 $M_{\odot }$, and a white dwarf cooling age of 100 $\pm $ 13 Myr. The MS lifetime of a 3.1 $M_{\odot }$ star, as derived from evolutionary models, is $\sim $300 Myr, which corresponds to a total age of $\sim $400 Myr for the white dwarf.

  • A direct comparison of the EB members with evolutionary tracks (Landin et al. 2006; Baraffe et al. 1998) places them on the PMS with an age of $\sim $60 Myr, although with large uncertainties. We note that HIP 96515 A does not share the kinematical properties of known young moving groups in the solar vicinity.

  • If HIP 96515 A&B are coeval, the system must be older than 300 Myr. The discrepancy between this age and the 60 Myr provided by PMS tracks might be explained by the different evolution of EBs, which normally exhibit fast rotation and stronger magnetic fields than non-magnetic stars. Since HIP 96515 A shows fast rotation and strong X-ray emission, it can depart from the predictions of theoretical models, which place the binary on the PMS.

  • If we assume an age of $\sim $400 Myr for the triple system, the comparison of HIP 96515 Aa&Ab with evolutionary tracks by Baraffe et al. (1998) shows that the models underestimate the stellar radii by 15% and 10%, for the primary and the secondary, respectively. This behavior has been previously observed in similar low-mass EBs.
Finally, we note that additional RV observations are needed to constrain more robustly the orbital and physical parameters of the eclipsing binary.

Acknowledgements
We are very grateful to I. Saviane for obtaining the NTT/EMMI spectrum presented here. N.H. is indebted to the Spanish Programa Juan de la Cierva. This research has been funded by Spanish grant MEC ESP2007-65475-C02-02 and MEC/Consolider-CSD2006-0070. N.H. and C.T. gratefully acknowledge support from ESO-DGDF 2007 program. G.C. acknowledges support from the Faculty of the European Space Astronomy Centre (ESAC). H.B. acknowledges the funding from the European Commission's Sixth Framework Program as a Marie Curie Outgoing International Fellow (MOIF-CT-2005-8389). L.P.R.V. and N.R.L. gratefully acknowledge financial support from the Brazilian agencies CAPES, CNPq and FAPEMIG. P.B. is a Cottrell Scholar of the Research Corporation for Science Advancement and is supported in part by the NSERC Canada and by the Fund FQRNT(Québec). This research has made use of the SIMBAD database, operated at CDS, Strasbourg, France, of NASA's Astrophysics Data System Abstract Service, and of the Washington Double Star Catalog maintained at the US Naval Observatory.

References

Footnotes

... companion[*]
Based on observations collected at the Paranal Observatory under programs 77.C-0483(A) and 81.C-0826(A).
...[*]
Marie Curie Outgoing International Fellow.
... ASAS[*]
http://www.astrouw.edu.pl/~gp/asas/asas.html
... dwarfs[*]
See http://www.astro.umontreal.ca/~bergeron/CoolingModels

All Tables

Table 1:   Main properties of HIP 96515 A.

Table 2:   Photometric data of the visual binary HIP 96515 A&B.

Table 3:   Spectroscopic data of the eclipsing binary HIP 96515 A.

Table 4:   Grid of solutions for 7 values of the mass ratio, q(=M2/M1), for HIP 96515 A.

Table 5:   Preliminary physical parameters for HIP 96515 Aa&Ab.

Table 6:   Relative astrometry of HIP 96515 A & B.

Table 7:   Best-fit parameters for HIP 96515 B.

All Figures

  \begin{figure}
\par\includegraphics[width=9cm,clip]{1514fig1.eps}
\end{figure} Figure 1:

Hipparcos and ASAS V-band light-curves of HIP 96515 A.

Open with DEXTER
In the text

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

The phase-folded lightcurves of HIP 96515 A from ASAS (+ symbols) and Hipparcos ($\times $ symbols), consisting of 818 and 96 observations, respectively, with our theoretical solution for q=0.9 obtained in Sect. 3.2. The Hipparcos observations were shifted by 0.6 mag, in order to separate the light curves. Note that ASAS observations were not corrected from the long term variation of the maxima, shown in Fig. 1.

Open with DEXTER
In the text

  \begin{figure}
\par\includegraphics[width=9cm,clip]{1514fig3.eps}
\end{figure} Figure 3:

The residuals O-C between the V light curves of HIP 96515 A from Hipparcos ($\times $ symbols) and ASAS (+ symbols) shown in Fig. 2, and our theoretical solution for q = 0.9 obtained in Sect. 3.2. The standard error (1 observation) is shown in the upper left part of each panel.

Open with DEXTER
In the text

  \begin{figure}
\par\includegraphics[width=9cm,clip]{1514fig4.eps}
\end{figure} Figure 4:

The residuals O-C between the V light curves of HIP 96515 A from ASAS (+ symbols) and Hipparcos (grey circles) shown in Fig. 3, folded with the period of (0.24965 $\pm $ 0.00001) day. Note that the pattern is present in both Hipparcos and ASAS data sets (two cycles are plotted).

Open with DEXTER
In the text

  \begin{figure}
\par\includegraphics[width=9cm,clip]{1514fig5.ps}
\end{figure} Figure 5:

VLT/NACO $K_{\rm s}$-band image of HIP 96515 A and its visual companion, HIP 96515 B, obtained in June 2008. The projected separation is of 378 AU at a distance of 44 pc. Note that the eclipsing binary (Aa & Ab) is not spatially resolved by NACO.

Open with DEXTER
In the text

  \begin{figure}
\par\includegraphics[width=7.5cm,clip]{1514fig6.eps}
\end{figure} Figure 6:

Differences in right ascention and declination between HIP 96515 A & B at different epochs (black circles). The expected position of a background object (grey circles and labels) at each epoch is also displayed. Top: all the astrometric data from 1910 to 2008. Bottom: only the NACO data, which provides the highest astrometric accuracy. The figure shows that the two objects are comoving, that is, they are a common proper motion pair.

Open with DEXTER
In the text

  \begin{figure}
\par\includegraphics[width=8.5cm,clip]{1514fig7.ps}
\end{figure} Figure 7:

Top: NTT/EMMI Optical spectrum of HIP 96515 B. We have marked all the identified He I lines. Bottom: fit of the optical spectrum using pure helium model atmospheres. The best fit parameters are displayed at the top right corner.

Open with DEXTER
In the text

  \begin{figure}
\par\includegraphics[width=9cm,clip]{1514fig8.eps}
\end{figure} Figure 8:

Stellar luminosity ( top panel) and surface gravity ( bottom panel) versus effective temperature. HIP 96515 Aa and Ab are represented by filled circles. The grey and black solid lines represent evolutionary tracks by Baraffe et al. (1998) and Landin et al. (2006), respectively, for 0.6 and 0.5 $M_{\odot }$ stellar masses. The dotted lines represent the uncertainties in the stellar masses for the ATON models. The comparison of HIP 96515 A with the two sets of evolutionary tracks provides an age of $\sim $60 Myr, that is, places the eclipsing binary members on the pre-main sequence, although with large uncertainties.

Open with DEXTER
In the text

  \begin{figure}
\par\includegraphics[width=8.5cm,clip]{1514fig9.eps}
\end{figure} Figure 9:

Mass-luminosity relation for all known eclipsing binaries with masses M < 0.7 $M_{\odot }$ (see Shkolnik et al. 2008, and references therein). HIP 96515 Aa and Ab are represented by filled circles. The dotted, dashed, and solid lines represent evolutionary tracks by Baraffe et al. (1998) for 50, 100, and 500 Myr, respectively.

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.