1 Draft version May 5, 2021 Typeset using LATEX twocolumn style in AASTeX63

The consistency of chemical clocks among coeval

Francisca Espinoza-Rojas,1 Julio Chaname,´ 1 Paula Jofre,´ 2 and Laia Casamiquela3

1Instituto de Astrof´ısica, Pontificia Universidad Cat´olica de Chile Av. Vicu˜naMackenna 4860, 782-0436 Macul, Santiago, Chile 2N´ucleo de Astronom´ıa,Facultad de Ingenier´ıay Ciencias, Universidad Diego Portales Ej´ercito 441, Santiago, Chile 3Laboratoire d’Astrophysique de Bordeaux, Univ. Bordeaux CNRS, B18N, all´eeGeoffroy Saint-Hilaire, 33615 Pessac, France

(Received XX; Revised YY; Accepted ZZ)

ABSTRACT The abundance ratios of some chemical species have been found to correlate with stellar age, leading to the possibility of using stellar atmospheric abundances as stellar age indicators. These chemical clocks have been calibrated with solar-twins, open clusters and red giants, but it remains to be seen whether they can be effective at identifying coeval stars in a field population that spans a broad parameter space (i.e., the promise of chemical tagging). Since the components of wide binaries are known to be stars of common origins, they constitute ideal laboratories for testing the usefulness of chemical clocks for the age dating of field stars. We determined the abundances of a new sample of 5 binaries and collected data for other 31 systems from the literature in order to test the applicability of chemical clocks. We recover the well known result that the components of wide binaries have more consistent chemistry than that of random pairs. However, we also show for the first time that abundance ratios designed as chemical clocks are even more consistent among the components of wide binaries than their [X/Fe] ratios. Not only that, but the special case of the pair HIP 34426/HIP 34407 may indicate that chemical clocks are consistent for coeval stars even when the individual abundances are not. If the assumption that chemical clocks are reliable age indicators is correct, this would constitute first quantitative, statistically significant evidence that the components of wide binaries in the Galactic field are indeed coeval, validating a large body of published work that relies on that to be the case. Moreover, our results provide strong evidence that chemical clocks indeed carry important information about stellar birthplaces and chemical evolution, and thus we propose that including them in chemical tagging efforts may facilitate the identification of nowadays dissolved stellar groups.

Keywords: stars: abundances- binaries - solar-type - solar neighbourhood

1. INTRODUCTION happens is explained in Nissen(2015). It is based on

arXiv:2105.01096v1 [astro-ph.SR] 3 May 2021 During recent , several studies have found in the fact that different nucleosynthesis processes occur at FGK low- stars that there are some abundance ra- specific stages of the stellar evolution depending on the tios like [Y/Mg], [Y/Al] and [Ba/Mg] that follow a lin- of the stars. One can observe, for example, that ear relation with stellar age (da Silva et al. 2012; Nissen the abundance of elements that are produced in core- 2015; Feltzing et al. 2017; Titarenko et al. 2019; Del- collapse supernovae from massive stars, like α-elements gado Mena et al. 2019). A way to understand why this (Nomoto et al. 2013), has a very different production rate with time compared with those produced in late AGB stages of intermediate mass stars, like neutron- 1 Released on XX, YY, ZZ capture elements (Karakas & Lattanzio 2014). Because Corresponding author: Francisca Espinoza Rojas these two processes enrich the interstellar medium (ISM) [email protected] 2 Espinoza-Rojas et al. at different timescales, we now see a dependence with because the consistency in the detailed chemistry of their stellar age when comparing the two families of elements. components (e.g. abundances and metallicity) imply a Following this idea of different chemical enrichment common origin for these stars. In addition, wide binaries rates for elements of different nucleosynthetic channels, span a wider range of age and [Fe/H] than open clus- Jofr´eet al.(2020) looked for more trends in abundance ters, thus making it possible to understand the possibil- ratios between different nucleosynthetic families using ity of chemical tagging over a more representative frac- a large sample of solar twins analysed by Bedell et al. tion of the Galactic stellar populations. Several studies (2018). They found that 55 different abundance ratios have already analyzed chemical abundances in wide bi- have similar chemical clock-behavior. All these ratios nary systems (e.g. Oh et al. 2017; Ram´ırezet al. 2019; involved neutron-capture elements (s-process produced Hawkins et al. 2020) and found that they are very ho- by AGB mostly). These abundance also have shown to mogeneous (within ∼ 0.1 dex and ∼ 0.05 dex for stars add crucial information in chemical evolution studies for with similar Teff), but also found that there are some metal-poor stars (see Sk´ulad´ottiret al. 2019). cases in which these systems seem to have statistically Although chemical clocks have been the subject of significant differences in chemical composition. If the many studies, their applications and limitations are chemical abundances follow a trend with condensation not yet fully understood. In particular, Feltzing et al. temperature, this can be typically attributed to plan- (2017) challenged the chemical clock applicability of etary formation or planet engulfment (Mel´endez et al. [Y/Mg] by showing that it depends on metallicity. Such 2009; Teske et al. 2015; Tucci Maia et al. 2019), en- dependencies have been further studied for [Y/Mg], hancing the abundances of refractory and/or Fe-peak [Y/Al] by Casali et al.(2020) and [Sr/Mg] by Nissen elements. Hawkins et al.(2020) used a sample of 25 et al.(2020) showing a significant effect. Casamiquela wide binaries such that their components were of similar et al.(2021) used a large sample of open clusters with spectral type and found that odd-z elements, Fe-peak, well-known ages, to show that offsets and dispersion are α-elements, and neutron-capture elements are very con- introduced in the age-abundance relations when reach- sistent, in contrast to what is found when random pairs ing out of the local bubble (d > 1 kpc). Therefore, it is are used. To test the applicability of chemical tagging, not straightforward to apply these relations outside the high precision abundances are key. solar-metallicity solar-vicinity. Such high precision abundance analysis can be done by performing line-by-line differential analyses using a Another test that could be performed to evaluate the with similar stellar parameters as reference, like applicability of chemical clocks would be to study how the for solar-twins. This helps to minimise the consistent they are between stars that were born to- effects of systematic errors in the derivation of the gether, with systems that allow covering a much greater abundances such as approximations in the model atmo- range of ages and metallicities. Wide binary systems spheres and/or inaccuracies in the atomic data, and thus (100AU < a < 1pc) meet these conditions, since their achieve high-precision measurements (see Jofr´eet al. possible formation scenarios, such as the association of 2019; Nissen & Gustafsson 2018, for more details). pairs of stars during the dissolution of their cluster of These levels of precision (∼ 0.01 dex Mel´endezet al. stars birth of a single age (Kouwenhoven et al. 2010; 2009) open the door to new ways of multidimensional Pe˜narrubia 2021), dynamic deployment of triple sys- studies of the chemical evolution of the Milky Way (e.g. tems (Elliott & Bayo 2016), turbulent fragmentation using evolutionary trees Jofr´eet al. 2017a; Jackson et al. (Lee et al. 2017) and gravitational attraction of nearby 2020). prestellar nuclei (Tokovinin 2017), imply that they are The aim of this work is to evaluate how consistent the formed by stars which have a common origin, being an various chemical clocks are among the components of open cluster analog of only two stars. As a consequence, wide binary systems, taking advantage of their expected wide binaries are expected to have similar chemical com- contemporary nature. Section2 explains the data used, positions and ages. To our knowledge, the only work and the process with which we derived the stellar pa- that studies the latter with a large sample of systems rameters and chemical abundances of our sample is de- is that of Kraus & Hillenbrand(2009), where they con- scribed in Section3. The results regarding chemical ho- firmed that young binary stars in Taurus-Auriaga asso- mogeneity are explained in Section4, and in Section5 ciation are coeval within ∼ 0.16 dex in log τ. we present a detailed analysis of the consistency found Originally proposed by Andrews et al.(2018, 2019), in chemical clocks. Finally, our conclusions are in Sec- wide binaries are an ideal sample for studying Galactic tion6. chemical tagging (Freeman & Bland-Hawthorn 2002), Chemical clocks in wide binaries 3

2. DATA For this study we used a sample of 36 wide binary systems in total. Among them, 5 were observed and an- alyzed by us, while the rest was compiled from different sources in the literature. This allowed us to have better statistics in our results.

2.1. Our sample of wide binaries Our 5 wide binaries have solar-type components and have been previously classified as bonafide binaries us- ing common proper motions and photometry (Chanam´e & Gould 2004; Gould & Chanam´e 2004). We checked their parallaxes ($) from and radial velocities from Gaia EDR3 (Gaia Collaboration et al. 2020), to confirm that the pairs belong to a wide binary system. We also com- puted the projected separation of each system, consider- ing the inverse of the mean parallax as distance and cal- Figure 1. One-dimensional spectra of our sample of wide culating their angular separation θ, as done in Andrews binaries, centered in the oxygen triplet lines at ∼777 nm. et al.(2018). Given the consistency in these parameters The primary component is shown in blue and the secondary we can say that they are bonafide wide binaries. in orange. These five systems were observed between October 1 2011 and January 2012 with the high-resolution Mag- For systems HAT-P (Saffe et al. 2017), HD 20782/HD ellan Inamori Kyocera Echelle (MIKE) spectrograph 20781 (Mack et al. 2014), HD 80606/HD 80607 (Liu on the 6.5 m Magellan Clay telescope, located at Las et al. 2018), WASP-94 (Teske et al. 2016), XO-2 Campanas Observatory, which has a resolution of R ∼ (Ram´ırezet al. 2015), HD 134439/HD 134440 (Reggiani 65, 000 for the red, and R ∼ 85, 000 for the blue or- & Mel´endez 2018), and ζ Ret (Saffe et al. 2016), we used ders. In addition, we considered a solar spectrum taken the available EWs published by the authors to derive from VESTA, also observed in that period. The spectra their abundances using the same methodology than for ˚ ˚ cover a range of 3380A to 9370A and have signal-to- our own observations (explained in Section3). noise (S/N) above 250. Data reduction was done using In almost all of these systems, except for the last two, the MIKE Carnegie-Python pipeline. it has been possible to probe that at least one of its Figure1 shows the oxygen triplet of our five wide bi- components hosts planets. Although the existence of naries, where the primary and secondary component are planets in HD 134439/HD 134440 and ζ Ret has not in blue and orange, respectively. Here we see that some been confirmed, they also present a slight difference in components have differences in shape and/or strength their chemical composition, of about 0.06 dex (Reggiani of their lines, which besides being caused by abundance & Mel´endez 2018) and 0.04 dex (Saffe et al. 2016) in differences can be due to line broadening (e.g. in HIP [Fe/H], respectively. This has been interpreted as possi- 58240/HIP 58241). ble signs of engulfment in HD 134439/HD 134440, and the existence of a debris disk around ζ2 Ret (Eiroa et al. 2.2. Wide binaries in Literature 2010). To improve the statistics of our analysis of chemical This last system is very peculiar. Nissen et al.(2020) clocks, we include 32 wide binary systems for which found it is depleted in yttrium compared to other stars high-resolution spectral analysis have been performed of similar ages. In addition to this, their isochrone and in the literature, and abundances or equivalent widths activity ages are significantly lower than the chemical 1 (EWs) are available. More information about these ages derived with [Y/Mg] (9.1 ± 0.5 Gyr for ζ Ret and 2 systems, like stellar parameters, and abun- 9.4 ± 0.5 Gyr for ζ Ret), which implies that its atmo- dances can be found in the tables included as online sphere behaves like that of a young star, despite being material. chemically and kinetically old. A plausible explanation for this if that they are blue straggler stars formed by short period binaries (Rocha-Pinto et al. 2002). 1 CN2012A-063 and CN2012A-064 CNTAC programs 4 Espinoza-Rojas et al.

Table 1. Signal-to-noise ratio, radial velocities, and stellar atmospheric parameters of our sample of five wide binaries derived with iSpec.

Star S/N RV Teff log g [M/H] vmic s (km/s) (K) (dex) (dex) (km s−1) (AU) HIP32871 344 0.25±0.26 5938 ± 82 4.27 ± 0.10 -0.24 ± 0.01 1.21 ± 0.04 HIP32865 342 0.0±0.24 5843 ± 72 4.39 ± 0.09 -0.23 ± 0.01 1.01 ± 0.04 1254.51 HIP34426 506 -12.15±0.28 5970 ± 93 4.13 ± 0.10 -0.55 ± 0.01 1.26 ± 0.06 HIP34407 409 -12.64±0.24 5974 ± 73 4.23 ± 0.10 -0.34 ± 0.01 1.26 ± 0.04 8218.91 HIP52792 332 29.63±0.42 6275 ± 94 4.19 ± 0.11 -0.43 ± 0.02 2.03 ± 0.05 2594.17 HIP52793 355 28.99±0.28 5994 ± 91 4.17 ± 0.11 -0.47 ± 0.02 1.33 ± 0.04 HIP58240 299 6.21±0.22 5848 ± 65 4.52 ± 0.08 -0.01 ± 0.02 1.41 ± 0.04 664.31 HIP58241 297 7.0±0.32 5901 ± 92 4.74 ± 0.06 0.03 ± 0.02 2.00 ± 0.08 HIP15304 310 31.54±0.24 6120 ± 42 4.26 ± 0.08 0.22 ± 0.09 1.62 ± 0.08 7374.10 HIP15310 288 32.18±0.24 5966 ± 66 4.50 ± 0.07 0.25 ± 0.08 1.35 ± 0.03 VESTA 255 -5.26±0.2 5799 ± 66 4.49 ± 0.09 -0.01 ± 0.07 1.07 ± 0.04

p 2 2 Note— s = θ × 1/$, θ = (αA − αB ) cos δA cos δB + (δA − δB ) [M/H] = [Fe/H] + log(0.638fα + 0.362), log fα = [α/H]

The abundances of the remaining 25 wide binaries The bottom panel shows a kiel diagram (Teff vs log g) were taken from Hawkins et al.(2020, hereafter H20), colored by metallicity, where the dashed lines point the where they derived parameters and abundances with location of the Sun. Each binary is connected with a spectral synthesis using the BACCHUS code (Masseron black line, which allows us to see that the components et al. 2016). We noted that their system called WB01 is have very similar colors, absolute magnitudes and stellar in fact HD 80606/HD 80607(Liu et al. 2018). In order parameters. to avoid counting it twice, we only kept its data from Most of the pairs are well placed in the expected main H20. sequence locus, except for HD 134439/HD 134440, which There are no studies that confirm or refute the exis- is located below. The components of this binary are very tence of planets in our systems or those of H20 (except metal-poor(−1.43 dex and −1.39 dex, respectively), so for WB01 and HIP 34426/HIP 34407 (Ram´ırez et al. they are much less affected by blanketing, and there- 2019)). Despite this, in H20 they mention that the fore, bluer than a star with the same mass but higher condensation temperature trend of the 5 systems with metallicity (i.e. solar). In Reggiani & Mel´endez(2018) ∆[Fe/H] > 0.1dex was studied and that, in some cases, they found that its chemical composition and kinemat- there is indeed a clear correlation indicating accretion ics are not consistent with low-α nor high-α halo stars, of rocky material, but they do not go further into the but instead suggest that it might have formed in a dwarf subject. similar to . This possible intergalactic ori- The advantage of including some binary systems gin makes it an interesting object to test chemical clocks. whose components present variations in their chemical patterns, is that it will allow us to test the consistency of chemical clocks even when there have been additional 3. CHEMICAL ABUNDANCES OF WIDE factors that could alter the chemical composition of the BINARIES stars, such as primordial differences (i.e. planet forma- 3.1. Our sample tion and engulfment) or stellar evolution. The top panel of Figure2 shows our entire sample of The analysis of the MIKE spectra was done using the stars in Gaia EDR3 color-magnitude diagram (CMD) spectroscopic framework iSpec (Blanco-Cuaresma et al. divided into 3 groups: our five wide binaries (green), 2014; Blanco-Cuaresma 2019), which is a wrapper ca- H20 systems (orange), and the rest of the systems from pable to run several state-of-the-art radiative transfer literature (blue). For reference we show in grey circles codes and model atmospheres. iSpec has also imple- the catalog of wide binaries from El-Badry & Rix(2018). mented several functions for basic spectra manipulation and computation of radial velocities. Chemical clocks in wide binaries 5

3.1.2. Stellar parameters Both abundances and stellar parameters were deter- mined from the observed EW, which where measured with iSpec fitting a Gaussian profile to each line. For the abundance determination, we selected the lines whose strength falls mainly in the linear part of the curve of growth (Gray 2008). In all our procedure we use the MARCS 1D models (Gustafsson et al. 2008) and LTE2 MOOG radiative transfer code (Kurucz 1993; Sbordone et al. 2004). The stellar atmospheric parameters were determined with iSpec using the classic spectroscopic method: excitation-ionization balance. This method requires no correlation between excitation potential and iron abun- dance derived from several absorption lines (excitation balance). In addition, it requires having same abun- dances derived form neutral and ionized iron lines (ion- ization balance). That way it is possible to constrain the effective temperature and surface gravity, respectively. Finally, the microturbulence velocity is obtained impos- ing a flat slope in the iron abundance as a function of the reduced equivalent width (REW) of each line. We used the extended selection of lines provided by iSpec that covers a range of 480 nm − 615 nm with 154 neu- tral and 11 singly-ionized iron lines. Atomic data was taken from the fifth version of the linelist of atomic data from Gaia-ESO Survey (GES) (Heiter et al. 2020) and we considered the solar abundance from Grevesse et al. Figure 2. Top: Gaia DR3 color-magnitude diagram of the stars used in this work. Background stars (grey) were taken (2007). Finally, we discarded outlier lines using a ro- from El-Badry & Rix(2018) catalog of wide binaries. Bot- bust regression of 0.9 with iSpec, to prevent any strong tom: Kiel diagram (Teff vs log g) colored by metallicity. The dispersion in the data. Our resulting stellar parameters dashed lines show the position of the Sun. . HD 134439/HD are listed in Table1. 13440 is separated from the rest of the sample and lies be- The system HIP 34426/HIP 34407 is the only one low the main sequence, given that it is very metal-poor and of our sample that has been spectroscopically analysed suffers much less blanketing in comparison with other more in the literature. Ram´ırez et al.(2019) determined metal-rich stars of the same mass. their stellar parameters using a differential approach with the q2 code3 and obtained T = 5925 ± 12 K, 3.1.1. Data preparation eff log g = 4.33 ± 0.03 dex and [Fe/H] = -0.367 ± 0.010 dex For each order its continuum was fitted by applying a for HIP 34407, and Teff = 6007±8 K, log g = 4.32±0.02 median-maximum filter to the fluxes in order to reduce dex and [Fe/H] = -0.544 ± 0.006 dex for HIP 34426. the effects of noise and ignore the deeper fluxes of ab- These values are in good agreement with our results sorption lines. Then, it was modeled with a third-degree within the uncertainties. Also, they estimated that their spline function every 5nm ignoring prefixed strong lines ages are 6.5 ± 0.8 Gyr and 6.6 ± 0.8 Gyr, respectively, like tellurics and strong absorption lines identified and assuming that both stars have identical internal metal- masked using a telluric line list from iSpec. After licities of about -0.28 dex. normalization, the orders were stacked using scombine task from IRAF using the median value of the fluxes. 3.2. Chemical abundances The one-dimensional spectra were corrected by their barycentric velocities with iSpec before determining the radial velocities from cross-correlation with a high S/N 2 Local Thermodynamic Equilibrium solar spectrum from NARVAL (370−1048 nm) provided 3 The Qoyllur-Quipu (q2) code(Ram´ırez et al. 2014) https:// by this framework. github.com/astroChasqui/q2 6 Espinoza-Rojas et al.

Chemical abundances were determined from EWs us- certainties. The final uncertainty is obtained by adding ing an extended selection of lines from different chemical in quadrature the element sensitivity and the standard species provided by iSpec. This list is adapted for the error. model atmospheres, radiative transfer code and linelist p 2 2 we use (see Section 3.1.2). σtot = SE + σsens (2) The linelist and the radiative transfer code include a treatment for hyperfine structure splitting (HFS) for Sc 3.3. Abundances of wide binaries from the literature I, V I Mn I, Co I, Cu I, Ba II, Eu II, La II, Pr II, Nd The abundances of systems collected from the litera- II, Sm II. That way we measured abundances of 23 dif- ture (except H20) were determined from their EW list, ferent elements from the following nucleosynthetic fam- following a similar method as ours (Section 3.2). In ad- ilies: carbon, α (O, Mg, Si, Ca, Ti), neutron-capture dition, the authors of the analyses of the seven systems (Sr, Y, Zr, Ba, La, Ce), odd-Z (Na, Al, Sc, V, Cu), taken from the literature reported the EWs of the ref- and iron-peak elements(Cr, Mn, Fe, Co, Ni, Zn). We erence star that they used (usually VESTA, except for use the same nucleosynthetic classification as in Jofr´e HD 134439/HD 134440 where they used HD 103095). et al.(2020). Oxygen abundances were determined us- We used this information to perform the differential ing the oxygen triplet at 777 nm, and non-LTE correc- analysis in the same way as for our systems. Something tions where extracted from M. Kovalev et al.(2018). important to recall is that each of these systems were Carbon abundances were determined using C I lines. observed within a different spectral range, therefore Once individual line absolute abundances are mea- their abundances were determined with different lists of sured, strictly line-by-line differential abundances can lines. be computed. For the element X, the differential abun- dance ∆X is defined as the mean value of the difference Our final sample of wide binaries consists of 36 sys- in abundance of each line with respect to the same lines tems (72 stars) with stellar parameters that span a range in the reference star: of ∼ 4946 − 6682 K in Teff, ∼ 3.79 − 4.83 dex in log g

N and ∼ −1.43 − 0.41 dex in [Fe/H]. 1 X ∆X = (AX,i − AX,ref) (1) N 4. CHEMICAL HOMOGENEITY OF WIDE i=1 BINARIES Where AX,i is the absolute abundance of the ith-line of the given star, AX,ref is the absolute abundance of The chemical homogeneity of stars in wide binary sys- the same line measured in the reference star, and N is tems has been previously studied (e.g. Andrews et al. the number of lines. Here, the reference star is the Sun, 2018; Hawkins et al. 2020), and these components have whose abundances are obtained from the spectrum of generally been found to exhibit very similar chemical VESTA. patterns, supporting evidence that they have a common origin. Uncertainties in element abundances√ were obtained using the standard error (SE = σ/ N) when more than three lines were available. This error estimation 4.1. Chemical homogeneity in our systems method is frequently used (e.g. Mel´endezet al. 2012; The abundance differences ∆X between the compo- Ram´ırezet al. 2015; Adibekyan et al. 2016), but it is nents of each of our systems are shown in Figure3. important to recall that in the case of elements with few In order to help to visualize the precision of the mea- lines it can underestimate the real dispersion, and it is surements, differential abundances with uncertainties recommended to use the median with the interquartile greater than 0.05 dex are colored in grey, whereas or- (q75 − q25) range as its uncertainty, because it is a more ange points depict differences with smaller uncertainties. robust estimator of the abundances. (see Jofr´eet al. Each panel also includes the difference in Teff, log g and 2018). [Fe/H] between the binary components. For elements with only one spectral line available we In some cases we were not able to measure all the el- used a conservative error of 0.03, which is the max- ements listed in Figure3 because the same spectral line imum error for the elements with more than 3 lines. was not correctly measured in one or both components Also, we take into account the effects of the uncertain- of that system. Specifically, this is the case of Zr, Sr I, ties in stellar parameters by computing the sensitivity Sc I, La II and Ce II, most of them n-capture elements σsens for each element. Following Casamiquela et al. which are known to be hard to measure since their lines (2020) this σsens is the dispersion in abundances ob- are weak and sometimes blended (Jofr´eet al. 2019). In tained when varying the stellar parameters by their un- general, ∆X is very close to zero (in average ∼ 0.1 dex), Chemical clocks in wide binaries 7

Figure 3. Differential chemical abundances of wide binary pairs with solar-twin components. Each panel includes their difference in Teff, log g and [Fe/H]. Grey points are abundances with uncertainties larger than 0.05 dex, whereas orange points depict uncertainties smaller than 0.05 dex. HIP 32871/HIP 32865 is a remarkable case of homogeneity, while HIP 34426/HIP 34407 has large abundance differences of 0.2 dex in average. bearing out that wide binaries systems are indeed chem- 4.2. Chemical homogeneity of the whole sample ically homogeneous, in agreement with the literature. Figure4 shows a comparison of the [X/Fe] for the 24 HIP 32871/HIP 32865 is a remarkable case of ho- chemical species, for the full dataset of wide binaries. mogeneity whereas HIP 34426/HIP 34407 has average The full sample analyzed in this paper is depicted as abundance differences of 0.2 dex in most of the elements, follows: our wide binaries (green), H20 (orange) and being the most inhomogeneous system in our sample. the rest of the systems from the literature (blue). Each This system was also studied in Ram´ırezet al.(2019), panel shows a different element, which is indicated at where they proposed and analysed several scenarios that its bottom right. [Fe/H] is included in the bottom right could be causing this difference in their chemical com- panel. The abundances are compared considering one position: (i) perhaps they are not binary from birth, (ii) component in the X axis and one component in Y axis, engulfment of planetary material rich in refractories, and expecting to follow a 1-1 line if the abundances were to (iii) inhomogeneities in the gas cloud where they were match exactly. formed. Despite the efforts, it has not yet been possible In order to quantitatively study the degree of ho- to confirm which of these scenarios led to this peculiar mogeneity of the abundances in wide binaries we com- system. pare them with the chemical behaviour of random pairs, Regarding how consistent wide binaries are in the which are not expected to follow the 1-1 line as wide different nucleosynthetic families, we find that the mean binaries, given that their components were not formed ∆X are ∼0.02 dex for neutron-capture (Sr, Y, Zr, Ba, together. We constructed random pairs selecting stars La, Ce), ∼0.003 dex for α (O, Mg, Si, Ca, Ti), ∼0.03 according to their position in the CMD (Figure2) as dex for odd-Z (Na, Al, Sc, V, Cu), and ∼0.03 dex for follows. We set a star as a reference and choose the two Fe-peak elements (Cr, Mn, Fe, Co, Ni, Zn). This is in stars (which are not its true companion) closest to it in good agreement with Hawkins et al.(2020). the diagram. This is done because our sample has stars that were analyzed with different lists of lines, so it is very likely that by only selecting the closest star, there 8 Espinoza-Rojas et al.

Figure 4. Comparison of the individual abundances [X/Fe] of the two components of wide binaries for all the elements measured in this work. [Fe/H] is in a separated panel at the bottom right of the figure. The systems are coloured as: binaries in our sample (green), H20 (orange) and other binaries from literature (blue). The distribution of random pairs is shown as grey contour lines. If there are less than 3 random pairs with [X/Fe], their distribution is not shown. HIP 34426/HIP 34407, the most inhomogeneous system of our sample, is marked with a black square. will not be an absolute match between the lines of the II and Co I, which span a wider range in [X/Fe], reach- components of the random systems. This could con- ing values above 0.5 dex. Regarding random pairs, they siderably reduce the number of these systems in each also tend to be concentrated within this range, except [X/Fe] and affect our statistics. This exercise is done in some cases (Mn I, Sc I, etc). for each of the 72 stars of the whole sample, creating a We further see that some elements show systematic total of 144 random systems. Grey contours show the differences among the three samples in the sense that distribution of the random pairs if at least three of them they seem to cover different abundance regimes, e.g. Mn have measurements of that [X/Fe]. We note H20 does I, Sc I, Ba II, Co I. The reason behind this might be not have abundances of O or Ce. That is why some of that each of these groups were analysed in different ways these panels, in addition of lacking orange points, lack (see Section 2.2). Different methodologies, even when grey contours. using the same prescription and input data, can lead In Figure4 we can see that wide binaries are generally to systematic differences in abundances. An extensive more tightly distributed along the 1-1 line than random discussion in this subject can be found in Jofr´eet al. pairs, as expected from previous works (Andrews et al. (2017b). 2018; Desidera et al. 2006). In general, abundances of According to the excitation-ionization balance, neu- almost all of the elements show solar abundance pat- tral and ionized abundances of an element should be terns and follow a similar distribution centered at zero, equal, but in practice this is not always the case given possibly because these stars belong to the solar neigh- non-LTE, 3D effects or HFS (see Jofr´e et al. 2015). bourhood. Exceptions to this are Ba II, Mn I, Cu I, La We have 3 elements with abundances for the neutral Chemical clocks in wide binaries 9 and ionized species: Sc, Cr, Ti. Previous studies have The linear regression that better fits the data is plot- found that Sc I and Cr I are strongly affected by non- ted as a blue line, and its parameters are included in each LTE (Zhang et al. 2008; Bergemann 2010; Battistini & panel. We leave the other sample of systems taken from Bensby 2015), while for Ti I the corrections are minor, literature (blue systems in Figure2) out of this plot, be- but it is still more sensitive to non-LTE (Bodaghee et al. cause the presence or absence of trend has already been 2003). It is for this reason that we decided to only use studied in their respective works, and we found their the Sc II, Cr II and Ti II measurements as the abun- slopes to be consistent with previous results. dances of each element for the rest of the following anal- In H20 they briefly discussed about this trend in the ysis. systems with ∆[Fe/H] ≥ 0.1dex, which we identify as These results and those found in previous studies re- WB02, WB05, WB09, WB16 and WB21. Here, we garding the similar chemical composition of wide bina- find that three of them, WB02, WB05, and WB21, ries (within 0.1 dex), serve as sufficient evidence that present a considerable trend, but also WB10, WB11, these systems are excellent laboratories for chemical tag- WB22 and WB24. Regarding our sample, we find that ging, studying effects such as planet presence, and also HIP 34426/HIP 34407 and HIP 58240/HIP 58241 also to test the consistency of chemical clocks, since they sup- show an important trend with condensation tempera- port the premise that these stars were born together. ture. Ram´ırezet al.(2019) studied this trend in HIP 34426/HIP 34407 and found a similar behaviour. The particular chemical inhomogeneity of this system causes 4.3. Condensation temperature trend it to have the highest slope among all the wide binaries The formation of planets can leave a subtle imprint in in this study, for which one would point to the forma- their host’s stellar atmosphere (∼ 0.01 dex) (Chambers tion or accretion of rocky planets as the main responsible 2010). While in general planets are formed from the for their nature (as discussed in Ram´ırezet al.(2019)). same material as their host star, rocky planets in par- Something similar could be proposed for the other sys- ticular would preferentially remove refractory content tems. Despite this, there is still no way to know with (higher Tc), so they may leave the stars slightly poor certainty which specific mechanism is causing this cor- in these elements (Mel´endezet al. 2009; Ram´ırezet al. relation between chemical difference and Tc just with 2009; Nissen 2015). Another scenarios in which the stel- the information that we handle. Probably, a systematic lar composition can be altered are planet engulfment or search for planets in this system would help understand- accretion of material with planet-like composition (e.g. ing if the slope is related to this effect. Saffe et al. 2017; Tucci Maia et al. 2019), thus given the nature of these cases and following the same logic as be- 5. CHEMICAL CLOCKS fore, they should also leave evidence in the abundance of refractory elements of the star In addition to common chemistry, that we have ev- Figure5 shows ∆[X/H] as a function of condensation idenced in the previous section, the expected common origins of the components of wide binaries imply that temperature, Tc (Lodders 2003), for H20 and the wide binaries in our sample that present a trend with an ab- they should have similar ages as well (e.g. Kraus & Hil- solute value of the slope greater that 3 × 10−5 dex/K lenbrand 2009). This coevality allows us to test the reli- with a significance greater than 2σ. This cut is based ability of the abundance ratios known as chemical clocks on the slopes that have been reported in the literature, (Section1) with which a trend has been confirmed or denied (Mack We follow Jofr´eet al.(2020), who used a sample of et al. 2014; Saffe et al. 2016; Teske et al. 2016; Liu et al. solar twins to present a large number of combinations 2018; Reggiani & Mel´endez 2018). The linear regression of chemical abundance ratios shown to have strong de- pendency with ages. These relations can be explained was calculated using the function linregress of SciPy library4. Since this function does not include the indi- with similar arguments as [Y/Mg], which is widely dis- vidual error of the elements to derive the uncertainty cussed in the literature (e.g. Nissen 2015; Delgado Mena of the fit, nor does it allow to make a weighted fit as et al. 2019; Casamiquela et al. 2021). The argument is suggested in Adibekyan et al.(2016), we performed the the very different timescales for the production of an linear regression for the abundances with errors smaller alpha-capture element such as Mg in comparison with than 0.05 dex (orange points in Figure5). a neutron-capture element such as Y. When combined, a tight relation of [Y/Mg] can be found for stars of the same metallicity but different ages. Jofr´eet al. 4 https://docs.scipy.org/doc/scipy/reference/generated/scipy. (2020) and Casamiquela et al.(2021) found that the stats.linregress.html clocks showing largest age dependency involved neutron- 10 Espinoza-Rojas et al.

Figure 5. Differential abundance as a function of elemental condensation temperature from the H20 systems and our sample of wide binaries that present a slope greater than 3 × 10−5 dex/K. The blue line shows the linear regression that better fits the data, which parameters are included in the legend of each system. Differential abundances are colored in the same way as in Figure3. capture elements with elements produced by other nu- sample. As well as in Figure4, we believe that this cleosynthesis channels, as expected given their produc- is mostly due to a systematic effect in the derived abun- tion scenarios. dances, as mentioned in Section 5.3. We recovered 42 of the abundance ratios5 found in Jofr´eet al.(2020). Figure6 is similar to Figure4 but 5.1. Consistency of chemical ages in wide binaries we plot chemical clocks instead of [X/Fe] abundance ra- In order to quantify the degree of correlation in our tios. In each axis we plot each component of the binary, data, or lack of thereof, we compute the Pearson cor- expecting the values to follow a 1-1 line. Furthermore, relation coefficient (PCC)6 for both wide binaries and we included the distribution of random pairs, like in Fig- random pairs for clocks and standard abundance ratios. ure4. Since H20 lacks of Ce and O abundances, there In Figure7 we plot the distributions of the PCC from are 8 chemical clocks for which we are unable to have Figures4 and6, showing wide binaries in orange and an estimate of a random pair correlation. random pairs in blue. The coefficients were derived only In addition, we highlight the system HIP 34426/HIP from those chemical clocks and [X/Fe] ratios with at 34407, which is the least homogeneous wide binary of least 10 random pairs, to ensure robust statistics. our sample (see individual results in second panel of Fig- The top panel of Figure7 shows that standard [X/Fe] ure3). This system is marked with a black square, and abundances are more consistent in wide binaries than despite being really inhomogeneous this does not seem random pairs, but the differences between distributions to have an effect on some of its chemical clocks, espe- are not so stark. On the other hand, in the bottom cially those that depend on Ba (e.g. [Ti/Ba], [Al/Ba], panel we can see that the PCC distributions of chemical [Ni/Ba]). Some chemical clocks span a wider range clocks are very different from each other, since random in abundance ratio than others, like [Cu/Ba], [Zn/Ba], pairs show a peak at ∼ 0.1 while in wide binaries this is [Na/Sr], and [O/Y], particularly because of the H20 at ∼ 0.8 (see Table2). Moreover, the areas covered by

√ 5 [X/Y] = [X/H] - [Y/H]; δ[X/Y] = δX2 + δY 2 6 Derived with SciPy: https://scipy.org Chemical clocks in wide binaries 11

Figure 6. Comparison between the components of wide binaries for 42 potential chemical clocks considered in this work. The different systems are coloured like in Figure4 and the distribution of random pairs is shown in grey curves. Same as in Figure4, HIP 34426/HIP 34407 is marked with a black square for easy visual identification. each distribution almost do not overlap, which means fected by different mechanisms (e.g. planet formation, that there is a very low chance that stars that do not planet engulfment, etc), Figure7 helps us to confirm necessarily have a common origin can have the same the idea that chemical clocks are sensitive to the age of chemical clocks, thus age. the stars regardless of the possible differences in their If chemical clocks were just reflecting abundance ra- compositions. tios, their 1-1 correlation should be comparable to the In a different scenario, if we simply consider true that correlation found for other individual abundances be- chemical clocks are capable of identifying contemporary tween coeval stars that were born together. However, systems, our results would confirm that the components our results imply that chemical clocks are even more of wide binaries are coeval, and it would be one of the consistent, therefore, they do not only trace chemical first works in literature to do an analysis with such a homogeneity among wide binaries, they must be carry- large and varied sample of them. ing age information as well. Moreover, since our sample included systems which composition has been subtly af- 12 Espinoza-Rojas et al.

All that said, however, the case of HIP 34426/HIP 34407 is puzzling and is worth discussing. From Fig- ure3 we see that abundances are very different (about 0.2 dex), with C and O the only exceptions. These stars were already analysed by Ram´ırezet al.(2019), who sug- gest the differences might be coming from stars not be- ing truly siblings but chance alignments, the birth cloud not being homogeneous, or engulfment of planetary ma- terial. After performing an analysis of the astrometry and the radial velocities of the system, they found that HIP 34426 shows a variation in radial velocity of ap- proximately 60 m/s, something that could indicate the presence of a planet such as a large Jupiter or the exis- tence of a third component of the system. Furthermore, the difference in the total velocity of the system posi- tions it as a physically unbound system, but given its age, Ram´ırezet al.(2019) propose that it was a binary at formation. If these stars were not siblings, but instead were (1) formed at separate sites with a difference of about 0.2 dex in metallicity, or (2) sites that suffered the possi- ble effects of chemical evolution (e.g. some considerable time span between the formation of both components) or space-phase (e.g. stellar migration), then, unless a large coincidence, they most likely would have significantly different ages and, therefore, different clocks. However, Figure 7. Distribution of the Pearson correlation coeffi- cient (PCC) from Figures4 (top) and6 (bottom), colored according to Ram´ırezet al.(2019) their ages are very in orange for wide binaries and in blue for random pairs. consistent, which most likely rules out these possibili- Statistical features of each plot are listed in Table2. ties and that of a chance alignment, since the proba- bility for two independent and coeval stars to become a wide binary by dynamical capture in the field is neg- Table 2. Median and MAD of the PCC distribu- ligible. Furthermore, we see in Figure6 that within tions of chemical clocks and [X/Fe], for wide binaries the 22 chemical clocks measured in this system, 14 of (WBs) and random pairs in Figure7. them show that the system is located almost perfectly on the one-to-one line despite the seemingly remarkable Median MAD inhomogeneity (for a typical wide binary system) in in- Chemical clocks WBs 0.81 0.06 dividual abundances [X/Fe]. random pairs 0.17 0.18 If we take the similarity of the chemical clocks as an in- [X/Fe] WBs 0.61 0.16 dication for the stars in this binary being coeval, then we random pairs 0.22 0.12 interpreted our results as strong evidence that the cloud material from which they formed was inhomogeneous to these levels, on a scale equivalent to the projected sep- aration of the components of this pair. This may have 5.2. The case of HIP 344026/HIP 34407 important implications for the idea of chemical tagging for Galactic archaeology. The case of binary systems like In previous sections we have checked that stars in this stands in stark contrast with the recent detailed wide binaries have homogeneous chemical compositions study of Kos et al.(2020), who determined chemical within the uncertainties. This confirms expectations abundances of about 300 GALAH stars distributed in and previous results (Andrews et al. 2018; Hawkins et al. the Orion , finding that the nebula is chemically 2020) since all known plausible formation scenarios in- homogeneous (typically within 0.1 dex in each element, dicate that they share common origins, hence they must except Li), at higher level of precision than the difference share their chemical composition, especially if they are observed in the HIP 34426/HIP 34407 system. on somewhat similar evolutionary states. Chemical clocks in wide binaries 13

Recently, Pe˜narrubia(2021) proposed that the forma- It is worth commenting that the stellar temperature tion of ultra-wide binaries (s > 0.1 pc) is via chance has a large impact in the line strength of neutral atomic entrapment of unrelated stars in the tidal streams of lines (this is why we adopted the excitation balance disrupting clusters, and found that most pairs formed method to determine temperature, see Section3). The in tidal streams can easily be disrupted by subtle grav- line strength can have a particular large implication in itational interactions with other objects. Here we pro- determining abundances which have strong effects due pose that HIP 34426/HIP 34407 might be part of a tidal to hyperfine structure splitting (HFS), such as Mn, Co stream. This scenario could explain some of the char- or Ba. The way different methods deal with HFS can acteristics that we observe in this system, such as its indeed be a major source of uncertainty for particular el- coevality and the astrometric differences that in prin- ements like Mn or Co (del Peloso et al. 2005; Jofr´eet al. ciple classify it as a non-physically bound system (see 2017b). For stars of different temperatures, abundances Ram´ırezet al. 2019). of Mn can be very different if different prescriptions for HFS are employed. We also refer to Casamiquela et al. 5.3. Systematic uncertainties and applicability of (2020) for the discussions about striking dependencies chemical clocks with Teff are seen (remarkably for Ba, Na, and Mn). The systematic difference in chemical clocks found here We are aware that different spectral methodologies is therefore a combination of line selection, HFS treat- can lead to significantly different abundance results ment and stellar temperature. (Hinkel et al. 2016; Jofr´eet al. 2017b, 2019; J¨onsson et al. 2018). Our analysis depends on data taken from There are some discussions in literature about the the literature, and therefore can be affected by system- metallicity dependency on chemical clocks (e.g. Feltz- atic uncertainties. ing et al. 2017; Delgado Mena et al. 2019). Casali et al. One of them is that, in spite that we used the same (2020) found that the slopes of [Y/Mg] and [Y/Al] ver- linelist for our analysis for the pairs for which we had sus stellar age change with metallicity (by using steps EWs, H20 derived their chemical abundances using a of ±0.2 dex in [Fe/H]), and related this to the star for- different selection of lines than ours. Jofr´eet al.(2019) mation history of the Milky Way. They concluded that extensively discusses the effects of line selections in re- the position of the stars in the Galactic plane is an- sulting abundances. other important parameter to take into account when In particular, while being consistent among all pairs in using chemical clocks to determine ages (e.g. whether it our full sample, we can see from Figure6 that the ratios belongs to the thin disk or to the thick disk), because [Zn/Ba], [Mg/Ba], [Mn/Ba], and [Cu/Ba] have values at different locations stars can have different metallic- very different for H20 than the rest of the sample: their ities. Most recently, Casamiquela et al.(2021) used a values show a tendency to negative values, that do not sample of red clump stars in 47 open clusters to study overlap with range of abundances found in Jofr´eet al. the abundance-age relation in a large spatial volume (2020). We note that Jofr´eet al.(2020) used a sample (6 kpc < R < 12 kpc), covering regions outside the of stars spanning very young to very old ages, leaving CG local bubble, unlike most studies on chemical clocks that little space for a systematic offset in these abundance use solar twins. By studying [Y/Mg] and [Y/Al] in more ratios due to age. detail, they find that the clusters located outside the so- To investigate this further, in Figure8 we plot the lar neighborhood appear to have a significantly greater same abundances as in Figure6 but with the pairs dispersion in abundances compared to the local ones, coloured by the average T of the system. It is evi- eff without finding any dependence on metallicity because dent from the plot that the cases of [X/Y] with large of the narrow range in [Fe/H] of their sample. They at- range along the 1-1 line, also show a gradient with stel- tribute this finding to galactic dynamics processes, such lar temperature, i.e. mass. as radial migration of stars or clusters, which are more Generally, stars with higher temperatures (∼ 6000 − likely to be sampled when covering a wider range in 6200 K) tend to produce ratios concentrated at lower ab- distance and can trace different chemical evolutionary solute values than cooler stars. We particularly remark histories. the cases where Ba is involved (e.g. [Ti/Ba], [Si/Ba], In general, the results of other studies indicate that [Na/Ba]), showing that there probably is some strong chemical clocks have several limitations caused by the effect depending on the stellar mass for the retrieval of chemical and dynamic nature of our Milky Way, so one abundances of this element. Moreover hotter stars seem must be very careful when using them to derive ages. to be more dispersed around the 1-1 line. We recall that Even so, their conclusions do not present problems or re- the hottest stars are also those analyzed by H20. 14 Espinoza-Rojas et al.

Figure 8. Comparison between the components of wide binaries for potential chemical clocks considered in this work. Colors depict the mean temperature of the system. Generally, stars with higher temperatures (∼ 6000−6200 K) tend to produce ratios concentrated at lower absolute values than cooler stars, and hotter stars seem to be more dispersed around the 1-1 line. strictions regarding chemical clocks being consistent in We can see that despite coming from independent coeval star systems. Little however has been discussed spectroscopic measurements, the chemical clocks of the in terms of the effects in deriving ages using chemical wide binaries fit nicely to the 1-1 line, and although we clocks when such abundances are hampered by system- are aware that there are systematic errors in our sample atic uncertainties such as those discussed here. More due to the inclusion of systems analyzed with different such analyses are needed in order to fully assess the ap- methods, these do not influence when plotting the age plicability of chemical clocks for wide Galactic archae- consistency of these systems. But, if they are to be im- ology studies. plemented in chemical tagging studies, it is necessary to follow a spectroscopic analysis as consistent as possible 5.4. Most consistent clocks in wide binaries in order to avoid any systematic uncertainty which could Based on our results we have selected three chemical make the identification of stellar groups very difficult. clocks that have the highest correlation coefficient for It is worth emphasizing that Figure9 shows the po- wide binaries, and that also trace the age consistency tential of these abundance ratios to represent the ages of HIP 34426/HIP 34407 system despite the chemical of systems that arise from a common origin, even when differences of its components. Figure9 shows [Sc/Ba], certain chemical inhomogeneities are present, but still, [Al/Ba] and [Ti/Ba] and the PCC values for wide bina- it remains to be seen how these clocks relate to the in- ries and random pairs. dividual ages of the components of wide binaries. Chemical clocks in wide binaries 15

Figure 9. Chemical clocks with the higher correlation in our sample of wide binaries: [Sc/Ba], [Al/Ba] and [Ti, Ba]. The distribution of random pairs is shown in grey curves and HIP 34426/HIP 34407 is marked with a black square, as in Figure6. These abundance ratios can trace the coevality of HIP 34426/HIP 34407, despite its chemical inhomogeneity.

6. CONCLUSIONS Also, we compared the degree of correlation of indi- In the present work we studied how consistent are cer- vidual abundances and chemical clocks, and found that tain group of abundance ratios, collectively known as chemical clocks are more consistent between components chemical clocks, among the components of wide bina- of wide binaries than individual abundances, even for ries. Since the latter are expected to share a common stars that might have gone through any process that al- formation site and to be coeval, the comparison of their ter its atmospheric composition. Therefore, we conclude chemical clocks provides a basic check on the concept that chemical clocks really carry more information than itself of a chemical clock and its overall validity. For just the plain run of individual chemical abundances, this, we spectroscopically analyzed a sample of five sys- information probably related to age and a common for- tems, four of which for the first time, doing a differential mation site. analysis with respect to the Sun and using the equiva- Finally, we propose that including them in the chem- lent widths method. We found that the components of ical tagging parameter space would be very beneficial these systems are chemically homogeneous (within 0.1 and allow better identification of members of disrupted dex), in good agreement with what has been found in stellar groups. other wide binaries. One remarkable exception of this is HIP 34426/HIP 34407, which presents differences of ∼ 0.2 dex in abundance of most of the elements mea- sured. We found their chemical clocks to be consistent, in agreement with the similarity in their ages that is ACKNOWLEDGMENTS reported in previous works. With this evidence we dis- This project is framed in the project “Gaia sibling carded some possible formation scenarios which could stars as probes of the chemodynamical evolution of explain its chemical difference, like stellar migration or the Milky Way disk”, funded by ECOS-ANID grant a chance alignment. number 180049. F.E-R. acknowledges Iv´anRam´ırez, In order to attempt to test chemical clocks with a big- and appreciates the invaluable guidance and support ger and more representative sample of stars, we enlarged of Marcelo Tucci-Maia, Claudia Aguilera and Danielle our sample by including 31 systems from the literature, de Brito. P.J. acknowledges financial support from finding that the components of these wide binaries fol- FONDECYT Iniciaci´onGrant Number 11170174 and low a 1-1 relation in all the abundance ratios studied, as Regular grant Number 1200703. J.C. acknowledges sup- what is expected thanks to what we know about the for- port from CONICYT project Basal AFB-170002, and mation scenarios of these systems. Our results indicate from the Agencia Nacional de Investigaci´ony Desarrollo that systematic uncertainties from different methodolo- (ANID), via Proyecto FONDECYT Regular 1191366. gies can affect the final abundances, and probably have LC acknowledge support from ”Programme National a significant effect on the different chemical clocks. In de Physique Stellaire” (PNPS) and from the ”Pro- order to better assess and quantify this effect, a more gramme National Cosmologie et ” (PNCG) of careful analysis of common stars but different methods CNRS/INSU. should be performed, as well as having the ages of all these stars. 16 Espinoza-Rojas et al.

APPENDIX

A. CATALOG

Table 3. Summary of the whole sample of wide binaries used in this work. Astrometric data is taken from Gaia DR3, except for radial velocities. The table is available as online material.

Column Units Description system System name Reference Paper reference where the data was taken from Instrument Spectrograph used to obtained the data Spectral range nm Spectral range of the instrument used star Name of the star Teff K Effective temperature eTeff K Effective temperature error logg dex Surface gravity elogg dex Surface gravity error [Fe/H] dex Metallicity e[Fe/H] dex Metallicity error vmic km/s Microturbulence velocity evmic km/s Microturbulence velocity error s AU Projected physical separation of the two stars ra deg ra error deg Right ascension error dec deg dec error deg Declination error parallax mas Parallax parallax error mas Parallax error pmra mas/yr in right ascension direction pmra error mas/yr Standard error of proper motion in right ascension direction pmdec mas/yr Proper motion in declination direction pmdec error mas/yr Standard error of proper motion in declination direction dr2 radial velocity km/s Spectroscopic barycentric radial velocity from Gaia DR2 dr2 radial velocity error km/s Standard error of spectroscopic barycentric radial velocity from Gaia DR2 phot g mean mag mag G-band mean magnitude (Vega scale) phot bp mean mag mag Integrated BP mean magnitude phot rp mean mag mag Integrated RP mean magnitude

REFERENCES

Adibekyan, V., Delgado-Mena, E., Figueira, P., et al. 2016, Andrews, J. J., Chanam´e,J., & Ag¨ueros,M. A. 2018, A&A, 591, A34, doi: 10.1051/0004-6361/201628453 MNRAS, 473, 5393, doi: 10.1093/mnras/stx2685

Andrews, J. J., Anguiano, B., Chanam´e,J., et al. 2019, Battistini, C., & Bensby, T. 2015, A&A, 577, A9, ApJ, 871, 42, doi: 10.3847/1538-4357/aaf502 doi: 10.1051/0004-6361/201425327 Chemical clocks in wide binaries 17

Bedell, M., Bean, J. L., Mel´endez,J., et al. 2018, ApJ, 865, Gustafsson, B., Edvardsson, B., Eriksson, K., et al. 2008, 68, doi: 10.3847/1538-4357/aad908 A&A, 486, 951, doi: 10.1051/0004-6361:200809724 Bergemann, M. 2010, arXiv e-prints, arXiv:1009.5693. Hawkins, K., Lucey, M., Ting, Y.-S., et al. 2020, MNRAS, https://arxiv.org/abs/1009.5693 492, 1164, doi: 10.1093/mnras/stz3132 Blanco-Cuaresma, S. 2019, MNRAS, 486, 2075, Heiter, U., Lind, K., Bergemann, M., et al. 2020, arXiv doi: 10.1093/mnras/stz549 e-prints, arXiv:2011.02049. Blanco-Cuaresma, S., Soubiran, C., Heiter, U., & Jofr´e,P. https://arxiv.org/abs/2011.02049 2014, A&A, 569, A111, Hinkel, N. R., Young, P. A., Pagano, M. D., et al. 2016, doi: 10.1051/0004-6361/201423945 ApJS, 226, 4, doi: 10.3847/0067-0049/226/1/4 Bodaghee, A., Santos, N. C., Israelian, G., & Mayor, M. Jackson, H., Jofre, P., Yaxley, K., et al. 2020, arXiv 2003, A&A, 404, 715, doi: 10.1051/0004-6361:20030543 e-prints, arXiv:2011.06453. Casali, G., Spina, L., Magrini, L., et al. 2020, A&A, 639, https://arxiv.org/abs/2011.06453 A127, doi: 10.1051/0004-6361/202038055 Jofr´e,P., Das, P., Bertranpetit, J., & Foley, R. 2017a, Casamiquela, L., Tarricq, Y., Soubiran, C., et al. 2020, MNRAS, 467, 1140, doi: 10.1093/mnras/stx075 A&A, 635, A8, doi: 10.1051/0004-6361/201936978 Jofr´e,P., Heiter, U., & Soubiran, C. 2019, ARA&A, 57, Casamiquela, L., Soubiran, C., Jofr´e,P., et al. 2021. 571, doi: 10.1146/annurev-astro-091918-104509 https://arxiv.org/abs/2103.14692 Jofr´e,P., Heiter, U., Tucci Maia, M., et al. 2018, Research Chambers, J. 2010, Terrestrial Planet Formation, ed. Notes of the American Astronomical Society, 2, 152, S. Seager, 297–317 doi: 10.3847/2515-5172/aadc61 Chanam´e,J., & Gould, A. 2004, ApJ, 601, 289, Jofr´e,P., Jackson, H., & Tucci Maia, M. 2020, A&A, 633, doi: 10.1086/380442 L9, doi: 10.1051/0004-6361/201937140 da Silva, R., Porto de Mello, G. F., Milone, A. C., et al. Jofr´e,P., Heiter, U., Soubiran, C., et al. 2015, A&A, 582, 2012, A&A, 542, A84, doi: 10.1051/0004-6361/201118751 A81, doi: 10.1051/0004-6361/201526604 del Peloso, E. F., Cunha, K., da Silva, L., & Porto de Jofr´e,P., Heiter, U., Worley, C. C., et al. 2017b, A&A, 601, Mello, G. F. 2005, A&A, 441, 1149, A38, doi: 10.1051/0004-6361/201629833 doi: 10.1051/0004-6361:20053077 J¨onsson,H., Allende Prieto, C., Holtzman, J. A., et al. Delgado Mena, E., Moya, A., Adibekyan, V., et al. 2019, 2018, AJ, 156, 126, doi: 10.3847/1538-3881/aad4f5 A&A, 624, A78, doi: 10.1051/0004-6361/201834783 Karakas, A. I., & Lattanzio, J. C. 2014, PASA, 31, e030, Desidera, S., Gratton, R. G., Lucatello, S., & Claudi, R. U. doi: 10.1017/pasa.2014.21 2006, A&A, 454, 581, doi: 10.1051/0004-6361:20064896 Kos, J., Bland-Hawthorn, J., Buder, S., et al. 2020, arXiv Eiroa, C., Fedele, D., Maldonado, J., et al. 2010, A&A, 518, L131, doi: 10.1051/0004-6361/201014594 e-prints, arXiv:2011.02485. https://arxiv.org/abs/2011.02485 El-Badry, K., & Rix, H.-W. 2018, MNRAS, 480, 4884, doi: 10.1093/mnras/sty2186 Kouwenhoven, M. B. N., Goodwin, S. P., Parker, R. J., Elliott, P., & Bayo, A. 2016, MNRAS, 459, 4499, et al. 2010, MNRAS, 404, 1835, doi: 10.1093/mnras/stw926 doi: 10.1111/j.1365-2966.2010.16399.x Feltzing, S., Howes, L. M., McMillan, P. J., & Stonkut˙e,E. Kraus, A. L., & Hillenbrand, L. A. 2009, ApJ, 704, 531, 2017, MNRAS, 465, L109, doi: 10.1093/mnrasl/slw209 doi: 10.1088/0004-637X/704/1/531 Freeman, K., & Bland-Hawthorn, J. 2002, ARA&A, 40, Kurucz, R. 1993, SYNTHE Spectrum Synthesis Programs 487, doi: 10.1146/annurev.astro.40.060401.093840 and Line Data. Kurucz CD-ROM No. 18. Cambridge, 18 Gaia Collaboration, Brown, A. G. A., Vallenari, A., et al. Lee, J.-E., Lee, S., Dunham, M. M., et al. 2017, Nature 2020, arXiv e-prints, arXiv:2012.01533. Astronomy, 1, 0172, doi: 10.1038/s41550-017-0172 https://arxiv.org/abs/2012.01533 Liu, F., Yong, D., Asplund, M., et al. 2018, A&A, 614, Gould, A., & Chanam´e,J. 2004, ApJS, 150, 455, A138, doi: 10.1051/0004-6361/201832701 doi: 10.1086/381147 Lodders, K. 2003, ApJ, 591, 1220, doi: 10.1086/375492 Gray, D. F. 2008, The Observation and Analysis of Stellar M. Kovalev, S. Brinkmann, M. Bergemann, & MPIA Photospheres IT-department. 2018, NLTE MPIA web server, [Online]. Grevesse, N., Asplund, M., & Sauval, A. J. 2007, SSRv, Available: http://nlte.mpia.de Max Planck Institute for 130, 105, doi: 10.1007/s11214-007-9173-7 Astronomy, Heidelberg. 18 Espinoza-Rojas et al.

Mack, Claude E., I., Schuler, S. C., Stassun, K. G., & Ram´ırez,I., Khanal, S., Aleo, P., et al. 2015, ApJ, 808, 13, Norris, J. 2014, ApJ, 787, 98, doi: 10.1088/0004-637X/808/1/13 doi: 10.1088/0004-637X/787/2/98 Reggiani, H., & Mel´endez,J. 2018, MNRAS, 475, 3502, Masseron, T., Merle, T., & Hawkins, K. 2016, BACCHUS: doi: 10.1093/mnras/sty104 Brussels Automatic Code for Characterizing High Rocha-Pinto, H. J., Castilho, B. V., & Maciel, W. J. 2002, accUracy Spectra. http://ascl.net/1605.004 A&A, 384, 912, doi: 10.1051/0004-6361:20011815 Mel´endez,J., Asplund, M., Gustafsson, B., & Yong, D. Saffe, C., Flores, M., Jaque Arancibia, M., Buccino, A., & 2009, ApJL, 704, L66, Jofr´e,E. 2016, A&A, 588, A81, doi: 10.1088/0004-637X/704/1/L66 doi: 10.1051/0004-6361/201528043 Mel´endez,J., Bergemann, M., Cohen, J. G., et al. 2012, Saffe, C., Jofr´e,E., Martioli, E., et al. 2017, A&A, 604, L4, A&A, 543, A29, doi: 10.1051/0004-6361/201117222 doi: 10.1051/0004-6361/201731430 Nissen, P. E. 2015, A&A, 579, A52, Sbordone, L., Bonifacio, P., Castelli, F., & Kurucz, R. L. doi: 10.1051/0004-6361/201526269 2004, Memorie della Societa Astronomica Italiana Nissen, P. E., Christensen-Dalsgaard, J., Mosumgaard, Supplementi, 5, 93. J. R., et al. 2020, A&A, 640, A81, https://arxiv.org/abs/astro-ph/0406268 ´ doi: 10.1051/0004-6361/202038300 Sk´ulad´ottir, A., Hansen, C. J., Salvadori, S., & Choplin, A. Nissen, P. E., & Gustafsson, B. 2018, A&A Rv, 26, 6, 2019, A&A, 631, A171, doi: 10.1007/s00159-018-0111-3 doi: 10.1051/0004-6361/201936125 Nomoto, K., Kobayashi, C., & Tominaga, N. 2013, Teske, J. K., Ghezzi, L., Cunha, K., et al. 2015, ApJL, 801, ARA&A, 51, 457, L10, doi: 10.1088/2041-8205/801/1/L10 doi: 10.1146/annurev-astro-082812-140956 Teske, J. K., Khanal, S., & Ram´ırez,I. 2016, ApJ, 819, 19, Oh, S., Price-Whelan, A. M., Hogg, D. W., Morton, T. D., doi: 10.3847/0004-637X/819/1/19 & Spergel, D. N. 2017, AJ, 153, 257, Titarenko, A., Recio-Blanco, A., de Laverny, P., Hayden, doi: 10.3847/1538-3881/aa6ffd M., & Guiglion, G. 2019, A&A, 622, A59, Pe˜narrubia, J. 2021, MNRAS, 501, 3670, doi: 10.1051/0004-6361/201833721 doi: 10.1093/mnras/staa3700 Tokovinin, A. 2017, MNRAS, 468, 3461, doi: 10.1093/mnras/stx707 Ram´ırez,I., Khanal, S., Lichon, S. J., et al. 2019, MNRAS, Tucci Maia, M., Mel´endez,J., Lorenzo-Oliveira, D., Spina, 490, 2448, doi: 10.1093/mnras/stz2709 L., & Jofr´e,P. 2019, arXiv e-prints, arXiv:1906.04195. Ram´ırez,I., Mel´endez,J., & Asplund, M. 2009, A&A, 508, https://arxiv.org/abs/1906.04195 L17, doi: 10.1051/0004-6361/200913038 Zhang, H. W., Gehren, T., & Zhao, G. 2008, A&A, 481, Ram´ırez,I., Mel´endez,J., Bean, J., et al. 2014, A&A, 572, 489, doi: 10.1051/0004-6361:20078910 A48, doi: 10.1051/0004-6361/201424244