Free Access
Issue
A&A
Volume 607, November 2017
Article Number A10
Number of page(s) 9
Section Stellar structure and evolution
DOI https://doi.org/10.1051/0004-6361/201730641
Published online 30 October 2017

© ESO, 2017

1. Introduction

Precise astrometry of the orbital motion of stars in binary systems allows us to obtain accurate, model-independent estimates of the masses of the individual components, relevant to provide tests of stellar evolution models (Guirado et al. 1997; Torres et al. 2012; Reid & Honma 2014). In the case of low- and very-low-mass stars, only a handful of stellar systems with dynamically-determined masses have been reported (Hillenbrand & White 2004; Stassun et al. 2004; Mathieu et al. 2007). In this context, the pre-main-sequence (PMS) stellar system AB Doradus, the main star of the AB Doradus moving group (AB Dor-MG), is of particular interest. AB Doradus is placed at a distance of ~15 pc and it is actually a quadruple system formed by two pairs of stars separated by 9′′, AB Dor A/C and AB Dor Ba/Bb (Close et al. 2005; Guirado et al. 2006).

AB Dor A (=HD 36705 =HIP 25647) is the main star of the system. Given its high-level of radio emission, during the early 90’s AB Dor A was monitored astrometrically with VLBI as part of a program to link the optical reference frame, defined by the positions of the stars observed by the ESA’s astrometric satellite Hipparcos (Lindegren & Kovalevsky 1995; Lestrade et al. 1995), and the celestial reference frame, defined by VLBI positions (Lestrade et al. 1995). The VLBI-Hipparcos combined data revealed the presence of a low-mass companion, 0.090 M, AB Dor C (Guirado et al. 1997), orbiting AB Dor A at an average angular distance of 0.2′′. Imaging this new low-mass companion was difficult and it had to wait until the VLT near-infrared observations carried out by Close et al. (2005, 2007) using commissioning time of the Simultaneous Differential Imager (SDI) at the VLT, which allowed independent photometry of AB Dor C, free from contamination of the much brighter AB Dor A.

The combination of precise measurements of the dynamical mass and photometry of AB Dor C provided relevant consequences in terms of calibration of the mass-luminosity relationship for young, low-mass objects: basically, stellar model predictions tend to underestimate the mass of the observed objects. Given the need for an accurate calibration of the models, any possible ambiguity or uncertainty affecting these few calibrating objects should be studied in detail. A vivid discussion about the nature, age, and characterization of this low-mass object followed from these results, including new observations made using the VLT chronograph (Boccaletti et al. 2008), VLT SINFONI (Thatte et al. 2007), Spitzer (Plavchan et al. 2009), and VLTI (Guirado et al. 2011).

Table 1

Journal of observations of AB Dor A.

thumbnail Fig. 1

Clean maps of the star AB Dor A at the three LBA epochs. In each map, the lowest contour level corresponds to three times the statistical root-mean-square (0.20, 0.05, and 0.18 mJy beam-1) with a scale factor between contiguous contours of 2\hbox{$\sqrt{2}$}. The peak flux densities in the images are, respectively, 11.1, 5.84, and 10.10 mJy beam-1. For image parameters see Table 2. In all the maps we have centered the position of AB Dor A at the origin. Absolute positions are presented in Table 3.

In this work we present a deeper analysis of the VLBI observations of the AB Doradus system reported in Azulay et al. (2015), based on a new correlation pass centered on AB Dor A/C. Our analysis improves the determination of the orbital parameters of the system and obtains precise values of the dynamical masses of both components. We also discuss the comparisons of the dynamical masses of this pair with PMS stellar evolution models.

Table 2

Circular Gaussian fits corresponding to the VLBI maps of the components of AB Dor A.

Table 3

Compilation of astrometric measurements for the AB Dor A system.

Table 4

Estimates of the astrometric and orbital parameters of AB Dor A.

2. Observations and data reduction

AB Dor A was observed in three epochs using the Long Baseline Array (LBA), the Australian VLBI Network, at 8.4 GHz (see Table 1). The observations were scheduled in phase-reference mode, using the ICRF-defining source BL Lac PKS 0516621 (about 3.6° away) as the phase calibrator. The sequence target-calibrator-target lasted about four minutes. Results from these observations were reported in Azulay et al. (2015) for AB Dor B, the companion to AB Dor A, that is 10′′ away from it. Since both sources lie within the primary beam of all telescopes, the AB Dor A data were obtained from a new correlator pass using the nominal AB Dor A position as phase center.

The data calibration was carried out using the program Astronomical Image Processing System (AIPS) of the National Radio Astronomy Observatory (NRAO) with standard routines. We followed the process already described in Azulay et al (2015), and outlined briefly here: (i) calibration of the visibility amplitudes using the gain curves of each antenna and the system temperatures of each station; (ii) correction of both the parallactic angle and the ionospheric contribution; and (iii) fringe-search on the calibrator PKS 0516621 to remove residual contributions to the phases, whose corrections were applied to the target source AB Dor A. The phase-referenced naturally-weighted images of AB Dor A were obtained using the Caltech imaging program DIFMAP (Shepherd et al. 1994) and are shown in Fig. 1.

thumbnail Fig. 2

Absolute orbits of the components AB Dor A and AB Dor C using the orbital elements in Table 4 (with aA and aC, respectively). Positions of the component A (triangles, Hipparcos data; circles, VLBI data) and C (squares) are marked. The center of mass of the system is placed at the origin.

thumbnail Fig. 3

Postfit residuals in right ascension (upper panel) and declination (lower panel) of the component AB Dor A. The weighted rms is 1.1 mas (see text).

thumbnail Fig. 4

Results of our MCMC exploration of the orbital-parameter space. All axis ranges are normalized to the values given in Table 4. Bottom left: cross-correlations between parameters. Top right: posterior probability distributions. Dotted lines mark the ranges used to derive the averages (and standard deviations) of the parameters, given in Table 4. We note the secondary peak in T0, which is separated from the first peak by the orbital period.

thumbnail Fig. 5

Comparison of AB Dor A with some PMS theoretical models (BCAH98, top left; MDKH04, top right; TDP12, bottom left; BHAC15, bottom right). Isochrones (dashed lines) correspond to 10 (top isochrone), 16, 25, 40, 50, and 100 Myr. Isomasses (solid lines) are for 0.75 (bottom isomass), 0.80, 0.85, 0.90, 0.95, and 1.00 M in the case of BCAH98 and TDP12, for 0.80 (bottom isomass), 0.85, 0.90, 0.95, and 1.00 M in the case of MDKH04, and for 0.80 (bottom isomass), 0.90, and 1.00 M in the case of BHAC15. We have considered the value αML = 1.90 for BCAH98 and the value αML = 1.74 for TDP12.

thumbnail Fig. 6

Comparison of AB Dor A with a set of PMS evolutionary tracks with [Fe/H] = 0 for different initial helium abundances (standard abundance, left panels; low abundance, right panels), and for different mixing length values αML (1.74, top panels; 1.00, bottom panels). The shaded area corresponds to the region between track with masses in 0.820.98 M, as given by the uncertainty on dynamical mass.

thumbnail Fig. 7

As Fig. 6 for [Fe/H] = +0.1.

3. Results

3.1. Maps of the AB Dor A/C system

The star AB Dor A is clearly detected in all three maps presented in Fig. 1. Although AB Dor A is known to be a strong radio emitter (White et al. 1988) and other VLBI observations have been previously reported (Guirado et al. 1997), these are (to our knowledge) the first VLBI images published of this object. To model these VLBI images, we made least-squares fits to the interferometric visibilities by using circular Gaussians (i.e., model fitting in the uv-plane; see Table 2). We notice that the use of elliptical Gaussians is not necessary to model the structure of AB Dor A, adding unnecessary free parameters to the analysis.

From the flux density displayed in Table 2 and the parallax given in Table 4, we calculated a mean value of the absolute radio luminosity for AB Dor A of LR = 3.16 × 1015 erg Hz-1 s-1, very similar to other radio stars of the AB Dor-MG. This value indicates that the radio emission from this star arises from non-thermal, mildly relativistic electrons, or non-thermal gyrosynchtron emission.

The companion AB Dor C, directly imaged by Close et al. (2005) with the NACO SDI camera, was not detected in any of the three epochs, with a noise floor of 0.20, 0.05, and 0.18 mJy, respectively. We made a specific search for the emission of AB Dor C by stacking the individual images of the system. Taking advantage of the predicted position of AB Dor C with respect to AB Dor A at all observing epochs (Guirado et al. 2006; see Sect. 3.2), we shifted each of the images so that the expected position of component C was at the origin of coordinates. In practice, we made a stacked image by convolving the three images with a common restoring beam (6.4 × 5.4 mas with a PA of 90°, which is nearly that of epoch 2010.816, the largest beam given in Table 1). Although there was no detection of AB Dor C, we set, thus, an upper bound to the radio emission of this object of 0.07 mJy.

3.2. Orbital parameters of AB Dor A/C

The new detections of AB Dor A in our VLBI observations have allowed us to determine new absolute positions of the star (i.e., positions of the star with respect to the ICRF defining source the PKS 0516621), which are shown in Table 3 along with other absolute and relative positions (AB Dor A AB Dor C) previously reported. Our new VLBI positions of AB Dor A necessarily brought us to a re-calculation of the orbital elements and masses of this system.

To carry out this revision, we used two fitting methods, a least-squares analysis and a Bayesian approach, that we describe in turn. First, we followed the least-squares fit method detailed in Azulay et al. (2015). We combined the absolute positions of AB Dor A and all the relative positions available to determine the astrometric and orbital parameters which are shown in Table 4. As stated in Azulay et al. (2015), a priori values of the orbital elements are essential to find convergence of the least-squares method; we used the orbital elements given in Guirado et al. (2006) as a priori values for the least-squares analysis, which greatly facilitated the convergence of our procedure. Diagrams of the absolute orbits are shown in Fig. 2. The weighted rms of the postfit residuals of AB Dor A (plotted in Fig. 3) is 1.1 mas; although some unmodelled effects are present, the shown uncertainties are sufficient to exclude the presence of further companions to this system with masses below 0.15 MJ at a distance of 1 AU from the central star.

In complement to the above approach, we used a Bayesian approach based on a Markov-chain Monte Carlo (MCMC) routine, which in turn uses the Metropolis-Hastings algorithm (Robert 2015). We explored the solution space corresponding to the range of parameters given in Table 4. The orbital parameters obtained agree within uncertainties with those resulting from the least-squares method, which are already reported in Table 4. We show in Fig. 4 the posterior probability distributions for all parameters, together with their cross-correlations. The coincidence of both methods shows the robustness of our orbital solution for the AB Dor A/C system.

thumbnail Fig. 8

Comparison of AB Dor C with PMS theoretical models of Chabrier et al. (2000) (top) and Baraffe et al. (2015) (bottom). Isochrones (dashed lines) correspond to 10 (top isochrone), 50, and 100 Myr, in the case of the DUSTY models, and to 10 (top isochrone), 16, 25, 40, 50, and 100 Myr, in the case of BHAC15. Isomasses (solid lines), in both cases, are for 0.04 (rightmost isomass), 0.05, 0.07, 0.09, and 0.10 M.

From the values of the mass ratio mC/mA, the period P (in years), and the parallax π (in arcsec) in Table 4, we could derive the sum of the masses of the system using the Kepler’s third law mA+mC=[1+(mCmA)-1]3(aA/π)3P2,\begin{eqnarray*} m_{\mathrm{A}} + m_{\mathrm{C}} = \bigg[1+\bigg(\frac{m_{\mathrm{C}}}{m_{\mathrm{A}}}\bigg)^{-1}\bigg]^{3} \frac{(a_{\mathrm{A}}/\pi)^{3}}{P^{2}} \mbox{,} \end{eqnarray*}which is 0.98 ± 0.08 M. Moreover, with the semimajor axis of the absolute orbit of AB Dor A, aA, and using the Kepler’s third law in the form mC3(mA+mC)2=(aA/π)3P2,\begin{eqnarray*} \frac{m_{\mathrm{C}}^{3}}{(m_{\mathrm{A}} + m_{\mathrm{C}})^{2}} = \frac{(a_{\mathrm{A}} / \pi)^{3}}{P^{2}} \mbox{,} \end{eqnarray*}we derived a value of 0.090 ± 0.008 M for the mass of AB Dor C. Since we did not detect the component AB Dor C in our VLBI observations, we could not determine the semimajor axis of the absolute orbit of this component, aC and, consequently, the mass of the component A; so the mass of AB Dor A follows from subtracting the value of the sum of the mass of both components and the value of the mass of AB Dor C. We obtained mA = 0.89 ± 0.08 M.

4. Discussion

The dynamical determination of the individual masses of AB Dor A and AB Dor C is an important result that can be used to test theoretical predicitons of PMS stellar evolution models. Both AB Dor A and C can be placed in H-R diagrams using the measured luminosity and the effective temperature, along with theoretical isomasses and isochrones. Different theoretical models are used depending on their availability according to the masses of the objects.

Regarding AB Dor A, in Fig. 5 we compared the measurements of the effective temperature and luminosity (5081 ± 50 K and 0.388 ± 0.008 L, respectively; Close et al. 2007) with the predictions of the PMS models of Baraffe et al. (1998; BCAH98), Montalbán et al. (2004; MDKH04), Tognelli et al. (2011, 2012; TDP12), and Baraffe et al. (2015; BHAC15). We refer to the cited papers for a complete discussion of the physics and calculations involved. The two first model comparisons (BCAH98 and MDKH04) are reported in Guirado et al. (2011) and are shown here again for completeness along with the other, newer models. Considering the isomasses displayed in Fig. 5, we see a general good agreement between the predicted and the dynamical mass, although a persistent underprediction is seen at all models (between 2% and 6%). Nevertheless, these differences are covered by the uncertainties shown for the measurements of effective temperature and luminosity. Regarding the isochrones, the models shown in Fig. 5 favor an age for AB Dor A of 4050 Myr, very similar to that predicted in Guirado et al. (2011).

We took advantage of the flexibility of the TDP12 models to investigate how theoretical predictions are affected by the adopted values of Y (initial helium mass fractional abundance) and Z (initial total metallicity) by showing the effect of varying these parameters on the tracks position in the H-R diagram. The quantities Y and Z are not directly obtained from the observations, but they are derived from the measured [Fe/H]. Therefore, we computed models for [Fe/H] = +0.0 and +0.1, the extreme values of the reported estimates of [Fe/H] for the AB Dor-MG, namely, [Fe/H] = 0.02 ± 0.02 (Barenfeld et al. 2013), 0.10 ± 0.03 (Biazzo et al. 2012), and 0.04 ± 0.05 (Viana Almeida et al. 2009). We analyzed the effect of the models of adopting different initial helium content in the star, by producing a set of stellar tracks with a low initial helium abundance. The results of the computations with the two [Fe/H] and Y values are plotted in Figs. 6 and 7, for stellar tracks in the mass range 0.700.95 M. For each figure, models with two values of the mixing length parameter, namely αML = 1.74 (solar calibrated mixing length parameter, top panels) and αML = 1.00 (bottom panel) are shown. We overplotted the 16, 25, 40, 50, and 100 Myr isochrones.

thumbnail Fig. 9

Comparison of AB Dor C with a set of PMS evolutionary tracks with [Fe/H] = 0 for different initial helium abundances (standard abundance, left panels; low abundance, right panels), and for different mixing length values αML (1.74, top panels; 1.00, bottom panels). The shaded area corresponds to the region between track with masses in 0.0820.099 M, as given by the uncertainty on dynamical mass.

thumbnail Fig. 10

As Fig. 9 for [Fe/H] = +0.1.

By comparing the models shown in panels of Figs. 6 and 7 it emerges that the variation of [Fe/H] produces a relevant effect on theoretical models. Indeed, an increase of 0.1 dex in [Fe/H] produces an increase of about 20% in Z, which significantly alters the position of the tracks in the HR-diagram making the stars cooler and fainter. Consequently, the derived mass gets larger if [Fe/H] = + 0.1 models are adopted. Moreover, at a fixed [Fe/H], a reduction of the initial helium content produces cooler and fainter stars. Even in this case the derived mass slightly increases, partially reducing the discrepancy between the measured and inferred stellar mass. Concerning the mixing length parameter, from the comparisons it seems to emerge that TDP12 solar calibrated models, i.e., αML = 1.74 are probably preferred (similarly to what found by Guirado et al. 2011). However, we remark that the uncertainties on both the luminosity and temperature do not allow us to reject the αML = 1.00 models.

From the comparisons shown in Figs. 6 and 7, we can conclude that AB Dor A is better reproduced by the set of values [Fe/H] = + 0.1 and αML = 1.74, with a derived mass of ~0.85 M and an age of about 4050 Myr. However, more precise measurements of the total luminosity and the effective temperature are required to better constrain the models.

Regarding AB Dor C, the mass derived for this object (0.090 ± 0.008 M) does not differ substantially from that reported in Guirado et al. (2006), therefore the comparisons reported in Close et al. (2007) and Boccaletti et al. (2008) still hold. For completeness we include here the H-R diagram (Fig. 8) using the recent BHAC15 models (along with those from Chabrier et al. 2000; DUSTY models). We used the effective temperature from Close et al. (2007) and the luminosity derived using the Ks magnitude given by Boccaletti et al. (2008) and derived using the procedure (distance and bolometric correction) given in Close et al. (2007). The models predict a compatible estimate with our dynamical mass, although slightly underestimating the result (by approximately 10%; although this comparison is limited by the uncertainties in the effective temperature). Both models suggest an age between 25 and 120 Myr, consistent with the values published in Close et al. (2007) and Boccaletti et al. (2008).

We performed the analysis using also the TDP12 models. Similarly to what done for the AB Dor A component, we compared the data with models computed with the two values of [Fe/H], namely + 0.0 and + 0.1, with the standard and low initial helium content. Even in this case we used models with two mixing length parameter values, namely 1.74 and 1.00. The results are shown in Figs. 9 and 10, were we plotted the tracks in the 0.050.10 M mass range and the isochrones of 25, 40, 50, 100 and 200 Myr. To place the star in the H-R diagram, we also used the effective temperature given by Close et al. (2007) and the luminosity obtained from the Ks magnitude given by Boccaletti et al. (2008).

As clearly visible, the models are fully compatible, within the observational uncertainties, with the data. In this case, the tracks are only marginally affected by the adopted chemical composition and mixing length parameter. Consequently, given the relatively large uncertainties, we can not state which of the plotted sets gets the best agreement. From the comparisons, we derived a central mass of about 0.07 M and a central age of 4050 Myr. We also tried to compare the TDP12 models with the data using the luminosity given by Close et al. (2007), which is about 0.1 dex lower than that derived using the Boccaletti et al. (2008) Ks magnitude. In this case, we obtained a mass of approximately 0.08 M and an age of about 80 Myr.

5. Conclusions

We revisited the 8.4 GHz observations presented in Azulay et al. (2015) by re-correlating the observed data centered in the PMS stellar system AB Dor A/C, thus complementing the previous astronomic monitoring by Guirado et al. (1997, 2006). The main star AB Dor A was clearly detected in all three epochs, with a flux density varying between 5 and 11 mJy, meanwhile the low-mass component AB Dor C was not detected. We have set an upper bound to the radio emission of this low-mass object of 0.07 mJy.

With these new epochs of observations of AB Dor A, we re-analyzed the orbital motion of the system following two different approaches (based on least-squares fits and MCMC analysis, respectively) which combined absolute positions of AB Dor A (from VLBI and Hipparcos) and near-infrared relative positions available (from VLT/NACO). From the estimated orbital parameters we determined values of the masses of 0.89 ± 0.08 M and 0.090 ± 0.008 M for AB Dor A and AB Dor C, respectively.

We compared these dynamical masses with the predictions of different PMS stellar evolution models. In general, we found a good agreement between the predicted and the dynamical masses, although a slight underestimation (already present in previous works) is predicted by all models. This data-theory comparison is limited by the uncertainties in the determination of the effective temperature and/or the luminosity. We have shown that the theoretical predictions could also be affected by the adopted chemical composition, that with a small variation may cause large differences in the derived masses.

Regarding the age, the models considered in this paper favor an age of 40–50 Myr for AB Dor A and 25–120 Myr for AB Dor C. This model-dependent age is on the younger side of the age estimates of several authors: 40–50 Myr (Zuckerman et al. 2004; López-Santiago et al. 2006; Guirado et al. 2011), 50–100 Myr (Nielsen et al. 2005; Janson et al. 2007; Boccaletti et al. 2008). However, it is barely coincident with (especially in the case of AB Dor A) the lower limit of the AB Dor nucleus star (110 Myr) established by Barenfeld et al. (2013) or the recent isochronical age of the AB Dor moving group (130–200 Myr, Bell et al. 2015).

New, more sensitive VLBI observations will be needed to detect the companion AB Dor C. Although the orbital motion of the AB Doradus system is remarkable well determined, knowledge of the submas-precise radio position of AB Dor C will dramatically improve the precision of the estimated dynamical masses.

Acknowledgments

R.A., J.C.G., J.M.M, and E.R. were partially supported by the Spanish MINECO projects AYA2012-38491-C02-01 and AYA2015-63939-C2-2-P and by the Generalitat Valenciana project PROMETEO/2009/104 and PROMETEOII/2014/057. E.T. was supported by the “PRA 2016 Università di Pisa”. E.T. also acknowledges the INFN iniziativa specifica TAsP. The Long Baseline Array is part of the Australia Telescope National Facility which is funded by the Commonwealth of Australia for operation as a National Facility managed by CSIRO.

References

  1. Azulay, R., Guirado, J. C., Marcaide, J. M., Martí-Vidal, I., & Arroyo-Torres, B. 2014, A&A, 561, A38 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  2. Azulay, R., Guirado, J. C., Marcaide, J. M., et al. 2015, A&A, 578, A16 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  3. Baraffe, I., Chabrier, G., Allard, F., & Hauschildt, P. H. 1998, A&A, 337, 403 [NASA ADS] [Google Scholar]
  4. Baraffe, I., Homeier, D., Allard, F., & Chabrier, G. 2015, A&A, 577, A42 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  5. Barenfeld, S. A., Bubar, E. J., Mamajek, E. E., & Young, P. A. 2013, ApJ, 766, 6 [NASA ADS] [CrossRef] [Google Scholar]
  6. Bell, C. P. M., Mamajek, E. E., & Naylor, T. 2015, MNRAS, 454, 593 [NASA ADS] [CrossRef] [Google Scholar]
  7. Biazzo, K., D’Orazi, V., Desidera, S., et al. 2012, MNRAS, 427, 2905 [NASA ADS] [CrossRef] [Google Scholar]
  8. Boccaletti, A., Chauvin, G., Baudoz, P., & Beuzit, J.-L. 2008, A&A, 482, 939 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  9. Chabrier, G., Baraffe, I., Allard, F., & Hauschildt, P. 2000, ApJ, 542, 464 [NASA ADS] [CrossRef] [Google Scholar]
  10. Close, L. M., Lenzen, R., Guirado, J. C., et al. 2005, Nature, 433, 286 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  11. Close, L. M., Thatte, N., Nielsen, E. L., et al. 2007, ApJ, 665, 736 [NASA ADS] [CrossRef] [Google Scholar]
  12. Guirado, J. C., Reynolds, J. E., Lestrade, J.-F., et al. 1997, ApJ, 490, 835 [NASA ADS] [CrossRef] [Google Scholar]
  13. Guirado, J. C., Martí-Vidal, I., Marcaide, J. M., et al. 2006, A&A, 446, 733 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  14. Guirado, J. C., Marcaide, J. M., Martí-Vidal, I., et al. 2011, A&A, 533, A106 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  15. Hillenbrand, L. A., & White, R. J. 2004, ApJ, 604, 741 [NASA ADS] [CrossRef] [Google Scholar]
  16. Janson, M., Brandner, W., Lenzen, R., et al. 2007, A&A, 462, 615 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  17. Lestrade, J.-F., Jones, D. L., Preston, R. A., et al. 1995, A&A, 304, 182 [NASA ADS] [Google Scholar]
  18. Lindegren, L., & Kovalevsky, J. 1995, A&A, 304, 189 [NASA ADS] [Google Scholar]
  19. López-Santiago, J., Montes, D., Crespo-Chacón, I., & Fernández-Figueroa, M. J. 2006, ApJ, 643, 1160 [NASA ADS] [CrossRef] [Google Scholar]
  20. Martí-Vidal, I., Marcaide, J. M., Quirrenbach, A., et al. 2011, A&A, 529, A115 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  21. Mathieu, R. D., Baraffe, I., Simon, M., Stassun, K. G., & White, R. 2007, Protostars and Planets V, 411, eds. B. Reipurth, D. Jewitt, & K. Keil (Tucson: University of Arizona Press) [Google Scholar]
  22. Mérand, A., Bordé, P., & Coudé du Foresto, V. 2006, A&A, 447, 783 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  23. Montalbán, J., D’Antona, F., Kupka, F., & Heiter, U. 2004, A&A, 416, 1081 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  24. Nielsen, E. L., Close, L. M., Guirado, J. C., et al. 2005, Astron. Nachr., 326, 1033 [NASA ADS] [CrossRef] [Google Scholar]
  25. Plavchan, P., Werner, M. W., Chen, C. H., et al. 2009, ApJ, 698, 1068 [NASA ADS] [CrossRef] [Google Scholar]
  26. Reid, M. J., & Honma, M. 2014, ARA&A, 52, 339 [NASA ADS] [CrossRef] [Google Scholar]
  27. Robert, C. P. 2015, ArXiv e-prints [arXiv:1504.018963] [Google Scholar]
  28. Shepherd, M. C., Pearson, T. J., & Taylor, G. B. 1994, BAAS, 26, 987 [NASA ADS] [Google Scholar]
  29. Stassun, K. G., Mathieu, R. D., Vaz, L. P. R., Stroud, N., & Vrba, F. J. 2004, ApJS, 151, 357 [NASA ADS] [CrossRef] [Google Scholar]
  30. Tatulli, E., Millour, F., Chelli, A., et al. 2007, A&A, 464, 29 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  31. Thatte, N., Abuter, R., Tecza, M., et al. 2007, MNRAS, 378, 1229 [NASA ADS] [CrossRef] [Google Scholar]
  32. Tognelli, E., Prada Moroni, P. G., & Degl’Innocenti, S. 2011, A&A, 533, A109 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  33. Tognelli, E., Degl’Innocenti, S., & Pra da Moroni, P. G. 2012, A&A, 548, A41 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  34. Torres, R. M., Loinard, L., Mioduszewski, A. J., et al. 2012, ApJ, 747, 18 [NASA ADS] [CrossRef] [Google Scholar]
  35. Viana Almeida, P., Santos, N. C., Melo, C., et al. 2009, A&A, 501, 965 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  36. White, G. L., Jauncey, D. L., Batty, M. J., Peters, W. L., & Gulkis, S. 1988, PASP, 100, 825 [NASA ADS] [CrossRef] [Google Scholar]
  37. Zuckerman, B., Song, I., & Bessell, M. S. 2004, ApJ, 613, L65 [NASA ADS] [CrossRef] [Google Scholar]

All Tables

Table 1

Journal of observations of AB Dor A.

Table 2

Circular Gaussian fits corresponding to the VLBI maps of the components of AB Dor A.

Table 3

Compilation of astrometric measurements for the AB Dor A system.

Table 4

Estimates of the astrometric and orbital parameters of AB Dor A.

All Figures

thumbnail Fig. 1

Clean maps of the star AB Dor A at the three LBA epochs. In each map, the lowest contour level corresponds to three times the statistical root-mean-square (0.20, 0.05, and 0.18 mJy beam-1) with a scale factor between contiguous contours of 2\hbox{$\sqrt{2}$}. The peak flux densities in the images are, respectively, 11.1, 5.84, and 10.10 mJy beam-1. For image parameters see Table 2. In all the maps we have centered the position of AB Dor A at the origin. Absolute positions are presented in Table 3.

In the text
thumbnail Fig. 2

Absolute orbits of the components AB Dor A and AB Dor C using the orbital elements in Table 4 (with aA and aC, respectively). Positions of the component A (triangles, Hipparcos data; circles, VLBI data) and C (squares) are marked. The center of mass of the system is placed at the origin.

In the text
thumbnail Fig. 3

Postfit residuals in right ascension (upper panel) and declination (lower panel) of the component AB Dor A. The weighted rms is 1.1 mas (see text).

In the text
thumbnail Fig. 4

Results of our MCMC exploration of the orbital-parameter space. All axis ranges are normalized to the values given in Table 4. Bottom left: cross-correlations between parameters. Top right: posterior probability distributions. Dotted lines mark the ranges used to derive the averages (and standard deviations) of the parameters, given in Table 4. We note the secondary peak in T0, which is separated from the first peak by the orbital period.

In the text
thumbnail Fig. 5

Comparison of AB Dor A with some PMS theoretical models (BCAH98, top left; MDKH04, top right; TDP12, bottom left; BHAC15, bottom right). Isochrones (dashed lines) correspond to 10 (top isochrone), 16, 25, 40, 50, and 100 Myr. Isomasses (solid lines) are for 0.75 (bottom isomass), 0.80, 0.85, 0.90, 0.95, and 1.00 M in the case of BCAH98 and TDP12, for 0.80 (bottom isomass), 0.85, 0.90, 0.95, and 1.00 M in the case of MDKH04, and for 0.80 (bottom isomass), 0.90, and 1.00 M in the case of BHAC15. We have considered the value αML = 1.90 for BCAH98 and the value αML = 1.74 for TDP12.

In the text
thumbnail Fig. 6

Comparison of AB Dor A with a set of PMS evolutionary tracks with [Fe/H] = 0 for different initial helium abundances (standard abundance, left panels; low abundance, right panels), and for different mixing length values αML (1.74, top panels; 1.00, bottom panels). The shaded area corresponds to the region between track with masses in 0.820.98 M, as given by the uncertainty on dynamical mass.

In the text
thumbnail Fig. 7

As Fig. 6 for [Fe/H] = +0.1.

In the text
thumbnail Fig. 8

Comparison of AB Dor C with PMS theoretical models of Chabrier et al. (2000) (top) and Baraffe et al. (2015) (bottom). Isochrones (dashed lines) correspond to 10 (top isochrone), 50, and 100 Myr, in the case of the DUSTY models, and to 10 (top isochrone), 16, 25, 40, 50, and 100 Myr, in the case of BHAC15. Isomasses (solid lines), in both cases, are for 0.04 (rightmost isomass), 0.05, 0.07, 0.09, and 0.10 M.

In the text
thumbnail Fig. 9

Comparison of AB Dor C with a set of PMS evolutionary tracks with [Fe/H] = 0 for different initial helium abundances (standard abundance, left panels; low abundance, right panels), and for different mixing length values αML (1.74, top panels; 1.00, bottom panels). The shaded area corresponds to the region between track with masses in 0.0820.099 M, as given by the uncertainty on dynamical mass.

In the text
thumbnail Fig. 10

As Fig. 9 for [Fe/H] = +0.1.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.