Issue 
A&A
Volume 561, January 2014



Article Number  A101  
Number of page(s)  4  
Section  Stellar structure and evolution  
DOI  https://doi.org/10.1051/00046361/201322923  
Published online  14 January 2014 
Research Note
Refined masses and distance of the young binary Haro 114 C^{⋆}
^{1}
UJFGrenoble 1/CNRSINSU, Institut de Planétologie et d’Astrophysique de
Grenoble (IPAG) UMR 5274,
38041
Grenoble,
France
email:
jeanbaptiste.lebouquin@obs.ujfgrenoble.fr
^{2}
European Southern Observatory, 85748
Garching by München,
Germany
^{3}
Lowell Observatory, 1400 West Mars Hill Road, Flagstaff, AZ
86001,
USA
^{4}
The CHARA Array of Georgia State University, Mount Wilson
Observatory, Mount
Wilson, CA
91023,
USA
Received: 28 October 2013
Accepted: 13 December 2013
Aims. We aim to refine the dynamical masses of the individual component of the lowmass premain sequence binary Haro 114C.
Methods. We combine the data of the preliminary orbit presented previously with new interferometric observations obtained with the four 8 m telescopes of the Very Large Telescope Interferometer.
Results. The derived masses are M_{a} = 0.905 ± 0.043 M_{⊙} and M_{b} = 0.308 ± 0.011 M_{⊙} for the primary and secondary components, respectively. This is about five times better than the uncertainties of the preliminary orbit. Moreover, the possibility of larger masses is now securely discarded. The new dynamical distance, d = 96 ± 9 pc, is smaller than the distance to the Ophiuchus core with a significance of 2.6σ. Fitting the spectral energy distribution yields apparent diameters of φ_{a} = 0.13 ± 0.01 mas and φ_{b} = 0.10 ± 0.01 mas (corresponding to R_{a} = 1.50 R_{⊙} and R_{b} = 1.13 R_{⊙}) and a visual extinction of A_{v} ≈ 1.75. Although the revised orbit has a nearly edgeon geometry, the system is unlikely to be a longperiod eclipsing binary.
Conclusions. The secondary in Haro 114C is one of the few lowmass, premain sequence stars with an accurately determined dynamical mass and distance.
Key words: binaries: close / stars: formation / techniques: interferometric
© ESO, 2014
1. Introduction
Haro 114C (HBC 644) is one of the few lowmass, premain sequence binaries whose orbit can be observed both with radial velocity and resolved astrometric measurements. This provides a way of measuring the dynamical masses and distance and of testing the theoretical evolutionary models. Consequently, Haro 114C attracted the attention of observers over the past decade. Reipurth et al. (2002) first identified the system as a singlelined spectroscopic binary. Simon & Prato (2004) detected the radial velocity of the secondary and determined the full doublelined spectroscopic orbit. Finally, Schaefer et al. (2008) spatially resolved the system for the first time using optical longbaseline interferometry and derived a preliminary threedimensional orbit.
The apparent magnitudes of Haro 114C in the nearinfrared are 8.0 ± 0.05 mag and 7.78 ± 0.03 mag in the H and K bands, respectively (2MASS Catalog). This is too faint for the existing suit of instruments installed at the focus of imaging interferometers, which operate with 1 m class telescopes (such as CHARA or VLTIAT). Observations from Schaefer et al. (2008) were obtained with the two 10 m telescopes of the Keck Interferometer (KeckI). However, the single KeckI baseline provides poor sampling in (u,v) coverage. This is especially true for Haro 114C, given its declination. Consequently, the preliminary threedimensional orbit suffers from a long tail of possible solutions in the probability distribution toward smaller inclinations and toward larger masses for the components: , and pc.
Fig. 1
Left and middle panels: PIONIER square visibilities and closure phase (dots) plotted together with the bestfit binary model (solid lines). The colors represent the four different pointings detailed in Table 1. The rightmost panel shows the reduced χ^{2} of the model with respect to the position of the secondary. The red cross represents the astrometry of the bestfit model, while the open circle represents the prediction from the published orbit at the time of observation. The gray line indicates the trace of the orbit, the gray circle its periastron passage, and the gray dashed line shows the line of nodes. 
2. Observations
Supplementary interferometric data were obtained with the PIONIER combiner (Le Bouquin et al. 2011) at the Very Large Telescope Interferometer (VLTI, Haguenauer et al. 2010). The use of the four 8 m unit telescopes (instead of the four 1.8 m auxiliary telescopes) was necessary to observe Haro 114C (H = 8.0 mag). The visible magnitude (V = 12.3 mag) was enough for the MACAO adaptive optics system to provide an adequate correction during the observations. Data were dispersed over three spectral channels across the H band. We obtained four calibrated points during the night 20130619 (ESO convention). Telescope UT4 was missing during the first point because of technical problems. To calibrate the fringe visibilities and the closure phase, we interleaved observations of unresolved single stars between those of Haro 114C. These reference stars were found with the SearchCal^{1} software in its faint mode (Bonneau et al. 2011). Table 1 summarizes the log of observations. Data were reduced and calibrated with the pndrs package (Le Bouquin et al. 2011).
Log of observations with the Uniform Disk Diameter (UDD) of the reference stars.
The calibrated square visibilities and closure phases are shown in Fig. 1 together with the bestfit binary model. The expected diameters of the individual components, <0.15 mas, are unresolved by the longest VLTI baselines and were set to zero in the fitting procedure. Thanks to the simultaneous combination of four telescopes, the binary was detected without ambiguity in the relative position of the two components. We performed a bootstrap analysis by randomly selecting subsets of the observations and adding noise to these datasets. The uncertainty is given by the dispersion of the bestfit parameters when fitting these datasets. The bestfit position of the faint component with respect to the bright component is (δEast = −2.66 mas, δNorth = + 0.74 mas). The 3σ error ellipse is (0.53 mas, 0.22 mas) with a position angle for the major axis of 130deg, measured east from north. The average Modified Julian Day of observation is 56 463.104.
Fig. 2
Left: bestfit orbital solution compared with radial velocities of the primary (blue) and of the secondary (red). Right: bestfit orbital solution compared with the pseudoastrometric uncertainty regions from KeckI (colors) and VLTI (black, smaller than the symbol size). In both plots the black circles are the predictions from the bestfit orbital solution. 
Fig. 3
Probability distributions for the most relevant parameters (inclination, distance, and masses of the components). 
3. Combined orbit
At first consideration, the orbit published by Schaefer et al. (2008) is largely inconsistent with the PIONIER astrometry (reduced ). The observed astrometry is in the vicinity of the periastron, but the date of observation is far from the expected periastron date (MJD = 56 616.1), which is well constrained by the radial velocities of the primary.
However, the KeckI observations are insensitive to the absolute orientation because of the lack of closurephase measurements: the published orbit can be flipped by Ω′ = Ω + 180deg. Doing so yields a better agreement and a discrepancy now oriented along the large uncertainties of the preliminary orbit. Starting from this flipped orbit, we ran a combined fit on all available data: radial velocities, KeckI square visibilities, and PIONIER closure phases and square visibilities. All orbital elements were let free to vary, but this new fit mainly updates the inclination. The bestfit orbital elements are summarized in Table 2 and compared with the preliminary solution from Schaefer et al. (2008).
Bestfit orbital elements and related physical parameters.
We performed a simple bootstrap analysis by adding noise to the data (radial velocities, visibilities, and closure phases) according to the uncertainties, and fitting these noisy datasets. For illustration purposes, we computed the binary separations at the times of the interferometer observations to build the pseudoastrometric uncertainty regions shown in Fig. 2. We discarded the VLTI observation when building the Keck pseudoastrometry, and vice versa. The bestfit orbital solution and its uncertainties were computed with the entire dataset. The probability distribution for the most relevant parameters (inclination, distance, and masses) are shown in Fig. 3. The preliminary orbit suffers from a long tail in the probability distribution toward smaller inclinations as a direct consequence of the almost linear arrangement of all observed (u,v) points. This is not true anymore with the addition of the PIONIER observations: the probability distributions can be approximated by normal distributions. The uncertainties presented Table 2 are the standard deviation of the probability distribution of the corresponding parameter in the bootstrap analysis. The relative uncertainties on the stellar masses are now ± 5% and ± 3.5% for the primary and the secondary component, respectively. This is about five times better than the uncertainties of the preliminary orbit. Moreover, the tail toward larger masses in the probability distributions is now securely discarded.
The revised orbital solution yields a dynamical distance of d = 96 ± 9 pc. Interestingly, this value is lower than the distance to the Ophiuchus core (119 ± 4 pc, Loinard et al. 2008) with a significance of 2.6σ.
4. Discussion
The secondary in Haro 114C is one of the few lowmass premain sequence stars with an accurately determined of its dynamical mass and distance. The remaining main limiting factor in making comparisons with the evolutionary tracks is the estimation of the effective temperature. Given the binary separation, its period and total brightness, GAIA should be able to accurately measure the flux ratio in the visible. Combined with the flux ratio in the nearinfrared, this provides key information for separating the spectral energy distribution of the secondary, and possibly for estimating its effective temperature with respect to that of the primary. It will also be possible to compare the observations with models using the V − K colors of each component in place of effective temperatures. This approach is potentially more precise (Torres et al. 2013).
Another possible solution lies in a direct measurement of the component diameters. However, they are out of reach of current optical interferometers with hectometric baseline lengths. Given the high inclination of the revised orbit, one might wonder whether Haro 114C might be a long period eclipsing binary. To test this idea, we fitted the spectral energy distribution with the sum of two atmospheric models (see Fig. 4). We used effective temperatures of 4400 ± 200 K for the primary and 3500 ± 200 K for the secondary (based on the spectral types from Schaefer et al. 2008). The interferometric flux ratio in the nearinfrared corresponds to a ratio of 0.75 ± 0.1 for the linear radii. We then adjusted the photometry by varying the apparent diameters and the visual extinction A_{V}, assuming that the latter is the same for the two components. The W4 photometric point at 22 μm was discarded from the fit because it is a clear outlier in the RayleighJeans regime. The best match is obtained for A_{V} = 1.8 ± 0.7 and for φ_{a} = 0.131 ± 0.015 mas and φ_{b} = 0.098 ± 0.011 mas using ATLAS9 models (Castelli & Kurucz 2003). Differences using NextGen models (Allard et al. 1997) are negligible. Unfortunately, Haro 114C cannot be an eclipsing binary because the closest approach between the two components in the apparent orbit is 0.32 mas. These diameters correspond to R_{a} = 1.50 R_{⊙} and R_{b} = 1.13 R_{⊙} in linear radius at the distance of the object.
Fig. 4
Bestfit synthetic fluxes (horizontal bars) to the spectralenergy distribution from NOMAD, 2MASS, and WISE (red circles). The corresponding synthetic spectra for the two individual components and their sum are overlaid (dashed, dotted, and solid gray lines). 
Finally, it would be very interesting to confirm the distance discrepancy between Haro 114C and the core of the Ophiuchus cloud. The dynamical distance is the physical parameter that is the least constrained by the existing dataset. This is because the distance linearly depends on the size of the apparent orbit (a_{app}), which is still degenerate at the 10% level with the orientation on
sky (Ω). A few additional, highquality astrometric observations at various epochs are needed to lift this degeneracy.
Acknowledgments
PIONIER is funded by the Université Joseph Fourier (UJF), the Institut de Planétologie et d’Astrophysique de Grenoble (IPAG), the Agence Nationale pour la Recherche (ANR06BLAN0421 and ANR10BLAN0505), and the Institut National des Science de l’Univers (INSU PNP and PNPS). The integrated optics beam combiner is the result of a collaboration between IPAG and CEALETI based on CNES R&T funding. The authors warmly thank everyone involved in the VLTI project. This work is based on observations made with the ESO telescopes. It made use of the Smithsonian/NASA Astrophysics Data System (ADS) and of the Centre de Donnees astronomiques de Strasbourg (CDS). All calculations and graphics were performed with the freeware Yorick. The work of L.P. was supported by NSF grant AST1009136.
References
 Allard, F., Hauschildt, P. H., Alexander, D. R., & Starrfield, S. 1997, ARA&A, 35, 137 [NASA ADS] [CrossRef] [Google Scholar]
 Bonneau, D., Delfosse, X., Mourard, D., et al. 2011, A&A, 535, A53 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Castelli, F., & Kurucz, R. L. 2003, IAU Symp. No 210, poster A20 [Google Scholar]
 Haguenauer, P., Alonso, J., Bourget, P., et al. 2010, in SPIE Conf. Ser., 7734 [Google Scholar]
 Le Bouquin, J.B., Berger, J.P., Lazareff, B., et al. 2011, A&A, 535, A67 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Loinard, L., Torres, R. M., Mioduszewski, A. J., & Rodríguez, L. F. 2008, ApJ, 675, L29 [NASA ADS] [CrossRef] [Google Scholar]
 Reipurth, B., Lindgren, H., Mayor, M., Mermilliod, J.C., & Cramer, N. 2002, AJ, 124, 2813 [NASA ADS] [CrossRef] [Google Scholar]
 Schaefer, G. H., Simon, M., Prato, L., & Barman, T. 2008, AJ, 135, 1659 [NASA ADS] [CrossRef] [Google Scholar]
 Simon, M., & Prato, L. 2004, ApJ, 613, L69 [NASA ADS] [CrossRef] [Google Scholar]
 Torres, G., RuízRodríguez, D., Badenas, M., et al. 2013, ApJ, 773, 40 [NASA ADS] [CrossRef] [Google Scholar]
All Tables
Log of observations with the Uniform Disk Diameter (UDD) of the reference stars.
All Figures
Fig. 1
Left and middle panels: PIONIER square visibilities and closure phase (dots) plotted together with the bestfit binary model (solid lines). The colors represent the four different pointings detailed in Table 1. The rightmost panel shows the reduced χ^{2} of the model with respect to the position of the secondary. The red cross represents the astrometry of the bestfit model, while the open circle represents the prediction from the published orbit at the time of observation. The gray line indicates the trace of the orbit, the gray circle its periastron passage, and the gray dashed line shows the line of nodes. 

In the text 
Fig. 2
Left: bestfit orbital solution compared with radial velocities of the primary (blue) and of the secondary (red). Right: bestfit orbital solution compared with the pseudoastrometric uncertainty regions from KeckI (colors) and VLTI (black, smaller than the symbol size). In both plots the black circles are the predictions from the bestfit orbital solution. 

In the text 
Fig. 3
Probability distributions for the most relevant parameters (inclination, distance, and masses of the components). 

In the text 
Fig. 4
Bestfit synthetic fluxes (horizontal bars) to the spectralenergy distribution from NOMAD, 2MASS, and WISE (red circles). The corresponding synthetic spectra for the two individual components and their sum are overlaid (dashed, dotted, and solid gray lines). 

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.