Issue |
A&A
Volume 646, February 2021
|
|
---|---|---|
Article Number | A78 | |
Number of page(s) | 22 | |
Section | Stellar structure and evolution | |
DOI | https://doi.org/10.1051/0004-6361/202039212 | |
Published online | 10 February 2021 |
Testing abundance-age relations beyond solar analogues with Kepler LEGACY stars⋆
1
Space sciences, Technologies and Astrophysics Research (STAR) Institute, Université de Liège, Quartier Agora, Allée du 6 Août 19c, Bât. B5C, 4000 Liège, Belgium
e-mail: tmorel@uliege.be
2
Université Côte d’Azur, Observatoire de la Côte d’Azur, CNRS, Laboratoire Lagrange, Bd de l’Observatoire, CS 34229, 06304 Nice Cedex 4, France
3
School of Physics and Astronomy, University of Birmingham, Edgbaston, Birmingham B15 2TT, UK
Received:
19
August
2020
Accepted:
6
November
2020
The prospects of using abundance ratios as stellar age indicators appear promising for solar analogues, but the usefulness of this technique for stars spanning a much wider parameter space remains to be established. We present abundances of 21 elements in a sample of 13 bright FG dwarfs drawn from the Kepler LEGACY sample to examine the applicability of the abundance-age relations to stars with properties strongly departing from solar. These stars have precise asteroseismic ages that can be compared to the abundance-based estimates. We analyse the well-known binary 16 Cyg AB for validation purposes and confirm the existence of a slight metal enhancement (∼0.02 dex) in the primary, which might arise from planetary formation and/or ingestion. We draw attention to systematic errors in some widely used catalogues of non-seismic parameters that may significantly bias asteroseismic inferences. In particular, we find evidence that the ASPCAP Teff scale used for the APOKASC catalogue is too cool for dwarfs and that the [Fe/H] values are underestimated by ∼0.1 dex. In addition, a new seismic analysis of the early F-type star KIC 9965715 relying on our spectroscopic constraints shows that the star is more massive and younger than previously thought. We compare seismic ages to those inferred from empirical abundance-age relations based on ages from PARSEC isochrones and abundances obtained in the framework of the HARPS-GTO programme. These calibrations depend on the stellar effective temperature, metallicity, and/or mass. We find that the seismic and abundance-based ages differ on average by 1.5–2 Gyr, while taking into account a dependency on one or two stellar parameters in the calibrations leads to a global improvement of up to ∼0.5 Gyr. However, even in that case we find that seismic ages are systematically larger by ∼0.7 Gyr. We argue that it may be ascribed to a variety of causes including the presence of small zero-point offsets between our abundances and those used to construct the calibrations or to the choice of the set of theoretical isochrones. The conclusions above are supported by the analysis of literature data for a larger number of Kepler targets. For this extended sample, we find that incorporating a Teff dependency largely corrects for the fact that the abundance-based ages are lower/larger with respect to the seismic estimates for the cooler/hotter stars. Although investigating age dating methods relying on abundance data is worth pursuing, we conclude that further work is needed to improve both their precision and accuracy for stars that are not solar analogues.
Key words: asteroseismology / stars: fundamental parameters / stars: abundances
© ESO 2021
1. Introduction
Stellar ages play a central role in various fields of astrophysics. Yet, it is one of the most difficult properties to accurately determine in field stars. Although traditionally used in the past to tag young or old stellar populations, [Fe/H] is now recognised not to be a reliable age proxy for stars in the Galactic discs (e.g., Edvardsson et al. 1993). However, several studies of solar twins and analogues spanning a narrow range of parameters have recently unveiled remarkable correlations between stellar ages (derived from isochrone fitting) and other abundance ratios (e.g., da Silva et al. 2012; Bedell et al. 2018; Nissen 2015, 2016; Nissen et al. 2017, 2020; Adibekyan et al. 2016; Spina et al. 2016, 2018; Tucci Maia et al. 2016; Jofré et al. 2020; Lin et al. 2020). Abundance ratios of elements produced through different nucleosynthesis channels (e.g., [Y/Mg]) are particularly sensitive to age because the relative amounts of ejecta released in the interstellar medium (ISM) strongly varied along the evolution of the Galaxy (e.g., Spina et al. 2016). The proposed relations between the abundance ratios and age generally extend over ∼10 Gyr and are relatively tight (scatter of ∼1 Gyr or even less for a given abundance ratio) for thin-disc stars.
Although this novel dating technique is potentially powerful, doubts remain as to whether it can be applied to different stellar populations and help in tackling some important issues, such as the formation history of our Galaxy. This question is particularly relevant and timely in view of the flood of stellar abundances that are already delivered by several medium-to-high-resolution spectroscopic surveys (e.g., Buder et al. 2020; Gilmore et al. 2012). A wealth of information into the processes that shaped our Galaxy is encoded in these data, but reconstructing the timeline of events is essential. The sensitivity of the abundance-age relations against the stellar properties remains largely unexplored. As a matter of fact, recent studies have warned that care should be exercised when applying the abundance-age relations derived from solar analogues to stars with characteristics (especially metallicity) significantly different from those of the Sun, not belonging to the Galactic thin disc, or even located outside the solar circle (e.g., Feltzing et al. 2017; Titarenko et al. 2019; Casali et al. 2020; Delgado Mena et al. 2019). The last point is a concern because old stars are believed to have radially migrated over large distances across the Galaxy (e.g., Sellwood & Binney 2002). Although these dependences cannot be ignored, they are currently neither well understood nor well quantified. Some abundance ratios appear to be relatively insensitive to some of these issues (e.g., [Y/Mg] compared to [Sr/Mg]; Nissen et al. 2020), but the dependency may be quantitatively different depending on the nucleosynthesis details of the pair of chemical elements involved.
Another line of progress is related to the fact that the calibrations are based on isochrone ages that fully rely on the predictions of evolutionary models and are of limited applicability for unevolved dwarfs. Even in favourable cases (e.g., sub-giants), the accuracy in the age determination of field stars is severely limited by systematic uncertainties (e.g., see discussion in Sahlholdt et al. 2019). In contrast, ages from asteroseismology can be derived with good confidence even for stars on or slightly off the zero-age main sequence (ZAMS), and are believed to be precise at the ∼10–15% level (e.g., Chaplin & Miglio 2013; Lebreton & Goupil 2014). In addition, isochrone ages may be affected by chromospheric activity for stars typically younger than 4–5 Gyr (Spina et al. 2020). Calibrating the abundance ratios against seismic ages would thus constitute a step forward, but the current samples of seismic targets with very precise abundances are still too limited in size. First attempts in that direction were made by Nissen et al. (2017, hereafter N17), who determined the abundances of 12 elements in ten solar-metallicity stars with an age of at most 7 Gyr observed by the Kepler satellite. Their results support the abundance-isochrone age relations for solar twins and analogues, although a larger scatter in the calibrations is noticeable. One of the reasons may be that their targets are significantly hotter and more evolved than the Sun.
In this paper, we present the abundances of 21 elements in a sample of 13 bright Kepler targets spanning a wide parameter range and with precise seismic ages extending up to ∼12 Gyr to further examine the performance of the abundance-age relations for stars with properties departing from solar. Our abundances probe the main nucleosynthesis production channels (iron-peak, α-, and neutron-capture) and also include lithium. These prime Kepler targets are amongst the main-sequence, solar-like stars with the best set of fundamental parameters (e.g., mass or age) known from asteroseismology. Our spectroscopic constraints complement the seismic data and will also aid further theoretical modelling.
2. Selection of targets and basic properties
Our targets were drawn from the so-called Kepler LEGACY sample made up of a total of 66 solar-like dwarfs and sub-giants (Lund et al. 2017; Silva Aguirre et al. 2017). They have all been the subject of intensive and long-duration (up to four years) observations by the Kepler satellite, and they are currently the stars with the best seismic parameters available. In particular, as discussed in Sect. 4, two independent studies derived the ages to 5–10% precision through a thorough modelling of these data. Because the amplitude of solar-like oscillations is proportional to stellar luminosity, this sample is biased towards stars that are hotter and more evolved than the Sun. However, stars at the end of the core-hydrogen burning phase that exhibit modes of mixed character were excluded (Lund et al. 2017; Silva Aguirre et al. 2017).
To adequately investigate the dependence of the abundance-age relations as a function of [Fe/H] (e.g., Feltzing et al. 2017; Skúladóttir et al. 2019; Casali et al. 2020), we first chose stars in Buchhave & Latham (2015) with metallicities deviating by more than 2σ (i.e. 0.2 dex) from solar. Our stars span a wide range of [Fe/H] values: from –0.95 to +0.41. Given the limited aperture size of the telescope, the final selection of the targets was mainly driven by signal-to-noise (S/N) considerations, and 13 stars with V < 10 mag could eventually be observed. Three near-solar-metallicity stars (KIC 6106415 and 16 Cyg AB) with high-precision abundance studies in the literature (N17 and references therein) were included in this sample for validation purposes. Except 16 Cyg B (Cochran et al. 1997), none of our targets are known to host planets, nor are they Kepler Objects of Interest (KOIs).
Very precise radial velocities (RVs) are obtained from our instrument cross-correlation (CCF) data (a G2 mask was used). A comparison with the line-of-sight (LOS) values of Lund et al. (2017) reveals clear RV changes in KIC 7871531 and KIC 12317678 for which we obtain values larger by ∼0.7 and 20.4 km s−1, respectively1. The variations for the former are small, but significant at the ∼6σ level. We therefore identify these two stars as single-lined (SB1) binaries. The binarity of our targets has recently been investigated on the basis of Apache Point Observatory Galactic Evolution Experiment (APOGEE) data. First, Price-Whelan et al. (2018) applied on the DR14 dataset a method optimised for sparse multi-epoch RV observations, but classified KIC 7871531 as single. On the other hand, El-Badry et al. (2018a) used another approach and fitted DR13 spectra to single out RV variables and stars with composite spectra (SB2). Three of our targets are flagged as single stars (KIC 3656476, KIC 6603624, and KIC 8694723), whereas three others are classified as SB2s (KIC 7871531, KIC 9965715, and KIC 12317678). No significant RV variations are apparent in our CCF data of KIC 9965715, which were acquired over two consecutive nights. The analysis of El-Badry et al. (2018a) confirms the binary nature of KIC 7871531 and KIC 12317678, although we do not detect the signature of the secondaries in our optical spectra. The APOGEE observations cover the near-IR H-band, and are therefore sensitive to cool companions. We computed the absolute magnitudes of the primaries in the V band, MV, assuming Gaia DR2 parallaxes (Gaia Collaboration 2018). Although the E(B − V) values are compatible with zero within the uncertainties, we also used reddening estimates from the 3D maps of Lallement et al. (2018)2. For the secondaries, we used the solar-metallicity isochrones of Spada et al. (2013) with a mixing-length parameter, α = 1.875, to estimate MV assuming the seismic masses and ages of the primaries discussed in Sect. 4, along with the mass ratios of El-Badry et al. (2018a). We conclude that the companions do not contribute to more than ∼10% to the flux in the optical. No visual companions have been detected in ten of our targets (including KIC 7871531 and KIC 12317678) from adaptive optics imaging in the visible (Schonhut-Stasik et al. 2017). If confirmed, the composite nature of the spectra is expected to bias our abundances by about 0.05 dex at most (El-Badry et al. 2018b). Leaving the visual binary 16 Cyg aside, to our knowledge there are no signs of binarity in the other targets. For KIC 8006161 in particular, no RV changes are detected in our data secured two days apart. This lack of variability is supported by the long-term CORAVEL monitoring of Halbwachs et al. (2018).
We used Gaia DR2 astrometric and proper motion data (Gaia Collaboration 2018), along with our RVs, to compute the space components of the targets relative to the local standard of rest (see Morel et al. 2003). Following Bensby et al. (2014), we then estimated the probability that the stars belong to a given component of the Galaxy. We find compelling evidence that KIC 8760414 displays kinematic properties consistent with a thick-disc membership (it is about seven times less likely to belong to the halo). The situation is much more ambiguous based on kinematical information alone for three other stars with an age above 8 Gyr. We find that they are about two times more likely to belong to the thin disc. However, given that kinematical properties do not allow a clean separation between the populations (e.g., Bensby et al. 2014), additional criteria – including the abundance pattern and the age – must be considered. As KIC 8760414, KIC 7970740 exhibits enhanced abundances of the α elements3 ([α/Fe]∼+0.2 dex; see Sect. 6) at a level more compatible with the high-α sequence defined by thick-disc stars in the [Fe/H]–[α/Fe] plane (e.g., Reddy et al. 2006). Its chemical properties therefore also support a thick-disc membership. We consider below that these two stars, which are the oldest and most metal-poor objects in our sample, are members of the thick disc. It can be noted that their ages (10.7 and 12.0 Gyr; Sect. 4) are consistent with this Galactic component having formed about 11 Gyr ago (Silva Aguirre et al. 2018; Miglio et al. 2021; Montalbán et al. 2020). The two other thick-disc candidates (KIC 3656476 and KIC 6603624) have large metallicities, but solar abundances of the α elements. This indicates that they do not belong to the population of metal-rich, α-enhanced stars (the so-called hαmr stars) claimed to be distinct by Adibekyan et al. (2011). It is likely that these metal-rich, solar [α/Fe] stars with ages of 8–9 Gyr migrated from the inner Galaxy (e.g., Miglio et al. 2021). The rest of the sample are clearly thin-disc members.
Our spectral analysis might be affected by chromospheric activity, although the impact is, as expected, much larger for (very) young objects (e.g., Flores et al. 2016; Yana Galarza et al. 2019). As recently shown by Spina et al. (2020), the effect on spectral lines can be significant for stars with log , where
is an activity proxy determined through measurements of the Ca II H+K emission-line fluxes. A homogeneous determination of the activity level is provided for 11 of our targets by Brewer et al. (2016). The log
values are found to tightly cluster in a range roughly solar (from –5.21 to –4.89), even for the only two targets that are much younger than the Sun (the early-F stars KIC 9965715 and KIC 12317678; see Sect. 4). A much higher activity level at the time of our observations cannot be ruled out, as illustrated by the strong activity cycle of KIC 8006161 (Kiefer et al. 2017; Karoff et al. 2018)4, but activity is not expected to be a major concern in our sample, which is largely dominated by stars with ages similar to that of the Sun or older. The basic properties of our targets are summarised in Table 1.
Basic properties of the targets.
3. Observations
The observations were secured with the échelle, fibre-fed SOPHIE spectrograph installed at the 1.93-m telescope of the Observatoire de Haute Provence (OHP, France) during the following period: 10–12 July 2018. The observations were carried out in the object+sky configuration. The spectra in high-efficiency (HE) mode have a resolving power, R, of about 40 000, and a complete wavelength coverage from about 3875 to 6945 Å. Two to three exposures were obtained for the faintest targets and co-added according to their S/N. The final S/N per pixel ranges from 213 to 350 at ∼5550 Å, with a median of 289.
For the solar spectrum to be used as reference, we averaged the six HE observations of asteroids (one of Ceres and five of Vesta) in the SOPHIE archives with an S/N above 100 (the mean S/N is 168) and a satisfactory blaze correction. Observations of a point-like source have to be preferred (e.g., Gray et al. 2000), while co-adding spectra from various reflecting bodies is not an issue (e.g., Bedell et al. 2014).
Initial processing steps were carried out with the SOPHIE reduction pipeline (DRS). All spectra were corrected for RV shifts based on the DRS CCF data. The spectrograph is simultaneously fed through two circular optical fibres, one illuminated by the target and the other one by the nearby sky. The spectrum obtained through the latter aperture was used for the sky subtraction of the stellar exposure. Finally, the mean spectra were split up into segments with a typical length of about 70 Å, and a wavelength coverage carefully chosen in such a way that the continuum level could be adequately fitted with low-order cubic spline or Legendre polynomials. To ensure the highest level of homogeneity, the same wavelength bounds and fitting functions were used for all stars. Standard tasks implemented in the IRAF5 software were used for the final reduction steps.
4. Asteroseismic data
We make use below of the age, mass, M, and surface gravity, log g, determined from asteroseismology. Two teams independently performed a seismic modelling of the oscillation frequencies and/or their separation ratios estimated from the full-length, short-cadence Kepler dataset (Lund et al. 2017). While the analysis of Creevey et al. (2017, hereafter C17) is based on a single pipeline (AMP: Asteroseismic Modeling Portal; Metcalfe et al. 2009), Silva Aguirre et al. (2017, hereafter SA17) followed another approach and used six different modelling procedures. As discussed by SA17, the results provided by the various pipelines are consistent despite the variety of codes employed or the different sensitivity of the seismic diagnostics to surface effects, for instance. We thus assume in the following the unweighted mean of all their individual values. It can be noted that C17 only used frequency ratios that are largely insensitive to the uncertain modelling of the near-surface layers. Stellar parameters of Kepler dwarfs and sub-giants determined through a detailed modelling of the full frequency dataset have to be preferred over those inferred from the global seismic quantities because they are more precise by a factor of a few. Notwithstanding, there is no evidence for the LEGACY stars for systematic differences significantly exceeding the uncertainty level (Serenelli et al. 2017).
There is an overall good agreement between the parameters of interest determined by C17 and SA17 (see Fig. 1), but a systematic difference for M and log g is apparent for the most evolved Kepler LEGACY stars with log g ≲ 4.2. It is beyond the scope of this paper to investigate in depth the origin of these slight discrepancies, which may stem from, for example, details of the fitting scheme and optimisation procedures, or choice of input physics. We simply note that the main differences between the models in SA17 and C17 are the exploration of the initial helium abundance and mixing length. Some of the pipelines of SA17 fixed the mixing length or restricted it, whereas others imposed a chemical enrichment law. Their C2kSMO method (Lebreton & Goupil 2014) is the closest to that employed by C17. An enrichment law was not imposed in C17, and high initial helium mass fractions, Yi, and low masses were occasionally obtained (e.g., KIC 8694723 and KIC 9965715 for which Yi = 0.309 and 0.310, respectively). As a test, we fixed Yi near to 0.257 and 0.267, and the new masses are revised upwards to 1.08 and 1.09 M⊙. In any case, this M − Yi correlation does not affect the age. Furthermore, the relatively small differences observed for our sample between the results obtained by C17 and SA17 have no impact on our conclusions. Given that there is also no obvious reason to prefer one set of values over the other, in the following we simply adopt the average values quoted in Table 2.
![]() |
Fig. 1. Comparison for Kepler LEGACY sample between the stellar mass (top panel), age (middle panel), and surface gravity (bottom panel) determined by C17 and SA17. The differences for a given parameter are the values of C17 minus those of SA17. The stars in our sample are shown as filled red circles. KIC 9965715 is shown for completeness as a cross. |
Summary of seismic parameters for the sample.
A different procedure was adopted for KIC 9965715 because the significantly different and more reliable spectroscopic constraints we obtain (as argued in Sect. 7.1) led us to redetermine the seismic parameters of this star. Two independent modelling approaches referred to in the following as ASTEC+ADIPLS and CLES+LOSC were performed. Full details are provided in Appendix A. The best-fit parameters obtained are shown in the bottom part of Table 2, and the straight mean values are adopted in the following. Compared to C17 and SA17, we find a mass revised upwards by ∼7% and that the star is ∼30% younger.
The comparison between the results of various pipelines suggests that the seismic ages used in our study are precise. This is supported by the fact that, for instance, inversion techniques provide an age for 16 Cyg AB (in the range 7.0–7.4 Gyr; Buldgen et al. 2016) fully compatible with the values adopted. However, we caution that the accuracy of seismic ages cannot be firmly evaluated, although they appear to be relatively robust against the choice of the input physics when detailed seismic information is available (e.g., Lebreton & Goupil 2014). It is reassuring that both C17 and SA17 closely reproduced the solar age to much better than ∼10%, but this level of accuracy may not be representative of the whole sample.
5. Abundance analysis
Except for log g (see above), the atmospheric parameters and abundances of 21 chemical elements were determined from a line-by-line differential analysis relative to the Sun, plane-parallel, 1D MARCS model atmospheres (Gustafsson et al. 2008), and the 2017 version of the line-analysis software MOOG originally developed by Sneden (1973).
The line list was taken from Reddy et al. (2003), as it was shown that the abundance analysis described below leads to an effective temperature, Teff, in good agreement with interferometric-based estimates (Kervella et al. 2017; Heiter et al. 2015) for the G2 and K1 main-sequence stars α Cen AB (Morel 2018) and the G0 sub-giant β Hyi (unpublished). In addition, [Fe/H] is within the range of the commonly accepted values for these three stars (e.g., Jofré et al. 2014).
The equivalent widths (EWs) were measured manually and assuming Gaussian profiles (multiple fits were used for well-resolved blends). Features with an unsatisfactory fit or significantly affected by telluric features based on the atlas of Hinkle et al. (2000) were discarded. Strong spectral features with RW = log (EW/λ) > –4.80 were also excluded provided that a sufficient number of lines are left after this operation. These lines are best avoided for different reasons; for example, they lie on the non-linear part of curve of growth or are formed in the upper atmosphere and particularly sensitive to chromospheric activity (Spina et al. 2020).
5.1. Determination of atmospheric parameters
A ‘constrained’ analysis, whereby the surface gravity is fixed to the very precise seismic value (Table 2), was enforced. The model parameters (Teff, microturbulence, [Fe/H], and mean abundance ratio of the α elements with respect to iron, [α/Fe]) were iteratively modified until the following conditions were simultaneously fulfilled: (1) the Fe I abundances exhibit no trend with RW; (2) the mean abundances derived from the Fe I and Fe II lines are identical; and (3) [Fe/H] and [α/Fe] are consistent with the model values. The α-element abundance of the model varied depending on [Fe/H] following Gustafsson et al. (2008). For instance, [α/Fe] = +0.2 for [Fe/H] = –0.4. Although excitation balance of iron is not formally fulfilled in our constrained analysis, it can be noted that the slope between the Fe I abundances and the lower excitation potential (LEP) deviates from zero by less than ∼2σ in our sample. The only clear exception is KIC 6603624 for which it appears that the adopted seismic log g is lower by ∼0.3 dex than the value that would be inferred from spectroscopy. For the solar analysis, Teff and log g were held fixed to 5777 K and 4.44, respectively, whereas the microturbulence, ξ, was left as a free parameter.
Following N17, we corrected the Fe I abundances from departures from local thermodynamic equilibrium (LTE). We made use of Spectrum Tools6 to interpolate the non-LTE corrections as a function of the stellar parameters. The calculations are described in Bergemann et al. (2012) and are available for a representative set of lines (about two thirds of the total). This was deemed unnecessary for the Fe II-based abundances because non-LTE effects are negligible. Even for Fe I, the differential corrections are small (less than 0.02 dex), and taking them into account has very little impact on the parameters derived through iron ionisation balance.
Freezing log g to the seismic value has a much more profound effect on our results. We find that relaxing this constraint would result in log g, Teff and [Fe/H] values larger on average by 0.11 dex, 50 K, and 0.03 dex, respectively. Our choice of a constrained analysis is motivated by the fact that severe discrepancies between the spectroscopic log g and the reference value are occasionally observed on a star-to-star basis (e.g., up to 0.5 dex in the extreme case of KIC 12317678). There is still no consensus among the community as to whether such an analysis has to be preferred for seismic targets (e.g., Doyle et al. 2017). A constrained analysis is undoubtedly more precise, but whether it is also more accurate is unclear. Our study does not shed new light on this issue. We simply note, as argued in Sect. 7.2, that a comparison with the handful of long-baseline interferometric measurements available lends some support to the cooler Teff scale resulting from the constrained analysis. This is consistent with the recent claim based on a larger sample of interferometric benchmark targets that a constrained analysis is more accurate, especially for F-type stars (Gent et al., in prep.).
5.2. Determination of chemical abundances
Our line list for elements other than iron was taken from Reddy et al. (2003). However, for some key elements (Mg, Al, and Zn) we substituted this line list with that of Meléndez et al. (2014) because not enough lines are included. This was also done for Co because not all lines in Reddy et al. (2003) have hyperfine structure (HFS) data available. HFS and isotopic splitting were taken into account for Sc, V, Mn, Co, and Cu using atomic data from the Kurucz database7 and assuming the Cu isotopic ratio of Rosman & Taylor (1998). The barium abundance is solely based on Ba IIλ5853 for which HFS and isotopic splitting can safely be neglected (e.g., Mashonkina et al. 1999). It is supported by our own tests on the Sun using the atomic data of Prochaska et al. (2000).
We do not further discuss the abundances derived for Si II and Cr II because they are based on a single feature, whereas those of the corresponding neutral species rely on up to six lines. However, despite the uncertainty plaguing the Si II- and Cr II-based abundances, ionisation balance is fulfilled for these two elements in the vast majority of cases (Fig. 2).
![]() |
Fig. 2. Comparison for Si (top panel) and Cr (bottom panel) between the mean abundances yielded by the neutral and singly-ionised ion. The mean difference is given in each panel (the number in brackets is the number of stars the calculation is based on). A representative error bar is also shown. |
The determination of the lithium abundance relied on a spectral synthesis of Li Iλ6708 (for details, see Morel et al. 2014). This line is confidently detected in seven targets; upper limits are provided for the others. The rotational and macroturbulent broadening values were taken from Brewer et al. (2016). For KIC 8694723 and KIC 8760414, which are not included in this study, we adopted the v sin i values of Bruntt et al. (2012) and computed the macroturbulence from the calibrations as a function of Teff of Bruntt et al. (2010). Tests using the broadening parameters from the works above, but also from Doyle et al. (2014), show that the exact choice of these quantities is irrelevant. For the Sun, we obtain an absolute lithium abundance, log ϵ⊙(Li) = +1.08, which is almost identical to the recommended photospheric solar value (e.g., Asplund et al. 2009). For the five stars with a firm detection in common with Bruntt et al. (2012) and Tucci Maia et al. (2019), there is a good agreement (within 0.1 dex) without any evidence for a systematic offset (see also Sect. 8). Our non-detection in five stars is also consistent with the findings of Bruntt et al. (2012) and Beck et al. (2017).
For a number of reasons, we did not attempt to correct our abundances for the combined effect of departures from LTE and time-dependent convection that would in principle require us to drop the rudimentary assumption of 1D. First, performing full 3D, non-LTE radiative transfer calculations is extremely tedious. As a result, theoretical predictions are currently only available for a few elements and very often do not span the whole parameter space of our sample (Amarsi et al. 2019a; Amarsi & Asplund 2017; Nordlander & Lind 2017; Gallagher et al. 2020; Bergemann et al. 2019; Mott et al. 2020). To our knowledge, spectral line formation under non-LTE has not even been investigated for some key elements in the context of our study (e.g., Y). Correcting for either non-LTE or 3D effects is not necessarily recommended because the corrections might be of opposite sign and counterbalance each other (e.g., Amarsi et al. 2019b). As a result, there is no guarantee that the corrected abundances are more accurate. Second, and more importantly, the abundance-age relationships that are the main focus of this paper are built on the assumptions of 1D and LTE (e.g., Delgado Mena et al. 2017). Should our abundance data be used in a wider context, we caution that care should be exercised for some elements that are known to be particularly sensitive to non-LTE and/or 3D effects, especially in the metal-poor regime (e.g., Mn; Bergemann et al. 2019).
6. Results
The atmospheric parameters and chemical abundances are given in Tables 1 and B.1, respectively. The random uncertainties were computed following common practice (see Morel 2018). For the abundances, the line-to-line scatter, σint, and the uncertainties arising from errors in the stellar parameters were added in quadrature. For iron, σint typically amounts to 0.04 dex. For the elements with a single diagnostic line, we assumed σint = 0.05 dex. The uncertainties are much larger for the three metal-poor, hottest stars because of the paucity and weakness of the diagnostic Fe I lines.
We explored to what extent the choice of another family of model atmospheres would affect our results. Namely, we repeated the analysis using ATLAS9 Kurucz models (Castelli & Kurucz 2003) for five representative stars spanning the parameter space of our sample in terms of Teff, log g, and [Fe/H]: KIC 7871531, KIC 8006161, KIC 8694723, KIC 8760414, and 16 Cyg A. We find in all cases very small Teff differences not exceeding 5 K. Although the use of Kurucz models leads to systematically lower microturbulent velocities (on average by 0.035 km s−1), this is paralleled in the Sun. As a net result, we find that all the abundance ratios differ by a negligible amount (virtually in all cases below 0.01 dex).
7. Validation of stellar parameters
7.1. Comparison with other spectroscopic studies
A comparison with the results of Buchhave & Latham (2015) obtained with the Stellar Parameters Classification (SPC) tool is particularly relevant because, with the exception of KIC 9965715 and 16 Cyg AB, their Teff and [Fe/H] values were adopted by C17 and SA17 to perform their seismic modelling. A constrained analysis was also enforced by adopting the seismic gravities of Chaplin et al. (2014) as priors. The results of Buchhave & Latham (2015) are based on data obtained with the TRES échelle spectrograph installed on the 1.5-m telescope of the Fred Lawrence Whipple Observatory (Mt. Hopkins, Arizona). The spectra have similar characteristics to ours (R ∼ 44 000 and covering 3800–9100 Å), but the S/N is modest.8 Despite this difference in terms of data quality, there is a very close agreement between the stellar parameters without any evidence for systematic discrepancies (Fig. 3).
![]() |
Fig. 3. Upper panels: differences (this study minus literature) between our Teff and [Fe/H] values and those adopted by C17 and SA17 for their seismic modelling. Bottom panels: same as upper panels, but for the SDSS and ASPCAP non-seismic constraints adopted by Serenelli et al. (2017). |
In sharp contrast, our Teff and [Fe/H] values for KIC 9965715 are discrepant with those of a different origin used by C17 and SA17 to model the Kepler data.9 We find that the star is less metal poor and much hotter, which is much more in line with an F2 V classification and agrees with the conclusions of Molenda-Żakowicz et al. (2013) and Brewer et al. (2016). On the other hand, Compton et al. (2018) reported problems when trying to model this star with the Teff adopted by C17 and SA17.
Overall, one can note a poorer level of agreement between our values and those quoted in the APOCASK10 catalogue (Serenelli et al. 2017), either derived from the analysis of near-IR APOGEE DR13 spectra with the APOGEE Stellar Parameters and Chemical Abundances Pipeline (ASPCAP) or, to a lesser extent, from Sloan Digital Sky Survey (SDSS) griz photometry. The ASPCAP Teff scale is indeed known to be too cool for dwarfs and sub-giants (e.g., Serenelli et al. 2017; Martinez et al. 2019), which likely explains the [Fe/H] discrepancy also seen in Fig. 3. Generally speaking, it is commonplace to find noticeable differences between [Fe/H] determinations based on optical or near-IR spectra, or even between the results of different pipelines applied to APOGEE data (Sarmento et al. 2020).
In Table 3, we present a global comparison with other literature studies. Molenda-Żakowicz et al. (2013) and Furlan et al. (2018) used various methods. For the former work, there are discrepancies amounting to almost 0.5 dex between the spectroscopic and seismic gravities. One can note that a cooler/hotter Teff scale is often associated with lower/larger surface gravities on average. This is a well-known consequence of the degeneracy in the determination of these two quantities from spectroscopy (e.g., Torres et al. 2012) that implies that the Teff scale is sensitive to the presence of systematic biases in the determination of log g. The Teff and log g differences indeed appear to be positively correlated when comparing our results to those of Molenda-Żakowicz et al. (2013) or Brewer et al. (2016). This is also true for the SPC and combined results of Furlan et al. (2018). We find a good agreement with our Teff scale when allowance is made for the fact that the spectroscopic gravities for the unconstrained studies above deviate from the very precise seismic estimates. Nonetheless, there is evidence that our metallicities are, on average, larger by ∼0.04 dex than those of previous studies. Although generally speaking there is a fair consistency between the various sets of parameters, a closer look on a star-to-star basis reveals that the agreement is very good for solar analogues, but it worsens considerably for hotter and/or more metal-poor stars. This is not surprising, because a robust determination is much more difficult to achieve.
Comparison with stellar parameters in the literature.
7.2. Comparison with interferometric Teff scale
Four stars in our sample have a nearly model-independent Teff determined from a combination of absolute flux and limb-darkened CHARA angular diameter measurements: KIC 6106415 and KIC 8006161 (Huber et al. 2012), as well as 16 Cyg AB (White et al. 2013). We do not discuss the CHARA observations of 16 Cyg AB by Boyajian et al. (2013) because they were carried out without the PAVO beam combiner whose use has been claimed to yield more reliable measurements for stars that are not well resolved (e.g., Casagrande et al. 2014; Karovicova et al. 2018). The Teff scale of Boyajian et al. (2013) and White et al. (2013) appear indeed discrepant, with the former being ∼100 K cooler.
To add a comparison point at higher effective temperatures, we analysed the Gaia benchmark β Vir (F9 V) following exactly the same steps as the main targets. The reference Teff (6083 ± 41 K) and log g (4.10 ± 0.02) values were taken from Heiter et al. (2015), and primarily rely on the interferometric and seismic observations of North et al. (2009) and Carrier et al. (2005), respectively. Our analysis was based on a high-quality spectrum (S/N ∼ 295) retrieved from the SOPHIE archives. Because no HE spectra are available, we had to rely on a high-resolution (HR) spectrum with R ∼ 75 000. However, our results are expected to be representative of those obtained with HE spectra. To ensure consistency, a HR asteroid spectrum was also necessary. We made use of the co-addition of the ten best-quality exposures in the archives, which results in a mean S/N of 330. For β Vir, we obtain the following: Teff = 6145 ± 35 K and [Fe/H] = +0.16 ± 0.03. The recommended metallicity value assuming LTE quoted by Jofré et al. (2014) is [Fe/H] = +0.21 ± 0.07. Incidentally, an unconstrained analysis yields a log g lower by 0.01 dex only. The Teff and [Fe/H] values are hence very similar.
Figure 4 shows a comparison between our effective temperatures and those derived from absolute flux and long-baseline interferometric measurements. The values obtained by N17 from iron ionisation balance are also added. Except for 16 Cyg AB, there is indication that the spectroscopic values (either ours or those from N17 for KIC 6106415) are systematically larger. There have been previous reports of systematically hotter Teff scales determined from spectroscopy (e.g., Boyajian et al. 2013; Brewer et al. 2016). Definitive conclusions can hardly be drawn considering the paucity of data points and the recognition that, as discussed above, there is some concern about the robustness of the interferometric data. In particular, we find that the interferometric radii of KIC 6106415 and KIC 8006161 tend to be larger than the seismic estimates of both C17 and SA17. The difference is particularly outstanding for KIC 6106415, and was ascribed by SA17 to the fact that it is poorly spatially resolved and possibly affected by calibration problems. The limb-darkened angular diameters, θLD, of 16 Cyg AB are larger than those of KIC 6106415 and KIC 8006161 (θLD ∼ 0.5 vs 0.3 mas), while β Vir is well resolved (θLD ∼ 1.45 mas). Order-of-magnitude estimates show that the Teff discrepancy observed in the two stars that are barely resolved may solely be explained by an overestimation of their angular diameters. Reddening should not be an issue given that these stars are all located within 50 pc and comfortably lie inside the Local Bubble. In any case, the cooler Teff scale obtained when fixing log g (see Fig. 4) seems to support the choice made in Sect. 5.1 of a constrained analysis (see also Gent et al., in prep.).
![]() |
Fig. 4. Comparison between our ionisation temperatures (constrained and unconstrained shown as filled and open black circles, respectively) and those in the literature derived from interferometric and absolute flux measurements (filled red squares). The blue crosses show the results of N17 determined from iron ionisation balance. The stars are ordered as a function of increasing Teff. |
8. Validation of chemical abundances
The bright, G1.5 V + G3 V binary system 16 Cyg AB is a target of great interest in the context of asteroseismology (e.g., Buldgen et al. 2016; Verma et al. 2014; Bazot et al. 2019; Farnir et al. 2020), and has been the subject of numerous high-precision abundance studies. Interestingly, there is evidence of a small metallicity difference between the two components, which might be attributable to the formation and/or ingestion of planets. The existence of a slight enhancement in metals (Δ[Fe/H] ∼ 0.03 dex) in the primary was first convincingly established by Laws & Gonzalez (2001) and is no longer disputed (Tucci Maia et al. 2019 and references therein). Detecting such a subtle difference represents a stringent test of the precision of our results. In Appendix C, we compare our results to those of N1711 and Tucci Maia et al. (2019), which are based on data of exceptional quality (R = 115 000 − 160 000 and S/N ∼ 800). A remarkable agreement is found in both cases, with abundances with respect to hydrogen differing by only ∼0.01 dex on average. Furthermore, contrary to previous studies that failed to reveal this peculiarity (Schuler et al. 2011; Takeda 2005), our results support the slightly more metal-rich nature of the primary. We obtain a weighted mean difference of the abundances with respect to hydrogen of ∼0.021 dex, which – albeit very small – is highly significant (see Appendix C for further details).
The G0 V star KIC 6106415 was studied by N17. We find an ionisation Teff lower by ∼35 K. A comparison between the abundance data is shown in Fig. 5, as a function of the 50% condensation temperature for a Solar System composition gas (Tc; Lodders 2003). Although there is a slight, systematic offset (∼0.02 dex) between the [X/H] values, the scatter is very small (∼0.015 dex). The zero-point offset affecting the absolute abundances may have several causes (e.g., continuum placement). The largest discrepancy (∼0.05 dex) is found for Zn, which has uncertain abundances in both studies (see also Nissen et al. 2020). The agreement is remarkable when the abundances with respect to iron are considered: ⟨Δ[X/Fe]⟩ (this study minus N17) = –0.007 ± 0.015 dex.
![]() |
Fig. 5. Upper panel: abundance pattern of KIC 6106415 with respect to hydrogen, [X/H], as a function of Tc. Our [X/H] results and those of N17 are shown as filled and open circles, respectively. A dotted, horizontal line is drawn at our [Fe/H] value. The solid line shows our weighted, linear fit of [X/H] as a function of Tc. The fit obtained by N17 is overplotted as a dashed line. Lower panel: differences (our values minus theirs) with respect to N17, δ[X/H]. The average δ[X/H] value is given (the number of elements in common is indicated in brackets). To guide the eye, a dotted line is drawn at δ[X/H] = 0. |
Although the tests described above suggest that our abundance results are precise at the 0.02–0.03 dex level for early G dwarfs, we warn the reader that this conclusion cannot be extended to the stars in our sample with properties significantly different from solar. It is because the effects of various physical phenomena (departures from LTE, convection inhomogeneities, atomic diffusion) no longer cancel out to the first order through our differential analysis.
For a more general validation of the elemental abundances, we consider the comprehensive studies of the Kepler targets undertaken by Bruntt et al. (2012) and Brewer et al. (2016). The surface gravity was held fixed to the seismic value in the former work. Both analyses are based on spectral synthesis performed with (semi-)automatic tools. The results are compared to ours in Fig. 6. Although some outliers are evident (e.g., KIC 12317678 when compared to Brewer et al. 2016), there is a reasonable agreement with only a slight offset in general for the abundance ratios with respect to hydrogen (of the order of 0.04 dex on average when considering the median differences). The most conspicuous difference is the lower carbon abundances we obtain compared to Bruntt et al. (2012), but their uncertainties are quite large for some stars (up to 0.15 dex). However, whether or not we exclude this element, we tend to find a slightly better agreement with Brewer et al. (2016).
![]() |
Fig. 6. Comparison between our abundance ratios and those of Bruntt et al. (2012) and Brewer et al. (2016). The abundance differences are those of this study minus those of the literature. The vertical dashed line connects the extreme values. The box covers the first to third quartile of the data, while the thick horizontal line inside the box shows the median. As recommended by Bruntt et al. (2012), we chose their Fe II-based abundances for [Fe/H] (but used the Mg I, Si I, Ti I, and Cr I abundances). For KIC 5184732, KIC 7970740, and KIC 8006161, we chose the values of Brewer et al. (2016) corresponding to the spectrum with the highest S/N. The number of stars the calculation is based on is indicated for each abundance ratio. The red symbols show the abundance ratios used as age indicators (Sect. 9). |
9. Discussion
9.1. Age inferences from abundance data for the LEGACY sample
Several linear relations directly linking isochrone ages and abundance ratios have been proposed for solar twins and analogues (e.g., Bedell et al. 2018; Nissen 2015; Tucci Maia et al. 2016). Such linear relationships may be an oversimplification (Spina et al. 2016, 2018). It is also becoming increasingly clear that other variables must be taken into account in order to apply them to stars with parameters strongly departing from the solar values (e.g., Feltzing et al. 2017). Recently, Delgado Mena et al. (2019, hereafter DM19) used a large sample of main-sequence FGK stars observed as part of the HARPS-Guaranteed Time Observations (HARPS-GTO) programme to empirically explore the dependency of the abundance-isochrone age relationships as a function of Teff, metallicity, and mass. They found that several abundance ratios, when combined with one or two stellar parameters, can be used to infer stellar ages to within 1.5–2 Gyr. They proposed three types of relationships: 1D (age vs abundance ratio), 2D (age vs abundance ratio and either Teff, [Fe/H], or M), and finally 3D (age vs abundance ratio and two of either Teff, [Fe/H], and M). The main rationale behind including a Teff and M dependency is to capture the effect of atomic diffusion, although it can be noted that using metal abundance ratios largely reduces its importance (Dotter et al. 2017). They concluded that 2D relations perform better compared to those in 1D (especially when [Fe/H] or M are folded in) and that 3D relations do not lead to a major additional improvement. Furthermore, it was found that 3D relations yield similar results regardless of the choice of the two stellar parameters. It should be noted that, contrary to some studies (e.g., Nissen 2016), thin- and thick-disc stars were not separated when constructing these calibrations. Thick-disc stars constitute about 20% of the sub-sample of 354 stars with the most precise isochrone ages (uncertainty below 1.5 Gyr) built for that purpose. Therefore, we apply the calibrations to our whole sample in the following. We discuss below the ages we obtain from these relationships and compare them to the seismic estimates.
Several abundance ratios that can be used as ‘chemical clocks’ were proposed by DM19, among which fourteen can be computed from our data. Although we refrain from discussing the time evolution of other abundance ratios in our small dataset, we briefly comment on the behaviour of the predominantly slow (s-) neutron-capture elements Sr and Y. As previously found (e.g., Spina et al. 2018), a smooth decline of [Sr/Fe] and [Y/Fe] with increasing age is seen. However, there is some evidence for an upturn for the two thick-disc stars with an age above ∼10 Gyr (Fig. 7). This is consistent with a picture in which there was a more vigorous production, during the formation of the thick disc, of s-elements from low-mass asymptotic giant branch (AGB) stars relative to iron from Type Ia supernovae (e.g., Battistini & Bensby 2016). This is apparently not seen in Fig. 7 for Ba, while the only Ce abundance at old ages is very uncertain. See DM19 for a discussion of the behaviours of the light and heavy s-elements.
![]() |
Fig. 7. Abundances of neutron-capture elements relative to iron as a function of seismic age. The thin- and thick-disc stars are shown as circles and squares, respectively. The solar age is shown as a vertical dotted line. |
In Fig. 8, we show how some abundance ratios selected by DM19 vary as a function of seismic age. The data are colour-coded as a function of Teff, [Fe/H], and M. Previous 1D relations proposed in the literature are also overplotted (DM19; Nissen 2016; Nissen et al. 2020; Tucci Maia et al. 2016; Spina et al. 2018; Bedell et al. 2018; Titarenko et al. 2019). We recall that, with the exception of DM19 and Titarenko et al. (2019), these relationships only apply to solar analogues and are not expected to provide a good fit to our data. Also, some of them may not be valid for old or thick-disc stars. The deviations with respect to 1D calibrations may be ascribed to stellar parameters differing from the solar values. A good example is [Y/Si], where the systematic deviations clearly noticeable for the old stars are largely removed when using 2D or 3D relations. A more general discussion is provided below.
![]() |
Fig. 8. Left panels: variation of abundance indicators, as a function of seismic age. The data are colour-coded as a function of Teff. The thin- and thick-disc stars are shown as circles and squares, respectively. The solar age is shown as a vertical dotted line. Middle and right panels: same as left panels, but colour-coded as a function of [Fe/H] and M, respectively. Linear 1D relations in the literature linking abundance ratios and ages are overplotted as dashed lines: Bedell et al. (2018, red), DM19 (blue), Nissen (2016) and Nissen et al. (2020, cyan), Spina et al. (2018, magenta), Titarenko et al. (2019, yellow), and Tucci Maia et al. (2016, green). |
We compare in Table 4 for the whole sample the seismic ages and those obtained for each type of relation (e.g., 1D) and abundance indicator (e.g., [Y/Mg]). Both the uncertainties in the abundance ratios and in the calibrations were propagated to the ages. The best abundance indicators proposed by DM19 differ depending on the parametrisation. In this respect, it is interesting to note that [Y/Al], which was often discussed in the literature in the context of solar analogues (e.g., Nissen 2016), is not among the ratios displaying the best 1D relation when considering the much more heterogeneous sample of DM19.
Unweighted mean age deviation with respect to seismic estimates (abundance-based minus seismic) for each abundance indicator and relation (the number in brackets is the number of stars the calculation is based on).
The young, F-star KIC 12317678 falls close to the upper Teff boundary of the relationships of DM19 and has abundance ratios that are particularly uncertain. However, we find that including it in our sample does not bias our results. For instance, removing it leads to mean deviations between the seismic- and abundance-based ages for the whole sample that differ very little (in the ± 0.2 Gyr range only depending on the type of relation). This is mainly because unphysical (negative) ages are often obtained – especially when they are based on the high Y abundance – and evidently ignored. Negative ages were also obtained in two cases for another young star, KIC 9965715.
As seen in Table 4, the standard deviation of the distributions of age differences is typically ∼1.5–2 Gyr, which is comparable to the intrinsic precision of the calibrations of DM19. To put these results in perspective, internal errors of the order of 0.5–1 Gyr are common for solar analogues (e.g., Tucci Maia et al. 2016; Spina et al. 2018). Generally speaking, determining ages from a single abundance indicator, whatever the relation used, is not recommended because it is prone to errors amounting to up to several Gyr. Some abundance-based ages (e.g., using [Ti/Fe]) appear to be less precise than others. There is tentative evidence that the age scatter decreases as the intrinsic quality of the 1D and 2D Teff relations increases (as parametrised by the goodness-of-fit measure, adj-R2; see DM19 for definition). It is not seen for the other relations, but the adj-R2 range is much smaller. The scatter between the abundance- and seismic-based ages is typically reduced by ∼20% when averaging the results of all abundance ratios. As may be expected, combining the results of many chemical clocks leads to more precise ages, but the gain is relatively modest.
Ages estimated from pairs of elements with significantly different condensation temperatures, Tc, may be more sensitive to star-to-star differences in their abundance-Tc trend (N17; Nissen & Gustafsson 2018). Namely, stars with the same age might have different volatile-to-refractory abundance ratios (e.g., Biazzo et al. 2015). In our case, one might thus expect the two abundance indicators involving the volatile element Zn ([Sr/Zn] and [Y/Zn]) to yield less precise ages. However, it is not borne out by our results (Table 4).
There is clear evidence that the calibrations provide systematically younger ages compared to asteroseismology by ∼0.7 Gyr on average, irrespective of the relationship used. To assess the robustness of this result, we computed the ages of KIC 6106415 and 16 Cyg AB using the abundance data of N17 and Tucci Maia et al. (2019). As seen in Fig. 9, the discrepancy persists. As an additional test, we redetermined the abundances of the metal-poor, G0 sub-giant KIC 8694723 adopting the line lists used to construct the calibrations of DM19. They were taken from Neves et al. (2009) for Mg, Al, Si, and Ti (with a few lines excluded, as described in Adibekyan et al. 2012), and from Delgado Mena et al. (2017) for Zn, Sr, and Y. We did not choose a solar analogue because non-LTE and 3D effects may be efficiently erased in that case through a differential analysis whatever the line list. On the contrary, the combined effect might be quantitatively different depending on the choice of the diagnostic lines when the reference star occupies another region of the parameter space. We obtain ages on average ∼0.8 Gyr older. The fact that it apparently solves the discrepancy mentioned above might be accidental: much more extensive tests are needed. We also used our atmospheric parameters, while systematic differences (e.g., in the Teff scale) are also likely to play a role. In any case, it shows that age deviations of the right order of magnitude might be ascribed to small (≲0.05 dex) zero-point abundance offsets. It is particularly true for Sr and Y because these two elements enter most calibrations. However, systematic differences between various sets of theoretical isochrones are well known. DM19 made use of PARSEC isochrones (Bressan et al. 2012) to infer the ages through Bayesian inference. Adopting another set of evolutionary models may be enlightening in this respect.
![]() |
Fig. 9. Mean ages of KIC 6106415 and 16 Cyg AB obtained using our abundance data (circles), those of N17 (squares), and those of Tucci Maia et al. (2019, triangles). The shaded horizontal stripes show the seismic ages, along with their ±1σ uncertainty. |
Figure 10 illustrates the age deviations, ΔA, as a function of Teff, [Fe/H], and M. The Sun is overplotted to ensure that its age is correctly reproduced and that no systematic age offsets are present. The size of the error bar in this particular case also provides an indication of the age uncertainty arising from the calibrations alone. There is no evidence in Fig. 10 that the thick-disc stars are outliers. We define a quantity, δA, which is a measure of the improvement brought about by using 2D or 3D relations instead of 1D ones. A positive value implies that the age is closer by δA Gyr to the reference seismic value when using a 2D or 3D parametrisation (and vice versa). The relative performance of the various relations is better appraised when examining the δA values for the four abundance indicators common to all relations: [Sr/Mg], [Sr/Ti], [Y/Mg], and [Y/Zn] (Fig. 11). The use of more sophisticated relations generally leads to a better agreement with the seismic ages (δA positive and reaching up to +0.5 Gyr), although the effect is small (see also Fig.13 of DM19).
![]() |
Fig. 10. Mean age deviation, ΔA, with respect to the seismic estimate (abundance-based minus seismic) when considering all the abundance indicators, as a function of Teff, [Fe/H], and M. The quantity δA is a measure of the improvement brought about by using 2D or 3D instead of 1D relations (see Sect. 9.1). The thin- and thick-disc stars are shown as filled and open circles, respectively. The Sun is shown as an open triangle. The rightmost panels show the breakdown of the ΔA and δA values. |
However, it is conceivable that a much more significant improvement is achieved in some regions of the parameter space. For instance, the age derived for the most metal-poor LEGACY star, KIC 8760414, is revised upwards by almost 4 Gyr using the relations making use of [Fe/H], and is no longer severely underestimated compared to the seismic result. DM19 cautioned that their relations are not robust for [Fe/H] ≲ –0.8. Whether the improvement observed for KIC 8760414 is only a fluke or indicates that the relations are still of some applicability for metal-poor stars must be investigated further. The improvement arising from the use of 2D or 3D relations might be buried in the noise for the other stars in our sample with parameters closer to solar.
9.2. Age inferences from abundance data for an extended sample
To further examine the performance of the abundance-age relations proposed by DM19, we considered an extended sample of Kepler dwarfs and sub-giants with a homogeneous determination in the literature of both the abundances and the seismic properties. Namely, we selected stars in common from Brewer et al. (2016) and Serenelli et al. (2017). Abundances of the key elements Sr and Y are not provided by Bruntt et al. (2012). We note that Brewer et al. (2016) empirically corrected their abundance dataset for trends as a function of Teff. In the following, we use the seismic results of Serenelli et al. (2017) based on SDSS data. We ignored stars with Teff taken from Brewer et al. (2016) below 5300 K because they fall outside the validity domain of the calibrations of DM19. Four stars exhibit anomalously high abundance ratios involving yttrium: we excluded KIC 9025370 because it is an SB2 (N17), and KIC 11026764 because significantly lower values were obtained by Metcalfe et al. (2010). We retained the remaining two stars: KIC 3733735 and KIC 9812850. No additional information about a general enhancement of the s-process elements12 is available, but, in any case, the ages inferred from the high Y abundances are all negative and therefore not considered further. We count a final total of 63 stars. It should be noted that they are on average younger (∼3 Gyr) and more evolved (log g values down to 3.5) than our LEGACY sample. The former feature likely explains the large proportion (∼20%) of unphysical (negative) ages obtained. Following Battistini & Bensby (2016), we simply flagged stars older than 8 Gyr as thick-disc members. There are only four stars fulfilling this criterion.
Three limitations must be kept in mind. First, only seven abundance ratios could be computed from the data of Brewer et al. (2016), because Zn and Sr are not available. Second, the stellar parameters of Serenelli et al. (2017) are based on global seismic quantities and are expected to be a factor of 2–3 less precise than those based on the modelling of the full set of frequencies or their separation ratios. Finally, we argue in Sect. 7.1 that for all but one LEGACY stars, our spectroscopic parameters are consistent with those adopted by C17 and SA17 for their seismic modelling. For the remaining star, KIC 9965715, we performed a new seismic analysis using our Teff and [Fe/H] values (Appendix A). In contrast, here there are systematic differences between the classical parameters adopted for the determinations of the chemical abundances and ages (see Fig. 12). For the stars selected, we find that the Teff scale of Brewer et al. (2016) is on average ∼120 K cooler than the SDSS one. For completeness, we also show the comparison with the ASPCAP Teff scale. It supports our previous conclusion that these Teff values are too low. In addition, the metallicities of Brewer et al. (2016) are ∼0.11 dex larger, which is also fully consistent with our finding that the ASPCAP values may be underestimated by this amount
![]() |
Fig. 12. Differences between stellar parameters of Brewer et al. (2016) and Serenelli et al. (2017). The mean values for the SDSS and ASPCAP Teff scales are given in each panel. We assumed the ASPCAP [M/H] values that are the overall scaled solar abundances. |
(Sect. 7.1). In contrast, there is a satisfactory agreement between the spectroscopic and seismic gravities, as indeed anticipated (see Brewer et al. 2015).
Figure 13 shows the abundance-age trends for our data and the sub-sample of 63 stars from Brewer et al. (2016). As shown in Fig. 14, the analysis of this much larger sample broadly confirms the conclusions drawn above for the LEGACY stars. Namely, the scatter of the deviations between the abundance- and seismic-based ages is of the order of 2 Gyr and an improvement of up to ∼0.5 Gyr occurs with 2D or 3D relations. Noteworthy is the fact that the 2D Teff and 3D relationships (the latter also involving Teff) quite efficiently erase the clear trend between the age difference and Teff. The lack of improvement when incorporating [Fe/H] as an additional variable might be due to the near-solar metallicity range. We find that negative δA values for the 2D [Fe/H] calibration are generally associated with the most evolved stars (log g ≲ 4.1). One may notice two sequences running almost parallel to the x-axis at ΔA∼ < SPSDOUBLEDOLLAR > −2 and +2 Gyr (best seen for the 2D mass relation) in the panels as a function of stellar mass in Fig. 14. If real, there is no clear explanation for this behaviour. Once again, the seismic ages appear overall systematically larger by ∼0.6 Gyr.
![]() |
Fig. 13. Abundance-age trends for our data (filled symbols) and those of Brewer et al. (2016, open symbols). The ages were derived by different means: detailed modelling of the oscillation frequencies and from global seismic quantities, respectively. The very high [Y/Al] abundance ratio (∼+1.35) obtained by Brewer et al. (2016) for KIC 3733735 is off scale. The thin- and thick-disc stars are shown as circles and squares, respectively. Linear 1D relations in the literature are overplotted as dashed lines: Bedell et al. (2018, red), DM19 (blue), Nissen (2016) and Nissen et al. (2020, cyan), Spina et al. (2018, magenta), Titarenko et al. (2019, yellow), and Tucci Maia et al. (2016, green). The solar age is shown as a vertical dotted line. |
10. Summary and conclusions
10.1. Implication for further seismic modelling
Our non-seismic parameters (especially Teff, [Fe/H], and [α/Fe]) complement the superb Kepler oscillation data obtained for these bright stars and will aid further seismic modelling. These are important observational constraints that are necessary to exploiting the full potential of asteroseismology (see e.g., Chaplin et al. 2014; Serenelli et al. 2017; Creevey et al. 2012; Valle et al. 2018). We classified two of our targets with an age exceeding 8 Gyr as thick-disc members mostly based on chemodynamical arguments. Their chemical pattern significantly different from the solar mixture (in particular an enhancement of the α elements) must be taken into account to achieve sensible seismic inferences (e.g., Ge et al. 2015; Li et al. 2020). Our results can also be useful in other more specific contexts, such as constraining the Galactic helium enrichment law (e.g., Verma et al. 2019) or improving the modelling of the outermost stellar layers (e.g., Compton et al. 2018). Two Li-detected stars, KIC 5184732 and KIC 8694723, also have a determination of their rotational period from Kepler light curves (García et al. 2014; Karoff et al. 2013), which paves the way for a detailed theoretical modelling of their interior.
We find evidence that a few of our targets are associated in a binary system with a large mass ratio. This suggests the existence in the Kepler LEGACY sample of a sizeable fraction of binaries. This is supported by the fact that N17 detected two SB2s in their sample of ten stars from snapshot observations only. It would not be surprising considering the ubiquity of companions in solar-like stars (e.g., Duquennoy & Mayor 1991). This aspect certainly deserves further investigation (i.e. a dedicated RV monitoring).
As guidance for further ensemble seismic modelling of Kepler main-sequence stars, our comparison with the results of (semi-)automatic methods suggests that the study of Buchhave & Latham (2015) is a valuable source of effective temperatures and metallicities. For the elemental abundances, the works of Bruntt et al. (2012) and Brewer et al. (2016) appear to be of similar quality for most elements. However, our results tend to be in closer agreement with those of the latter study. We confirm that the ASPCAP effective temperatures used for the APOKASC catalogue are systematically too cool for dwarfs. This is not completely unexpected given that the pipeline is naturally optimised for red giants (García Pérez et al. 2016). The consequence is that a grid-based modelling that makes use of the SDSS photometric Teff scale may be more reliable for relatively unevolved Kepler targets, as already claimed by Serenelli et al. (2017). However, one can see, in Fig. 12, the existence of a systematic discrepancy between the SDSS Teff values and those of Brewer et al. (2016). We find evidence in our sample that the APOKASC metallicities for dwarfs (solely based on ASPCAP) are typically underestimated by ∼0.1 dex. A similar conclusion is reached when considering the data of Brewer et al. (2016). Correcting for such a bias would lead to a fractional improvement of a few percent in the stellar properties quoted in the APOKASC catalogue (e.g., ∼5% in age; Serenelli et al. 2017). As illustrated by our new seismic analysis of KIC 9965715, even larger changes can be expected if the quality of the classical parameters is mediocre.
10.2. Stellar ages from chemical clocks beyond solar analogues
We find that the quality of stellar ages empirically determined from abundance data is currently limited by the precision of the calibrations, which is typically 1.5–2 Gyr when venturing away from stars closely resembling our Sun. Constructing more complex relationships that take into account a dependency on one or two stellar parameters does improve the situation, if only slightly. For instance, incorporating Teff as an additional variable globally improves the agreement between the abundance- and seismic-based ages for the APOKASC stars studied by Brewer et al. (2016). The calibrations presently available appear significantly less precise than what can be achieved for stars akin to the Sun (precision often below 1 Gyr), but some abundance ratios may hopefully prove more suitable than others outside the somewhat restricted domain defined by solar analogues in the future. A concern is that applying relationships constructed from a specific abundance dataset (such as that of Delgado Mena et al. 2017 used by DM19) is prone to the presence of study-to-study, zero-point abundance offsets, which are virtually unavoidable (e.g., Jofré et al. 2017). More generally, the possible existence of systematic differences between isochrone and seismic ages must be investigated further to improve the accuracy of the relationships (see Berger et al. 2020 for an example of comparison).
Constructing empirical abundance-age calibrations allowing one to accurately infer the age of stars spanning a wide range in spectral type, evolutionary status, and metallicity – and even born in different parts of the Galactic discs characterised by different star-formation histories – appears to be a daunting endeavour. A large number of stars with solar-like oscillations will be detected
by the TESS satellite, but this sample will be largely dominated by evolved sub-giants (Campante et al. 2016). The full accomplishment of this goal might thus be dependent on the release of seismic ages for thousands of bright main-sequence stars by the PLATO space missions (Rauer & Heras 2018), complemented by a high-quality determination of their surface chemical composition.
The LOS values of Lund et al. (2017) for these two stars are based on spectroscopic data obtained in early June 2014: see Kepler Community Follow-up Observing Program (CFOP); website https://exofop.ipac.caltech.edu/cfop.php.
See online tool at https://stilism.obspm.fr/.
We note that our spectroscopic observations of KIC 8006161 were accidentally obtained during an activity minimum (Karoff et al. 2018).
Available at http://kurucz.harvard.edu/linelists.html.
In the 30–110 range per resolution element at 5110 Å based on information about the observing runs in Furlan et al. (2018).
The adopted values are the preliminary Kepler Asteroseismology Science Consortium (KASC) estimates available at the time of the analysis (Silva Aguirre, priv. comm.) and are not taken from Pinsonneault et al. (2012, 2014), as erroneously quoted in C17 and SA17.
It should be noted that non-LTE corrections were applied to most elements by N17. Unfortunately, their LTE abundances cannot be computed because the mean differential corrections are not provided on a star-to-star basis. However, as discussed by N17, they are exceedingly small for solar analogues and probably irrelevant to the purpose of our discussion. Their uncertainties for elements other than iron also refer to [X/Fe], but we regard them in the following as being representative of those for [X/H].
For the very young (∼0.5 Gyr) star KIC 3733735, only [Y/Al] is very high, which points instead to a problem with the Al abundance.
Acknowledgments
We would like to thank the anonymous referee for her/his careful reading of the manuscript and thoughtful suggestions. TM acknowledges financial support from Belspo for contract PRODEX PLATO mission development, and wishes to thank Werner Verschueren (Belspo) for allowing the funding of the mission at OHP. JM, AM, and EW acknowledge support from the ERC Consolidator Grant funding scheme (project ASTEROCHRONOMETRY, https://www.asterochronometry.eu, G.A. n. 772293). This work has made use of data from the European Space Agency (ESA) mission Gaia (https://www.cosmos.esa.int/gaia), processed by the Gaia Data Processing and Analysis Consortium (DPAC, https://www.cosmos.esa.int/web/gaia/dpac/consortium). Funding for the DPAC has been provided by national institutions, in particular the institutions participating in the Gaia Multilateral Agreement. This article made use of AIMS, a software for fitting stellar pulsation data, developed in the context of the SPACEINN network, funded by the European Commission’s Seventh Framework Programme. This research has also made use of NASA’s Astrophysics Data System Bibliographic Services and the SIMBAD database operated at CDS, Strasbourg (France).
References
- Adelberger, E. G., García, A., Robertson, R. G. H., et al. 2011, Rev. Mod. Phys., 83, 195 [NASA ADS] [CrossRef] [Google Scholar]
- Adibekyan, V. Z., Santos, N. C., Sousa, S. G., & Israelian, G. 2011, A&A, 535, L11 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Adibekyan, V. Z., Sousa, S. G., Santos, N. C., et al. 2012, A&A, 545, A32 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Adibekyan, V., Delgado-Mena, E., Figueira, P., et al. 2016, A&A, 592, A87 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Amarsi, A. M., & Asplund, M. 2017, MNRAS, 464, 264 [NASA ADS] [CrossRef] [Google Scholar]
- Amarsi, A. M., Nissen, P. E., & Skúladóttir, Á. 2019a, A&A, 630, A104 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Amarsi, A. M., Barklem, P. S., Collet, R., Grevesse, N., & Asplund, M. 2019b, A&A, 624, A111 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Angulo, C., Arnould, M., Rayet, M., et al. 1999, Nucl. Phys., 656, 3 [Google Scholar]
- Asplund, M., Grevesse, N., Sauval, A. J., & Scott, P. 2009, ARA&A, 47, 481 [NASA ADS] [CrossRef] [Google Scholar]
- Ball, W. H., & Gizon, L. 2014, A&A, 568, A123 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Battistini, C., & Bensby, T. 2016, A&A, 586, A49 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Bazot, M., Benomar, O., Christensen-Dalsgaard, J., et al. 2019, A&A, 623, A125 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Beck, P. G., do Nascimento, J. D., Duarte, T., et al. 2017, A&A, 602, A63 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Bedell, M., Meléndez, J., Bean, J. L., et al. 2014, ApJ, 795, 23 [NASA ADS] [CrossRef] [Google Scholar]
- Bedell, M., Bean, J. L., Meléndez, J., et al. 2018, ApJ, 865, 68 [NASA ADS] [CrossRef] [Google Scholar]
- Bensby, T., Feltzing, S., & Oey, M. S. 2014, A&A, 562, A71 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Bergemann, M., Lind, K., Collet, R., Magic, Z., & Asplund, M. 2012, MNRAS, 427, 27 [NASA ADS] [CrossRef] [Google Scholar]
- Bergemann, M., Gallagher, A. J., Eitner, P., et al. 2019, A&A, 631, A80 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Berger, T. A., Huber, D., van Saders, J. L., et al. 2020, AJ, 159, 280 [Google Scholar]
- Biazzo, K., Gratton, R., Desidera, S., et al. 2015, A&A, 583, A135 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Böhm-Vitense, E. 1958, ZAp, 46, 108 [NASA ADS] [Google Scholar]
- Boyajian, T. S., von Braun, K., van Belle, G., et al. 2013, ApJ, 771, 40 [NASA ADS] [CrossRef] [Google Scholar]
- Bressan, A., Marigo, P., Girardi, L., et al. 2012, MNRAS, 427, 127 [NASA ADS] [CrossRef] [Google Scholar]
- Brewer, J. M., Fischer, D. A., Basu, S., Valenti, J. A., & Piskunov, N. 2015, ApJ, 805, 126 [NASA ADS] [CrossRef] [Google Scholar]
- Brewer, J. M., Fischer, D. A., Valenti, J. A., & Piskunov, N. 2016, ApJS, 225, 32 [Google Scholar]
- Bruntt, H., Bedding, T. R., Quirion, P.-O., et al. 2010, MNRAS, 405, 1907 [NASA ADS] [Google Scholar]
- Bruntt, H., Basu, S., Smalley, B., et al. 2012, MNRAS, 423, 122 [NASA ADS] [CrossRef] [Google Scholar]
- Buchhave, L. A., & Latham, D. W. 2015, ApJ, 808, 187 [NASA ADS] [CrossRef] [Google Scholar]
- Buder, S., Sharma, S., Kos, J., et al. 2020, MNRAS, submitted [arXiv:2011.02505] [Google Scholar]
- Buldgen, G., Reese, D. R., & Dupret, M. A. 2016, A&A, 585, A109 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Campante, T. L., Schofield, M., Kuszlewicz, J. S., et al. 2016, ApJ, 830, 138 [NASA ADS] [CrossRef] [Google Scholar]
- Carrier, F., Eggenberger, P., D’Alessandro, A., & Weber, L. 2005, New A, 10, 315 [Google Scholar]
- Casagrande, L., Portinari, L., Glass, I. S., et al. 2014, MNRAS, 439, 2060 [Google Scholar]
- Casali, G., Spina, L., Magrini, L., et al. 2020, A&A, 639, A127 [CrossRef] [EDP Sciences] [Google Scholar]
- Castelli, F., & Kurucz, R. L. 2003, in Modelling of Stellar Atmospheres, eds. N. Piskunov, W. W. Weiss, & D. F. Gray, 210, A20 [Google Scholar]
- Chaplin, W. J., & Miglio, A. 2013, ARA&A, 51, 353 [Google Scholar]
- Chaplin, W. J., Basu, S., Huber, D., et al. 2014, ApJS, 210, 1 [Google Scholar]
- Christensen-Dalsgaard, J. 2008a, Ap&SS, 316, 113 [Google Scholar]
- Christensen-Dalsgaard, J. 2008b, Ap&SS, 316, 13 [NASA ADS] [CrossRef] [Google Scholar]
- Cochran, W. D., Hatzes, A. P., Butler, R. P., & Marcy, G. W. 1997, ApJ, 483, 457 [Google Scholar]
- Compton, D. L., Bedding, T. R., Ball, W. H., et al. 2018, MNRAS, 479, 4416 [Google Scholar]
- Cox, J. P., & Giuli, R. T. 1968, Principles of Stellar Structure [Google Scholar]
- Creevey, O. L., Doǧan, G., Frasca, A., et al. 2012, A&A, 537, A111 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Creevey, O. L., Metcalfe, T. S., Schultheis, M., et al. 2017, A&A, 601, A67 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- da Silva, R., Porto de Mello, G. F., Milone, A. C., et al. 2012, A&A, 542, A84 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Deal, M., Richard, O., & Vauclair, S. 2015, A&A, 584, A105 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Delgado Mena, E., Tsantaki, M., Adibekyan, V. Z., et al. 2017, A&A, 606, A94 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Delgado Mena, E., Moya, A., Adibekyan, V., et al. 2019, A&A, 624, A78 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Dotter, A., Conroy, C., Cargile, P., & Asplund, M. 2017, ApJ, 840, 99 [Google Scholar]
- Doyle, A. P., Davies, G. R., Smalley, B., Chaplin, W. J., & Elsworth, Y. 2014, MNRAS, 444, 3592 [Google Scholar]
- Doyle, A. P., Smalley, B., Faedi, F., Pollacco, D., & Gómez Maqueo Chew, Y. 2017, MNRAS, 469, 4850 [Google Scholar]
- Duquennoy, A., & Mayor, M. 1991, A&A, 248, 485 [NASA ADS] [Google Scholar]
- Edvardsson, B., Andersen, J., Gustafsson, B., et al. 1993, A&A, 275, 101 [NASA ADS] [Google Scholar]
- El-Badry, K., Ting, Y.-S., Rix, H.-W., et al. 2018a, MNRAS, 476, 528 [Google Scholar]
- El-Badry, K., Rix, H.-W., Ting, Y.-S., et al. 2018b, MNRAS, 473, 5043 [Google Scholar]
- Farnir, M., Dupret, M. A., Buldgen, G., et al. 2020, A&A, 644, A37 [EDP Sciences] [Google Scholar]
- Feltzing, S., Howes, L. M., McMillan, P. J., & Stonkutė, E. 2017, MNRAS, 465, L109 [Google Scholar]
- Ferguson, J. W., Alexander, D. R., Allard, F., et al. 2005, ApJ, 623, 585 [Google Scholar]
- Flores, M., González, J. F., Jaque Arancibia, M., Buccino, A., & Saffe, C. 2016, A&A, 589, A135 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Foreman-Mackey, D., Hogg, D. W., Lang, D., & Goodman, J. 2013, PASP, 125, 306 [Google Scholar]
- Formicola, A., Imbriani, G., Costantini, H., et al. 2004, Phys. Lett. B, 591, 61 [Google Scholar]
- Furlan, E., Ciardi, D. R., Cochran, W. D., et al. 2018, ApJ, 861, 149 [Google Scholar]
- Gabriel, M., Noels, A., Montalbán, J., & Miglio, A. 2014, A&A, 569, A63 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Gaia Collaboration (Brown, A. G. A., et al.) 2018, A&A, 616, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Gallagher, A. J., Bergemann, M., Collet, R., et al. 2020, A&A, 634, A55 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- García, R. A., Ceillier, T., Salabert, D., et al. 2014, A&A, 572, A34 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- García Pérez, A. E., Allende Prieto, C., Holtzman, J. A., et al. 2016, AJ, 151, 144 [Google Scholar]
- Ge, Z. S., Bi, S. L., Li, T. D., et al. 2015, MNRAS, 447, 680 [Google Scholar]
- Gilmore, G., Randich, S., Asplund, M., et al. 2012, The Messenger, 147, 25 [NASA ADS] [Google Scholar]
- Gray, D. F., Tycner, C., & Brown, K. 2000, PASP, 112, 328 [Google Scholar]
- Grevesse, N., & Sauval, A. J. 1998, Space Sci. Rev., 85, 161 [Google Scholar]
- Gustafsson, B., Edvardsson, B., Eriksson, K., et al. 2008, A&A, 486, 951 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Halbwachs, J. L., Mayor, M., & Udry, S. 2018, A&A, 619, A81 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Heiter, U., Jofré, P., Gustafsson, B., et al. 2015, A&A, 582, A49 [NASA ADS] [CrossRef] [EDP Sciences] [PubMed] [Google Scholar]
- Hidalgo, S. L., Pietrinferni, A., Cassisi, S., et al. 2018, ApJ, 856, 125 [Google Scholar]
- Hinkle, K., Wallace, L., Valenti, J., & Harmer, D. 2000, Visible and Near Infrared Atlas of the Arcturus Spectrum 3727–9300 A [Google Scholar]
- Huber, D., Ireland, M. J., Bedding, T. R., et al. 2012, ApJ, 760, 32 [Google Scholar]
- Iglesias, C. A., & Rogers, F. J. 1996, ApJ, 464, 943 [NASA ADS] [CrossRef] [Google Scholar]
- Irwin, A. W. 2012, FreeEOS: Equation of State for Stellar Interiors Calculations [Google Scholar]
- Jofré, P., Heiter, U., Soubiran, C., et al. 2014, A&A, 564, A133 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Jofré, P., Heiter, U., Worley, C. C., et al. 2017, A&A, 601, A38 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Jofré, P., Jackson, H., & Tucci Maia, M. 2020, A&A, 633, L9 [CrossRef] [EDP Sciences] [Google Scholar]
- Karoff, C., Metcalfe, T. S., Chaplin, W. J., et al. 2013, MNRAS, 433, 3227 [Google Scholar]
- Karoff, C., Metcalfe, T. S., Santos, Â. R. G., et al. 2018, ApJ, 852, 46 [Google Scholar]
- Karovicova, I., White, T. R., Nordlander, T., et al. 2018, MNRAS, 475, L81 [Google Scholar]
- Kervella, P., Bigot, L., Gallenne, A., & Thévenin, F. 2017, A&A, 597, A137 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Kiefer, R., Schad, A., Davies, G., & Roth, M. 2017, A&A, 598, A77 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Komatsu, E., Smith, K. M., Dunkley, J., et al. 2011, ApJS, 192, 18 [NASA ADS] [CrossRef] [Google Scholar]
- Krishna Swamy, K. S. 1966, ApJ, 146, 731 [Google Scholar]
- Lallement, R., Capitanio, L., Ruiz-Dern, L., et al. 2018, A&A, 616, A132 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Laws, C., & Gonzalez, G. 2001, ApJ, 553, 405 [Google Scholar]
- Lebreton, Y., & Goupil, M. J. 2014, A&A, 569, A21 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Li, T., Bedding, T. R., Christensen-Dalsgaard, J., et al. 2020, MNRAS, 495, 3431 [Google Scholar]
- Lin, J., Asplund, M., Ting, Y.-S., et al. 2020, MNRAS, 491, 2043 [Google Scholar]
- Lodders, K. 2003, ApJ, 591, 1220 [Google Scholar]
- Lund, M. N., & Reese, D. R. 2018, in Asteroseismology and Exoplanets: Listening to the Stars and Searching for New Worlds, eds. T. L. Campante, N. C. Santos, & M. J. P. F. G. Monteiro, 49, 149 [Google Scholar]
- Lund, M. N., Silva Aguirre, V., Davies, G. R., et al. 2017, ApJ, 835, 172 [Google Scholar]
- Mamajek, E. E., Prsa, A., Torres, G., et al. 2015, ArXiv e-prints [arXiv:1510.07674] [Google Scholar]
- Martinez, C. F., Cunha, K., Ghezzi, L., & Smith, V. V. 2019, ApJ, 875, 29 [Google Scholar]
- Mashonkina, L., Gehren, T., & Bikmaev, I. 1999, A&A, 343, 519 [NASA ADS] [Google Scholar]
- Meléndez, J., Ramírez, I., Karakas, A. I., et al. 2014, ApJ, 791, 14 [Google Scholar]
- Metcalfe, T. S., Creevey, O. L., & Christensen-Dalsgaard, J. 2009, ApJ, 699, 373 [Google Scholar]
- Metcalfe, T. S., Monteiro, M. J. P. F. G., Thompson, M. J., et al. 2010, ApJ, 723, 1583 [Google Scholar]
- Miglio, A., Chiappini, C., Mackereth, T., et al. 2021, A&A, 645, A85 [EDP Sciences] [Google Scholar]
- Molenda-Żakowicz, J., Sousa, S. G., Frasca, A., et al. 2013, MNRAS, 434, 1422 [Google Scholar]
- Montalbán, J., Mackereth, J. T., Miglio, A., et al. 2020, ArXiv e-prints [arXiv:2006.01783] [Google Scholar]
- Morel, T. 2018, A&A, 615, A172 [CrossRef] [EDP Sciences] [Google Scholar]
- Morel, T., Micela, G., Favata, F., Katz, D., & Pillitteri, I. 2003, A&A, 412, 495 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Morel, T., Miglio, A., Lagarde, N., et al. 2014, A&A, 564, A119 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Mott, A., Steffen, M., Caffau, E., & Strassmeier, K. G. 2020, A&A, 638, A58 [EDP Sciences] [Google Scholar]
- Neves, V., Santos, N. C., Sousa, S. G., Correia, A. C. M., & Israelian, G. 2009, A&A, 497, 563 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Nissen, P. E. 2015, A&A, 579, A52 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Nissen, P. E. 2016, A&A, 593, A65 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Nissen, P. E., & Gustafsson, B. 2018, A&ARv., 26, 6 [Google Scholar]
- Nissen, P. E., Silva Aguirre, V., Christensen-Dalsgaard, J., et al. 2017, A&A, 608, A112 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Nissen, P. E., Christensen-Dalsgaard, J., Mosumgaard, J. R., et al. 2020, A&A, 640, A81 [CrossRef] [EDP Sciences] [Google Scholar]
- Nordlander, T., & Lind, K. 2017, A&A, 607, A75 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- North, J. R., Davis, J., Robertson, J. G., et al. 2009, MNRAS, 393, 245 [Google Scholar]
- Pietrinferni, A., Cassisi, S., Salaris, M., & Castelli, F. 2004, ApJ, 612, 168 [NASA ADS] [CrossRef] [Google Scholar]
- Pinsonneault, M. H., An, D., Molenda-Żakowicz, J., et al. 2012, ApJS, 199, 30 [Google Scholar]
- Pinsonneault, M. H., Elsworth, Y., Epstein, C., et al. 2014, ApJS, 215, 19 [Google Scholar]
- Price-Whelan, A. M., Hogg, D. W., Rix, H.-W., et al. 2018, AJ, 156, 18 [Google Scholar]
- Prochaska, J. X., Naumov, S. O., Carney, B. W., McWilliam, A., & Wolfe, A. M. 2000, AJ, 120, 2513 [NASA ADS] [CrossRef] [Google Scholar]
- Ramírez, I., Meléndez, J., Cornejo, D., Roederer, I. U., & Fish, J. R. 2011, ApJ, 740, 76 [Google Scholar]
- Ramírez, I., Khanal, S., Aleo, P., et al. 2015, ApJ, 808, 13 [Google Scholar]
- Rauer, H., & Heras, A. M. 2018, Space Missions for Exoplanet Science: PLATO, 86 [Google Scholar]
- Reddy, B. E., Tomkin, J., Lambert, D. L., & Allende Prieto, C. 2003, MNRAS, 340, 304 [Google Scholar]
- Reddy, B. E., Lambert, D. L., & Allende Prieto, C. 2006, MNRAS, 367, 1329 [Google Scholar]
- Reese, D. R., Chaplin, W. J., Davies, G. R., et al. 2016, A&A, 592, A14 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Rendle, B. M., Buldgen, G., Miglio, A., et al. 2019, MNRAS, 484, 771 [Google Scholar]
- Rogers, F. J., & Nayfonov, A. 2002, ApJ, 576, 1064 [Google Scholar]
- Rosman, K. J. R., & Taylor, P. D. P. 1998, J. Phys. Chem. Ref. Data, 27, 1275 [Google Scholar]
- Sahlholdt, C. L., Feltzing, S., Lindegren, L., & Church, R. P. 2019, MNRAS, 482, 895 [Google Scholar]
- Sarmento, P., Delgado Mena, E., Rojas-Ayala, B., & Blanco-Cuaresma, S. 2020, A&A, 636, A85 [EDP Sciences] [Google Scholar]
- Schonhut-Stasik, J. S., Baranec, C., Huber, D., et al. 2017, ApJ, 847, 97 [Google Scholar]
- Schuler, S. C., Cunha, K., Smith, V. V., et al. 2011, ApJ, 737, L32 [Google Scholar]
- Scuflaire, R., Théado, S., Montalbán, J., et al. 2008a, Ap&SS, 316, 83 [Google Scholar]
- Scuflaire, R., Montalbán, J., Théado, S., et al. 2008b, Ap&SS, 316, 149 [Google Scholar]
- Sellwood, J. A., & Binney, J. J. 2002, MNRAS, 336, 785 [Google Scholar]
- Serenelli, A., Johnson, J., Huber, D., et al. 2017, ApJS, 233, 23 [Google Scholar]
- Silva Aguirre, V., Lund, M. N., Antia, H. M., et al. 2017, ApJ, 835, 173 [Google Scholar]
- Silva Aguirre, V., Bojsen-Hansen, M., Slumstrup, D., et al. 2018, MNRAS, 475, 5487 [NASA ADS] [Google Scholar]
- Skúladóttir, Á., Hansen, C. J., Salvadori, S., & Choplin, A. 2019, A&A, 631, A171 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Sneden, C. A. 1973, Ph.D. Thesis, The University of Texas AT Austin, USA [Google Scholar]
- Spada, F., Demarque, P., Kim, Y.-C., & Sills, A. 2013, ApJ, 776, 87 [Google Scholar]
- Spina, L., Meléndez, J., Karakas, A. I., et al. 2016, A&A, 593, A125 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Spina, L., Meléndez, J., Karakas, A. I., et al. 2018, MNRAS, 474, 2580 [NASA ADS] [Google Scholar]
- Spina, L., Nordlander, T., Casey, A. R., et al. 2020, ApJ, 895, 52 [Google Scholar]
- Takeda, Y. 2005, PASJ, 57, 83 [NASA ADS] [Google Scholar]
- Thoul, A. A., Bahcall, J. N., & Loeb, A. 1994, ApJ, 421, 828 [NASA ADS] [CrossRef] [Google Scholar]
- Titarenko, A., Recio-Blanco, A., de Laverny, P., Hayden, M., & Guiglion, G. 2019, A&A, 622, A59 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Torres, G., Fischer, D. A., Sozzetti, A., et al. 2012, ApJ, 757, 161 [Google Scholar]
- Tucci Maia, M., Ramírez, I., Meléndez, J., et al. 2016, A&A, 590, A32 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Tucci Maia, M., Meléndez, J., Lorenzo-Oliveira, D., Spina, L., & Jofré, P. 2019, A&A, 628, A126 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Valle, G., Dell’Omodarme, M., Prada Moroni, P. G., & Degl’Innocenti, S. 2018, A&A, 620, A168 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Verma, K., Faria, J. P., Antia, H. M., et al. 2014, ApJ, 790, 138 [Google Scholar]
- Verma, K., Raodeo, K., Basu, S., et al. 2019, MNRAS, 483, 4678 [Google Scholar]
- White, T. R., Huber, D., Maestro, V., et al. 2013, MNRAS, 433, 1262 [Google Scholar]
- Yana Galarza, J., Meléndez, J., Lorenzo-Oliveira, D., et al. 2019, MNRAS, 490, L86 [Google Scholar]
Appendix A: Seismic analysis of KIC 9965715
A.1. First approach (ASTEC+ADIPLS)
We approached the new interpretation of the data from KIC 9965715 in two steps. The first step consisted of using a dense grid of BASTI stellar models (Pietrinferni et al. 2004; Hidalgo et al. 2018) to match the Teff, ⟨Δν⟩, and νmax. The global seismic quantities were taken from Lund et al. (2017). These models use a chemical enrichment law, and we refer to the relevant papers for details of the physics. We used a simple Bayesian approach in 2D (mass, age) on grids of different metallicities with mild priors that have little influence on our results: mass is restricted to [0.8,1.3] M⊙, age is restricted to [0,10] Gyr, and we used a prior on the evolution state defined as the relative amount of time spent in the specific phase of the star’s life (main-sequence, sub-giant, and giant). The uncertainties are obtained by marginalising over the mass and age, and they are defined as the 68% confidence interval centred on the median values. We also recorded the model observables, and so we can define an optimal surface gravity log g = 4.280 ± 0.004, where we also accounted for the different metallicity value.
In a second step, we used the fitted parameters from the first step as initial starting points for a detailed stellar modelling. The use of the global seismic quantities brings important constraints on certain stellar parameters such as the surface gravity or density. However, for the best precision and improved accuracy (scaling relations are still subject to some systematic errors) on other parameters such as mass and age, we must use the information contained in the detailed seismic data. We adopted a similar approach to that presented in C17 by using the frequency ratios (r01, r02) in the optimisation. These values were calculated directly from the individual frequencies from Lund et al. (2017), and we derived the covariance matrix C by performing simple Monte Carlo-like simulations.
In the optimisation approach, we calculated the following likelihood:
where
Here, x are the observational constraints and xM are the corresponding observables from the model. We additionally included the constraint of the mean large and small frequency separations, ⟨Δν⟩ and ⟨δν⟩, by calculating this value within the range of the observed frequencies. We also included the metallicity and Teff constraint.
For this second step of stellar modelling, we used the Aarhus STellar Evolution Code (ASTEC) and the Adiabatic Pulsation code (ADIPLS) from Christensen-Dalsgaard (2008a,b). The setup of the physics is as follows: we assumed a non-rotating, non-magnetic star. The opacities and equation-of-state tables were taken from the OPAL collaboration (Iglesias & Rogers 1996; Rogers & Nayfonov 2002). Nuclear reaction rates were taken from Angulo et al. (1999), and we included the values obtained by the LUNA collaboration for the 14N(p,γ)15O reaction (Formicola et al. 2004). We used the solar mixture of Grevesse & Sauval (1998). Apart from the fixed physics described above, the only other parameters that control the evolution of the star are the mass, M, and the initial metallicity and helium mass fraction, Zi and Yi, where Xi + Yi + Zi = 1, and X refers to the hydrogen mass fraction. We also needed to set the mixing-length parameter α, defined as the ratio of the mean-free path of a fluid element, and the pressure scale height (see Böhm-Vitense 1958). The final parameter that controls the current structure of the model is the age.
We approached our optimisation of the parameters by building several small 3D grids (M, Zi, and age) while fixing (Yi, α) at discrete values. We show some of these results in Table A.1, where the first set of (Yi, α) is equivalent to the values used in the BASTI library. The parameters are defined as the maximum likelihood ones, and the uncertainties were calculated as the 68% confidence interval with the highest probability, as shown by the dashed lines in Fig. A.1. We find that the mass varies between 1.05 and 1.13 M⊙, with a strong dependence on Yi. The optimal age remains in the range of 2.7–3.3 Gyr. We adopt the values from the third column, and, to account for the other parameters, we add the differences to the uncertainties in quadrature to obtain M = 1.07 ± 0.04 M⊙ and an age of 2.8 ± 0.4 Gyr.
![]() |
Fig. A.1. Marginalised distributions of mass and age for the set (Yi, α) = (0.287,1.94), scaled so the maximum is equal to unity. The adopted stellar parameter (maximum likelihood) and the lower and upper values (68% confidence interval) are shown by the dashed lines. |
Stellar parameters for KIC 9965715 using the first approach.
A.2. Second approach (CLES+LOSC)
The stellar parameters of KIC 9965715 were also estimated by using the open-source code AIMS (Asteroseismic Inference on a Massive Scale; Reese et al. 2016; Lund et al. 2018; Rendle et al. 2019) that implements a Bayesian inference approach. AIMS evaluates the posterior distributions of stellar parameters using a Markov chain Monte Carlo (MCMC) ensemble sampler (Foreman-Mackey et al. 2013), and it selects stellar models that best fit observational data by interpolating (evolutionary tracks and frequencies) in a pre-computed grid.
The grid of stellar models at the base of the procedure was computed using the code CLES (Scuflaire et al. 2008a; Gabriel et al. 2014), and the oscillation frequencies of radial (ℓ = 0) and non-radial (ℓ = 1, 2) modes for each stellar model using the adiabatic oscillation code LOSC (Scuflaire et al. 2008b). The stellar model’s computation follows the evolution from the pre-main sequence up to the point where the hydrogen content at the stellar centre, Xc, is below 10−6. The prescriptions for the input physics are the following: nuclear reaction rates of Adelberger et al. (2011) or from NACRE (Angulo et al. 1999) if those of the former are not available. The equation of state is FreeEoS (Irwin 2012). We adopted the solar metal mixture from Asplund et al. (2009) and opacity tables were built combining OPAL (Iglesias & Rogers 1996) opacity values at high temperature with those of Wichita State University for the low temperature domain (Ferguson et al. 2005). Atmospheric boundary conditions were provided by the Krishna-Swamy T(τ) law (Krishna Swamy 1966), and the corresponding atmosphere structure was added on the top of the interior model. Convection was treated with the ‘mixing length’ formalism (Cox & Giuli 1968). The corresponding α parameter was kept fixed for all the grids, and it was derived from the solar calibration (α = 2.338) with the same physics and including microscopic diffusion of chemical elements (Thoul et al. 1994).
We assumed extra mixing at the boundaries of convective regions introducing core-overshooting (ov = 0.1 and 0.2 Hp – pressure scale height), and undershooting below the convective envelope (∼0.2 Hp). In both cases, the chemical mixing is instantaneous and the temperature gradient in the mixed region is the radiative one. Finally, microscopic diffusion of chemical elements is included in one of the grids.
The grids are parametrised by mass (from 0.9 to 1.8 M⊙ with a step of 0.01 M⊙) and initial [Fe/H] (from −0.45 to 0.00, with a step of 0.05 dex). We assumed an enrichment law Δ Y/Δ Z provided by the solar calibration and a primordial He abundance, Yp, of 0.2485 (Komatsu et al. 2011).
The stellar parameters of KIC 9965715 were derived using these grids and the following observational constraints: Teff and [Fe/H], as classical constraints, as well as the average νmax and the individual frequencies of 40 modes (ℓ = 0, 1, and 2) of Lund et al. (2017) as seismic data. We also performed a second set of optimisations using the values of frequency separation ratios (r01, r10, r02) and the frequency of the lowest order radial mode as seismic information. In the former case, we used the two-terms surface effect correction from Ball & Gizon (2014).
The analysis of the results allowed us to rule out models with an extra mixing as large as 0.2 Hp. On the other hand, the overshooting below the convective envelope has no effect on the stellar parameters derived. From ‘Bayesian evidence’ values, models including microscopic diffusion of helium and metals are clearly preferred over those without it.
The results of the different optimisations we performed indicate an initial chemical composition close to the solar one. So, stellar parameters based on non-diffusion grids tend to deviate from the observed surface chemical composition. A good match of surface composition can be obtained at the cost of a mediocre fit of the seismic properties. The stellar parameters based on frequency separation ratios are in good agreement with those derived using the individual frequencies. The parameters for each optimisation are defined from the posterior distributions, as the maximum likelihood ones, and the uncertainties are calculated as the 68% confidence interval with the highest probability. We summarise these results in Table A.2 in two sets of stellar parameters, depending on whether they are based on the model grids with or without microscopic diffusion. We adopt the values of the first column and inflate the uncertainties to take into account that the target could be up to 1 Gyr older if the constraint on the surface chemical composition is relaxed.
Stellar parameters for KIC 9965715 using the second approach.
Appendix B: Abundance results
Abundance results for the stars in our sample.
Appendix C: Comparison with literature results for 16 Cyg AB
The detailed abundance pattern of the components in the well-studied binary system 16 Cyg AB has recently been investigated with high precision by a number of studies. These are the two brightest solar analogues in the Kepler field. Noteworthy is the fact that the primary is slightly more metal rich than the secondary. The following results were obtained as a function of the increasing significance level: Δ[Fe/H] (A–B) = +0.025 ± 0.009 (Laws & Gonzalez 2001), +0.031 ± 0.010 (N17), and +0.041 ± 0.004 dex (Tucci Maia et al. 2019). The origin of this difference in metal content is unclear, but it might be related to planetary formation and possibly subsequent engulfment. A Jupiter-mass planet is actually known to orbit 16 Cyg B along a very eccentric (e ∼ 0.6) orbit (Cochran et al. 1997).
In addition to the analysis described in Sect. 5, we performed a differential analysis of 16 Cyg A with respect to 16 Cyg B. The parameters of 16 Cyg B in Table 1 were adopted. Our analysis procedures strictly followed those employed relative to solar for the other targets. Such an approach for binary components with very similar parameters is expected to increase the precision of the results, and, consequently, to better reveal subtle chemical differences (e.g., Ramírez et al. 2015). We indeed find that it leads to a better precision with respect to the solar analysis. For instance, the line-to-line scatter in the iron abundances is divided by a factor of about two (here a mere ∼0.01 dex). The results are given in Table C.1.
With respect to N17 and Tucci Maia et al. (2019), we obtain a slightly cooler Teff scale by ∼20 K on average. Figure C.1 shows a comparison between the abundance patterns of the two components. As can be seen, there is an overall remarkable agreement with the literature data, with average differences at the 0.01-dex level. On an element-to-element basis, the differences are systematically below 0.05 dex. Considering the weighted (by the inverse variance) average of all the abundances with respect to hydrogen, [X/H], we obtain ⟨Δ[X/H]⟩ (A–B) = +0.0209 ± 0.0031 dex. This is compelling on statistical grounds (at a 6.7σ significance level) and supports the slight enhancement in metals of the primary. This conclusion is strengthened by the fact that Δ[X/H] (A–B) is systematically positive for the 20 elements (Fig. C.1). The stellar parameters of the two components are so close that, for example, differential diffusion effects are not expected to modify this conclusion (Deal et al. 2015). The detection of such a small metallicity difference between the two components indicates that a precision at the 0.02–0.03 dex level is realistically achieved for the abundances of the solar analogues in our sample (once again, it does not pertain to stars outside of this category). It is an interesting conclusion in the context of large-scale spectroscopic surveys, as it shows that such subtle abundance differences can be revealed even with data of limited quality – especially resolving power – provided that adequate analysis strategies are implemented. However, we are unable to find evidence for the increase of Δ[X/H] (A–B) as a function of Tc reported by Tucci Maia et al. (2019), and, to a lesser extent, by Laws & Gonzalez (2001) and N17. In contrast, the Δ[X/H] (A–B)-Tc behaviour we find is completely flat (–0.21 ± 0.79 × 10−5 dex K−1). This might be ascribed to our lower data quality (but see, e.g., Ramírez et al. 2011).
![]() |
Fig. C.1. Abundance patterns with respect to hydrogen, [X/H], and differences with respect to literature values (this study minus literature), δ[X/H], as a function of Tc. The left and right panels show the comparison with N17 and Tucci Maia et al. (2019), respectively. We ignored the molecular-based C abundances of Tucci Maia et al. (2019). The results for 16 Cyg A, 16 Cyg B, and 16 Cyg A with respect to 16 Cyg B are shown in the top, middle, and bottom panels, respectively. Our [X/H] results and those in the literature are shown as filled and open circles, respectively. A dotted, horizontal line is drawn at our [Fe/H] value. The solid lines show the weighted, linear fit of [X/H] as a function of Tc. The fits obtained in the literature for the full set of elements are overplotted as dashed lines (the trends for 16 Cyg AB with respect to the Sun are not available in Tucci Maia et al. 2019). The average δ[X/H] values are given in the relevant panels (the number of elements in common is indicated in brackets). To guide the eye, a dotted line is drawn at δ[X/H] = 0. |
All Tables
Unweighted mean age deviation with respect to seismic estimates (abundance-based minus seismic) for each abundance indicator and relation (the number in brackets is the number of stars the calculation is based on).
All Figures
![]() |
Fig. 1. Comparison for Kepler LEGACY sample between the stellar mass (top panel), age (middle panel), and surface gravity (bottom panel) determined by C17 and SA17. The differences for a given parameter are the values of C17 minus those of SA17. The stars in our sample are shown as filled red circles. KIC 9965715 is shown for completeness as a cross. |
In the text |
![]() |
Fig. 2. Comparison for Si (top panel) and Cr (bottom panel) between the mean abundances yielded by the neutral and singly-ionised ion. The mean difference is given in each panel (the number in brackets is the number of stars the calculation is based on). A representative error bar is also shown. |
In the text |
![]() |
Fig. 3. Upper panels: differences (this study minus literature) between our Teff and [Fe/H] values and those adopted by C17 and SA17 for their seismic modelling. Bottom panels: same as upper panels, but for the SDSS and ASPCAP non-seismic constraints adopted by Serenelli et al. (2017). |
In the text |
![]() |
Fig. 4. Comparison between our ionisation temperatures (constrained and unconstrained shown as filled and open black circles, respectively) and those in the literature derived from interferometric and absolute flux measurements (filled red squares). The blue crosses show the results of N17 determined from iron ionisation balance. The stars are ordered as a function of increasing Teff. |
In the text |
![]() |
Fig. 5. Upper panel: abundance pattern of KIC 6106415 with respect to hydrogen, [X/H], as a function of Tc. Our [X/H] results and those of N17 are shown as filled and open circles, respectively. A dotted, horizontal line is drawn at our [Fe/H] value. The solid line shows our weighted, linear fit of [X/H] as a function of Tc. The fit obtained by N17 is overplotted as a dashed line. Lower panel: differences (our values minus theirs) with respect to N17, δ[X/H]. The average δ[X/H] value is given (the number of elements in common is indicated in brackets). To guide the eye, a dotted line is drawn at δ[X/H] = 0. |
In the text |
![]() |
Fig. 6. Comparison between our abundance ratios and those of Bruntt et al. (2012) and Brewer et al. (2016). The abundance differences are those of this study minus those of the literature. The vertical dashed line connects the extreme values. The box covers the first to third quartile of the data, while the thick horizontal line inside the box shows the median. As recommended by Bruntt et al. (2012), we chose their Fe II-based abundances for [Fe/H] (but used the Mg I, Si I, Ti I, and Cr I abundances). For KIC 5184732, KIC 7970740, and KIC 8006161, we chose the values of Brewer et al. (2016) corresponding to the spectrum with the highest S/N. The number of stars the calculation is based on is indicated for each abundance ratio. The red symbols show the abundance ratios used as age indicators (Sect. 9). |
In the text |
![]() |
Fig. 7. Abundances of neutron-capture elements relative to iron as a function of seismic age. The thin- and thick-disc stars are shown as circles and squares, respectively. The solar age is shown as a vertical dotted line. |
In the text |
![]() |
Fig. 8. Left panels: variation of abundance indicators, as a function of seismic age. The data are colour-coded as a function of Teff. The thin- and thick-disc stars are shown as circles and squares, respectively. The solar age is shown as a vertical dotted line. Middle and right panels: same as left panels, but colour-coded as a function of [Fe/H] and M, respectively. Linear 1D relations in the literature linking abundance ratios and ages are overplotted as dashed lines: Bedell et al. (2018, red), DM19 (blue), Nissen (2016) and Nissen et al. (2020, cyan), Spina et al. (2018, magenta), Titarenko et al. (2019, yellow), and Tucci Maia et al. (2016, green). |
In the text |
![]() |
Fig. 9. Mean ages of KIC 6106415 and 16 Cyg AB obtained using our abundance data (circles), those of N17 (squares), and those of Tucci Maia et al. (2019, triangles). The shaded horizontal stripes show the seismic ages, along with their ±1σ uncertainty. |
In the text |
![]() |
Fig. 10. Mean age deviation, ΔA, with respect to the seismic estimate (abundance-based minus seismic) when considering all the abundance indicators, as a function of Teff, [Fe/H], and M. The quantity δA is a measure of the improvement brought about by using 2D or 3D instead of 1D relations (see Sect. 9.1). The thin- and thick-disc stars are shown as filled and open circles, respectively. The Sun is shown as an open triangle. The rightmost panels show the breakdown of the ΔA and δA values. |
In the text |
![]() |
Fig. 11. Same as Fig. 10, but only for the four abundance indicators common to all the relations. |
In the text |
![]() |
Fig. 12. Differences between stellar parameters of Brewer et al. (2016) and Serenelli et al. (2017). The mean values for the SDSS and ASPCAP Teff scales are given in each panel. We assumed the ASPCAP [M/H] values that are the overall scaled solar abundances. |
In the text |
![]() |
Fig. 13. Abundance-age trends for our data (filled symbols) and those of Brewer et al. (2016, open symbols). The ages were derived by different means: detailed modelling of the oscillation frequencies and from global seismic quantities, respectively. The very high [Y/Al] abundance ratio (∼+1.35) obtained by Brewer et al. (2016) for KIC 3733735 is off scale. The thin- and thick-disc stars are shown as circles and squares, respectively. Linear 1D relations in the literature are overplotted as dashed lines: Bedell et al. (2018, red), DM19 (blue), Nissen (2016) and Nissen et al. (2020, cyan), Spina et al. (2018, magenta), Titarenko et al. (2019, yellow), and Tucci Maia et al. (2016, green). The solar age is shown as a vertical dotted line. |
In the text |
![]() |
Fig. 14. Same as Fig. 10, but for the extended sample. |
In the text |
![]() |
Fig. A.1. Marginalised distributions of mass and age for the set (Yi, α) = (0.287,1.94), scaled so the maximum is equal to unity. The adopted stellar parameter (maximum likelihood) and the lower and upper values (68% confidence interval) are shown by the dashed lines. |
In the text |
![]() |
Fig. C.1. Abundance patterns with respect to hydrogen, [X/H], and differences with respect to literature values (this study minus literature), δ[X/H], as a function of Tc. The left and right panels show the comparison with N17 and Tucci Maia et al. (2019), respectively. We ignored the molecular-based C abundances of Tucci Maia et al. (2019). The results for 16 Cyg A, 16 Cyg B, and 16 Cyg A with respect to 16 Cyg B are shown in the top, middle, and bottom panels, respectively. Our [X/H] results and those in the literature are shown as filled and open circles, respectively. A dotted, horizontal line is drawn at our [Fe/H] value. The solid lines show the weighted, linear fit of [X/H] as a function of Tc. The fits obtained in the literature for the full set of elements are overplotted as dashed lines (the trends for 16 Cyg AB with respect to the Sun are not available in Tucci Maia et al. 2019). The average δ[X/H] values are given in the relevant panels (the number of elements in common is indicated in brackets). To guide the eye, a dotted line is drawn at δ[X/H] = 0. |
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.