Issue 
A&A
Volume 541, May 2012



Article Number  A29  
Number of page(s)  9  
Section  Stellar structure and evolution  
DOI  https://doi.org/10.1051/00046361/201118707  
Published online  24 April 2012 
Orbits and masses in the multiple system LHS 1070 ^{⋆}
^{1} MaxPlanckInstitut für Astronomie, Königstuhl 17, 69117 Heidelberg, Germany
email: koehler@mpia.de
^{2} Landessternwarte, Zentrum für Astronomie der Universität Heidelberg, Königstuhl, 69117 Heidelberg, Germany
email: r.koehler@lsw.uniheidelberg.de
^{3} UniversitätsSternwarte München, LudwigMaximiliansUniversität, Scheinerstr. 1, 81679 München, Germany
Received: 22 December 2011
Accepted: 15 March 2012
Aims. We present a study of the orbits of the triple system LHS 1070, with the aim to determine individual masses of its components.
Methods. Sixteen new relative astrometric positions of the three components in the K band were obtained with NACO at the VLT, Omega CASS at the 3.5 m telescope on Calar Alto, and other highspatialresolution instruments. We combined them with data from the literature and fit orbit models to the dataset. We derived an improved fit for the orbit of LHS 1070 B and C around each other, and an estimate for the orbit of B and C around A.
Results. The orbits are nearly coplanar, with a misalignment angle of less than 10°. The masses of the three components are M_{A} = 0.13...0.16 M_{⊙}, M_{B} = 0.077 ± 0.005 M_{⊙}, and M_{C} = 0.071 ± 0.004 M_{⊙}. Therefore, LHS 1070 C is certainly, and LHS 1070 B probably a brown dwarf. Comparison with theoretical isochrones shows that LHS 1070 A is either fainter or more massive than expected. One possible explanation would be that it is a binary. However, the close companion reported previously could not be confirmed.
Key words: stars: lowmass / brown dwarfs / stars: fundamental parameters / stars: individual: LHS 1070 / binaries: close / celestial mechanics
© ESO, 2012
1. Introduction
LHS 1070 (other common names are GJ 2005, LP 88164, 2MASS J002444192708242) is a nearby highpropermotion star located in the south galactic pole region. It was observed by Leinert et al. (1994) as part of their nearinfrared speckle survey for duplicity of nearby southern M dwarves. The group detected in 1993 two companions located 1.1′′ and 1.3′′ north of the primary by using the SHARP camera mounted at the NTT (La Silla, Chile) and applying Speckle imaging. In a consecutive measurement the companions exhibited an orbital motion around each other and around the primary while following the proper motion of the system. A fourth component D was identified by Henry et al. (1999) with the fine guidance sensors onboard the Hubble Space Telescope (HST). A separation of only about 50 mas from component A was reported. Speckle interferometric and adaptive optics measurements, however, were not able to confirm LHS 1070 D (Leinert et al. 2001; Seifahrt et al. 2008). This might be due to the small separation, but it is also possible that the detection is spurious (Henry, priv. comm.).
When comparing the colors of the components, the companions appear redder than the primary (Leinert et al. 1994). First estimates already identified the two companions as stars close to the hydrogen burning limit that divides brown dwarves from mainsequence stars. This mass range is characterised by a strong decrease of the effective temperature and the onset of dust formation in the atmospheres. Photometric and spectroscopic measurements with the HST in the visual indeed were represented by model atmospheres containing dust (Leinert et al. 2000). Also the low masses could be confirmed. The spectral classifications of the companions are M8.5V and M99.5V. The primary has an earlier spectral type of M5.56V.
Reiners et al. (2007) found for the components B and C a similarily high rotational velocity of vsini = 16 km s^{1}, which is twice that of the primary. The authors conclude that the temperaturedependent magnetic braking was acting for 1 Gyr on all three components. The main difference between the otherwise very similar components B and C is the higher activity of component B. This might be related to the higher temperature or the higher magnetic flux. In the HST spectra Hα emission was found towards LHS 1070 A and B, but not towards component C (Leinert et al. 2000). Recently, photometric activity of the two more massive stars was identified (Almeida et al. 2011). Component B showed a brightness increase in the visual that is among the largest ever observed in a flare star.
The very low mass companions of LHS 1070 provide the possibility to refine the models of (pre)mainsequence stars close to the transition region to brown dwarves. It is thus very important to derive the exact masses of these stars. The small separation of the components B and C, i.e. their short orbital period is very well suited to allow even for this low mass stars a mass determination by fitting the orbital elements. The dynamical mass has the advantage of being independent from theoretical models. It is thus a precious probe to test evolutionary models.
The main uncertainty for the mass determination is the distance to the object. The trigonometric parallax of LHS 1070 derived by van Altena et al. (1995) is 135.3 ± 12.1 mas. This value was later refined to 129.47 ± 2.48 mas, placing the system at a distance of 7.72 ± 0.15 pc (Costa et al. 2005). The proper motion of LHS 1070 is 653.7 ± 0.3 mas/yr, corresponding to 23.9 km s^{1}, and directed along the postion angle 348.3° ± 0.44°. In addition, radial velocity measurements found that LHS 1070 is approaching with 36.4 km s^{1} to the Sun (Basri & Marcy 1995). Due to its kinematic properties, LHS 1070 can be associated with the old disk population.
A first successful fit of the close orbit of component C around component B was presented by Leinert et al. (2001). Based on Speckle interferometric measurements a semimajor axis of a = 446 ± 29 mas and a period of P = 16.1 ± 1.4 yr were derived. When taking the refined distance of Costa et al. (2005) into account, the corresponding dynamical mass is 0.157 ± 0.042 M_{⊙}. A consistent mass of 0.157 ± 0.009 M_{⊙} was found by Seifahrt et al. (2008), who derived a wider (a = 461.9 ± 0.7 mas) orbit, but with P = 17.0 yr. Due to the better coverage the authors presented also a first reliable fit of the wide orbit of the components B and C around the primary. The dynamical mass of the whole system was determined as 0.272 ± 0.017 M_{⊙}. The orbits were found to be coplanar.
In this study, we present new measurements taken both with Speckle interferometric techniques and adaptive optics in the Kband (Sect. 2). The data is carefully selected and calibrated on a case by case basis with the aim to reach the highest degree of consistency possible. These new data are not only used to further refine the fits of the close (Sect. 3) and the wide orbit (Sect. 4). We utilise the primary component as astrometric reference, which enables the determination of the mass ratio of the components B and C. With this knowledge the individual masses of all three components can be derived. In Sect. 5 we discuss the implications of our findings with respect to the stability of the system and its evolutionary stage. We conclude in Sect. 6.
Fig. 1 Images of LHS 1070 obtained with NACO in December 2002, December 2004, October 2006, and November 2008. North is up, and East to the left. The separation between the two components B and C changes from about 330 mas in 2002 to about 240 mas in 2008 (cf. Table 1). 
Astrometric measurements of LHS 1070 BC.
2. Observations and data reduction
Our group monitored LHS 1070 since the discovery of its two companions, using a number of different telescopes and instruments. Table 1 gives a journal of observations and the measured relative positions. The following subsections describe the individual instruments and data reduction procedures.
2.1. ESO NTT/SHARP
The SHARP I camera (System for High Angular Resolution Pictures) of the MaxPlanckInstitut for Extraterrestrial Physics (Hofmann et al. 1992) was already used for many of the observations presented in Leinert et al. (2001). We used it again during an observing campaign in June and July 2001 at the European Southern Observatory (ESO) 3.5 m New Technology Telescope (NTT) on La Silla, Chile. LHS 1070 was observed in the Kband at 2.2 μm in four different nights. SHARP is a camera for speckle interferometry, which means that one observation consists of several hundred frames with short integration times (0.5 s). To reduce the data, we used our speckle program (Köhler et al. 2000). The same observing and data reduction strategy was already employed for the data published in Leinert et al. (2001) and a number of other multiplicity studies.
2.2. ESO 3.6 m/ADONIS/SHARP II
In January and December 2001, LHS 1070 was observed with the adaptive optics (AO) system ADONIS (Rousset & Beuzit 1999) and the SHARP II camera (Hofmann et al. 1995) at the ESO 3.6 m telescope on La Silla, Chile. We used the Kband filter of this instrument, which has a central wavelength of 2.177 μm. The observing strategy was “AOassisted speckle interferometry”, meaning we took many frames with short integration times. In January, we recorded 1000 frames with an integration time of 0.5 s each. In December, we took 240 frames with a longer integration time of 3 s each. For the data reduction, we used the same programs and algorithms as for the speckle data.
2.3. Calar Alto 3.5 m/Omega Cass
On several occasions between 2000 and 2006, LHS 1070 was observed with the Omega Cass camera (Lenzen et al. 1998) at the 3.5 mtelescope on Calar Alto, Spain. The camera is equipped with a 1024 × 1024pixel detector, but we used only a subarray of 128 × 128 pixels to enable the fast readout mode required for Speckleinterferometry. The observations were done in the Kband and with the highest resolution optics of the instrument (pixel scale ~95 mas/pixel). Our speckle program was used again to reduce the data.
2.4. ESO VLT/NACO
The majority of new observations were taken with NAOS/CONICA (NACO for short), the adaptive optics, nearinfrared camera at the ESO Very Large Telescope (VLT) on Cerro Paranal, Chile (Rousset et al. 2003; Lenzen et al. 2003). LHS 1070 was observed in the course of several programs (see Table 1). To ensure a consistent data set, only imaging observations in the K_{s} photometric band were used for the orbit determination.
The NACO images were sky subtracted with a median sky image, and bad pixels were replaced by the median of the closest good neighbors. Finally, the images were visually inspected for any artifacts or residuals. Figure 1 shows an example of the results. The Starfinder program (Diolaiti et al. 2000) was used to measure the positions of the stars. The positions in several images taken during one observation were averaged, and their standard deviation used to estimate the errors.
2.5. Plate scale and orientation
For a heterogenous data set like this, it is crucial to calibrate the absolute pixel scale and orientation of each observation. To this end, we took images of fields in the Orion Trapezium during each observing campaign, and reduced them in the same way as the images of the science targets. The measured positions of the cluster stars were compared with the coordinates given in McCaughrean & Stauffer (1994). The mean pixel scale and orientation were computed from a global fit of all star positions. The scatter of values derived from star pairs were used to estimate the errors (see Sect. 2.6).
The errors of the calibration are usually comparable to or larger than the errors of the measured positions of the science target, indicating the importance of a proper astrometric calibration. For this reason, we decided to use in this work only data where images of the Orion Trapezium were taken within a few days. We have two observations of LHS 1070 that were taken in June and July, when Orion was not observable. These observations are not used here, in order to ensure a consistent calibration of all the data. For the same reason, we do no use data from the literature that could not be recalibrated by our group.
2.6. Astrometric error estimates
Independent of instrument and telescope used, our observation resulted in at least four images or speckle cubes for each epoch. The final result was obtained by averaging the relative positions measured in individual images. As estimate for the errors, we use the standard deviation.
The error of the astrometric calibration was estimated in a similar way. We computed the pixel scale by dividing the separation of star pairs in arcseconds by the measured separation in pixels. Similarly, the orientation is the difference of the position angle on the sky and on the detector. The standard deviation of the results for the many star pairs in our calibration field was used as estimate for the error.
An additional source of error might be the motion of the telescope between the observation of the calibration field and the target. This should not cause much uncertainty for telescopes on a parallactic mount (i.e. the Calar Alto 3.5 m and the ESO 3.6 m). However, telescopes on an altazimuth mount (ESO NTT and VLT) have to correct for image rotation, which might introduce an error on the order of 0.1° (W. Brandner, priv. comm.). Additionally, the opening and closing of the AOloop, in connection with the active optics of the primary mirror, might induce small changes in the image scale. Our calibrations show that there are small changes in pixel scale and orientation on the time scale of weeks or months, but – to the best of our knowledge – there has been no study of shortterm variations.
This unknown error source might explain why we seem to underestimate our errors, which results in orbit fits with χ^{2} > 1 (cf. Tables 2 and 3).
3. The orbit of LHS 1070 BC
We estimated the orbital parameters of the BC pair by fitting orbit models to all observations listed in Table 1. We followed the procedure described in Köhler et al. (2008): a gridsearch in eccentricity e, period P, and time of periastron T_{0}. At each grid point, the ThieleInnes elements were determined by a linear fit to the observational data using Singular Value Decomposition. From the ThieleInnes elements, the semimajor axis a, the angle between node and periastron ω, the position angle of the line of nodes Ω, and the inclination i were computed.
Since the orbit of LHS 1070 BC is already quite wellknown (Leinert et al. 2001; Seifahrt et al. 2008), only a small range of parameter values had to be scanned: 200 points within 0.018 ≤ e ≤ 0.028, 200 points within 17.21 yr ≤ P < 17.31 yr, and initially 200 points for T_{0} distributed over one orbital period. After the initial scan over T_{0}, the best estimate for T_{0} was improved by rescanning a narrower range in T_{0} centered on the minimum found in the coarser scan. This grid refinement was repeated until the step size was less than one day.
Parameters of the best orbital solution for the pair BC.
We improved the results of the gridsearch with a LevenbergMarquardt χ^{2} minimization algorithm (Press et al. 1992) that fits for all 7 parameters simultaneously. The simpleminded approach would be to use the orbital elements with the minimum χ^{2} found with the gridsearch. However, initial test runs showed that the algorithm does not converge on the global minimum. For the same reason, we did not use one of the previously published orbit solutions as starting point (Leinert et al. 2001; Seifahrt et al. 2008). To make sure we find the globally minimum χ^{2}, we decided to use all orbits resulting from the gridsearch as starting points that had . The number 9 was chosen arbitrarily, to avoid starting from obviously bad orbits. The orbit with the globally minimum χ^{2} found in by the LevenbergMarquardt fit is shown in Fig. 2, and its elements are listed in Table 2. To convert the semimajor axis from mas to AU, we used the distance of 7.72 ± 0.15 pc (Costa et al. 2005).
Fig. 2 The orbit of component C around component B. The observed positions are marked by their error ellipses and lines connecting the observed and calculated position at the time of the observations. The observations with NACO are marked by crosses. Their errors are too small to be discernible. The dashdotted line indicates the line of nodes, the dashed line the periastron, and the arrow shows the direction of the orbital motion. 
The reduced χ^{2} of 3.5 is higher than expected for a good fit, which indicates that we underestimate our astrometric errors. The mean absolute difference between observed and predicted positions is also larger than typical errors of our measurements. This might be caused by unknown errors in the calibration, as mentioned in Sect. 2.6. However, the residuals show no pattern indicating a systematic error, and the errors given in Table 2 are our best guess. Therefore, we decided to accept the fit as it is.
Errors of the orbital elements were determined by studying the χ^{2} function around its minimum. Since we are interested in the confidence interval for each parameter taken separately, we have to perturb one parameter (for example T_{0}) away from the minimum, and optimize all the other parameters. Any perturbation of a parameter will of course lead to a larger χ^{2}. The range in T_{0} within which defines the 68% confidence interval for T_{0}. This interval is usually not symmetric around the T_{0} of the best fit, therefore we list in Table 2 separate limits for positive and negative perturbations. It should be noted that these limits describe the parameter range that contain 68% of the probability distribution, which is equivalent to the commonly used 1σerrors. However, the errors are not normally distributed, therefore a 2σ interval will not contain 95% of the probability distribution.
Parameters of the best solution for the orbit of BC around A.
Estimating the error of the mass required a special procedure. The mass itself is computed using Kepler’s third law (M = a^{3}/P^{2}). The semimajor axis a and the period P are usually strongly correlated. To obtain a realistic estimate for the mass error, we did not use the naive way of error propagation. Instead, we considered a set of orbital elements where the semimajor axis was replaced by the mass. This is possible because Kepler’s third law gives an unambiguous relation between the two sets of elements. With the mass being one of the orbital elements, we can treat it as one of the independent fit parameters and determine its error with the method described in the previous paragraph.
With an orbit derived from astrometric measurements, there always remains the ambiguity which of the two nodes is the ascending node (defined as the node where the companion is receding from the observer). Fortunately, Seifahrt et al. (2008) measured the relative radial velocities of LHS 1070 B and C on 2006 October 9. They found it to be negative, i.e. LHS 1070 C was approaching us^{1}. Therefore, the position angle of the ascending node is 14.65° (cf. Fig. 2), and the position angle of the descending node is 194.65°. Our orbit model predicts a relative radial velocity of −4.8 ± 0.1 km s^{1} on 2006 Oct. 9, in agreement with the measurement of Seifahrt et al. (2008).
4. The orbit of LHS 1070 BC around A
The orbit of components B and C around each other allows us to determine only the combined mass of B and C. To compute the individual masses, we need to know the mass ratio q, which can be computed if the position of the center of mass (CM) of B and C is known. Unfortunately, we cannot observe the CM directly. However, we know that the CM of B and C is in orbit around component A^{2}, and that B and C are in orbit around their CM. The CM is always on the line between B and C, and its distance from B is the constant fraction q/(1 + q) of the separation of B and C.
We follow the method that was used by Köhler et al. (2008) to derive masses in the triple system T Tauri. The position of the CM of B and C is described in two ways: first, it is on a Keplerian orbit around A, which is described by 7 orbital elements. Second, the position of the CM can be computed from the observed positions of B and C, and the mass ratio (which is treated as a free parameter). Standard error propagation is used to obtain an error estimate for this position. To compute χ^{2}, we compare the position of the CM from the orbit around A with the positions derived from the observations. Our model has therefore 8 free parameters, the 7 elements which describe the orbit of the CM of B+C around A, and the parameter f = q/(1 + q). The parameter f is often called fractional mass (Heintz 1978), since it is the secondary star’s fraction of the total mass in a binary. It is useful in our case because it also describes the fractional offset of the CM from B, i.e. the separation between B and the CM divided by the separation between B and C. For a gridsearch, f is better suited than q, because f is confined to the range 0 to 1, while q is a number between 0 and infinity.
The fitting procedure is similar to that used for the orbit of BC, except that the gridsearch is carried out in 4 dimensions: eccentricity e, period P, time of periastron T_{0}, and the fractional mass f. Singular Value Decomposition was used to fit the ThieleInnes constants, which give the remaining orbital elements. It is worth noting that the orbital elements in this fit describe the orbit of the ABC binary, only the fractional mass f refers to the pair B+C.
For the four dimensional grid search, we used 100 points each for f, e, and P. The grid ranged from 0 to 0.99 in f, 0 to 0.99 in e, and 30 to 300 years in P. The grid in T_{0} started with 100 points distributed uniformly over one orbital period. Similar to the fit for the orbit of BC, the grid in T_{0} was refined until the grid spacing was less than one day.
The results of the gridsearch were improved with a LevenbergMarquardt χ^{2} minimization algorithm (Press et al. 1992). As starting points, we used all orbits resulting from the gridsearch that had . The orbit with the globally minimum χ^{2} found by the LevenbergMarquardt fit is shown in Figs. 3 and 4. Table 3 lists its elements, as well as the mass ratio q, and the individual masses of all three components. As with the inner orbit of B and C, χ^{2} is somewhat larger then expected for a good fit.
Fig. 3 The orbit of the center of mass of components B and C around component A. The observed positions are marked by their error ellipses and lines connecting the observed and calculated position at the time of the observations. Note that these positions are computed from the observed positions of components B and C, and the mass ratio q, which is a free parameter of the model fit. The dashdotted line indicates the line of nodes, the dashed line the periastron, and the arrow shows the direction of the orbital motion. The two orbits shown are the model with the globally minimum χ^{2}, and an almost circular orbit with a longer period of about 80 years. 
Fig. 5 χ^{2} as function of P and e for the orbit of BC around A. The cross at P = 44.4 yr, e = 0.52 marks the minimum. The contour lines encircle the 68% confidence region (which is a series of mostly unconnected patches), the 95%, and the 99.7% confidence region (corresponding to 1σ, 2σ, and 3σ in the case of normally distributed errors). The dashed line marks the stability limit for the system, only orbits below this line are longterm stable (see Sect. 5.1). 
Errors for the parameters were again estimated by analyzing the χ^{2} function around the minimum. Since only a relatively small fraction of the orbit has been observed so far, the uncertainties for the orbital elements are much larger than for the orbit of B and C around each other. Figure 5 shows χ^{2} as function of period and eccentricity. Orbits with periods of more than 100 years are within the 95% confidence region, and periods of more than 200 years are in the 99.7% confidence region. The error of the system mass was derived in the same way as with the orbit of LHS 1070 B and C around each other, i.e. by treating M as an independent fit parameter.
The ambiguity of the ascending node is resolved with the help of Seifahrt et al. (2008). They measured in 2006 a positive relative radial velocity between the barycentre of LHS 1070 B and C and component A. Therefore, the binary BC was receding from the observer, and the position angle of the ascending node must be 26.8° (cf. Fig. 3). Our orbit model predicts a relative radial velocity of 3.1 ± 0.1 km s^{1} on 2006 Oct. 9, in agreement with the measurement of Seifahrt et al. (2008).
Finally, the mass of component A can be computed from the system mass of the triple and the mass of the binary B+C (derived in Sect. 3). The largest contribution to the error of both masses is the uncertainty in the distance to the system. Since this contribution to the errors of M_{ABC} and M_{BC} is correlated, it would not be correct to add the errors in quadrature. Instead, we computed masses from the semimajor axes in mas (which results in the unusual mass unit of mas^{3}/years^{2}), subtracted M_{BC} from M_{ABC}, and converted the resulting M_{A} into solar masses by multiplying it with the distance cubed. This way, the error of the distance enters the calculation only once, resulting in the correct estimate for the error of M_{A}.
5. Discussion
5.1. Is the LHS 1070 system stable?
Our best fit for the orbit of the binary BC around A is only a factor 2.5 larger than the orbit of BC itself. The outer orbit has also a rather large eccentricity. As a result, the distance from B or C to A becomes smaller than the distance to its binary companion during the periastron passage of the outer orbit (Fig. 6). This raises the question whether a triple system like this can be stable over timescales comparable to its age.
Fig. 6 Distances between the three stars in the system as function of time. Here we assume that the stars stay on the Keplerian orbits derived in Sects. 3 and 4, i.e. we do not simulate the threebody problem posed by the triple system. 
Fig. 7 Results of the fits for the outer orbit as function of their period. Shown are all model orbits within the 99% confidence region (3σ). Unstable orbits (according to the criterion by Donnison & Mikulskis 1995) are indicated by gray symbols, stable orbits are black. The orbit with the minimal χ^{2} is marked by an asterisk. 
Donnison & Mikulskis (1995) presented a criterion for the longterm stability of coplanar triple systems: The periastron distance q_{2} = a_{2}(1 − e_{2}) of the outer orbit should be at least a factor of 3.3 larger than the semimajor axis a_{1} of the inner orbit. In noncoplanar systems, increasing the inclination decreases the region of stability, i.e. the outer orbit has to be even larger to be stable (Donnison 2009). In Fig. 5, the limit for coplanar orbits (q_{2} > 3.3 a_{1}) is indicated by the dashed line, only orbits below the line are stable. The orbit with the minimum χ^{2} is clearly in the unstable regime. However, the 2σ and 3σ confidence regions extend well into the regime of stable orbits. It is therefore possible that the true orbit has a longer period of at least 80 years (cf. Fig. 3).
How does this uncertainty about the orbit influence our conclusions? Figure 7 shows the parameters of the outer orbit that we are most interested in. Shown are all orbits within the 99% confidence region in our grid of possible solutions. Unstable orbits (according to the criterion by Donnison & Mikulskis 1995) are indicated by gray symbols, stable orbits are black. The top panel shows the angle between the planes of the inner and outer orbits. This angle is 12.5° for the orbit with the minimum χ^{2}, and less than 10° for all stable orbits. We find a few orbital solutions with angles larger than 100°, but all these orbits are unstable. We therefore conclude that the inner and outer orbit are almost coplanar. Therefore, the stability criterion by Donnison & Mikulskis (1995) is sufficient in our case, although it does not apply to inclined orbits.
The second panel of Fig. 7 shows the masses of LHS 1070 A and B. Orbits with short periods result in unrealistically high masses for A, which also indicates that the true period is longer. The mass of component B lies between 0.076 and 0.077 M_{⊙}, almost independent of the orbital period. The mass of component C was not plotted, since it can easily be computed as 0.149 M_{⊙} − M_{B}.
Finally, the bottom panel of Fig. 7 shows the mass ratio between LHS 1070 C and B. All stable orbits result in a mass ratio of 0.923, which is the same as the mass ratio of the orbit with the minimum χ^{2}.
We conclude that it is possible to find stable orbital solutions that are compatible with the astrometric measurements. In these stable configurations, the orbit of B and C around each other, and the orbit of B+C around A are almost coplanar, with angles of less than 10°. The mass of component A lies between 0.13 and 0.16 M_{⊙} for stable orbits.
As a final note on the stability of the system, we point out that the stability criterion used here does not take into account resonances. For example, a period of the outer orbit of 52 years would be in a 3:1 resonance with the inner orbit. This orbit would be within the 1σ region around the orbit model with minimum χ^{2}. However, the uncertainties of our orbit fit are too large to provide a reasonable starting point for a search for stable resonant configurations. Furthermore, it would not significantly change our conclusions about coplanarity and the masses of the three components if the period of the outer orbit was between about 50 and 80 years.
5.2. Comparison with theoretical models
The goal of dynamical mass determinations is to test theoretical predictions for the mass and luminosity of the stars. Figure 8 shows a massluminositydiagram with the isochrones of Baraffe et al. (1998) and our results for the components of LHS 1070 (using the Vmagnitudes given by Leinert et al. 2000). Additionally, it shows the empirical MassLuminosityRelation by Henry et al. (1999), which is similar to the theoretical 500Myrisochrone.Within the error bars, both LHS 1070 B and C are compatible with an age between 500 and 800 Myr.
On the other hand, LHS 1070 A is either 3^{mag} too faint or 0.07 M_{⊙} too heavy compared to the theoretical models and the empirical relation. This is also true if we take into account that stable orbit configurations result in lower masses (indicated by the dotted line in Fig. 8). In the Kband, LHS 1070 A is about 1^{mag} too faint compared to the model isochrones. This excludes the possibility that the underluminosity is due to extinction by circumstellar or foreground material (as unlikely as it may be for a star of this age and distance), since then the effect in the infrared should be smaller, only about 0.3^{mag}.
Fig. 8 Massluminosity diagram with the components of LHS 1070, the models of Baraffe et al. (1998), and the empirical MassLuminosityRelation by Henry et al. (1999). Shown are the theoretical isochrones for ages between 500 Myr and 5 Gyr. The dotted horizontal line indicates the mass range for A resulting from stable orbit configurations. The points labeled “Aa photometry” and ”D photometry” indicate the magnitudes if A is split into two components (Henry et al. 1999), with the masses adjusted to put the points onto the 630Myr isochrone. 
The simplest explanation would be that A is indeed a binary as reported by Henry et al. (1999). In the following, we will call the components Aa and D, and continue to use A for the (unresolved) binary system. We computed the absolute Vmagnitudes of Aa and D from the combined apparent magnitude V = 15.35^{mag} (Leinert et al. 2000) and the magnitude difference ΔV = 2.46^{mag} (Henry et al. 1999), using the distance of 7.72 pc (Costa et al. 2005). By finding the intersection with the 630Myr isochrone from the models of Baraffe et al. (1998), we can estimate the mass of Aa and D to be 0.10 M_{⊙} and 0.08 M_{⊙}, resp. The combined mass of A resulting from photometry is therefore about 0.18 M_{⊙}, in reasonable agreement with the dynamical mass resulting from the orbit determinations.
In an attempt to find a possible close companion to LHS 1070 A, we analyzed all available highresolution data of the system. The resolution required to detect component D can only be reached with a 8 mclass telescope, preferentially at short wavelengths, e.g. Jband. This means we have to rely on NACO data. There are two datasets taken with NACO in the Jband on June 27. and December 12., 2003, and two datasets in the narrowband filter NB_2.17 in a mode that allows speckle postprocessing of the data. In all four datasets, we can exclude an companion with a flux ratio brighter than 0.2, i.e. a magnitude difference of 1.75^{mag}. The resolution of the four observations is 30–40 mas. With the masses of Aa and D derived above, the magnitude difference should be about 0.8^{mag}. Assuming a semimajor axis of 50 mas and a system mass of 0.17 M_{⊙} results in an orbital period of about 0.6 years. It is possible that all our observations were carried out at times when the companion was too close to the primary to be visible. However, given the large number of attempt to detect LHS 1070 D this appears unlikely.
We also looked into the possibility to detect a periodic astrometric signal of LHS 1070 D in the residuals of the ABC separation. If LHS 1070 A is an unresolved binary, then an astrometric shift of its center of light with the period of its orbit can occur. The magnitude of this shift depends on the position of the center of mass and the position of the center of light, hence the mass ratio and the flux ratio. The mass ratio of our suspected AaD binary is 0.8, while the flux ratio computed from the models of Baraffe et al. (1998) is 0.48. The expected motion of the center of light is therefore about 0.12 times the size of the orbit of D around Aa (6 mas for an orbit of 50 mas). Unfortunately, our sampling is not very dense, with 38 observations in 15.25 years. The Nyquist frequency ν_{N} = N/(2T) corresponds to a period of 0.8 years. Since our sampling is rather uneven, the Nyquist frequency is not a sharp limit for the detectability of a periodic signal. Nonetheless, due to the sparse sampling, it is unlikely that we can detect periods much shorter than 0.8 years. In fact, neither classical nor generalized LombScargleperiodograms (Scargle 1982; Zechmeister & Kürster 2009) show any significant peaks. Even adding an artificial signal with a period of 0.6 years and a large amplitude of 100 mas does not result in a significant detection, although peaks at 1/2 and 1/4 of the frequency appear. In summary, we conclude that an astrometric detection of LHS 1070 D is not possible with the available data.
Another possible explanation could be that LHS 1070 A is a binary too close to be resolved by the observations available so far. This hypothetical companion would not be identical with component D. However, at the moment the most likely explanation might be that we overestimate the system mass. This would not be too surprising, given the small observational coverage of the orbit.
6. Summary and conclusions
We present new relative positions of LHS 1070 A, B, and C, collected in the years 2000 to 2008. They were used to derive an improved model for the orbit of B and C around each other, and an estimate for the orbit of B and C around A. The orbit of B and C is welldetermined by now, with an orbital period of 17.24 ± 0.01 years and a system mass of 0.149 ± 0.009 M_{⊙} (cf. Table 2).
The observations of LHS 1070 B relative to A span only a range in position angle of 36°. The orbit of B and C around A derived from these observations is therefore not so wellconstrained. The orbital solution minimizing χ^{2} has a rather short period of about 44 years, but it results in an unstable configuration of the triple system. It is more likely that the true orbital period is in the range 80–200 years. Orbital solutions in this range of periods are still within the 95% and 99.7% confidence region for the fit to the astrometric observations (which corresponds to 2σ and 3σ for the case of normally distributed measurement errors). Despite this uncertainly, we can constrain the system mass M_{A} + M_{B} + M_{C} to 0.28 to 0.31 M_{⊙} if we accept only stable orbit configurations.
The outer orbit also yields the mass ratio of B and C, which is quite wellconstrained to 0.92 ± 0.01. Taken together, these results yield individual masses of M_{A} = 0.13...0.16 M_{⊙}, M_{B} = 0.077 ± 0.005 M_{⊙} and M_{C} = 0.071 ± 0.004 M_{⊙}.
Placing the three stars in a massluminosity diagram and comparing with theoretical isochrones shows that B and C are coeval within the measurement errors, with an age between 500 and 800 Myr. On the other hand, LHS 1070 A appears to be too faint for its mass, or too massive for its luminosity. One possible explanation could be that LHS 1070 A is itself a binary. There has been one report of the discovery of a close companion (Henry et al. 1999), but this discovery could not be confirmed, despite the large number of observations collected. The detection might have been caused by a glitch in the data (Henry, priv. comm.). To confirm or disprove the binarity of LHS 1070 A, observations sensitive to companions at very small separations would be helpful, either spectroscopic or interferometric. For the time being, the most likely explanation might be that our orbit fit overestimates the system mass.
Our results for the masses of the three stars show that LHS 1070 C is almost certainly a brown dwarf. LHS 1070 B is very close to the hydrogenburning masslimit, possibly also a brown dwarf. While LHS 1070 A is clearly above the hydrogenburning limit, it is still a very low mass star. The very low mass triple system LHS 1070 thus has the very interesting property to contain one component above, one below, and one at the hydrogen burning limit. This combination makes it an ideal candidate for testing the change of atmosphere properties in this mass regime. Furthermore, it allows some conclusions about its formation. One theory for the formation of brown dwarfs is that they got ejected out of a multiple system before they could accrete enough mass to start hydrogen burning (Reipurth & Clarke 2001). It is difficult how a system like LHS 1070 can survive such an event. Even if it managed to remain bound, one would expect more eccentric and inclined orbits than found in this work. It is more likely that LHS 1070 formed in the same way as normal stars, where a number of triple and higherorder multiple systems are known.
It is not entirely clear whether Seifahrt et al. (2008) give the velocity of C minus B, but they report Ω = 14.5°, which is consistent with our interpretation
We did not detect the companion D reported by Henry et al. (1999). Therefore, by “A” we mean the suspected close binary composed of A and D (if it exists) in the nomenclature of Henry et al. (1999).
Acknowledgments
We thank the many observers and support astronomers who carried out the many separate observations necessary for this work. We are also grateful for the helpful comments by the editor.
References
 Almeida, L. A., Jablonski, F., & Martioli, E. 2011, A&A, 525, A84 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Baraffe, I., Chabrier, G., Allard, F., & Hauschildt, P. H. 1998, A&A, 337, 403 [NASA ADS] [Google Scholar]
 Basri, G., & Marcy, G. W. 1995, AJ, 109, 762 [NASA ADS] [CrossRef] [Google Scholar]
 Costa, E., Méndez, R. A., Jao, W.C., et al. 2005, AJ, 130, 337 [NASA ADS] [CrossRef] [Google Scholar]
 Diolaiti, E., Bendinelli, O., Bonaccini, D., et al. 2000, A&AS, 147, 335 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Donnison, J. R. 2009, Planet. Space Sci., 57, 771 [NASA ADS] [CrossRef] [Google Scholar]
 Donnison, J. R., & Mikulskis, D. F. 1995, MNRAS, 272, 1 [NASA ADS] [Google Scholar]
 Heintz, W. D. 1978, Double Stars, Geophysics and Astrophysics Monographs No. 15 (Dordrecht, Holland: D. Reidel Publishing Company) [Google Scholar]
 Henry, T. J., Franz, O. G., Wasserman, L. H., et al. 1999, ApJ, 512, 864 [NASA ADS] [CrossRef] [Google Scholar]
 Hofmann, R., Blietz, M., Duhoux, P., et al. 1992, in Progress in Telescope and Instrumentation Technologies, ESO Conference and Workshop Proc. No. 42, ed. M.H. Ulrich (ESO Garching), 617 [Google Scholar]
 Hofmann, R., Brandl, B., Eckart, A., Eisenhauer, F., & TacconiGarman, L. E. 1995, in SPIE 2475, ed. A. M. Fowler, 192 [Google Scholar]
 Köhler, R., Kunkel, M., Leinert, C., & Zinnecker, H. 2000, A&A, 356, 541 [NASA ADS] [Google Scholar]
 Köhler, R., Ratzka, T., Herbst, T. M., & Kasper, M. 2008, A&A, 482, 929 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Leinert, C., Weitzel, N., Richichi, A., Eckart, A., & TacconiGarman, L. E. 1994, A&A, 291, L47 [NASA ADS] [Google Scholar]
 Leinert, C., Allard, F., Richichi, A., & Hauschildt, P. H. 2000, A&A, 353, 691 [NASA ADS] [Google Scholar]
 Leinert, C., Jahreiß, H., Woitas, J., et al. 2001, A&A, 367, 183 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Lenzen, R., Bizenberger, P., Salm, N., & Storz, C. 1998, in SPIE Conf. Ser. 3354, ed. A. M. Fowler, 493 [Google Scholar]
 Lenzen, R., Hartung, M., Brandner, W., et al. 2003, in Instrument Design and Performance for Optical/Infrared Groundbased Telescopes, ed. M. Iye, & A. F. M. Moorwood, SPIE Proc., 4841, 944 [Google Scholar]
 McCaughrean, M. J., & Stauffer, J. R. 1994, AJ, 108, 1382 [NASA ADS] [CrossRef] [Google Scholar]
 Press, W. H., Teukolsky, S. A., Vetterling, W. T., & Flannery, B. P. 1992, Numerical Recipes in C, 2nd edn. (Cambridge, UK: Cambridge University Press) [Google Scholar]
 Reiners, A., Seifahrt, A., Käufl, H. U., Siebenmorgen, R., & Smette, A. 2007, A&A, 471, L5 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Reipurth, B., & Clarke, C. 2001, AJ, 122, 432 [NASA ADS] [CrossRef] [Google Scholar]
 Rousset, G., & Beuzit, J.L. 1999, The COMEON/ADONIS systems, ed. F. Roddier, 171 [Google Scholar]
 Rousset, G., Lacombe, F., Puget, P., et al. 2003, in Adaptive Optical System Technologies II, ed. P. L. Wizinowich, & D. Bonaccini, SPIE Proc. 4839, 140 [Google Scholar]
 Scargle, J. D. 1982, ApJ, 263, 835 [NASA ADS] [CrossRef] [Google Scholar]
 Seifahrt, A., Röll, T., Neuhäuser, R., et al. 2008, A&A, 484, 429 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 van Altena, W. F., Lee, J. T., & Hoffleit, E. D. 1995, The general catalogue of trigonometric [stellar] parallaxes, ed. W. F. van Altena, J. T. Lee, & E. D. Hoffleit [Google Scholar]
 Zechmeister, M., & Kürster, M. 2009, A&A, 496, 577 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
All Tables
All Figures
Fig. 1 Images of LHS 1070 obtained with NACO in December 2002, December 2004, October 2006, and November 2008. North is up, and East to the left. The separation between the two components B and C changes from about 330 mas in 2002 to about 240 mas in 2008 (cf. Table 1). 

In the text 
Fig. 2 The orbit of component C around component B. The observed positions are marked by their error ellipses and lines connecting the observed and calculated position at the time of the observations. The observations with NACO are marked by crosses. Their errors are too small to be discernible. The dashdotted line indicates the line of nodes, the dashed line the periastron, and the arrow shows the direction of the orbital motion. 

In the text 
Fig. 3 The orbit of the center of mass of components B and C around component A. The observed positions are marked by their error ellipses and lines connecting the observed and calculated position at the time of the observations. Note that these positions are computed from the observed positions of components B and C, and the mass ratio q, which is a free parameter of the model fit. The dashdotted line indicates the line of nodes, the dashed line the periastron, and the arrow shows the direction of the orbital motion. The two orbits shown are the model with the globally minimum χ^{2}, and an almost circular orbit with a longer period of about 80 years. 

In the text 
Fig. 4 Enlarged section of the orbits in Fig. 3, showing the part covered by observations. 

In the text 
Fig. 5 χ^{2} as function of P and e for the orbit of BC around A. The cross at P = 44.4 yr, e = 0.52 marks the minimum. The contour lines encircle the 68% confidence region (which is a series of mostly unconnected patches), the 95%, and the 99.7% confidence region (corresponding to 1σ, 2σ, and 3σ in the case of normally distributed errors). The dashed line marks the stability limit for the system, only orbits below this line are longterm stable (see Sect. 5.1). 

In the text 
Fig. 6 Distances between the three stars in the system as function of time. Here we assume that the stars stay on the Keplerian orbits derived in Sects. 3 and 4, i.e. we do not simulate the threebody problem posed by the triple system. 

In the text 
Fig. 7 Results of the fits for the outer orbit as function of their period. Shown are all model orbits within the 99% confidence region (3σ). Unstable orbits (according to the criterion by Donnison & Mikulskis 1995) are indicated by gray symbols, stable orbits are black. The orbit with the minimal χ^{2} is marked by an asterisk. 

In the text 
Fig. 8 Massluminosity diagram with the components of LHS 1070, the models of Baraffe et al. (1998), and the empirical MassLuminosityRelation by Henry et al. (1999). Shown are the theoretical isochrones for ages between 500 Myr and 5 Gyr. The dotted horizontal line indicates the mass range for A resulting from stable orbit configurations. The points labeled “Aa photometry” and ”D photometry” indicate the magnitudes if A is split into two components (Henry et al. 1999), with the masses adjusted to put the points onto the 630Myr isochrone. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext 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 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.