Äîêóìåíò âçÿò èç êýøà ïîèñêîâîé ìàøèíû. Àäðåñ îðèãèíàëüíîãî äîêóìåíòà : http://www.sao.ru/hq/balega/PUBL_BAL/PUB_2008/aa_478_187.pdf Äàòà èçìåíåíèÿ: Wed Feb 24 16:01:53 2010 Äàòà èíäåêñèðîâàíèÿ: Tue Oct 2 11:27:26 2012 Êîäèðîâêà: Ïîèñêîâûå ñëîâà: arp 220

A&A 478, 187­191 (2008) DOI: 10.1051/0004-6361:20078594
c ESO 2008

Astronomy & Astrophysics

A methodology for studying physical and dynamical properties of multiple stars. Application to the system of red dwarfs Gl 22
J. A. Docobo1 , V. S. Tamazian1 ,Y. Y.Balega2 , M. Andrade1 , D. Schertl3 , G. Weigelt3 ,P.Campo1 , and M. Palacios4
1

2 3 4

Observatorio AstronÑmico RamÑn MarÌa Aller, Universidade de Santiago de Compostela, Avenida das Ciencias s/n., Santiago de Compostela, Spain e-mail: [oadoco;oatamaz;oandrade;oapedro]@usc.es Special Astrophysical Observatory, N. Arkhyz, Karachai-Cherkasia 369167, Russia e-mail: balega@sao.ru Max-Planck-Institut fÝr Radioastronomie, Auf dem HÝgel 69, 53121 Bonn, Germany e-mail: [ds;weigelt]@mpifr-bonn.mpg.de Grupo de MecÀnica Espacial, Universidad de Zaragoza, Spain e-mail: mpala@unizar.es

Received 31 August 2007 / Accepted 30 October 2007
ABSTRACT

Aims. The main aim of this study is the elaboration of a methodology for studying physical and dynamical properties of multiple stars and its application to the hierarchical triple system Gl 22. A careful evaluation of the components' masses and a study of the system's overall stability and long-term dynamical evolution were also pursued. Methods. New NIR speckle interferometric observations with the 6 m telescope of the Special Astrophysical Observatory (Russia) in the K photometric band have been carried out. We have made use of the method for orbit calculation reported by Docobo (1985). Results. An original methodology was elaborated and applied to evaluate the most probable elements of the outer orbit. Due to the almost definitive orbit of the inner pair, which just covered a full revolution, the motion of B relative to MC A has been carefully calculated. The position of MC A was estimated on the basis of differential photometry and empirical mass-luminosity relationships. A weak sinusoidal pattern in the apparent motion of the component B was noticed. Conclusions. Our methodology was successfully applied to the triple system Gl 22. The newly calculated outer orbit exhibits a moderate eccentricity (e = 0.29), which differs from the previously known circular solutions. Both orbits are coplanar and corevolving. This already known suggestion is now based on a much larger set of observational data, including a significant number of speckle measurements. Gl 22 is most likely a dynamically stable system, at least on the time scale of 10 Myr. The sinusoidal pattern in the motion of the B component could be caused by a fourth, unseen, very low-mass object with a mass of 0.015 M (16 MJ ) on a circular orbit around B with a period of 15 yr and semimajor axis 0. 35.
Key words. astrometry ­ stars: binaries: visual ­ stars: individual: Gl 22 ­ stars: fundamental parameters ­ stars: late-type ­ stellar

dynamics

1. Introduction
Visual companion B, at a distance of 2. 3 from the nearby bright star Gl 22 (HIP 2552; ADS 440), was discovered by Espin & Milburn (1926). Using the photographic plates obtained at Leander McCormick Observatory, Alden (1947) reported the primary component A to be an astrometric binary (hereafter components Aa and Ab) itself. On the basis of red/infrared spectra (6300 to 9000 å), spectral types M2.0V and M3.0V were assigned to the combined spectrum of the Aa+Ab and B components, respectively (Kirkpatrick et al. 1991). Throughout this paper, pairs Aa-Ab (with its mass center MC A ) and MC A - B will be referred to as the inner and the outer orbits, respectively (with the corresponding definition of the inner and outer orbits). Hershey (1973) published first orbits for the outer (a = 4. 1, P = 320 yr, e = 0.05) and inner ( = 0. 125, P = 15.95 yr, e = 0.0) pairs and reached the conclusion that the two orbits are coplanar and the stars are co-revolving. At the same time, he obtained estimates for the total mass of the Aa-Ab system (0.53 M ) and the B component (0.12 M ).

Prior to this, on the basis of a photographic plate collection obtained at Sproul Observatory, Heintz (1955) measured an absolute parallax of 0. 102 for Gl 22. Furthermore, with a larger set of plates, Heintz (1993) recalculated the orbits of Hershey, but their circular characters remained the same. Using near-infrared (NIR) speckle interferometry techniques, McCarthy et al. (1991) presented the first two directly resolved observations of the inner pair. The improved inner and outer orbits, both being circular, were obtained by SÆderhjelm (1999) and Lampens & Strigachev (2001), respectively. With our most recent measurements presented in this paper, the inner pair has now covered a full revolution. They demonstrate very small residuals with regard to our most recently calculated orbit (Docobo & Tamazian 2006) which is likely to be close to definitive. While the outer pair covered an arc of about 80 , it seemed reasonable to calculate its most probable orbital elements using an original methodology, described below. Taking into account new orbital elements of both pairs, we studied the system's dynamical stability by the application of several well-known stability criteria. Furthermore, we analysed

Article published by EDP Sciences and available at http://www.aanda.org or http://dx.doi.org/10.1051/0004-6361:20078594


188

J. A. Docobo et al.: A methodology for studying multiple stars

3.1. The actual orbit of the inner pair

Ab Aa

Our first orbit for the Aa-Ab pair was reported by Docobo et al. (2002) and then improved by Docobo & Tamazian (2006, but also see Docobo et al. 2006). The most recent speckle observations performed with the 6 m telescope in the K band led to observational residuals of +0 6 and -0. 003 with regard to the . improved orbit. It is worth noting that with this measurement, a full revolution has now been completed (see Fig. 2).
3.2. Is it possible to calculate the outer orbit with a reasonable accuracy?

B
Fig. 1. Diffraction-limited image of Gl 22 in the K band (the components Aa, Ab and B are marked). Field of view 5. 5, north is up, and east is to the left. Table 1. Astrometry and differential photometry in the K band. Intens. ratio (mas) (deg) Ab/Aa 2005.803 451 ± 9 45.5 ± 1.5 0.143 ± 0.005 B/(Aa+Ab) 2005.803 3954 ± 20 173.7 ± 0.5 0.235 ± 0.005 Comp. Epoch

the long-term evolution of the orbital elements on a time scale of 10 Myr, and also its long-term stability.

Regarding the outer orbit, to date, the secondary covers an arc of 80 with respect to MCA . However, we will calculate its most probable preliminary orbital elements according to the methodology that follows. First, luminosities are obtained on the basis of differential photometry in the K band (Woitas et al. 2003), to which the latest differential photometry data from our Table 1 were added. The empirical masses of each component are then derived on the basis of infrared mass-luminosity relations given by Delfosse et al. (2000) and Henry & McCarthy (1993) by using weighted means of the above-mentioned differential photometry data. The obtained estimates are given in Table 2. Furthermore, we obtain the most probable parallax by imposing the following condition: the empirical mass sum of the outer pair ( MAa+Ab = 0.515 ± 0.031 M ) must be equal to the dynamical mass determined through a direct application of Kepler's third law. This condition leads to a parallax of 0. 102 ± 0.001, which is exactly the same as that obtained by Heintz on the basis of the Sproul plates in 1955. It coincides with the Hipparcos parallax (0. 09874 ± 0. 00337) within a 1 margin.
3.3. New determination of the outer orbit

2. Observations
The speckle interferograms of Gl 22 were recorded on 21 October 2005 with the 6 m telescope of the Special Astrophysical Observatory (SAO, Russia) with exposure times of 590 ms each. The detector of the speckle camera was a Rockwell HAWAII array (only one 512 â 512 quadrant was used) with a pixel size corresponding to 28.7 mas on the sky. A K -band filter with a central wavelength of 2142 nm and a bandwidth of 293 nm was used. The direct image is shown in Fig. 1. The data set for Gl 22 consists of 523 speckle interferograms. To compensate the atmospheric speckle transfer function, 568 speckle interferograms of the unresolved star 2MASS 00312573+6738205 were used. The seeing (FWHM)was 3. 0. The modulus of the Fourier transform of the object (visibility) was obtained with the speckle interferometry method (Labeyrie 1970). Images with a diffraction-limited resolution of 75 mas were reconstructed from the data set using the bispectrum speckle interferometry method (Weigelt 1977; Weigelt & Wirnitzer 1983; Lohmann et al. 1983; Hofmann & Weigelt 1986). The epoch of the observation in fractional Besselian year, the measured separations, position angles, and intensity ratios of the components are listed in Table 1. Although there is a wealth of visual, photographic and CCD measurements of the outer pair, it must be kept in mind that they give the position of B with respect to the light centre of the inner pair Aa-Ab (LC A ), whereas speckle data give the position of B directly with respect to Aa. Therefore, all measurements must be transformed to positions with respect to MCA . In our case, this is possible due to the knowledge of 1) the brightness and mass of either component of the inner pair; 2) the precise inner orbit, allowing us to determine the position of each component with respect to both the MCA and LCA at any given time. Once these transformations bring the measurements to mutually comparable positions, we can proceed with the orbit calculation by using Docobo's method (Docobo 1985). Since the parallax is also known, for each generated orbit, we calculate the corresponding total mass. The latter will be a key differentiating factor for selecting the best orbit, since the rms of the O-C residuals in the and of generated orbits vary within a narrow margin (see Table 5) and thus do not allow us to choose the orbit with certainty. Therefore, we should search for an orbit with a dynamical mass close to the total mass of the system 0.692 ± 0.034 M given in Table 2. With this aim, let us select three base points (listed in Table 4) situated at the beginning, the middle, and the end of the measured arc. As a first approximation, let us give a step of 20 to the variable V (the difference between the eccentric anomalies corresponding to the base points 3 and 1) appearing in Table 5. It is seen that, indeed, the generated orbits have very similar rms of residuals but that the total mass varies in a relatively wide range.

3. The orbits
The set of available NIR speckle measurements of both pairs has been greatly widened due to observations with the 3.5 m telescope reported by Woitas et al. (2003), where data on differential photometry in the K band were also given.


J. A. Docobo et al.: A methodology for studying multiple stars

189

0.4

Gl 22 Aa

__ 1989.778 __ 2005.803

0.2

0

-0.2

E
-0.4

N
-0.6 -0.4 -0.2 0 0.2 0.4 0.6

Fig. 2. The apparent orbit of the Aa-Ab pair (the scale on both axes is in arcseconds). Each measurement (represented by a dot) is connected to its predicted position by an O-C line, and the first and last measurements are marked to show the orbit's closure. The dashed line passing through the primary star is the line of nodes, and the arrow shows the direction of the orbital motion.

Table 2. Luminosity and mass estimates. Component Aa Ab B Total mass M
K

Table 4. The base points. No. 1 2 3 ( ) 98.05 154.24 174.28 () 2.603 4.000 4.010 t 1924.250 1976.706 2005.8029

6.42 8.44 7.88

Mass (M ) 0.377 ± 0.030 0.138 ± 0.007 0.177 ± 0.014 0.692 ± 0.034

Table 3. Orbital elements and system mass. Orbital element P (yr) T e a( ) i ( ) ( ) ( ) Outer MC A - B pair 223.
+7 2 3-9..8 5.3 2.4 + 1859.4- + 0.293- + 3.322-

Inner Aa- Ab pair 15.64 ± 0.20 2000.76 ± 0.20 0.174 ± 0.003 0.511 ± 0.005 44.6 ± 1.5 106.8 ± 5.0 175.1 ± 1.0

Apart from having a moderate eccentricity (contrary to all previous solutions), the orbit has an inclination and node very similar to those of the inner pair. Assigning indices 1 and 2 to the parameters of the inner and outer orbits, respectively and hereafter in this paper, the mutual inclination can readily be obtained through a well-known expression cos I = cos i1 cos i2 + sin i1 sin i2 cos(2 - 1 ) yielding I = 2 7, which confirms that the orbits, in perfect ac. cordance with Hershey's suggestion, are practically coplanar. Finally, in Table 8 we present the ephemerides for the position of B relative to the MCA and LCA , and its position relative to the main component Aa. Notice that the positions predicted in Table 8 regarding LCA are useful to test those measurements that cannot resolve the inner pair.

0.044 0.025 0.040 0.060 0.5 0.3 2.0 3.8 2.7 1.3

+ 47.3-

+ 146.3- + 174.9-

4. Stability and dynamical evolution of the system
4.1. Stability criteria for hierarchical triple systems

Yet, we notice that the masses we are interested in appear in the range 90 < V < 110 . Therefore, we select this range of V but lower its step to 2 in the next approximation. After a second approximation, we find a total mass of 0.69 M for V = 105 deg, which corresponds to an orbit with a period of 223.4 yr. Due to the total mass uncertainty, its values between 0.66 and 0.73 M are also valid. Hence, the orbits leading to the masses in this range should be taken into account in order to calculate the standard errors of the orbital elements. A further improvement of the rms in and leads us to the final orbit with a dynamical mass of 0.693+0..069 M , shown in -0 062 Fig. 3 (its elements are listed in Table 3).

The motion in this hierarchical three-body system is a combination of two perturbed motions: the inner pair Aa-Ab and the B component with regard to MCA . In order to determine whether or not such a system is stable for a given set of initial orbital elements, several numerical criteria have been proposed by Harrington (1972, hereafter HA), Graziani & Black (1981, hereafter GB), Eggleton & Kiseleva (1995, hereafter EK), and Aarseth & Mardling (1999, hereafter AM). Both stability and critical parameters are calculated for each of the criteria mentioned above. The system is stable when its stability parameter is larger than the critical one (or smaller for


190 Table 5. Approximations. V ( ) 90 110 130 103 104 105 106 107
3

J. A. Docobo et al.: A methodology for studying multiple stars

P (yr) 336.9 198.0 130.7 234. 229. 223. 217. 212. 9 0 4 9 7

T

e

a ()

i ( )

( ) 162.0 2.1 31.9 172. 173. 174. 176. 177. 2 4 7 0 4

( ) 292.1 316.4 278.7 150. 148. 146. 144. 142. 0 1 2 3 4

rms

rms

Total mass (M ) 0.60 0.77 10.05 0. 0. 0. 0. 0. 66 67 69 70 71

1928.8 1870.7 1902.1 1854. 1857. 1859. 1861. 1864. 77 17 54 86 15

First approximation 0.046 4.175 52.2 0.411 3.181 46.9 0.914 5.669 71.2 Second approximation 0.249 3.387 47.8 0.271 3.347 47.5 0.294 3.310 47.3 0.317 3.277 47.1 0.340 3.247 47.0

0.96 1.32 1.32 0. 0. 0. 0. 0. 96 97 97 97 97

0.11 0.15 0.14 0. 0. 0. 0. 0. 11 11 11 11 10

Gl 22

WDS 00321+6715

H T

2

Micrometric measurement Photographic " Hipparcos " Tycho " CCD " Speckle (LC) " Speckle " Three-body model orbit Four-body " "

Wobble motion

1

H T

H T
0

-1

E
-2 -4

N
-3 -2 -1 0 1 2 3

Fig. 3. Left: the apparent orbit of component B relative to the mass centre of Aa­Ab (the scale on both axes is in arcseconds). Each measurement (see legend) is connected to its predicted position by an O­C line. The dashed line passing through the primary star is the line of nodes, and the arrow shows the direction of the orbital motion. Right: zoom of the orbit showing the wobble motion of the B component due to the possible new object in the area of measurements with the highest weight. Table 7. The range of periodic variations of orbital elements at the end of 10 Myr.

Table 6. Stability margin. HA GB EK AM Values -0.151 +0.405 -0.240 +0.226 Stability ? yes ? yes

the GB criterion). Otherwise, nothing can be stated regarding its stability. Thus, the stability margin of a triple system can be defined (Orlov & Petrova 2000) as p - pc p = ± , pc where p represents the stability parameter and pc the critical parameter ("+" for HA, EK, and AM criteria; "-" for GB criterion). The larger p , the larger the stability margin of the system will be according to the given criterion. Thus, if p < 0, then the system could be dynamically unstable, i.e. a violation of hierarchy or an escape of the distant component would be possible. In principle, the small ratio of the (outer to inner) semimajor axis (equal to 6.48) is a sign of the system's instability. With the orbital elements given in Table 3, the stability margin for each criterion is calculated, and results are given in Table 6. According to the GB and AM criteria, the system would be stable, whereas those of HA and EK do not allow one to reach the same conclusion. In any case, the obtained values are small enough to clearly decide whether or not the system will be stable. Indeed, it is very close to the stability boundary.

Orbital element P (yr) T e a (AU) i ( ) ( ) ( )

Outer pair 4.7 secular 0.023 0.457 1.2 secular 1.6

Inner pair 0.46 secular 0.095 0.097 4.5 secular 6.1

4.2. Stability and long-term evolution

To investigate the long-term evolution of the system in a precise manner, a more detailed analysis is necessary. To integrate the hierarchical (nearly-Keplerian) three-body system, we used an implicit Runge-Kutta integrator. All integrations conserve energy and angular momentum to better than 10-9 and 10-4 , respectively. We explore the orbital elements and parameter space in 10 Myr integrations to demonstrate that neither semimajor axes nor eccentricities exhibit secular variations. Indeed, as can be seen in Table 7, all orbital elements undergo relatively small periodical variations, except for arguments of periastron that advance secularly (1 and 2 cover a full revolution in about 104 and 5 â 104 yr respectively) and times of periastron passage. Moreover, the mutual inclination remains approximately constant with a maximum periodical variation as small as I = 0 2 for the considered interval. .


J. A. Docobo et al.: A methodology for studying multiple stars Table 8. Ephemerides. Epoch 2006. 2008. 2010. 2012. 2014. 2016. 0 0 0 0 0 0 B - MC ( ) 174.2 175.5 177.2 179.0 181.1 183.0
A

191

orbit ( ) 3.983 3.991 3.991 3.967 3.905 3.800

B- ( ) 174.9 176.3 177.8 179.3 180.8 182.4

LC A ( ) 4.016 3.982 3.944 3.900 3.852 3.798

B-Aa ( ) ( ) 173.5 3.949 174.7 4.001 176.5 4.042 178.8 4.038 181.4 3.961 183.6 3.802

5. The possible existence of a fourth body
From a visual inspection of Fig. 3, a weak sinusoidal pattern in the apparent motion of component B has been noticed. This pattern can be attributed to either a very unusual distribution of observational residuals or an unseen fourth body in the system. The latter situation is supported by the fact that such apparent motion has been obtained by considering a collection of different techniques, as can be seen in Fig. 3. In this case, star B would consist of the components Ba and Bb, whose mass centre will hereafter be defined as MC B . Assuming Bb to be a very low-mass object of 0.015 M (16 M J ) on a circular orbit around Ba with a period of 15 yr, semimajor axis 0. 35, and coplanar with two other orbits, the observational residuals of the outer orbit are improved. In such a case, the component Ba would be moving relative to the MCB on an orbit with a semimajor axis of 0. 03. The blue line in Fig. 3 corresponds to the motion of B relative to MCA (or that of MCB relative to MCA , if Bb really does exist), whereas the red line shows the motion of Ba affected by the putative component Bb. Similar to all visual, photographic, and CCD observations, a single speckle measurement marked as "speckle (LC)" was initially showing the position of B relative to the LCA . For the calculation of orbits, we reduced all such measurements to MCA .

like Gl 22, with an almost definitive inner orbit and a reasonable arc covered by the wide pair. We have also carried out numerical studies of the stability of the Gl 22 system. The application of four well-known empirical criteria allowed us to conclude that Gl 22 is very close to the boundary of stability. However, the application of an implicit Runge-Kutta integrator extended over 10 Myr showed no significant secular changes in major semi-axes, eccentricities or inclinations. A weak sinusoidal pattern in the apparent motion of component B has been noticed. If we assume the existence of a very low-mass object "Bb" with a mass of 0.015 M (16 M J ) on a circular orbit with a period of 15 yr, the observational residuals of the outer orbit are improved. If confirmed, this would be the first example of a "speckle-astrometric" binary ever discovered. However, a further careful investigation with more extended astrometric and radial velocity data is needed to verify this trend.
Acknowledgements. This paper was financed by research project AYA200407003 and PGIDIT06 PXIB243031PR of Spanish Ministerio de EducaciÑn y Ciencia and Xunta de Galicia respectively. This research made use of the Washington Double Star Catalogs maintained at the US Naval Observatory and SIDONIe Database maintained at the Nice-CÒte d'Azur Observatory (France).

References
Aarseth, S. J., & Mardling, R. A. 1999, in The Dynamics of Small Bodies in the Solar System, a Major Key to Solar System Studies, ed. B. A. Steves, & A. E. Roy (Dordrecht: Kluwer Academic Publishers), 385 Alden, H. L. 1947, AJ, 52, 138 Delfosse, X., Forveille, T., SÈgransan, D., et al. 2000, A&A, 364, 217 Docobo, J. A. 1985, Celest. Mech., 36, 143 Docobo, J. A., & Tamazian, V. S. 2006, IAU Comm. 26, 160 Docobo, J. A., Tamazian, V. S., Woitas, J., & Leinert, Ch. 2002, IAU Comm. 26 Inf. Circ., 147 Docobo, J. A., Tamazian, V. S., Balega, Y. Y., & Melikian, N. D. 2006, AJ, 132, 994 Eggleton, P., & Kiseleva, L. 1995, ApJ, 455, 640 Espin, T. E., & Milburn, W. 1926, MNRAS, 86, 131 Graziani, F., & Black, D. C. 1981, ApJ, 251, 337 Harrington, R. S. 1972, Celest. Mech., 6, 322 Hartkopf, W. I., & Mason, B. D. 2007, Sixth Catalog of Orbits of Visual Binary Stars, http://ad.usno.navy.mil/wds/orb6.html Heintz, W. D. 1955, AN, 282, 221 Heintz, W. D. 1993, PASP, 105, 44 Henry, T. J., & McCarthy, D. W. 1993, AJ, 106, 773 Hershey, J. L. 1973, AJ, 78, 935 Hofmann, K.-H., Koch, R., & Weigelt, G. 1986, Mitt. Astron. Gesellschaft, 67, 430 Kirkpatrick, J. D., Henry, T. J., & McCarthy, D. W., Jr. 1991, ApJS, 77, 417 Labeyrie, A. 1970, A&A, 6, 85 Lampens, P., & Strigachev, A. 2001, A&A, 368, 572 Lohmann, A. W., Weigelt, G., & Wirnitzer, B. 1983, App. Opt., 22, 4028 McCarthy, Jr., D. W., Henry, T. J., McLeod, B., & Christou, J. C. 1991, AJ, 101, 214 Orlov, V. V., & Petrova, A. V. 2000, Astron. Lett., 26, 250 Pettersen, B. R. 1975, A&A, 41, 113 SÆderhjelm, S. 1999, A&A, 341, 121 Weigelt, G. 1977, Opt. Comm., 21, 59 Weigelt, G., & Wirnitzer, B. 1983, Opt. Lett., 8, 389 Woitas, J., Tamazian, V. S., Docobo, J. A., & Leinert, Ch. 2003, A&A, 406, 293

6. Conclusions
A methodology for studying the physical and dynamical properties of multiple stars is presented. On this basis, an attempt is made to obtain maximal information on the astrophysical and dynamical properties of the hierarchical triple system of red dwarfs Gl 22. Recent speckle measurements and differential photometry in the K band allowed us to calculate a better orbital solution for the outer orbit in the hierarchical triple system of the red dwarfs Gl 22 and to confirm that the orbit of the inner pair Aa-Ab is almost definitive. The obtained orbital elements, parallax, and individual masses represent a robust system of compatible parameters for this interesting system of low-mass, main-sequence stars exhibiting a low flare activity. Our results confirm Hershey's suggestion that two orbits are coplanar and co-revolving (Hershey 1973). A newly developed methodology can be applied to the study of triple systems