A&A 484, 413-418 (2008)
DOI: 10.1051/0004-6361:20078381
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 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 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 -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 2 to be a very close binary of
0.16 arcsec separation, having a flux ratio of
0.91, thus having almost equal brightness and indistinguishable spectral types within the errors. We show that the two tentative components of Cha H
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
100 K, corresponding to a spectral type of M6
1 and a surface gravity of
,
both from a comparison with GAIA model atmospheres. Furthermore, we obtained an optical extinction of
4.3 mag from this comparison.
Conclusions. We derive masses of 0.110
(
0.070
)
and
0.124
(
0.077
)
for the two components of Cha H
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 2
Neuhäuser & Comerón (1998) found Cha H 2, also called ISO-ChaI 111, to be a member of the Cha I star-forming cloud with an age of
2 Myr. Neuhäuser & Comerón (1999) then classified Cha H
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
2.
Neuhäuser et al. (2002) searched for companions around Cha H
1-12
with the Hubble Space Telescope and presented the first indication for binarity of
Cha H
2 due to elongation of the Point Spread Function (PSF) in filters R and H
.
The rotational period of 3.21
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
2 signs of accretion in H
and Pa
,
making it the only object out of nine to show clear signs of accretion
at a rate of
/yr and an inclination of
65 to
75
.
Observations with XMM-Newton by Stelzer et al. (2004)
revealed X-ray emission from Cha H
2.
Recently, Apai et al. (2005) succeeded in using new infrared spectra
to conclude that the disk around Cha H
2 is still closer to
the flat disk model, even though a prominent enstatite silicate feature
could be found at 9.3
m.
Here, we present evidence for the binarity of Cha H 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.
We observed Cha H 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.
![]() |
Figure 1:
NACO Ks-band image of the low-mass binary Cha H![]() |
Open with DEXTER |
Table 1: VLT/NACO observations log.
Table 2:
Proper motions of Cha H 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 (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 2 is clearly resolved in a double object with an apparent separation of
160 mas at a position angle of
40
,
see Fig. 1. The southwestern component appears slightly brighter, hence we call
this component A (or SW), the fainter northeastern component B (or NE).
![]() |
Figure 2:
Observed separation between the two components of Cha H![]() |
Open with DEXTER |
![]() |
Figure 3:
Observed position angle between the two components of Cha H![]() |
Open with DEXTER |
To check for common proper motion of the two tentative components of Cha H 2, we used the proper motion (here PM) of Cha H
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
50
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 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 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
and 1.5
(Fig. 2); 2.2
and 2.3
(Fig. 3); 1.4
and 1.0
(Fig. 4); 2.1
and 1.9
(Fig. 5)
that either one of the tentative components of Cha H
2 is a non-moving background object.
Viz, for both components the background hypothesis is rejected by
3.7
,
respectively.
For two objects slightly above the sub-stellar mass limit (see below) and the given separation (at
168 pc), the orbital period is
110 yrs, so that the maximum change in separation due to orbital motion (for circular edge-on orbit) is
12 mas/yr (
/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
significance, we detect orbital motion as deviation of the separation in 2005 and 2006 in contrast to the 2007 measurement, with 1.9
and 3.9
significance. From the negligible change in PA during our
2 years of epoch difference, we can exclude a circular pole on orbit by 2.8
,
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
and 2.3
for component A and B not being a non-moving background star, respectively.
![]() |
Figure 4:
Observed separation between the two components of Cha H![]() |
Open with DEXTER |
![]() |
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.
![]() |
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, ![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
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 2 from the Two Micron All Sky Survey (2MASS) catalogue of J = 12.210
0.024 mag and K = 10.675
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.
![]() |
Figure 7:
Plot of the true and expected position of Cha H![]() ![]() ![]() |
Open with DEXTER |
We combined all available spectra of Cha H 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
m to 2.45
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 = 1 is free of dust. Since little or no dust forms at temperatures above
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
100 K,
a surface gravity
= 4.0
+0.75-0.5 and a visual extinction of
= 4.3
0.25 mag, see Fig. 6 and the online material for a color version.
We note that our extinction value is in between the
of 7.9 mag measured in a slope in the
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
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
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
Cha 10 in Luhman (2004).
We also obtained low-S/N (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 2 system.
Table 5:
Brightnesses of the Cha H 2 system.
![]() |
Figure 8:
Observed K-band NACO spectra of the Cha H![]() ![]() |
Open with DEXTER |
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
of 3.0
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
of -1.21
0.14 and -1.25
0.14, and radii of 0.92
-0.14 +0.19 and 0.88
for Cha H
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
0.005
(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
0.17 days (Joergens et al. 2003) and the
v sin i of 12.8
1.2 km s-1 (Joergens & Guenther 2001)
that Cha H
2 A and B have inclinations of
36.1 to 143.9
and 40.4 to 139.6
,
respectively,
consistent with
65 to 75
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
and 0.281
for A and B,
thus leaving Cha H
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
and 0.110
for Cha H
2 A and B for an age of 2 Myr, a temperature of
3025 K, a
of
3.6, and the luminosity values as derived above of log
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
2 of 0.14
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 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
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 27 AU at
168 pc,
so that further investigation of the interaction between the two components in
Cha H
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 H2 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.
![]() |
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, ![]() ![]() ![]() ![]() |
Open with DEXTER |