A&A 484, 413-418 (2008)
DOI: 10.1051/0004-6361:20078381

Confirmation of the binary status of Chamaeleon H$\alpha $ 2 - a very young low-mass binary in Chamaeleon[*],[*]

T. O. B. Schmidt1 - R. Neuhäuser1 - N. Vogt2,3 - A. Seifahrt1 - T. Roell1 - A. Bedalov1


1 - Astrophysikalisches Institut und Universitäts-Sternwarte, Universität Jena, Schillergäßchen 2-3, 07745 Jena, Germany
2 - Departamento de Física y Astronomía, Universidad de Valparaíso, Avenida Gran Bretaña 1111, Valparaíso, Chile
3 - Instituto de Astronomía, Universidad Catolica del Norte, Avda. Angamos 0610, Antofagasta, Chile

Received 30 July 2007 / Accepted 29 January 2008

Abstract
Context. Neuhäuser & Comerón (1998, Science, 282, 83; 1999, A&A, 350, 612) presented direct imaging evidence, as well as first spectra, of several young stellar and sub-stellar M6- to M8-type objects in the Cha I dark cloud. One of these objects is Cha H$\alpha $ 2, classified as brown dwarf candidate in several publications and suggested as possible binary in Neuhäuser et al. (2002, A&A, 384, 999).
Aims. We have searched around Cha H$\alpha $ 2 for close and faint companions with adaptive optics imaging.
Methods. Two epochs of direct imaging data were taken with the Very Large Telescope (VLT) Adaptive Optics instrument NACO in February 2006 and March 2007 in $K{\rm s}$-band together with a Hipparcos binary for astrometric calibration. Moreover, we took a J-band image in March 2007 to get color information. We retrieved an earlier image from 2005 from the European Southern Observatory (ESO) Science Archive Facility, increasing the available time coverage. After confirmation of common proper motion, we deduce physical parameters of the objects by spectroscopy, like temperature and mass.
Results. We find Cha H$\alpha $ 2 to be a very close binary of $\sim$0.16 arcsec separation, having a flux ratio of $\sim$0.91, thus having almost equal brightness and indistinguishable spectral types within the errors. We show that the two tentative components of Cha H$\alpha $ 2 form a common proper motion pair, and that neither component is a non-moving background object. We even find evidence for orbital motion. A combined spectrum of both stars spanning optical and near-infrared parts of the spectral energy distribution yields a temperature of 3000 $\pm$ 100 K, corresponding to a spectral type of M6 $\pm$ 1 and a surface gravity of $\log{g}= 4.0^{+0.75}_{-0.5}$, both from a comparison with GAIA model atmospheres. Furthermore, we obtained an optical extinction of $A_{\rm V}$ $\simeq$ 4.3 mag from this comparison.
Conclusions. We derive masses of $\sim$0.110 $M_{\odot}$ ($\geq$0.070 $M_{\odot}$) and $\sim$0.124 $M_{\odot}$ ($\geq$0.077 $M_{\odot}$) for the two components of Cha H$\alpha $ 2, i.e., probably low-mass stars, but one component could possibly be a brown dwarf.

Key words: stars: low-mass, brown dwarfs - stars: binaries: close - stars: binaries: visual - stars: individual: Cha H$\alpha $ 2

1 Introduction

Neuhäuser & Comerón (1998) found Cha H$\alpha $ 2, also called ISO-ChaI 111, to be a member of the Cha I star-forming cloud with an age of $\sim$2 Myr. Neuhäuser & Comerón (1999) then classified Cha H$\alpha $ 2 as a candidate brown dwarf of spectral type M6.5, using medium-resolution optical spectroscopy. Comerón et al. (2000) presented evidence for near- to mid-infrared excess indicating a disk. Natta & Testi (2001) and Apai et al. (2002) argued that the spectrum can be explained with an optically flat dust disk, but not with a flared disk, due to lack of any prominent silicate feature previously expected for the disk of Cha H$\alpha $ 2. Neuhäuser et al. (2002) searched for companions around Cha H$\alpha $ 1-12 with the Hubble Space Telescope and presented the first indication for binarity of Cha H$\alpha $ 2 due to elongation of the Point Spread Function (PSF) in filters R and H$\alpha $.

The rotational period of 3.21 $\pm$ 0.17 days, found in photometric data by Joergens et al. (2003), confirmed the spectroscopic period of 2.9 +1.4-1.0 days already proposed in Joergens & Guenther (2001). Natta et al. (2004) observed nine stars in the Cha I dark cloud and found only at Cha H$\alpha $ 2 signs of accretion in H $\alpha $ and Pa $\beta$, making it the only object out of nine to show clear signs of accretion at a rate of $\sim$ $10^{-10}~M_{\odot}$/yr and an inclination of $\sim$65 to 75$^{\circ}$. Observations with XMM-Newton by Stelzer et al. (2004) revealed X-ray emission from Cha H$\alpha $ 2. Recently, Apai et al. (2005) succeeded in using new infrared spectra to conclude that the disk around Cha H$\alpha $ 2 is still closer to the flat disk model, even though a prominent enstatite silicate feature could be found at 9.3 $\mu $m.

Here, we present evidence for the binarity of Cha H$\alpha $ 2. In Sects. 2 and 3, we present the observations, data reduction, and astrometric results. The available photometric properties of both components will be discussed in Sect. 4. The spectroscopy of the system is presented in Sect. 5. We end with conclusions in Sect. 6.

2 Observations with VLT/NACO

We observed Cha H$\alpha $ 2 in two epochs in February 2006 and in March 2007, see Table 1 for the observations log. Included in Table 1 are the data found in the archive from March 2005.

  \begin{figure}
\par\includegraphics[width=7.4cm,clip]{8381fig1.eps} \end{figure} Figure 1: NACO Ks-band image of the low-mass binary Cha H$\alpha $ 2 from 25 March 2005. The southwestern component appears slightly brighter, hence we call this component A (or SW), the fainter northeastern component B (or NE).
Open with DEXTER

Table 1: VLT/NACO observations log.

Table 2: Proper motions of Cha H$\alpha $ 2.

We observed with the European Southern Observatory (ESO) Very Large Telescope (VLT) instrument Naos-Conica (NACO, Lenzen et al. 2003; Rousset et al. 2003). DIT, NDIT, NINT, and filter bands are listed in Table 1. In all cases, we used the S13 camera ($\sim$13 mas/pixel pixel scale) and the double-correlated read-out mode.

For the raw data reduction, we subtracted a mean dark from all the science frames and flatfield frames, then divided by the mean dark-subtracted flatfield and subtracted the mean background using eclipse / jitter.

In all three images, Cha H$\alpha $ 2 is clearly resolved in a double object with an apparent separation of $\sim$160 mas at a position angle of $\sim$40$^{\circ}$, see Fig. 1. The southwestern component appears slightly brighter, hence we call this component A (or SW), the fainter northeastern component B (or NE).


  \begin{figure}
\par\includegraphics[width=8cm,clip]{8381fig2.eps} \end{figure} Figure 2: Observed separation between the two components of Cha H$\alpha $ 2. The archival data point from 2005 and our two measurements from 2006 and 2007 are shown. The short-dashed lines enclose the area for constant separation. The dash-dotted line is the change expected if the B component is a non-moving background star. The opening cone enclosed by the solid lines are its estimated errors. The waves of this cone show the differential parallactic motion, which has to be taken into account if the other component is a non-moving background star. The opening long-dashed cone indicates the amplitude of possible orbital motion.
Open with DEXTER


  \begin{figure}
\par\includegraphics[width=7.8cm,clip]{8381fig3.eps} \end{figure} Figure 3: Observed position angle between the two components of Cha H$\alpha $ 2 (measured from the brighter component north over east to south). See Fig. 2 for details.
Open with DEXTER

3 Astrometry

To check for common proper motion of the two tentative components of Cha H$\alpha $ 2, we used the proper motion (here PM) of Cha H$\alpha $ 2 published in the literature (USNO-B and PSSPMC); we calculated additional values by comparing SuperCOSMOS Sky Survey (here SSS) data from 1985 with images from VLT/FORS1 (epoch 1999, published by Comeron et al. 2000, obtained directly from F. Comerón, priv. com.) and New Technology Telescope instrument Son of ISAAC (NTT/SofI) from 2000, observed by us (10 $\times$ 50 $\times$ 1.2 s integrations), reduced in the same way as the NACO data with eclipse. See Table 2 for the results (originally from Gaedke 2005, but revised here by us). We use the mean proper motion for checking, whether the two objects show common proper motion below.

Table 3: Astrometric calibration and astrometric results for Cha H$\alpha $ 2 AB.

We calibrated the NACO data using the wide binary star HIP 73357 for our two measurements in 2006 and 2007, and the closer binary HIP 73111 for which we found archived data taken in the same night as the Cha H$\alpha $ 2 observations in 2005, resulting in the astrometric calibration summarized in Table 3. The error bars of pixel scale and orientation include, possible orbital motion of the respective binary since its measurement by Hipparcos - i.e., a maximum change in separation due to possible orbital motion for circular edge-on orbit and a maximum change in position angle for circular pole-on orbit - as well as the uncertainties in parallax, total mass, position measurements in our images and separation at the epoch of Hipparcos of the respective binary system.

To determine the positions of both components we constructed a reference PSF from the binary itself by using the undeformed remote sides of both objects shifted to the same peak position. Thus, we obtained a clean reference PSF for each single image. With IDL/starfinder, we scaled and shifted the reference PSF simultaneously to both components in each of our individual images by minimizing the residuals. As a result, we obtained positions (see Table 3) and relative photometry (see Table 4), including realistic error estimates for each object by averaging the results of all single images taken within each epoch.

The development of separation and position angle with time is shown in Figs. 2-5. We can exclude by 0.9$\sigma $ and 1.5$\sigma $ (Fig. 2); 2.2$\sigma $ and 2.3$\sigma $ (Fig. 3); 1.4$\sigma $ and 1.0$\sigma $ (Fig. 4); 2.1$\sigma $ and 1.9$\sigma $ (Fig. 5) that either one of the tentative components of Cha H$\alpha $ 2 is a non-moving background object. Viz, for both components the background hypothesis is rejected by $\geq$3.7$\sigma $, respectively. For two objects slightly above the sub-stellar mass limit (see below) and the given separation (at $\sim$168 pc), the orbital period is $\sim$110 yrs, so that the maximum change in separation due to orbital motion (for circular edge-on orbit) is $\sim$12 mas/yr ($\sim$$3^{\circ}$/yr in PA for pole-on orbit). While orbital motion could be detected in PA only from one point in 2006 in contrast to 2007, with 1.9$\sigma $ significance, we detect orbital motion as deviation of the separation in 2005 and 2006 in contrast to the 2007 measurement, with 1.9$\sigma $ and 3.9$\sigma $ significance. From the negligible change in PA during our $\sim$2 years of epoch difference, we can exclude a circular pole on orbit by 2.8$\sigma $, however, we cannot rule out a highly eccentric pole-on orbit.

In Fig. 7, the data from Hubble Space Telescope by Neuhäuser et al. (2002) are also included. Only the plot for the assumption of component B being a non-moving background star is shown here. The Hubble Space Telescope data hence yield another 2.2$\sigma $ and 2.3$\sigma $ for component A and B not being a non-moving background star, respectively.


  \begin{figure}
\par\includegraphics[width=8cm,clip]{8381fig4.eps} \end{figure} Figure 4: Observed separation between the two components of Cha H$\alpha $ 2. Same as Fig. 2, but the opening cone enclosed by the solid lines is the change expected if the A component is a non-moving background star.
Open with DEXTER


  \begin{figure}
\par\includegraphics[width=8cm,clip]{8381fig5.eps} \end{figure} Figure 5: Same as Fig. 3, but position angle between the two components measured from the B component and assuming the A component is a non-moving background star.
Open with DEXTER

While the difference seen in PA and the separation between the different observations are consistent with common proper motion and some orbital motion, a possible difference in proper motion between both objects of up to a few mas/yr cannot be excluded from the data. Such a difference in proper motion is typical for the velocity dispersion in star forming regions like Cha I (Ducourant et al. 2005), so that we cannot yet exclude that both objects are independent members of Cha I, but not orbiting each other. However, it is extremely unlikely that two objects with almost the same temperatures, J and K magnitudes, gravities, radii, proper motions, and extinctions are located within less than 170 mas. Even if this would be the case, age, and distance would be the same as assumed below, hence also the mass estimation.


  \begin{figure}
\par\includegraphics[angle=90,width=16.6cm,clip]{8381fig6.eps} \end{figure} Figure 6: Bottom: spectrum of both components from the optical part to the near-infrared combined from data from Comerón et al. (2000) and data from the ESO Paranal and La Silla instruments ISAAC and SofI retrieved from the ESO Science Archive Facility. Top: the best fit to the spectrum: GAIA models computed with the PHOENIX code by Brott & Hauschildt (2005) for a temperature of 3000 K, $\log{g}$ of 4.0 and a visual extinction  $A_{\rm V}\simeq 4.3$ mag. Note that the missing flux depression in J-band (between 1.2 and 1.3 micron) as well as in the peak of the H-band (around 1.65 micron) is due to missing FeH opacities in the spectral models. Note the strong H$\alpha $ emission in Cha H$\alpha $ 2 at 0.6568 $\mu $m.
Open with DEXTER

4 Photometry

As described in the last section, we obtained from the PSF fitting of both components also the relative flux ratio, see Table 4. We considered these flux ratios in the astrometric proper motion analysis since the background hypotheses is based on propermotions obtained by centroid measurements of the unresolved system. Hence, if one component were a non-moving background object, the centroid (and thus the proper motion of the combined object) would only shift by about half the relative change in separation between both objects.

Using the unresolved photometry of Cha H$\alpha $ 2 from the Two Micron All Sky Survey (2MASS) catalogue of J = 12.210 $\pm$ 0.024 mag and K = 10.675 $\pm$ 0.021 mag, we obtain the photometry of each component using our measured flux ratios, see Table 5. We see no evidence for photometric variability.


  \begin{figure}
\par\includegraphics[width=8cm,clip]{8381fig7.eps} \end{figure} Figure 7: Plot of the true and expected position of Cha H$\alpha $ 2 A if Cha H$\alpha $ 2 B is a non-moving background star. Shown are the NACO measurements from 2005 to 2007, as in Figs. 2-5, and, in addition, the measurement from 2000 using the Hubble Space Telescope in R-band by Neuhäuser et al. (2002). As in Fig. 2, the line showing parallactic wobbles is the expected position of component A if component B is a non-moving background star, in this case from epoch 2000.8 until 2007.2, from left to right. The opening dashed cone are its errors. The HST data point obtained at epoch 2000.8 is significantly different (2.3$\sigma $) from the expectation, if background.
Open with DEXTER

5 Spectroscopy

We combined all available spectra of Cha H$\alpha $ 2 from the archive and an optical spectrum obtained directly from F. Comerón (priv. com.) published in Comerón et al. (2000) to create a combined spectrum with a nearly complete coverage from 0.64 $\mu $m to 2.45 $\mu $m. The data were reduced in a standard manner, doing flat field correction, sky subtraction, spectrum extraction, wavelength calibration, and standard correction with G stars. The spectra are not flux-calibrated. We found the relative offsets from the overlaps of the spectra.

We fitted this spectrum with GAIA COND models computed with the PHOENIX code by Brott & Hauschildt (2005). COND represents the case where dust forms, as in the similar DUSTY models, but rains out completely from the photosphere into deeper layers and, hence, the region where $\tau$ = 1 is free of dust. Since little or no dust forms at temperatures above $\sim$2400 K both models with different treatment of dust formation give the same result at the given temperatures. Our best fit results give an effective temperature of 3000 $\pm$ 100 K, a surface gravity $\log{g}$ = 4.0 +0.75-0.5 and a visual extinction of $A_{\rm V}$ = 4.3 $\pm$ 0.25 mag, see Fig. 6 and the online material for a color version. We note that our extinction value is in between the $A_{\rm V}$ of 7.9 mag measured in a slope in the $A_{\rm V}$ map from Kainulainen et al. (2006) and the first ever measured value of 3.55 mag by Neuhäuser & Comerón (1998). Given the quality of our fit, we find here further evidence for the physical binarity of Cha H$\alpha $ 2 being a double of almost identical mid- to late-M-type objects. If one of the components would be a background (or foreground) object, both the spectral type and the extinction of the components would differ and the unresolved spectrum of the binary would not be consistent with one spectral type over the whole spectral range.

The derived temperature can be converted to a spectral type of M6 $\pm$ 1 with the empirical temperature scale for young M-type objects from Luhman (1999). Moreover, our combined spectrum in Fig. 6 resembles very much the spectrum of the M5.75 object $\epsilon$ Cha 10 in Luhman (2004).

We also obtained low-S/N ($\sim$15) K-band spectra of both objects with VLT/NACO, see Fig. 8. They are indistinguishable, giving further evidence that neither of the two objects are background or foreground stars, which would exhibit a significantly different K-band spectrum.

Table 4: Flux ratios of the Cha H$\alpha $ 2 system.

Table 5: Brightnesses of the Cha H$\alpha $ 2 system.


  \begin{figure}
\par\includegraphics[angle=90,width=7.8cm,clip]{8381fig8.eps} \end{figure} Figure 8: Observed K-band NACO spectra of the Cha H$\alpha $ 2 system. In each panel, one of the components of the system is shown ( top in each panel) compared to the higher S/N spectrum of both objects combined from Fig. 6 ( bottom in each panel). Top: component A (southwestern component). Bottom: component B (northeastern component). Both spectra are dereddened using the visual extinction of $A_{\rm V}\simeq 4.3$ mag derived from the spectra shown in Fig. 6.
Open with DEXTER

6 Conclusions

From the photometry, the visual extinction, an extinction law to convert the visual to a near-infrared extinction by Rieke & Lebofsky (1985), a bolometric correction BC$_{\rm K}$ of 3.0 $\pm$ 0.1 mag from Golimowski et al. (2004) and the distance to a particular member of the Cha I cloud of 168 +28-24 from Bertout et al. (1999) (calculated by HIPPARCOS measurements of 4 members of Cha I and resulting in an average distance of 168 +14-12 to the cloud), we find luminosities log $(L/L_{\odot})$ of -1.21 $\pm$ 0.14 and -1.25 $\pm$ 0.14, and radii of 0.92  -0.14 +0.19 and 0.88 $_{~-0.13}^{~+0.18}~R_{\odot}$ for Cha H$\alpha $ 2 A (SW) and B (NE), respectively. This makes both components more massive than the more massive component of the eclipsing binary brown dwarf found by Stassun et al. (2006) in Orion of M = 0.054 $\pm$ 0.005 $M_{\odot}$ (assuming similar age, which is justified).

It is, a priori, not known to which of both components the photometric rotation period refers. However, since both stars are very young and have been formed in similar environmental conditions, we may assume that they rotate at the same speed. In this case, we can conclude from the photometric period of 3.21 $\pm$ 0.17 days (Joergens et al. 2003) and the v sin i of 12.8 $\pm$ 1.2 km s-1 (Joergens & Guenther 2001) that Cha H$\alpha $ 2 A and B have inclinations of 36.1 to 143.9 $\deg$ and 40.4 to 139.6 $\deg$, respectively, consistent with $\sim$65 to 75 $\deg$ by Natta et al. (2004). Assuming an alignment of the orbital axis and the rotational axis of both components, this is also consistent with our finding of astrometric orbital motion in the form of a change in separation rather than in a change of the PA, i.e., with an orbital plane being more edge-on than pole-on.

Given the luminosities, temperatures, radii, and gravities, the masses of both components are 0.308 $^{~+1.430}_{~-0.231}~M_{\odot}$ and 0.281 $^{~+1,302}_{~-0.210}~M_{\odot}$ for A and B, thus leaving Cha H$\alpha $ 2 AB above the stellar - brown dwarf boundary. However, component B could still be a brown dwarf, according to this error estimation. Although we arrive at masses that could possibly be higher than the sun's (given our large error bars), this can be ruled out due to the fact that such very young objects could be of approximately same size as the sun, but would not have a similar surface gravity at the same time. Using theoretical models by Baraffe et al. (1998), we find a best fit for masses of 0.124 $M_{\odot}$ and 0.110 $M_{\odot}$ for Cha H$\alpha $ 2 A and B for an age of 2 Myr, a temperature of $\sim$3025 K, a $\log~g$ of $\sim$3.6, and the luminosity values as derived above of log $(L/L_{\odot})$ of -1.21 and -1.25, but we would like to mention that the determination of ages and masses from evolutionary models is uncertain until at least 10 Myr (Chabrier et al. 2005). These mass values are in good agreement with the mass for the unresolved Cha H$\alpha $ 2 of 0.14 $M_{\odot}$ from Mohanty et al. (2005). While there are currently almost one hundred systems known of even lower total mass (Siegler 2007), a review of the field and the existence of very low mass binaries at different ages can be found in Burgasser et al. (2007).

After submission of our manuscript Ahmic et al. (2007) reported the possible binarity of Cha H$\alpha $ 2 found in their data from 2005, included as 1st epoch in our common proper motion analysis. We note that separation, position angle, and Ks-band magnitudes measured coincide well within the 1$\sigma $ error bars of each other. Ahmic et al. (2007) could not show common proper motion because they had only one epoch available.

The projected separation between A and B corresponds to $\sim$27 AU at $\sim$168 pc, so that further investigation of the interaction between the two components in Cha H$\alpha $ 2 and its disk (with silicate feature) would be very important.

Acknowledgements
T.O.B.S. acknowledges support from Evangelisches Studienwerk e.V. Villigst.
Moreover he would like to thank F. Comerón for providing FORS1 images as well as an electronic optical spectrum of Cha H$\alpha $ 2 from his publication Comerón et al. (2000), P. Hauschildt for the electronic GAIA models, J. Kainulainen for the electronic extinction map of the Cha I cloud for accurate reading of the extinction values, an anonymous referee for helpful comments and Amy Mednick for language editing.
N.V. acknowledges support by FONDECYT grant 1061199.
A.B. and R.N. would like to thank DFG for financial support in projects NE 515/ 13-1 and 13-2.
This publication makes use of data products from the Two Micron All Sky Survey, which is a joint project of the University of Massachusetts and the Infrared Processing and Analysis Center/California Institute of Technology, funded by the National Aeronautics and Space Administration and the National Science Foundation.
We use imaging data from the SuperCOSMOS Sky Survey, prepared and hosted by Wide Field Astronomy Unit, Institute for Astronomy, University of Edinburgh, which is funded by the UK Particle Physics and Astronomy Research Council.
This research has made use of the VizieR catalogue access tool and the Simbad database, both operated at the Observatoire Strasbourg.

References

 

  
Online Material


  \begin{figure}
\par\includegraphics[angle=90,width=16.6cm,clip]{8381Onl6.eps} \end{figure} Figure 6: Spectrum of both components from the optical part to the near-infrared combined from data from Comerón et al. (2000) (purple), and data from the ESO Paranal and La Silla instruments ISAAC (blue) and SofI (black: blue arm spectra, red: red arm spectra) retrieved from the ESO Science Archive Facility. The best fit to the spectrum, a GAIA model computed with the PHOENIX code by Brott & Hauschildt (2005) for a temperature of 3000 K, $\log~g$ of 4.0 and a visual extinction of 4.3 mag is superimposed in green. Note that the missing flux depression in J-band (between 1.2 and 1.3 micron) as well as in the peak of the H-band (around 1.65 micron) is due to missing FeH opacities in the spectral models. Note the strong H$\alpha $ emission in Cha H$\alpha $ 2 at 0.6568 $\mu $m.
Open with DEXTER



Copyright ESO 2008