Issue 
A&A
Volume 579, July 2015



Article Number  A59  
Number of page(s)  16  
Section  Stellar structure and evolution  
DOI  https://doi.org/10.1051/00046361/201425568  
Published online  29 June 2015 
Gridbased estimates of stellar ages in binary systems
SCEPtER: Stellar CharactEristics Pisa Estimation gRid^{⋆}
^{1} INAF–Osservatorio Astronomico di Collurania, via Maggini, 64100 Teramo, Italy
email: valle@df.unipi.it
^{2} INFN, Sezione di Pisa, Largo Pontecorvo 3, 56127 Pisa, Italy
^{3} Dipartimento di Fisica “Enrico Fermi”, Università di Pisa, Largo Pontecorvo 3, 56127 Pisa, Italy
Received: 22 December 2014
Accepted: 28 April 2015
Aims. We investigate the performance of gridbased techniques in estimating the age of stars in detached eclipsing binary systems. We evaluate the precision of the estimates due to the uncertainty in the observational constraints – masses, radii, effective temperatures, and [Fe/H] – and the systematic bias caused by the uncertainty in convective core overshooting, element diffusion, mixinglength value, and initial helium content.
Methods. We adopted the SCEPtER grid, which includes stars with mass in the range [0.8; 1.6] M_{⊙} and evolutionary stages from the zeroage main sequence to the central hydrogen depletion. Age estimates have been obtained by a generalisation of the maximum likelihood technique described in our previous work.
Results. We showed that the typical 1σ random error in age estimates – due only to the uncertainty affecting the observational constraints – is about ± 7%, which is nearly independent of the masses of the two stars. However, such an error strongly depends on the evolutionary phase and becomes larger and asymmetric for stars near the zeroage main sequence where it ranges from about + 90% to −25%. The systematic bias due to the including convective core overshooting – for mild and strong overshooting scenarios – is about 50% and 120% of the error due to observational uncertainties. A variation of ± 1 in the heliumtometal enrichment ratio ΔY/ ΔZ accounts for about ± 150% of the random error. The neglect of microscopic diffusion accounts for a bias of about 60% of the error due to observational uncertainties. We also introduced a statistical test of the expected difference in the recovered age of two coeval stars in a binary system. We find that random fluctuations within the current observational uncertainties can lead genuine coeval binary components to appear to be noncoeval with a difference in age as high as 60%.
Key words: binaries: eclipsing / methods: statistical / stars: evolution / stars: lowmass
Appendix A is available in electronic form at http://www.aanda.org
© ESO, 2015
1. Introduction
The possibility of obtaining direct measurements of stellar mass and radius makes eclipsing binary systems the objects of fundamental importance (see e.g. Andersen 1991; Torres et al. 2010). Nowadays, for stars in doublelined systems, these quantities can be determined with precisions of a few percentage points or better (see e.g. Clausen et al. 2008; Pavlovski et al. 2014). Thus, detached eclipsing binaries offer a unique test bed for evolutionary models of single stars, giving the opportunity of obtaining accurate estimates of their ages.
However, this age estimation is usually performed by fitting the observed data with isochrones or stellar evolutionary tracks (see among many others Southworth et al. 2011; Torres et al. 2012; Vos et al. 2012; Pavlovski et al. 2014). Recently, different techniques have been established for single stars, based upon a maximum likelihood estimation over a grid of precomputed models (e.g. Gai et al. 2011; Basu et al. 2012; Gennaro et al. 2012; Mathur et al. 2012; Prada Moroni et al. 2012; Valle et al. 2014, 2015). Such methods are very fast and flexible since they can simultaneously take all the available observables into account, reaching high precision in the estimates and selfconsistently evaluating their errors.
Despite these appealing features, very few efforts can be found in the literature that address the performance of these techniques for binary systems and that investigate the possible sources of their internal bias. Two studies that address, among other topics, the problem of binary ages are by Gennaro et al. (2012) and Schneider et al. (2014). The former is a Bayesian investigation of the age reconstruction of premain sequence stars, which also considers binary systems. The authors adopted nine synthetic test systems, subjected each of them to 100 Monte Carlo perturbations, and evaluated the fraction of fake noncoeval estimates at 1σ level. The results indicated a fraction of erroneous noncoevality from 0% to 95% with a mean of 13%. Given the restricted sample of explored systems, Gennaro et al. (2012) recognised the importance of a full analysis of the expected differences in age estimates in different ranges of mass, metallicity, and evolutionary phases. The last was a study by Schneider et al. (2014), who tested a Bayesian technique for age estimation and, st other tasks, compared the inferred stellar ages for the primary and secondary stars of eclipsing binaries in the mass range between 4.5 and 28 M_{⊙}. In a sample of 18 binary systems, they found agreement at 1σ level in the estimated ages of 17 of them. They also evaluated the mean difference in age of 14 systems (systems with mass ratio higher than 0.97 were excluded) to be 0.9 ± 2.3 Myr at the 95% confidence level.
The aim of this paper is to theoretically analyse the accuracy of gridbased age estimates for detached eclipsing binary systems. We also study the bias on these estimates due to the uncertainties in some physical mechanisms adopted in the models (convective core overshooting and element microscopic diffusion efficiencies, mixinglength value), and in the stellar chemical composition (mainly the initial helium content).
Finally, we introduce a statistical test of the expected difference in the reconstructed ages of two coeval binary components simply due to the observational uncertainties. Such a test is of interest because the recovered noncoevality of stars in binary systems – i.e. the inability to fit both components with a single isochrone – is often used to claim some deficiency in current generation of stellar models. This sometimes leads to introducing and/or calibrating some physical processes in evolutionary codes, such as convective core overshooting, or to varying the external convection efficiency through the mixinglength parameter (see e.g. Andersen 1991; Pols et al. 1997; Ribas et al. 2000; Torres et al. 2006, 2010, 2014; Morales et al. 2009; Clausen et al. 2009, 2010, and references therein).
We restrict our analysis to central hydrogenburning stars with mass in the range from 0.8 to 1.6 M_{⊙}. A similar investigation has already been performed (Valle et al. 2015) for the age determination of isolated single stars for which the main observables adopted in the reconstruction were different from the present ones. In that case neither the mass nor the radius were available, while in the present study we do not rely on the asteroseismic observables. As discussed in the following, the performances and the systematic biases of gridbased recovery techniques sensitively depend on the specific set of observables actually adopted in the recovery procedure.
The structure of the paper is the following. In Sect. 2 we discuss the method and the grids used in the estimation process. The main results are presented in Sects. 3 and 4. In Sect. 5 we present a statistical test of the expected differences between the singlestar age estimates. Some concluding remarks can be found in Sect. 6. Appendix A describes the sampling strategy followed to build the synthetic dataset.
2. Gridbased recovery technique
The basic technique, derived from Basu et al. (2012) and suitable for isolated star estimation, is described in Valle et al. (2014, 2015), hereafter V14 and V15. Age estimates are based upon a modified SCEPtER (Stellar CharactEristics Pisa Estimation gRid) scheme, adapted for binary stars. The code and the standard grid developed for this work are available in the R package SCEPtERbinary^{1} on CRAN.
The adopted implementation assumes that and are detached binary system stars for which the following vectors of observed quantities are available: , , . We let be the nominal uncertainty in the observed quantities. For each point j on the estimation grid of stellar models, we define q^{j} ≡ { T_{eff,j}, [Fe/H] _{j},M_{j},R_{j} }. Let ℒ^{1,2}_{j} be the singlestar likelihood functions defined as (1)where (2)The singlestar likelihood functions are independently evaluated for the two stars for each grid point within 3σ of all the variables from ; let ℒ^{1,2}_{max} be the two maximum values obtained in this step. The singlestar ages are estimated by averaging the corresponding quantity of all the models with likelihood greater than . Informative priors can be inserted as a multiplicative factor in Eq. (1), as a weight attached to the grid points.
In the following we assume, implicitly or explicitly, that the stars in the binary system are coeval. This assumption justifies the adoption of a single age as representative of the whole system. A first possible way to estimate the age of the binary system is simply to take the mean of the ages of the single components, computed independently. This method has a computational complexity of O(n_{1} + n_{2}), where n_{1,2} are the sizes of the 3σ gridpoint samples.
Fig. 1 Left: age relative error in dependence on the mass of the primary star with its 1σ envelope (red solid line). The blue long dashed line marks the relative errors median. Right: same as the left panel, but in dependence on the relative age of the primary star. 

Open with DEXTER 
A second approach explicitly assumes in the likelihood computation that the stars are coeval and computes the joint likelihood function only for the couples of models in the two 3σ boxes with ages within 10 Myr. The joint likelihood is computed as the product of the singlestar likelihood functions. Let be the maximum value obtained in this step. The jointstar estimated age is obtained by averaging the corresponding quantity of all the couples of models with likelihood greater than . This technique has a computational complexity of O(n_{1} × n_{2}). Readers interested in the technical solutions implemented are referred to the code provided in the abovementioned SCEPtERbinary R package. This estimation method is assumed as our standard in the following.
The described technique is similar to those previously employed in the literature in the framework of isochrone χ^{2} fitting for binary age estimations (see, amongst many, Pols et al. 1997; Lastennet & VallsGabaud 2002). The main differences are that the quoted χ^{2} approaches adopt a timeconsuming numerical functional minimisation as the best model estimates and use χ^{2} profiles to infer the errors on the estimated parameters, often fixing the value of the masses to their observed values. However these methods can lead to possibly severe error underestimation (see e.g. Basu et al. 2010; Quirion et al. 2010). The technique adopted in this and similar gridbased approaches (e.g. Gai et al. 2011; Basu et al. 2012) computes as best estimate a local mean of the bestmatching models, as described above. The error on the estimates are then computed by a Monte Carlo simulation. This approach allows a confidence interval to be provided for the age estimate and observational parameter correlations to be considered that are not addressed in χ^{2} profiles (see e.g. Lastennet & VallsGabaud 2002; Basu et al. 2010).
Although in the current paper we are not interested in obtaining a statistical confidence interval for the age of observed systems, the adopted method can be easily used to this purpose. To address this point, we reside on a generation of a synthetic sample of n binary systems, starting from the observed values and following, for each star, a multivariate normal distribution with vector of mean and covariance matrix Σ. For the Monte Carlo simulations, a value of n = 10 000 can be adopted since it provides a fair balance between computation time and accuracy of the results^{2}. The median of the age of the n systems, obtained with one of the two methods described above, is taken as the best estimate of the true values; the 16th and 84th quantiles of the n values are adopted as a 1σ confidence interval.
2.1. Standard stellar model grid
The standard estimation grid of stellar models is obtained using FRANEC stellar evolution code (Degl’Innocenti et al. 2008; Tognelli et al. 2011), in the same configuration as adopted to compute the Pisa Stellar Evolution Data Base^{3} for lowmass stars (Dell’Omodarme et al. 2012; Dell’Omodarme & Valle 2013).
The grid consists of 141 680 points (110 points for 1288 evolutionary tracks), corresponding to evolutionary stages from the zeroage main sequence (ZAMS) to the central hydrogen depletion^{4}. Models are computed for masses in the range [0.80; 1.60] M_{⊙} with a step of 0.01 M_{⊙}. The initial metallicity [Fe/H] is assumed in the range [−0.55; 0.55] with a step of 0.05 dex. The solar scaled heavyelement mixture by Asplund et al. (2009) is adopted. The initial helium abundance is obtained using the linear relation with the primordial abundance Y_{p} = 0.2485 from WMAP (Cyburt et al. 2004; Steigman 2006; Peimbert et al. 2007a,b), and assuming ΔY/ ΔZ = 2 (Pagel & Portinari 1998; Jimenez et al. 2003; Gennaro et al. 2010). The models are computed by assuming the solarscaled mixinglength parameter α_{ml} = 1.74. The convective core overshooting extension is not considered. Further details on the input adopted in the computations are available in Valle et al. (2009, 2013a,b).
SCEPtER median (p_{50}) and 1σ envelope boundaries (p_{16} and p_{84}) for age relative error as a function of the masses, relative ages, and mass ratio of the stars.
3. Binary system age estimates: internal accuracy
To evaluate the internal accuracy of gridbased age estimates of binary systems, we analysed the ideal case in which stellar models are in perfect agreement with stars. This is the most favourable case where the error in the estimated ages arises only from the uncertainties affecting the observational constraints used in the recovery procedure.
To achieve such an ideal situation, the first ageestimation test was performed on a synthetic dataset obtained by sampling N = 50 000 artificial binary systems from the same standard estimation grid of stellar models used in the recovery procedure. To simulate the effect of observational uncertainties, we added a Gaussian noise in all the observed quantities. We assumed as standard deviations 100 K in T_{eff}, 0.1 dex in [Fe/H], 1% in mass, and 0.5% in radius. Since the observationally inferred values of a given physical quantity for the two binary components are correlated, correlation coefficients should be used in the covariance matrix Σ whenever a realistic noise has to be simulated. In our standard scenario, we assumed a correlation of 0.95 between the primary and secondary effective temperatures, 0.95 between the metallicities, 0.8 between the masses, and no correlation between the radii. Motivations of these choices and an analysis of the impact on the results of different assumptions of correlations amongst stellar quantities, radii included, are presented in Sect. 3.1.
Figure 1 shows the Monte Carlo relative error in estimated age as a function of the mass and the relative age of the primary star. In this and in the following analogous figures, a positive value of the age relative error corresponds to an overestimated age. The relative age is defined as the ratio between the age of the star and the age of the same star at central hydrogen exhaustion (the age is conventionally set to 0 at the ZAMS position). It shows the position of the median relative error and the position of the 1σ envelope of agerelative error as a function of mass and relative age. The envelope is obtained as in V15, computing the 16th and 84th quantiles of the relative errors over a moving window in mass and relative age^{5}.
Table 1 reports the position of the envelope as a function of the mass of the two stars, of their relative age, and of the mass ratio q of the system^{6}. Table 1 and Fig. 1 show that the typical relative age error as a function of the primary star mass is about 7%. The envelope is larger for models near the ZAMS, where the envelope is highly asymmetric and – for the primary star – it ranges from −23% to 93%, with a clear edgeeffect distortion (see the extensive discussion in V14 and V15). For relative ages of the primary star larger than 0.4 the precision is always better than 12%. This accuracy is about four times better than is attainable by adopting asteroseismic constraints (i.e. the average large frequency spacing Δν and the frequency of maximum oscillation power ν_{max}), without the knowledge of mass and radius of the stars (V15 and references therein). Moreover, in V15 we showed that the systematic in age estimates due to different stellar codes was of the same order of magnitude as the random errors. Since the random component is significantly smaller for binary systems, we expect the systematic bias due to the adoption of different grids of stellar models to play a major role in the age determinations.
Several resulting features may depend in principle on the technique employed to obtain the system age and on the adopted perturbation strategy, so we devoted Sects. 3.1 and 3.2 to discussing the influence of these possible bias sources. A detailed discussion of the sampling strategy adopted to build the synthetic dataset and an analysis of its effect on the obtained results can be found in Appendix A.
3.1. Effect of correlations
To simulate the observational uncertainties (see Sect. 2), we added a Gaussian perturbation to the four observables (i.e. M, R, T_{eff}, [Fe/H]) of the stars in the previously generated synthetic dataset of artificial binary systems. In the standard scenario they are sampled from a multivariate normal distribution assuming correlation of ρ = 0.95 between the two effective temperatures, ρ = 0.95 between the two metallicities, ρ = 0.8 between the two masses, and zero between the two radii. These correlations arise from the fact that the analysis of data for binary stars usually determines some function of the observables at higher accuracy than the observables themselves.
The strongest correlation is expected between the effective temperatures of the two stars, since the difference in or the ratio of the effective temperatures in an eclipsing system is usually determined with relative accuracy that is about three times greater than that of the individual effective temperatures themselves (see e.g. Claret 2003; Southworth & Clausen 2007; Southworth 2013; Torres et al. 2014). A direct computation, by means of error propagation under normality assumption, on the data presented in the cited literature showed that the correlation of the two effective temperatures is reproduced well when assuming ρ = 0.95.
Regarding [Fe/H], we adopted a high correlation coefficient (ρ = 0.95) since usually only a mean value for the entire system is presented in the analyses. For mass determination, the mass ratio q is usually determined with relative accuracy of about 1.6 higher than the mean one of the masses themselves (see e.g. Popper et al. 1986; Lacy et al. 2008; Torres et al. 2009; Brogaard et al. 2011; Vos et al. 2012; Sandquist et al. 2013). Assuming normality, this corresponds to a correlation of ρ = 0.8 between the two star masses. However, even significantly lower values for the correlation (ρ ≈ 0.15) resulted from literature examination (Hełminiak et al. 2009; Meibom et al. 2009).
The situation is much more complex in the case of radii determination. The light curve analysis poses independent constraints on the sum of radii r_{s} and on their ratio k, causing a dependence between the two evaluated stellar radii; however, the magnitude and the sign of the correlation between the estimated radii depend on the actual values of r_{s} and k and of their errors. Direct computations from the values quoted in the literature showed correlations varying from ρ = 0.8 to ρ = −0.5 (Popper et al. 1986; Grundahl et al. 2008; Torres et al. 2009; Brogaard et al. 2011; Vos et al. 2012). Luckily, such high variability is not a serious problem because, as shown below, the actual magnitude of the correlation between radii has very little impact on the results. Thus, for the reference scenario we assumed uncorrelated radii.
Fig. 2 Upper row: 1σ envelope of the age relative error in dependence on the mass of the primary star (left panel), and on its relative age (right panel). The solid black line corresponds to accounting for correlations in covariance matrix (see text) and assumes that the stars are coeval; the red solid line assumes the same covariance matrix without assumption of coeval stars in the recovery (see Sect. 2); the dashed line assumes coeval stars and a diagonal covariance matrix. Lower row: same as the upper row, but for the secondary star. 

Open with DEXTER 
To quantify the impact on the age estimate uncertainty of accounting for the aforementioned correlations in the procedure followed to build the synthetic datasets, we applied the SCEPtER pipeline on datasets produced by relying on different assumptions about correlations. As a first test, we provided a dataset of artificial binary stars with uncorrelated observables; i.e., the Gaussian noise was added assuming a diagonal covariance matrix. Although such a choice is quite extreme and not realistic, the comparison between the outcomes of the age estimate procedure on this dataset with those on the standard one will quantify the maximum effect of the correlations. Figure 2 shows the impact of the correlation between observables on the relative errors in the estimated age of the system. The neglect of correlation on the corresponding observables of the two stars causes a mean relative shrinkage of the envelope of about 20%. Direct simulations (not shown) proved that the shrinkage is due to dropping the correlation between the two masses, while the correlations between the effective temperatures and between the metallicities account for modest variations.
As a second test, we built two synthetic datasets by assuming the same correlations of the standard case, with the exception of the correlation coefficient between the radii of the binary components, namely ρ = ± 0.7. The comparison of the performances of the recovery procedure on these datasets with those on the standard one will assess the impact of the correlation amongst the radii. Table 2 reports, as a function of the mass of the primary star, the results obtained accounting for correlations between the radii. The mean variations on the envelope boundaries are about 0.2%, so it is safe to neglect the correlation among radii.
3.2. Recovery algorithms
Figure 2 also shows the impact of the adopted agerecovery technique. The simple average of the individually independent age estimates of the two stars appears to be as good as the joint estimate in the case of a primary star of low mass (M ≤ 1.1M_{⊙}), while it shows a notable larger variance in binary systems composed of a highmass primary and a lowmass secondary. For these combinations of objects, the secondary stars is sampled in the first stages of the evolution, since the maximum time is set by the faster evolution of the more massive star. The age estimate of a star in the early evolutionary stages is subject to larger relative uncertainty, and this uncertainty is propagated in the final mean. Such an effect is not present in the joint evaluation, since in this case the age is mainly determined by the more evolved star. In fact the age spread of the models in the 3σ box around the primary star is generally lower than the spread the secondary star, since the primary is sampled at a higher relative age than the secondary. Therefore the joint likelihood estimate is limited by the range of age spanned by the primary star.
Median (p_{50}) and 1σ envelope boundaries (p_{16} and p_{84}) for age relative error as a function of the mass of the primary star, assuming correlations ρ = ±0.7 between the radii.
The increase in the variance in age estimates when coevality is not explicitly assumed is shown further in Fig. 3. The top row of the figure shows, as a function of the masses M_{1} and M_{2} of the binary system components, the position of the 2D lower (left panel) and upper (right panel) 1σ envelopes^{7}. To summarise the results, the figure reports the mean value (expressed as percent) of the envelope boundaries in three mass ranges: M_{1} < 1.1 M_{⊙} and M_{2} < 1.1 M_{⊙}; M_{1} > 1.1 M_{⊙} and M_{2} > 1.1 M_{⊙}; M_{1} > 1.1 M_{⊙} and M_{2} < 1.1 M_{⊙}. The lower row of the figure displays the ratio of the boundaries for joint likelihood and independent age estimates. It is apparent that the two ageestimation techniques give nearly equivalent results as long as the system stars are close to each other in mass (q ≈ 1), while for unbalanced systems, the ratio of the boundary is as low as 0.2 for M_{1} = 1.6 M_{⊙}, M_{2} = 0.8 M_{⊙}, which is, an increase in the envelope width of a factor of five.
Fig. 3 Upper row: left lower boundary of the 1σ 2D relative error envelope for joint likelihood estimates as a function of the mass of the binary system stars. The percentages on the plot refer to mean values in different ranges of mass (see text). Right: same as the left panel, but for the upper boundary of the 2D relative error envelope. Lower row: left ratio of lower boundary of the 1σ 2D relative error envelope due to different recovery techniques. The ratio is computed by dividing the position of the joint likelihood estimates boundary by the corresponding value obtained averaging the independently estimated ages. Right: same as the left panel, but for the upper boundary of the 2D relative error envelope. 

Open with DEXTER 
3.3. Contribution of the secondary star to age estimates
In Sect. 3.2 we discussed the dominant impact of the primary star on joint likelihood age determination. It is therefore interesting to explore the contribution of the secondary star to the whole system age determination. We performed this test by comparing the envelopes for the joint likelihood age estimate and for the primary singlestar estimates. The results are presented in Fig. 4. Computed as a function of the primary mass, the joint estimate’s relative error envelope shows negligible differences for both the upper and the lower boundaries with respect to the envelope computed only for primary stars. The contribution of the secondary star is only visible for young systems, where the age of the primary star is less constrained. In fact, the righthand panel of the figure shows that the relative age envelope that takes the secondary contribution into account is less prone to age underestimation at a relative age lower than 0.15. In summary, the secondary star plays a minor role with respect to the primary one because its age is generally estimated with a larger uncertainty being closer to the ZAMS than the more massive companion of the same age.
Fig. 4 Left: 1σ envelope of the agerelative error in dependence on the mass of the primary star (solid line), compared with the same quantity computed only with the data of the primary star (dashed line). Right: same as the left panel, but as a function of relative age of the primary star. 

Open with DEXTER 
Fig. 5 Left: standard 1σ envelope of the binary system age relative error in dependence on the mass of the primary star (solid black), compared with the same quantity obtained by sampling from grids with different core overshooting efficiency β = 0.2 (dashed red) and β = 0.4 (dot dashed green). Middle: same as the left panel, but in dependence on the mass of the secondary star. Right: same as the left panel, but in dependence on the system mass ratio q. 

Open with DEXTER 
4. Stellar model uncertainty propagation
The accuracy of age estimates of real binary systems obviously depends on the reliability of the adopted grid of stellar models. In this section, following V14 and V15, we provide an analysis of the impact of the various uncertainty sources that affect the stellar computations. We focus here on some issues considered in V15: the convective core overshooting efficiency, the initial helium content, and the neglect of the microscopic diffusion. We constructed several artificial grids characterised by the variation in the single input under examination up to its extreme allowed values, while keeping all the others fixed to their reference values. We then built the corresponding synthetic datasets by sampling artificial binary systems from these nonstandard grids of stellar models and by adding a Gaussian noise in all the observed quantities to simulate observational uncertainties following the same prescriptions for covariance matrix described in the standard case. Finally, the ages are estimated by means of the SCEPtER pipeline by relying on the standard grid of stellar models.
Although the use of a solarcalibrated mixinglength value for stars that differ from the Sun for mass, composition, or evolutionary phase could be inappropriate in the present paper; unlike V15, we do not consider this source of uncertainty. Such a choice is motivated by the fact that the dependence of the mixinglength value on the mass, metallicity, level of activity, and the evolutionary phase of the star is still under study and that no definitive conclusion has been reached (see, among many, Clausen et al. 2009; Deheuvels & Michel 2011; Bonaca et al. 2012; Mathur et al. 2012; Tanner et al. 2014). In V15 we studied the impact on age estimates due to a change in the mixinglength value by adopting synthetic grids with nonsolar mixinglength values, which did not vary inside a given grid regardless of the mass and metallicity of the star. Since in the present work, we consider systems of stars with different masses, hence in different evolutionary stages, we feel that the approach of V15 might not be appropriate. Besides the effect of a change in the mixing length of the sampling grid, a differential effect of the possible evolution of the mixinglength value with the mass of the stars could be important. Nevertheless, we performed a numerical experiment, by sampling from grids with α_{ml} ± 0.24 with respect to our solarcalibrated value (α_{ml} = 1.74), which showed a bias ranging from about ∓ 4.0% to ∓ 6.0%, the lowest values for mass ratio q near 1.0. These simulations suggest that the bias due to the different choices of the mixinglength values is certainly important but not dominant. For a more rigorous investigation, better knowledge of the dependence of the external convection efficiency on stellar parameters, such as mass, chemical composition, evolutionary phase, and stellar activity, is mandatory.
The observables adopted here in the age reconstruction are different from those in V15, which made use of classical (i.e. the effective temperature of stars and their metallicity [Fe/H]) and asteroseismic (i.e. the average large frequency spacing Δν and the frequency of maximum oscillation power ν_{max}) observable constraints. We therefore expect that the impact of the uncertainty sources is different from what is found in our previous work. In particular, the availability of precise stellar masses is paramount since it severely restricts the number of stellar models that can actually contribute to the recovery procedure.
SCEPtER bias on the recovered age due to the change in the stellar code input.
SCEPtER bias on the recovered age due to the change in the stellar code input.
4.1. Convective core overshooting
As in V15 we parametrised the extension of the extramixing region beyond the canonical border, as defined by the Schwarzschild criterion, in terms of the pressure scale height H_{p}: l_{ov} = βH_{p}, where β is a free parameter. To quantify the impact of taking the convective core overshooting into account, we computed – only for models more massive than 1.1 M_{⊙}− two additional grids with values of β = 0.2 and 0.4, the latter representing a possible maximum value (see e.g. the discussion in Valle et al. 2009). We then extended these two grids of stellar models at masses lower than 1.1 M_{⊙} with the standard one, which neglect core overshooting.
Fig. 6 Left: relative age bias due to the mild core overshooting as a function of the masses of the binary system. Artificial stars are sampled from the grid with mild core overshooting (β = 0.2) and their age is estimated on the standard grid. Right: same as the left panel, but for the strong overshooting scenario (β = 0.4). 

Open with DEXTER 
From these nonstandard grids we built two synthetic data sets (for β = 0.2 and β = 0.4) by sampling N = 50 000 artificial binary systems each and adding a Gaussian noise. The age estimate is then performed by adopting our standard grid of models (without overshooting) in the recovery procedure. The results are shown in Tables 3 and 4 and in Figs. 5 and 6. To evaluate the magnitude of the obtained biases, they can be compared with the standard 1σ envelope due to observational uncertainties in Table 1. It turns out that the median bias due to the mild (i.e. β = 0.2) core overshooting scenario – computed as function of the mass ratio q – is about 50% of the half width of the standard envelope uncertainties. The corresponding bias for the strong (β = 0.4) core overshooting scenario is about 120% of the standard envelope half width. In both cases the relative importance of the bias is generally greater for low values of q; for the strong overshooting scenario, in particular, the bias reaches values of about 160% of the standard envelope half width at q = 0.5, while it is only 60% at q = 1.0. This trend is expected since part of the systems with q near 1.0 are composed of stars with both masses lower than 1.1 M_{⊙}, for which core overshooting is not taken into account so that the discussed effect is partially masked.
Figure 6 shows that the dominant effect for the bias is due to the mass of the primary star. As an example, the mean bias for mild overshooting scenario (left panel in the figure) is −0.2% for a primary star that is less massive than 1.1 M_{⊙}, while it is −3.7% for a massive primary coupled to a light secondary star. The shift due to having a massive secondary star is −6.0%, so it only accounts for an additional −2.3% change. The contribution of the secondary star is even less important in the strong overshooting scenario (right panel in the figure), when it accounts for about onefifth of the primary star one.
The results presented in this section have particular importance since binary systems are often adopted to calibrate the efficiency of the convective core overshooting by isochrone fine tuning (see e.g. Claret 2007; Lacy et al. 2008; Clausen et al. 2010). Although our results do not directly allow estimation of the errors in the recovered best convective coreovershooting parameter that is consistent with the observations, we found that the bias induced by a mild convective core overshooting scenario on age estimate is about onehalf of the 1σ expected error in age estimate. This implies a difficulty unambiguously identifying the overshooting effect, which can be masked by random fluctuations.
The present results suggest that a mild overshooting scenario is hardly distinguishable from a no overshooting one. This result agrees with the finding by Claret (2007) in the lower mass range considered in that paper. Although a statistical calibration – obtained from several systems – can be meaningful, the large errors that propagate into the final calibrated overshooting parameter shed some doubt on the calibration of single systems. Owing to the widespread use of binary systems for calibration purposes, we feel that further research is needed to statistically quantify the errors in the calibrated parameters.
4.2. Initial helium abundance
We quantified the impact on the age estimate of the current uncertainty in the initial helium abundance following V15. We computed two additional grids of stellar models with the same metallicity values Z as in the standard grid, but by changing the heliumtometal enrichment ratio ΔY/ ΔZ to values of 1 and 3. Then, we built two synthetic datasets, each of N = 50 000 artificial binary systems, by sampling the objects from these two nonstandard grids and adding Gaussian noise. The age of the objects is reconstructed using the SCEPtER pipeline relying on the standard grid with ΔY/ ΔZ = 2.
Fig. 7 Left: standard 1σ envelope of the binary system age relative error in dependence on the mass of the primary star (solid black), compared with the same quantity obtained by sampling from the grid with a different initial helium abundance computed assuming ΔY/ ΔZ = 1 (dashed red) and ΔY/ ΔZ = 3 (dot dashed green). Middle: same as the left panel, but in dependence on the mass of the secondary star. Right: same as the left panel, but in dependence on the system mass ratio q. 

Open with DEXTER 
Fig. 8 Left: relativeage bias due to adopting a low initial helium abundance as a function of the masses of the binary system. Artificial stars are sampled from the grid with ΔY/ ΔZ = 1, and their age is estimated on the standard grid (ΔY/ ΔZ = 2). Right: same as the left panel, but sampling from a grid with ΔY/ ΔZ = 3. 

Open with DEXTER 
Fig. 9 Left: standard 1σ envelope of the binarysystem agerelative error in dependence on the mass of the primary star (solid black), compared with the same quantity obtained by sampling from the standard grid and reconstruction on a grid of stellar models computed without microscopic diffusion (dashed red). Middle: same as the left panel, but in dependence on the mass of the secondary star. Right: same as the left panel, but in dependence on the system’s mass ratio q. 

Open with DEXTER 
The results are presented in Tables 3 and 4 and in Figs. 7 and 8. The mean impact of the initial helium abundance uncertainty is higher than for the convective core overshooting. The median bias for ΔY/ ΔZ = 1 – computed as a function of the mass ratio q (Table 4) – is about the −140% of the halfwidth of the standard 1σ envelope due to observational uncertainties reported in Table 1. The corresponding bias for ΔY/ ΔZ = 3 is about 160% of the standard envelope half width in Table 1. Figure 8 shows that both the changes in the bias due to the primary and secondary stars are modest, so, in the examined mass range, the initial helium induced bias can be considered almost constant at about ± 10% for a change in ΔY/ ΔZ by ± 1.
The results presented in this section on the initial helium uncertainty impact are very different from those reported in V15, where this uncertainty source have been found to be nearly negligible. This striking difference results from the observational constraints adopted in the estimation. The analysis in V15 in fact showed that the little influence of the helium content on age estimates came from two opposite effects that nearly cancel each other. The first obvious effect was due to the impact of the helium content on the evolutionary time scale, so that heliumrich models evolve faster. A second effect was due to the placement on the standard recovery grid of the nonstandard helium models; heliumrich models mimicked more massive standard models, leading to a mass overestimation bias and therefore to reduced age estimates. The two effects compensate for each other.
This balancing effect cannot occur here since the mass is adopted as an observational constraint, and this leads to the strong influence of the initial helium content on the final age estimates.
4.3. Element microscopic diffusion
We evaluated the impact on the age estimates of neglecting microscopic diffusion in the stellar models adopted in the recovery procedure following, as in V15, a different approach with respect to the previous cases. Here we preferred to build the synthetic data set by sampling N = 50 000 artificial binary systems from stellar models that take diffusion into account and by adding Gaussian noise. We then reconstructed their age by means of the SCEPtER pipeline but adopting the nonstandard grid of models, which neglect diffusion.
Fig. 10 Relativeage bias due to neglecting elements diffusion as a function of the masses of the binary system. Artificial stars are sampled from the standard grid, and their ages are estimated on a grid computed without microscopic diffusion. 

Open with DEXTER 
The results are presented in Tables 3 and 4 and in Figs. 9 and 10. The bias due to the microscopic diffusion has a mean value of 3.6% (see Table 4), it is similar in absolute value to the one due to mild overshooting, and it ranges from about onehalf to onethird of the bias due to initial helium change. Figure 10 confirms an expected result that the diffusion is most important in systems composed of two lowmass stars. For primary and secondary masses lower than 1.1 M_{⊙}, the bias in age estimates is 5.6%, while it drops to 3.7% for more massive primary stars. The trend can be understood since the evolutionary time scale of these stars is faster than for the diffusion processes, which are therefore quite inefficient (see left panel in Fig. 9).
As already discussed in Sect. 4.2, this result confirms that the bias in age grid estimates, owing to the different input in the evolutionary stellar codes, and strongly depends on the assumed observational constraints. In fact V15 has shown that – if adopting asteroseismic constraints – the microscopic diffusioninduced bias dominated the other examined biases, while it is not the case here.
5. The coevality problem: a statistical approach
In the literature there are several cases of real eclipsing binary systems whose components are estimated to be noncoeval, because they cannot be fitted simultaneously by a single isochrone (e.g. Clausen et al. 2009; Vos et al. 2012; Sandquist et al. 2013). Such an occurrence is often used to claim that the current generation of stellar models present some weaknesses and that additional physical processes should be included or recalibrated. Given its considerable implications, it is worth discussing whether the estimates of noncoevality amongst the two components of real binary systems are statistically reliable or are simply a fluctuation. In this regard it is of particular interest to quantify the expected difference in the estimated ages of two genuine coeval stars caused simply by the observational uncertainties. To do that, we analysed our previously described synthetic datasets of coeval binary stars.
As detailed in Sect. 2, the artificial binary systems sampled from the grid of stellar models are subject to random perturbations to account for the observational errors. Therefore the ages of the single stellar members, which are equal by construction in the sampling stage, are generally estimated as different. It is interesting to estimate how large the expected differences are in age owing only to these random perturbations, since this effect also interests real world observations.
A detailed analysis of the problem has to deal with many technical aspects^{8} and it is therefore outside the aims of the present work. A devoted paper addresses the statistical aspects of this question presenting more detailed results (Valle et al., in prep.). Nevertheless, we present some basic findings here that were obtained with the standard sampling described in Sect. 2.
We define A_{1} and A_{2} as the estimated ages of the two members, with A_{1}>A_{2}. We focused our attention on the statistics W = (A_{1}−A_{2}) /A_{1}, which then ranges from 0, when the stars are correctly estimated coeval, and 1, when one star is estimated to be much younger than the other. To develop a statistical test based on W we have to estimate how large W can become only because of the random fluctuations in age estimates A_{1} and A_{2}. Once the distribution of W – by means of a Monte Carlo simulation has been empirically estimated, it is possible to choose a critical value by identifying the range of values of W that are too extreme to be compatible with the coevality hypothesis. Usually the critical value is chosen as the 95th or the 99th quantile of the distribution of the statistics under consideration. The choice of the critical value to adopt (in the following, we refer to this as W_{1−α}) defines the “level” α of the test (α thus being the probability that an observed value of W is larger than W_{1−α}). The set of values higher than the critical values are said to lie in the “rejection region”.
We present the rejection regions at level α = 0.05, hence the critical values W_{0.95} identifying them, of the null hypothesis that the stars are coeval. Reconstructed values of W higher than the critical value W_{0.95} lead to the rejection of the null hypothesis, implying that standard stellar models cannot account for the coevality of the stars with the assumed input or parameters.
A first obvious complication comes from the value of W_{0.95} depending on the mass of the two stars, on their metallicity, and on the (observationally unknown) relative age. Here we focus on the impact of the binary star masses.
To perform the criticalvalue estimation, we generated a sample of N = 50 000 artificial binary systems, perturbed their observables as explained in Sect. 2, and then independently estimated the ages of the two stars. From these estimates, we computed the statistics W. Then we binned the W values according to the primary and secondary stellar masses, as for constructing the 2D envelope in Sect. 3.1. In each of these bins, we computed the empirical W 95th quantile, which approximates the required critical value. A detailed discussion of the dependence of such a critical value on the actual uncertainties affecting the observational constraints and on the evolutionary phase of the two binary members will be presented in a forthcoming paper (Valle et al., in prep.).
Fig. 11 Critical values W_{0.95} of the statistic W for the expected differences in the ages of the two stars, only due to the observational uncertainties (see text). 

Open with DEXTER 
The results of the analysis are displayed in Fig. 11. The figure shows that genuine coeval stars can be reconstructed as noncoeval with a sizeable relative age difference only because of the current uncertainty in the observational constraints. While the critical value W_{0.95} is about 0.25 for systems of nearly equal masses, it can be greater than 0.6 for unbalanced binary masses.
This result should be carefully considered when determining the stellar ages of a binary system. Thus, a determination of two discrepant stellar ages, which leads to a value of W<W_{0.95}, does not allow a statistically grounded rejection of the hypothesis that the stars are coeval since the difference might only be due to a random fluctuation in the observational constraints.
As discussed in Sect. 4.1, one should be aware of this effect whenever the extent of convective core overshooting is calibrated on eclipsing binary by requiring that both stellar components are simultaneously fitted by a single isochrone, i.e. by imposing their coevality. If the ages estimated by means of stellar models without core overshooting give a value of W<W_{0.95}, one cannot conclude on statistical grounds that these models are not able to simultaneously fit both stars and that a convective core overshooting must be taken into account. Therefore, the reliability of the calibration method of the convective core overshooting consisting in the isochrones fine tuning to the observation has to be cautiously evaluated.
In contrast, when standard stellar models provide ages of the binary components such that W>W_{0.95}, the possibility that noncoevality is the result of observational uncertainties should be considered as negligible. As an example, binary systems in which at least one star is suspected of supporting a strong surface activity, generally show a large discrepancy among the ages of the two components, such as V636 Cent and EF Aqr (Clausen et al. 2009; Vos et al. 2012). In this system, the active secondary stars are estimated to be much older than the primary, far above the limit due to random fluctuations. The coevality is recovered adopting a significantly lower mixinglength value for the secondary stars. This is an expected behaviour since it is recognised that in active stars the dynamoinduced magnetic fields can suppress convection and produce starspots, causing differences in temperature and radii with respect to standard stellar models, which can be mimicked by a lower mixinglength value (Gabriel 1969; Cox et al. 1981; Clausen et al. 1999). It is, however, clear that further research is mandatory to determine the expected errors in the calibrated mixing length in a sound statistical way.
6. Conclusions
We performed a detailed theoretical investigation of the age estimation of detached doublelined eclipsing binaries by means of a maximum likelihood gridbased technique. We analysed the impact of the current uncertainty on several input of stellar model computations, such as the chemical composition, the efficiency of convective core overshooting, and the efficiency of microscopic processes.
We adopted the gridbased pipeline SCEPtER, which is extensively described in Valle et al. (2014, 2015). As observational constraints we assumed the stellar effective temperature, the metallicity [Fe/H], the mass, and the radius of the stars. The grid of stellar models, computed for the evolutionary phases from ZAMS to the central hydrogen depletion covers the mass range [0.8; 1.6] M_{⊙}.
We computed the statistical errors arising from the uncertainties in observational quantities. We found an overall relative uncertainty of about ± 7% in age estimates of the system, which is nearly independent of the masses of the stars and of the mass ratio of the binary systems. However, as described in Valle et al. (2015), the relative error in the age estimates sensitively depends on the evolutionary phase and becomes larger for models near those ZAMS for which it can reach values of about 90% for the upper boundary, with an envelope half width of about 50%. The large uncertainty and the bias toward higher ages found near the ZAMS are mainly an edge effect distortion, since the estimated age cannot be negative. The relative error envelope quickly shrinks as the stars evolve from the ZAMS, and for primary stars of relative age 0.3, the half width of the error envelope is about 15%.
We studied the impact on the age estimates of different choices of the sampling algorithms, of the correlation between the effective temperatures, the masses, the radii, and [Fe/H] of the two binary components, and of the method of estimating the joint age of the stars. It resulted that the first source of uncertainty plays a minor role since it accounts for a variation in the relative error half width of about 1%. The neglect of correlations between the observables of the two stars was shown to reduce the age error envelope width by about onefifth. The age estimates of the binary system obtained by averaging the ages of the two stars computed independently – without explicitly assuming coevality – lead to negligible differences only for stars of nearly equal masses. Whenever an unbalanced system is considered (M_{1}> 1.1M_{⊙} and M_{2} ≤ 1.1M_{⊙}), the envelope width of the age relative errors for the average of independent age estimates is larger by a factor of two, but it can be as high as a factor of five for systems with q ≈ 0.5.
We evaluated the systematic biases on age estimates owing to the uncertainties in initial helium content, in the convective core overshooting, and in the microscopic diffusion efficiency. We found that the bias due to a mild (β = 0.2) and strong (β = 0.4) convective core overshooting scenarios are about 50% to 120% of the half width of the standard age relative error envelope due only to the observational uncertainties. The bias due to the initial helium content was explored assuming an uncertainty of ± 1 on the heliumtometal enrichment ratio ΔY/ ΔZ. The effect is about −140% and 160% of the 1σ random error, so slightly more than that of the strong overshooting scenario. Finally, the bias due to the neglect of the element microscopic diffusion was about 60% of the half width of the standard agerelative error envelope.
The comparison of these results to those reported in Valle et al. (2015), which were computed by adopting asteroseismic constraints and without the knowledge of stellar mass and radius, showed relevant differences. In fact, Valle et al. (2015) showed that the heliuminduced bias on age estimates was negligible, while the bias due to the microscopic diffusion was dominant. These differences arise from the corresponding bias on mass estimates that could occur in Valle et al. (2015), which led to some compensations for the age estimates. These mechanisms cannot occur in the present study, since the stellar mass is assumed as an observational constraint. An important point is therefore that the bias in the gridbased age estimates – due to adopting different input in a stellar evolutionary code – strongly depends on the assumed observational constraints and should be evaluated casebycase whenever a new grid is adopted for estimates.
We also introduced a statistical test on the expected differences in the estimated ages of two genuine coeval binary components due simply to the observational uncertainties affecting the observables used in the recovery procedure. We found that the relative difference between the reconstructed ages of two coeval binary members depends on the mass ratio q. For systems with q ≈ 1, the relative difference in the inferred ages can be about 20%, while it is higher than 60% for system with q ≲ 0.8. This result should be carefully taken into account before using the apparent noncoevality of stars in binary systems – i.e. the inability to find a single isochrone that fits both components – to claim the weakness of current generation of stellar models and to calibrate the efficiency of some poorly known physical process, such as convective core overshooting, by isochrone fine tuning with the observations.
The joint density was obtained by a bidimensional kernel density estimation using the function kde2d in the R library MASS. Details on the technique can be found in Venables & Ripley (2002).
A LOESS (LOcal regrESSion) smoother is a nonparametric, locally weighted polynomial regression technique that is often used to show the underlying trend of scattered data (see e.g. Feigelson & Babu 2012; Venables & Ripley 2002).
Acknowledgments
We thank our anonymous referee for the stimulating comments that helped us to clarify and improve the paper. This work has been supported by PRINMIUR 20102011 (Chemical and dynamical evolution of the Milky Way and Local Group galaxies, P.I. F. Matteucci), and PRININAF 2012 (The M4 Core Project with Hubble Space Telescope, P.I. L. Bedin).
References
 Andersen, J. 1991, A&ARv, 3, 91 [NASA ADS] [CrossRef] [Google Scholar]
 Asplund, M., Grevesse, N., Sauval, A. J., & Scott, P. 2009, ARA&A, 47, 481 [NASA ADS] [CrossRef] [Google Scholar]
 Basu, S., Chaplin, W. J., & Elsworth, Y. 2010, ApJ, 710, 1596 [NASA ADS] [CrossRef] [Google Scholar]
 Basu, S., Verner, G. A., Chaplin, W. J., & Elsworth, Y. 2012, ApJ, 746, 76 [NASA ADS] [CrossRef] [Google Scholar]
 Bonaca, A., Tanner, J. D., Basu, S., et al. 2012, ApJ, 755, L12 [NASA ADS] [CrossRef] [Google Scholar]
 Brogaard, K., Bruntt, H., Grundahl, F., et al. 2011, A&A, 525, A2 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Claret, A. 2003, A&A, 399, 1115 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Claret, A. 2007, A&A, 475, 1019 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Clausen, J. V., Helt, B. E., & Olsen, E. H. 1999, in Stellar Structure: Theory and Test of Connective Energy Transport, eds. A. Gimenez, E. F. Guinan, & B. Montesinos, ASP Conf. Ser., 173, 321 [Google Scholar]
 Clausen, J. V., Torres, G., Bruntt, H., et al. 2008, A&A, 487, 1095 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Clausen, J. V., Bruntt, H., Claret, A., et al. 2009, A&A, 502, 253 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Clausen, J. V., Frandsen, S., Bruntt, H., et al. 2010, A&A, 516, A42 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Cox, A. N., Hodson, S. W., & Shaviv, G. 1981, ApJ, 245, L37 [NASA ADS] [CrossRef] [Google Scholar]
 Cyburt, R. H., Fields, B. D., & Olive, K. A. 2004, Phys. Rev. D, 69, 123519 [NASA ADS] [CrossRef] [Google Scholar]
 Degl’Innocenti, S., Prada Moroni, P. G., Marconi, M., & Ruoppo, A. 2008, Ap&SS, 316, 25 [NASA ADS] [CrossRef] [Google Scholar]
 Deheuvels, S., & Michel, E. 2011, A&A, 535, A91 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Dell’Omodarme, M., & Valle, G. 2013, The R Journal, 5, 108 [Google Scholar]
 Dell’Omodarme, M., Valle, G., Degl’Innocenti, S., & Prada Moroni, P. G. 2012, A&A, 540, A26 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Feigelson, E. D., & Babu, G. J. 2012, Modern Statistical Methods for Astronomy with R applications (Cambridge University Press) [Google Scholar]
 Gabriel, M. 1969, in LowLuminosity Stars, ed. S. S. Kumar, 267 [Google Scholar]
 Gai, N., Basu, S., Chaplin, W. J., & Elsworth, Y. 2011, ApJ, 730, 63 [NASA ADS] [CrossRef] [Google Scholar]
 Gennaro, M., Prada Moroni, P. G., & Degl’Innocenti, S. 2010, A&A, 518, A13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Gennaro, M., Prada Moroni, P. G., & Tognelli, E. 2012, MNRAS, 420, 986 [NASA ADS] [CrossRef] [Google Scholar]
 Grundahl, F., Clausen, J. V., Hardis, S., & Frandsen, S. 2008, A&A, 492, 171 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Hełminiak, K. G., Konacki, M., Ratajczak, M., & Muterspaugh, M. W. 2009, MNRAS, 400, 969 [NASA ADS] [CrossRef] [Google Scholar]
 Jimenez, R., Flynn, C., MacDonald, J., & Gibson, B. K. 2003, Science, 299, 1552 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Lacy, C. H. S., Torres, G., & Claret, A. 2008, AJ, 135, 1757 [NASA ADS] [CrossRef] [Google Scholar]
 Lastennet, E., & VallsGabaud, D. 2002, A&A, 396, 551 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mathur, S., Metcalfe, T. S., Woitaszek, M., et al. 2012, ApJ, 749, 152 [NASA ADS] [CrossRef] [Google Scholar]
 Meibom, S., Grundahl, F., Clausen, J. V., et al. 2009, AJ, 137, 5086 [NASA ADS] [CrossRef] [Google Scholar]
 Morales, J. C., Torres, G., Marschall, L. A., & Brehm, W. 2009, ApJ, 707, 671 [NASA ADS] [CrossRef] [Google Scholar]
 Pagel, B. E. J., & Portinari, L. 1998, MNRAS, 298, 747 [NASA ADS] [CrossRef] [Google Scholar]
 Pavlovski, K., Southworth, J., Kolbas, V., & Smalley, B. 2014, MNRAS, 438, 590 [NASA ADS] [CrossRef] [Google Scholar]
 Peimbert, M., Luridiana, V., & Peimbert, A. 2007a, ApJ, 666, 636 [NASA ADS] [CrossRef] [Google Scholar]
 Peimbert, M., Luridiana, V., Peimbert, A., & Carigi, L. 2007b, in From Stars to Galaxies: Building the Pieces to Build Up the Universe, eds. A. Vallenari, R. Tantalo, L. Portinari, & A. Moretti, ASP Conf. Ser., 374, 81 [Google Scholar]
 Pols, O. R., Tout, C. A., Schroder, K.P., Eggleton, P. P., & Manners, J. 1997, MNRAS, 289, 869 [NASA ADS] [CrossRef] [Google Scholar]
 Popper, D. M., Lacy, C. H., Frueh, M. L., & Turner, A. E. 1986, AJ, 91, 383 [NASA ADS] [CrossRef] [Google Scholar]
 Prada Moroni, P. G., Gennaro, M., Bono, G., et al. 2012, ApJ, 749, 108 [NASA ADS] [CrossRef] [Google Scholar]
 Quirion, P.O., ChristensenDalsgaard, J., & Arentoft, T. 2010, ApJ, 725, 2176 [NASA ADS] [CrossRef] [Google Scholar]
 Ribas, I., Jordi, C., & Giménez, Á. 2000, MNRAS, 318, L55 [NASA ADS] [CrossRef] [Google Scholar]
 Sandquist, E. L., Shetrone, M., Serio, A. W., & Orosz, J. 2013, AJ, 146, 40 [NASA ADS] [CrossRef] [Google Scholar]
 Schneider, F. R. N., Langer, N., de Koter, A., et al. 2014, A&A, 570, A66 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Southworth, J. 2013, A&A, 557, A119 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Southworth, J., & Clausen, J. V. 2007, A&A, 461, 1077 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Southworth, J., Pavlovski, K., Tamajo, E., et al. 2011, MNRAS, 414, 3740 [NASA ADS] [CrossRef] [Google Scholar]
 Steigman, G. 2006, Int. J. Mod. Phys. E, 15, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Tanner, J. D., Basu, S., & Demarque, P. 2014, ApJ, 785, L13 [NASA ADS] [CrossRef] [Google Scholar]
 Tognelli, E., Prada Moroni, P. G., & Degl’Innocenti, S. 2011, A&A, 533, A109 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Torres, G., Lacy, C. H., Marschall, L. A., Sheets, H. A., & Mader, J. A. 2006, ApJ, 640, 1018 [NASA ADS] [CrossRef] [Google Scholar]
 Torres, G., Sandberg Lacy, C. H., & Claret, A. 2009, AJ, 138, 1622 [NASA ADS] [CrossRef] [Google Scholar]
 Torres, G., Andersen, J., & Giménez, A. 2010, A&ARv, 18, 67 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Torres, G., Fischer, D. A., Sozzetti, A., et al. 2012, ApJ, 757, 161 [NASA ADS] [CrossRef] [Google Scholar]
 Torres, G., Vaz, L. P. R., Sandberg Lacy, C. H., & Claret, A. 2014, AJ, 147, 36 [NASA ADS] [CrossRef] [Google Scholar]
 Valle, G., Marconi, M., Degl’Innocenti, S., & Prada Moroni, P. G. 2009, A&A, 507, 1541 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Valle, G., Dell’Omodarme, M., Prada Moroni, P. G., & Degl’Innocenti, S. 2013a, A&A, 549, A50 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Valle, G., Dell’Omodarme, M., Prada Moroni, P. G., & Degl’Innocenti, S. 2013b, A&A, 554, A68 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Valle, G., Dell’Omodarme, M., Prada Moroni, P. G., & Degl’Innocenti, S. 2014, A&A, 561, A125 (V14) [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Valle, G., Dell’Omodarme, M., Prada Moroni, P. G., & Degl’Innocenti, S. 2015, A&A, 575, A12 (V15) [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Venables, W., & Ripley, B. 2002, Modern applied statistics with S, Statistics and computing (Springer) [Google Scholar]
 Vos, J., Clausen, J. V., Jørgensen, U. G., et al. 2012, A&A, 540, A64 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
Online material
Appendix A: Synthetic dataset: sampling strategy
In constructing the sample adopted in the error estimation, we did not try to match any observed distribution of the mass of the binary stars or of the mass ratio q. Since we primarily aimed to investigate the performances of gridbased age estimates of binary stars, we paid more attention when sampling the whole range in mass, metallicity, and evolutionary phases covered by the available grid of models.
In detail, the datasets of artificial stars are sampled by adopting the following scheme. First, a star is randomly sampled from the grid. No priors on mass, age, and metallicity are assumed in the sampling procedure. Then, a second star is coupled to the first one with two requirements: it must share with the first star both the same initial [Fe/H] and the age, with a tolerance of 10 Myr. The couple of stars is reordered to have the most luminous star as a first member.
Fig. A.1 Left: joint density of mass and relative ages in the sample of primary stars. The colours correspond to the different densities of probability. Right: same as the left panel for secondary stars. The colour scales of the two panels are different. 

Open with DEXTER 
Fig. A.2 Joint density of mass ratio q and primary relative age in the sample of binary stars. The solid line displays the LOESSsmoothed trend of relative age versus mass ratio q. 

Open with DEXTER 
To show which models are selected as primary and secondary stars, we estimated the joint density of mass and relative age for primary and secondary stars^{9}. From the joint density it is possible to verify if there are ranges of mass or relative age that are preferably sampled. The results are displayed in Fig. A.1. It is apparent that the density functions for the primary and secondary stars are different and that they are both not uniform. In particular, the sample of primary stars is biased towards high masses and high relative ages since these models have higher intrinsic luminosity. The distribution of secondary stars is more diffuse, since several young lowmass models are present in the sample. Figure A.2 shows the joint density of the mass ratio q and the primary relative age. The solid line in the figure is a LOESS smoother^{10} of the relative age versus q. It is apparent that for binary systems with near equal masses, the sampling returns more evolved stars. This bias comes from later evolutionary phases needing more points to adequately follow the rapid evolution. The median time step in this grid region is about 13 Myr. Once a model has been sampled in this grid region, several points in the same evolutionary track will pass the constraint on maximum age differences of 10 Myr from the first sampled model. Therefore there is a high probability that a point in the same track is coupled to the first selected stars, leading to a mass ratio of 1.0. This phenomenon is not important in the reconstruction phase, since all these points have basically the same age, as discussed in V15.
The effect of the sampling could be seen in Table 1, which shows a small shrink in the age relative error envelope at q greater than 0.9. This is a direct consequence of the relative error envelope narrowing at high relative age. To show that this artificial trend could be removed with different sampling strategies, in Table 1 we also present two additional series of results. In fact, we first explored the impact of a different sampling. For this test, we sampled the set of primary stars, then we coupled each of those objects to the secondary star by sampling only from the models respecting the constraint in metallicity and age, which are also less massive than the primary object. The results are in the rows labelled “mass ratio q (alternative sampling)” in Table 1. As a second test, to see the impact of the adhoc removal of the trend of primary relative ages versus q shown in Fig. A.2, we checked an alternative approach that rejects some objects in the subset q> 0.9, based on the standard method. In fact, since this is the region that is more populated, it was possible to draw from it a sample of size equal to that of the adjacent subset 0.8 <q ≤ 0.9, with the requirement that the new sample’s primary relative age matches the distribution of the subset 0.8 <q ≤ 0.9. The corresponding results are in the rows “mass
ratio q (rejection step)” in Table 1. It appears that the shrinkage in the upper q region is reduced in both two cases.
In contrast, the shrinkage of the error envelope at low mass ratio q (see Table 1) is due to an edge effect since there are only models in this region with mass 1.6 M_{⊙} coupled with 0.8 M_{⊙} models. In fact the lower lefthand region of Fig. A.2 is not populated because to produce a coeval combination, the more massive model cannot be too young. In fact, in this case the secondary is still in the premain sequence phase, which is not considered in the grid. Figure A.2 shows that the minimum relative age for these massive models is about 0.2. Therefore the most critical region for age estimation – that of low relative age – is avoided, and the estimate relative error envelope is narrower.
All Tables
SCEPtER median (p_{50}) and 1σ envelope boundaries (p_{16} and p_{84}) for age relative error as a function of the masses, relative ages, and mass ratio of the stars.
Median (p_{50}) and 1σ envelope boundaries (p_{16} and p_{84}) for age relative error as a function of the mass of the primary star, assuming correlations ρ = ±0.7 between the radii.
All Figures
Fig. 1 Left: age relative error in dependence on the mass of the primary star with its 1σ envelope (red solid line). The blue long dashed line marks the relative errors median. Right: same as the left panel, but in dependence on the relative age of the primary star. 

Open with DEXTER  
In the text 
Fig. 2 Upper row: 1σ envelope of the age relative error in dependence on the mass of the primary star (left panel), and on its relative age (right panel). The solid black line corresponds to accounting for correlations in covariance matrix (see text) and assumes that the stars are coeval; the red solid line assumes the same covariance matrix without assumption of coeval stars in the recovery (see Sect. 2); the dashed line assumes coeval stars and a diagonal covariance matrix. Lower row: same as the upper row, but for the secondary star. 

Open with DEXTER  
In the text 
Fig. 3 Upper row: left lower boundary of the 1σ 2D relative error envelope for joint likelihood estimates as a function of the mass of the binary system stars. The percentages on the plot refer to mean values in different ranges of mass (see text). Right: same as the left panel, but for the upper boundary of the 2D relative error envelope. Lower row: left ratio of lower boundary of the 1σ 2D relative error envelope due to different recovery techniques. The ratio is computed by dividing the position of the joint likelihood estimates boundary by the corresponding value obtained averaging the independently estimated ages. Right: same as the left panel, but for the upper boundary of the 2D relative error envelope. 

Open with DEXTER  
In the text 
Fig. 4 Left: 1σ envelope of the agerelative error in dependence on the mass of the primary star (solid line), compared with the same quantity computed only with the data of the primary star (dashed line). Right: same as the left panel, but as a function of relative age of the primary star. 

Open with DEXTER  
In the text 
Fig. 5 Left: standard 1σ envelope of the binary system age relative error in dependence on the mass of the primary star (solid black), compared with the same quantity obtained by sampling from grids with different core overshooting efficiency β = 0.2 (dashed red) and β = 0.4 (dot dashed green). Middle: same as the left panel, but in dependence on the mass of the secondary star. Right: same as the left panel, but in dependence on the system mass ratio q. 

Open with DEXTER  
In the text 
Fig. 6 Left: relative age bias due to the mild core overshooting as a function of the masses of the binary system. Artificial stars are sampled from the grid with mild core overshooting (β = 0.2) and their age is estimated on the standard grid. Right: same as the left panel, but for the strong overshooting scenario (β = 0.4). 

Open with DEXTER  
In the text 
Fig. 7 Left: standard 1σ envelope of the binary system age relative error in dependence on the mass of the primary star (solid black), compared with the same quantity obtained by sampling from the grid with a different initial helium abundance computed assuming ΔY/ ΔZ = 1 (dashed red) and ΔY/ ΔZ = 3 (dot dashed green). Middle: same as the left panel, but in dependence on the mass of the secondary star. Right: same as the left panel, but in dependence on the system mass ratio q. 

Open with DEXTER  
In the text 
Fig. 8 Left: relativeage bias due to adopting a low initial helium abundance as a function of the masses of the binary system. Artificial stars are sampled from the grid with ΔY/ ΔZ = 1, and their age is estimated on the standard grid (ΔY/ ΔZ = 2). Right: same as the left panel, but sampling from a grid with ΔY/ ΔZ = 3. 

Open with DEXTER  
In the text 
Fig. 9 Left: standard 1σ envelope of the binarysystem agerelative error in dependence on the mass of the primary star (solid black), compared with the same quantity obtained by sampling from the standard grid and reconstruction on a grid of stellar models computed without microscopic diffusion (dashed red). Middle: same as the left panel, but in dependence on the mass of the secondary star. Right: same as the left panel, but in dependence on the system’s mass ratio q. 

Open with DEXTER  
In the text 
Fig. 10 Relativeage bias due to neglecting elements diffusion as a function of the masses of the binary system. Artificial stars are sampled from the standard grid, and their ages are estimated on a grid computed without microscopic diffusion. 

Open with DEXTER  
In the text 
Fig. 11 Critical values W_{0.95} of the statistic W for the expected differences in the ages of the two stars, only due to the observational uncertainties (see text). 

Open with DEXTER  
In the text 
Fig. A.1 Left: joint density of mass and relative ages in the sample of primary stars. The colours correspond to the different densities of probability. Right: same as the left panel for secondary stars. The colour scales of the two panels are different. 

Open with DEXTER  
In the text 
Fig. A.2 Joint density of mass ratio q and primary relative age in the sample of binary stars. The solid line displays the LOESSsmoothed trend of relative age versus mass ratio q. 

Open with DEXTER  
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.